This Python 3.3 module implements a wrapper for map projections.
CHANGELOG:
NOTE:
All lengths are measured in meters and all angles are measured in radians unless indicated otherwise. By ‘ellipsoid’ below, i mean an oblate ellipsoid of revolution.
Bases: builtins.object
Represents a map projection of a given ellipsoid.
INSTANCE ATTRIBUTES:
EXAMPLES:
>>> from rhealpix_dggs.ellipsoids import WGS84_ELLIPSOID
>>> f = Proj(ellipsoid=WGS84_ELLIPSOID, proj='rhealpix', north_square=1, south_square=0)
>>> print(my_round(f(0, 30), 15))
(0.0, 3748655.1150495014)
>>> f = Proj(ellipsoid=WGS84_ELLIPSOID, proj='cea')
>>> print(my_round(f(0, 30), 15))
(0.0, 3180183.485774971)
NOTES:
When accessing a homemade map projection assume that it can be called via a function g(a, e), where a is the major radius of the ellipsoid to be projected and e is its eccentricity. The output of g should be a function object of the form f(u, v, radians=False, inverse=False). For example, see the healpix() function in pj_healpix.py.