[PROJ] [PROJ 6.0.0] Different behaviors for creating a transformation of type EPSG:9616

[Didier Richard]

1. - Shame on me, it is working !

2.- Tested on IGNF:NTFLAMB2E.NGF84 to IGNF:ETRS89LCC.EVRF2000
but the final type is PJ_TYPE_UNKNOWN, not PJ_TYPE_OTHER_COORDINATE_OPERATION :

    obj = proj_create_crs_to_crs(c, "IGNF:NTFLAMB2E.NGF84", "IGNF:ETRS89LCC.EVRF2000", NULL);
    if (obj != NULL) {
        fprintf(stderr,"%s is of type %d\n", proj_get_name(obj), proj_get_type(obj));
        proj_destroy(obj);
    }

outputs :

(null) is of type 0

$ projinfo -o PROJ -s IGNF:NTFLAMB2E.NGF84 -t IGNF:ETRS89LCC.EVRF2000
Candidate operations found: 3
-------------------------------------
Operation n°1:

unknown id, Inverse of LAMBERT II ETENDU + Transformation from NGF-LALLEMAND to EVRF2000 (UELN-95/98)(EUROPEAN VERTICAL REFERENCE FRAME) (ballpark vertical transformation) + NTF geographiques Paris (gr) vers NTF GEOGRAPHIQUES GREENWICH (DMS) + NOUVELLE TRIANGULATION DE LA FRANCE (NTF) vers WGS 84 + Inverse of ETRS 89 vers WGS 84 + ETRS89 LAMBERT CONFORMAL CONIC, unknown accuracy, FRANCE METROPOLITAINE, has ballpark transformation

PROJ string:
+proj=pipeline +step +inv +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +ellps=clrk80ign +pm=paris +step +proj=push +v_3 +step +proj=cart +ellps=clrk80ign +step +proj=helmert +x=-168 +y=-60 +z=320 +step +inv +proj=cart +ellps=GRS80 +step +proj=pop +v_3 +step +proj=lcc +lat_0=52 +lon_0=10 +lat_1=35 +lat_2=65 +x_0=4000000 +y_0=2800000 +ellps=GRS80

-------------------------------------
Operation n°2:

unknown id, Inverse of LAMBERT II ETENDU + Transformation from NGF-LALLEMAND to EVRF2000 (UELN-95/98)(EUROPEAN VERTICAL REFERENCE FRAME) (ballpark vertical transformation) + NTF geographiques Paris (gr) vers NTF GEOGRAPHIQUES GREENWICH (DMS) + NOUVELLE TRIANGULATION DE LA FRANCE (NTF) vers WGS 84 + Inverse of ETRS 89 vers WGS 84 + ETRS89 LAMBERT CONFORMAL CONIC, unknown accuracy, FRANCE METROPOLITAINE, has ballpark transformation

PROJ string:
+proj=pipeline +step +inv +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +ellps=clrk80ign +pm=paris +step +proj=push +v_3 +step +proj=cart +ellps=clrk80ign +step +proj=helmert +x=-168 +y=-60 +z=320 +step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3 +step +proj=lcc +lat_0=52 +lon_0=10 +lat_1=35 +lat_2=65 +x_0=4000000 +y_0=2800000 +ellps=GRS80

-------------------------------------
Operation n°3:

unknown id, Inverse of LAMBERT II ETENDU + Transformation from NGF-LALLEMAND to EVRF2000 (UELN-95/98)(EUROPEAN VERTICAL REFERENCE FRAME) (ballpark vertical transformation) + Inverse of Transformation from ETRS89 geographiques (dms) to NTF geographiques Paris (gr) altered to use prime meridian of ETRS89 geographiques (dms) + Ballpark geographic offset from NTF geographiques Paris (gr) altered to use prime meridian of ETRS89 geographiques (dms) to ETRS89 geographiques (dms) + ETRS89 LAMBERT CONFORMAL CONIC, unknown accuracy, World, has ballpark transformation

PROJ string:
+proj=pipeline +step +inv +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +ellps=clrk80ign +pm=paris +step +proj=lcc +lat_0=52 +lon_0=10 +lat_1=35 +lat_2=65 +x_0=4000000 +y_0=2800000 +ellps=GRS80


--
RICHARD Didier - Chef du Centre de Compétences Technologies des Systèmes d'Information
http://fr.linkedin.com/pub/didier-richard/98/2a3/a8/ - https://www.osgeo.org/member/didier/
IGN/Direction des Sciences et Technologies de l'Information/ENSG Géomatique
6/8 avenue Blaise Pascal - BP Champs-sur-Marne - 77455 MARNE-LA-VALLÉE CEDEX 2
Tél : +33 (0) 1 43 98 83 23


________________________________________
De : Even Rouault [even.rouault at spatialys.com]
Envoyé : mardi 30 avril 2019 16:19
À : Didier Richard
Cc : proj at lists.osgeo.org
Objet : Re: [PROJ] [PROJ 6.0.0] Different behaviors for creating a transformation of type EPSG:9616

On mardi 30 avril 2019 14:02:50 CEST Didier Richard wrote:
> 1.- Does not work either with
> "urn:ogc:def:coordinateOperation:IGNF::SG1250". I guess because of the
> authority ?

no just a typo. Missing T before SG1250

> 3.- Ok, I thought that operations located in other_transformation table were
> of type PJ_TYPE_OTHER_COORDINATE_OPERATION. Do you have an example of an
> operation for which the type is PJ_TYPE_OTHER_COORDINATE_OPERATION ? I
> found nothing in the test/ directory checking this.

For example, if you would do a transformation from a CompoundCRS to another
CompoundCRS (or a Geographic3D CRS) with both horizontal and vertical
transformations, there is nothing in ISO:19111 to model the composite of both.
So there's an internal PROJ C++ class, PROJBasedOperation, whose type would be
exposed as PJ_TYPE_OTHER_COORDINATE_OPERATION

--
Spatialys - Geospatial professional services
http://www.spatialys.com