El proyecto projectRaster no logra cambiar el crs cuando se aplica a un objeto de tipo lista en R.
Quiero apilar 6 rásteres en una lista llamada allrasters, pero primero debo corregir inconsistencias en crs y extensión. Aquí está mi intento de código para establecer el segundo ráster de la lista con el crs del tercer ráster de la lista:
“`{r}
projectRaster(allrasters[[2]], crs=crs(allrasters[[3]]))
<pre><code><br />Sin embargo, cuando ejecuto este código y verifico, allrasters[[2]] sigue siendo proj.merc y nada ha cambiado…
Información del ráster:
</code></pre>
crs(allrasters[[2]])
Argumentos CRS:
+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0
+x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext
+no_defs
crs(allrasters[[3]])
Argumentos CRS:
+proj=aea +lat_0=0 +lon_0=-120 +lat_1=34 +lat_2=40.5
+x_0=0 +y_0=-4000000 +datum=NAD83 +units=m +no_defs
“`
davy.ai
Has intentado utilizar la función
projectRaster
para establecer el CRS del segundo ráster enallrasters
al CRS del tercer ráster. Sin embargo, parece que esta función no ha cambiado el CRS del ráster como se esperaba.Antes de intentar apilar los rásteres, es importante asegurarse de que tengan la misma extensión y CRS. Para establecer el CRS del segundo ráster para que coincida con el CRS del tercer ráster, puedes utilizar la función
spTransform
del paquetergdal
. Aquí tienes un ejemplo de código:“`{r}
library(rgdal)
establecer el CRS del segundo ráster para que coincida con el CRS del tercer ráster
allrasters[[2]] <- spTransform(allrasters[[2]], CRSobj = crs(allrasters[[3]]))
comprobar que el CRS se ha actualizado
crs(allrasters[[2]])
“`
Este código debería actualizar el CRS del segundo ráster para que coincida con el del tercer ráster. Luego puedes repetir este proceso para los otros rásteres en la lista, según sea necesario, para asegurarte de que todos tengan el mismo CRS antes de intentar apilarlos.