Más

Guarde la capa vectorial y actualice el archivo de proyecto .qgs (PyQGIS)

Guarde la capa vectorial y actualice el archivo de proyecto .qgs (PyQGIS)


Necesita ayuda para resolver el siguiente caso de uso:

  1. QGIS comenzó con un nuevo proyecto (sin capas)
  2. Se crea una nueva capa mediante programación y se completan los datos:

    vl = QgsVectorLayer ("Punto", "puntos_temporales", "memoria") pr = vl.dataProvider () vl.startEditing () <… poner datos en vl> vl.commitChanges () QgsMapLayerRegistry.instance (). addMapLayer (vl)
  3. La nueva capa se guarda como shapefile:

    QgsVectorFileWriter.writeAsVectorFormat (vl, "c:  temp  my_shapes.shp", "CP1250", Ninguno, "ESRI Shapefile")
  4. Entonces supongo que necesito actualizar mi capa para cambiar el tipo de datos de "memoria" a "ogr" con la ruta al shapefile del paso # 3.

  5. Luego actualice el archivo del proyecto con esta información; por lo que la próxima vez que se cargue este proyecto, la capa leerá datos de my_shapes.shp.

"En cierto modo" sé cómo hacer estos pasos n. ° 4 y n. ° 5 usando la interfaz de usuario de QGIS:

  • guardar la capa como shapefile y establecer la casilla de verificación "Agregar capa guardada al mapa".
  • eliminar la capa original
  • Guardar proyecto

¿Debo implementar estos pasos programáticamente o hay una mejor manera?


¿Necesitas la capa de memoria?

Eche un vistazo al ejemplo "directamente desde características" en el libro de recetas de PyQGIS para ver cómo crear un Shapefile persistente a partir de características, para evitar crear una capa de memoria.


Abrir datos¶

Como parte de un ecosistema de software de código abierto, QGIS se basa en diferentes bibliotecas que, combinadas con sus propios proveedores, ofrecen capacidades para leer y, a menudo, escribir muchos formatos:

Los formatos de datos vectoriales incluyen formatos ESRI (Shapefile, Geodatabase ...), formatos de archivo MapInfo y MicroStation, AutoCAD DWG / DXF, GeoPackage, GeoJSON, GRASS, GPX, KML, valores separados por comas y muchos más ... Lea la lista completa de vectores OGR compatibles formatos.

Los formatos de datos ráster incluyen ArcInfo Binary Grid, ArcInfo ASCII Grid, JPEG, GeoTIFF, ERDAS IMAGINE, MBTiles, R o rásteres Idrisi, ASCII Gridded XYZ, GDAL Virtual, SRTM, Sentinel Data y muchos más ... Lea la lista completa de formatos compatibles con ráster .

Los formatos de base de datos incluyen PostgreSQL / PostGIS, SQLite / SpatiaLite, Oracle, DB2 o MSSQL Spatial, MySQL ...

El soporte de los servicios de datos web (WM (T) S, WFS, WCS, CSW, ArcGIS Servers…) también lo manejan los proveedores de QGIS (consulte QGIS como OGC Data Client).

También puede leer archivos compatibles de carpetas archivadas y usar formatos nativos de QGIS, como capas virtuales y de memoria.

A la fecha de este documento, los proveedores nativos GDAL / OGR y QGIS admiten más de 80 formatos vectoriales y 140 ráster.

No todos los formatos enumerados pueden funcionar en QGIS por varias razones. Por ejemplo, algunos requieren bibliotecas propietarias externas, o es posible que la instalación GDAL / OGR de su sistema operativo no se haya creado para admitir el formato que desea utilizar. Para ver la lista de formatos disponibles, ejecute la línea de comando ogrinfo --formats (para vector) y gdalinfo --formats (para ráster), o marque Configuración ‣ Opciones ‣ Menú GDAL (para ráster) en QGIS.

En QGIS, dependiendo del formato de datos, existen diferentes herramientas para abrirlo, principalmente disponibles en el menú Capa ‣ Agregar capa ‣ o desde la barra de herramientas Administrar capas (habilitada a través del menú Ver ‣ Barras de herramientas). Sin embargo, todas estas herramientas apuntan a un cuadro de diálogo único, el cuadro de diálogo Administrador de fuente de datos que puede abrir directamente con el Botón Abrir Administrador de fuentes de datos disponible en la barra de herramientas del Administrador de fuentes de datos o presionando Ctrl + L. De hecho, el cuadro de diálogo Administrador de fuentes de datos ofrece una interfaz unificada para abrir datos basados ​​en archivos vectoriales o ráster, así como bases de datos o servicios web compatibles con QGIS. Puede configurarse modal o no con el Diálogo del administrador de fuente de datos no modal en Configuración ‣ Opciones ‣ Menú general.

Cuadro de diálogo Administrador de fuente de datos de QGIS ¶

Además de este punto de entrada principal, también tiene el complemento DB Manager que ofrece capacidades avanzadas para analizar y manipular bases de datos conectadas. Se expone más información sobre las capacidades de DB Manager en el complemento de DB Manager.

También hay muchas otras herramientas, complementos nativos o de terceros, que lo ayudan a abrir formatos de datos dedicados.

Este capítulo describirá solo las herramientas proporcionadas por defecto en QGIS para cargar datos. Se centrará principalmente en el cuadro de diálogo Administrador de fuentes de datos, pero más que describir cada pestaña, también explorará las herramientas según el proveedor de datos o las especificidades del formato.


Vuelva a ejecutar el modelo 01 del tutorial para crear archivos de verificación

Una vez que esté listo, copie y guarde el Create_Empties.fvc archivo en el TUFLOWFV corre carpeta como trap_steady_04.fvc y abra este archivo en Notepad ++. Este archivo ahora solo contiene el comando para definir el formato GIS y la proyección del modelo. Editaremos este archivo para crear un nuevo archivo de ejecución.

Actualizar archivo FVC

1. Comente el último comando ya que no necesitamos generar archivos vacíos nuevamente.

Nota: Si ha completado las primeras 6 secciones de este tutorial y ha creado y ejecutado con éxito trap_steady_01.fvc de la carpeta TUFLOWFV runs, vaya al paso 5.

2. Ve a la carpeta Modelo_completo TUFLOWFV corre, copie el archivo FVC del tutorial original trap_steady_01.fvc a través de nuestra carpeta de trabajo el TUFLOWFV carpeta de ejecuciones.

3. Copie los datos de los límites de entrada estableQ_01.csv y estableWL_01.csv de Modelo_completo TUFLOWFV bc_dbase a nuestras carpetas de trabajo del modelo . TUFLOWFV bc_dbase

4. Copie el archivo de malla trap_steady_01.2dm de Complete_Model TUFLOWFV model geo a nuestra carpeta de malla de trabajo . TUFLOWFV model geo

5. Copiar a lo ancho TODAS comandos de trap_steady_01.fvc a la recién creada trap_steady_04.fvc archivo, debajo de las líneas de código existentes.

6. Cambie el tipo de salida a xmdf para facilitar la visualización al usar TUFLOW Viewer en QGIS:

7. Además, al final del archivo, agregue:

8. Ahora estamos listos para volver a ejecutar nuestra simulación con los archivos de verificación GIS habilitados. Ejecutaremos TUFLOW FV usando un archivo por lotes de Windows como se detalla en la siguiente sección.

Ejecutando desde un archivo por lotes.

El método preferido para ejecutar TUFLOW FV es mediante un archivo por lotes de Windows. Vamos a crear un archivo por lotes en el . TUFLOWFV corre carpeta llamada _run.bat. Para hacer esto, haga clic derecho dentro de la carpeta de ejecuciones, creando un nuevo archivo de texto de la siguiente manera:


Cambie el nombre del archivo de texto recién creado a _run.bat y ábralo en Notepad ++ (o su editor de texto preferido).


1. Una vez en un editor de texto, agregue este comando para establecer la ruta al ejecutable FV:

donde ruta es la ruta a su ejecutable TUFLOW FV.

2. Establezca la cantidad de subprocesos de CPU utilizados en la simulación según las especificaciones de su PC:

3. Agregue el siguiente comando para ejecutar el archivo FVC usando el ejecutable FV anterior:


4. Ejecute el modelo haciendo doble clic en el archivo por lotes. Si desea saber más sobre cómo ejecutar TUFLOW FV desde un lote, consulte nuestra página Wiki Ejecutando TUFLOW FV desde un archivo por lotes.

Inspeccione los archivos de verificación

Con GIS habilitado y usando el Escribir archivos de verificación == ..cheque comando, la versión TUFLOW FV 2019 puede exportar archivos GIS que le ayudarán a revisar lo que TUFLOW FV ha utilizado para la entrada. Comenzaremos investigando los tres archivos de verificación generados para trap_stready_04.fvc que se muestran a continuación.

1. Vaya a la carpeta TUFLOWFV run log, abra trap_steady_04.qgs archivo arrastrando este archivo a QGIS o simplemente haciendo doble clic en él.

Ahora todos los archivos de verificación se abren en las "Capas de salida". (Este modelo no tiene ninguna capa GIS de entrada, por lo que el grupo "Capas de entrada" está vacío).

2. Utilice el Aplicar estilo TUFLOW a capas abiertas para comprobar los estilos de las capas de comprobación.

Nota: Si no puede ver la malla en este punto, seleccione la capa "trap_steady_04_mesh_check (R)", vaya a ver & gt "acercar a la capa".

3. Haga clic derecho en el _mesh_check_R capa y abra la tabla de atributos.

