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

[Jorge Tornero]

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



-- 
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <http://lists.osgeo.org/pipermail/qgis-es/attachments/20181113/a49d3e0b/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/20181113/a49d3e0b/attachment-0001.jpg>