Census block maps in SAS

0

I've heard many customers ask "Can SAS do/use/handle census block maps?"

My answer was "of course!" ... but I never really had any examples to show/prove that.  Well, now I do! :)

First, you need to locate & download the shape files for the census blocks you're interested in plotting.  For the 2010 U.S. census, here is a link to one location where you can download the files from the census web site.  Each state has a zip file, with the FIPS code for that state in the name.  For example, North Carolina is FIPS code 37, and the North Carolina zip file is named 'tabblock2010_37_pophu.zip'.

Once you download and extract the zip file, you can import the shape file into a SAS data set using Proc Mapimport, such as ...

proc mapimport  datafile="tabblock2010_37_pophu.shp"  out=tabblock2010_37_pophu;  id blockid10;  run;

The output SAS data set (tabblock2010_37_pophu) contains both the x/y coordinates for the borders of the census blocks, and also the 2010 population of each block (stored in the variable named pop10).

Census blocks are very small/granular areas, therefore if you try to plot the entire state and shade each area, your computer screen just won't really have enough pixels for you to discern each census block.  I would recommend subsetting to show just one county.  For example, here is a map of Wake County, NC (where the SAS headquarters is located):

Wake county census block map, shaded by population

 

But what if you *do* want to show the entire state population at the block level???

The approach I took was to calculate the center coordinate of each census block (using the %centroid macro), and then plot a single annotated 'point' at that center, and make the color of the point represent how many people live in that census block.  You could set the color variable manually in a data step using if/else statements, but in this example I decided to try out a new (upcoming) feature in SAS 9.4 where I stack/overlay alpha-transparent gray points in the same location (the higher the population, the more points I stack, and the darker the point will appear).  I think you will like the results!

Be sure to click each of the maps above, to see them full-size, and the full-size map will allow you to hover your mouse over the cities and see the city names!

And if you'd like to see the exact SAS code used to create the above examples, here are links for the choro map, and the dot map .  This is a great way to apply the skills you learned in the SAS Mapping Class!  :)

 

Share

About Author

Robert Allison

The Graph Guy!

Robert has worked at SAS for over 20 years, and is perhaps the foremost expert in creating custom graphs using SAS/GRAPH. His educational background is in Computer Science, and he holds a BS, MS, and PhD from NC State University. He is the author of several conference papers, has won a few graphic competitions, and has written a book (SAS/GRAPH: Beyond the Basics).

Leave A Reply

Back to Top