...
- In the New values column, replace the “NODATA” value with a value of “0”.
- For ‘Output raster’, rename the raster from “Reclass_Flow1” to “Reclassify”.
- At the top of the Geoprocessing window, click the Environments tab on the right.
- Under the 'Processing Extent' section, for 'Extent', select Same As layer: DEMft.
- Under the 'Raster Analysis' section, for ‘Cell Size’, select Same as layer DEMft.
- For ‘Snap Raster’, select the DEMft raster.
- At the top of the Geoprocessing window, click the Parameters tab on the left.
- Ensure your Geoprocessing pane appears as shown below and click Run.
You will notice that the output raster is a rectangle that contains the geometry of encompasses the subbasin – , but does not match it exactly. follow its curved boundaries. In the next section (, 'Calculating Euclidean Distances) distances', you will assign specify a “mask” that will define a more clip the output raster to the specific geometry . The of the subbasin. In this case, the null values in this raster prevent this the input raster prevented a mask from being an option - , so we you will retroactively clip the reclassified raster.
- At the top left of the Geoprocessing pane, click the Back button.
- In the search box, type "clip".
- Click the Clip Raster tool.
- For ‘Input Raster’, select the Reclassify layer.
- Use the For ‘Output Extent’ drop-down menu to , select the Subbasin_StatePlane layer.
- Click Use Input Features for Clipping Geometry.
If this box was not checked, the output would be limited to the rectangular extent encompassing the subbasin, which is the same extent you already have.
- For ‘Output raster’, rename the raster from Reclassify_Clip to “FlowlineReclass”.
- Ensure that your Geoprocessing pane appears as shown below and clickRun.
Notice that the extent of the output file is now in the shape of the subbasin.
- Open the attribute table for the “Reclassify” Reclassify layer and note the ‘Count’ Count of cells with a value of 0.
- Now open the attribute table for the “FlowlineReclass” FlowlineReclass layer and note the ‘Count’ Count of cells with a value of 0.
Notice that, although the visual geometry of the two layers has changed (one is a rectangle, and one is the shape of the subbasin), the total number of cells has not changed. ArcGIS has not automatically updated the attribute table for the clipped raster.
- At the top left of the Geoprocessing pane, click the Back button.
- In the search box, type "build raster".
- Click the Build Raster Attribute Table tool.
- Drag in the “FlowlineReclass” For 'Input Raster', select the FlowlineReclass layer and click Run.
- Now Once again, open the “FlowlineReclass” attribute table for the FlowlineReclass layer again and note that the number of cells that have been assigned a 0 have decreased to reflect the clipped geometry.
- Remove “Reclassify”In the Contents pane, remove the Reclassify layer.
Now all cells that were previously assigned a value of NoData within the subbasin geometry have been assigned a value of 0. This will allow you to perform uniform mathematical operations across the entire subbasin in the future.
...
- At the top left of the Geoprocessing pane, click the Back button.
- In the search box, type "euclidean".
- Click the Euclidean Distance tool.
- For ‘Input raster or feature source data’, select the Flowlines_StatePlane layer.
- For ‘Output distance raster’, rename the raster from “EucDist_Flow1” to “FlowlineDistance”.
- Click the Environments… buttonAt the top of the Geoprocessing window, click the Environments tab on the right.
- For Processing 'Extent, use the ‘Extent’ drop-down menu to select Same as layer DEMft.Use the ‘Snap Raster’ drop-down menu to select the DEMft raster', select Same As layer: DEMft.
- For ‘Cell Size’, select Same as layer DEMft.
- For ‘Snap Raster’, select the DEMft raster.
In order to output a raster with the desired geometry (in this case the same as DEMft or subbasin) of the subbasin, you must apply a mask. Without a mask, the output raster will be a rectangle that contains containing the subbasin geometry., as seen before.
- For 'Mask’, Use the ‘Mask’ drop-down menu to select the Subbasin_StatePlane layer.
Go back to Parameters, ensure your Geoprocessing pane appears as shown below and click Run
Right-click FlowlineDistance layer and click Symbology.
Select Quantile for Method and 5 for classes.
You should now see cells that vary in color the farther away they are from the nearest flowline.
Performing map algebra
Finally, you are ready to perform the reconditioning step. You can think of this as digging a trench with a funnel-shaped cross-section along each flowline to ensure that water will flow in the intended direction. The math behind this step is explained in the diagram below and all units are in feet, since the State Plane Texas South Central project stores XY coordinates in feet and since you also recalculated the Z elevations in Lab 2 in feet.
At the top left of the Geoprocessing pane, click the Back button.
In the search box, type "calculator".
Click the Raster Calculator tool.
When entering expressions in the raster calculator, your syntax must be exact, so though you could technically type in the equation below, we’d recommend entering everything by clicking the mouse, so you do not get any syntax errors (at least until you’ve learned the proper syntax).
By double-clicking the layer names and single-clicking the calculator buttons, enter the following equation: "DEMft" - 30 * "FlowlineReclass" - 0.02 * (1500 - "FlowlineDistance") * ("FlowlineDistance" < 1500).
For ‘Output raster’, rename the raster from “demft_raster” to “DEMRecon”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
You should notice that the DEMRecon looks very similar to the original DEMft. To better understand the difference between the two, you will actually calculate the difference between their corresponding values.
Click the Tools tab.
Double-click the Raster Calculator tool again.
Using the same clicking technique as before, enter the following equation: "DEMft” – “DEMRecon”.
For ‘Output raster’, rename the raster to “DEMDiff”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
The DEMDiff layer displays the amount of earth that has been removed through the reconditioning process.
Part 3: Analyzing Hydrologic Terrain
...
The Fill tool fills sinks in the DEM raster to remove small imperfections in the data. If cells with higher elevations surround a cell, the water gets trapped in that cell and cannot flow. The Fill tool modifies the elevation of the original DEM to eliminate these problems.
At the top left of the Geoprocessing pane, click the Back button.
In the search box, type "fill".
Click the Fill tool.
For ‘Input surface raster’, select the DEMRecon layer.
For ‘Output surface raster’, rename the raster from “Fill_DEMReco1” to “FIL”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Calculating flow direction
Now you will calculate the cardinal direction in which water from each cell is expected to flow based on the elevations of surrounding cells. Each direction is assigned one of 8 integer values between 1 and 128 as indicated in the illustration below.
At the top left of the Geoprocessing pane, click the Back button.
In the search box, type "flow".
Click the Flow Direction tool.
For ‘Input surface raster’, select the FIL layer.
For ‘Output flow direction raster’, rename the raster from “FlowDir_FiIL1” to “FDR”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
The result may appear a bit strange at first, due to the color scheme assigned to each cardinal direction.
Symbolize the FDR layer using Unique Values with a smooth black and white gradient for the Color Scheme.
You will notice the layer now looks similar to a hillshade.
Calculating flow accumulation
Finally, you will calculate flow accumulation in each cell as illustrated below. Flow accumulation is dependent on the flow direction you calculated previously.
At the top left of the Geoprocessing pane, click the Back button.
In the search box, type "flow".
Click the Flow Accumulation tool.
For ‘Input flow direction raster’, select the FDR layer.
For ‘Output accumulation raster’, rename the raster from “FlowAcc_FDR1” to “FAC”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Turn off the Flowlines_StatePlane layer to better see the FAC layer.
Right-click the FAC layer and click the Symbology.
For ‘Symbology’, use drop-down menu to selectClassify.
In the ‘Method’ drop-down menu to select Equal Interval.
Use the ‘Classes’ drop-down menu to select8 classes.
For ‘Class Breaks’, double-click the numbers in the Upper Value column and typethe following list: 100, 300, 1000, 3000, 10000, 30000, 100000. (Leave the final max value, 351408 as is.)
Use the ‘Color Ramp’ drop-down menu to select the multipart color scheme.
Visually trace along the streams until you find the darkest path exiting the watershed. Notice that, according to the model you have just generated, the Buffalo-San Jacinto basin drains north into Lake Houston, instead of south into Galveston Bay. Look for the exact point on the map where this divergence is created. It is circled in the image below. A problem is that the topography in that area is extremely flat and, rather than a single well-defined channel, there are wide bodies of water connecting the bayous to Galveston Bay. Calculating the flow accumulation more accurately would require further editing the base raster so that the wide flat bodies of water are not treated as narrow flowlines. This example also illustrates how you must always double-check the results of your models to ensure they make sense. For the purposes of learning the remaining geoprocessing tools in this lab, you will continue with the existing flow accumulation model, as if it was correct.
Zoom into the northern outlet point, as indicated on the image above, so that you can see the individual pixels.
Creating and editing feature classes
You will now create a new point feature class to store the outlet point you just identified.
Click the Catalog tab.
Right-click the HydrologyLab geodatabase and selectNew > Feature Class….
For ‘Feature Class Name’, type “Outlet”.
Use the ‘Geometry Type’ drop-down menu to selectPoint.
For Coordinate System, selectCurrent Map.
ClickRun.
On the Standard toolbar, click the Edit tab and click the Create button.
On the Create Features pane, select Outlet layer and click Point. Click ->.
In the Map Display, click the outlet pixel to draw a point there.
On the Editor toolbar, click Save.
When asked if you want to save your edits, clickYes.
Close the Editor toolbar and the Create Features tab.
Delineating watersheds
Now you will delineate the watershed that flows to the outlet point.
Click Tools from Analysis tab.
In the Hydrology toolset, double-click the Watershed tool.
For ‘Input flow direction raster’, select the FDR layer.
For ‘Input raster or feature pour point data’ select the Outlet layer from the drop-down menu.
For ‘Output raster’, rename the raster from “Watersh_FDR1” to “Watershed”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Notice that this is very close, but does not exactly match the boundaries of the subbasin you were originally provided from NHDPlus.
In the Contents pane, turn off all layers beneath the Watershed layer, so that only the cells within the newly defined subbasin are visible.
Delineating streams
Now that you have delineated the subbasin boundaries, you will generate flowlines. The number of flowlines that are generated will depend upon the flow accumulation threshold you select in separating a flowline cell from a non-flowline cell. For the purposes of this lab, you will select a threshold of 3000.
Click Tools from Analysis tab.
In the Map Algebra toolset, double-click the Raster Calculator tool.
Using the same clicking technique as before, enter the following equation: ("FAC" > 3000) & ("Watershed" > 0).
For ‘Output raster’, rename the raster from “fac_rasterca” to “Streams”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Again you should have a binary raster, where flowline cells are assigned a value of 1 and all other cells are assigned a value of 0.
...
Now you will separate your entire raster of flowlines into individual flowline segments as illustrated below.
Click Tools from Analysis tab.
In the Hydrology toolset, double-click the Stream Link tool.
For ‘Input stream raster’, select the Streams layer.
For ‘Input flow direction raster’, select the FDR layer.
For ‘Output raster’, rename “StreamL_Stre1” to “StreamLinks”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Symbolize the StreamLinks layer using unique values.
In the Contents pane, collapse the StreamLinks symbology.
Delineating catchment basins
Now you will delineate the catchment basin associated with each stream link. These catchment basins are similar to the subwatersheds you downloaded from NHDPlus.
Click Tools from Analysis tab.
In the Hydrology toolset, double-click the Watershed tool.
For ‘Input D8 flow direction raster’, select the FDR layer.
For ‘Input raster or feature pour point data’ select the StreamLinks layer.
For ‘Output raster’, rename “Watersh_FDR2” to “Catchments”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Symbolize the Catchments layer using unique values.
In the Contents pane, collapse the Catchments symbology.
Converting rasters to features
Now you will convert the Catchments raster into a polygon feature class.
Click Tools from Analysis tab.
Double-click the Conversion Tools toolbox à From Raster toolset à Raster to Polygon tool.
For ‘Input raster’, select the Catchments layer.
For ‘Output polyline features’, rename “RatsterT_Catchme1” to “CatchmentPoly”.
UncheckSimplify polylines.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Symbolize the CatchmentPoly layer using unique values with ID as the value field.
In the Contents pane, collapse the CatchmentPoly symbology.
Now you will also convert the raster StreamLinks into a polyline feature class. While you could use the Raster to Polyline tool, you will instead use a tool in the Hydrology toolset, designed specifically for converting raster stream links into polyline features. The difference between the two methods is illustrated below.
Click Tools from Analysis tab.
Collapse the Conversion Tools toolbox.
In the Hydrology toolset, double-click the Stream to Feature tool.
For ‘Input stream raster’, select the StreamLinks layer.
For ‘Input flow direction raster’, select the FDR layer.
For ‘Output polyline features’, rename the feature class from “StreamT_StreamL1” to “DrainageLines”.
UncheckSimplify polylines.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Calculating Strahler stream order
Now you will use the Strahler method to assign a flow hierarchy to the stream links as illustrated below.
Click Tools from Analysis tab.
In the Hydrology toolset, double-click the Stream Order tool.
For ‘Input stream raster’, select the Streams layer.
For ‘Input flow direction raster’, select the FDR layer.
For ‘Output raster’, rename the raster from “StreamO_Stre1” to “StrahlerOrder”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
In the Contents pane, turn off all layers except for the new StrahlerOrder layer, so that its symbology is more apparent.
Calculating zonal statistics
Now you would like to append this Strahler order to the attribute table of the DrainageLines vector layer, so that you may symbolize that layer using graduated symbols according to the Strahler order number. First, you will need to create a table indicating the Strahler order of each stream.
Click Tools from Analysis tab.
In the Spatial Analyst Tools toolbox, double-click the Zonal toolbox > Zonal Statistics as Table.
For ‘Input raster or feature zone data’, select the StreamLinks layer from the drop-down menu.
Use the ‘Zone field’ drop-down menu to select the Value field.
For ‘Input value raster’, select the StrahlerOrder layer from the drop-down menu.
For ‘Output table’, rename the table from “ZonalSt_StreamL1” to “OrderTable”.
Ensure your Geoprocessing pane appears as shown below and clickRun.
Adding a new field to an attribute table
Now you will create a new field in the DrainageLines attribute table to store the Strahler order number.
Open the DrainageLines attribute table.
At the top of the DrainageLines attribute table, click the Add Field button.
For ‘Name:’, type “StrahlerOrder”.
Use the ‘Type:’ drop-down menu to select Short Integer.
Ensure your ‘Add Field’ window appears as shown below and click Save on top of the main tab.
Scroll all the way to the right and notice the new StrahlerOrder field has been added to the end of the table. All of the values are currently null, but will later be populated with the Strahler Order integers from the OrderTable.
...
You are now ready to join the OrderTable to the DrainageLines layer based on their common identification fields.
Click the Table Options button and select Joins and Relates > Add Join….
For ‘Input Join Field’, select the grid_code field.
For ‘Join Table’, select the OrderTable table.
For ‘Output Join Field’, select the VALUE field.
Ensure that your Geoprocessing pane appears as shown below and clickRun.
Scroll across the table and notice that most of the statistics (Min, Max, Mean, Majority, Minority, Median) are the same and all contain the Strahler order integer, so it doesn’t matter which field you copy, though you will choose the MIN field for this lab.
In the attribute table, right-click on the StrahlerOrder field header and select Calculate Field….
Scroll down the list of ‘Fields:’ and double-click the OrderTable.MINORITY field and click Run.
Scroll down the newly populated StrahlerOrder field and notice the values of 1, 2, and 3 corresponding to Strahler order value. You no longer need join.
Again, click the Table Options button and select Joins and Relates > Remove Join(s) > Remove All Joins.
Symbolize the DrainageLines layer using graduated symbols based on the StrahlerOrder field.
In the Contents pane, turn on only the Outlet, DrainageLines, and CatchmentPoly, layers and turn off all other layers.
FOR MAP LAYOUT TO BE TURNED IN
Create an 8.5 x 11 layout showing the vector catchment basins along with the stream links symbolized using graduated symbols according to their Strahler order and the outlet point that you used.
Deliverables
An exported map of DEMRecon layer turned on.
Create an 8.5 x 11 layout showing the vector catchment basins along with the stream links symbolized using graduated symbols according to their Strahler order and the outlet point that you used.