Как совместить данные АгроАтласа и WorldClimate?

[bzborisov]

Всем доброго дня, есть проблема. Существуют растровые климатические данные с АгроАтласа (http://www.agroatlas.ru/ru/content/Climatic_maps/) и WorldClimate (http://www.worldclim.org/). Данные очень хорошие, но хочется совместить их, т.е. сделать идентичными. Самая большая проблема это то, что данные АгроАтласа не просто в формате RST, а имеют жесткую проекцию

Проекция: "Равновеликая Альберса на СССР": ref. system : Alber's Equal Area Conic for Russia projection : Alber's Equal Area Conic datum : Pulkovo 1942 (S42) delta WGS84 : 28 -130 -95 ellipsoid : Krasovsky

Самое большее, что получается, это задать эту же проекцию для фрейма ArcGIS, конвертировать RST в GRID, но самое главное не получается провести команду ArcToolBox => Управление данными => Растр => Обработка растра => Изменить разрешение. Программа говорит, что ячейки растра слишком маленькие или слишком большие. Пространственное разрешение данных АгроАтласа 10 км, WorldClimate - 1 км, но даже если разбить растры АгроАтласа на 1 км ячейки, ArcGIS 10.1 все равно пишет эту ошибку. Я понимаю, что это происходит из-за большой разницы в конфигурации самих ячеек, но что делать?

[bolotoved]

Видимо проблема в том, что вы пытаетесь перепроецировать данные выходящие за 180 мередиан.
Во-первых, вам нужно определиться в какой проекции вы будете работать, т.е. в какой проекции находятся наиболее ценные для вас данные. Остальные данные нужно перепроецировать в эту выбранную проекцию.
Проекция, которую избрали авторы агроатласа удобна для работы на территории РФ: можно считать площади (равновеликая) и нет разрыва по 180 (не отрезана Чукотка).
Можно перепроецировать данные WorldClim (которые находятся в географической системе координат) в равновеликую Альберса. При этом, при перепроецировании в Альберса, нужно указать охват - совпадающий с охватом Агроатласа, растр на весь мир в Альберса перепроецировать не удастся (потому как Антарктида устремится в бесконечность). Для таких манипуляций в ArcGIS наверняка есть инструменты и знающий народ подскажет.

Но эта операция может быть выполнена и с помощью открытой утилиты GDAL:

Код: Выделить всё

gdalwarp исходный.tif перепроецированный.tif -s_srs "+proj=longlat +ellps=WGS84" -t_srs "параметры необходимой выходной проекции" -te координаты углов необходимого фрейма -tr вертикальное и горизонтальное разрешение -r способ интерполяции

Чтобы выяснить как правильно записать проекцию данных Агроатласа на языке proj4 (который понимает GDAL) нужно для растра Агроатласа выполнить команду:

Код: Выделить всё

gdalinfo агроатлас.tif

Что такое GDAL и как с ним работать.

[bzborisov]

bolotoved » Сегодня, 01:53
Видимо проблема в том, что вы пытаетесь перепроецировать данные выходящие за 180 мередиан.

О черт!!! Как я забыл!!!! Ну я тупоголовый Спасибо большое, попробую отрезать Чукотку, она в принципе мне не особо нужна в расчетах.