Difference between revisions of "Tutorial M11 001 GIS Inputs QGIS GPKG"

From Tuflow
Jump to navigation Jump to search
Line 162: Line 162:
  
 
<br>
 
<br>
[[File:Animation M11 001 GIS 07.gif]]<br>
+
[[File:Animation_M11_GPKG_001_GIS_07.gif]]<br>
 
<br>
 
<br>
 
</ol>
 
</ol>

Revision as of 13:51, 1 August 2023

Page Under Construction

Introduction

QGIS is used to create, import and view input data.

TUFLOW Project Re-Configuration

Re-configure the TUFLOW project to use and save empty files to the correct folder:

  1. Go to Plugins > TUFLOW > Editing > Configure / Create TUFLOW Project.
  2. Change the folder which contains TUFLOW to the Module_11 folder. This is the only change required, the model folder structure and template empty files are supplied from previous model.
    Note: The 'results' and 'check' folder get automatically created when the TUFLOW model is run.
  3. Click 'OK' and save the workspace by selecting Project > Save As.

    M11 GPKG ConfigureProject.PNG

  4. Set the workspace projection to EPSG:32760, see Set the Projection.


Method

Set up the GeoPackage Database:

  1. Navigate to the Module_11\TUFLOW\model\gis folder. Save a copy of M02_001.gpkg as M11_001.gpkg.
  2. Within the QGIS Browser Panel, right click on 'Favorites' and select 'Add a Directory...'.
  3. Navigate to the Module_11 folder and select it.

Navigate to the Module_11\Tutorial_Data folder:

  1. Copy and save the xs folder into the Module_11\TUFLOW\model folder.
  2. Within the QGIS Browser Panel, navigate to the Module_11\TUFLOW\model\xs folder. Drag and drop the following layers into the M11_001.gpkg:
    • 1d_xs_M11_bridges_002_L
    • 1d_xs_M11_channels_001_L
    • 1d_xs_M11_weirs_002_L

Within the QGIS Browser Panel, navigate to the Module_11\Tutorial_Data folder. Drag and drop the following layers into the M11_001.gpkg:

  • 1d_bc_M11_001_P
  • 1d_nwk_M03_culverts_001_L
  • 1d_nwk_M11_channels_001_L
  • 1d_wll_M11_001_L
  • 2d_bc_M03_culverts_001_P
  • 2d_bc_M11_hx_001_L
  • 2d_code_M11_null_creek_001_R
  • 2d_zsh_M11_hx_001_L
  • 2d_zsh_M11_hx_001_P

1D Cross Section

A 1D open channel requires cross sections for processing of the channel geometry. For each cross section there is an individual comma separated csv file:

  1. Load in the 1d_xs_M11_channels_001_L layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. The cross sections extend further than the 2d_code layer and DEM, here the model transitions from a 1D/2D model to a 1D only model. This method can be used where cross section data is available but topography data is not.
  5. Use the 'Identify Features' tool to inspect a cross section line. Each line has a unique name in the 'Source' attribute field that links to a cross section csv contained within the xs folder. Type attribute 'XZ' specifies distance vs elevation cross section:

    Animation M11 GPKG 001 GIS 01.gif

  6. Open the TUFLOW Viewer, the 1d_xs file automatically appears with a time series result 'XZ' showing. Select a 1d_xs line and view the cross section:

    Animation M11 001 GIS 02.gif

1D Channels

A 1d_nwk layer is used to model open channels. It contains information on inverts and to obtain data from the cross section the 1d_nwk and 1d_xs must be snapped:

  1. Load in the 1d_nwk_M11_channels_001_L layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. There are three 1d_nwk types used:
    • S - standard open channel.
    • W - weir.
    • X - connecting channel used for side channels entering the main channel.
    See Table 5-1 of the 2018 TUFLOW Manual for more information on the 1d_nwk types.
  5. Use the 'Identify Features' tool to inspect the attributes of the three 1d_nwk types:

    Animation M11 GPKG 001 GIS 03.gif

1D Boundary Conditions

The 1D model extents both upstream and downstream of the 2D active domain:

  1. Load in the 1d_bc_M11_001_P layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. Use the 'Identify Features' tool to inspect the attributes of the three boundaries. A QT type (flow vs time) hydrograph is used to apply flow at the upstream boundary and a HQ (water level vs flow) boundary is used to remove flow from the downstream boundary based on a water level vs flow curve:

    Animation M11 GPKG 001 GIS 09.gif

1D Culverts

1D culverts are included in a separate 1d_nwk layer. 1D culverts are described in detail in Module 3:

  1. Load in the 1d_nwk_M03_culverts_001_L layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. There are two 1d_nwk types used:
    • C - circular culvert.
    • R - rectangular culvert.
  5. Two culverts within the open channel link upstream 1D channel to the downstream 1D channel underneath the roads, the culvert to the east is embedded into the 2D domain directly. Use the 'Identify Features' tool to inspect the attributes:

    Animation M11 GPKG 001 GIS 04.gif

1D / 2D Links

The culvert on the east is connected to the 2D domain using 2d_bc layer SX type connection. It removes water out of the upstream cell(s) and deposits water back at the downstream cell(s):

  1. Load in the 2d_bc_M03_culverts_001_P layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. The 2d_bc SX points are snapped to the 1d_nwk culvert:

    Animation M11 GPKG 001 GIS 05.gif

The 1D open channel is linked to the 2D using a 2d_bc layer HX and CN type connections:

  1. Load in the 2d_bc_M11_hx_001_L layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. The 2d_bc HX type lines run along the top of bank. In HX type, the H indicates a head (water level) boundary and the X indicates a water level value coming from the 1D model. Once the water level in the 1D exceeds the water level in the 2D boundary cell, water can enter from 1D channel to 2D domain and vice versa. For water to flow in and out of floodplain overbank areas, HX lines extending into the floodplain at the upstream and downstream are required. Without them, the flow is assumed to be contained within bank which may not be the case.
  5. The CN lines connect the HX lines to the 1D channel. Where a CN line is connected to the HX line, the water level from the 1D nodes is transferred to the HX line. In between 1D nodes a linear interpolation of water level is applied.

    Animation M11 GPKG 001 GIS 06.gif

2D Breaklines

The water level is transferred between 1D to 2D through HX cells. It is important the 2D cell elevation reflects the correct level that water can spill out from one domain to another, thick breakline is required. 2D breaklines are described in detail in Module 2:

    2D Breaklines all together.jpg

  1. Load in the 2d_zsh_M11_hx_001_L and 2d_zsh_M11_hx_001_P layers from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layers, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layers (hold Ctrl to select multiple).
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. The 2d_zsh line layer is used to define breaklines along the bank of the channel:
    • No ‘Z’ attribute is specified, the elevations are taken from the 2d_zsh point layer.
    • The breakline has the same Shape_Width as the cell size, this is a thick breakline modifying whole cells.
    • The Shape_Option of 'MAX' is specified. The Zpt elevations only changes when the Z shape elevation at the Zpt is higher. The purpose of the Z shape is to enforce the highest elevations along the bank of the channel.
  5. The 2d_zsh point layer is used to apply elevation to the breaklines:
    • The values in the 'Z' attribute are elevations assigned to each point and are interpolated along the breakline.
    • The points must be snapped to the vertices of the 2d_zsh line to instruct TUFLOW that these layers are linked. If no points are in the layer, the elevation on the line object is used, i.e. the elevation along the line does not vary.

    Animation M11 GPKG 001 GIS 07.gif

Deactivate 2D Cells

The area of the 1D open channel defined by the 1D/2D HX boundaries must be removed from the active 2D cells:

  1. Load in the 2d_code_M11_null_creek_001_R layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. The 2d_code layer is used to deactivate the area representing the 1D open channel with attribute of '0'. All cells within the polygon are removed from the 2D computations.

    Animation M11 001 GIS 08.gif

