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 <[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
--
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
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
