скрипт-slice'инг исходных климатических 2D-полей
-
- Гуру
- Сообщения: 968
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
скрипт-slice'инг исходных климатических 2D-полей
Здравствуйте!
Подскажите, пожалуйста, каким корректным методом и каким инструментом быстро (около 800 NetCDF файлов) организовать скрипт-slice'инг исходных климатических 2D-полей (температура, солёность Мирового океана) в географических координатах разрешением 0.25 град. долготы X 0.16 град. широты в поля с таким же шагом и такими же координатами узлов сетки по долготе, а вот по широте - с другим шагом (0.25 град.). В R или Python, м.б., кто это делал? Какими пакетами, библиотеками?
Подскажите, пожалуйста, каким корректным методом и каким инструментом быстро (около 800 NetCDF файлов) организовать скрипт-slice'инг исходных климатических 2D-полей (температура, солёность Мирового океана) в географических координатах разрешением 0.25 град. долготы X 0.16 град. широты в поля с таким же шагом и такими же координатами узлов сетки по долготе, а вот по широте - с другим шагом (0.25 град.). В R или Python, м.б., кто это делал? Какими пакетами, библиотеками?
-
- Гуру
- Сообщения: 4130
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1096
- Ваше звание: программист
- Откуда: Казань
Re: скрипт-slice'инг исходных климатических 2D-полей
я делал подобное, для этого на R или Python пишется генератор набор команд для утилит gdal, которые решают задачу. Один из подобных наборов у меня содержит около 45000 вызовов разных утилит, работает пару часов.
Сначала много gdal_translate для перевода всех нужных тайлов в GeoTIFF, потом gdalbuildvrt для сборки из них покрытия, потом gdalwarp для вырезания нужного куска в нужной проекции
Сначала много gdal_translate для перевода всех нужных тайлов в GeoTIFF, потом gdalbuildvrt для сборки из них покрытия, потом gdalwarp для вырезания нужного куска в нужной проекции
-
- Гуру
- Сообщения: 968
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: скрипт-slice'инг исходных климатических 2D-полей
Уважаемый gamm, спасибо за декомпозицию дорожной карты! Попробую через GDAL... И именно через GeoTIFF'ы!
-
- Гуру
- Сообщения: 4130
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1096
- Ваше звание: программист
- Откуда: Казань
Re: скрипт-slice'инг исходных климатических 2D-полей
там главное правильно написать имена, для Модиса в hdf4 это выглядит так:
еще помню были утилиты командной строки специально для NetCDF, но gdal вроде его нынче понимает сам.
RTFM
с именами вроде так тыц
Код: Выделить всё
gdal_translate -of GTiff HDF4_EOS:EOS_GRID:"MCD12Q1.A2015001.h24v02.061.2022166050508.hdf":MCD12Q1:LC_Type1 test01.tif
RTFM
с именами вроде так тыц
Код: Выделить всё
gdal_translate -a_srs EPSG:4326 NETCDF:File_Name.nc:Band_Name -of ‘Gtiff’ Output_FileName.geotiff
-
- Гуру
- Сообщения: 968
- Зарегистрирован: 22 май 2010, 20:20
- Репутация: 154
Re: скрипт-slice'инг исходных климатических 2D-полей
gamm, премного благодарен Вам!
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость