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]]
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]]
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