Rásteres de longitud, latitud mediante la librería raster de lenguaje R

En un post anterior se refirió la posibilidad de producir rásteres de longitud, latitud empleando la librería rgdal. No obstante, un método más expedito es posible empleando el método ‘init’ de la librería raster. El código es el siguiente:

library(raster)
setwd('proyectosR')
r <- raster('tiznados_canoa.tif')
coor<-coordinates(r)
x<-function(x){coor[,1]}
y<-function(x){coor[,2]}
lon <- init(r, x, filename='lon.tif', overwrite=TRUE)
lat <- init(r, y, filename='lat.tif', overwrite=TRUE)

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(lon)
image(lat)

long_lat2

y en QGIS de esta otra (pseudocolor con 5 clases):

qgis

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

Una respuesta a Rásteres de longitud, latitud mediante la librería raster de lenguaje R

  1. Pingback: Rásteres de longitud, latitud con GDAL/PyQGIS | 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