Versions Compared

Key

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

...

  1. On the Standard toolbar, click the Insert tab and click New Map.
    Image Removed
    Image Added

  2. At the top of the Contents pane, rename the map to "Lab2".

...

  1. In the Catalog pane, expand the Databases section and the HydrologyLab.gdb geodatabase.
  2. Drag the Watersheds feature class into the Lab2 map view.
  3. In the Contents pane, double-click the Watersheds layer.
  4. In the ‘Layer Properties’ window, click the Source tab.
  5. Scroll down and expand the Spatial Reference section.

Notice that the layer is in a geographic coordinate system called NAD 1983, which stands for North American Datum 1983. Because the data has a geographic coordinate system, the coordinates are stored in degrees, which indicate the three-dimensional location of the data on Earth's

...

spheroid. Though the data itself is stored in a geographic coordinate system, your computer monitor is flat, so, even though no projection has been defined, the data must be displayed in a particular projection. Whenever ArcGIS displays data in a geographic coordinate system, it uses a pseudo plate carrée projection, where one degree of latitude by one degree of longitude is represented as a square, rather than a curved trapezoid. In other words, all lines of latitude and longitude are evenly spaced. This type of projection results in stretching in the east-west direction, which increases the farther north or south from the equator you are mapping.

  1. Close the ‘LayerProperties’ window.

Working with geographic coordinate systems is fine for creating purely visual maps, as you did in Lab 1 (though the visual distortion can be disorienting and misleading), but, in this lab, you will be calculating areas, distances, and overlaps between features. Such calculations require the three-dimensional coordinates to be projected down onto a two-dimensional plane, so that the coordinates are stored in linear units, such as feet or meters, rather than degrees. In order to facilitate measurements of distance and area, you will now project the Watersheds layer into the State Plane Texas South Central projection which is best suited to mapping the greater Houston region.

  1. In the Analysis tab, click the Tools button to open the Geoprocessing pane.
  2. In the 'Find Tools' search box, type "project".
  3. Click the Under Toolboxes, double-click the Data Management Tools toolbox → Projections and Transformations toolset → Project tool.
  4. For ‘Input Dataset or Feature Class’, use the drop-down menu to select the Watersheds layer.
  5. For ‘Output Dataset or Feature Class’, rename “Watersheds Watersheds_Project” Projectto “Watersheds_StatePlane”, since that is the name of the projection you will be using.
  6. Next to the ‘Output Coordinate System’ box, click the Set Select Coordinate System button.
  7. Double-click Projected Coordinate Systems → State PlaneNAD 1983 (2011) (US Feet).
  8. Select NAD 1983 StatePlane Texas S Central FIPS 4204 (US Feet) and click OK.

Because both the input and output coordinate systems are based on the NAD 1983 geographic coordinate system, no geographic transformation is required.

  1. Ensure your ‘Project’ window appears as shown below and click Run.
    Image Removed
    Image Added

    Close Geoprocessing Pane.

Now that you have the correctly projected layer, you no longer need the original NAD 1983 layer.

  1. RightIn the Contents pane, right-click the original Watersheds layer and select Remove.
  2. Double-click the new Watersheds_StatePlane layer.
  3. Scroll down and expand the Spatial Reference section.

Notice that the layer is now in a projected coordinate system, NAD

...

1983

...

StatePlane

...

Texas

...

S Central FIPS  4204 (US Feet).

  1. Close the ‘Layer Properties’ window.

You may have noticed that the visual appearance of the watersheds did not change in your Map Display, even though you projected them. That is because the data frame takes on the projection of the first layer added to it. Since you first added the original unprojected Watersheds layer into the data frame, which was in NAD 1983, the data frame still displays the data in NAD 1983 (or

...

psuedo plate carrée). Currently, the projected Watersheds_StatePlane layer is being projected-on-the-fly back into NAD 1983 for visual purposes

...

.

Move the cursor around the screen and notice that the coordinates

...

at the bottom

...

of the

...

map view are shown in decimal degrees. This is another clue that the data frame is still using a geographic coordinate system; however, you would like the data frame to display data using the local Houston projection.

  1. At the top of the Contents pane,  rightdouble-click the Lab 2 map title and click Propertiesto open the 'Map Properties' window.

  2. Click the Coordinate System Systems tab.

While you could search for or navigate to the State Plane Texas South Central

...

projection, as you did before, in this case, you know that the same coordinate system is already used by the Watersheds_StatePlane layer. In such an instance, it is often easier to import the coordinate system from another known layer, especially if you are not familiar with the hierarchy of the coordinate system folders

  1. Scroll to the top of the 'XY Coordinate Systems Available' list.

  2. Expand Layers.

Notice that all of the coordinate systems used by layers currently in the map are displayed. You can expand the coordinate systems to determine exactly which layers are stored in which coordinate systems.

  1. On the ‘Coordinate System’ tab toolbar, click the Add Coordinate System button and select Import….
    Image Removed

  2. Double-click the HyrdologyLab geodatabase.
  3. Select the Watersheds_StatePlane feature class and click Add.
    Notice that the familiar State Plane Texas South Central projection is now listed

    Click NAD 1983 StatePlane Texas S Central FIPS  4204 (US Feet)

    .

  4. Click OK.

Notice that the watershed boundaries are now more compact in the east-west direction, as expected, because the local projection results in less distortion than the

...

pseudo plate carrée projection used to represent geographic coordinate systems.

Part 2: Downloading DEM Data

...

Now you are ready to download digital elevation model (DEM) data for the Buffalo-San Jacinto subbasin. Some local government agencies, such as the Houston-Galveston Area Council        Council (H-GAC) contract to have LiDAR data collected, which provides high resolution data containing both the elevation of the bare land and the heights of features in the built environment. Though we are lucky to have this high-quality data available in this particular region, for projects anywhere in the U.S., the best available DEM data generally comes from the National Elevation Dataset (NED) produced by the USGS. More information regarding NED data can be found at ned.usgs.gov. You will download NED data from The National Map ViewerDownload Client.

  1. In a web browser, go to https://viewerapps.nationalmap.gov/basicdownloader/.

