Hi gdal-dev list, Not sure if this is a bug so posting here first...
I am having to reverse X and Y coordinates in a coordinate transformations (in at least 3 different dst_crs) otherwise I am getting -inf values. Python bindings with GDAL 3.0.1, Python 3.7.2 64-bit on Windows 10. For a point near the dateline in the Arctic at (175W, 80N): from osgeo import osr in_crs = osr.SpatialReference() in_crs.ImportFromEPSG(4326) out_crs = osr.SpatialReference() out_crs.ImportFromEPSG(3347) transformation = osr.CoordinateTransformation(in_crs, out_crs) print(transformation.TransformPoint(-175, 80, 0)) # lon, lat print(transformation.TransformPoint(80, -175, 0)) # lat, lon Gives : (inf, inf, inf) (4983944.437216174, 5754720.297247784, 0.0) # correct https://epsg.io/transform#s_srs=4326&t_srs=3347&x=-175.0000000&y=80.0000000 Regards, Benjamin Deschamps Environment and Climate Change Canada / Government of Canada
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/gdal-dev
