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

Reply via email to