データ変換を使用してQGISでラスターファイルを再投影する方法


11

tiffラスターファイルを投影A(ETRS89 / UTMゾーン32)から投影B(DHDN /ガウスクルーガーゾーン2)に再投影したい。2つの投影には異なるデータムがあるため、データム変換を使用する必要があります(私の場合はNTv2 BETa2007)。別のソフトウェアで使用したいので、プロジェクションBで新しいラスターファイルを作成します。そのため、オンザフライでの再投影は必要ありません。

[ラスター]> [変換](または同様の方法-英語のメニューがわからない)を試しましたが、データ変換をここで定義する方法がわかりません。

データム変換を使用してQgisでラスターファイルを永続的に再投影するにはどうすればよいですか?

回答:


13

Gdalwarpは、再投影するツールです。QgisのRaster-> Projektionen-> TransformierenまたはOSGEO4Wのスタンドアロンにあります。

基本的なコマンドは

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は同じフォルダにあるか、または絶対パスを使用する必要があります。

Qgisで、読み込まれたtifとそのCRSを選択し、target-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

1
アンドレ、返信ありがとうございます。今後のプロジェクトでは、300 tiffをUTMからガウスクルーガーに変換する必要があります。今日、私はそれをテストするガウスクルーガーのティフしか持っていません。'gdalwarp -s_srs "+ proj = tmerc + lat_0 = 0 + lon_0 = 6 + k = 1 + x_0 = 2500000 + y_0 = 0 + ellps = bessel + datum = potsdam + units = m + no_defs + wktext" -t_srsを試しました"+ proj = utm + zone = 32 + ellps = GRS80 + units = 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 'ですが、130 mのシフトがあります。したがって、gdalwarpはデータム変換を行わないと思います。
イェンス

4
それがもう1つの理由です。utm文字列に+ towgs84 = 0,0,0,0,0,0,0を追加する必要があります。+ DHDN文字列に+ nadgridsを記述する方が良いでしょう。標準のEPSGコードを使用すると、約2メートルのオフセットが得られます。
AndreJ

1
新しいパラメーターを使用すると、正常に動作します。Andreに感謝します。
イェンス

2
最終的なコマンドを投稿できますか?
クリフ
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.