Versions Compared

Old Version 17

changes.mady.by.user Jean Aroom

Saved on

compared with

New Version Current

changes.mady.by.user Jean Aroom

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
Comment: Migrated to Confluence 5.3

...

  1. By double-clicking the 'Raster' layer names and 'Tools' operations and typing the numbers, enter the following equation: "DEMft" - 30 * "FlowlineReclass" - 0.02 * (1500 - "FlowlineDistance") * ("FlowlineDistance" < 1500).
  2. For ‘Output raster’, rename the raster from “demft_raster” to “DEMRecon”.
  3. Ensure your Geoprocessing pane appears as shown below and click Run.



Notice that the DEMRecon layer 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.

...

  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the search box, type "flow".
  3. Click the Flow Direction tool.
  4. For ‘Input surface raster’, select the FIL layer.
  5. For ‘Output flow direction raster’, rename the raster from “FlowDir_FiIL1” to “FDR”.
  6. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

The result may appear a bit strange at first, due to the color scheme assigned to each cardinal direction.

 

 

  1. 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 more similar to a hillshade or 3D terrain visualization.

 

Calculating flow accumulation

Finally, you will calculate flow accumulation in each cell from the cells upstream, as illustrated below. Flow accumulation is dependent on the flow direction you calculated previously.

 

Image Modified

 

  1. Return to the Geoprocessing pane.
  2. At the top left of the Geoprocessing pane, click the Back button.
  3. In the search box, type "flow".
  4. Click the Flow Accumulation tool.
  5. For ‘Input flow direction raster’, select the FDR layer.
  6. For ‘Output accumulation raster’, rename the raster from “FlowAcc_FDR1” to “FAC”.
  7. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

  2. In the Contents pane, uncheck Turn off the Flowlines_StatePlane layer to better see the FAC layer.
  3. Right-click the FAC layer and click the select Symbology.
  4. For ‘Symbology’, use drop-down menu to select‘Primary symbology’, select Classify.
  5. In the ‘Method’ drop-down menu to For ‘Method’, select Equal Interval.
  6. Use the ‘Classes’ drop-down menu to select8 classes.
  7. For ‘Classes’, select 8.
  8. In the 'Classes' tableFor ‘Class Breaks’, double-click the numbers in the Upper Value column  and column and typethe following list: 100, 300, 1000, 3000, 10000, 30000, 100000. (Leave the final max value, 351408 as is.)
  9. Use the ‘Color Ramp’ drop-down menu to select the multipart color schemeFor ‘Color scheme’, scroll down to the very bottom and select the Yellow-Green-Blue (Continuous) color scheme, as shown below.

    Image Modified


Image Modified

 

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 southeast 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.

Image Removed

Zoom into the northern outlet point, as indicated on the image above, so that you can see the individual pixels.

  1. Zoom in to the southeastern outlet point, as indicated below, so that you can see the individual pixels.

Image Modified

Creating and editing feature classes

You will now create a new point feature class to store the outlet point you just identified.Click

  1. At the bottom of the Symbology pane, click the Catalog tab.
  2. Right-click the

...

  1. HydrologyLab.gdb geodatabase and select New > Feature Class

...

  1. .
  2. For

...

  1. ‘Name’, type “Outlet

...

Use the ‘Geometry Type’ drop-down menu to selectPoint.

For Coordinate System, selectCurrent Map.

ClickRun.

...

  1. ”.
  2. For ‘Feature Class Type’, select Point.
  3. At the bottom, click Next.
  4. Click Next again.
  5. For 'Current XY', ensure NAD 1983 StatePlane Texas S Central FIPS 4204 (US Feet) is selected.
  6. Click Finish.
  7. In the Catalog pane, drag the newly created Outlet feature class into the Lab3 map view.
  8. On the ribbon, click the Edit tab and click the Create button.

...


  1. Image Added

  2. On the Create Features pane, select the Outlet layer and click Point

