01-18-2017 04:57 PM
Statistics Canada has released the files for the new census boundaries and I'm trying to create a map for Alberta.
The data they provide is for all of Canada, so if I map all of Canada the map appears correct, but if I map only Alberta it's rotated. I'd like to rotate it a centered version, more similar to what you'd see in Google Maps - with the top boundary parallel to the top of the page.
Is there a straightforward option to do this? GPROJECT? My Google/Search is failing me today
I've attached a data step version to create the file required for mapping, in a text file (zipped because it was too big otherwise) and the code to create the map is below.
Other features that would be nice, that I'll work on next:
proc gmap map=alberta_CSD2016; id csdname ; choro csdname; run; quit;
01-19-2017 12:58 PM
Looking at the data in your map data set, it appears to be projected.
It is a similar projection to our MAPSGFK.CANADA data set. If you do all of Canada it is on the same angle as when you display Alberta only. You need to find what projection it is in to "unproject" it with Proc GPROJECT and the PROJ4 options (From= and To=).
Here is the code:
data alberta; set mapsgfk.canada; /*Uncomment for UNPROJECTED*/ /*x=long; y=lat;*/ if id1="CA-48" then output; run; proc gmap map=alberta data=alberta; id id; choro id; run;
Here is MAPSGFK.CANADA projected and unprojected: