скрипт-slice'инг исходных климатических 2D-полей

Вопросы общего характера по ГИС и дистанционному зондированию, не связанные с конкретным ПО.
Ответить
nickleb
Гуру
Сообщения: 968
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

скрипт-slice'инг исходных климатических 2D-полей

Сообщение nickleb » 14 ноя 2024, 11:03

Здравствуйте!
Подскажите, пожалуйста, каким корректным методом и каким инструментом быстро (около 800 NetCDF файлов) организовать скрипт-slice'инг исходных климатических 2D-полей (температура, солёность Мирового океана) в географических координатах разрешением 0.25 град. долготы X 0.16 град. широты в поля с таким же шагом и такими же координатами узлов сетки по долготе, а вот по широте - с другим шагом (0.25 град.). В R или Python, м.б., кто это делал? Какими пакетами, библиотеками?

gamm
Гуру
Сообщения: 4130
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1096
Ваше звание: программист
Откуда: Казань

Re: скрипт-slice'инг исходных климатических 2D-полей

Сообщение gamm » 14 ноя 2024, 15:27

я делал подобное, для этого на R или Python пишется генератор набор команд для утилит gdal, которые решают задачу. Один из подобных наборов у меня содержит около 45000 вызовов разных утилит, работает пару часов.
Сначала много gdal_translate для перевода всех нужных тайлов в GeoTIFF, потом gdalbuildvrt для сборки из них покрытия, потом gdalwarp для вырезания нужного куска в нужной проекции

nickleb
Гуру
Сообщения: 968
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Re: скрипт-slice'инг исходных климатических 2D-полей

Сообщение nickleb » 14 ноя 2024, 15:58

Уважаемый gamm, спасибо за декомпозицию дорожной карты! Попробую через GDAL... И именно через GeoTIFF'ы!

gamm
Гуру
Сообщения: 4130
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 1096
Ваше звание: программист
Откуда: Казань

Re: скрипт-slice'инг исходных климатических 2D-полей

Сообщение gamm » 14 ноя 2024, 17:33

там главное правильно написать имена, для Модиса в hdf4 это выглядит так:

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

gdal_translate -of GTiff HDF4_EOS:EOS_GRID:"MCD12Q1.A2015001.h24v02.061.2022166050508.hdf":MCD12Q1:LC_Type1 test01.tif 
еще помню были утилиты командной строки специально для NetCDF, но gdal вроде его нынче понимает сам.
RTFM
с именами вроде так тыц

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

 gdal_translate -a_srs EPSG:4326 NETCDF:File_Name.nc:Band_Name  -of ‘Gtiff’ Output_FileName.geotiff

nickleb
Гуру
Сообщения: 968
Зарегистрирован: 22 май 2010, 20:20
Репутация: 154

Re: скрипт-slice'инг исходных климатических 2D-полей

Сообщение nickleb » 14 ноя 2024, 17:50

gamm, премного благодарен Вам!

Ответить

Вернуться в «Общие вопросы»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость