[R-es] Crear raster a partir de tabla de datos (ERA: Pregunta sobre rLandsat)

[Francisco Rodriguez Sanchez]

Hola Ángel,

Prueba a pasarle el argumento `resolution = 30` cuando creas el raster. 
También puedes especificar el CRS al crear el raster (en ese caso no 
habría que reproyectar los datos a longlat previamente):

r1 <- raster(nrows=1386, ncols=1649,
               xmn=xmn, xmx=xmx,
               ymn=ymn, ymx=ymx,
resolution = 30 ,
crs = "+proj=utm +zone=14 +datum=WGS84")


Suerte

Paco


El 20/02/2020 a las 11:01, Marcelino De La Cruz Rot escribió:
> 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
>
-- 
Dr Francisco Rodriguez-Sanchez
https://www.linkedin.com/in/frodriguezsanchez/


	[[alternative HTML version deleted]]