La siguiente información se puede verificar a partir de los atributos:

  • IDX2: ID de celda 2D,
  • IDX3: ID de celda 3D (si el modo es 2D, IDX2 será igual a IDX3),
  • CTRDX: Coordenada X del centro de la celda,
  • CTRDY: Coordenada Y del centro de la celda,
  • Estera: ID del material de la celda,
  • ZB: Elevación del lecho del centro de la celda,
  • BRough: Rugosidad del lecho.


4. Desmarque la _ns_check_L check file en el Panel de capas para mostrar solo el archivo de verificación de malla y el archivo de verificación BC. Todas las caras de la celda límite están representadas por un segmento de línea. Selecciona el _bc_check_L comprobar archivo en el Panel de capas.


Utilice la herramienta Identificar entidad y haga clic en una de las caras de la celda para abrir una Ventana de atributos de entidad.

  • IDENTIFICACIÓN: ID de cadena de nodo asociado con el límite
  • TIPO: Tipo de condición de contorno. En este modelo "Q" para límite de flujo o "WL" para límite de nivel de agua.
  • Subtipo: Subtipo del límite
  • Fuente: Archivo de origen de límite vinculado con el límite



5. Desmarque la _bc_check_L y revisa el _ns_check_L en el panel de capas para revisar la ubicación y la dirección de las cadenas de nodos. Si las cadenas de nudos están configuradas correctamente, las flechas deben apuntar en la dirección del flujo.

Visualización de resultados de modelado con TUFLOW Viewer

Esta sección del tutorial nos mostrará cómo usar TUFLOW Viewer en QGIS para ver nuestros resultados XMDF.

1. Haga clic en el Visor TUFLOW (TuPlot) para abrir la ventana de TUFLOW Viewer. Vaya a Archivo & gt & gt & gtLoad Results - Map Outputs.

2. En la ventana emergente, navegue hasta TUFLOWFV resultados y abrir trap_steady_04.xmdf.sup

3. En TUFLOW Viewer, intente:

  • Hacer clic Alternar renderizado de malla para mostrar y ocultar la malla del modelo
  • Cambie el vector y el escalar para trazar en el Tipo de resultado panel
  • Use el control deslizante / botones de paso de tiempo para avanzar a través del tiempo

4. Haga clic con el botón derecho en el tipo de resultado y cambie el estilo de trazado en las pestañas Contorno / Vector / Malla.


CÉSPED¶

Configurar GRASS no es muy diferente de configurar SAGA. Primero, se debe definir la ruta a la carpeta GRASS, pero solo si está ejecutando Windows. Además, un intérprete de shell (normalmente msys.exe , que se puede encontrar en la mayoría de las distribuciones de GRASS para Windows) tiene que ser definido y su ruta configurada también.

De forma predeterminada, SEXTANTE intenta configurar su conector GRASS para usar la distribución GRASS que se envía junto con QGIS. Esto debería funcionar sin problemas en la mayoría de los sistemas, pero si tiene problemas, es posible que deba hacerlo manualmente. Además, si desea utilizar una versión de GRASS diferente, puede cambiar esa configuración y señalar la carpeta donde se guarda esa otra versión. GRASS 6.4 es necesario para que los algoritmos funcionen correctamente.

Si está ejecutando Linux, solo debe asegurarse de que GRASS esté correctamente instalado y de que pueda ejecutarse sin problemas desde una consola.

Los algoritmos de GRASS utilizan una región para los cálculos. Esta región se puede definir manualmente usando valores similares a los que se encuentran en la configuración SAGA, o automáticamente, tomando la extensión mínima que cubre todas las capas de entrada utilizadas para ejecutar el algoritmo cada vez. Si este es el comportamiento que prefiere, marque la casilla Usar región de cobertura mínima opción en los parámetros de configuración de GRASS.

GRASS incluye archivos de ayuda que describen cada algoritmo. Si configura el Carpeta de ayuda de GRASS , SEXTANTE los abrirá cuando utilice el [Mostrar ayuda] en la ventana de parámetros del algoritmo.

El último parámetro que debe configurarse está relacionado con el conjunto de mapas. Se necesita un conjunto de mapas para ejecutar GRASS, y SEXTANTE crea uno temporal para cada ejecución. Tienes que decirle a SEXTANTE si los datos con los que estás trabajando utilizan coordenadas geográficas (lat / lon) o proyectadas.

No se necesita configuración adicional para ejecutar algoritmos GDAL, ya que ya está incorporado a QGIS y SEXTANTE puede inferir su configuración a partir de él.


Métodos realizados en SAGA GIS para análisis hidrológico

Configuración básica de SAGA GIS

Antes de comenzar a configurar SAGA GIS para su uso, se recomienda que adquiera los datos que planea usar para su proyecto, tenga en cuenta dónde se encuentran sus datos en una carpeta y asegúrese de que todos los datos que está utilizando estén ubicados en ese carpeta.

Abriendo SAGA GIS

