Ресамплировать изображения GeoTIFF в том же разрешении?

21

У меня есть два геофизических изображения, которые я хотел бы изменить в том же разрешении. Если бы я использовал ArcGIS, это было бы легко сделать. Но используя QGIS, я не знаю, как это сделать. Некоторая помощь будет очень признательна. Благодарю.

Гил Гонсалвес
источник

Ответы:

29

В QGIS это тоже легко, хотя и немного менее очевидно. Есть несколько способов сделать это:

  • Растровый калькулятор - просто используйте растровый калькулятор, и вы можете установить разрешение и экстент там, чтобы они соответствовали другому растру, выбрав растровую полосу, которой вы хотите соответствовать, в списке Растровые полосы и затем нажав кнопку «Текущий экстент слоя». Поля столбцов и строк позволят вам установить разрешение. Однако этот метод не дает вам контроля над методом повторной выборки.
  • Используя GDAL_Warp - этот инструмент позволяет вам установить выходное разрешение, указав ширину и высоту выходного растра или указав ключ -tr (см. Документацию). Вы можете перейти к инструменту GDAL_warp, выбрав Raster-> Projection-> Warp (я говорил, что это не было очевидно с точки зрения повторной выборки!).

    • (v2.x) Если вы хотите использовать ключ -tr, заполните все поля для входного растра и выходных данных и т. д. (ваши значения исходного и целевого SRS, вероятно, будут одинаковыми в этом случае - хотя не обязательно, если ты тоже проектируешь) Затем нажмите маленький значок карандаша внизу и отредактируйте автоматически сгенерированную командную строку gdal-warp, чтобы включить в нее ключ -tr. Gdal_wrap позволяет вам указать алгоритм, который вы хотите использовать для повторной выборки, и поэтому он немного менее туп, чем использование растрового калькулятора.
    • (v3.x) Ключ -tr включается с помощью разрешения Выходной файл в целевых ячейках с географической привязкой . Например, чтобы уменьшить частоту дискретизации 1 м до 2 м, в этом поле можно ввести 2. Однако нет возможности передать два разных аргумента для неквадратных пикселей. Скажите, что ваш целевой размер пикселя 0.3125,0.25, то есть xresесть 0.3125и yresесть 0.25. Если вы сейчас передадите значение 0.3125в это поле, оно будет установлено -tr 0.3125 0.3125в команде. Чтобы снять это ограничение, просто скопируйте код, вставьте в командную строку, отредактируйте флаг -tr и запустите. Например:

      gdalwarp -t_srs EPSG:4326 -tr 0.3125 0.25 -r near -te 71.40625 24.875 84.21875 34.375 -te_srs EPSG:4326 -of GTiff foo.tiff bar.tiff

      (в зависимости от вашей установки и переменных среды вам также может понадобиться указать путь к gdalwarp).

MappaGnosis
источник
1
2-й вариант работает только в QGIS v. 2.x, редактирование команд GDAL недоступно в v. 3x questions.qgis.org/issues/15090
reima
Правда, к сожалению! Однако вы можете использовать gdal_warp из командной строки. С установленным QGIS у вас будет доступ к этому, поэтому, хотя неудобно, второй вариант все еще доступен.
MappaGnosis
3

Я обычно использую инструмент Raster Alignment.

Он локализован на вкладке растра и может быть использован для установки растров в той же степени. Он также может выполнять некоторые операции при повторной выборке с новым разрешением (например, среднее, минимальное и максимальное значения и т. Д.).

Вот ссылка с некоторой информацией об этом: https://docs.qgis.org/2.18/en/docs/user_manual/working_with_raster/raster_analysis.html?highlight=raster%20alignment#id3

Я использую версию 2.18.13. Я не знаю, доступен ли инструмент в более старых версиях.

Габриэль Гуарилья Перес
источник
1
Похоже, это не работает для пересэмплирования - уменьшения разрешения растра.
Якуб Сисак GeoGraphics
странно, у меня сработало в прошлый раз, когда я использовал. Это может быть немного глючит, хотя.
Габриэль Гуарилья Перес
Я использовал Warp (reproject), и он работал как шарм с агрегатом и другими опциями ...
Jakub Sisak GeoGraphics