GEOB 270

Lab 3: Planning for a Tsunami/Storm Surge

Spatial Analysis, Tables, Editing

Introduction

In this lab, you will be working more with spatial and tabular datasets, both raster and vector. You will be reviewing from lab 2 spatial and raster data, coordinates systems, projections, and then be introduced to a variety of basic GIS Analysis tools.

You will be given a scenario for which you will perform GIS analyses. You will identify:

  • the problem to be solved;
  • the data you need;
  • the geographic analysis to be performed;
  • the analytical operation or GIS tools you need to access in ArcGIS.

Scenario

You have been hired by the City of Vancouver GIS department. After geological and climate change weather events that have transpired recently around the world, your supervisor wants you to prepare a map that highlights the areas in the city that are more at risk from a tsunami or a storm surge. (Risk of tsunami in Vancouver is small due to Vancouver Island, however, there is still a risk, and the analysis in this lab could be applied to any coastal city in the world, or to sea level rise due to climate change, or storm surge driving coastal waters inland).

Although the height, speed and inland reach of a tsunami depend on many factors, we will simplify for this exercise that the immediate danger zone from a 10-metre wave is mainly to the low lying areas at, or below, 10 metres that lie up to 1 kilometre from the shoreline.

The municipality is interested in the following information:

  • total area threatened, listed by land use type;
  • length of the roads that might be affected, broken down by road type;

(highway, major road, local road, etc.);

  • number of educational and health care facilities within the danger zone, if any;
  • locations to put up tsunami warning signage

Learning Objectives:

Part 1 -Review Data for Geographic Analysis

Part 2 -Perform Basic Geographic Analysis to determine areas for possible Tsunami

a)perform buffer proximity analysis

b)reclassify raster layers

c)convert raster to vector data files

d)combine vector data layers with polygon overlay tool intersect

Part 3 - Performing Geographic Analysis to Extract Vancouver Data affected by possible Tsunami

a)combine vector data layers with the polygon overlay tool intersect

b)perform a proximity analysis using select by location

c)extract datasets with the polygon overly tool clip

Part 4 - Calculate statistics (areas, length) of Vancouver landuse and roads affected by a potential tsunami

a)create summary tables by area of landuse

b)create lists of facilities affected

c)create summary tables of road infrastructure affected

Part 5 - Add layer of potential signage points

a)Learn how to create a new feature class, explaining the different types (point, multipoint, etc…)

b)Introduce basic editing of features and tables (change values on individual table cells, modification/creation/deletion of features)

c)Introduce the concept of snapping parameters for more accurate positioning of new feature.

Due: 2 week lab, due at the beginning of your lab time in 2 weeks. Any labs submitted after the start of your lab time, will be marked late. Answers to lab questions will be accepted in a .jpg (jpeg) or pdf for map images and pdf for text answers. Submit answers to questions via the course web site:

Part 1 -Review Data for Geographic Analysis

After downloading your data for this lab, launch ArcCatalog as a standalone program and preview the contents of C:\temp\lab3\

The dataset downloaded contains a geodatabase of shapefiles and a grid dataset. Geodatabases provide a file management structure for the geo-spatial data in ArcGIS. With whatever GIS software you are using, it is always important for project management to carefully manage and store your data. Geodatabases offer structural, performance and data management advantages over a collection of shapefiles in a directory. In this lab, we have given you the data in a geodatabase. In the next lab, you will create a geodatabase.

Launch ArcCatalog and view the contents of the geodatabase called lab3.

QUESTION 1(1)

What is the coordinate system, projection, linear units, and datum for your geodatabase?

Close ArcCatalog.

Preparing a Geodatabase for Analysis

Launch ArcMap. In the Getting Started window select the Blank Map from My Templates and set your default geodatabase to C:\temp\lab3\analysis.gdb. Click on OK to begin. Save your map as c:\temp\lab3\mylab3.mxd.

Setting your default geodatabase to C:\temp\lab3\analysis.gdb will help you avoid losing your data. It will insure that the output of any analytical operation is stored in analysis.gdb as a feature class or a raster dataset, unless another location is chosen.

Part 2 -Perform Basic Geographic Analysis to determine areas for possible Tsunami

a)perform bufferproximity analysis

b)reclassify raster layers

c)convert raster to vector data files

d)combine vector data layers with polygon overlay tool intersect

Identifying Metro Vancouver’s danger zone

The first condition of our analysis isto find areas within 1 km of the shoreline. To do this, we will use the Buffer tool.

  • Add the shoreline feature class to your empty map.
  • Open the ArcToolbox window and expand the Analysis Tools by clicking on the ‘+’ sign next to the toolbox.
  • Expand the Proximity toolset and double-click on the Buffer tool.
  • To create a 1 km buffer around the shoreline, enter the following parameters and press OK:
  • Input features: shoreline
  • Output feature class: C:\temp\lab3\analysis.gdb\shore1km
  • Linear unit: 1 kilometers
  • Dissolve type: ALL

Every tool dialog box has a Help section that can be toggled on/off by using the Show Help > or < Hide Help buttons in the bottom-right of the tool’s dialog window. As you enter the different parameters, the help content changes to explain what each parameter means.

  • Click OK.

After a moment, the new shore1km layer is added to the Table of Contents (TOC).

NOTE

Notice that the buffer extends on both sides of the shoreline.

What does 1km represent; each side of the shoreline or the total width (use the measure tool)?

Notice as well that the buffer extends into the mountains, think about why this is, and whether it affect our analysis. This is an example of how you need to be aware of what is happening during GIS analysis and not just blindly click on buttons.

Identifying areas at 10-metre elevation or below

The second condition in our analysis is that the land affected must be at or below 10 metres.

  • Add the Metro Vancouver Digital Elevation Model (dem92g) to your map. Open its attribute table. Notice the number of pixels that exist at or below 10m

QUESTION 2 (1)

What is the resolution of the DEM?

(Hint: remember that units are important when expressing resolution)

QUESTION 3(2)

How many pixels have values of 10 metres or less in the Metro VancouverDEM? How much area in km2 does that represent?

(Hint: use the Statistics command)

For our analysis, we don’t really need the variety of elevation values. What we would like to do is to renumber all the low values as ‘1’ and the higher values to ‘NoData’. To do this simplification, we will use the Reclassify tool:

  • Make sure the Spatial Analyst extension is checked on under Customize >

Extensions…

  • Bring up the ArcToolbox window and select Spatial Analyst Tools > Reclass >Reclassify
  • In the Reclassify dialog window, enter the following parameters and press OK:
  • Input raster: dem92g
  • Reclass field: VALUE
  • Click on the Classify button
  • Under Method:Equal Interval(this will allow us to change thenumber of categories)
  • Classes:2
  • Under Break Values: type10for the first value and click OK
  • You should now have only 3 rows in the reclassification table. For the second row, change the New Value from ‘2’ to ‘NoData’
  • Set the output raster to C:\temp\lab3\analysis.gdb\lowlands. Click OK.

The lowlands raster layer is now added to the TOC. It contains values of ‘1’ for the areas under threat and nothing should be visible for the other areas.

Combining the 1km buffer & the<10 m areas using Intersect

The areas we are searching for are those that are within a) the 1 km buffer and b) below 10 m in elevation. We need to combine these layers to isolate the areas common to both (i.e. the areas susceptible to the tsunami).

You can easily identify their intersection cartographically. To do this:

  • Turn off all layers except shore1km and lowlands. Make sure the buffer is above the lowlands raster.
  • Open the properties for shore1km and under the Display tab enter a Transparency value of 30% and click OK.

The common areas between both layers should be easy to spot. However, in order to use this information in our analysis, we need to create a data layer containing these intersections. We first need to convert the lowlandsraster into vector:

  • In ArcToolbox, select Conversion ToolsFrom RasterRaster to polygon. Enter the following parameters and click OK. The new vector layer is added to the TOC:
  • Input raster: lowlands
  • Field: Value
  • Output polygon features: C:\temp\lab3\analysis.gdb\lowlands_vector
  • Ensure that Simplify polygons is NOT checked

If you have not done so yet, this would be a good time to save your map document. Make sure to save it to C:\temp\lab3\ and NOT inside the analysis.gdb. The geodatabase should only contain actual data files.

To create a vector layer that satisfies both conditions (area susceptible to a tsunami: within 1km buffer, and <10 m elevation) we will use the overlay tool called Intersect:

  • In ArcToolbox, select Analysis toolsOverlayIntersect
  • Enter the following parameters and click OK:
  • Features: lowlands_vector,shore1km
  • Output Feature Class: C:\temp\lab3\analysis.gdb\Danger
  • Accept all other default values
  • After the new layer has been added to the TOC, turn off the lowlands_vector

shore1km layer

  • Zoom in to Danger.

Part 3 - Performing Geographic Analysis to Extract Vancouver Data affected by possible Tsunami

a)combine vector data layers with the polygon overlay tool intersect

b)perform a proximity analysis usingselect by location

c)extract datasets with the polygon overly tool clip

Extract City of Vancouver data in area covered by a possible tsunami

Spatial datasets can be very large. They normally contain more data than are needed for the study at hand. Carrying around extra data can slow down performance. It is a good practice to trim datasets to the minimum extent necessary, although you should always keep a backup copy of your original data on a USB or DVD somewhere just in case. In our case, our interest lies in the City of Vancouver proper and not in the Metro Vancouver area, formerly known as GVRD. We will use a series of overlay tools to narrow down our focus.

  • Add the Vancouvermask layer to your map and zoom to its extent.
  • Using the Intersect tool, overlay the Vancouvermask with Danger and store the new layer as C:\temp\lab3\analysis.gdb\Vancouver_Danger
  • Once the new layer appears in the TOC, turn off all other layers, and add the following layers:
  • landuse
  • roads
  • education(educational institutions stored within feature dataset)
  • healthcare(health careinstitutions stored within feature dataset)
  • Zoom to the full extent.

Notice that the new data layers encompass the whole of Metro Vancouver. We will extract a Vancouver subset using a couple of different tools.

Extract Educational & Health Care institutions using the proximity analysis toolSelect by Location

One of the advantages of GIS is that you can select features in one layer based on their location relative to features in a different layer. Since we have a polygon or mask for Vancouver and the educational & health care layers are points, the Select by Location tool is ideal for our purpose.

  • In the ArcMap menu bar, click on SelectionSelect by Location
  • Input the following parameters before pressing OK
  • Selection method: select features from
  • Under Target layers put a check next to educationhealthcare
  • Source layer: Vancouvermask
  • Spatial selection method for target layer feature(s): intersect the sourcelayer feature

Many points should be selected from those two layers. We will create an independent layer for each that we can use in our analysis.

  • Right-click on education> Data > Export Data
  • Export Selected features and save the new layer as:
  • C:\temp\lab3\analysis.gdb\Vancouver_education
  • Repeat with Healthcare and name Vancouver_health
  • Once the new layers have been added to the TOC,
  • remove the original Metro Vancouver education and health care layers
  • clear the selection
  • save your map document

Extract Land use and Road layers using Overlay tools

It is sometimes necessary to select features in one layer that do not respect the same boundaries in other layers. For example, if you use Select by Location to extract the roads for Vancouver, you may end up with lines extending outside of the Vancouver mask. If you try this, you’ll see that roads in Burnaby and North Vancouver would also be selected. (If you do try this, make sure to clear all selection before performing the steps below).

If we want to impose a boundary on a second layer, we can use the Clip overlay tool:

  • In Arctoolbox, click on Analysis toolsExtractClip
  • To extract roads use the following parameters:
  • Input features: roads
  • Clip features: Vancouvermask
  • Output feature class: C:\temp\lab3\analysis.gdb\Vancouver_roads
  • Leave the other values at their default settings and click Ok

When using the Clip tool, pay special attention to the order in which you input the layers, because it affects the output. Look at the help in the tool for a graphic explanation.

  • Repeat the clipping with landuse to extract the land use for Vancouver, naming the resulting layer Vancouver_landuse.Always be careful where you are saving your files.
  • Change the Symbology of the Vancouver_landuse to represent the ‘Category’ value to better understand the landscape.
  • Remove the landuse and roads layers and zoom into the Vancouver_landuse extent.

Part 4 - Calculating statistics (areas, length) of Vancouver landuse and roads affected by a potential tsunami

a)create summary tables by area of landuse

b)create lists lists of facilities affected

c)create summary tables of road infrastructure affected

Now that we have delineated the danger zone in Vancouver and also extracted the necessary data for the city, we can continue using the overlay tools to determine the land use areas, schools and health care facilities and the parts of the road network that might be affected by a tsunami.

  • Using the appropriate overlay tool, create a layer with the different landuse types that might be affected by the tsunami.
  • Look at the table for the resultant landuse layer affected by the tsunami. Create a Summary Table of the total areas for each land category

(Hint: review Summary Tables in ArcGIS help to get a tabulated list of areas: You are summarizing the field CATEGORY by AREA – SUM)

QUESTION 4 (4)

What are the total areas for each land use category? What are the units of these areas?

What are the total areas for each land use category in Hectares?

Explain the procedure you used in ArcMap to get your answers.

QUESTION 5(1)

What percentage of the City of Vancouver’s total area is in danger?Explain what method you used.

QUESTION 6 (2)

By using the appropriate selection method, create a layer containing the health care and educational facilities within the City of Vancouver danger zone.

List the resulting facilities, if any. Check your results visibly.

Explain how you came up with your answer.

QUESTION 7 (4)

By using the appropriate overlay tool, create a data layer with the roadinfrastructure that might be affected by the tsunami. What are the lengths of damage, to the closest metre, for each road category.

(Hint: the field that describes the type of road is name ‘Road_Type’. At the end of this lab is a file which explains the Road_type code.

Part 5 - Add layer of potential signage points

a)Learn how to create a new feature class, explaining the different types (point, multipoint, etc…)

b)Introduce basic editing of features and tables (change values on individual table cells, modification/creation/deletion of features)

c)Introduce the concept of snapping parameters for more accurate positioning of new features

Now that the danger zone has been outlined, you need to create a layer with the possible locations for appropriate signage warning of the susceptibility of the area to tsunamis.

Creating a new feature class

The new data layers produced so far have been the result of analytical tools. It is sometimes necessary to create an empty canvas to locate new features.