Location of north pole in rotated map
pole = org2pol(origin)
pole = org2pol(origin,units)
pole = org2pol(origin)
returns the location of the North Pole in terms of the coordinate system after transformation based on the inputorigin
. Theorigin
is a three-element vector of the form[latitude longitude orientation]
, wherelatitude
andlongitude
are the coordinates that the new center (origin) had in the untransformed system, andorientation
is the azimuth of the true North Pole from the new origin point in the transformed system. The outputpole
is a three-element vector of the form[latitude longitude meridian]
, which gives the latitude and longitude point in terms of the original untransformed system of the new location of the true North Pole. The meridian is the longitude from the original system upon which the new system is centered.
pole = org2pol(origin,units)
allows the specification of the angular units of theorigin
vector, whereunits
is any valid angle unit. The default is'degrees'
.
When developing transverse or oblique projections, transformed coordinate systems are required. One way to define these systems is to establish the point at which, in terms of the original (untransformed) system, the (transformed) true North Pole will lie.
Perhaps you want to make (30ºN,0º) the new origin. Where does the North Pole end up in terms of the original coordinate system?
pole = org2pol([30 0 0]) pole = 60.0000 0 0
This makes sense: pull a point 30º down to the origin, and the North Pole is pulled down 30º. A little less obvious example is the following:
pole = org2pol([5 40 30]) pole = 59.6245 80.0750 40.0000