[gdal-dev] Problems with several dozen projections
Javier Jimenez Shaw
j1 at jimenezshaw.com
Wed Jul 30 06:07:12 PDT 2025
When you say "I was running a test of all available projections in PROJ
9.6.2 using GDAL 3.11.3", you mean all the available Coordinate Reference
Systems (CRS).
Some of those CRS have a projection that is not included in PROJ. For
instance, the first in the list, EPSG:2299 says
"Error when exporting to PROJ string: Unsupported conversion method:
Lambert Conic Conformal (West Orientated)"
On Wed, 30 Jul 2025 at 12:13, Ray at Daylon via gdal-dev <
gdal-dev at lists.osgeo.org> wrote:
> I was running a test of all available projections in PROJ 9.6.2 using
> GDAL 3.11.3 (both current as of this writing), on Windows.
>
> I got errors trying to project from lat/lon and back with the 55
> projections listed in the attached CSV file. If I leave off the
> SetAxisMappingStrategy calls, many more projections fail.
> No deprecated or celestial body projections were tested.
>
> When a transformation failed to go from lat/lon to the projected CRS,
> the reason was usually that the transformation wound up with an infinity.
>
>
> My test program is essentially this:
>
> LatLonCRS = OGRSpatialReference(EPSG:4326)
> LatLonCRS.SetAxisMappingStrategy(OAMS_TRADITIONAL_GIS_ORDER)
>
> Projections = proj_get_crs_info_list_from_database()
>
> For each Projection
> Wkt2 = proj_as_wkt(Projection, PJ_WKT2_2019)
>
> ProjCRS = OGRSpatialReference(Wkt2)
> ProjCRS.SetAxisMappingStrategy(OAMS_TRADITIONAL_GIS_ORDER)
>
> Xfm1 = OGRCreateCoordTransform(LatLonCRS, ProjCRS)
> Xfm2 = OGRCreateCoordTransform(ProjCRS, LatLonCRS)
>
> ProjCRS.GetAreaOfUse(West, South, East, North)
> // Get center of bbox
> X = (West + East) / 2
> Y = (North + South) / 2
>
> Xfm1.Transform(X, Y) // Go from lat/lon to Proj CRS
> Xfm2.Transform(X, Y) // and back to lat/lon
>
> Next projection
>
> There were various errors logged by PROJ:
>
> Cannot find coordinate operations from `EPSG:2299' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:2301' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:2303' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:2304' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:2305' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:2306' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:2307' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:2963' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:3144' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5017' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5800' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5801' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5802' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5803' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5808' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5809' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5810' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5811' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5812' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5813' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5814' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5815' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5816' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:5817' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:6715' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:10642' to `EPSG:4326'
> Cannot find coordinate operations from `EPSG:10740' to `EPSG:4326'
> No inverse operation
> Point outside of projection domain
> PROJ: proj_crs_get_geodetic_crs: CRS has no geodetic CRS
> PROJ: omerc: Invalid value for lat_1: lat_1 should be different from 0
> PROJ: pipeline: Pipeline: Bad step definition: inv (Invalid value for an
> argument)
> PROJ: pipeline: Pipeline: Bad step definition: proj=omerc (Invalid value
> for an argument)
> PROJ: pipeline: Pipeline: Inverse operation for wag7 is not available
>
> Are any of these bugs in PROJ or GDAL?
>
> If so, are there any workarounds?
>
> I was also wondering, should it be necessary to take a projection's axis
> mapping into account? Isn't that a detail that GDAL or PROJ can abstract
> away? Is SetAxisMappingStrategy() really needed?
>
> My thanks in advance for any help,
> Ray Gardener
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20250730/7c8ab6c0/attachment.htm>
More information about the gdal-dev
mailing list