Я работаю над программой, которая будет использовать GDAL для обработки различных выбранных пользователем глобальных наборов растровых данных с различным разрешением в качестве основы для распределения прогнозируемых выбросов загрязнения в масштабе страны в квадраты сетки.
ПодходЯ пытался использовать это следующим образом:
- Во-первых, обрезать глобальные наборы растровых данных, используя файл формы с многоугольником для данной страны.
- Далее, повторить выборкунаборы данных, чтобы все они имели одинаковое разрешение (размер ячейки) и одинаковые углы XLL / YLL и одинаковое количество строк / столбцов.
Все переделенные наборы данных будут грубее, чем исходные глобальные прокси-данные (большие размеры ячеек), и я хочу использовать метод суммирования.
Мне удалосьвыполнить шаг 1, используя деформацию GDAP, но пока не удалось заставить шаг 2 работать.Я понимаю, что правильный подход к агрегации состоит в том, чтобы использовать метод пересчета «-r среднее», а затем умножать результаты на отношение количества квадратов в исходных и пересчитанных сетках.Но мне не ясно, как создавать сетки с передискретизацией, которые имеют одинаковые строки / столбцы и углы xll / yll.
Может ли кто-нибудь предложить лучшие методы для этого с помощью gdal_warp, gdal_translate или аналогичных?
Примечание: я вызываю эти утилиты из другого приложения - поэтому действительно нужно решение, которое может быть вызвано автоматически (то есть из командной строки) без вмешательства пользователя.Полученные сетки будут отображаться на экране в приложении Windows, написанном на Delphi и использующем MapWinGIS для отображения карт.
Заранее благодарим за любые советы!