proj.db: Alter GDA94 transform pipeline accuracies

Prefer to use the conformal+distortion grid to transform between GDA94
and GDA2020, rather than the conformal-only grid or the 7-param
transform.

Additionally alter the GDA94 -> WGS84 transformations to go via
GDA2020 and also use the conformal+distortion grid

Author: craigds

data/sql/koordinates_overrides.sql

@@ -0,0 +1,51 @@
+-- This file is hand generated
+--
+-- When transforming from GDA2020 to WGS84,
+-- use the null transform rather than grid or 7-param transforms
+UPDATE
+    helmert_transformation_table
+SET
+    accuracy = 0.0299
+WHERE
+    auth_name = 'EPSG'
+    AND code = 8450;
+
+--
+-- When transforming from GDA94 to GDA2020,
+-- use the distortion+conformal grid in preference to the conformal-only grid.
+UPDATE
+    grid_transformation
+SET
+    accuracy = 0.0099
+WHERE
+    auth_name = 'EPSG'
+    AND code = 8447;
+
+/*
+ These three transforms from GDA94 to WGS84 are defined:
+ * EPSG:9688, GDA94 to WGS 84 (2), 3.0 m, Australia including Lord Howe Island, Macquarie Island, Ashmore and Cartier Islands, Christmas Island, Cocos (Keeling) Islands, Norfolk Island. All onshore and offshore.
+ * EPSG:1150, GDA94 to WGS 84 (1), 3.0 m, Australia including Lord Howe Island, Macquarie Island, Ashmore and Cartier Islands, Christmas Island, Cocos (Keeling) Islands, Norfolk Island. All onshore and offshore.
+ * DERIVED_FROM(EPSG):9689, GDA94 to WGS 84 (3), 3.0 m, Australia - Australian Capital Territory; New South Wales; Northern Territory; Queensland; South Australia; Tasmania; Western Australia; Victoria.
+ 
+ Of these, the last one is what we want: When projecting from GDA94 to WGS84, assume WGS84 equivalence with GDA2020.
+ Then prioritise the transformation using the conformal+distortion grid followed by a null transform.
+ 
+ So we prioritise that transform by improving its accuracy slightly:
+ */
+UPDATE
+    grid_transformation
+SET
+    accuracy = 2.99
+WHERE
+    auth_name = 'EPSG'
+    AND code = 9689;
+
+-- and this one makes no sense; GDA94 to WGS84 shouldn't be a null transform.
+-- setting the accuracy to null prevents PROJ from choosing this transform over another one.
+update
+    helmert_transformation_table
+set
+    accuracy = NULL
+where
+    auth_name = 'EPSG'
+    and code = 1150;

data/sql_filelist.cmake

@@ -63,6 +63,7 @@ list(APPEND SQL_FILES
   "${SQL_DIR}/wgs84_realizations_concatenated_operations.sql"
   "${SQL_DIR}/customizations.sql"
   "${SQL_DIR}/nkg_post_customizations.sql"
+  "${SQL_DIR}/koordinates_overrides.sql"
 )
 
 if (PROJ_DB_EXTRA_VALIDATION)

test/cli/test_projinfo.yaml

@@ -1403,13 +1403,13 @@ tests:
     unknown id, AGD66 to WGS 84 (22) + WGS 84 to WGS 84 (G1762), 6 m, Papua New Guinea - Papuan fold and thrust belt.
     unknown id, AGD66 to WGS 84 (23) + WGS 84 to WGS 84 (G1762), 6 m, Papua New Guinea - North Fly area (between 5°04'S and 6°36'S and west of 141°32'E).
     unknown id, AGD66 to WGS 84 (12) + WGS 84 to WGS 84 (G1762), 5 m, Australia - Australian Capital Territory.
-    unknown id, AGD66 to GDA94 (11) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 0.75 m, Australia - Australian Capital Territory; New South Wales; Northern Territory; Queensland; South Australia; Tasmania; Western Australia; Victoria., time-dependent operation
-    unknown id, AGD66 to GDA94 (12) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 3.25 m, Australia - Australian Capital Territory; New South Wales; Northern Territory; Queensland; South Australia; Tasmania; Western Australia; Victoria., time-dependent operation
+    unknown id, AGD66 to GDA94 (11) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 0.7099 m, Australia - Australian Capital Territory; New South Wales; Northern Territory; Queensland; South Australia; Tasmania; Western Australia; Victoria., time-dependent operation
+    unknown id, AGD66 to GDA94 (12) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 3.2099 m, Australia - Australian Capital Territory; New South Wales; Northern Territory; Queensland; South Australia; Tasmania; Western Australia; Victoria., time-dependent operation
     unknown id, AGD66 to WGS 84 (17) + WGS 84 to WGS 84 (G1762), 4.9 m, Australia - Australian Capital Territory; New South Wales; Northern Territory; Queensland; South Australia; Tasmania; Western Australia; Victoria.
-    unknown id, AGD66 to GDA94 (9) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 1.25 m, Australia - Northern Territory., time-dependent operation
-    unknown id, AGD66 to GDA94 (4) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 1.25 m, Australia - New South Wales and Victoria., time-dependent operation
-    unknown id, AGD66 to GDA94 (8) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 1.25 m, Australia - Tasmania including islands - onshore., time-dependent operation
-    unknown id, AGD66 to GDA94 (19) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 0.75 m, Australia - Australian Capital Territory., time-dependent operation
+    unknown id, AGD66 to GDA94 (9) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 1.2099 m, Australia - Northern Territory., time-dependent operation
+    unknown id, AGD66 to GDA94 (4) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 1.2099 m, Australia - New South Wales and Victoria., time-dependent operation
+    unknown id, AGD66 to GDA94 (8) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 1.2099 m, Australia - Tasmania including islands - onshore., time-dependent operation
+    unknown id, AGD66 to GDA94 (19) + GDA94 to GDA2020 (2) + Conversion from GDA2020 (geog2D) to GDA2020 (geocentric) + GDA2020 to WGS 84 (G1762) (1) + Conversion from WGS 84 (G1762) (geocentric) to WGS 84 (G1762) (geog2D), 0.7099 m, Australia - Australian Capital Territory., time-dependent operation
     unknown id, Ballpark geographic offset from AGD66 to WGS 84 (G1762), unknown accuracy, World, has ballpark transformation
 - args: -s EPSG:31467 -t ETRS89 --spatial-test intersects --grid-check none --bbox 8,48,9,49 --summary
   out: |