Как перепроектировать растровый файл в QGIS с преобразованием данных?

11

Я хочу перепроектировать файл растровых изображений из проекции A (ETRS89 / UTM Zone 32) в проекцию B (DHDN / Gauss-Kruger zone 2). Две проекции имеют разные данные, поэтому я должен использовать преобразование данных (NTv2 BETa2007 в моем случае). Я хочу создать новый растровый файл с проекцией B, потому что я хочу использовать его в другом программном обеспечении. Поэтому репроекция на лету - это не то, что мне нужно.

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

Как я могу сделать постоянное перепроецирование растрового файла в Qgis с преобразованием данных?

Jens
источник

Ответы:

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, и вставьте следующее в целевой 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
источник
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 + unit = m + no_defs + nadgrids = D: \ gis_spielwiese \ raster_reproject \ BETA2007.gsb + wktext" D: \ gis_spielwiese \ raster_reproject \ 4611.tif D: \ gse raster_reproject \ 4611utm_test7.tif ', но есть смещение 130 м. Так что я полагаю, что gdalwarp не выполняет преобразование данных.
Йенс
4
Это еще одна причина: вы должны добавить + towgs84 = 0,0,0,0,0,0,0 к строке utm. И лучше написать + nadgrids в строке DHDN. Со стандартными кодами EPSG вы получаете смещение около 2 метров.
AndreJ
1
С новым параметром теперь все работает нормально. Большое спасибо, Андре.
Йенс
2
Сможете ли вы опубликовать свою последнюю команду?
Клифф