Álgebra de mapas con lenguaje R usando la librería rgdal

En el artículo anterior se consideró el álgebra de mapas con lenguaje R sin el concurso de rgdal. A veces es necesario intentar otras vías; sobre todo cuando no se encuentra la manera de instalar una librería como rgdal la cual no viene incluida en la oferta del r-base. Sin embargo, en Debian y Debian derivadas como Ubuntu he encontrado un procedimiento para instalar sin problemas rgdal; algo que todavía no he probado en ambiente Windows.

Como ya tengo instalada la librería rgdal vamos a probar que tantas facilidades nos permite a la hora de realizar el álgebra de mapas en el intérprete de R. Para ello se va a considerar el mismo ejemplo de la determinación del NDVI del artículo anterior. Por tanto, lanzamos R, se carga la librería rgdal y luego nos cambiamos al directorio donde se tienen los mapas para cargar las bandas objeto de cálculo:

setwd('proyectoR')
> library(rgdal)
Loading required package: sp
rgdal: version: 0.8-16, (SVN revision 498)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.9.0, released 2011/12/29
Path to GDAL shared files: /usr/share/gdal/1.9
Loaded PROJ.4 runtime: Rel. 4.7.1, 23 September 2009, [PJ_VERSION: 470]
Path to PROJ.4 shared files: (autodetected)
>b3<-readGDAL('b3_c.tif') #lectura directa del formato GTiff
b3_c.tif has GDAL driver GTiff 
and has 256 rows and 280 columns
>b4<-readGDAL('b4_c.tif') #lectura directa del formato GTiff
b4_c.tif has GDAL driver GTiff 
and has 256 rows and 280 columns
> names(b3); names(b4)
[1] "band1"
[1] "band1"
>NDVI<-(b4$band1-b3$band1)/(b4$band1+b3$band1) #cálculo del NDVI
>ndvi<-b3 #clonado del objeto b3 en ndvi con formato SpatialGridDataFrame
>ndvi$band1<-NDVI #copiado de los atributos de NVDI a ndvi
writeGDAL(ndvi,"ndvi.tif",drivername="GTiff")

El mapa resultante (*.tif) se obtiene en tiempo bastante breve, con su georreferenciación incorporada y no hay que usar a QGIS como intermediario para ningún cambio de formato. La visualización en QGIS es la siguiente:

rgdal

A pesar de lo complicado que pueda lucir este procedimiento se puede automatizar en lenguaje de script y luce más compacto que Python.

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

3 respuestas a Álgebra de mapas con lenguaje R usando la librería rgdal

  1. Pingback: Cómo crear un SpatialGridDataFrame en lenguaje R | El Blog de José Guerrero

  2. Pingback: Álgebra de mapas con lenguaje R usando l...

  3. Pingback: Cómo hacer álgebra de mapas con la librería raster de lenguaje R | El Blog de José Guerrero

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