QGIS'te bir raster dosyasını datum dönüşümü ile nasıl yeniden kopyalayabilirim?

11

Bir tiff tarama dosyasını projeksiyon A'dan (ETRS89 / UTM Bölge 32) projeksiyon B'ye (DHDN / Gauss-Kruger bölge 2) yeniden yansıtmak istiyorum. İki projeksiyonun farklı bir referans noktası vardır, bu nedenle bir referans dönüşümü kullanmak zorundayım (benim durumumda NTv2 BETa2007). B projeksiyonu ile yeni bir tarama dosyası oluşturmak istiyorum çünkü başka bir yazılımda kullanmak istiyorum. Bu yüzden anında yeniden projeksiyon ihtiyacım olan şey değil.

Raster> Transform (veya benzer - İngilizce menüsünü bilmiyorum) denedim ama burada bir datum dönüşümü tanımlamak nasıl göremiyorum.

Qgis'te datum dönüşümü olan bir raster dosyasının kalıcı olarak yeniden projeksiyonunu nasıl yapabilirim?

qgis  raster  coordinate-system  datum 
— Jens
kaynak

Yanıtlar:

13

Gdalwarp yeniden projelendirme aracıdır, Raster-> Projektionen-> Transformieren veya OSGEO4W'da bağımsız olarak Qgis'de bulabilirsiniz.

Temel komut

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 aynı klasörde olmalı veya mutlak yol kullanmalıdır.

Qgis'te, yüklenen tif'i girdi olarak ve CRS'sini seçin ve aşağıdakileri hedef-CRS'ye yapıştırın:

+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
kaynak
1
Andre, cevabın için teşekkürler. Gelecekteki bir projede 300 tiff'i UTM'den Gauss-Kruger'a dönüştürmem gerekiyor. Bugün sadece Gauss-Kruger test etmek için tiff'lerim var. 'Gdalwarp -s_srs "denedim + proj = tmerc + lat_0 = 0 + lon_0 = 6 + k = 1 + x_0 = 2500000 + y_0 = 0 + ellps = bessel + datum = potsdam + birimleri = m + no_defs + wktext" -t_srs "+ proj = utm + bölge = 32 + ellps = GRS80 + birim = m + no_defs + nadgrids = D: \ gis_spielwiese \ raster_reproject \ BETA2007.gsb + wktext" D: \ gis_spielwiese \ raster_reproject \ 4611.tif D: \ gis_spielwies raster_reproject \ 4611utm_test7.tif 'ama 130 metrelik bir kayma var. Bu yüzden sanırım gdalwarp veri dönüşümü yapmadı.
— Jens
4
Bu başka bir neden: utm dizesine + towgs84 = 0,0,0,0,0,0,0,0 eklemelisiniz. Ve + nadgrid'leri DHDN dizesine yazmak daha iyidir. Standart EPSG kodlarıyla yaklaşık 2 metre ofset elde edersiniz.
— AndreJ
1
Yeni parametre ile şimdi iyi çalışıyor. Çok teşekkürler Andre.
— Jens
2
Son komutunuzu gönderebilir misiniz?
— Cliff
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.