Страница 1 из 1

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

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

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

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

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

Добавлено: 14 ноя 2024, 15:58
nickleb
Уважаемый gamm, спасибо за декомпозицию дорожной карты! Попробую через GDAL... И именно через GeoTIFF'ы!

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

Добавлено: 14 ноя 2024, 17:33
gamm
там главное правильно написать имена, для Модиса в 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

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

Добавлено: 14 ноя 2024, 17:50
nickleb
gamm, премного благодарен Вам!