Comment reprojeter un fichier raster dans QGIS avec une transformation de datum?

11

Je veux reprojeter un fichier raster tiff de la projection A (ETRS89 / UTM Zone 32) à la projection B (DHDN / Gauss-Kruger zone 2). Les deux projections ont un datum différent donc je dois utiliser une transformation de datum (NTv2 BETa2007 dans mon cas). Je veux créer un nouveau fichier raster avec la projection B parce que je veux l'utiliser dans un autre logiciel. Donc, la reprojection à la volée n'est pas ce dont j'ai besoin.

J'ai essayé Raster> Transform (ou similaire - je ne connais pas le menu anglais) mais je ne vois pas comment définir une transformation de datum ici.

Comment faire une reprojection permanente d'un fichier raster dans Qgis avec une transformation de datum?

Jens
la source

Réponses:

13

Gdalwarp est l'outil à reprojeter, vous le trouverez dans Qgis sous Raster-> Projektionen-> Transformieren ou autonome dans OSGEO4W.

La commande de base est

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 doit se trouver dans le même dossier ou utiliser un chemin absolu vers celui-ci.

Dans Qgis, sélectionnez le tif chargé en entrée et son CRS et collez ce qui suit dans le 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
la source
1
Andre, merci pour ta réponse. Dans un futur projet, je dois transformer 300 tiff d'UTM en Gauss-Kruger. Aujourd'hui, je n'ai que des tiffs Gauss-Kruger pour le tester. J'ai essayé '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_spielwies raster_reproject \ 4611utm_test7.tif 'mais il y a un décalage de 130 m. Donc je suppose que gdalwarp ne fait pas de transformation de datum.
Jens
4
C'est une autre raison: vous devez ajouter + towgs84 = 0,0,0,0,0,0,0 à la chaîne utm. Et il vaut mieux écrire les + nadgrids dans la chaîne DHDN. Avec les codes EPSG standard, vous obtenez un décalage d'environ 2 mètres.
AndreJ
1
Avec le nouveau paramètre, cela fonctionne bien maintenant. Merci beaucoup André.
Jens
2
Pourriez-vous publier votre commande finale?
Cliff