First, you will constrain your downloads to your area of interest.

  1. In the 'Area of Interest' section, use the 'Map Extent/Geometry' drop-down menu button to select Selectable Polygon.
  2. Use the 'Select...' drop-down menu to select 8-digit HU (Subbasin).
  3. At the top right of the map, in the 'Find address or place' search box, type “Houston” and press Enter.
    Image Added

After zooming into Houston, you can see the blue boundaries of the individual subbasins and each subbasin is labeled with its HUC-8 number on the map.

    select the data products you are interested in viewing and downloading.
  1. On the left side bar, check the Elevation Products (3DEP) section to expand it.
  2. Within the Product Search Filter, check 1/3 arc-second DEM for the subcategory and select ArcGrid for the File Format.
  3. On the left side bar, scroll below the Elevation Products section and check the Hydrography (NHDPlus HR, NHD, WBD) section to expand it.
  4. Within the Product Search Filter, select National Hydrography Dataset (NHD) for the subcategory, select HU-8 Subbasin for the Data Extent and select FileGDB10.1 for the File Format.
    Now you constrain your downloads to your area of interest.
    In the search bar just above the map on the right, typeHouston” and click Go.
    After zooming into Houston, you can see the blue boundaries of the individual subbasins and each subbasin is labeled with its HUC-8 number on the map.
    Image RemovedNear the center of the map, click within  within Subbasin 12040104 to  to select it.

    Notice that just above the map, your search is now constrained to Polygon: huc8 12040104.
  5. At the top of the left side bar, click Find Products.
    Image Removed
  6. Under Available Products, to the right of Elevation Products (3DEP), click results.
    You are now provided with a listing of all the DEM tiles covering the area of Subbasin 12040104.
    Image Removed
  7. In turn, click Download for USGS NED n30w096, n30w095, and n31w096.
    Image Removed
  8. Navigate to the location where the three zipped folders have been downloaded.
  9. Select and copy all three zipped folders.

Your data search will now be constrained to the polygon for subbasin 12040104.

Next, you will select the data products you are interested in viewing and downloading.

  1. On the left side bar, under the 'Data' section, check the Elevation Products (3DEP) section to expand it.
  2. If necessary, within the 'Subcategories' section, check1/3 arc-second DEM.
  3. Under the 'File Formats' section, select All, which will make the ArcGrid file format available (rather than only the GeoTIFF and IMG file formats)
  4. Scroll back to the top of the left side bar and click the Search Products button.

    Image Added

You are now provided with a listing of all the DEM tiles covering the area of subbasin 12040104. Although 12 results are listed, there are only 3 areas of coverage provided in 4 different file formats: ArcGrid, GeoTIFF, GridFloat, and IMG. You will download only the 3 files in the ArcGrid format 

  1. In turn, click Download ZIP for USGS NED 1/3 arc-second 1 x 1 degree ArcGrid 2019 for n30w096, n30w095, and n31w096.

    Image Added

  2. Once all files have finished downloading, navigate to the location where the three zipped folders have been downloaded. The downloads may take a few minutes to complete.
  3. Select and copy all three zipped folders.
  4. Using Using Windows Explorer, navigate to your HydrologyLab folder.
  5. Paste all three zipped folders directly inside your HydrologyLab folder. Do NOT paste them inside the  HydrologyLab.gdb geodatabase folder.
  6. Select all three zipped folders.
  7. Right-click any of the three selected folders and select 7-Zip → Extract to “*\”, which will create one unzipped folder for the contents of each zipped folder., which will create one unzipped folder for the contents of each zipped folder. (If you are on a personal computer without 7-Zip installed, then right-click each folder in turn and select Extract All. Leave the default location, which is the same location as the original zipped folder, and click Extract.)
  8. Return to ArcGIS Pro.

Since you just added new files to your folder, you will need to refresh it in order for them to appear in the Catalog pane.

  1. At the bottom of the Geoprocessing pane, click the Catalog tab.
  2. Double-click Folders > HydrologyLab.
  3. If you do not see your newly download folders, right-click the 
  4. Click the Project tab on the Catalog pane.
  5. Right-click your HydrologyLab folder and select Refresh.
  6. Expand n30w095 and the other all three USGS folders to preview their contents.

Each raster file with a name such as grdn30w095_13 corresponds to a 1x1 degree tile. The file name contains “grd” for grid, followed by the latitude and longitude of the top left corner of the tile.

Ading Adding raster data in ArcMap

  1. Drag the grdn30w095_13 raster into the Map Display Lab2 map view.

You will be asked if you would like to create pyramids. Pyramids cache the raster at multiple reduced resolutions, resulting in an increased file size, but better rendering performance in your

...

map view. It is normally a good idea to create pyramids, but since you will be processing the corresponding rasters into a new file in the next step, you will not create pyramids at this time.

  1. In the 'Build pyramids' window, click No.
  2. Repeat Steps 1 and 2 the previous two steps with the grdn30w096_13 and grdn31w096_13 rasters.

Notice that the 1x1 degree tiles appear as angled rectangles, because they are being displayed in the State Plane Texas South Central projection. If the data frame were still in the NAD 1983 geographic coordinate system, the tiles would look like squares. Now you will look up the native coordinate system of the raster files.

  1. In the Contents pane, double-click the grdn31w096_13 layer to open the 'Layer Properties' window.Under
  2. In the Source tab, scroll down and notice the spatial reference is listed as GCS_North_American_1983 expand the 'Spatial Reference' section.

Notice the Geographic Coordinate System is NAD 1983 and that no projected coordinate system is listed. Since this coordinate system is not the same as the one currently used by the data frame, the raster layer is being projected-on-the-fly into the State Plane Texas South Central projection.

Mosaicking raster files

The edges where the three tiles meet are currently visible, because the minimum and maximum elevation is different in each raster, causing the same value to be represented with a different shade of gray in each raster. To solve that visual problem and also simplify future processing steps, you will mosaic the three rasters into a single raster. Before creating a mosaic, you will need to look up a few pieces of some information from the original rasters.

  1. Scroll back to the top of the ‘Layer Properties’ window and expand the 'Raster Information' section.

