Click here to see the SAS code.
Click here to see the example.
-----
See my blog for more information about this graph!
Here's how I got the elevation data ...
http://www.ngdc.noaa.gov/mgg/gdas/gd_designagrid.html
Your Grid Id: bobsas2
Grid Database: ETOPO2 - 2-minute global relief
x - Windows
Lat/Lon (left defaults - lat(up/down) = 33N-38N , lon (across)=85W -> 73W
Cell size - 2-minute
Cell params - (4-byte floating point, tenths of meters precision)
Format: XYZ (lon,dat,depth), No Header, Space delimited, omit empty grid cells
Clicked 'Design-A-Grid' button.
Clicked 'Next-->' button
Clicked 'Compress and Retrieve' button
Clicked 'Retrieve' button
Clicked 'save'
Saved bobsas2.zip to my U:
Then, after saving it, I unzipped it, which created a directory structure...
There's a lot of 'junk', but the file I wanted was ...
/u/realliso/bobsas2_data/bobsas2/bobsas2.xyz
(I renamed that to rise.xyz in the current directory)
-----
For each of the elevation long/lat coordinates, I created an
annotate dot (pie), with the elevation determining what the
color/shade should be.
I also created annotated text labels & dots for some city names,
for points-of-reference.
I combined the annotate data sets with a US county map,
and 'gproject'd them, clipping along the boundaries for
the area of interest.
I then separated the map from the annotate, and then used
gmap to draw the map outline (using an 'empty' fill pattern),
and then annotating the elevation color behind the map,
and the city labels on top of the map.
I set a macro variable with the amount of sea level rise,
and then change the color for the elevations between 0 and
the rise level.
Back to Samples Index