<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>Buenos días:</p>
<p>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
<a class="moz-txt-link-freetext" href="https://www.rrim.jp/en/">https://www.rrim.jp/en/</a> . 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<br>
</p>
<p>Los principios del método son explicados aquí:
<a class="moz-txt-link-freetext" href="http://www.isprs.org/proceedings/XXXVII/congress/2_pdf/11_ThS-6/08.pdf">http://www.isprs.org/proceedings/XXXVII/congress/2_pdf/11_ThS-6/08.pdf</a>.
El RRIM es un procedimiento patentado (en Europa no) aunque existe
un permiso para su utilización académica, etc.<br>
</p>
<p>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.</p>
<p>Requerimientos:</p>
<p>1) QGIS y SAGA instalados en tu sistema</p>
<p>2) Un MDT para efectuar el proceso (descargado del IGN, por
ejemplo)</p>
<p>Flujo de trabajo:</p>
<p>1) Importamos el MDT en SAGA. En mi caso, dado que las descargas
del IGN son en formato .asc:</p>
<p>Geoprocessing -> File -> Grid -> Import -> ESRI
Arc/Info Grid<br>
</p>
<p>2) Generamos el mapa de pendientes a partir del MDT:</p>
<p>Geoprocessing -> Terrain Analysis -> Morphometry >
Slope, Aspect, Curvature</p>
<p>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</p>
<p>3) Generamos los mapas de Positive/negative openess:</p>
<p>Geoprocessing -> Terrain Analysis -> Ligthing ->
Topographic Openess</p>
<p>De nuevo, como Grid y Elevation de entrada seleccionamos las
correspondientes al MDT que hemos cargado, y generaremos tanto la
positive como la negative<br>
</p>
<p>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.</p>
<p>Por cierto, el proceso de generar la Topographic Openess puede
durar un buen rato.<br>
</p>
<p>4) Exportamos los resultados para poder cargarlos en QGIS:</p>
<p>Geoprocessing -> File -> GDAL/OGR -> GDAL: Export raster
to GeoTiff<br>
</p>
<p>Seleccionamos la grid correspondiente y exportamos cada imagen
por separado.<br>
</p>
<p>5) Cargamos los ráster generados con SAGA en QGIS: Pendiente,
Openess Positiva y Openess Negativa.<br>
</p>
<p>6) Creamos la capa RVI (Ridge Valley Index) mediante la
calculadora ráster. Recordad que RVI=(OpenessPos-OpenessNeg)/2</p>
<p>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)</p>
<p>8) Coloreamos de rojo la capa de pendientes, estableciendo el
gradiente de colores de la capa en "De blanco a negro"<br>
</p>
<p>9) Ya tenemos el RRIM generado!!!</p>
<p><br>
</p>
<p>Advertencias:</p>
<p>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</p>
<p> 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.</p>
<p>Espero que os pueda servir de inspiración y, si es de vuestro
interés, se pueda ir mejorando</p>
<p><br>
</p>
<p>Un saludo</p>
<p><br>
</p>
<p>Jorge<br>
</p>
<p><br>
</p>
<p><br>
</p>
<div class="moz-signature">-- <br>
<img src="cid:part1.E711C120.2F2F1DCC@ieo.es" border="0"></div>
</body>
</html>