Hola muy interesante,una duda: Cuando se dice Cargamos la capa RVI sobre la de Pendientes y le aplicamos como modo de mezcla Multiplicar,
A qué se refiere o como se hace el modo de mezcla multiplicar? Disculpen mi ignorancia. Saludos Carlos El dom., 24 feb. 2019 3:18 a. m., JORGE TORNEO NUNEZ <[email protected]> escribió: > Buenas tardes, Miguel: > > Me alegro que la pequeña aportación vaya siendo de utilidad (aunque la > aportación de verdad la hicieron los inventores del RRIM). > Efectivamente, es posible, com dices, hacerlo directamente a través de > QGIS pero lamentablemente el QGIS del que disponía en ese momento y el > sistema en conjunto no tenía disponibles los algoritmos de SAGA para > hacerlo sin salirse de QGIS. La verdad es que las imágenes que se > generan son casi siempre espectaculares, digo casi siempre porque > probablemente dependan mucho de la resolución de los MDT de entrada. > Nosotros hicimos unas pruebas con batimetrías del golfo de Cádiz y el > resultado no era tan espectacular. Pero cuando los datos o el relieve > se prestan a ello, salen cosas muy chulas. > > Un saludo > > Jorge Tornero > > Miguel Sevilla-Callejo <[email protected]> escribió: > > > 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 > > <[email protected]> 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 > >> [email protected] > >> https://lists.osgeo.org/mailman/listinfo/qgis-es > _______________________________________________ > QGIS-es mailing list > [email protected] > https://lists.osgeo.org/mailman/listinfo/qgis-es >
_______________________________________________ QGIS-es mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/qgis-es