Para ejecutar sus análisis hidrológicos, asegúrese de que SAGA GIS esté instalado correctamente en su computadora y sea funcional. Si no tiene SAGA GIS, puede descargarlo en la página de inicio de SAGA GIS.

  • Si está instalado y es funcional, uno simplemente tiene que ubicar y hacer clic en el icono de SAGA GIS, esto abrirá la ventana inicial de SAGA GIS.
  • Una vez que se haya abierto SAGA GIS, su interfaz de usuario se mostrará como se muestra a continuación.

Figura 1. Captura de pantalla de cómo se ve la interfaz gráfica de usuario de SAGA GIS al inicio

Importando DEM

Esta sección le enseñará cómo incorporar adecuadamente un DEM procesado en SAGA GIS para que pueda realizar un análisis hidrológico de manera eficaz.

  • En la interfaz gráfica de usuario (gui) de SAGA GIS, primero deberá navegar a la Módulos menú desplegable ubicado dentro de la barra de herramientas del menú principal en la parte superior de la pantalla.
  • En el menú Módulos navegue hasta Archivo y luego a GDAL / OGR y abre el GDAL: Importar ráster mando.
  • Una vez que se haya abierto el cuadro de diálogo Importar ráster, verá que la primera fila tiene la palabra Archivo en la columna de la izquierda y nada en la columna de la derecha. Debe hacer clic en la columna vacía, lo que hará que el cuadro muestre un botón de puntos suspensivos (.). Luego, hará clic en este botón, que abrirá un explorador de archivos y navegará hasta su archivo DEM. Cuando haya localizado su archivo de modelo de elevación, haga clic en Abierto y luego, en el cuadro de diálogo Importar ráster, haga clic en Okey. Nota: Al descargar al importar un DEM en SAGA, puede estar en cualquier formato ráster (p. Ej., .Rst, .GRID, .tiff) pero se prefiere el formato tiff.
  • Después de haber iniciado el proceso de importación del archivo DEM, puede ver el progreso del módulo en la ventana de mensajes ubicada en la parte inferior de la interfaz gráfica de usuario de SAGA. Cuando aparece la ventana de mensaje La ejecución del módulo se realizó correctamente, su DEM se ha importado correctamente a SAGA.
  • Un aspecto que es importante recordar al importar cualquier archivo a SAGA, incluido el archivo ráster DEM, es que el archivo importado aparecerá en la pestaña de datos del espacio de trabajo. sin embargo, para ver visualmente sus datos, debe hacer doble clic en el nombre del archivo que abrirá una nueva ventana titulada Agregar capa al mapa seleccionado, luego seleccione nuevo y haga clic en OK.
  • Ahora, el DEM se ha importado de forma completa y precisa a SAGA GIS y se mostrará como en la Figura 2.

Figura 2. Modelo de elevación digital para el lago Okanagan, Columbia Británica, Canadá (Fuente: NRCan, 2013)

Guardar datos en SAGA GIS

Para que un usuario guarde los datos importados o los datos que se generan mediante el uso de una herramienta de análisis, debe guardar cada archivo ráster individualmente. Esto se debe a que los archivos no se guardan automáticamente en ninguna parte de la computadora y al guardar todo su trabajo como un proyecto SAGA (.sprj) no se guardan realmente los archivos ráster individuales y, si los archivos no se guardan, no se pueden volver a manipular o procesar. si el documento se cierra o se bloquea. Para guardar archivos ráster correctamente en SAGA GIS:

  • Navegue hasta el archivo ráster de la lista de archivos en el Datos pestaña de la espacio de trabajo.
  • Ahora haga clic derecho en el nombre del archivo y desplácese hacia abajo hasta Guardar como. .
  • Al hacer clic en Guardar como se abrirá una nueva ventana que le permite especificar dónde desea que se guarde el archivo, cómo desea nombrarlo y en qué formato debe estar. (Nota: la extensión .sgrd es el formato de cuadrícula nativo para SAGA GIS y generalmente funciona mejor).
  • Una vez especificada toda la información en este Guardar cuadrícula ventana, haga clic en Ahorrar. Sus datos ráster ahora se han guardado correctamente.

Proyectar los datos correctamente para realizar el análisis

En SAGA GIS, tendrá que comenzar sus análisis proyectando primero su archivo DEM en un sistema de proyección Universal Transverse Mercator (UTM) porque SAGA funciona mejor con este sistema de referencia de coordenadas en lo que respecta a las herramientas de hidrología. Si su DEM ya estaba en un sistema UTM y en la zona correcta para la región de estudio, puede omitir este paso; de lo contrario, debe seguir los siguientes pasos para convertir la proyección original del archivo DEM en un sistema de zona UTM que SAGA puede trabajar en.

