Recortar archivos raster con extensión en metros y grados.
Necesito cortar un raster basado en otro más pequeño para realizar algunas operaciones, sin embargo, tengo un problema y es que la extensión está expresada de manera diferente en ambos (metros y grados), por lo que no puedo proceder, hay una forma de resolverlo.
El raster en el que quiero recortar lo obtuve apilando varias capas en un solo rasterstack mediante muestreo, también he intentado usar solo la capa de referencia de muestreo, pero obtengo el mismo error.
La capa que voy a cortar es el producto de una predicción RK.
Mi código a continuación:
Mgpred <- stack(raster("REALIZACION1_mapasrkRF.tif"))
Covs <- stack(read.RDS("covs1km")
extent(Mgpred)
#class : Extent
#xmin : -14388.19
#xmax : 1009756
#ymin : 1113870
#ymax : 2126756
extent(Covs)
#class : Extent
#xmin : -91.04483
#xmax : -90.89741
#ymin : 14.72128
#ymax : 14.86761
e <- extent(Covs)
cr <- crop(Mgpred, e)
#Error in .local(x, y, …) : las extensiones no se superponen.
davy.ai
Hay algunas opciones para resolver el problema de las diferentes expresiones de extensión en los dos raster. Una opción es convertir un raster a la misma expresión de extensión que el otro. En este caso, parece que el raster basado en la predicción RK se expresa en metros mientras que el otro raster se expresa en grados. Por lo tanto, podríamos convertir la extensión del raster de predicción RK a grados usando la función
projectRaster
del paqueteraster
:Ahora que ambos raster están en la misma expresión de extensión, podemos cortarlos usando la función
crop
.Otra opción es usar la función
resample
para hacer que ambos raster tengan la misma extensión y resolución, y luego cortarlos:Cualquiera de estas opciones debería permitirle cortar correctamente el raster basado en la predicción RK a la extensión del raster más pequeño, para que pueda realizar las operaciones deseadas.