GRIB2-файл -прочитать и записать в мЕньший по размеру домейн

[nickleb]

Здравствуйте. GRIB2-файл -прочитать и записать в мЕньший по размеру домейн (все band'-"ы" надо из исходного файла), вероятно, только через GDAL? "Дружелюбные" wrapper'-"ы" для этой операции в свободных ГИС, м.б., есть? Подскажите, пожалуйста, - кто с этим плотно работал. Пробовал в R через package rNOMADS - с ходу не пошло...

[gamm]

Подскажите, пожалуйста

ничего не понял - я эти сетки когда-то (лет ...дцать назад) "руками" читал (на С++), там вроде обычная сетка без особых затей, даже не HDF. Если формат посмотреть, то наверняка можно, используя двоичный ввод/вывод, на R собрать читалку, я так гигантские SURFER Binary Grid читал по частям, типа тайлами ...

[nickleb]

Подскажите, пожалуйста

ничего не понял - я эти сетки когда-то (лет ...дцать назад) "руками" читал (на С++), там вроде обычная сетка без особых затей, даже не HDF. Если формат посмотреть, то наверняка можно, используя двоичный ввод/вывод, на R собрать читалку, я так гигантские SURFER Binary Grid читал по частям, типа тайлами ...

gamm, здравствуйте! надо из большого файла формата GRIB2 с сонмом характеристик разного рода в регулярной сеточной области по всему земному шару выбрать в широтной зоне, скажем, ревущих сороковых subset всех band'ов этого исходного файла и записать это в меньший по размеру файл такого же формата GRIB2... и всё...

[gamm]

прочитать еще можно (наверное), если попотеть, а вот записать ... я вообще не уверен, что фрагмент туда можно записать, а не весь, чтобы его кто-то потом прочитал. Вообще, поскольку вся эта канитель исходно была сделана под Унихом, и там все работает лучше, чем порты под Уиндовс, может есть смысл поставить виртуалку с Убунту, и там все конвертировать.

[nickleb]

"... — Понятно… Разрешите хотя бы…"
www.cpc.ncep.noaa.gov/products/wesley/w ... _grib.html
"— А вот это попробуйте! Вам поручена эта операция, так что действуйте..."
[ "Бриллиантовая рука" ]

[bolotoved]

Помучившись с GDAL, пришел к выводу, что для этого специфического формата лучше использовать спец. утилиты. Я своих целей добился через wgrib и wgrib2. Но в grib сохранять я не пытался, т.ч. не знаю, умеет ли он это делать.

[nickleb]

Помучившись с GDAL, пришел к выводу, что для этого специфического формата лучше использовать спец. утилиты. Я своих целей добился через wgrib и wgrib2. Но в grib сохранять я не пытался, т.ч. не знаю, умеет ли он это делать.

Уважаемый bolotoved, спасибо Вам за ценную информацию!
"Опыт бесценен, плохо только, что за него приходится платить собственной молодостью" [С. Харви]

[nickleb]

Материал по ссылке
https://rda.ucar.edu/datasets/ds083.2/s ... stall+.txt
пригодился для того, чтобы в Windows 10 после инсталляции wgrib2 сделать желаемую операцию...
Вот так в командной строке из "большого" grib2 (имя Source.grb2) выбираем subset "маленького" grib2 в широтной зоне от 50 ю.ш. до 40 ю.ш. Этот subset записывается в файл с именем from_Source_Fi_50S_40S.grb2:

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

wgrib2 Source.grb2 -set_grib_type j -small_grib 0:360 -50:-40 from_Source_Fi_50S_40S.grb2