Pasos para proyectar archivos en SAGA GIS:

  • Busque su código de proj.4 específico en línea para su área de interés en Spatialreference.org. Nota: Asegúrese de encontrar un código para su DEM que esté en NAD83 (CSRS) UTM y también asegúrese de que lo está proyectando en la zona UTM correcta para su región de estudio.
  • Una vez que tenga el código proj.4, regrese a SAGA GIS y abra el Módulos pestaña en la ventana del espacio de trabajo.
  • Navegue hasta el conjunto de herramientas titulado Proyección - Proj.4 y expandirlo haciendo clic en el pequeño botón con el signo más.

  • Ahora baja a la herramienta titulada Proj.4 (argumentos de la línea de comandos, cuadrícula).
  • Ahora verá un cuadro de diálogo emergente con dos campos principales, el primero llamado Parámetros de origen y el segundo llamado Parámetros de destino.
  • En la segunda fila de Parámetros de origen, denominada Sistema de cuadrícula, deberá usar el menú desplegable para seleccionar el sistema de cuadrícula bajo el cual se enumera su conjunto de datos de elevación (según lo que está escrito en la sección de datos de la ventana del espacio de trabajo.

  • A continuación, debe seleccionar su archivo DEM usando el menú desplegable de la fila titulada fuente.
  • Ingresar el archivo DEM en el cuadro de diálogo cambiará automáticamente la entrada de la primera fila llamada Parámetros de proyección de fuente a la proyección adecuada.
  • Ahora, en los parámetros de destino, escriba el código de proyección Proj.4 recuperado del sitio web de referencia espacial.
  • Deje la casilla Crear cuadrículas X / Y sin marcar y continúe con la siguiente fila.
  • En el archivo de Interpolación, seleccione cualquier interpolación de superficie que prefiera de la lista desplegable (Vecino natural se usa para este tutorial y generalmente es una buena idea dependiendo de su proyecto).
  • En el último campo del cuadro de diálogo, use el menú desplegable una vez más y seleccione usuario definido.
  • El siguiente paso en el proceso de proyección es hacer clic Okey.
  • Esto volverá a abrir una nueva ventana que le permite manipular la cuadrícula definida por el usuario para su proyección.
  • En esta ventana de usuario, las primeras cuatro filas ya tendrán entradas correctas basadas en sus parámetros de destino, pero deberá ingresar parámetros definidos en el campo de tamaño de celda.
  • Asegúrese de que esté definido el tamaño de celda adecuado; de lo contrario, su archivo DEM recién proyectado no funcionará correctamente para el análisis hidrológico. Su entrada de tamaño de celda debe ser el mismo que el tamaño de celda en el DEM original. Para este tutorial sobre Lago okanagan delimitación de la cuenca hidrográfica, el tamaño de celda se estableció en 30.

  • Después de esto, lo único que queda por hacer es hacer clic Okey.
  • Ahora ha vuelto a proyectar con éxito el DEM original en las coordenadas UTM adecuadas para su análisis hidrológico.
  • Esta herramienta de proyección proporciona al archivo un sistema de referencia de coordenadas interno (CRS); sin embargo, para aplicar un CRS al archivo ráster real, debe utilizar el Establecer sistema de referencia de coordenadas.
  • La herramienta Establecer sistema de referencia de coordenadas se encuentra debajo Módulo | Proyecciones | Establecer sistema de referencia de coordenadas.

  • En la ventana de herramientas, deberá ingresar el código EPSG para la proyección que se aplicó al DEM usando la herramienta de proyección Proj.4. Este código EPSG también se puede encontrar en el sitio web de referencia espacial.
  • Una vez que se haya definido el código EPSG, el campo Parámetros de Proj4 cambiará automáticamente según el código de proyección.
  • A continuación, baje al Cuadrículas campo bajo el Objetos de datos sección y haga clic en la columna vacía para mostrar el botón de exploración .
  • Al hacer clic en el botón Examinar, accederá a una ventana emergente interna en la que deberá seleccionar su archivo DEM. que ya proyectaste y transfiéralo a la columna de la derecha, que debe estar vacía.
  • Hacer clic Okey.
  • Hacer clic Okey también en el cuadro de diálogo principal de la herramienta Establecer sistema de referencia de coordenadas.
  • Finalmente, su aplicación de proyección de su DEM está completa y su archivo ráster está listo para análisis hidrológicos.

Llenar los fregaderos de un DEM

Es absolutamente vital para las personas que realizan análisis hidrológicos llenar los sumideros contenidos en los modelos digitales de elevación. Al completar estas características de sumidero, esencialmente está limitando el número de artefactos geométricos, errores de reflexión especular del agua y otros errores en el modelo de elevación causados ​​por una variedad de aspectos. Para llenar los lavabos, siga los siguientes pasos:

  • Ve a la Módulos menú en la barra de herramientas principal.
  • A continuación, navegue hasta el Análisis de terreno menú, amplíelo y luego coloque el mouse sobre Preprocesamiento.
  • Una vez que se haya abierto el menú de preprocesamiento junto al análisis del terreno, deberá ubicar la herramienta de llenado de sumideros titulada Llenar fregaderos (Wang Lui).

  • Una vez que haya abierto el Llenar fregaderos (Wang Lui) pestaña, selecciona su sistema de cuadrícula (bajo el cual se guardan sus datos DEM) a través del menú desplegable al lado de la fila titulada Sistema de red.
  • Después de seleccionar su sistema de cuadrícula, vaya a la fila titulada DEM y seleccione su archivo DEM a través del menú desplegable.
  • Una vez que haya hecho esto, asegúrese de que los campos DEM lleno,Direcciones de flujo, y Cuencas hidrográficas a través de sus menús desplegables tienen [crear] seleccionado.

  • El siguiente paso es hacer clic Okey, ya que generalmente no es necesario editar la fila de Pendiente mínima opcional.
  • Después de todos estos pasos, tendrá un nuevo archivo DEM en el que se llenan todos los lavabos, lo que hace que su DEM sea hidrológicamente correcto y esto se muestra en la Figura 3 para el Lago okanagan DEM.

Figura 3. Modelo de elevación digital para el lago Okanagan, Columbia Británica, Canadá después de que se hayan llenado todos los sumideros (Proyección: NAD83 (CSRS) UTM Zona 10N).

Mapas de orientación y mapas de pendientes

Los mapas de aspecto y pendientes se crearon en este tutorial del proyecto porque estos mapas son una excelente manera de visualizar la información DEM y le permiten ver cómo se ve la topografía en el área como si estuviera volando sobre la región en un avión. Para realizar estos análisis, siga los siguientes procedimientos:

  • Vaya a la ventana del espacio de trabajo de la GUI de SAGA GIS y desplácese hacia abajo hasta el conjunto de herramientas titulado Análisis del terreno: morfometría

  • Una vez allí, expanda el conjunto de herramientas usando el pequeño botón con el signo más al lado del nombre del conjunto de herramientas.
  • Haga doble clic Pendiente, aspecto, curvatura, y esto abrirá una ventana para el cuadro de diálogo de la herramienta de análisis.
  • Seleccione su sistema de cuadrícula, bajo el cual se encuentra su archivo de datos DEM, con el menú desplegable a la derecha de la ventana en la fila etiquetada Sistema de red.
  • Luego baja a la fila Elevación y vaya a la derecha y seleccione sus datos de elevación, que es su archivo DEM, a través del menú desplegable.
  • Para crear sus mapas de pendiente y aspecto, asegúrese de seleccionar [crear] en el menú desplegable de la derecha para esas respectivas filas.
  • Asegúrese de que las filas tituladas Curvatura, Planificar la curvatura, y Curvatura del perfil tengo [no establecido] seleccionados en sus respectivos menús desplegables.
  • En la fila titulada Método puede elegir el algoritmo que desee, sin embargo, se encontró que el Ajuste 2 grados Polynom (Zevenbergen & amp Thorne 1987) parece funcionar muy bien y ese es el empleado para la región del lago Okanagan.

  • Después de completar estos pasos, haga clic en Okey en la parte superior derecha de la ventana.
  • Esto creará tus mapas.

Figura 4. Mapa de aspecto del área de estudio del lago Okanagan.

Figura 5. Mapa de pendientes de la región de estudio del lago Okanagan, Columbia Británica, Canadá.

Conversión de un polígono vectorial en ráster

En esta sección del tutorial se ejecutarán los procedimientos necesarios para tomar un archivo que está en formato de polígono vectorial y convertirlo en formato ráster. Este paso es muy importante ya que se requiere para una de las entradas para la siguiente parte del análisis hidrológico, donde evaluamos la delimitación de la cuenca. Al crear una cuenca hidrográfica, se requiere que se especifique un área objetivo en función de la cual se construirá la cuenca hidrográfica, y esta área debe ser una característica hidrológica como un lago o río. Debido al hecho de que las características hidrológicas generalmente no se encuentran en formato ráster, es importante seguir los procedimientos de esta sección para convertir el polígono en ráster para un análisis SAGA GIS más preciso.

Para convertir vector a ráster:

  • Primero navega hasta el Módulos pestaña de la barra de herramientas del menú principal.
  • Abra la pestaña Módulos haciendo clic en ella y luego navegue hacia abajo en la lista desplegable para Red.
  • En el menú Cuadrícula que debería haberse abierto después de navegar a la Red selección, necesitará encontrar el Gridding conjunto de herramientas.
  • En este conjunto de herramientas, baje la lista hasta Formas a la cuadrícula herramienta y ábrala.

  • El cuadro de diálogo de la herramienta se abrirá ahora con muchas opciones de entrada.
  • Las dos primeras filas de la herramienta están etiquetadas Formas y Atributo y estas opciones le dan al usuario los campos en los que se debe especificar el archivo de polígono vectorial que necesita ser convertido. Entonces, en la primera fila, deberá usar el menú desplegable y seleccionar el archivo de polígono vectorial que debe convertirse. Para este proyecto, se seleccionó el polígono vectorial del lago Okanagan.
  • En la segunda fila, deberá seleccionar el atributo del archivo vectorial que determinará los límites del contorno del archivo. Nota: debido a que el archivo se está convirtiendo en este caso para la delimitación de la cuenca hidrográfica, se puede elegir cualquier atributo, y en este proyecto se utilizó el campo ID de objeto..
  • La siguiente fila titulada Método para valores múltiples tendrá que seleccionar el Significar opción del menú desplegable. Esto le dice al algoritmo de conversión que tome el promedio si se encuentran dos o más valores iguales en el archivo vectorial.
  • El campo llamado Método para líneas se puede especificar como desee. Simplemente significa que el contorno de su polígono vectorial será grueso o delgado. Para este proyecto se empleó la opción delgada.
  • Tipo de cuadrícula de destino preferido El campo variará según el tipo de archivo de datos vectoriales que tenga. En este proyecto, el polígono vectorial del lago Okanagan era un vector de 8 bytes, por lo que Punto flotante (8 bytes) se eligió la opción.
  • El ultimo campo Cuadrícula de destino debe configurarse en Red para que SAGA sepa que es necesario crear una nueva red.
  • La herramienta Formas a cuadrícula se puede ver a continuación para ayudarlo a comprender mejor los procedimientos.

  • Una vez que haya ingresado todos los parámetros requeridos para la herramienta, haga clic en Okey y su característica hidrológica de polígono vectorial se convertirá en un ráster. A continuación, se muestra el ráster del lago Okanagan después de convertirlo a partir de un polígono vectorial.

Crear un mapa de cuencas hidrográficas

En esta sección se analizará el análisis hidrológico real de la generación de mapas de cuencas hidrográficas para una región. En este tutorial específico el Lago Okanagan La cuenca hidrográfica se genera únicamente, pero los procedimientos se aplicarían a cualquier región de estudio. Para hacer un mapa de cuencas hidrográficas:

  • Vaya a la ventana del espacio de trabajo Módulos Tab y desplácese hacia abajo para Análisis del terreno - Hidrología y expandirlo.
  • Desplácese hacia abajo hasta Área de pendiente ascendente y abra la herramienta haciendo doble clic.

  • Seleccione su sistema de cuadrícula a través del menú desplegable en el lado derecho de la ventana en el Sistema de red fila .
  • Baja a la fila titulada Área objetivo y seleccione su área de interés, que es la característica hidrológica vectorial convertida en la que se basará la cuenca. En este tutorial, el área objetivo es el Límite del lago Okanagan archivo de trama.
  • Baja al Elevación fila y seleccione su DEM lleno datos del menú desplegable en el lado derecho de la ventana.
  • Asegúrese de que la fila titulada Rutas de fregadero posee [no establecido] seleccionado a través del menú desplegable. Este campo se especifica como [no establecido] porque el DEM lleno ya tiene en cuenta las rutas de sumidero de los datos.
  • Baja a la fila titulada Área de pendiente ascendente y asegúrate [crear] se selecciona en el menú desplegable. Este es el campo que le dice al algoritmo que genere el mapa de la cuenca hidrográfica, que también se conoce como mapa de área de pendiente ascendente.
  • Hundirse Opciones y a la fila titulada Método y asegúrate Determinista 8 se selecciona en el menú desplegable. Nota: En este proyecto se elige el algoritmo Deterministic 8 (D8) porque se considera que es la mejor opción al generar mapas de cuencas hidrográficas. Este algoritmo D8 asume que cada celda del DEM lleno tiene un flujo de agua unidireccional, mientras que los otros algoritmos asumen que una celda puede tener agua fluyendo en múltiples direcciones. Para hacer un mapa de delineación de cuencas hidrográficas más preciso y comprensible, el algoritmo D8 es el mejor.
  • Ahora el siguiente campo se titula Convergencia y para esta fila se aconseja que por favor, no hagas caso el campo y deje el valor predeterminado. Es importante dejar el campo intacto porque este factor de convergencia solo se aplica si está utilizando el dirección de flujo múltiple algoritmo, pero debido a que hemos utilizado el algoritmo D8, este campo no se aplicará a la salida.

  • Una vez que haya realizado todos los pasos, haga clic en Okey en la parte superior derecha de la ventana.
  • Esto crea su mapa de cuencas hidrográficas.

Figura 6. Cuenca del lago Okanagan, Columbia Británica, Canadá.

Ahora, como se ve en la Figura 6, se ha generado un excelente mapa de delineación de cuencas hidrográficas para la región de estudio, sin embargo, también puede ver que el mapa de cuencas hidrográficas se proyecta sobre un fondo blanco, lo que hace que la información sea ambigua y no proporcione un significado geoespacial. Para que los hidrólogos y / o investigadores en el mundo real utilicen este mapa con éxito, definitivamente necesitará colocar este mapa sobre otro mapa (como un DEM, mapa base de la región de estudio o mapa de aspecto) que proporcione contexto geoespacial a la información. para que los investigadores puedan ver qué áreas de la región de estudio contribuyen al flujo de agua al lago Okanagan. En este proyecto, se le dio un contexto geoespacial a la cuenca hidrográfica colocándola sobre el mapa de aspecto porque este mapa describe claramente el paisaje de la región de estudio y ayuda a identificar las áreas que contribuyen con el flujo de agua al área objetivo.

Figura 7. La cuenca del lago Okanagan, Columbia Británica, Canadá se muestra envuelta sobre el mapa de aspecto para que los datos sean más significativos.

En la Figura 7 puede ver el mapa de cuencas hidrográficas desplegado sobre el mapa de aspecto. Para cubrir la cuenca hidrográfica sobre otros mapas, primero debe convertir el archivo de cuenca hidrográfica ráster en un polígono vectorial y hacer esto:

  • Primero deberá abrir la herramienta llamada Vectorización de clases de cuadrícula navegando por los menús de la barra de herramientas principal siguiendo los siguientes pasos: Módulos | Formas | Cuadrícula | Vectorización | Vectorización de clases de cuadrícula.

  • Una vez que se haya abierto el cuadro de diálogo de la herramienta, deberá seleccionar una vez más su sistema de red de la lista desplegable de la fila respectiva.
  • Luego baja al Red campo e ingrese el Área de pendiente ascendente mapa que se creó para el mapa de cuencas hidrográficas.
  • En el campo titulado Polígonos puede asegurarse de que la columna de selección se establezca como [crear] para que el algoritmo sepa que es necesario generar un polígono vectorial.
  • El siguiente campo llamado La selección de clase debe establecerse como Todas las clases y el identificador de clase debe establecerse como 1. Este campo especificará cuántas clases desea que se generen en el polígono del vector de salida.
  • El último campo de entrada llamado Clase vectorizada como. debe establecerse en un solo objeto (multi) polígono porque queremos que se cree un solo polígono. La otra opción de este campo crearía un archivo vectorial del área de pendiente ascendente pero con muchos polígonos para construirlo, y usando la única opción obtendrá un polígono que se necesita para cubrirlo sobre el mapa de aspecto.
  • Ahora haga clic en Okey y se generará un nuevo archivo vectorial del mapa de la cuenca.

Los pasos anteriores describen los pasos necesarios para convertir un ráster en un polígono vectorial para que el mapa de la cuenca hidrográfica se pueda colocar sobre el mapa de aspecto. Una vez que el mapa de la cuenca hidrográfica se ha vectorizado, se pueden tomar más medidas para colocarlo sobre otros mapas para el contexto geoespacial. Para hacer esto:

  • Ve a la Datos pestaña de la ventana del espacio de trabajo y busque el polígono de vector de área de pendiente ascendente expediente.
  • Haga clic en el nombre del archivo para seleccionarlo y cuando esté resaltado vaya a la ventana de propiedades que se encuentra en el lado derecho de la pantalla.
  • En la ventana de propiedades, asegúrese de estar en el Ajustes y luego navegue hasta el campo titulado Transparencia [%].
  • En este campo puede ingresar cualquier valor que prefiera en función de qué tan transparente le gustaría que se vea su mapa de cuencas hidrográficas cuando se muestra sobre el mapa de aspecto. Para el área de estudio del lago Okanagan, una transparencia de 40 de seleccionados.
  • Una vez que se haya especificado el porcentaje de transparencia, haga clic en Solicitar y entonces Ahorrar en la parte inferior de la ventana de propiedades.
  • Ahora vuelve al Datos pestaña de la ventana del espacio de trabajo y haga clic con el botón derecho en el archivo de vector de cuenca.
  • Selecciona el Add layer to Map opción.
  • En el Map Selection Window that pops up select the Upslope area mapa.
  • Finally click OK and now your Watershed map is fully completed and has great geospatial information along with its hydrological information just like Figure 7.

Creating a Wetness Index Map

This section will run through the analysis of defining a region's wetness, which describes how susceptible specific areas in a study region are to becoming saturated if high inputs of precipitation occur in a relatively short amount of time. To run a wetness index analysis conduct the following:

  • In the workspace window under the Modules tab scroll down toTerrain Analysis - Hydrology' and expand it using the small plus sign button.
  • Scroll down to SAGA Wetness Index tool and double click it to open the dialogue box to the analysis tool dialogue box.

  • Select your grid system through the drop down menu on the right side of the window on the Grid System row .
  • Go down to the Elevation row and select your Filled DEM data from the drop down menu in the empty column beside the row title.
  • Ensure Catchment area, Catchment slope, Modified catchment area y Wetness index all have [create] in each of the row's drop down menu of the window.
  • In the row titled t make sure the row value is set to ten (10), since that is the default for the tool and it causes the wetness index values to being floating point numbers with a precision of ten points.

  • After you have done this click Okay at the top right of the window.
  • This creates your Wetness Index Map.

Figure 8. Map displaying the Wetness Index of the Lake Okanagan Region where red shows highest wetness and blue shows driest areas.