Notice that there is 1 band, meaning that each pixel only stores a single value: the height above sea level in meters. Aerial imagery has 3 bands to store the 3 RGB values. The format is listed as GRID, which stands for an Esri Grid, which is a raster file format native to Esri software. The pixel type and depth is 32-bit floating point-bit floating point, which indicates that the cells can store decimal data.

  1. Close the ‘Layer Layer Properties’ window.
  2. At the bottom of the Catalog pane, click the Geoprocessing tab.
  3. At the top left of the Geoprocessing pane, click the Back buttonOpen Tools. Click Toolboxes on the geoprocessing pane.
  4. In the search box, type "mosaic to new raster".
  5. Click the  Data Management Tools toolbox, click the Raster toolset → Raster Dataset toolset → Mosaic to New Raster tool.
  6. For ‘Input Rasters’, select the three 1x1 degree raster layers.For ‘Output Location’, click the Browse button1x1 degree raster layers.
  7. For ‘Output Location’, click the Browse button.
  8. In the bottom right corner of the 'Output Location' window, use the drop-down menu to select File Geodatabases (instead of Folders), if necessary.
  9. Select the HydrologyLab  HydrologyLab.gdb geodatabase and click Add OK.
  10. For ‘Raster Dataset Name with Extension’, typeDEMMosaic. No extension is necessary when storing the raster in a file geodatabase.
  11. Use the ‘Pixel Type’ drop-down menu to select 32 bit float, since that was the same type stored in the original rasters.
  12. For ‘Number of Bands’, type1”. (The text appears on the right side of the field.)
  13. Ensure your ‘Mosaic To New Raster’ pane appears as shown below and click Run.
    Image Removed
    Image Added

The mosaic may take a couple minutes to process. When it is complete, notice there are no longer visual seams between the tiles in the mosaicked raster. Now that you have a single mosaic, you no longer need the three originals tiles.

  1. RemoveIn the Contents pane, remove all three original raster files from the Contents panelayers.

Projecting raster files

Now you need to project your mosaic into the State Plane Texas South Central projection, so that you can properly calculate spatial statistics based on the data it contains.

  1. In At the top left of the Geoprocessing pane, click Data Management Tools -> Projections and Transformations toolset -> the Raster toolset Project Raster tool. click the Back button.
  2. In the search box, type "project raster".
  3. Click the Project Raster tool.
  4. Use the ‘Input Raster’ drop-down menu to select the DEMMosaic layerFor ‘Input Raster’, select in the DEMMosaic layer in the drop-down menu.
  5. For ‘Output Raster Dataset’, rename the raster from “DEMMosaic_ProjectRaster” to “DEMStatePlane” DEMMosaic_ProjectRaster to “DEMStatePlane
  6. Use the For ‘Output Coordinate System’ , drop-down menu to select Watersheds_StateplaneStatePlane, which will import the same coordinate system used by that layer.

You should now see the NAD_1983_StatePlane_Texas_South_Central_FIPS_4204_Feet displayed in for the Output Coordinate System.

  1. Use the ‘Resampling Technique’ drop-down menu to select Cubic Convolution.

For cell size, you may normally want to keep the original resolution, which was approximately 30 m, but, in this case, you will reduce the resolution to expedite processing times during this lab. The cell size is always specified in the same units as the projection, which in this case is feet, so you will select a cell size of 300 feet.

  1. Leave ‘Output Cell Size’ blank, but­ type “300” 300for the ‘X’ and ‘Y’. Make sure that your screen is the same as the screenshot below. Click
  2. Ensure your Geoprocessing pane appears as shown and click Run.
    Image Removed
    Image Added

Remember that the data frame was already displaying all layers in State Plane Texas South Central, so you should not notice much of a difference between the two layers, other than that the cell size has increased, meaning the resolution has decreased.

  1. Remove the DEMMosaic layer from the Contents pane.

...

Now you are ready to clip the DEM mosaic to the Buffalo-San Jacinto subbasin.

  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the search box, type "clip raster".
  3. Click the Clip Raster toolIn ArcToolbox, back in the Raster toolset, double-click the Raster Processing toolset and click the Clip (Data Management) tool.
  4. For ‘Input Raster’, drag in the DEMStatePlane layer from the Contents pane select DEMStatePlane.
  5. For ‘Output Extent’, drag in the select Watersheds_StatePlane layer.
  6. Check Use Input Features for Clipping Geometry.

Checking that box ensures that the raster is limited to the actual shape of the watersheds, rather than a rectangle covering the same extent.

  1. For ‘Output Raster Dataset’, rename the raster from “DEMStatePlaneDEMStatePlane_Clip” to “DEMSubbasin”Clip” to “DEMSubbasin.
  2. Ensure your ‘Clip’ window Geoprocessing pane appears exactly as shown below and click Run.
    Image Removed
    Image Added

  3. In the Contents pane, remove Remove the DEMStatePlane layer from the Contents pane.
  4. Turn off the Watersheds_StatePlane layer to view the DEMSubbasin layer.

...

Right now, the elevations in the DEM are in units of meters. In order to convert these units into feet, you will use the raster calculator.

  1. In At the top left of the Geoprocessing pane, click the Spatial Analyst Tools toolbox, click the Map Algebra toolset, click the Raster Calculator tool the Back button.
  2. In the search box, type "raster calculator".
  3. Click the Raster Calculator tool.
  4. In the list of ‘Rasters’, double-click the DEMSubbasin layer.
  5. In the calculatorlist of 'Tools', double-click the  the * button symbol.
  6. In the equation box, type “3 3.281”281, which is the conversion factor from meters to feet.
  7. For ‘Output raster’, rename the raster from “rastercalc” to “DEMft”DEMft.
  8. Ensure your ‘Raster Calculator’ window appears as shown below and click Run.
    Image Removed
    Image Added

Visually, the meters and feet layers should be identical, but, if you look at the layers in the Contents pane, you will notice that the original layer goes from elevations of -8 to 62 meters and the newly calculated layer goes from elevations of -26 to 204 feet.

  1. Remove the DEMSubbasin layer from the Contents pane.

...

