Cómo hacer álgebra de mapas con la librería raster de lenguaje R

En artículos anteriores contemplamos la manera de hacer álgebra de mapas con lenguaje R usando rgdal o sin ella. Con el concurso de la librería raster también es posible. Para instalarla desde el intérprete, como superusuario, se hace uso de install.packages(), se escoge el mirror más cercano para finalmente seleccionar raster en el menú de opciones.

Una vez instalada, nos movemos con setwd al espacio de trabajo que contiene los rásteres, los cargamos y operamos con ellos. En el ejemplo que doy a continuación se combinan los procedimientos de rgdal y raster sólo con el fin de mostrar como pueden ser interconvertidos.

library(rgdal)
library(raster)
setwd('proyectosR')
b3<-readGDAL('b3_c.tif') #lectura con rgdal
b3_r<-raster(b3) #conversion a formato de lib raster
b4_r<-raster('b4_c.tif') #lectura directa con raster
ndvi<-(b4_r-b3_r)/(b4_r+b3_r) #calculo del ndvi
writeRaster(ndvi,'ndvi.tif',drivername='GTiff') #escritura directa con raster

El ráster obtenido, visualizado con QGIS, es el siguiente:

ndvi_R

Esta entrada fue publicada en Lenguaje R, SIG, Software Libre. Guarda el enlace permanente.

2 respuestas a Cómo hacer álgebra de mapas con la librería raster de lenguaje R

  1. baltazar dijo:

    hola jose muchas gracias por tu aportes.
    podrias subir por favor mas informacion sobre este paquete de R llamado raster.
    saludos.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. 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 )

Google+ photo

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

Conectando a %s