Lab 1: Digital Terrain Data Models

Due by October 8 before class.

Introduction

 

This lab has three separate tasks that will let you explore different digital terrain data models. You will download a DEM from an Internet website. Task 2 involves the use of points and a DEM to learn how ArcGIS handles 3D objects. In Task 3, you will convert a DEM to a TIN data model. The exercise will be done using ArcGIS. This is your opportunity to explore the software and to get acquainted with the interface. Use the help system to learn about the different software options and functionality. Spend some time browsing the websites that you visit in this lab. Type the answers to all questions at the end of each task, attach your map, and submit them to the instructor by the due date.

 

Before you start the first lab, you need to know where to save your lab data. You are only allowed to put data in the C:\Users folder on the computers in CH469 or CH475. You can use your odin ID or your name to create a working folder in C:\Users so that your data won’t get mixed with other students’. If you want to access your working folder from other computers on campus, please follow the instructions on this page. Then create a folder for lab1.

 

 

Task 1. Download and Process DEM and DLG from the Internet

 

Instructions

Task 1 guides you through the process of downloading a DEM and converting it to a grid. The DEM is in SDTS format, and the SDTS files are in the TAR.GZ zipped format.

 

  1. Go to the OIS Data Depot website http://data.geocomm.com/. (OIS Data Depot is one of three websites for downloading USGS data.) Please note that there are many other DEM online sources, for example, http://ned.usgs.gov/.

 

  1. Use the following links to download the Portland, OR DEM 24K: USGS DEMs > Download DEM Data Here> Click Oregon on the map> Select Washington > Click DEM 24K> Select the bluish download button of Portland, OR > Set up an account with GeoCommunity (make sure you provide a valid email account for verification) > Click 1646292.DEM.SDTS.TAR.GZ (l0 meter) to download > Save the file to your lab1 folder. This DEM has a 10-meter resolution. Also save the 1646292.DEM.SDTS.TXT file in your lab1 folder.

 

  1. Now we need to unzip (decompress) the gz file. Double-click on the gz file and check “Select the program from a list”. Click OK and then click on the “Browse” button. Navigate to “C:\Program Files\7-zip” folder and select “7zFM.exe”. The 7-Zip File Manager window will appear. Click on the Extract icon, specify your output folder (e.g., C:\Users\Myfolder\Lab1\) and click OK. This will unzip the gz file to your lab1 folder. (Note: You can unzip gz files using WinZip as well, but we don’t have WinZip installed on our lab computers.) You will find a new file with the .TAR file extension. Now, create a folder call elevation under lab1 folder and move the .TAR file to the elevation folder and repeat the procedures above to unzip the data (Use the “edit/select all” option to extract all files). When done, there are several SDTS data descriptive files (DDF) in the folder.

 

  1. Now you will download a DLG and a master data dictionary for the same area. Go to the USGS Geographic Data Download website (http://edc.usgs.gov/geodata/) (Note: USGS has a new web portals for GIS data. The link is http://eros.usgs.gov.) Click on the 1:24K DLG STDS Format, and select State. Click Oregon, and then Portland. Click on hydrography and then 903780.HY.sdts.tar.gz. Save the file in your folder. Go back to the USGS Geographic Data Download web page. Click on 1:24K as a single tar file under SDTS DLGs require Master Data Dictionary. Save the file (00MASTERDD_LRG.SDTS.tar.gz) in your folder.

 

  1. At this point, you should have 903780.HY.sdts.tar.gz and 00MASTERDD_LRG.SDTS.tar.gz, both zipped files, in your lab1 folder. Create a new folder named hydrography and another new folder named masterdd. Move 903780.HY.sdts.tar.gz to the hydrography folder and 00MASTERDD_LRG.SDTS.tar.gz to the masterdd folder. Now you can unzip the two tar files.

 

  1. Start ArcCatalog. Click Show/Hide ArcToolbox Windows to open the ArcToolbox window in ArcCatalog. Double-click the Import From SDTS tool in the Coverage Tools/Conversion/To Coverage toolset. In the Import From SDTS dialog, use the browse button to navigate to the data files in the elevation folder. The data files all have the 7974 prefix. Double-click any of the files. The Input SDTS Transfer File Prefix in the dialog should list 7974. Change the output grid name to pdxdem and save it in your lab1 folder (not the elevation folder). Click OK to execute the conversion.

 

  1. Double-click the Import From SDTS tool in the Coverage Tools/Conversion/To Coverage toolset. Use the browse button to navigate to the data files in the hydrography folder. The data files all have the HY01 prefix. Double-click any of the files. Save the coverage as pdxhydro in your lab1 folder (not the elevation folder).

 

  1. Launch ArcMap. Add pdxdem and the arc layer of pdxhydro to the data frame. Change the color of pdxhydro.arc to blue. The two data sets should register spatially.

 

Questions

  1. What is the unit of the x and y coordinates in pdxdem and what is the unit of elevation values in pdxdem? Where did you find the information?
  2. What is the elevation range (in meters) in pdxdem?
  3. What is the coordinate system that both pdxdem and pdxhydro are based on?
  4. What other layers besides arc are contained in the pdxhydro coverage?

 

 

Task 2. Extracting Elevation Data to a Point Layer

 

Instructions

You will use the DEM you just created to assign elevation data to a RLIS fire station layer.

 

1.      Add I:\Students\data\GIS\RLIS\2008_Aug\ESRISHAPEFILES\PLACES\fire_sta.lyr to the same data frame you created in Task1. Click the selection tab of the left panel of ArcMap, check the fire stations layer and uncheck all other layers. Use the Select Features tool to select only the fire stations within the extent of the DEM.

 

2.      Right-click on the fire_sta layer and select Data/Export Data. Check the data frame option button so that the output uses the same coordinate system as that of the data frame (i.e., DEM’s coordinate system. Type in the output shapefile’s path and name (e.g., C:\users\odin\firestations.shp) and click OK to finish the feature extraction process. Make sure you have read/write access privilege to the path you specified in the output field. When done, open the attribute table of firestations and make sure it has 13 records. (Alternative method: This method requires the use of “project” tool to reproject the fire station layer to the projection of the DEM. Click Show/Hide ArcToolbox Windows to open the ArcToolbox window in ArcMap. Double-click the Select tool in the Analysis Tools/Extract toolset. Click on the down arrow of the input features field and select fire_sta from the dropdown list. Set output feature class as firestation shapefile in your lab1 folder. Click OK to continue. When using the PROJECT tool in Arctoolbox, you need to specify two transformation methods: NAD_1983_To_HARN_WA_OR and then NAD_1927_To_NAD_1983_6.)

 

3.      Start ArcScene. Add firestations and pdxdem to the data frame in ArcScene. (Do you see a flat map of DEM and points that you can spin around when you drag the cursor?) Double-click on pdxdem layer to open the layer properties window. Click on the Base Heights tab and click Obtain heights for layer from surface. You should see pdxdem listed in the field of that option. Click on the Symbology tab and set Histogram Equalize as the stretch method. Click OK to continue. Double-click on the firestations layer to open the layer properties window. Click on the Base Heights tab and click Obtain heights for layer from surface. Use pdxdem as the surface layer. Click OK. Now you should see a 3D rendering of the DEM and the firestations points.

 

4.      Open the ArcToolBox window in ArcScene and make sure that the Spatial Analyst extension is enabled. Double-click the Extract Values to Points tool in the Spatial Analyst Tools/Extraction toolset. Select firestations as the input point feature and pdxdem as the input raster. Set output as firest_z shapefile in your lab1 folder. Click OK to continue. Open and compare the attribute tables of firestations and firest_z. Answer Question 1 below now.

 

5.      Go to View/Toolbars on the ArcScene pulldown menu and open the 3D Analyst toolbar. Select “Features to 3D” from the 3D Analyst/Convert menu. Set firest_z as the input feature and use “RASTERVALU” as the input feature attribute option. Save the output as firest_3d in your lab1 folder.

 

6.      The firest_3d layer is automatically added to the data frame. Is it displayed as a flat map of points? Double-click on the firest_3d layer. Click on the Base Heights tab and check its setting. Open and compare the attribute tables of firestations and firest_3d. Answer Question 2 below now.

 

7.      Double-click the Add XY coordinates tool in the Data Management Tools/Features toolset in ArcToolBox. Set firestations as input features and click OK. Use the same tool to process the firest_3d layer. Open and compare the attribute tables of firestations and firest_3d. Answer Question 3 below now.

 

8.      (Optional step) Repeat step 5 and use the raster or TIN surface option (specify pdxdem as the raster) to convert firestations to a new 3D shapefile. What happened? Explain.

 

Questions

  1. What’s the difference between the attributes of firestations and firest_z?
  2. What’s the difference between the attributes of firestations and firest_3d?
  3. The firestations layer is a 2D featureclass and firest_3d is a 3D featureclass. How can you tell if a featureclass is 3D (or 2.5D) in ArcGIS?
  4. How many fire stations have an elevation higher than 50 meters? What’s the highest elevation and lowest elevation (in meters) among the 13 fire stations?

 

 

Task 3. Build and Display a TIN

 

Instructions

You will use the pdxdem to construct a TIN. You will also display different features of the TIN.

 

1.      Double-click the Raster to TIN tool in the 3D Analyst Tools/Conversion toolset. Select pdxdem for the input raster, specify pdxtin for the output TIN, and change the Z Tolerance value to 10. Click OK to run the command.

 

2.      You can view pdxtin in a variety of ways in ArcMap or ArcScene. In ArcMap, select Properties from the context menu (i.e., right-click menu) of pdxtin. Click the Symbology tab. Click the Add button below the Show frame. An Add Renderer scroll list appears with choices related to the display of edges, faces, or nodes that make up pdxtin. Click "Faces with the same symbol" in the list, click Add, and click Dismiss. Uncheck all the boxes in the Show frame except Faces. Make sure that the box to show hillshade illumination effect in 2-D display is checked. Click OK on the Layer Properties dialog. With its faces in the same symbol, pdxtin can be used as a background in the same way as a shaded relief map for displaying map features such as streams, vegetation, and so on.

 

Questions

1.      The default Z tolerance in the Convert Raster to TIN dialog is 125.4. What happens when you change the tolerance to 10?

2.      How many nodes and triangles are in pdxtin?

3.      Make a map showing the TINs you created with different Z tolerance values.