Now you will determine which portions of land may be affected by a 15’ storm surge. Obviously, complex inundation models will take more variables into account, but, in this instance, you will simply highlight all the areas of land with an elevation of 15’ or less.

  1. Again, go back to In the Raster Calculator tool. Delete , delete the previous expression "DEMSubbasin" * 3.281.
  2. In the list of 'Rasters’, double-click the DEMft layer.
  3. In the calculatorlist of 'Tools', double-click the the <= button symbol.
  4. In the equation box, type “15” 15.
  5. For ‘Output raster’, rename the raster from “demftdemft_raster” raster to “Flood15ft”Flood15ft.
  6. Ensure your ‘Raster Calculator’ window Geoprocessing pane appears as shown below and clickRun.
    Image Removed
  7. Close Geoprocessing pane.

  8. Image Added

  9. In the Contents pane, rightRight-click the Flood15ft layer . Click and select Symbology.
  10. In the table of valuesSymbology pane, right-click the 0 value and click Remove.
    Image Removed
    Image Added

...


  1. Click the rectangle symbol to the left of the 1 value and select Blue.

Image Removed

11. Save your map document.


  1. Image Added

Now the cells containing an elevation of 15 feet or less are highlighted on top of the elevation rastertop of the elevation raster.

  1. In the Contents pane, uncheck the Watersheds_StatePlane layer and the Topographic basmap.
  2. Save your project.

FOR MAP LAYOUT TO BE TURNED IN

Create an 8.5 x 11 layout showing the DEM clipped to the subbasin with areas less than or equal to 15 feet in elevation highlighted.

Image Removed

...

Image Added

 

Symbolizing raster files

Now you will create a new map document from the one you are currently using.

  1. Click On the ribbon, click the Insert menu tab and select click the New Map…Map button.
  2. Double-click “Map1” in the contents pane and type “Lab2Topo” as name. Click Save.
  3. In the Contents pane, rename the map “Lab2Topo”.
  4. Return to Expand Databases in the Catalog pane. Click
  5. Expand the HydrologyLab.gdb geodatabase.
  6. Drag the DEMft and Watersheds_StateplaneStatePlane layers into the Lab2Topo map displayview.
  7. Turn off the In the Contents pane, uncheck the Watersheds_Stateplane layer.
  8. Right-click the DEMft layer and click the select Symbology.
  9. Scroll For ‘Color scheme’, scroll down the ‘Color Ramp’ drop-down menu to the very bottom and select the rainbow ramp Multipart Color Scheme or another continuous color scheme of your choice.
    Image Removed
    Image Added

Patterns within the data are now easier to see, especially at the lower elevations.

...

In addition to raster DEM data, it is sometimes useful to be able to represent elevation using vector contour lines. Contour lines at any regular intervals or discrete values can be created from DEM data.

  1. Open Return to the Geoprocessing pane by clicking Tools from Analysis tab.
  2. At the top left of the Geoprocessing pane, click the Back button.
  3. In the search box, type "contour".
  4. Click the Contour (Spatial Analyst Tools) In the Spatial Analyst Tools toolbox, click the Surface toolset then click the Contour tool.
  5. For ‘Input raster’, select the DEMft layer.
  6. For the ‘Output polyline features’, rename the feature class from “ContourContour_DEMft1” to “Contours10ft”DEMft1” to “Contours10ft.
  7. For ‘Contour interval’, type “10” 10.

The base contour is the lowest contour that will be shown. Since the lowest elevation is -26 ft, you will set your base contour to -20 ft.

  1. For ‘Base contour’, type  “-20”20.

Since the x,y XY coordinates are in the units of your State Plane Texas South Central projection, which is feet, and you have also used the Raster Calculator to convert the z Z units stored within the raster cells to feet, you do not need a custom Z factor.

  1. For ‘Z factor’, leave the default value of 1.
  2. Ensure your ‘Contour’ window Geoprocessing pane appears as shown below and clickRun.
    Image Removed
    Image Added

  3. Turn off the DEMft layer to better see the contours.

...

  1. Right-click the Contours10ft layer

...

  1. and select Symbology.

...

  1. In the Symbology pane, use the 'Primary symbology' drop-down menu

...

  1. to select Unique values.

...

  1. Use the

...

  1. ‘Field 1’ drop-down menu to select the Contour field.

...

  1. At the bottom of the ‘Color

...

  1. scheme’ drop-down menu

...

  1. check Show All and select

...

  1.  one of the

...

  1. continuous color

...

  1. schemes, as shown below. 

...


  1. Image Added

Now you can tell which contours are highest and lowest.14. Turn

  1. In the Contents pane, turn off and collapse the Contours10ft layer.

...

  1. Turn back on the DEMft layer.

Generating a hillshade

Now you will create a hillshade raster, which provides a shaded relief of the terrain based on a certain sun angle. It stores a value between 0 and 255 indicating the extent to which the cell would be shaded from the sun.

  1. Open Return to the Geoprocessing pane by clicking Tools from Analysis tab.
  2. In the Surface toolset, double-click the Hillshade tool.
  3. For ‘Input raster’, select the DEMft layer.
  4. For ‘Output raster’, rename the raster from “HillSha_DEMf1” to “Hillshade”.

The azimuth and altitude refer to sun angles. For now, you will stick with the default values.

As previously mentioned, the Z factor is used when the horizontal units of the projection and units of the elevation measurements stored in the raster cells are not the same. In this case, both are measured in feet, so a value of 1 is technically correct. Typically, hillshades are used to create realistic three-dimensional representations of mountains, which can be much easier for an audience to interpret than a flat DEM or contour lines. Unfortunately, the opposite is true in Houston, because the area is so flat.  It would be difficult to see any changes in elevation using a hillshade, since few shadows would be cast by the terrain. To compensate visually for the flat terrain, you will exaggerate the vertical elevations.

  1. For ‘Z factor’, type “20”.
  2. Ensure your ‘Hillshade’ window appears as shown and click Run.
    Image Removed

Typically hillshades are shown beneath transparent layers conveying other information, just to give the map a realistic appearance.

  1. In the Contents pane, drag the Hillshade layer beneath the DEMft layer.
  2. In the DEMft layer, in the feature-dependent Appearance tab, in the Effects group, slide the transparency to 60%.

Now the gradual changes in elevation are apparent from the coloring and are complimented by a realistic display of the terrain created with the hillshade.

  1. Turn on the Watersheds_StatePlane layer.

10. Symbolize the Watersheds_StatePlane layer with a hollow fill and thick black outline.

As a final touch, you will add the flowlines you downloaded in Lab1 to your Map Display.

11. Click the Catalog tab.

12. In the HydrologyLab geodatabase, drag the Flowlines feature class into the Map Display.

13. In the Contents pane, drag the Flowlines layer above the DEMft layer, but beneath the Watersheds_StatePlane layer.

14. Symbolize the Flowlines layer with a thin blue line.

15. Save your map document.

 

 

FOR MAP LAYOUT TO BE TURNED IN

Create an 8.5 x 11 layout showing transparent elevation in graduated colors on top of a hillshade, with watershed boundaries and flowlines visible.

Image Removed

 

 

Calculating zonal statistics

Now you will calculate elevation statistics by watershed.

...

  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the search box, type "hillshade".
  3. Click the Hillshade (Spatial Analyst Tools) tool.
  4. In the Surface toolset, double-click the Hillshade tool.
  5. For ‘Input raster’, select the DEMft layer.
  6. For ‘Output raster’, rename the raster from “HillSha_DEMf1” to “Hillshade”.

The azimuth and altitude refer to sun angles. For now, you will stick with the default values.

As previously mentioned, the Z factor is used when the horizontal units of the projection and units of the elevation measurements stored in the raster cells are not the same. In this case, both are measured in feet, so a value of 1 is technically correct. Typically, hillshades are used to create realistic three-dimensional representations of mountains, which can be much easier for an audience to interpret than a flat DEM or contour lines. Unfortunately, the opposite is true in Houston, because the area is so flat.  It would be difficult to see any changes in elevation using a hillshade, since few shadows would be cast by the terrain itself. To compensate visually for the flat terrain, you will exaggerate the vertical elevations.

  1. For ‘Z factor’, type20”.
  2. Ensure your ‘Hillshade’ window appears as shown and click Run.

    Image Added

Typically, hillshades are shown beneath transparent layers conveying other information, just to give the map a realistic appearance.

  1. In the Contents pane, drag the Hillshade layer beneath the DEMft layer, but above the Topographic basemap.
  2. In the Contents pane, select the DEMft layer.
  3. In the ribbon, click the contextual Appearance tab.
  4. In the Effects group, slide the transparency to 60%.

Now the gradual changes in elevation are apparent from the coloring and are complimented by a realistic display of the terrain created with the hillshade.

  1. Turn on the Watersheds_StatePlane layer.
  2. Symbolize the Watersheds_StatePlane layer with no transparency, a hollow fill, and thick black outline.

As a final touch, you will add the flowlines you downloaded in Lab 1 to your Map Display.

  1. Click the Catalog tab.
  2. In the HydrologyLab geodatabase, drag the Flowlines feature class into the Lap2Topo map view.
  3. In the Contents pane, drag the Flowlines layer above the DEMft layer, but beneath the Watersheds_StatePlane layer.
  4. Symbolize the Flowlines layer with a thin blue line.
  5. Save your project.

FOR MAP LAYOUT TO BE TURNED IN

Create an 8.5 x 11 layout showing transparent elevation in graduated colors on top of a hillshade, with watershed boundaries and flowlines visible.

Image Added

Calculating zonal statistics

Now you will calculate elevation statistics by watershed.

  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 "zonal statistics".
  4. Click the Zonal Statistics as Table (Spatial Analyst Tools) tool.
  5. Use the ‘Input raster or feature zone data’ drop-down menu to select the Watersheds_StatePlane layer.
  6. For ‘Zone field’, select the ‘HU_10_NAME’ field.
  7. For ‘Input value raster, select the DEMft layer.
  8. For ‘Output table’, rename the table from “ZonalSt_Watersh1” to “WatershedElevation”.
  9. Ensure your Geoprocessing pane appears as shown and click Run.

    Image Added

  10. At the bottom of the Contents pane, right-click the WatershedElevation table and select Open.

This table tells you the statistics regarding all the elevation values within each watershed zone.

  1. In the Geoprocessing pane, search for "table to excel’.
  2. Click the Table To Excel tool.
  3. For 'Input Table', select WatershedElevation.
  4. For ‘Output Excel File’, click the Browse button.
  5. Locate and double-click the HydrologyLab folder to save the file within the folder.
  6. For ‘Name’, type “WatershedElevation”.
  7. Click Save.
  8. In the Geoprocessing pane, click Run.
  9. Close the WatershedElevation table.

Now you will open the exported text file in Excel.

  1. Using File Explorer, open your WatershedElevation.xls spreadsheet in Microsoft Excel.
  2. Delete the OBJECTID, ZONE_CODE, COUNT, AREA, and SUM fields.
  3. Highlight the minimum and maximum values in each column.

Continue formatting the table until you are satisfied with its appearance.

  1. Click the File menu and select Save As.
  2. Navigate to your HydrologyLab folder.
  3. For ‘File name:’, typeWatershedElevation”.
  4. Use the ‘Save as type:’ drop-down menu to select Excel Workbook (*.xlsx).
  5. Click Save and close Excel.

FOR TABLE TO BE TURNED IN

Create a table highlighting the minimum and maximum values for the minimum, maximum, range, mean, and standard deviation of all elevation values within each watershed.

Image Added

Part 3: Downloading rainfall data

Downloading data from Climate Data Online.

Now you will download rain gauge station data created by the National Climatic Data Center (NCDC) using the Climate Data Online (CDO) interface.

  1. In a web browser, go to www.ncdc.noaa.gov/cdo-web/.
  2. Click the Mapping Tool tab.

    Image Added

You will search for data using the NHD hydrologic units. Previously, you had been working with the Buffalo-San Jacinto subbasin (HUC = 12040104). In this case, you will step up two levels to the Galveston Bay-San Jacinto subregion (HUC = 1204). 

  1. On the Surface Maps tab, click Normals.

    Image Added

  2. In the left sidebar, on the Layers tab, uncheck Daily Climate Normals, and check Annual Climate Normals.
  3. To the right of Annual Climate Normals, click the Map Tools button.

    Image Added

  4. In the new ‘ANNUAL CLIMATE NORMALS TOOLS’ window, click Location.
  5. Use the drop-down menu to select USGS HUC.
  6. Use the ‘Select a HUC type’ drop-down menu to select Subregions (4-digit).
  7. Use the ‘Select a HUC’ drop-down menu to select Galveston Bay-San Jacinto.

    Image Added

  8. Click Zoom to location.
  9. The left sidebar switches to the Results tab. Click Get All Location Data.

    Image Added


  10. For Step 1, select Custom Annual/Seasonal Normals CSV for the output format and click CONTINUE.

    Image Added


  11. For ‘Station Detail & Data Flag Options’, check Station Name, Geographic Location, and Include Data Flags to include those variables the data table.

    Image Added

  12. For ‘Select data types for custom output’, expand the Annual Precipitation category.
  13. Check Long-term averages of annual precipitation totals (ANN-PRCP-NORMAL).

    Image Added


  14. At the bottom of the window, click CONTINUE.
  15. Type your email address twice and click SUBMIT ORDER.

Check your email. You should receive two emails a couple minutes apart, although it may take a few hours to receive the second email. The first one indicates that your data request was submitted and the second one includes the requested data.

  1. In your email, click the Download link to download the requested CSV file.

    Image Added

Excel

  1. Navigate to the location where the CSV file was stored.
  2. Double-click the CSV file to open it using Excel.

The first column contains the unique station identification code and the second column contains the station name. Next are the latitude, longitude, and elevation of the stations. The annual precipitation field contains long-term averages of annual precipitation totals in hundredths of inches. More information is available on the Datasets portion of the CDO website.

 

Info

As of March 7, 2021, the CSV download includes ALL climate variables, even though only ANN-PRCP-NORMAL was selected for download. In this case, the annual precipitation field, ANN-PRCP-NORMAL, will not appear in column F, after the ELEVATION field. Instead, it appears in column BP. In this case, you will need to copy and paste the data from column BP to column F.


Before opening this table in ArcGIS, you must reformat some of the field names, which cannot have special characters and are recommended to be 13 characters or less.

  1. Rename “ANN-PRCP-NORMAL” to “ANNPRCP_HI”, for annual precipitation in hundredths of inches.
  2. Along the top of the worksheet, drag across the column letters to select columns A through F.
  3. Copy the selected columns.
  4. At the bottom left of the worksheet, click the New sheet button (plus) to create a new worksheet. In cell A1, paste your previously copied columns. 
  5. Delete the previous sheet. 
  6. Along the top of the worksheet, drag across the column letters to select columns A through F.
  7. Hover your mouse between columns E and F until the cursor changes to two outward facing arrows and double-click to auto-size the column widths.
  8. At the bottom left of the worksheet, rename the worksheet from Sheet1 to “PrecipStations”.
  9. Click the File menu and select Save As.
  10. Navigate to your HydrologyLab folder.
  11. For ‘File name:’, type "PrecipStations".
  12. Use the ‘Save as type:’ drop-down menu to select Excel Workbook.
  13. Click Save.
  14. Close Excel.

Displaying XY data

Now you are ready to start a new map and display the tabular rain gage data you just downloaded.

  1. Return to ArcGIS Pro.
  2. On the ribbon, click the Insert tab and click the New Map button.
  3. In the Contents pane, rename the map “Lab2Precip”.
  4. Return to the Catalog pane.
  5. Drag Watersheds_StatePlane into the Lab2Precip map view.
  6. Right-click the HydrologyLab folder and select Refresh.
  7. Expand the PrecipStations.xlsx Excel file to see the individual worksheets it contains.
  8. Drag the PrecipStations$ worksheet into the Lab2Precip map view.
  9. In the Contents pane, right-click the PrecipStations$ table and select Display XY Data.
  10. For ‘X Field:’, select the LONGITUDE field.
  11. For ‘Y Field:’, select the LATITUDE field.
  12. For 'Coordinate System', click the Select coordinate system button.

Because the coordinates are in the form of latitude and longitude in decimal degrees, you know you will need to select a geographic coordinate system, rather than a projected coordinate system.  While the data could theoretically be in any geographic coordinate system, you will select the North American Datum 1983, commonly abbreviated NAD 83, because this is coordinate system of the data provided on the NCDC website.

  1. Scroll towards the top of the 'XY Coordinate Systems Available'.
  2. Ensure that Geographic Coordinate Systems is already expanded, then expand North America > USA and Territories.
  3. Select NAD 1983 and click OK.
  4. Ensure that your window matches that below and click OK.

    Image Added


The points should now appear on top of the watersheds, though they also extend beyond the watersheds in the Buffalo-San Jacinto subbasin, since we downloaded them for the entire Galveston Bay-San Jacinto subregion.

  1. In the Contents pane, right-click the PrecipStations$ table and select Remove

...

  1. In the Contents pane, right-click the WatershedElevation table and select Open.

This table tells you the statistics regarding all the elevation values within each watershed zone.

  1. In the Geoprocessing pane, search for ‘Table to Excel’.

10. For the Input Table, select ‘WatershedElevation’ from the drop-down menu.

11. For ‘Output Excel File’, click the Browse button.

If you save this table inside your file geodatabase, you will not be able to open it outside of ArcGIS, so you will instead save it as a text file outside of your file geodatabase.

12. Save the table within the HydrologyLab folder.

13. For ‘Name’, name as “WatershedElevation”.

14. In the ‘Table to Excel’ window, click Run.

15. Close the Table.

Now you will open the exported text file in Excel.

  1. Open WatershedElevation.xls in Microsoft Excel.
  2. Delete the OBJECTID, ZONE_CODE, COUNT, AREA, and SUM fields.
  3. Highlight the minimum and maximum values in each column.

Continue formatting the table until you are satisfied with its appearance.

  1. Click the File menu and select Save As.
  2. Navigate to your HydrologyLab folder.
  3. For ‘File name:’, type “WatershedElevation”.
  4. Use the ‘Save as type:’ drop-down menu to select Excel Workbook (*.xlsx).
  5. Click Save and close Excel.

FOR TABLE TO BE TURNED IN

Create a table highlighting the minimum and maximum values for the minimum, maximum, range, mean, and standard deviation of all elevation values within each watershed.

Image Removed

 

 

Part 3: Downloading rainfall data

Downloading data from Climate Data Online.

Now you will download rain gauge station data created by the National Climatic Data Center (NCDC) using the Climate Data Online (CDO) interface.

...

Excel

  1. Navigate to the location where the CSV file was stored.
  2. Double-click the CSV file to open it using Excel.

