[QGIS-es] Generación de RRIM (Red Relief Image Maps)
carlos h
carloslda256 en gmail.com
Dom Feb 24 06:31:59 PST 2019
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 <jorge.tornero en ieo.es>
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 <msevilla00 en gmail.com> 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
> > <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
> _______________________________________________
> 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/20190224/b4bda159/attachment.html>
Más información sobre la lista de distribución QGIS-es