Расшифровать Grib-файл

Полезные ссылки:
Метеоцентр.Азия - наш сайт с высокодетализированными прогнозами погоды по пунктам Pоссии и мира
Облегчённая версия Метеоклуба (для смартфонов)

Расшифровать Grib-файл

Сейчас в Метеоклубе:
Участников - 1 [ slava31 ]
Максимальное одновременное количество посетителей: 308 [2 Ноя 2013 22:42]
Гостей - 296 / Участников - 12

 - Начало - Ответить - Статистика - Pегистрация - Поиск -

МЕТЕОКЛУБ : независимое сообщество любителей метеорологии (Европа и Азия) : ФОРУМ О ПОГОДЕ И ПРИРОДЕ / Компьютерная техника и интернет в метеорологии / Расшифровать Grib-файл
<< . 1 . 2 .
Автор Сообщение
lalex
Участник

Москва+Крым
# Дата: 2 Jun 2011 21:02 - Поправил: lalex


Да нет, вроде бы оба Grib2. Раскодировал, понял - который меньше это одноградусный ;)
Вы, я так понимаю, Grib-1 и одноградусную модель предпочитаете?

Corvus
Автор сайта

###Corvus###
# Дата: 2 Jun 2011 21:24


Вы, я так понимаю, Grib-1 и одноградусную модель предпочитаете?
lalex

На Метеоцентре используется 0.5-градусная ГФС в формате ГРИБ1.

lalex
Участник

Москва+Крым
# Дата: 3 Jun 2011 14:23 - Поправил: lalex


Ясно, и с inventory тоже вроде начало проясняться. А насчет выбора конкретных серверов: каталог rotating обновляется быстрее, но в нем заблаговременность только 48 часов, вроде бы?

Corvus
Автор сайта

###Corvus###
# Дата: 3 Jun 2011 16:02


каталог rotating обновляется быстрее, но в нем заблаговременность только 48 часов, вроде бы?
lalex

Нет, почему же.

180 часов.

http://nomad3.ncep.noaa.gov/pub/gfs/rotating/

lalex
Участник

Москва+Крым
# Дата: 3 Jun 2011 19:45


Хм. В этом каталоге вроде бы одноградусные грибы, матрица 360х181 (на пробу распаковал). Размер 32 Мб.
Полуградусные, я так понял, в каталоге http://nomad3.ncep.noaa.gov/pub/gfs/rotating-0.5/ (если на Nomad3). И вот там уже - только по 48 часов (файлы притом 56 Мб и сразу Grib-2).

Corvus
Автор сайта

###Corvus###
# Дата: 3 Jun 2011 20:01


Полуградусные, я так понял, в каталоге http://nomad3.ncep.noaa.gov/pub/gfs/rotating-0.5/ (если на Nomad3). И вот там уже - только по 48 часов (файлы притом 56 Мб и сразу Grib-2).
lalex

Да, точно, я был неправ.

Можно брать 0.5-градусную на 180 ч вот тут
http://nomad1.ncep.noaa.gov/cgi-bin/ftp2u_gfs0.5.s h

Я там беру для Метеоцентра данные ГФС. Но система не шибко надёжная, от 1 до 5 дней в месяц данные не обновляются.

В общем, лучше берите СМС - там всё просто и стабильно :)

lalex
Участник

Москва+Крым
# Дата: 3 Jun 2011 21:52 - Поправил: lalex


"от 1 до 5 дней в месяц данные не обновляются"

Я как раз хотел об этом спросить... Наверно, по этой причине на Метеоцентре висит прогноз за 01.06?
С другой стороны, на основном сервере Nomads вроде бы все давно уже обновилось, по третье июня включительно, прогон за 12 засветился недавно. И заблаговременность полноценная - 192 часа... Ну в общем к этому всему надо пристреляться самостоятельно, как я понял.
http://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/ prod/

насчет СМС внимательно читаю все отзывы, слышал, думаю. Как и насчет других моделей. Но это отдельная тема, пока что.

Corvus
Автор сайта

###Corvus###
# Дата: 3 Jun 2011 21:57


