<div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Bom dia colegas listeiros, como vão?</div><div><br></div><div>Estou com dificuldades em criar um objeto raster a partir de um shapefile ou a partir de um dataframe, mas mantendo os valores associados aos pontos após a transformação.</div><div>Tentei as seguintes abordagens:</div><div><br></div><div><div>r1 = rasterize(shapefile, raster_vazio)</div><div>r1[!<a href="http://is.na">is.na</a>(r1)] = value</div><div>r1 = mask(merge(r1, raster_vazio), raster_vazio, filename = label, format = "GTiff", overwrite = T)</div><div><br></div><div>r2 = rasterFromXYZ(data.frame(coordenadas_pontos, dados_alvo))</div></div><div dir="ltr"><br></div><div dir="ltr"><div dir="ltr">modelo = gstat(formula = 

shapefile$dados_alvo ~ 1,</div><div dir="ltr">                 locations = shapefile,</div><div dir="ltr">                 nmax = 1,</div><div dir="ltr">                 set = list(idp = 0)</div><div dir="ltr">                 )</div><div dir="ltr"><br></div><div dir="ltr">r3 = raster(predict(modelo, as(raster_vazio, 'SpatialPixels')))</div><div dir="ltr"><br></div><div>Dessas abordagens, r1 me retorna um raster vazio; r2 está dando erro, mesmo eu passando um dataframe de grade regular (um ponto a cada 20m) e em r3, o modelo deveria sempre atribuir ao pixel exatamente o valor do ponto mais próximo,</div></div>mas isso não está ocorrendo.</div><div>Caso algum colega possa me elucidar uma nova solução ou apontar qual o erro estou cometendo numa das abordagens listadas, seria de grande ajuda.</div><div><br></div><div>Desde já, agradeço pelo auxilio de todos! </div><div dir="ltr"><br clear="all"><div><div dir="ltr" class="gmail_signature"><div dir="ltr"><div>Yury Duarte<br></div>Engenheiro Agrônomo - ESALQ/USP<br></div></div></div></div></div></div>