[gdal-dev] Problems with several dozen projections
Javier Jimenez Shaw
j1 at jimenezshaw.com
Thu Jul 31 00:49:49 PDT 2025
Looking a bit more detail the list I see different reasons
A few non implemented projections in PROJ:
"Lambert Conic Conformal (West Orientated)"
"Bonne (South Orientated)"
Engineering CRSs
Some are height systems. The transformation make no sense
Others are not defined from projected to geographic. Those projections do
not have an inverse function in PROJ.
Some ESRI methods may not be implemented.
On Wed, 30 Jul 2025 at 15:07, Javier Jimenez Shaw <j1 at jimenezshaw.com>
wrote:
> 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/20250731/2e803fe7/attachment.htm>
More information about the gdal-dev
mailing list