[QGIS-es] Generación de RRIM (Red Relief Image Maps)

Sab Feb 23 04:17:43 PST 2019

Hola Jorge,

Muchas gracias por tu aportación a la construcción de del RRIM con
herramientas libres.

Me resultó muy interesante cuando lo planteaste en el grupo de Telegram
peor no había tenido tiempo de mirarlo hasta estos días en los que estaba
mostrándole el potencial de QGIS a un amigo mexicano.

Te comento que en realidad no hace falta salirse de QGIS para lanzar las
herramientas de SAGA: Desde Processing Toolbox puedes hacer los cálculos de
los índices de apertura topográfica corriendo los algoritmos de SAGA.

De este modo también es posible crear un modelo (Processing Modeler) que
nos permita automatizar las salidas para esta interesante representación
topográfica. Y es esto lo que hice y os adjunto.

Aunque el resultado considero que es aún muy mejorable, para que podáis
recrear la imagen que os muestro, de los volcanes de la sierra del
Popocatepetl en las inmediaciones de Ciudad de México, os adjunto
igualmente los archivos de estilo de las dos capas.

Espero que sea de utilidad.

Un saludo

--
*Miguel Sevilla-Callejo*
Doctor en Geografía

PD. Los datos del modelo digital de elevaciones son los del SRTM de la NASA
(existe un plug-in que te permite automatizar su descarga).

On Tue, 13 Nov 2018 at 13:23, Jorge Tornero <jorge.tornero en ieo.es> wrote:

> Buenos días:
>
> Como ya comenté en el grupo de Telegram, una compañera de trabajo me
> preguntó cómo hacer un Red Relief Image Map, descrito en
> https://www.rrim.jp/en/ . Estos mapas son una forma novedosa de
> representar el relieve, y básicamente, hasta lo que alcanzo a comprender,
> se trata de combinar información de pendientes de un MDT con la información
> de convexidad y concavidad del terreno (lo que los autores de este método
> llaman Positive Openess y Negative Openess), asignando el color mediante la
> pendiente y el brillo mediante la combinación de concavidad y convexidad. A
> esta útima combinación la denominan Ridge-Valley Index y es el resultado de
> (PositiveOpeness - NegativeOpeness)/2
>
> Los principios del método son explicados aquí:
> http://www.isprs.org/proceedings/XXXVII/congress/2_pdf/11_ThS-6/08.pdf.
> El RRIM es un procedimiento patentado (en Europa no) aunque existe un
> permiso para su utilización académica, etc.
>
> Su implementación en QGIS puro, hasta lo que yo sé, no es posible, ya que
> QGIS carece de herramientas propias para calcular la Openess. Sin embargo,
> SAGA sí que las posee por lo que mi flujo de trabajo hace uso de SAGA.
> Desafortunadamente, en mi instalación de QGIS no cuento con la integración
> de SAGA, por lo que lo he utilizado directamente, exportando la información
> necesaria a QGIS.
>
> Requerimientos:
>
> 1) QGIS y SAGA instalados en tu sistema
>
> 2) Un MDT para efectuar el proceso (descargado del IGN, por ejemplo)
>
> Flujo de trabajo:
>
> 1) Importamos el MDT en SAGA. En  mi caso, dado que las descargas del IGN
> son en formato .asc:
>
> Geoprocessing -> File -> Grid -> Import -> ESRI Arc/Info Grid
>
> 2) Generamos el mapa de pendientes a partir del MDT:
>
> Geoprocessing -> Terrain Analysis -> Morphometry > Slope, Aspect, Curvature
>
> Como Grid y Elevation de entrada seleccionamos las correspondientes al MDT
> que hemos cargado. Y en cuanto a datos de salida, sólo sería necesario
> generar Slope
>
> 3) Generamos los mapas de Positive/negative openess:
>
> Geoprocessing -> Terrain Analysis -> Ligthing -> Topographic Openess
>
> De nuevo, como Grid y Elevation de entrada seleccionamos las
> correspondientes al MDT que hemos cargado, y generaremos tanto la positive
> como la negative
>
> Los parámetros de generación de openess por defecto funcionan bien, aunque
> sería posible modificarlos y observar el efecto en el mapa generado.
>
> Por cierto, el proceso de generar la Topographic Openess puede durar un
> buen rato.
>
> 4) Exportamos los resultados para poder cargarlos en QGIS:
>
> Geoprocessing -> File -> GDAL/OGR -> GDAL: Export raster to GeoTiff
>
> Seleccionamos la grid correspondiente y exportamos cada imagen por
> separado.
>
> 5) Cargamos los ráster generados con SAGA en QGIS: Pendiente, Openess
> Positiva y Openess Negativa.
>
> 6) Creamos la capa RVI (Ridge Valley Index) mediante la calculadora
> ráster. Recordad que RVI=(OpenessPos-OpenessNeg)/2
>
> 7) Cargamos la capa RVI sobre la de Pendientes y le aplicamos como modo de
> mezcla Multiplicar, estableciendo los valores de máximo y mínimo de capa de
> manera adecuada (aquí yo he tenido que ir probando, aunque lo más adecuado
> ha sido, casi siempre, Estirar a min-max y cargar valores min-max)
>
> 8) Coloreamos de rojo la capa de pendientes, estableciendo el gradiente de
> colores de la capa en "De blanco a negro"
>
> 9) Ya tenemos el RRIM generado!!!
>
>
> Advertencias:
>
> Este flujo de trabajo ha sido hecho para intentar reproducir el RRIM por
> alguien con bastantes lagunas en el manejo avanzado de GIS, por lo que
> puede ser considerado como un punto de partida más que otra cosa. Así, por
> ejemplo, con respecto al manejo de los valores mínimos y máximos de las
> capas, he tenido problemas para trabajar con una imagen del estrecho, ya
> que incluye un sector grande con NaNs
>
>  Aunque el modo de mezcla que más me resulta apropiado es el de
> Multiplicar, se pueden obtener diferentes resultados con otros modos de
> mezcla. Además, es posible que lo que especifica el método no se consiga
> estrictamente con modos de mezcla, sino por ejemplo creando una imagen
> multicanal en la que el canal rojo sea la pendiente y el canal de brillo el
> RVI. Crear semejante imagen en QGIS por ahora se me escapa y puede que esté
> diciendo incluso alguna barbaridad.
>
> Espero que os pueda servir de inspiración y, si es de vuestro interés, se
> pueda ir mejorando
>
>
> Un saludo
>
>
> Jorge
>
>
>
> --
> _______________________________________________
> QGIS-es mailing list
> QGIS-es en lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/qgis-es
>
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <http://lists.osgeo.org/pipermail/qgis-es/attachments/20190223/50754a92/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: Firma-Correo-Jorge-Tornero.jpg
Type: image/jpeg
Size: 83670 bytes
Desc: no disponible
URL: <http://lists.osgeo.org/pipermail/qgis-es/attachments/20190223/50754a92/attachment-0002.jpg>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: rrim_popocatepetl_20190223.jpg
Type: image/jpeg
Size: 143891 bytes
Desc: no disponible
URL: <http://lists.osgeo.org/pipermail/qgis-es/attachments/20190223/50754a92/attachment-0003.jpg>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: rrim1.model3
Type: application/octet-stream
Size: 9832 bytes
Desc: no disponible
URL: <http://lists.osgeo.org/pipermail/qgis-es/attachments/20190223/50754a92/attachment-0003.obj>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: rrim_rvi.qml
Type: application/octet-stream
Size: 1521 bytes
Desc: no disponible
URL: <http://lists.osgeo.org/pipermail/qgis-es/attachments/20190223/50754a92/attachment-0004.obj>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: rrim_slope.qml
Type: application/octet-stream
Size: 2690 bytes
Desc: no disponible
URL: <http://lists.osgeo.org/pipermail/qgis-es/attachments/20190223/50754a92/attachment-0005.obj>