Package pygeodesy :: Module osgr
[frames] | no frames]

Module osgr

Ordinance Survey Grid References (OSGR) references.

Classes Osgr and OSGRError and functions parseOSGR and toOsgr.

Pure Python implementation of OS Grid Reference functions using an ellipsoidal earth model, transcoded from JavaScript originals by (C) Chris Veness 2005-2016 published under the same MIT Licence**, see OS National Grid and Module osgridref.

OSGR provides geocoordinate references for UK mapping purposes, converted in 2015 to work with WGS84 datum by default or OSGB36 as option.

See Guide, Proposed Changes, Confirmation and Ordnance Survey National Grid.

See also Karney 'Transverse Mercator with an accuracy of a few nanometers', 2011 (building on Krüger 'Konforme Abbildung des Erdellipsoids in der Ebene', 1912), Seidel 'Die Mathematik der Gauß-Krueger-Abbildung', 2006 and Transverse Mercator: Redfearn series.


Version: 21.06.09

Classes
  OSGRError
Ordinance Survey Grid References (OSGR) parse or other Osgr issue.
  Osgr
Ordinance Survey Grid References (OSGR) coordinate.
Functions
 
parseOSGR(strOSGR, Osgr=<class 'pygeodesy.osgr.Osgr'>, name='')
Parse a string representing an OSGR grid reference, consisting of "[grid] easting northing".
 
toOsgr(latlon, lon=None, datum=Datum(name='WGS84', ellipsoid=Ellipsoids.WGS84, transform=Tran..., Osgr=<class 'pygeodesy.osgr.Osgr'>, name='', **Osgr_kwds)
Convert a lat-/longitude point to an OSGR coordinate.
Variables
  __all__ = _ALL_LAZY.osgr
Function Details

parseOSGR (strOSGR, Osgr=<class 'pygeodesy.osgr.Osgr'>, name='')

 

Parse a string representing an OSGR grid reference, consisting of "[grid] easting northing".

Accepts standard OS Grid References like 'SU 387 148', with or without whitespace separators, from 2- up to 10-digit references (1 m × 1 m square), or all-numeric, comma-separated references in meters, for example '438700,114800'.

Arguments:
  • strOSGR - An OSGR coordinate (str).
  • Osgr - Optional class to return the OSGR coordinate (Osgr) or None.
  • name - Optional Osgr name (str).
Returns:
The OSGR coordinate (Osgr) or an EasNor2Tuple(easting, northing) if Osgr is None.
Raises:

Example:

>>> g = parseOSGR('TG 51409 13177')
>>> str(g)  # TG 51409 13177
>>> g = parseOSGR('TG5140913177')
>>> str(g)  # TG 51409 13177
>>> g = parseOSGR('TG51409 13177')
>>> str(g)  # TG 51409 13177
>>> g = parseOSGR('651409,313177')
>>> str(g)  # TG 51409 13177
>>> g.toStr(prec=0)  # 651409,313177

toOsgr (latlon, lon=None, datum=Datum(name='WGS84', ellipsoid=Ellipsoids.WGS84, transform=Tran..., Osgr=<class 'pygeodesy.osgr.Osgr'>, name='', **Osgr_kwds)

 

Convert a lat-/longitude point to an OSGR coordinate.

Arguments:
  • latlon - Latitude (degrees) or an (ellipsoidal) geodetic LatLon point.
  • lon - Optional longitude in degrees (scalar or None).
  • datum - Optional datum to convert lat, lon from (Datum, Ellipsoid, Ellipsoid2 or a_f2Tuple).
  • Osgr - Optional class to return the OSGR coordinate (Osgr) or None.
  • name - Optional Osgr name (str).
  • Osgr_kwds - Optional, additional Osgr keyword arguments, ignored if Osgr=None.
Returns:
The OSGR coordinate (Osgr) or an EasNor2Tuple(easting, northing) if Osgr is None.
Raises:
  • OSGRError - Invalid latlon or lon.
  • TypeError - Non-ellipsoidal latlon or invalid datum or conversion failed.

Example:

>>> p = LatLon(52.65798, 1.71605)
>>> r = toOsgr(p)  # TG 51409 13177
>>> # for conversion of (historical) OSGB36 lat-/longitude:
>>> r = toOsgr(52.65757, 1.71791, datum=Datums.OSGB36)