[PROJ] [PROJ 6.1.0] Tests failed on vertical grid shift for France mainland onshore

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

[PROJ] [PROJ 6.1.0] Tests failed on vertical grid shift for France mainland onshore

Didier Richard
Hi all,

While running make check after a fresh install, the following error appeared :

...
============================================
Running ../../test/cli/testprojinfo using ../../src/projinfo:
============================================
diff testprojinfo_out with testprojinfo_out.dist
--- testprojinfo_out    2019-05-18 10:52:56.203561622 +0200
+++ ../../test/cli/testprojinfo_out.dist        2019-05-06 07:12:09.000000000 +0200
@@ -762,18 +762,18 @@
 -------------------------------------
 Operation n°1:
 
-INVERSE(DERIVED_FROM(EPSG)):8371, Inverse of RGF93 to NGF IGN69 height (2), 0.02 m, France - mainland onshore
+INVERSE(DERIVED_FROM(EPSG)):8885, Inverse of RGF93 to NGF IGN69 height (3), 0.01 m, France - mainland onshore
 
 PROJ string:
-+proj=vgridshift +grids=RAF09.gtx +multiplier=1
++proj=vgridshift +grids=RAF18.gtx +multiplier=1
 
 -------------------------------------
 Operation n°2:
 
-INVERSE(DERIVED_FROM(EPSG)):8885, Inverse of RGF93 to NGF IGN69 height (3), 0.01 m, France - mainland onshore
+INVERSE(DERIVED_FROM(EPSG)):8371, Inverse of RGF93 to NGF IGN69 height (2), 0.02 m, France - mainland onshore
 
 PROJ string:
-+proj=vgridshift +grids=RAF18.gtx +multiplier=1
++proj=vgridshift +grids=RAF09.gtx +multiplier=1
 
 -------------------------------------
 Operation n°3:

PROBLEMS HAVE OCCURRED
test file testprojinfo_out saved


RAF18.gtx is related to France mainland onshore while RAF09.gtx is related to French Antilla.

Looking at sql files :

sql/grid_transformation.sql:INSERT INTO "grid_transformation" VALUES('EPSG','8371','RGF93 to NGF IGN69 height (2)',NULL,NULL,'EPSG','1073','Geographic3D to GravityRelatedHeight (IGN2009)','EPSG','4965','EPSG','5720','EPSG','1326',0.02,'EPSG','8666','Geoid (height correction) model file','RAF09.mnt',NULL,NULL,NULL,NULL,NULL,NULL,'IGN Fra 09',0);
sql/grid_transformation.sql:INSERT INTO "grid_transformation" VALUES('EPSG','8885','RGF93 to NGF IGN69 height (3)',NULL,NULL,'EPSG','1073','Geographic3D to GravityRelatedHeight (IGN2009)','EPSG','4965','EPSG','5720','EPSG','1326',0.01,'EPSG','8666','Geoid (height correction) model file','RAF18.tac',NULL,NULL,NULL,NULL,NULL,NULL,'IGN Fra 18',0);

sql/supersession.sql:INSERT INTO "supersession" VALUES('grid_transformation','EPSG','8371','grid_transformation','EPSG','8885','EPSG');

-- Remove supression entry of EPSG:8371 ("RGF93 to NGF IGN69 height (2)", using RAF09.mnt) by
-- EPSG:8885 ("RGF93 to NGF IGN69 height (3)", using RAF18.tac), since
-- we have only RAF09.mnt in proj-datumgrid-europe
-- Our code to remove superseded operations should probably take into account grid availability
DELETE FROM "supersession" WHERE superseded_table_name = 'grid_transformation' AND
                                 superseded_auth_name = 'EPSG' AND
                                 superseded_code = '8371' AND
                                 replacement_table_name = 'grid_transformation' AND
                                 replacement_auth_name = 'EPSG' AND
                                 replacement_code = '8885' AND
                                 source = 'EPSG';

There seams to have no difference between these lines in 6.0.0 and 6.10 sql files.
Any idea ?

Best regards,
--
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

_______________________________________________
PROJ mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/proj
Reply | Threaded
Open this post in threaded view
|

Re: [PROJ 6.1.0] Tests failed on vertical grid shift for France mainland onshore

Even Rouault-2
Hi Didier,

I assume that you have run the tests with the RAF09.gtx file present and not
the RAF18.gtx one ?
If so the failure is expected. 'make check' will run fine with just the
default proj-datumgrid is installed, but the presence of other grids might
cause changes in results. If the new grid RAF18.gtx from proj-datumgrid-
europe-1.3 is present, you'll get also a slightly different result. Nothing to
worry about.

> RAF18.gtx is related to France mainland onshore while RAF09.gtx is related
> to French Antilla.

Hum, not according to the EPSG dataset... Both are registered against
EPSG:1326 area of use 'France - mainland onshore'. Maybe you are confusing
RAF09.gtx (the former version of RAF18.gtx) with RGAF09 which is indeed for
French Antilla.

> -- Remove supression entry of EPSG:8371 ("RGF93 to NGF IGN69 height (2)",
> using RAF09.mnt) by -- EPSG:8885 ("RGF93 to NGF IGN69 height (3)", using
> RAF18.tac), since -- we have only RAF09.mnt in proj-datumgrid-europe
> -- Our code to remove superseded operations should probably take into
> account grid availability DELETE FROM "supersession" WHERE
> superseded_table_name = 'grid_transformation' AND superseded_auth_name =
> 'EPSG' AND superseded_code = '8371' AND
>                                  replacement_table_name =
> 'grid_transformation' AND replacement_auth_name = 'EPSG' AND
> replacement_code = '8885' AND
>                                  source = 'EPSG';
>

The above entry dates back indeed to PROJ 6.0 / proj-datumgrid-europe-1.2 when
there was only the RAF09.gtx grid available, and not RAF18. The removal of
this deprecation entry doesn't hurt by itself, although it is no longer be
necessary. So I've just removed it in master.

Even

--
Spatialys - Geospatial professional services
http://www.spatialys.com
_______________________________________________
PROJ mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/proj