The first column contains the unique station identification code and the second column contains the station name. Next are the elevation, latitude, and longitude of the stations. The annual precipitation field contains long-term averages of annual precipitation totals in hundredths of inches. More information is available on the Data Set Documentation and Samples portion of the CDO website. The annual precipitation field "ANN-PRCP-NORMAL" may not be in the F column, it might be in the BP column, in which case you will need to move it there using the copy-paste function. 

Before opening this table in ArcGIS, you must reformat some of the field names, which cannot have special characters and must be 13 characters or less.

  1. Rename “ANN-PRCP-NORMAL” to “ANNPRCPHI”, for annual precipitation in hundredths of inches.
  2. Along the top of the worksheet, drag across the column letters to select columns A through H.
  3. Copy the columns. Along the bottom of the worksheet, click the "plus" icon to create a new sheet. In the A1 cell, paste your previously copied columns. 
  4. Delete the previous sheet. 
  5. Along the top of the worksheet, drag across the column letters to select columns A through H.
  6. Hover your mouse between columns G and H until the cursor changes to two outward facing arrows and double-click to auto-size the column widths.
  7. At the bottom left of the worksheet, rename the worksheet “PrecipStations”.
  8. Click the File menu and select Save As.

10. Navigate to your Lab 1 folder.

11. For ‘File name:’, type PrecipStations.

12. Use the ‘Save as type:’ drop-down menu to select Excel 97-2003 Workbook.

13. Click Save.

14. Close Excel.

Displaying XY data

Now you are ready to start a new map document and display the tabular rain gage data you just downloaded.

...

Exporting XY data

Since the points appear to be in reasonable locations (rather than in another country or the middle of the ocean), you will want to export them to a new feature class in your ElevationRainfall geodatabase.  Exporting to a feature class will allow you to reuse this points layer in other future map documents without having to go through the display XY data process each time.

    Right-click the PrecipStations_Layer layer and select Data > Export Features.
    Because you last exported a text file outside of your geodatabase, notice the output feature class is now set to a shapefile in your Lab 2 folder, rather than a feature class within your ElevationRainfall geodatabase.
  1. Next to ‘Output feature class’, click the Browse button.
  2. Navigate to the HydrologyLab geodatabase.
  3. Double-click your HydrologyLab geodatabase to save your feature class in it.
  4. For ‘Name:’, type “PrecipStations_Features” and click Save.
    Image RemovedClick Run.
    Since you are now using a permanent feature class, you may remove your temporary Events layer and the corresponding Excel table.
  5. Right-click the PrecipStations_Layer layer and select Remove.
  6. Right-click the PrecipStations table and select Remove.

Projecting vector data

Repeating the technique you learned earlier in this lab to project vector data, project the PrecipStations_XYTableToPoint layer into the State Plane Texas South Central projection.Save the resulting feature class and name it “PrecipStations_StatePlane”. Remove the original PrecipStations_XYTableToPoint layer from the Contents pane.

...

  1. In the Contents pane, Ctrl-select the PrecipStations_StatePlane and Watersheds_StatePlane layers.
  2. Right-click either selected layer and selectZoom To Layers.
  3. Open Geoprocessing pane by clicking Tools from Analysis tab.
  4. Return to the Geoprocessing pane.
  5. At the top left of the Geoprocessing pane, click the Back button.
  6. In the search box, type "thiessen".
  7. Click the Create Thiessen Polygons toolDouble-click the Analysis toolbox  then the Proximity toolset then the Create Thiessen Polygons tool.

Before populating the variables in this tool, you will change an Environment setting so that the Thiessen polygons are calculated for the entire region that you just zoomed to.

  1. At the top of the Geoprocessing window, click the Environments button tab on the right.
  2. Under the 'Processing Extent' section, for 'Extent', select Current Display Extent.
  3. At the top of the Geoprocessing window, return to the Parameters tab on the leftFor Extent, use the drop-down menu to select Current Display Extent and go back to Parameters.
  4. For ‘Input Features’,  drag in the  select the PrecipStations_StatePlane layer.
  5. For ‘Output Feature Class’, rename the feature class from “PrecipStationsPrecipStations_StatePlane_Cr” to “PrecipThiessen”Cr” to “PrecipThiessen.
  6. Use the For ‘Output Fields’ drop-down menu to select All select All fields.
  7. Ensure your ‘Create Thiessen Polygons’ window appears as shown below and click Run.

    Image AddedImage Removed

You will notice that polygons now fill the entire Map Display indicating which areas are closest to which rain gages.

  1. Open the PrecipThiessen layer attribute table.

Notice that all of the fields that you originally downloaded from CDO are still included, because you selected to output all fields when running the Create Thiessen Polygons tool. If you do not see all of the same fields, re-run the tool and this time output all fields.

  1. Close the Table PrecipThiessen attribute table.

Intersecting two polygon layers

In order to determine which portions of the resulting polygons overlap with which watersheds, you will now perform an intersect operation between the two layers. The result will allow you to calculate weighted averages of the precipitation in each watershed.

  1. At the top left of the Geoprocessing pane, click the Back button.
  2. In the search box, type "intersect".
  3. Click the Intersect toolIn the Analysis Tools toolbox, click the Overlay toolset then the Intersect tool.
  4. For ‘Input Features’, drag in the select the PrecipThiessen and the Watersheds_StatePlane layers.
  5. For ‘Output Feature Class’, rename it from “PrecipThiessenPrecipThiessen_Intersect” Intersect to “ThiessenWatershedIntersect”ThiessenWatershedIntersect.
  6. Ensure your ‘Intersect’ window appears as shown below and click Run.
    Image Removed
    Image Added

  7. Remove In the Contents pane, remove the PrecipThiessen layer from the Contents pane.
  8. Zoom to the ThiessenWatershedIntersect layer.