...

  1.  button.

    Image Modified

  2. In the

...

  1. map view, click near the center of the outlet pixel to draw a point there.

    Image Modified

  2. On the

...

  1. ribbon, click the Save button.

    Image Added

  2. When asked if you want to save

...

  1. all edits, click

...

  1.  Yes.

Delineating watersheds

Now you will delineate the watershed that flows to the outlet point.

 

Image Modified

 Click Tools from Analysis tab

  1. At the bottom of the Create Features pane, click the Geoprocessing tab.
  2. At the top of the Geoprocessing pane, click the Back button.
  3. In the

...

  1. search box, type "watershed".
  2. Click the Watershed (Spatial Analyst Tools) tool.
  3. For ‘Input D8 flow direction raster’, select the FDR

...

  1. raster.
  2. For ‘Input raster or feature pour point data’ select the Outlet layer

...

  1. .
  2. For ‘Output raster’, rename the raster from “Watersh_FDR1” to “Watershed”.
  3. Ensure your Geoprocessing pane appears as shown below and clickRun.

...


  1. Image Added

  2. In the Contents pane, right-click the Watershed layer and select Zoom To Layer.

You can toggle the Subbasin_StatePlane layer on and off. Notice that the watershed boundary you have created

 

Notice that this is very close, but does not exactly match the boundaries of the subbasin you were originally provided from NHDPlus.

 

Image Added

 

  1. 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

  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the

...

  1. search box, type "calculator".
  2. Click the Raster Calculator (Spatial Analyst Tools) tool.
  3. Using the same clicking technique as before, enter the following equation: ("FAC" > 3000) & ("Watershed" > 0).
  4. For ‘Output raster’, rename the raster from “fac_rasterca” to “Streams”.
  5. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

 


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.

 

Image Modified

Now you will separate your entire raster of flowlines into individual flowline segments as illustrated below.

 

Image Modified

...


  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the

...

  1. search box, type "stream".
  2. Click the Stream Link tool.
  3. For ‘Input stream raster’, select the Streams

...

  1.  raster.
  2. For ‘Input flow direction raster’, select the FDR

...

  1.  raster.
  2. For ‘Output raster’, rename “StreamL_Stre1” to “StreamLinks”.
  3. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

  2. Symbolize the StreamLinks layer using unique values.
  3. 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.

...

  1. Return to the Geoprocessing pane.
  2. At the top left of the Geoprocessing pane, click the Back button.
  3. In the search box, type "watershed".
  4. Click the Watershed (Spatial Analyst Tools) tool.
  5. For ‘Input D8 flow direction raster’, select the FDR

...

  1.  raster.
  2. For ‘Input raster or feature pour point data’ select the StreamLinks layer.
  3. For ‘Output raster’, rename “Watersh_FDR2” to “

...

  1. CatchmentsRaster”.
  2. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

  2. Symbolize the

...

  1. CatchmentsRaster layer using unique values.

...

  1. In the Contents pane, collapse the

...

  1. CatchmentsRaster symbology.

Converting rasters to features

Now you will convert the Catchments CatchmentsRaster raster into a polygon feature class.

Click Tools from Analysis tab.

...

  1. Return to the Geoprocessing pane.
  2. At the top left of the Geoprocessing pane, click the Back button.
  3. In the search box, type "raster to polygon".
  4. Click the Raster to Polygon tool.
  5. For ‘Input raster’, select the

...

  1. CatchmentsRaster layer.
  2. For ‘Output polyline features’, rename “RatsterT_Catchme1” to “

...

  1. CatchmentsPoly”.
  2. Uncheck Simplify

...

  1. polygons.
  2. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

  2. Symbolize the CatchmentPoly layer using unique values with ID as the value field. Make sure to click the Add all values button.

...


  1. Image Added

  2. 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.

 

Image Modified

...

Click Tools from Analysis tab.

...

 

  1. Return to the Geoprocessing pane.
  2. At the top left of the Geoprocessing pane, click the Back button.
  3. In the

