[R-es] Pregunta sobre rLandsat

[Angel Cervantes]

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 <r-help-es-bounces using r-project.org> en nombre de Francisco Rodr�guez <fjroar using hotmail.com>
Enviado: mi�rcoles, 19 de febrero de 2020 04:40 p. m.
Para: r-help-es using r-project.org <r-help-es using r-project.org>
Asunto: [R-es] Pregunta sobre rLandsat



	[[alternative HTML version deleted]]