Suma de archivos Raster en R

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)





De este link se pueden descargar los datos del presente ejercicio.

2 comentarios sobre “Suma de archivos Raster en R

Responder

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión /  Cambiar )

Google photo

Estás comentando usando tu cuenta de Google. Cerrar sesión /  Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión /  Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión /  Cambiar )

Conectando a %s