GRIB2 import

[dime1]

Здравствуйте,
довольно срочно нужно научиться импортировать такие гриб-файлы в QGIS 2.4.
Был когда-то опыт импорта netcdf в ArcGIS, потому имею представление, что скорее всего проблема в том, что нет стандартизированной структуры этих файла.
Пробовал zygrib, но он сказал, что у меня не grib-файл.

Пример во вложении.

[Игорь Белов]

Программа Degrib показывает две записи в структуре файла. Так что это, видимо, контейнер данных GRIB. Однако сами записи не читаются: Degrib жалуется на неизвестную координатную систему.
GrADS требует файл описания данных с расширением ctl.

[Александр Мурый]

Откуда был получен этот GRIB-файл?

Вот вывод gdalinfo для файла:

Спойлер

Driver: GRIB/GRIdded Binary (.grb)
Files: h05_20140717_0900_24_rom.grb
h05_20140717_0900_24_rom.grb.aux.xml
Size is 1900, 900
Coordinate System is:
PROJCS["unnamed",
GEOGCS["Coordinate System imported from GRIB file",
DATUM["unknown",
SPHEROID["Spheroid imported from GRIB file",6378140,298.2529810614887]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Geostationary_Satellite"],
PARAMETER["central_meridian",0],
PARAMETER["satellite_height",35785831],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0]]
Origin = (-5571679.196145262569189,5574935.774084444157779)
Pixel Size = (5861.840290526315584,-12374.996168888888860)
Corner Coordinates:
Upper Left (-5571679.196, 5574935.774)
Lower Left (-5571679.196,-5562560.778)
Upper Right ( 5565817.356, 5574935.774)
Lower Right ( 5565817.356,-5562560.778)
Center ( -2930.920, 6187.498) ( 0d 1'34.78"W, 0d 3'21.45"N)
Band 1 Block=1900x1 Type=Float64, ColorInterp=Undefined
Description = 0.000[-] undefined ()
Min=0.025 Max=59.275
Minimum=0.025, Maximum=59.275, Mean=3.040, StdDev=5.280
NoData Value=9999
Metadata:
GRIB_COMMENT=Estimated precipitation [kg/(m^2)]
GRIB_ELEMENT=
GRIB_FORECAST_SECONDS=0 sec
GRIB_PDS_PDTN=30
GRIB_PDS_TEMPLATE_NUMBERS=1 0 8 54 1 0 0 0 0 0 255 255 255 255 255
GRIB_REF_TIME=1405587600 sec UTC
GRIB_SHORT_NAME=0 undefined
GRIB_UNIT=[kg/(m^2)]
GRIB_VALID_TIME=1405587600 sec UTC
STATISTICS_MAXIMUM=59.275001525879
STATISTICS_MEAN=3.0404158207925
STATISTICS_MINIMUM=0.025000035762787
STATISTICS_STDDEV=5.2801207938831
Band 2 Block=1900x1 Type=Float64, ColorInterp=Undefined
Description = 0.000[-] undefined ()
Min=73.000 Max=73.000
Minimum=73.000, Maximum=73.000, Mean=73.000, StdDev=0.000
NoData Value=9999
Metadata:
GRIB_COMMENT=(prodType 3, cat 0, subcat 192) [-]
GRIB_ELEMENT=unknown
GRIB_FORECAST_SECONDS=0 sec
GRIB_REF_TIME=1405587600 sec UTC
GRIB_SHORT_NAME=0 undefined
GRIB_UNIT=[-]
GRIB_VALID_TIME=1405587600 sec UTC
STATISTICS_MAXIMUM=73
STATISTICS_MEAN=73
STATISTICS_MINIMUM=73
STATISTICS_STDDEV=0

Из вывода gdalinfo видно, что не указан датум и сфероид; канала в растре два (причём только первый содержит полезную информацию, Estimated precipitation [kg/(m^2)]).

Простая команда

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

gdalwarp h05_20140717_0900_24_rom.grb h05_20140717_0900_24_rom.tif

делает двухканальный геотифф, который можно скормить QGIS. Правда, лично мне не очень ясно, как использовать такие данные в проекции "Geostationary_Satellite".

qgis_grib2.png (66.4 КБ) 4998 просмотров

[dime1]

Откуда был получен этот GRIB-файл?

Это EUMETSATовские продукты H-SAF H05 (accumulated precipitation). По-моему они просто использовали WGS_84. Т.е. получается, что если бы они напрямую в структуре файла указывали: PROJCS, GEOGCS, DATUM, SPHEROID, то проблем с чтение .grib и переводом в другой формат не возникало бы?
Еще у них на FTP нашёл декодеры информации (во вложении).

Спасибо за помощь)
Мне уже не к спеху, тренинг завершился. Теперь буду понемногу сам разбираться.