The resulting layer integrates all of the boundaries from both the Thiessen polygons and the watersheds, limited to the extent of their overlap.

  1. Open the ThiessenWatershedIntersect layer attribute table.

    Notice that the original 8 watersheds have now been divided into 42 41 sections indicating which areas of each watershed are closest to each rain gage. Let Pk denote the annual precipitation associated with each rain gage and Aik denote the area of the intersected polygon associated with rain gage k and watershed i. The area weighted precipitation associated with each watershed is
    Image Removed
    Image Added
    You will add a new field to the table to calculate the elements of the numerator of the equation.|


  2. At the top left of the attribute table, click the Click Add Field… button on top , which will open up the Fields view of the table display.



  3. For ‘Name:’, type “APProd”At the bottom of the Fields table, rename the new field from Field to “APProd.
  4. Use the ‘Type:’, ‘Data Type’ drop-down menu to select Double and click.
  5. On the Ribbon, on the Fields tab, click the Save button.

    Image Added


  6. Close the Fields view of the table.
  7. Scroll to the far right of the ThiessenWatershedIntersect table. Save on the Fields tab.
    Image Removed
  8. Right-click the APProd field name and select Calculate  Calculate Field.
  9. Using the fields and buttons or by typing, enter “!ANNPRCPHI! * !Shape_Area!”.
  10. In the 'Fields' list, double-click ANNPRCP_HI.
  11. Underneath the 'Helpers' list, click the * button.
  12. In the 'Fields' list, double-click Shape_Area.
  13. Ensure your ‘Calculate Field’ window Ensure your ‘Field Calculator’ appears as shown below and click Run OK.

    Image Removed

    Image Added

The APProd field now contains the numerator values in the equation. You are now ready to summarize the calculated statistics by watershed.

  1. Right-click the HU_10_NAME field name and select Summarize…. field name and select Summarize.

  2. For ‘Output Table’, rename the table from ThiessenWatershedInteract_S to “WatershedPrecip”.
  3. For 'Statistics Field(s)'For statistics fields, use the drop-down menu to select the Shape_Area field and select Sum.

  4. Select the APProd field and select Sum.
  5. For ‘Output Table’, rename the table from “ThiessenWatershedInteract_St” to “WatershedPrecip”.
  6.  'Field' and the Sum 'Statistic Type'.

  7. For the second statistics field, select the APProd 'Field' and the Sum 'Statistic Type'.
  8. Ensure the 'Case field' is the HU_10_NAME field, which will summarize the statistics by watershed.
  9. Ensure your 'Summary Statistics" Ensure your ‘Summarize’ window appears as shown below and click Run OK.
    Image Removed
    Image Added

The resulting table gives the numerator and denominator in the equation for each watershed.

  1. Open the WatershedPrecip table.

Repeating the techniques you just learned, add a new field to the WatershedPrecip table called “Precip” of type double. Use the field calculator to evaluate [Sum_APProd]/[Sum_Shape_Area]. The result is the precipitation for each subwatershed. Export the table to Excel and format it. Also include the mean annual precipitation over the entire watershed.

...

Now you will calculate precipitation for each watershed using a different interpolation method.

  1. Close the WatershedPrecip table and the ThiessenWatershedIntersect attribute table.
  2. Turn off the ThiessenWatershedIntersect layer.
  3. Right-click the Watersheds_StatePlane layer and select Zoom to Layer.
  4. Open ToolboxAt the top left of the Geoprocessing pane, click the Back button.
  5. In the search box, type "spline".
  6. Click the Spline (Spatial Analyst Tools toolbox, click the Interpolation toolset then the Spline ) tool.
  7. On At the top of the window, click the Environments.the Geoprocessing window, click the Environments tab on the right.
  8. Under the 'Processing Extent' section, for 'Extent', select For Extent, use the drop-down menu to select Current Display Extent.
  9. Use the ‘Mask’ drop-down menu to select the Watersheds_StatePlane layer.
  10. Under the 'Raster Analysis' section, for ‘Mask', select the Watersheds_StatePlane layer, which will clip the resulting raster.
  11. At the top of the Geoprocessing window, return to the Parameters tab on the left.
  12. Go back to Parameters. For ‘Input point features’, drag in the select the PrecipStations_StatePlane layer.
  13. Use the For ‘Z value field’ drop-down menu to select the ANNPRCPHI ANNPRCP_HI field that contains the values you wish to interpolate.
  14. For ‘Output raster’, rename the exported raster from “Spline_Preci1” to “AnnPrecip” Spline_Preci1 to “AnnPrecip.
  15. Use the For ‘Spline type’ drop-down menu to select TENSIONselect Tension.
  16.  Ensure your ‘Spline’ window Geoprocessing pane appears as shown below, and click Run.
    Image Removed
    Image Added

The result is a solid surface estimating the rainfall at each cell, based on the data collected at each rain gage. Turn the

...

ThiessenWatershedIntersect polygons back on

...

and give them a hollow fill. Also symbolize the Watersheds_StatePlane layer and put it on top of the Thiessen layer. Symbolize the rain gages as you desire.

...

 

 

FOR MAP LAYOUT TO BE TURNED IN

Create an 8.5 x 11 layout showing the Thiessen polygon boundaries, interpolated rainfall layer and rain gage locations.

TOP MAP WAS THE ORIGINAL MAP, BOTTOM IS MY SCREENSHOT WHICH HAS SOMEWHAT DIFFERENT DATA. NOT SURE IF IT MATTERS BUT THOUGHT IT WORTHWHILE TO AT LEAST CHECK

Image Removed

Image Removed

Image Added

Using the techniques you learned earlier in this lab, use the " Zonal Statistics as Table " tool to export a table containing the total annual rainfall for each watershed, as calculated from the AnnPrecip spline interpolation layer and format the table in Excel.

...

Create a table containing the total annual precipitation for each watershed, as calculated using the spline interpolation method, along with the total annual precipitation over the entire subbasin. 

Deliverables

  1. Create an 8.5 x 11 layout showing the DEM clipped to the subbasin with areas less than or equal to 15 feet in elevation highlighted.
  2. Create an 8.5 x 11 layout showing transparent elevation in graduated colors on top of a hillshade, with watershed boundaries and flowlines visible.
  3. Create a table highlighting the minimum and maximum values for the minimum, maximum, range, mean, and standard deviation of all elevation values within each watershed.
  4. Create a table containing the weighted mean annual precipitation for each watershed, as calculated using Thiessen polygons, along with the total mean annual precipitation over the entire subbasin.
  5. Create an 8.5 x 11 layout showing the Thiessen polygon boundaries, interpolated rainfall layer, and rain gage locations.
  6. Create a table containing the total annual precipitation for each watershed, as calculated using the spline interpolation method, along with the total annual precipitation over the entire subbasin.