Más

¿Obtiene valores de pendiente incorrectos de un ASTER DEM en QGIS?

¿Obtiene valores de pendiente incorrectos de un ASTER DEM en QGIS?


Tengo una capa DEM del satélite ASTER y quiero obtener la información de la pendiente de ella.

Utilizo análisis de terreno> pendiente, pero la salida es muy incorrecta. ¿Porqué es eso?

Aquí está el DEM:

Y aquí está la "Pendiente" que está totalmente equivocada. Debería estirar los valores de 0 a 89,9097, pero no es así.


Los DEM de Aster están en EPSG 4326 (lat lon WGS84). Según la documentación GDAL DEM:

Para proyecciones LatLong cerca del ecuador, donde las unidades de latitud y las unidades de longitud son similares, las unidades de elevación (z) se pueden convertir para que sean compatibles usando scale = 370400 (si la elevación está en pies) o scale = 111120 (si la elevación está metros). Para ubicaciones que no estén cerca del ecuador, sería mejor reproyectar su cuadrícula usandogdalwarpantes de usargdaldem.

En esta respuesta para un problema con sombreado (la escala y el factor Z no tienen efecto en el análisis de sombreado en QGIS), el usuario @Luke proporciona una forma de calcular la escala correcta ("relación de unidades verticales a horizontales") multiplicando el coseno de la latitud multiplicado por 111320. La escala afecta los cálculos de pendiente.

En la imagen de abajo pude producir pendientes en grados de un DEM en EPSG 4326 en latitudes que oscilan entre 19 ° 45N y 20 ° N. La escala (-s) se estableció en104550.67.

Si quieres trabajar con un DEM en un CRS que defina las unidades de medida de las coordenadas y las elevaciones en metros debes reproyectar el DEM usando raster / proyecciones / warp (reproyectar). Por ejemplo, reproyección de EPSG 4326 (lat lon wgs84) a EPSG 32164 (UTM WGS84 14N). Vea abajo:

En este caso, al calcular la pendiente, la "relación de unidades verticales a horizontales" (escala) se establecería en 1.


Probablemente porque la proyección no estaba en metros.

Pendiente con DEM proyectada en metros (UTM / 31N: EPSG 25831):

Pendiente con el mismo DEM proyectado en grados (WGS 84: EPSG 4326) yz = 1:

Pendiente con el mismo DEM proyectado en grados (WGS 84: EPSG 4326) yz = 111320 (en el ecuador una unidad horizontal, 1 °, correspondería aproximadamente a 111320 m):


Ver el vídeo: Ráster de pendiente a partir de un Modelos Digital de Elevación DEM con QGIS