Obtener Modelo Digital de Elevación de la cuenca del río Biabo

Antiguamente el acceso a datos topográficos o de relieve en una zona de estudio era escasos o inexistentes. Actualmente las nuevas tecnologías están permitiendo que tengamos información detallada de un terreno. Dependiendo del nivel de detalle de los datos, podemos encontrar diferentes tipos de información en formato raster que asocia la altitud de un terreno con la información de pixeles de una imagen; a este tipo de información se le conoce como modelos digitales de elevación (MDE) o en sus siglas en inglés Digital Model Elevation (DEM).

Cuando se observa un modelo digital de elevación (DEM) en un mapa, no se ve una matriz de celdas. En lugar de eso, se ve una capa simbolizada con una rampa de color. El rango de colores nos muestra la variación que existe en terreno desde sus partes más altas a las más bajas.

Estos (DEM) según su resolución pueden provenir de diferentes fuentes como son la conversión de archivos vectoriales de puntos o curvas de nivel, captura de altitud de datos usando drones, o recopilación de información satelital. De esta última fuente Google Earth Engine nos permite obtener información de diferentes fuentes, entre las más conocidas podemos encontrar: ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer), STRM (stands for Shuttle Radar Topography Mission), GTOPO30, etc.

Para el estudio realizado se usó la información ASTER DEM la cual proveniente de la colección Terra ASTER, se centra en la corrección de las anomalías de la morfología del territorio y la presencia de masas de agua. Los cuales tienen un ancho de pixel de 30 metros de resolución. 

En el siguiente código fuente usando en Google Earth Engine y el área de la cuenca del río Biabo, se pudo obtener la información de los niveles de altitud de la cuenca y los valores de las pendientes.

//Codigo fuente en GGE para extrar valores de altitud y pendientes usando Aster DEM

//Autor: Mg.Sc. Abel Carmona Arteaga

// definir geometria de trabajo (crear una región o poligono que encierre nuestra cuenca en formato shapefile)

var area_trabajo = ee.FeatureCollection("users/abeljamiro/Biabo");

//el valor de pixel es 30m

var ALOSDSMDEM = ee.Image('NASA/NASADEM_HGT/001');

var elevation = ALOSDSMDEM.select('elevation').clip(area_trabajo);

var pendiente = ee.Terrain.slope(elevation).clip(area_trabajo);

var corte1 = elevation.clip(area_trabajo)

var corte2 = pendiente.clip(area_trabajo)

Map.addLayer( pendiente, {min: 0, max: 100}, 'pendiente');

Map.addLayer( elevation, {max:4000,min:0,palette: ['ce7e45','df923d','f1b555','fcd163','207401']}, 'elevation');

//Exportación de DEM y pendientes de ALOS DSM

Export.image.toDrive({

 image: corte1.select("elevation"),

 description: 'Altitud_ALOS_DSM_30m',

 scale: 30,

 region: area_trabajo})

   Export.image.toDrive({

 image: corte2.select("pendiente"),

 description: 'Pendiente_ALOS_DSM_30m',

 scale: 30,

 region: area_trabajo})

   Map.addLayer( area_trabajo,{palette: 'red'} ,'Cuenca')