Jak zmienić projekt pliku rastrowego w QGIS z transformacją bazy danych?

11

Chcę ponownie rzutować plik rastrowy Tiff z projekcji A (strefa 32 ETRS89 / UTM) na projekcję B (strefa 2 DHDN / Gaussa-Krugera). Te dwie projekcje mają inny układ odniesienia, więc muszę użyć transformacji układu odniesienia (w moim przypadku NTv2 BETa2007). Chcę utworzyć nowy plik rastrowy z projekcją B, ponieważ chcę go użyć w innym oprogramowaniu. Dlatego nie jestem w stanie potrzebować ponownej projekcji w locie.

Próbowałem Rastra> Przekształć (lub podobnie - nie znam menu w języku angielskim), ale nie widzę tutaj, jak zdefiniować transformację układu odniesienia.

W jaki sposób mogę wykonać trwałe odrzucenie pliku rastrowego w Qgis z transformacją układu odniesienia?

Jens
źródło

Odpowiedzi:

13

Gdalwarp to narzędzie do ponownej projekcji, znajdziesz je w Qgis pod Raster-> Projektionen-> Transformieren lub samodzielne w OSGEO4W.

Podstawowe polecenie to

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 powinien znajdować się w tym samym folderze lub podać bezwzględną ścieżkę do niego.

W Qgis wybierz załadowany tif jako dane wejściowe i jego CRS i wklej następujące elementy w docelowym CRS:

+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
źródło
1
Andre, dzięki za odpowiedź. W przyszłym projekcie muszę przekształcić 300 tiff z UTM do Gaussa-Krugera. Dziś mam tylko tiffy Gaussa-Krugera do przetestowania. Próbowałem 'gdalwarp -s_srs "+ proj = tmerc + lat_0 = 0 + lon_0 = 6 + k = 1 + x_0 = 2500000 + y_0 = 0 + ellps = bessel + datum = poczdam + jednostki = m + no_defs + wktext" -t_srs "+ proj = utm + zone = 32 + ellps = GRS80 + jednostki = m + no_defs + nadgrids = D: \ gis_spielwiese \ raster_reproject \ BETA2007.gsb + wktext" D: \ gis_spielwiese \ raster_reproject \ 4611.tif D: \ gis_spielse raster_reproject \ 4611utm_test7.tif ', ale przesunięcie wynosi 130 m. Więc przypuszczam, że gdalwarp nie dokonuje transformacji układu odniesienia.
Jens
4
To kolejny powód: musisz dodać + towgs84 = 0,0,0,0,0,0,0,0 do ciągu utm. I lepiej napisać + nadgrids w ciągu DHDN. Standardowe kody EPSG zapewniają przesunięcie o około 2 metry.
AndreJ
1
Dzięki nowemu parametrowi działa teraz dobrze. Wielkie dzięki Andre.
Jens
2
Czy byłbyś w stanie opublikować swoje ostateczne polecenie?
Cliff