Water Level Lines

Water Level Line (WLL) are used to display 1D results in map outputs. WLL do not change 1D hydraulic computations and are only used to display 1D results in plan (2D) formats. For more information on WLL, see Section 9.5.2 of the 2018 TUFLOW Manual:

  1. Load in the 1d_wll_M11_001_L layer from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layer, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layer.
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. The 1d_wll layer must be digitised perpendicular to the flow direction and from left to right looking downstream. An attribute of '2.5' is specified for all lines as a minimum distance interval along the WLL to generate triangulation points.

    Animation M11 001 GIS 10.gif

Reporting Locations

Reporting locations allow for the plotting of time series results across both the 1D and 2D sections of the model. Points are treated as water level output and lines as flow output.
Create a reporting location layer:

  1. Click on the ‘Import Empty File’ symbol from the TUFLOW Plugin toolbar.
  2. Select '0d_rl' from the Empty Type list.
  3. In the Run ID box write 'M11_001'.
  4. Tick the 'Line' and 'Point' Geometry Type.
  5. In the Spatial Database Options, select 'All to one'. Navigate to the Module_11/TUFLOW/model/gis folder and select M11_001.gpkg. Click Save and hit 'OK'.
  6. The 0d_rl (reporting location) files appear in the Layers panel.

Copy the line and point from the 2d_po layer:

  1. Load in the 2d_po_M01_001_L and 2d_po_M01_001_P layers from M11_001.gpkg, either:
    • Within the QGIS Browser Panel, navigate to M11_001.gpkg and double click the layers, or
    • In File Explorer, drag and drop the M11_001.gpkg into the QGIS workspace and select the layers (hold Ctrl to select multiple).
  2. Use the 'Apply TUFLOW Styles to Open Layers'.
  3. Use the 'Apply GPKG Layer Names'.
  4. Right click on the 0d_rl_M11_001_L and toggle on Editing.
  5. Copy the line from the 2d_po line layer into the 0d_rl line layer.
  6. Open the attribute table and enter Name as 'Loc1_L'.
  7. Turn off editing to save the edits.

    Animation M11 001 GIS 11e.gif

  8. Right click on the 0d_rl_M11_001_P and toggle on Editing.
  9. Copy the point from the 2d_po point layer into the 0d_rl point layer.
  10. Open the attribute table and enter Name as 'Loc1_P'.
  11. Turn off editing to save the edits.

    Animation M11 001 GIS 12e.gif

  12. To ensure the point reporting functionality is correct, it needs to be snapped to the start/end of a 1d_nwk line.
  13. Right click on the 0d_rl_M11_001_P and toggle on Editing.
  14. Use the 'Vertex Tool' to move the point.
  15. Enable snapping and use the 'Edit Advanced Configuration' tool to select 1d_nwk_M11_channels_001_L.
  16. Use the 'Vertex Tool' to snap the point to the start/end of the 1d_nwk line.
  17. Turn off editing to save the edits.

    Animation M11 001 GIS 13a.gif


Conclusion

  • A 1d_xs layer was inspected using the TUFLOW Viewer, cross sections were plotted for the 1D open channel.
  • A 1d_nwk layer was inspected, three types of 1D channel are used and connected to the 1D cross sections.
  • Boundary conditions were defined in the 1D using a 1d_bc layer.
  • A 1d_nwk layer was inspected containing information on 1D culverts.
  • The 1D open channel model was linked to the 2D using the 2d_bc 'HX' type.
  • Breaklines were introduced to enforce the top of bank elevations using 2d_zsh points and lines.
  • A 2d_code layer with zero attribute deactivating cells from the 2D domain was inspected.
  • A 1d_wll layer was inspected, lines were digitised perpendicular to flow and left to right looking downstream.
  • A 0d_rl layer was created to report flow and water level time series for the 1D and 2D domain.


Up
Go-up.png Back to Module 11 Main Page