...

  1. search box, type "stream".
  2. Click the Stream to Feature tool.
  3. For ‘Input stream raster’, select the StreamLinks

...

  1.  raster.
  2. For ‘Input flow direction raster’, select the FDR

...

  1.  raster.
  2. For ‘Output polyline features’, rename the feature class from “StreamT_StreamL1” to “DrainageLines”.
  3. Uncheck Simplify polylines.
  4. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

...


Image Added

Calculating Strahler stream order

Now you will use the Strahler method to assign a flow hierarchy to the stream links as illustrated below.

Image Removed

...

 

Image Added


  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the

...

  1. search box, type "stream".
  2. Click the Stream Order tool.
  3. For ‘Input stream raster’, select the Streams

...

  1. raster.
  2. For ‘Input flow direction raster’, select the FDR

...

  1. raster.
  2. For ‘Output raster’, rename the raster from “StreamO_Stre1” to “StrahlerOrder”.
  3. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

  2. In the Contents pane, turn off all layers except for the new StrahlerOrder layer, so that its symbology is more apparent.

    Image Modified

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

  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the

...

  1. search box, type "zonal statistics".
  2. Click the Zonal Statistics as Table (Spatial Analyst) tool.
  3. For ‘Input raster or feature zone data’, select the StreamLinks layer

...

  1. .

...

  1. For ‘Zone field’

...

  1. , select the Value field.
  2. For ‘Input value raster’, select the StrahlerOrder

...

  1.  raster.
  2. For ‘Output table’, rename the table from “ZonalSt_StreamL1” to “OrderTable”.
  3. Ensure your Geoprocessing pane appears as shown below and click Run.

...


  1. Image Added

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.

  1. Open the DrainageLines attribute table.
  2. At the top of the DrainageLines attribute table, click the Add Field button.

...


  1. Image Added

  2. In the bottom row of the table, for ‘ Field Name

...

  1. :’, type

...

  1. “StrahlerOrder”.
  2. Use the

...

  1. ‘Data Type:’ drop-down menu to select Short

...

  1.  for short integer.
  2. Ensure your

...

  1. Fields table view appears as shown below

...

  1. .

    Image Added

  2. On the ribbon, under the Fields tab, click the Save button.

    Image Added

  3. Close the Fields table view.

Notice

 

Scroll all the way to the right and notice the new StrahlerOrder field has been added to the end of the DrainageLines attribute 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

  1. At the top right of the DrainageLines attribute table, click the Table Options button and select Joins and Relates > Add

...

  1. Join.

...


  1. Image Added

  2. For ‘Input Join Field’, select the grid_code field.
  3. For ‘Join Table’, select the OrderTable table.
  4. For ‘Output Join Field’, select the

...

  1. Value field.
  2. Ensure that your

...

  1. Add Join window appears as shown below and click

...

  1.  OK.

...


  1. Image Added

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.

  1. In the attribute table, right-click on the StrahlerOrder field header and select Calculate

...

  1. Field.
  2. Scroll down the list of

...

  1. ‘Fields’ and double-click the OrderTable.

...

  1. MIN field and click

...

  1. OK.
  2. Scroll down the newly populated StrahlerOrder field and notice the values of 1, 2, and 3 corresponding to Strahler order value.

Since you have copied the Strahler order into a field in the original DrainageLines attribute table, you You no longer need the tabular join.

  1. Again, click the Table Options button and select Joins and Relates

...

  1. > Remove All Joins.
  2. When asked if you are sure you want to remove all joins, click Yes.
  3. Symbolize the DrainageLines layer using graduated symbols based on the StrahlerOrder field.
  4. In the Contents pane, turn on only the Outlet, DrainageLines, and CatchmentPoly, layers and turn offall other layers.

...


  1. Image Added

 

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

...

  1. Create an 8.5 x 11 layout showing the DEMRecon layer

...

  1. .
  2. 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.