[R-es] Pregunta sobre rLandsat

[Marcelino De La Cruz Rot]

Hola Ángel:


Yo creo que tendrías que establecer el sistema de coordenadas de 
referencia de tu objeto raster antes de salvarlo como GTiff. Algo así:

crs(r1) <-"+proj=utm +zone=14 +datum=WGS84"

Saludos,

Marcelino.



El 20/02/2020 a las 1:42, Angel Cervantes escribió:
> 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]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es


-- 
Marcelino de la Cruz Rot
Depto. de Biología y Geología
Física y Química Inorgánica
Universidad Rey Juan Carlos
Móstoles España