Як відтворити растровий файл у QGIS з перетворенням даних?


11

Я хочу перепрофілювати растровий файл tiff від проекції A (ETRS89 / UTM Zone 32) до проекції B (DHDN / зона Гаусса-Крюгера 2). Дві проекції мають різну дату, тому мені доведеться використовувати перетворення дат (у моєму випадку NTv2 BETa2007). Я хочу створити новий растровий файл з проекцією B, оскільки я хочу використовувати його в іншому програмному забезпеченні. Тож перепроектування на ходу - це не те, що мені потрібно.

Я спробував Raster> Transform (або подібне - я не знаю англійського меню), але не можу зрозуміти, як тут визначити перетворення дат.

Як я можу зробити постійну репроекцію растрового файлу в Qgis з перетворенням даної?

Відповіді:


13

Gdalwarp - це інструмент для відтворення проектів, ви знайдете його у Qgis у розділі Растр-> Проекція-> Трансформація або автономне в 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 та вставте наступне у цільовий-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 TIF від UTM до Гаусса-Крюгера. Сьогодні у мене є лише прихильники Гаусса-Крюгера, щоб перевірити це. Я спробував 'gdalwarp -s_srs "+ proj = tmerc + lat_0 = 0 + lon_0 = 6 + k = 1 + x_0 = 2500000 + y_0 = 0 + ellps = bessel + date = potsdam + units = m + no_defs + wktext" -t_srs "+ proj = utm + зона = 32 + ellps = GRS80 + одиниці = m + no_defs + nadgrids = D: \ gis_spielwiese \ raster_reproject \ BETA2007.gsb + wktext" D: \ gis_spielwiese \ raster_reproject \ 4611.tif D: \ gis_spiel raster_reproject \ 4611utm_test7.tif ', але є зсув 130 м. Тому я припускаю, що gdalwarp не перетворює дату.
Єнс

4
Це ще одна причина: вам потрібно додати + towgs84 = 0,0,0,0,0,0,0 до рядка utm. І краще написати + надрешітки в рядку DHDN. За допомогою стандартних кодів EPSG ви отримуєте зсув близько 2 метрів.
AndreJ

1
З новим параметром він працює нормально. Дякую Андре.
Єнс

2
Чи зможете ви опублікувати остаточну команду?
Кліф
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.