R: generación modelo digital de elevaciones de RD uniendo dos imágenes separadas

Script de R para generar un modelo digital de elevaciones de RD uniendo dos imágenes separadas en EPSG:4326 (lat/lon WGS84). El ámbito del país se obtiene del shapefile “provincias” de ONE, y el DEM proviene del SRTM, descargado con la función getData() del paquete raster de R, que a la vez se obtiene desde servidor FTP del Consorcio para la Información Espacial CGIAR.

#CARGA DE PAQUETES
library(raster)
library(rgdal)
 
#FIJACIÓN DIRECTORIO
setwd(tempdir()) #FIJAR EL DIRECTORIO DE TRABAJO
getwd() #MUESTRA DÓNDE ESTÁ EL DIRECTORIO DE TRABAJO
 
#ÁMBITO RD A PARTIR DE PROVINCIAS DE ONE
temp <- tempfile() #ASIGNAR EL NOMBRE DE UN ARCHIVO TEMPORAL AL OBJETO temp
download.file("http://geografiafisica.org/r/divisionrd/provincias.zip",temp) #DESCARGA EL ARCHIVO ZIP Y LO NOMBRA COMO TEMPORAL
unzip(temp) #DESCOMPRIME EN EL DIRECTORIO DE TRABAJO
unlink(temp) #DESVINCULA DEL ARCHIVO TEMPORAL Y LO BORRA
provincias <- shapefile('PROVCenso2010.shp') #LECTURA DE LAS PROVINCIAS DESCARGADAS
provincias <- spTransform(provincias,CRS('+init=epsg:4326')) #REPROYECCIÓN A LAT/LON WGS84
ambitord <- extent(provincias) #GENERACIÓN DE UN OBJETO DE EXTENSIÓN QUE ABARCA EL ÁMBITO DEL PAÍS (ENVELOPE) DEL OBJETO provincias
 
#DESCARGAR DEM
getData('SRTM',lat=19,lon=-70) #DESCARGA DEM ORIENTAL AL DIRECTORIO DE TRABAJO 
getData('SRTM',lat=19,lon=-71) #DESCARGA DEM OCCIDENTAL AL DIRECTORIO DE TRABAJO 
demO <- raster(readGDAL("srtm_23_09.tif")) #LECTURA DESDE DIRECTORIO Y CONVERSIÓN A RASTER
demE <- raster(readGDAL("srtm_22_09.tif")) #LECTURA DESDE DIRECTORIO Y CONVERSIÓN A RASTER
dem <- raster::merge(demE,demO) #UNIÓN
plot(dem) #MUESTRA EL DEM. NOTAR QUE ESTE ABARCA UNA PARTE MARINA AL SUDOESTE DEL PAÍS
dem <- crop(dem,ambitord) #RECORTANDO AL ÁMBITO (ENVELOPE) DE LAS PROVINCIAS
plot(dem) #MUESTRA EL DEM RECORTADO

Created by Pretty R at inside-R.org

Dr. José Ramón Martínez Batlle (Ph.D)

Un pensamiento en “R: generación modelo digital de elevaciones de RD uniendo dos imágenes separadas

  1. Pingback: Estadísticos descriptivos de la altura por municipios de República Dominicana | Geografía Física – República Dominicana – Dr. José Ramón Martínez Batlle

Deja un comentario

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *