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