Откуда был получен этот 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 КБ) 4995 просмотров