org2pol

Location of north pole in rotated map

Syntax

pole = org2pol(origin)
pole = org2pol(origin,units)

Description

pole = org2pol(origin)returns the location of the North Pole in terms of the coordinate system after transformation based on the inputorigin. Theoriginis a three-element vector of the form[latitude longitude orientation], wherelatitudeandlongitudeare the coordinates that the new center (origin) had in the untransformed system, andorientationis the azimuth of the true North Pole from the new origin point in the transformed system. The outputpoleis 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 theoriginvector, whereunitsis 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.

Examples

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

|