Bagaimana cara memproyeksi ulang file raster di QGIS dengan transformasi datum?

11

Saya ingin memproyeksi ulang file tiff raster dari proyeksi A (ETRS89 / UTM Zone 32) ke proyeksi B (DHDN / Gauss-Kruger zone 2). Kedua proyeksi memiliki datum yang berbeda jadi saya harus menggunakan transformasi datum (NTv2 BETa2007 dalam kasus saya). Saya ingin membuat file raster baru dengan proyeksi B karena saya ingin menggunakannya dalam perangkat lunak lain. Jadi proyeksi ulang dengan cepat bukan yang saya butuhkan.

Saya telah mencoba Raster> Transform (atau serupa - saya tidak tahu menu bahasa Inggris) tapi saya tidak bisa melihat cara mendefinisikan transformasi datum di sini.

Bagaimana saya bisa melakukan proyeksi ulang permanen file raster di Qgis dengan transformasi datum?

qgis raster coordinate-system datum Jens
sumber

Jawaban:

13

Gdalwarp adalah alat untuk memproyeksi ulang, Anda menemukannya di Qgis di bawah Raster-> Projektionen-> Transformieren atau standalone di OSGEO4W.

Perintah dasar adalah

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 harus berada di folder yang sama, atau gunakan path absolut untuk itu.

Di Qgis, pilih tif yang dimasukkan sebagai input, dan CRS-nya, dan rekatkan yang berikut ini di 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
AndreJ
sumber
1
Andre, terima kasih atas balasan Anda. Dalam proyek masa depan saya harus mengubah 300 tiff dari UTM ke Gauss-Kruger. Hari ini saya hanya punya tiffs Gauss-Kruger untuk mengujinya. Saya mencoba '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 + unit = m + no_defs + nadgrids = D: \ gis_spielwiese \ raster_reproject \ BETA2007.gsb + wktext" D: \ gis_spielwiese \ raster_reproject \ gis raster_reproject \ 4611utm_test7.tif 'tetapi ada pergeseran 130 m. Jadi saya kira gdalwarp tidak melakukan transformasi datum.
Jens
4
Itulah alasan lain: Anda harus menambahkan + towgs84 = 0,0,0,0,0,0,0 ke string utm. Dan lebih baik untuk menulis + nadgrids di string DHDN. Dengan kode EPSG standar Anda mendapatkan offset sekitar 2 meter.
AndreJ
1
Dengan parameter baru itu berfungsi dengan baik sekarang. Terima kasih banyak, Andre.
Jens
2
Apakah Anda dapat memposting perintah terakhir Anda?
Tebing