Main Content


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 input origin. The origin is a three-element vector of the form [latitude longitude orientation], where latitude and longitude are the coordinates that the new center (origin) had in the untransformed system, and orientation is the azimuth of the true North Pole from the new origin point in the transformed system. The output pole 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 the origin vector, where units 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

See Also