pymgp.crstrans.ellcurvature#

pymgp.crstrans.ellcurvature(lat, height, unit='rad', ellipsoid='current')[source]#

Multiplication factors to convert latitude and longitude angles to meters.

Parameters:

latarray_like or float: Latitude in radians or degrees.
heightarray_like or float: Height in meters.
unit{‘rad’, ‘deg’}, default = ‘rad’: Unit of latitude (input) and for multiplication factors (output).
ellipsoidstr or list of floats, default = ‘current’: Text string with the name of the ellipse or a list [a, 1/f] with the semi-major axis a and inverse flattening 1/f. The current ellipsoid can be set by setell. Default for the current ellips is ‘WGS-84’.

Returns:

flat, flonfloat: Conversion factors for latitude and longitude in unit.

Examples

>>> flat, flon = ellcurvature(52, 0, unit='deg')
>>> print(flat, flon)
111267.44260909088 68677.7753788433

>>> print(f'In Delft one arcsec is {flat/3600:.3f} m in latitude and {flon/3600 :.3f} m in longitude.')
In Delft one arcsec is 30.908 m in latitude and 19.077 m in longitude.