насчет СМС внимательно читаю все отзывы, слышал, думаю
lalex

Основные копья ломаются относительно кол-ва конвективных осадков по СМС (оно часто завышено). Если брать из СМС обычные осадки, то будет примерно то же самое, что и по ГФС.

Правда, при этом 30% гроз и ливней будут просто-напросто неспрогнозированы (это присуще ГФС), но для многих метеоклубовцев это меньшее зло, чем перестраховочные прогнозы СМС вкупе с 90% предупреждённостью опасных конв. явлений (у ГФС эта предупреждённость 70%).

gasha
Участник

Россия, Карелия, Петрозаводск
# Дата: 14 Jun 2011 21:37 - Поправил: gasha


Отсюда скачиваю 65 grib2 файлов (0-192 часа). Затем в батнике командой copy /B склеиваю их в один большой гриб(gfs.t00z.mastergrb2f). По размеру он как исходные 65 грибов (gfs.t00z.mastergrb2f00 и т.д.).

После идёт следующий батник:

perl g2ctl.pl gfs.t00z.mastergrb2f >gfs.t00z.mastergrb2f.ctl
gribmap -i gfs.t00z.mastergrb2f.ctl

Так вот в получившемся ctl файле 61 момент времени, а не 65.

dtype grib2
ydef 41 linear 50.000000 0.5
xdef 101 linear 0.000000 0.500000
tdef 61 linear 00Z14jun2011 3hr

Версия g2ctl.pl у меня 0.0.8.2

wgrib2 версии 0.1.9.1
gribmap (This is gribmap for GrADS Version 2.0.a9.oga.1)

Куда пропадают 4 момента времени (183, 186, 189, 192 часа)?

Вот пример скрипта (tmp2m.gs):

V=tmp2m
V2=pressfc
hour=0

'set mpdset hires'
'set map auto'
'set mproj scaled'
'set lat 50 70'
'set lon 0 50'

'open gfs.t00z.mastergrb2f.ctl'

timex=1

while( timex <= 65)
'set t ' timex

'value='V'-273.16'
'value2='V2''

'set gxout shaded'
'd value'
'run cbarn.gs'
'set gxout contour'
'set ccolor 0'
'd value2'

'q time'
times = sublin(result,1)
timeword = subwrd(times,3)
'draw title 'V' / 'V2', valid: 'timeword''

'printim 'V'_'hour'.png white x800 y800'

timex = timex + 1
hour=hour+3
'clear'
endwhile

Генерируемые файлы от tmp2m_0.png до tmp2m_180.png содержат информация, с ними всё хорошо. А вот оставшиеся с tmp2m_183.png по tmp2m_192.png пустые.
grads пишет:

Data Request Warning: Request is completely outside file limits
Data Request Warning: Request is completely outside file limits
Cannot contour grid - all undefined values
Cannot plot color bar: No shading information
Cannot contour grid - all undefined values

и таких 4 пары сообщений.
Что доказывает, что g2ctl.pl проигнорировал 4 момента времени (183-192). Где в нём чего надо исправить, чтобы заработало. Сижу в никаю в этот скрипт, пока не нашёл.

P.S. Вот изменённый скрипт meteogram_gfs.gs, данные для построения беруться из большого гриба) И понятное дело, что строиться метеограмма по 61 моменту, а не по 65 как хотелось бы. Так как после g2ctl.pl получается не совсем правильный ctl

P.P.S. ОТБОЙ ТРЕВОГИ!
Всё хорошо)) Ошибка была при скачивании грибов: как раз последние 4 из 65 качались с параметрами 61-ого. Теперь все проблемы исчезли))

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 5 Май 2014 13:23


Как декодировать GRIB1 или GRIB2 без использования сторонних программ? Реализовать, например на С++ , в интернете всё облазил, не нашёл схему декодирования...

Corvus
Автор сайта

###Corvus###
# Дата: 5 Май 2014 17:39 - Поправил: Corvus


Mesocyclon

Зачем реализовывать это всё с нуля, если есть прекрасные консольные программы (под Винду и Юникс) wgrib и wgrib2?

