<div dir="auto">Hola muy interesante,una duda:<div dir="auto"><br></div><div dir="auto">Cuando se dice <span style="font-family:sans-serif;font-size:14px">Cargamos la capa RVI sobre la de Pendientes y le aplicamos como modo de mezcla Multiplicar,</span><span style="font-family:sans-serif;font-size:14px"> </span></div><div dir="auto"><span style="font-family:sans-serif;font-size:14px"><br></span></div><div dir="auto"><span style="font-family:sans-serif;font-size:14px">A qué se refiere o como se hace el modo de mezcla multiplicar?</span></div><div dir="auto"><span style="font-family:sans-serif;font-size:14px"><br></span></div><div dir="auto"><font face="sans-serif"><span style="font-size:14px">Disculpen mi ignorancia.</span></font></div><div dir="auto"><font face="sans-serif"><span style="font-size:14px"><br></span></font></div><div dir="auto"><font face="sans-serif"><span style="font-size:14px">Saludos</span></font></div><div dir="auto"><font face="sans-serif"><span style="font-size:14px"><br></span></font></div><div dir="auto"><font face="sans-serif"><span style="font-size:14px">Carlos</span></font></div><div dir="auto"><font face="sans-serif"><span style="font-size:14px"><br></span></font></div><div dir="auto"><font face="sans-serif"><span style="font-size:14px"><br></span></font></div><div dir="auto"><br></div><div dir="auto"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El dom., 24 feb. 2019 3:18 a. m., JORGE TORNEO NUNEZ <<a href="mailto:jorge.tornero@ieo.es">jorge.tornero@ieo.es</a>> escribió:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"> Buenas tardes, Miguel:<br>
<br>
Me alegro que la pequeña aportación vaya siendo de utilidad (aunque la <br>
aportación de verdad la hicieron los inventores del RRIM). <br>
Efectivamente, es posible, com dices, hacerlo directamente a través de <br>
QGIS pero lamentablemente el QGIS del que disponía en ese momento y el <br>
sistema en conjunto no tenía disponibles los algoritmos de SAGA para <br>
hacerlo sin salirse de QGIS. La verdad es que las imágenes que se <br>
generan son casi siempre espectaculares, digo casi siempre porque <br>
probablemente dependan mucho de la resolución de los MDT de entrada. <br>
Nosotros hicimos unas pruebas con batimetrías del golfo de Cádiz y el <br>
resultado no era tan espectacular. Pero cuando los datos o el relieve <br>
se prestan a ello, salen cosas muy chulas.<br>
<br>
Un saludo<br>
<br>
Jorge Tornero<br>
<br>
Miguel Sevilla-Callejo <<a href="mailto:msevilla00@gmail.com" target="_blank" rel="noreferrer">msevilla00@gmail.com</a>> escribió:<br>
<br>
> Hola Jorge,<br>
> <br>
> Muchas gracias por tu aportación a la construcción de del RRIM <br>
> con herramientas libres.<br>
> <br>
> Me resultó muy interesante cuando lo planteaste en el grupo de <br>
> Telegram peor no había tenido tiempo de mirarlo hasta estos días en <br>
> los que estaba mostrándole el potencial de QGIS a un amigo mexicano.<br>
> <br>
> Te comento que en realidad no hace falta salirse de QGIS para <br>
> lanzar las herramientas de SAGA: Desde Processing Toolbox puedes <br>
> hacer los cálculos de los índices de apertura topográfica corriendo <br>
> los algoritmos de SAGA.<br>
> <br>
> De este modo también es posible crear un modelo (Processing <br>
> Modeler) que nos permita automatizar las salidas para esta <br>
> interesante representación topográfica. Y es esto lo que hice y os <br>
> adjunto.<br>
> <br>
> Aunque el resultado considero que es aún muy mejorable, para que <br>
> podáis recrear la imagen que os muestro, de los volcanes de la <br>
> sierra del Popocatepetl en las inmediaciones de Ciudad de México, os <br>
> adjunto igualmente los archivos de estilo de las dos capas.<br>
> <br>
> Espero que sea de utilidad.<br>
> <br>
> Un saludo<br>
> <br>
> --<br>
> MIGUEL SEVILLA-CALLEJO<br>
> Doctor en Geografía<br>
> <br>
> PD. Los datos del modelo digital de elevaciones son los <br>
> del SRTM de la NASA (existe un plug-in que te permite automatizar su <br>
> descarga).<br>
><br>
> On Tue, 13 Nov 2018 at 13:23, Jorge Tornero <br>
> <<a href="mailto:jorge.tornero@ieo.es" target="_blank" rel="noreferrer">jorge.tornero@ieo.es</a>> wrote:<br>
><br>
>> Buenos días:<br>
>><br>
>> Como ya comenté en el grupo de Telegram, una compañera de <br>
>> trabajo me preguntó cómo hacer un Red Relief Image Map, descrito en <br>
>> <a href="https://www.rrim.jp/en/" rel="noreferrer noreferrer" target="_blank">https://www.rrim.jp/en/</a> . Estos mapas son una forma novedosa de <br>
>> representar el relieve, y básicamente, hasta lo que alcanzo a <br>
>> comprender, se trata de combinar información de pendientes de un <br>
>> MDT con la información de convexidad y concavidad del terreno (lo <br>
>> que los autores de este método llaman Positive Openess y Negative <br>
>> Openess), asignando el color mediante la pendiente y el brillo <br>
>> mediante la combinación de concavidad y convexidad. A esta útima <br>
>> combinación la denominan Ridge-Valley Index y es el resultado de <br>
>> (PositiveOpeness - NegativeOpeness)/2<br>
>><br>
>> Los principios del método son explicados aquí: <br>
>> <a href="http://www.isprs.org/proceedings/XXXVII/congress/2_pdf/11_ThS-6/08.pdf" rel="noreferrer noreferrer" target="_blank">http://www.isprs.org/proceedings/XXXVII/congress/2_pdf/11_ThS-6/08.pdf</a>. El <br>
>> RRIM es un procedimiento patentado (en Europa no) aunque existe un <br>
>> permiso para su utilización académica, etc.<br>
>><br>
>> Su implementación en QGIS puro, hasta lo que yo sé, no es <br>
>> posible, ya que QGIS carece de herramientas propias para calcular <br>
>> la Openess. Sin embargo, SAGA sí que las posee por lo que mi flujo <br>
>> de trabajo hace uso de SAGA. Desafortunadamente, en mi instalación <br>
>> de QGIS no cuento con la integración de SAGA, por lo que lo he <br>
>> utilizado directamente, exportando la información necesaria a QGIS.<br>
>><br>
>> Requerimientos:<br>
>><br>
>> 1) QGIS y SAGA instalados en tu sistema<br>
>><br>
>> 2) Un MDT para efectuar el proceso (descargado del IGN, por ejemplo)<br>
>><br>
>> Flujo de trabajo:<br>
>><br>
>> 1) Importamos el MDT en SAGA. En mi caso, dado que las <br>
>> descargas del IGN son en formato .asc:<br>
>><br>
>> Geoprocessing -> File -> Grid -> Import -> ESRI Arc/Info Grid<br>
>><br>
>> 2) Generamos el mapa de pendientes a partir del MDT:<br>
>><br>
>> Geoprocessing -> Terrain Analysis -> Morphometry > Slope, <br>
>> Aspect, Curvature<br>
>><br>
>> Como Grid y Elevation de entrada seleccionamos las <br>
>> correspondientes al MDT que hemos cargado. Y en cuanto a datos de <br>
>> salida, sólo sería necesario generar Slope<br>
>><br>
>> 3) Generamos los mapas de Positive/negative openess:<br>
>><br>
>> Geoprocessing -> Terrain Analysis -> Ligthing -> Topographic Openess<br>
>><br>
>> De nuevo, como Grid y Elevation de entrada seleccionamos las <br>
>> correspondientes al MDT que hemos cargado, y generaremos tanto la <br>
>> positive como la negative<br>
>><br>
>> Los parámetros de generación de openess por defecto funcionan <br>
>> bien, aunque sería posible modificarlos y observar el efecto en el <br>
>> mapa generado.<br>
>><br>
>> Por cierto, el proceso de generar la Topographic Openess puede <br>
>> durar un buen rato.<br>
>><br>
>> 4) Exportamos los resultados para poder cargarlos en QGIS:<br>
>><br>
>> Geoprocessing -> File -> GDAL/OGR -> GDAL: Export raster to GeoTiff<br>
>><br>
>> Seleccionamos la grid correspondiente y exportamos cada imagen <br>
>> por separado.<br>
>><br>
>> 5) Cargamos los ráster generados con SAGA en QGIS: Pendiente, <br>
>> Openess Positiva y Openess Negativa.<br>
>><br>
>> 6) Creamos la capa RVI (Ridge Valley Index) mediante la <br>
>> calculadora ráster. Recordad que RVI=(OpenessPos-OpenessNeg)/2<br>
>><br>
>> 7) Cargamos la capa RVI sobre la de Pendientes y le aplicamos <br>
>> como modo de mezcla Multiplicar, estableciendo los valores de <br>
>> máximo y mínimo de capa de manera adecuada (aquí yo he tenido que <br>
>> ir probando, aunque lo más adecuado ha sido, casi siempre, Estirar <br>
>> a min-max y cargar valores min-max)<br>
>><br>
>> 8) Coloreamos de rojo la capa de pendientes, estableciendo el <br>
>> gradiente de colores de la capa en "De blanco a negro"<br>
>><br>
>> 9) Ya tenemos el RRIM generado!!!<br>
>><br>
>> <br>
>><br>
>> Advertencias:<br>
>><br>
>> Este flujo de trabajo ha sido hecho para intentar reproducir <br>
>> el RRIM por alguien con bastantes lagunas en el manejo avanzado de <br>
>> GIS, por lo que puede ser considerado como un punto de partida más <br>
>> que otra cosa. Así, por ejemplo, con respecto al manejo de los <br>
>> valores mínimos y máximos de las capas, he tenido problemas para <br>
>> trabajar con una imagen del estrecho, ya que incluye un sector <br>
>> grande con NaNs<br>
>><br>
>> Aunque el modo de mezcla que más me resulta apropiado es el <br>
>> de Multiplicar, se pueden obtener diferentes resultados con otros <br>
>> modos de mezcla. Además, es posible que lo que especifica el método <br>
>> no se consiga estrictamente con modos de mezcla, sino por ejemplo <br>
>> creando una imagen multicanal en la que el canal rojo sea la <br>
>> pendiente y el canal de brillo el RVI. Crear semejante imagen en <br>
>> QGIS por ahora se me escapa y puede que esté diciendo incluso <br>
>> alguna barbaridad.<br>
>><br>
>> Espero que os pueda servir de inspiración y, si es de vuestro <br>
>> interés, se pueda ir mejorando<br>
>><br>
>> <br>
>><br>
>> Un saludo<br>
>><br>
>> <br>
>><br>
>> Jorge<br>
>><br>
>> <br>
>><br>
>> <br>
>> --<br>
>><br>
>> _______________________________________________<br>
>> QGIS-es mailing list<br>
>> <a href="mailto:QGIS-es@lists.osgeo.org" target="_blank" rel="noreferrer">QGIS-es@lists.osgeo.org</a><br>
>> <a href="https://lists.osgeo.org/mailman/listinfo/qgis-es" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-es</a><br>
_______________________________________________<br>
QGIS-es mailing list<br>
<a href="mailto:QGIS-es@lists.osgeo.org" target="_blank" rel="noreferrer">QGIS-es@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/qgis-es" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/qgis-es</a><br>
</blockquote></div>