def _to3zBll(lat, lon, cmoff=True):
'''(INTERNAL) Return zone, Band and lat- and (central) longitude in degrees.
@arg lat: Latitude (C{degrees}).
@arg lon: Longitude (C{degrees}).
@kwarg cmoff: Offset B{C{lon}} from zone's central meridian.
@return: 4-Tuple (zone, Band, lat, lon).
'''
z, lat, lon = _to3zll(lat, lon) # in .utmupsBase
if _UTM_LAT_MIN > lat or lat >= _UTM_LAT_MAX: # [-80, 84) like Veness
t = ' '.join((_outside_, _UTM_, _range_, '[%s,' % (_UTM_LAT_MIN,),
'%s)' % (_UTM_LAT_MAX,)))
raise RangeError(lat=degDMS(lat), txt=t)
B = _Bands[int(lat + 80) >> 3]
x = lon - _cmlon(z) # z before Norway/Svaldbard
if abs(x) > _UTM_ZONE_OFF_MAX:
t = ' '.join((_outside_, _UTM_, _zone_, str(z), 'by', degDMS(x, prec=6)))
raise RangeError(lon=degDMS(lon), txt=t)
if B == 'X': # and 0 <= int(lon) < 42: z = int(lon + 183) // 6 + 1
x = {32: 9, 34: 21, 36: 33}.get(z, None)
if x: # Svalbard
z += 1 if lon >= x else -1
elif B == 'V' and z == 31 and lon >= 3:
z += 1 # SouthWestern Norway
if cmoff: # lon off central meridian
lon -= _cmlon(z) # z after Norway/Svaldbard
return Zone(z), Band(B), Lat(lat), Lon(lon)
def upsZoneBand5(lat, lon, strict=True):
'''Return the UTM/UPS zone number, (polar) Band letter, pole and
clipped lat- and longitude for a given location.
@arg lat: Latitude in degrees (C{scalar} or C{str}).
@arg lon: Longitude in degrees (C{scalar} or C{str}).
@kwarg strict: Restrict B{C{lat}} to UPS ranges (C{bool}).
@return: A L{UtmUpsLatLon5Tuple}C{(zone, band, hemipole,
lat, lon)} where C{hemipole} is the C{'N'|'S'} pole,
the UPS projection top/center.
@raise RangeError: If B{C{strict}} and B{C{lat}} in the UTM
and not the UPS range or if B{C{lat}} or
B{C{lon}} outside the valid range and
L{rangerrors} set to C{True}.
@raise ValueError: Invalid B{C{lat}} or B{C{lon}}.
'''
z, lat, lon = _to3zll(*parseDMS2(lat, lon))
if lat < _UPS_LAT_MIN: # includes 30' overlap
z, B, p = _UPS_ZONE, _Band(lat, lon), _S_
elif lat > _UPS_LAT_MAX: # includes 30' overlap
z, B, p = _UPS_ZONE, _Band(lat, lon), _N_
elif strict:
t = ' '.join((_inside_, _UTM_, _range_, '[%s,' % (_UPS_LAT_MIN,),
'%s]' % (_UPS_LAT_MAX,)))
raise RangeError(lat=degDMS(lat), txt=t)
else:
B, p = NN, _hemi(lat)
return UtmUpsLatLon5Tuple(z, B, p, lat, lon, Error=UPSError)
def _clipped_(angle, limit, units):
'''(INTERNAL) Helper for C{clipDegrees} and C{clipRadians}.
'''
c = min(limit, max(-limit, angle))
if c != angle and _rangerrors:
t = NN(fstr(angle, prec=6, ints=True), 'beyond',
copysign(limit, angle), units)
raise RangeError(t, txt=None)
return c