Hola a todos, quisiera pedirles su ayuda. Estoy tratando de crear un raster a
partir de una tabla de datos que contiene coordenadas UTM. Determino los
renglones y columnas que debe de tener el raster pero cuando visualizo en un
SIG me dice que no tiene referencia espacial y el tama�o de pixel es de 29.9 y
no de 30 como originalmente era y eso mete ruido al modelo que estoy generando.
Alguien me podr�a ayudar?
xmn=min(table_2[,1]); xmx=max(table_2[,1])
ymn=min(table_2[,2]); ymx=max(table_2[,2])
b <- data.frame(cbind(table_2$x, table_2$y))
coordinates(b) = ~X1 + X2
proj4string(b) <-CRS("+proj=utm +zone=14 +datum=WGS84")
b1 <- spTransform(b,CRS("+proj=longlat"))
data3 <- table_2
data3$long <- b1$X1
data3$lat <- b1$X2
r1 <- raster(nrows=1386, ncols=1649,
xmn=xmn, xmx=xmx,
ymn=ymn, ymx=ymx )
ras1 <- rasterize(table_2[,1:2], r1, field = table_2[,14])
writeRaster(ras1, filename = "Results\\ras2.tif",format= "GTiff",
datatype="FLT4S",overwrite=T)
________________________________
De: R-help-es <[email protected]> en nombre de Francisco
Rodr�guez <[email protected]>
Enviado: mi�rcoles, 19 de febrero de 2020 04:40 p. m.
Para: [email protected] <[email protected]>
Asunto: [R-es] Pregunta sobre rLandsat
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
