EXERCISE 8 – 3D ANALYST and VIEWSHEDS
- or -
CAN I SEE MY HOUSE FROM HERE?

Environmental Resources 372:362

Intermediate Environmental Geomatics


Today we’ll be exploring 3D analysis in ArcGIS using 3D Analyst and Spatial Analyst. In addition to this lab and your notes, check out ESRI’s 3D Analyst tutorial to learn more about this extension. It’s available online from ESRI or you can find a copy in the class folder http://downloads.esri.com/support/documentation/ao_/787Using_3D_Analyst_Tutorial.pdf
\\ad-rsc\data\teach\intgeo\ClassWork\3DAnalystTutorial.pdf


I. TINs and 3D Analyst

In ArcMap, turn on the Spatial Analyst and 3D analyst extensions (View | Toolbars). You will see a toolbar appear but if you click on the pull-down menu you will see that nothing is highlighted. You must enable the extensions (Tools | Extensions). When you’re done using the 3D analyst extension, turn it off.

Last week we worked with a raster elevation format, the DEM. Rasters aren't the only way to store elevation data. The triangulated irregular network (TIN) is a way to store elevation data (or other surface data) as a set of triangles derived from irregularly spaced sample points.

In ArcMap, add the DEM of part of Hopewell Township (\\ad-rsc\data\teach\intgeo\classwork\hopewell\penndem_sp). To create a TIN from a raster data set containing elevation, activate the 3D Analyst | Convert | Raster to Tin command. Keep the default Z tolerance and be sure you're saving the TIN to your directory. Once the TIN is created, zoom in to the hilly area in the northwest part. Notice how the surface is composed of triangles. The TIN display defaults to hillshaded in order to give you an idea of the 3D-ness, but it's not really 3D, is it? Don't you feel kind of cheated? Let's fix that.

To get to something approaching real 3D, we need to launch ArcScene, ArcGIS's 3D scene viewer. To do this, click on the ArcScene icon to the right of the 3D Analyst toolbar. The ArcScene viewer looks a lot like the ArcMap interface, so you should feel relatively comfortable with the controls. Copy the TIN in ArcMap and paste the layer into ArcScene. Now that looks a bit more like it, huh? If New Jersey wasn't so flat it would look even better. Luckily we can exaggerate the vertical relief (View | Scene Properties | General). Try a vertical exaggeration of 10.

Let's see what we can do with the controls. The pan and zoom functions work similar to those in ArcMap. Since we're in 3D, we can also rotate the scene using the navigate control. Go slow, you might get dizzy. One of the super nifty neato things we can do is fly-through the scene using the fly-through control, which should be next to the navigate tool in ArcScene’s display tool bar. Once this control is activated, the position of the mouse determines the direction the camera moves, the left mouse button increases the fly-through speed, and the right mouse button decreases the speed.

Now that you've spent some time pretending you're a bird, let's try displaying something more exciting than just the TIN. Add an aerial photo (\\ad-rsc\data\teach\intgeo\classwork\hopewell\771.jpg) to the scene. Now, we need to tell ArcScene how to drape the photo over the TIN. Under the aerial photo’s Layer Properties | Base Heights, check the box that tells ArcScene to use the TIN as the data layer to obtain the base heights from. Turn off the TIN and explore the draped aerial photo a bit. Set the extent of the scene to be equal to the extent of the aerial photo by right-clicking on ‘Scene layers’ in the TOC and selecting Properties | Extent. Now only data overlapping the aerial photo will displayed (theoretically).

We can also drape thematic data over the TIN. Add \\ad-rsc\data\teach\intgeo\classwork\hopewell\pennlulc to the scene. Set the symbology to use the field "level1", the transparency to 60%, and tell ArcScene what layer to use for the base heights. As you can see, ArcScene isn't perfect, but given the right circumstance you may find it useful.

Some other applications of ArcScene you should be aware of but we won't demo:

Viewshed Analysis -

Ok, so flying through the landscape is cool and all, but what about real 3-D analysis? Let's go back to ArcMap. In 3-D Analyst, we can create topographical profiles of lines, draw lines of sight from points, and find the steepest path downhill from a specific point. One of the more widely used 3-D analyses is the creation of viewsheds. A viewshed is simply the land area visible from any given point. You can perform viewshed analysis in either 3-D Analyst or Spatial Analyst.

We need two data layers in order perform a viewshed analysis: a TIN or raster representing elevation and a point or line layer representing the vantage points from which we want to determine the viewshed. We'll use the TIN we created earlier as our elevation data. The vantage points are represented in \\ad-rsc\data\teach\intgeo\classwork\hopewell\viewpoints. They are prospective overlooks on a trail along the ridge in the north of the DEM. Create the viewshed (Surface Analysis | Viewshed under the 3D or Spatial Analyst menu), accepting the default cell size and saving the viewshed to your directory. The output is a raster data layer, automatically symbolized to show what is visible and not visible from the points. Open the attribute table for the viewshed and you'll see it's a little more complicated than that. Cells with a value of 0 are not visible from any of the points, cells with a value of 1 are visible from one point, and those with a value of 2 are visible from all points.

III. Creating Topographical Profiles
3-D Analyst also allows us to examine the profile of data in 2D. In ArcMap, add the Hopwell Township DEM (penndem_sp) that we used at the beginning of lab. This is the only file we need.

First, use the Interpolate Line tool in 3D Analyst to draw a line between a point of high elevation and a point of lower elevation. When you are finished, double-click to stop digitizing. Click the Create Profile Graph tool. When the graph pops up, double click on it to open up the Properties window. Here you can change the title, subtitle, and colors, and add another data series or a trend line. Under the Appearance tab, change the title and add units to the Y-axis label (e.g., “Z (feet)”). Let’s add a trend line, too. Under the Series tab, select New Function from the Add dropdown menu at the bottom of the window. Function type should be Curve Fit. Change the color of the line so it’s not the same color as the profile line. When you’re happy with the way your graph looks, click OK.


Assignment

  1. Make a map using the viewshed results and the landcover found at \\ad-rsc\data\teach\intgeo\classwork\hopewell\pennlulc.shp. Your map should emphasize the Level 1 land use categories visible from either point. Although you can't tell precisely which areas are viewable from which point, you can make a general guess. Explain on your map which point should be the overlook point if the goal is to emphasize the rural nature of the area.
  2. Make a map of your profile graph. The map should show the DEM, the interpolated line (Hint: double-click on the line in Data View to change its color and width so it’s visible against the DEM), and the profile graph.

Assignment due at the beginning of lab on Monday, April 6th.