Cómo producir rásteres de longitud, latitud mediante lenguaje R

Un usuario de un foro sobre temas SIG solicita un procedimiento para producir rásteres de longitud y latitud para ser empleados en una fórmula de regresión que genera valores de precipitación. Se señala que en ArcGis no es posible hacerlo de una manera inmediata; aunque Idrisi si cuenta con un módulo para tal efecto. Como mi interés es señalar que este tipo de tareas es factible de realizar mediante el empleo de software libre, inmediatamente visualicé que ello era posible, sin programar, mediante el uso de lenguaje R.

Para ello:

setwd('proyectosR') #establece directorio de trabajo
library(rgdal)
tiz<-readGDAL('tiznados_canoa.tif') #carga raster de interes
coor<-coordinates(tiz) #extrae coordenadas x,y del raster
long<-coor[,1] #extrae coordenada x
lat<-coor[,2] #extrae coordenada y
tiz_long<-tiz #raster de longitud
tiz_long$band1<-long #sustitucion de values por longitud
tiz_lat<-tiz #raster de latitud
tiz_lat$band1<-lat #sustitucion de values por latitud
writeGDAL(tiz_long,'tiz_long.tif',drivername='GTiff') #raster de longitud
writeGDAL(tiz_lat,'tiz_lat.tif',drivername='GTiff') #raster de latitud

Los rásteres de longitud y latitud (formato SpatialGridDataFrame) visualizados en el intérprete de R lucen de esta manera:

par(mfcol=c(1,2))
image(tiz_long)
image(tiz_lat)

long_lat2

Ambos rásteres en formato GTiff, conjuntamente con el DEM, se cargaron en QGIS. El plugin Value Tool de QGIS permitió corroborar que los valores del ráster de longitud solo aumentaban en sentido horizontal mientras que los de latitud solo disminuían en sentido vertical; ambos según las dimensiones de celda del ráster (es decir, tal como era de esperar). A continuación, una imagen que resume el proceso:

long_lat

donde los valores que señala Value Tool son los que corresponden a la posición del cursor en la imagen.

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

Una respuesta a Cómo producir rásteres de longitud, latitud mediante lenguaje R

  1. Pingback: Rásteres de longitud, latitud mediante 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