Click here to see the SAS code used to draw the map.
Click here to see the example.
---------------------------------------------------------------
We (SAS) don't have a world map with the prime meridian on the right-hand-side
(and Europe/Asia on the left), but that's the layout that lends itself best to
plotting certain data ...
I start with maps.world, and rename long & lat to x & y (I need to work with
unprojected values to do what I'm doing).
I use "proc gproject" with "project=none" to chop the map in two (and close the
resulting polygons), so that I have a left & right side.
I then add an offset to the x values of the right side of the map,
so they wrap-around and are on the left side), and also added 5000 to
the segment so that the countries that were "split" would have
new/unique segment numbers.
And I then put the 2 halves together into 1 data set again.
When I "proc gproject" the map using the cylindri projection,
I have to use "norangecheck" because otherwise it would think the coordinates were
out of range, because of the offset I used to put the right 1/2 on the left.
Back to Samples Index