Due at 8am on Monday, 12 February 2018
This assignment is worth 10 points: 5 points for your map and 5 points for your data, awarded according to the evaluation rubric. Your map must present your results accurately and clearly, and adhere to cartographic design principles. Your data must use correct units, data types, and file format.
In this homework, we'll analyze some potential impacts of extreme (as in, maybe Greenland melts...) sea-level rise in Santa Barbara County. To reduce computational requirements, we've restricted our region of interest (ROI) to the downtown Santa Barbara area, rather than the entire county's coastline.
Your assignment is to estimate the impact of ten scenarios of sea level rise: 1m, 2m, . . . , 10m increases in mean sea level. You will compute three metrics for each sea level rise scenario, as shown in Table 1.
sea level rise
We projected all available data into the California Teale Albers (NAD 1983 datum) projection, and they're available in the HW2 geodatabase:
Work through Law & Collins Chapters 15 (Querying data) and 16 (Selecting features by location) on attribute-based and location-based queries.
HW2.7zand extract the contents into your
HW2.mxdand add all the data from your HW2 geodatabase. Organize your map by creating a New Group Layer. Rename it to "Basemap", then move Streets, Cities, County, and California into it.
ParcelMetadata.pdf. For the parcel data, there are several property value attributes. We will use
NET_AV: the net assessed value (in USD) for the parcel.
For the elevation data, the only attribute is
and its value represents the sea level rise (m) required to
flood the land at that elevation. The data also uses an
interval range; that is, the value marks the upper end of a
range. For example, interpret
GRIDCODE = 4 as 3
GRIDCODE ≤ 4m.
The Join Data dialog (Figure 1) is invoked by
right-clicking on the Elevation layer and selecting Joins
and Relates → Join.... Choose the option to "Join data
from another layer based on spatial location". Select Parcels
as the layer to join to the active layer (Elevation),
select Sum as your aggregate function, and name the
spatial join output as
Figure 1: Join Data dialog
Join_Output.GRIDCODEas the join keys.
Figure 2: Field Calculator dialog
Click the "Python" radiobutton and the "Show
Codeblock" checkbox. Type in (or just cut-n-paste) the
following Python code. For each calculation, substitute
for whichever field you're using.
Pre-Logic Script Code:
total = 0 def accumulate(increment): global total if total: total += increment else: total = increment return totalExpression:
Join_Output.Shape_Areafield, which is in m2.
Join_Output.Sum_NET_AVfield, which is in USD.
Design a map that communicates your results. Consider these questions:
HW2.png(as in HW1, use 300 dpi and 8-bit color palette options).
HW2.dbfusing the dBase (DBF) file format. Make sure that your "Field Map" has only 4 fields (S, A, P, and L).