[Qgis-user] Using a NTv2 grid file for transforming between different CRSs.

Thomas Wunder foren.thomas at thomas-wunder.de
Mon Jul 27 02:02:09 PDT 2020


Hello,
I've created a new project (CRS=EPSG:25832 ETRS89/UTM32N)  using QGIS 3.10.8 
(on Debian 11, PROJ version 7.1.0) which contains several multi polygon layers 
(land parcels). Those layers were imported from shape files containing either 
EPSG:25832 (as is the project) or legacy EPSG:31468 (DHDN/3° Gauss-Krueger 
Zone 4) data.  QGIS gave me a couple of transformation options via Project -> 
Properties -> CRS -> Datum Transformations but none of those is accurate 
enough for my purposes. I need sub centimeter accuracy.

Therefore I'd like to go for an NTv2 grid based approach using a gsb file. I've 
put the grid file to /usr/share/proj/kanu_ntv2_oberfranken.gsb and added 
several rows to /usr/share/qgis/resources/srs.db like so (don't remember all 
the variations):

INSERT INTO tbl_datum_transform VALUES (99999, 1000997, 4314, 4326, 9615, 
'kanu_ntv2_mittelfranken.gsb', NULL, NULL, NULL, NULL, NULL, NULL, 'Für ALKIS 
Bayern, amtliche Gitterdatei des LDBV', 'Für Genauigkeitsanforderungen im cm 
Bereich', 1, 0, NULL);

or

INSERT INTO tbl_datum_transform VALUES (99999, 1000998, 4314, 4258, 9615, 
'kanu_ntv2_mittelfranken.gsb', NULL, NULL, NULL, NULL, NULL, NULL, 'Für ALKIS 
Bayern, amtliche Gitterdatei des LDBV', 'Für Genauigkeitsanforderungen im cm 
Bereich', 1, 0, 3339);

No matter what I add there QGIS never lists a new transformation option in 
Project -> Properties -> CRS -> Datum Transformations -> '+' when selecting 
Source = EPSG:31468 and Destination = EPSG:25832. 

I only get
+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 
+y_0=0 +ellps=bessel +step +proj=hgridshift +grids=de_adv_BETA2007.tif +step 
+proj=utm +zone=32 +ellps=GRS80

+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 
+y_0=0 +ellps=bessel +step +proj=push +v_3 +step +proj=cart +ellps=bessel 
+step +proj=helmert +x=598.1 +y=73.7 +z=418.2 +rx=0.202 +ry=0.045 +rz=-2.455 
+s=6.7 +convention=position_vector +step +inv +proj=cart +ellps=GRS80 +step 
+proj=pop +v_3 +step +proj=utm +zone=32 +ellps=GRS80

+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 
+y_0=0 +ellps=bessel +step +proj=push +v_3 +step +proj=cart +ellps=bessel 
+step +proj=helmert +x=597.1 +y=71.4 +z=412.1 +rx=0.894 +ry=0.068 +rz=-1.563 
+s=7.58 +convention=position_vector +step +inv +proj=cart +ellps=GRS80 +step 
+proj=pop +v_3 +step +proj=utm +zone=32 +ellps=GRS80

+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 
+y_0=0 +ellps=bessel +step +proj=push +v_3 +step +proj=cart +ellps=bessel 
+step +proj=helmert +x=584.8 +y=67 +z=400.3 +rx=0.105 +ry=0.013 +rz=-2.378 
+s=10.29 +convention=position_vector +step +inv +proj=cart +ellps=GRS80 +step 
+proj=pop +v_3 +step +proj=utm +zone=32 +ellps=GRS80

+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 
+y_0=0 +ellps=bessel +step +proj=push +v_3 +step +proj=cart +ellps=bessel 
+step +proj=helmert +x=590.5 +y=69.5 +z=411.6 +rx=-0.796 +ry=-0.052 +rz=-3.601 
+s=8.3 +convention=position_vector +step +inv +proj=cart +ellps=GRS80 +step 
+proj=pop +v_3 +step +proj=utm +zone=32 +ellps=GRS80

I would expect there to be something like +nadgrids=kanu_ntv2_oberfranken.gsb.

What could be going wrong here?
Am I missing something? 

Thank you in advance and best regards!
 Thomas





More information about the Qgis-user mailing list