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 en r-project.org> en nombre de Francisco Rodr?guez <fjroar en hotmail.com> > Enviado: mi?rcoles, 19 de febrero de 2020 04:40 p. m. > Para: r-help-es en r-project.org <r-help-es en r-project.org> > Asunto: [R-es] Pregunta sobre rLandsat > > > > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > R-help-es en 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
Hola ?ngel con todos mis respetos, pero la pregunta que hab?a planteado yo es la
siguiente (donde en principio espero a ver si alguien me responde, entiendo que
su problema es otro relacionado con rLandsat que convendr?a que tuviese otro
motivo distinto, no s? si la lista de correos la ha liado o ha habido una
confunsi?n, pero mejor separar los 2 temas):
Hola grupo, estoy siguiendo una gu?a de la librer?a rLandsat que me la he
descargado de:
devtools::install_github("socialcopsdev/rLandsat")
Y tras hacer los siguiente (obviamente tengo me he registrado previamente en la
api correspondiente):
product_id = c("LC08_L1TP_145049_20180301_20180308_01_T1",
"LC08_L1TP_145049_20170330_20170414_01_T1",
"LC08_L1TP_145049_20140407_20170424_01_T1")
result_order = espa_order(product_id,
product = "rs")
Me da el siguiente error:
Invalid product selection. Please check the 'product_available'
dataframe
He probado con todas las opciones para product
No s? si me he equivocado en algo, hay un bug o hay alg?n sitio m?s recomendable
para atacar, desde R, estas im?genes
Un saludo y cualquier ayuda es bienvenida, gracias de antemano
[[alternative HTML version deleted]]
Francisco Rodriguez Sanchez
2020-Feb-20 11:19 UTC
[R-es] Crear raster a partir de tabla de datos (ERA: Pregunta sobre rLandsat)
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 en r-project.org> en nombre de
Francisco Rodr?guez <fjroar en hotmail.com>
>> Enviado: mi?rcoles, 19 de febrero de 2020 04:40 p. m.
>> Para: r-help-es en r-project.org <r-help-es en r-project.org>
>> Asunto: [R-es] Pregunta sobre rLandsat
>>
>>
>>
>> [[alternative HTML version deleted]]
>>
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es en 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]]
Hola Francisco,
Es posible que haya una errata en el código?
result_order = espa_order(product_id,
product = "rs")
(product = "rs" en lugar de product = "sr")
Si este paquete no te convence, hay otros para descargar landsat, como
RGISTOOLS:
https://www.rdocumentation.org/packages/RGISTools/versions/1.0.0/topics/lsDownload
Suerte
Paco
El 20/02/2020 a las 11:15, Francisco Rodríguez escribió:> Hola ?ngel con todos mis respetos, pero la pregunta que hab?a planteado yo
es la siguiente (donde en principio espero a ver si alguien me responde,
entiendo que su problema es otro relacionado con rLandsat que convendr?a que
tuviese otro motivo distinto, no s? si la lista de correos la ha liado o ha
habido una confunsi?n, pero mejor separar los 2 temas):
>
> Hola grupo, estoy siguiendo una gu?a de la librer?a rLandsat que me la he
descargado de:
>
> devtools::install_github("socialcopsdev/rLandsat")
>
> Y tras hacer los siguiente (obviamente tengo me he registrado previamente
en la api correspondiente):
>
> product_id = c("LC08_L1TP_145049_20180301_20180308_01_T1",
> "LC08_L1TP_145049_20170330_20170414_01_T1",
> "LC08_L1TP_145049_20140407_20170424_01_T1")
>
> result_order = espa_order(product_id,
> product = "rs")
>
> Me da el siguiente error:
>
> Invalid product selection. Please check the 'product_available'
dataframe
>
> He probado con todas las opciones para product
>
> No s? si me he equivocado en algo, hay un bug o hay alg?n sitio m?s
recomendable para atacar, desde R, estas im?genes
>
> Un saludo y cualquier ayuda es bienvenida, gracias de antemano
>
>
> [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en 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]]