def toUtmUps8(latlon,
lon=None,
datum=None,
falsed=True,
Utm=Utm,
Ups=Ups,
pole='',
name='',
**cmoff):
'''Convert a lat-/longitude point to a UTM or UPS coordinate.
@arg latlon: Latitude (C{degrees}) or an (ellipsoidal)
geodetic C{LatLon} point.
@kwarg lon: Optional longitude (C{degrees}) or C{None}.
@kwarg datum: Optional datum to use this UTM coordinate,
overriding B{C{latlon}}'s datum (C{Datum}).
@kwarg falsed: False both easting and northing (C{bool}).
@kwarg Utm: Optional class to return the UTM coordinate (L{Utm})
or C{None}.
@kwarg Ups: Optional class to return the UPS coordinate (L{Ups})
or C{None}.
@kwarg pole: Optional top/center of UPS (stereographic)
projection (C{str}, C{'N[orth]'} or C{'S[outh]'}).
@kwarg name: Optional name (C{str}).
@kwarg cmoff: DEPRECATED, use B{C{falsed}}. Offset longitude
from zone's central meridian, for UTM only (C{bool}).
@return: The UTM or UPS coordinate (B{C{Utm}} respectively B{C{Ups}})
or a L{UtmUps8Tuple}C{(zone, hemipole, easting, northing,
band, datum, convergence, scale)} if B{C{Utm}} respectively
B{C{Ups}} is C{None} or B{C{cmoff}} is C{False}.
@raise RangeError: If B{C{lat}} outside the valid UTM or UPS bands
or if B{C{lat}} or B{C{lon}} outside the valid
range and L{rangerrors} set to C{True}.
@raise TypeError: If B{C{latlon}} is not ellipsoidal or B{C{lon}}
value is missing.
@raise UTMUPSError: UTM or UPS validation failed.
@raise ValueError: Invalid B{C{lat}} or B{C{lon}}.
@see: Functions L{toUtm8} and L{toUps8}.
'''
lat, lon, d, name = _to4lldn(latlon, lon, datum, name)
z, B, p, lat, lon = utmupsZoneBand5(lat, lon)
f = falsed and cmoff.get('cmoff', True)
if z == _UPS_ZONE:
u = toUps8(lat,
lon,
datum=d,
falsed=f,
Ups=Ups,
pole=pole or p,
name=name)
else:
u = toUtm8(lat, lon, datum=d, falsed=f, Utm=Utm, name=name)
return u
def _to7zBlldfn(latlon, lon, datum, falsed, name, zone, Error, **cmoff):
'''(INTERNAL) Determine 7-tuple (zone, band, lat, lon, datum,
falsed, name) for L{toEtm8} and L{toUtm8}.
'''
f = falsed and _xkwds_get(cmoff, cmoff=True) # DEPRECATED
lat, lon, d, name = _to4lldn(latlon, lon, datum, name)
z, B, lat, lon = _to3zBll(lat, lon, cmoff=f)
if zone: # re-zone for ETM/UTM
r, _, _ = _to3zBhp(zone, B)
if r != z:
if not _UTM_ZONE_MIN <= r <= _UTM_ZONE_MAX:
raise Error(zone=zone)
if f: # re-offset from central meridian
lon += _cmlon(z) - _cmlon(r)
z = r
return z, B, lat, lon, d, f, name
def toUps8(latlon, lon=None, datum=None, Ups=Ups, pole='',
falsed=True, strict=True, name=''):
'''Convert a lat-/longitude point to a UPS coordinate.
@param latlon: Latitude (C{degrees}) or an (ellipsoidal)
geodetic C{LatLon} point.
@keyword lon: Optional longitude (C{degrees}) or C{None}
if B{C{latlon}} is a C{LatLon}.
@keyword datum: Optional datum for this UPS coordinate,
overriding B{C{latlon}}'s datum (C{Datum}).
@keyword Ups: Optional (sub-)class to return the UPS
coordinate (L{Ups}) or C{None}.
@keyword pole: Optional top/center of (stereographic) projection
(C{str}, C{'N[orth]'} or C{'S[outh]'}).
@keyword falsed: False both easting and northing (C{bool}).
@keyword strict: Restrict B{C{lat}} to UPS ranges (C{bool}).
@keyword name: Optional B{C{Ups}} name (C{str}).
@return: The UPS coordinate (B{C{Ups}}) or a
L{UtmUps8Tuple}C{(zone, hemipole, easting, northing,
band, datum, convergence, scale)} if B{C{Ups}} is
C{None}. The C{hemipole} is the C{'N'|'S'} pole,
the UPS projection top/center.
@raise RangeError: If B{C{strict}} and B{C{lat}} outside the
valid UPS bands or if B{C{lat}} or B{C{lon}}
outside the valid range and L{rangerrors}
set to C{True}.
@raise TypeError: If B{C{latlon}} is not ellipsoidal.
@raise ValueError: If B{C{lon}} value is missing or if B{C{latlon}}
is invalid.
@see: Karney's C++ class U{UPS
<https://GeographicLib.SourceForge.io/html/classGeographicLib_1_1UPS.html>}.
'''
lat, lon, d, name = _to4lldn(latlon, lon, datum, name)
z, B, p, lat, lon = upsZoneBand5(lat, lon, strict=strict) # PYCHOK UtmUpsLatLon5Tuple
E = d.ellipsoid
p = str(pole or p)[:1].upper()
N = p == 'N' # is north
a = lat if N else -lat
A = abs(a - 90) < _TOL # at pole
t = tan(radians(a))
T = E.es_taupf(t)
r = hypot1(T) + abs(T)
if T >= 0:
r = 0 if A else 1 / r
k0 = getattr(Ups, '_scale0', _K0) # Ups is class or None
r *= 2 * k0 * E.a / E.es_c
k = k0 if A else _scale(E, r, t)
c = lon # [-180, 180) from .upsZoneBand5
x, y = sincos2d(c)
x *= r
y *= r
if N:
y = -y
else:
c = -c
if falsed:
x += _Falsing
y += _Falsing
if Ups is None:
r = UtmUps8Tuple(z, p, x, y, B, d, c, k)
else:
if z != _UPS_ZONE and not strict:
z = _UPS_ZONE # ignore UTM zone
r = Ups(z, p, x, y, band=B, datum=d, falsed=falsed,
convergence=c, scale=k)
r._hemisphere = _hemi(lat)
if isinstance(latlon, _LLEB) and d is latlon.datum:
r._latlon_to(latlon, falsed) # XXX weakref(latlon)?
return _xnamed(r, name)