Наверняка и их исходные коды есть в инете (на C), их можно встроить в любые свои приложения (хотя, ИМХО, гораздо проще запускать wgrib и wgrib2 в виде отдельных задач по планировщику, а потом своими программами "кушать" уже расшифрованные ГРИБы).

Corvus
Автор сайта

###Corvus###
# Дата: 5 Май 2014 17:40


Кроме того, в инете есть подробное описание кода ГРИБ. Я 12 лет назад писал на Дельфи самопальную программку для декодирования ГРИБ.

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 5 Май 2014 18:36 - Поправил: qwerty


Скачал гриб ukmet

Отсюда
http://nomads.ncep.noaa.gov/pub/data/nccf/com/ukme t/prod/ukmet.201405/

Создал файл ukm.ctl

Пытался открыть в Градсе в итоге он мне выдал ошибку Can't open description file

В чём дело?

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 5 Май 2014 19:18


Corvus

Прочитал что в формате GRIB2 использован алгоритм сжатия JPEG2000, спрашивается зачем было использовать видеоформат.

http://www.gdal.org/frmt_grib.html

либо я не так понял?

Corvus
Автор сайта

###Corvus###
# Дата: 6 Май 2014 11:41


Mesocyclon

Не забивайте себе голову лишней информацией, лучше изучите программы ВГРИБ и ВГРИБ2, там и исходные коды есть

http://www.cpc.ncep.noaa.gov/products/wesley/wgrib .html
http://www.cpc.ncep.noaa.gov/products/wesley/wgrib 2/

TornadoF5
Участник

Харьков, Украина. (Игорь)
# Дата: 6 Май 2014 13:36


Corvus
Проверьте пожалуйста почту, я Вам письмо отправлял дня 2 назад :)

Corvus
Автор сайта

###Corvus###
# Дата: 6 Май 2014 16:24


TornadoF5

Я видел, при наличии времени отвечу на него.

met2
Участник

Минск
# Дата: 8 Май 2014 15:15


Corvus
TornadoF5

ну или degrib. она еще проще

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 17 Jul 2014 14:57 - Поправил: qwerty


Почему при открытие ctl файла выскакивает ошибка Can't open description file ?

qwerty
Участник

Владимирская область, г.Киржач
# Дата: 16 Авг 2014 15:13


Как открыть grib файл?

Mesocyclon
Участник

Троицк, Ключевка, Челябинская область
# Дата: 16 Авг 2014 16:55


qwerty

попробуйте этой программой

www.giss.nasa.gov/tools/panoply/download_win.html

Meteoman16
Участник

Белогорск, Крым
# Дата: 17 Янв 2015 14:10


чем можно открыть GRIB2? уже пытался разными программами, не получается,может как-то можно конвертировать в GRIB1, его я могу открыть

matrix145
Участник

Курган
# Дата: 17 Янв 2015 14:30


Meteoman16

Файл grb2 легко открывается программой panoplay, выше уже кидали ссылку.
http://www.giss.nasa.gov/tools/panoply/download_wi n.html

Meteoman16
Участник

Белогорск, Крым
# Дата: 17 Янв 2015 14:58


matrix145 скачал, однако открыть не удалось, тормозит

matrix145
Участник

Курган
# Дата: 17 Янв 2015 15:19


matrix145 скачал, однако открыть не удалось, тормозит

Я хорошо открываю через него файлы, даже большого размера.

Meteoman16
Участник

Белогорск, Крым
# Дата: 17 Янв 2015 19:12


matrix145

я так понимаю, что работает он с JAva. У меня открывается окно после загрузки пустое, так и стоит

Corvus
Автор сайта

###Corvus###
# Дата: 18 Янв 2015 01:28


Meteoman16

Нужно установить или обновить Java на Вашем компьютере.
В Файрфоксе это делается через закладку "Плагины".

<< . 1 . 2 .
Ваш ответ

          Отменить *Что это?

 » Логин  » Пароль 
 
 


Поддержка: miniBB forum software © 2001-2024