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)
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