The Data menu provides options for bathymetry interpolation as well as accessing online bathymetry data. Being able to view bathymetry underneath the grid helps determine the resolution required in the deep and main channel regions compared to shallow regions, which can often have lower resolutions. In many cases, bathymetry data needs to be interpolated to cover the whole domain. This section describes features unique to Grid+ whereby a user may interpolate between two cross-sections along even over a sinusoidal domain. The Grid Bathymetry Interpolation and Download Online Bathymetry options are shown in Figure 1 and their use is described below.
| Anchor | ||||
|---|---|---|---|---|
|
Figure 1. Data menu options.
...
Simple Interpolation of Grid Bathymetry
...
In cases where cross-section data is available, it can be viewed in Grid+ with files in ldb format. For example, a file with name XS.ldb contains cross-sections with the following information (Figure 2).
- The first row is the name of the cross-section.
- The second row has a contains the number of data points (e.g 11), and a the number of columns (e.g 3) of for each data point. It This includes coordinates (X, Y) and Z (elevation)
- Data points start from the third row to end.
| Anchor | ||||
|---|---|---|---|---|
|
Figure 2. Cross-section data format.
From the Import button in the main toolbar, select Import Cross-Sections option, as shown in Figure 3. The Import Cross-Sections form will be displayed. Browse to the cross-section file, then click the Open button (Figure Data Menu#Data Menu#Figure 4), as a result, cross-sections are displayed as shown in Figure 5.
| Anchor | ||||
|---|---|---|---|---|
|
Figure 3. Import cross-sections option.
| Anchor | ||||
|---|---|---|---|---|
|
Figure 4. Import cross-section file.
The Coordinate System form will pop up. We need to select proper projection for the imported cross sections as shown in Data Menu#Data Menu#Figure 5. Click the OK button and the cross-sections will be displayed as shown in Data Menu#Data Menu#Figure 6.
Anchor Figure 5 Figure 5
Figure 4 Import cross-section file.
Figure 5 5. Define projection.
Anchor Figure 6 Figure 6
Figure 6. Cross-sections loaded.
To display the elevation of the cross-sections by color ramp, press Alt+C from the keyboard, and the Specify Elevation will pop up. Enter the Min and Max values for the elevation range (this range is based on the Z column value in the cross-section file). In this case, enter 0, -15 (Figure 6Data Menu#Data Menu#Figure 7). Then click the OK button, and the cross-section will be displayed as shown in Figure 7 Data Menu#Data Menu#Figure 8.
Anchor Figure 67 Figure 67
Figure 6 7. Cross-sections loaded.
Anchor Figure 78 Figure 78
Figure 7 8. Cross-sections displayed in color ramp.
Grid Bathymetry Interpolation for Dense Cross-Sections
For the river segment, the available In cases where a grid that covers the cross-sections are dense, so it is possible is available we can load the grid to interpolate grid bathymetry directly. First, draw a grid for the selected river segment, it is advised to create a fine grid to advance the interpolation results. Then, select the grid layer in the Layer Control panel, and onto it.
Go to the Data menu, then select Grid Bathymetry Interpolation option from the Data menu as shown in Figure 8 Data Menu#Data Menu#Figure 9. The Topographic Interpolation form will be displayed pop up as shown in Figure 9 Data Menu#Data Menu#Figure 10.
- Grid Layer : is the name of the selected grid layer that we are going to interpolate its for which we will interpolate the bathymetry.
- Using Visible Data: is the number of data points. In this case, there are 2234 points from in the cross-section file.
- Interpolation Options: There are two options, including Update Missing Only (update elevation values for the missing grid nodes only, in case other nodes already have elevation values) and Overwrite All Data (create new elevation values for all grid nodes and overwrite the nodes that have existing elevation values).
- Interpolation Method: Grid+ currently uses only one interpolation method, the Inverse Distance Weighting (IDW) is used. The IDW Options frame allows changing its of the interpolation properties.
Click the Interpolate button to
...
proceed with the interpolation process. When it is complete, click the OK button to close the form.
After the interpolation process is donecomplete, RMC on the grid layer and select Show Properties / Bottom Elevation to display the bottom elevation of the grid after its grid nodes are assigned elevation values, as shown in Figure 10 Data Menu#Data Menu#Figure 11.
| Info | ||
|---|---|---|
| ||
The bottom elevation values for the grid are for grid nodes, not for the centroid of grid cells. |
Anchor Figure 89 Figure 89
Figure 8 9. Grid Bathymetry Interpolation.
Anchor Figure 910 Figure 910
Figure 9 10. Topographic Interpolation form.
Anchor Figure 1011 Figure 1011
Figure 10 11. Grid with bottom elevation.
...
Select the Add a new polyline button from the main toolbar, then start drawing a polygon that covers the grid as shown in Figure 11 Data Menu#Data Menu#Figure 12. Change the cursor to select object mode (press the S key or select the Select Object in the main toolbar). Select the overlay layer (the layer that contains the polygon), then LMC on the polygon (it will be highlighted by changing color to red), then RMC to display options, and select the Interpolate Grid Elevations option as shown in Figure 12 Data Menu#Data Menu#Figure 13. A message will be shown as in Figure 13 Data Menu#Data Menu#Figure 14 when the interpolation process is complete. Click the OK button to close the message.
RMC on the grid layer and select Show Properties / Bottom Elevation to display the bottom elevation of the grid after its grid nodes are assigned elevation values, as shown in 2127790207 Data Menu#Data Menu#Figure 15.
Anchor Figure 1112 Figure 1112
Figure 11 12. Draw a polygon that covers the grid.
Anchor Figure 1213 Figure 1213
Figure 12 13. Interpolation by using a polygon.
Anchor Figure 1314 Figure 1314
Figure 13 14. Grid Bathymetry Interpolation message.
Anchor Figure 1415 Figure 1415
Figure 14 15. Grid with bottom elevation(2).
...
Other Interpolation Methods
In reality, the available cross-sections in some river segments are few ( as the example in Figure 15). Therefore, Grid+ support interpolates grid bathymetry by blocks and then extends to the entire grid. Firstly, draw a fine grid to cover the cross-section as shown in Figure 16. Then draw a polygon to cover the grid cells near the available cross-section
...
Figure 15 Sparse cross-sections in the river segment.
...
Figure 16 Grid to cover the cross-sections.
...
the second case, we will look at there are few cross-sections, and they are far apart from each other. The Inverse Distance Weighting (IDW) interpolation method used in the previous introduction is not an appropriate approach. Instead, interpolation using an orthogonal curvilinear grid will be used in this case.
Suppose that there are two cross-sections loaded as shown in Data Menu#Data Menu#Figure 16. Next, generate a grid that covers the two cross-sections. Note that the grid should have high resolution (grid cells).
The grid is displayed as shown in Data Menu#Data Menu#Figure 17. We can see that only some grid cells in two ends of the grid cover the cross sections. We should interpolate for the bathymetry of those cells first.
| Anchor | ||||
|---|---|---|---|---|
|
Figure 16. Cross-sections loaded and located in the river segment.
Anchor Figure 17 Figure 17
Figure 17. High-resolution grid covers the cross-sections.
Click Add a polyline button from the Main Toolbar to draw a polygon to cover the grid cells near the cross-sections, as shown in Data Menu#Data Menu#Figure 18. Then, change the cursor to select object mode (press the S key). Select the overlay layer (the layer that contains the new polygon), then LMC on the polygon to select it (it will be highlighted by changing its color to red), then RMC to display options, and select the Interpolate Grid Elevations option as shown in Data Menu#Data Menu#Figure 19. A message will be shown as in Data Menu#Data Menu#Figure 20 when the interpolation process is complete. Click the OK button to close the message. Apply similar steps to interpolate the bottom elevation of grid cells near the second cross-section.
Anchor Figure 18 Figure 18
Figure 18. Draw two polygons covering grid cells and cross-section.
| Anchor | ||||
|---|---|---|---|---|
|
Figure 19. Interpolate grid cells near the cross-section.
Anchor Figure 20 Figure 20
Figure 20. Grid Bathymetry Interpolation message.
To view the bottom elevation of cells that are interpolated, RMC on the grid layer and select Show Properties / Bottom Elevation as shown in Data Menu#Data Menu#Figure 21. The grid bottom elevation will be displayed as shown in Data Menu#Data Menu#Figure 22.
Next, we will interpolate the bottom elevation for the entire grid based on the cells with the existing bottom elevation in two ends of the grid. Select the grid layer in the Control Panel and set the cursor in selection mode. Select the first grid node and hold the Shift key to select the second node (select two grid nodes to define the grid block we want to interpolate). In this case, we select two grid nodes to cover the entire grid. The grid block will be highlighted by changing the color to red, then RMC to display options and select the Interpolate Elevations in I Direction option as shown in Data Menu#Data Menu#Figure 23. After the interpolation process is done, the grid bottom elevation is displayed as shown in Data Menu#Data Menu#Figure 24.
Anchor Figure 21 Figure 21
Figure 21. Show the bottom elevation of grid cells.
Anchor Figure 22 Figure 22
Figure 22. Grids with bottom elevation at two cross-sections.
Anchor Figure 23 Figure 23
Figure 23. Interpolate the grid block in I direction.
Anchor Figure 24 Figure 24
Figure 24. Grid with bottom elevation.
Extrapolation of Bathymetry
Besides this interpolation feature, Grid+ has an extrapolation feature. In this case, some cells in the two ends of the grid are missing the bottom elevation. To fill in their bottom elevation, select the grid layer in the Control Panel and set the cursor in selection mode. Select a grid block that covers both cells missing bottom elevation and cells having bottom elevation. The grid block will be highlighted by changing the color to red. Then RMC to the grid cells to display options and select the Extend Elevations in I Direction option as shown in Data Menu#Data Menu#Figure 25. The bottom elevation values will be extended to the cells that are missing values. As a result, the bottom elevation for those cells is assigned as shown in Data Menu#Data Menu#Figure 26.
Anchor Figure 25 Figure 25
Figure 25. Extend elevation in I-direction.
Anchor Figure 26 Figure 26
Figure 26. Grid cells after extension bottom elevation.
Exporting to a Scatter Data File
After interpolating the bathymetry of the grid cells, Grid+ allows exporting the bathymetry to a scatter data file. To do this, RMC on grid layer, then select Create Scatter Data from Elevation option as shown in Data Menu#Data Menu#Figure 27. This will extract the grid node’s elevation and then create an XYZ scatter layer in the Layer Control. It is named “XYZ from Grid”. RMC on XYZ from Grid layer, select Export Data Point ( as shown in Data Menu#Data Menu#Figure 28) and browse to the location to save the file, enter a file name, then click OK button to export the layer to a file.
Anchor Figure 27 Figure 27
Figure 27. Create scatter data from elevation.
Anchor Figure 28 Figure 28
Figure 28. Export bottom elevation to data file.
Download Online Bathymetry
After selecting the Download Online Bathymetry option, the Download Online Data form will be displayed as shown in 2127790207 Figure 29. Users are able to select data set using the drop-down Data Set. The currently available open-source sites are:
- GEBCO 2021: This is 2021 GEBCO’s gridded bathymetric data set, the GEBCO_2021 Grid, is a global terrain model for ocean and land, providing elevation data, in meters, on a 15 arc-second interval grid. The data information can be found on GEBCO's website GEBCO 2021.
- GEBCO 2014: This is 2014 GEBCO’s gridded bathymetric data set. The data information can be found on GEBCO's website GEBCO 2014.
- GEBCO 08: a continuous terrain model for ocean and land with a spatial resolution of 30 arc-seconds released by the General Bathymetric Chart of the Oceans (GEBCO) in 2009. Bathymetric data is available to download, and the user manual can be found on the website GEBCO 08.
Depending on the data set selected, different data information is provided in the Spatial Information section. The Data Extraction Limits is auto-set to the whole model domain, but the user can manually select the specific area by checking to Using Polygon File and browser to the closed polygon that contains the area inside.
Click on Download button and Grid+ will download and apply the bathymetry data of the selected Data Set to the specific defined area.
Anchor Figure 1529 Figure 1529
Figure 15 29. Download online data for bathymetry.
...