def elevation2(lat, lon, timeout=2.0):
'''Get the geoid elevation at an C{NAD83} to C{NAVD88} location.
@arg lat: Latitude (C{degrees}).
@arg lon: Longitude (C{degrees}).
@kwarg timeout: Optional, query timeout (seconds).
@return: An L{Elevation2Tuple}C{(elevation, data_source)}
or (C{None, "error"}) in case of errors.
@raise ValueError: Invalid B{C{timeout}}.
@note: The returned C{elevation} is C{None} if B{C{lat}} or B{C{lon}}
is invalid or outside the C{Conterminous US (CONUS)},
if conversion failed or if the query timed out. The
C{error} is the C{HTTP-, IO-, SSL-, Type-, URL-} or
C{ValueError} as a string (C{str}).
@see: U{USGS National Map<https://NationalMap.gov/epqs>},
the U{FAQ<https://www.USGS.gov/faqs/what-are-projection-
horizontal-and-vertical-datum-units-and-resolution-3dep-standard-dems>},
U{geoid.py<https://Gist.GitHub.com/pyRobShrk>}, module
L{geoids}, classes L{GeoidG2012B}, L{GeoidKarney} and
L{GeoidPGM}.
'''
try:
x = _qURL(
'https://NED.USGS.gov/epqs/pqs.php', # 'https://NationalMap.gov/epqs/pqs.php'
(
'%s=%.6F' % (_x_, Lon(lon)),
'%s=%.6F' % (_y_, Lat(lat)),
'%s=%s' % (_units_, 'Meters'), # 'Feet', capitalized
'output=xml'), # json, case_sensitive
timeout=Scalar(timeout, name=_timeout_))
if x[:6] == '<?xml ':
e = _xml('Elevation', x)
try:
e = float(e)
if -1000000 < e < 1000000:
return Elevation2Tuple(e, _xml('Data_Source', x))
e = 'non-CONUS %.2F' % (e, )
except (TypeError, ValueError):
pass
else:
e = 'no %s "%s"' % (
_XML_,
clips(x, limit=128, white=_SPACE_),
)
except (HTTPError, IOError, TypeError, ValueError) as x:
e = repr(x)
e = _error(elevation2, lat, lon, e)
return Elevation2Tuple(None, e)
def geoidHeight2(lat, lon, model=0, timeout=2.0):
'''Get the C{NAVD88} geoid height at an C{NAD83} location.
@arg lat: Latitude (C{degrees}).
@arg lon: Longitude (C{degrees}).
@kwarg model: Optional, geoid model ID (C{int}).
@kwarg timeout: Optional, query timeout (seconds).
@return: An L{GeoidHeight2Tuple}C{(height, model_name)}
or C{(None, "error"}) in case of errors.
@raise ValueError: Invalid B{C{timeout}}.
@note: The returned C{height} is C{None} if B{C{lat}} or B{C{lon}} is
invalid or outside the C{Conterminous US (CONUS)}, if the
B{C{model}} was invalid, if conversion failed or if the
query timed out. The C{error} is the C{HTTP-, IO-, SSL-,
Type-, URL-} or C{ValueError} as a string (C{str}).
@see: U{NOAA National Geodetic Survey
<https://www.NGS.NOAA.gov/INFO/geodesy.shtml>},
U{Geoid<https://www.NGS.NOAA.gov/web_services/geoid.shtml>},
U{USGS10mElev.py<https://Gist.GitHub.com/pyRobShrk>}, module
L{geoids}, classes L{GeoidG2012B}, L{GeoidKarney} and
L{GeoidPGM}.
'''
try:
j = _qURL('https://Geodesy.NOAA.gov/api/geoid/ght',
('lat=%.6F' % (Lat(lat), ), 'lon=%.6F' %
(Lon(lon), ), 'model=%s' % (model, ) if model else ''),
timeout=Scalar(timeout, name='timeout')) # PYCHOK 5
if j[:1] == '{' and j[-1:] == '}' and j.find('"error":') > 0:
d = _json(j)
if isinstance(d.get('error', 'N/A'), float):
h = d.get('geoidHeight', None)
if h is not None:
m = _str(d.get('geoidModel', 'N/A'))
return GeoidHeight2Tuple(h, m)
e = 'geoidHeight'
else:
e = 'JSON'
e = 'no %s "%s"' % (e, clips(j, limit=256, white=' '))
except (HTTPError, IOError, TypeError, ValueError) as x:
e = repr(x)
e = _error(geoidHeight2, lat, lon, e)
return GeoidHeight2Tuple(None, e)
def clipStr(bstr, limit=50, white=NN):
'''DEPRECATED, use function L{clips}.
'''
from pygeodesy.basics import clips
return clips(bstr, limit=limit, white=white)