Wie projiziere ich eine Rasterdatei in QGIS mit Bezugstransformation neu?

11

Ich möchte eine TIFF-Rasterdatei von Projektion A (ETRS89 / UTM Zone 32) auf Projektion B (DHDN / Gauss-Kruger Zone 2) neu projizieren. Die beiden Projektionen haben ein unterschiedliches Datum, daher muss ich eine Bezugstransformation verwenden (in meinem Fall NTv2 BETa2007). Ich möchte eine neue Rasterdatei mit Projektion B erstellen, da ich sie in einer anderen Software verwenden möchte. Eine schnelle Projektion ist also nicht das, was ich brauche.

Ich habe versucht, Raster> Transformieren (oder ähnliches - ich kenne das englische Menü nicht), aber ich kann hier nicht sehen, wie eine Bezugstransformation definiert wird.

Wie kann ich eine permanente Reprojektion einer Rasterdatei in Qgis mit einer Bezugstransformation durchführen?

Jens
quelle

Antworten:

13

Gdalwarp ist das Tool zum Neuprojektieren. Sie finden es in Qgis unter Raster-> Projektionen-> Transformieren oder eigenständig in OSGEO4W.

Grundbefehl ist

gdalwarp -s_srs EPSG:25832 -t_srs "+proj=tmerc +lat_0=0 +lon_0=6 +x_0=2500000 +y_0=0 +k=1.000000 +ellps=bessel +units=m +nadgrids=./BETA2007.gsb +wktext" input.tif output.tif

BETA2007.gsb sollte sich im selben Ordner befinden oder den absoluten Pfad dazu verwenden.

Wählen Sie in Qgis das geladene tif als Eingabe und dessen CRS aus und fügen Sie Folgendes in das Ziel-CRS ein:

+proj=tmerc +lat_0=0 +lon_0=6 +x_0=2500000 +y_0=0 +k=1.000000 +ellps=bessel +units=m +nadgrids=D:\path\to\your\BETA2007.gsb +wktext
AndreJ
quelle
1
Andre, danke für deine Antwort. In einem zukünftigen Projekt muss ich 300 tiff von UTM auf Gauss-Kruger umstellen. Heute habe ich nur Gauß-Krüger-Tiffs, um es zu testen. Ich habe versucht 'gdalwarp -s_srs "+ proj = tmerc + lat_0 = 0 + lon_0 = 6 + k = 1 + x_0 = 2500000 + y_0 = 0 + ellps = bessel + datum = potsdam + unit = m + no_defs + wktext" -t_srs "+ proj = utm + zone = 32 + ellps = GRS80 + Einheiten = m + no_defs + nadgrids = D: \ gis_spielwiese \ raster_reproject \ BETA2007.gsb + wktext" D: \ gis_spielwiese \ raster_reproject \ 4611.tif D: \ gis_spielwiese \ raster_reproject \ 4611utm_test7.tif ', aber es gibt eine Verschiebung von 130 m. Ich nehme also an, dass gdalwarp keine Bezugstransformation durchführt.
Jens
4
Das ist ein weiterer Grund: Sie müssen + utgs84 = 0,0,0,0,0,0,0 zur utm-Zeichenfolge hinzufügen. Und es ist besser, die + nadgrids in den DHDN-String zu schreiben. Mit Standard-EPSG-Codes erhalten Sie einen Versatz von ca. 2 Metern.
AndreJ
1
Mit dem neuen Parameter funktioniert es jetzt einwandfrei. Vielen Dank Andre.
Jens
2
Könnten Sie Ihren letzten Befehl veröffentlichen?
Cliff