Dear All Happy New Year!
I'm trying to use the GetAreaOfUse() method on a osgeo.osr.SpatialReference to check the validity of coordinates. With the following for the Irish National Grid (EPSG:29902) I get the WGS84 bounds - west, east, south, north - of the spatial reference system to be -10.56, -5.93, 51.39, 55.43: >>> from osgeo import osr >>> srs = osr.SpatialReference() >>> srs.ImportFromEPSG(29902) 0 >>> aou = srs.GetAreaOfUse() >>> aou.west_lon_degree, aou.east_lon_degree, aou.south_lat_degree, >>> aou.north_lat_degree (-10.56, -5.93, 51.39, 55.43) However, at https://epsg.io/29902 the WGS84<https://epsg.io/29902%20the%20WGS84> bounds are listed as -10.56, -5.34, 51.39, 55.43. How is the GetAreaOfUse() method implemented? And is there a bug in the values stored for the EPSG:29902 projection? Or is there another reason why the values differ to those listed at epsg.io. Many thanks, best wishes, Simon. Simon Wright JBA Risk Management Limited -- e: simon.wri...@jbarisk.com d:+44 (0)1756 587258 t: +44 (0)1756 799919 www.jbarisk.com All JBA Risk Management's email messages contain confidential information and are intended only for the individual(s) named. If you are not the named addressee you should not disseminate, distribute or copy this e-mail. Please notify the sender immediately by email if you have received this email by mistake and delete this email from your system. JBA Risk Management Limited is registered in England, company number 07732946, 1 Broughton Park, Old Lane North, Broughton, Skipton, North Yorkshire, BD23 3FD, England.
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev