Wednesday, 7 October 2015

coordinate system - correct projection (PROJ.4 string) to rotate map


I want to project a map as being rotated by x degrees (as opposed to the standard "north is up" orientation).



So far I'm able to get some arbitrary rotation using these projections:


Two Point Equidistant as used in this Gist


+proj=tpeqd +lat_1=35 +lat_2=35 +lon_1=-80 +lon_2=-122  +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs

Azimuthal as described here: Rotating ~90° using Two-Point Equidistant projection with Proj4


+proj=aeqd +ellps=sphere +lat_0=90 +lon_0=-90

However, I'm having trouble actually converting a value in degrees (like 45°) to a corresponding PROJ.4 string for a specific map. Can anyone point me in the right direction?




No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...