Is this exception “osgeo::proj::io::FactoryException” a bug? in 9.4.0 about proj

[rouault]

When I compile and run the qgis project, I got this exception.

This exception might be emitted by PROJ with such an unsual semi-major axis, but it should be caught, at least in the C function that QGIS uses: which is it?
With only the content of this report, it is not clear if there's actually something to change in PROJ or in QGIS
Might be worth to open it rather as a QGIS issue with the description of the related QGIS scenario

[loftyscholar]

Thanks for your reply.
This unsual semi-major axis is located in ellipsoid table in proj.db, and its code is "107911" and its name is "Ananke_2000_IAU_IAG".

[rouault]

@loftyscholar
ok, I see indeed that's an issue when trying to do something like projinfo -s "+proj=longlat +a=10000 +no_defs +type=crs" -t ESRI:104911 which should be a no-op operation, because there are several celestial bodies that match the 10000 radius. But couldn't you just avoid using a deprecated PROJ.4 string like "+proj=longlat +a=10000 +no_defs +type=crs" and use the ESRI:104911 code in QGIS. That should avoid the issue.

[rouault]

or use the WKT2 representation:

GEOGCRS["GCS_Ananke_2000",
    DATUM["D_Ananke_2000",
        ELLIPSOID["Ananke_2000_IAU_IAG",10000,0,
            LENGTHUNIT["metre",1]]],
    PRIMEM["Reference_Meridian",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    USAGE[
        SCOPE["Not known."],
        AREA["World."],
        BBOX[-90,-180,90,180]],
    ID["ESRI",104911]]

Alternatively you may define the PROJ_IGNORE_CELESTIAL_BODY=YES environment variable to avoid an error like

bin/projinfo -s  "+proj=longlat +R=10000 +no_defs +type=crs" -t ESRI:104911
createOperations() failed with: Source and target ellipsoid do not belong to the same celestial body (Non-Earth body vs Ananke). You may override this check by setting the PROJ_IGNORE_CELESTIAL_BODY environment variable to YES.

[loftyscholar]

deprecated
Ok, I'll have a try soon.

[rouault]

Improvements in PROJ queued in #4149

With those, projinfo -s "+proj=longlat +R=10000 +no_defs +type=crs" -t ESRI:104911 now works out of the box

[loftyscholar]

Improvements in PROJ queued in #4149

With those, projinfo -s "+proj=longlat +R=10000 +no_defs +type=crs" -t ESRI:104911 now works out of the box

Thanks for your reply, have all the changes been commited to the proj source code? Where should I get the changes?

[rouault]

have all the changes been commited to the proj source code?

This has just been committed in the master branch, and will be soon in the 9.4 branch as well. This will be in the 9.4.1 release