Suma de archivos raster en R
El objetivo de este ejercicio es realizar la suma de 12 archivos raster de manera óptima en R, estos archivos son la precipitación para Colombia en los 12 meses del año, la idea es que el archivo final me represente la precipitación acumulada de este país.
Hacemos uso de las librerias «raster» y «rgdal»
require(raster)
require(rgdal)
Hacemos un listado de los archivos raster (estos deben de estar sin más archivos en una sola carpeta).
listado = list.files(«E:/Blogger/Post_3/_datos/_prec_mensual/», pattern = «asc», full.names = T) #usamos el pattern para que solo escoga los archivos tif y fullnames para que tome toda la dirección de donde se encuentran los archivos.
Luego le decimos a R que esos archivos son de tipo Raster
datos = lapply(listado, FUN=raster)
Hacemos un stack, es decir poner todos los archivos unidos.
datos_raster_stack = stack(datos)
Ahora si hacemos la suma entre los archivos.
suma_raster = sum(datos_raster_stack)
Y por último, escribimos el archivo resultante donde queramos.
writeRaster(suma_raster, «E:/Blogger/Post_3/_datos/_prec_acum/sum_prec.tif»)
plot(suma_raster)
plot(suma_raster)
De este link se pueden descargar los datos del presente ejercicio.
SIG con R - Tutoriales 
Que tan sencillo parece….Super Fabio. Muy pedagógico.
Me gustaMe gusta
Muy claro, muy bien explicado Fabio
Me gustaMe gusta