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;
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:
Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.
If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.