ECMWF | Advancing global NWP through international collaboration后使用快捷导航没有帐号？
只需一步, 快速开始
只需一步，快速开始
查看: 9690|回复: 26
关于欧洲中心GRIB数据格式的解码问题
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
之前国家气象局下发的是每个预报时次一个文件，文件大小都是一样的，49M左右，分辨率为0.25x0.25,361x281个数据，假设文件名为ec_old.bin，用wgrib -V ec_old.bin还可以理解，解码如下：
rec 1:153480:date
T kpds5=130 kpds6=100 kpds7=100 levels=(0,100) grid=255 100 mb 12hr fcst:
&&T=Temperature [K]
&&timerange 0 P1 12 P2 0 TimeU 1&&nx 361 ny 281 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat
60.000000 to -10.000000 by 0.250000&&nxny 101441
& && && & long 60.000000 to 150.000000 by 0.250000, (361 x 281) scan 0 mode 128 bdsgrid 1
&&min/max data 184.643 229.439&&num bits 12&&BDS_Ref 184.643&&DecScale 0 BinScale -6
rec 2:305760:date
FAL kpds5=243 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc 12hr fcst:
&&FAL=Forecast albedo [(0 - 1)]
&&timerange 0 P1 12 P2 0 TimeU 1&&nx 361 ny 281 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat&&60.000000 to -10.000000 by 0.250000&&nxny 101441
& && && & long 60.000000 to 150.000000 by 0.250000, (361 x 281) scan 0 mode 128 bdsgrid 1
&&min/max data 0.06 0.864688&&num bits 8&&BDS_Ref 0.06&&DecScale 0 BinScale -8
但最近国家局下发的格式变了，每个时次的文件大小不一致了，但文件命名方式没有改变，假设文件名为ec_new.bin，用wgrib -V ec_new.bin解码如下：
rec 1:178560:date
W kpds5=135 kpds6=100 kpds7=20 levels=(0,20) grid=255 20 mb 27hr fcst:
&&W=Vertical velocity [Pa s**-1]
&&timerange 0 P1 27 P2 0 TimeU 1&&nx 78915 ny 1 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat&&60.000000 to -10.000000 by 0.250000&&nxny 78915
& && && & long 60.000000 to 150.000000 by 0.250000, (78915 x 1) scan 0 mode 128 bdsgrid 1
&&min/max data -0..0330269&&num bits 15&&BDS_Ref -0.158023&&DecScale 0 BinScale -18
rec 2:326640:date
W kpds5=135 kpds6=100 kpds7=200 levels=(0,200) grid=255 200 mb 27hr fcst:
&&W=Vertical velocity [Pa s**-1]
&&timerange 0 P1 27 P2 0 TimeU 1&&nx 78654 ny 1 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat&&60.000000 to -10.000000 by 0.250000&&nxny 78654
& && && & long 60.000000 to 150.000000 by 0.250000, (78654 x 1) scan 0 mode 128 bdsgrid 1
&&min/max data -3.67&&num bits 16&&BDS_Ref -3.42821&&DecScale 0 BinScale -13
rec 12:2703960:date
TCW kpds5=136 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc 27hr fcst:
&&TCW=Total column water [kg m**-2]
&&timerange 0 P1 27 P2 0 TimeU 1&&nx 1414042 ny 1 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat&&90.000000 to -90.000000 by 0.125000&&nxny 1414042
& && && & long 0.000000 to -0.125000 by 0.125000, (1414042 x 1) scan 0 mode 128 bdsgrid 1
&&min/max data 0..199&&num bits 12&&BDS_Ref 0.230094&&DecScale 0 BinScale -5
解码显示，还是lat&&60.000000 to -10.000000 by 0.250000 long 60.000000 to 150.000000 by 0.250000,但格点数有时 78915，有时78654，反正所有的要素都不一样，而且地面的解码lat&&90.000000 to -90.000000 by 0.125000&&long 0.000000 to -0.125000 by 0.125000，
请教各位高手，请问这到底是什么意思呢？
密码修改失败请联系qq：
在线时间 小时
就是说，即使取出78915个数据，但无法对应到 361 x 281个格点上啊？因为它说是：
lat 60.000000 to -10.000000 by 0.250000&&nxny 101441
long 60.000000 to 150.000000 by 0.250000,
密码修改失败请联系qq：
在线时间 小时
我用以下命令：wgrib W_NAFP_C_ECMF_35_P_C1D00001.bin -V -d all -bin -nh -o&&001.dat & 001.txt解码没有问题：
rec 11::date
CAPE kpds5=59 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc 6hr fcst:
&&CAPE=Convective available potential energy [J kg**-1]
&&timerange 0 P1 6 P2 0 TimeU 1&&nx 2880 ny 1441 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat&&90.000000 to -90.000000 by 0.125000&&nxny 4150080
& && && & long 0.000000 to -0.125000 by 0.125000, (2880 x 1441) scan 0 mode 128 bdsgrid 1
&&min/max data 0 11464.8&&num bits 16&&BDS_Ref 0&&DecScale 0 BinScale -2
也许是你的命令不对，也许你的wgrib不对。地面 long 0.000000 to -0.125000 by 0.125000表示经度范围是0度到359.875（360-0.125）度。
密码修改失败请联系qq：
在线时间 小时
我刚才又试以一下，还是原来的情况，就是
rec 2:326640:date
W kpds5=135 kpds6=100 kpds7=200 levels=(0,200) grid=255 200 mb 27hr fcst:
&&W=Vertical velocity [Pa s**-1]
&&timerange 0 P1 27 P2 0 TimeU 1&&nx 78654 ny 1 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat&&60.000000 to -10.000000 by 0.250000&&nxny 78654
& && && & long 60.000000 to 150.000000 by 0.250000, (78654 x 1) scan 0 mode 128 bdsgrid 1
&&min/max data -3.67&&num bits 16&&BDS_Ref -3.42821&&DecScale 0 BinScale -13
78654个格点和 lat 60.000000 to -10.000000 by 0.250000，long 60.000000 to 150.000000 by 0.250000对应不上，南北271个格点，东西361个数据，一共应该是271x361=101441个才对啊？
那请问您的wgrib版本是多少，我的是在http://www.cpc.ncep.noaa.gov/products/wesley/wgrib.html下载的，不过是编译好的版本，不是源代码的，v1.0.8.9m1 (8-22-05),另外，您的ec文件所有的文件都是一样大小吗？我这里得到的每个预报时次都不一样，有的90多M,有的100多M
密码修改失败请联系qq：
在线时间 小时
本帖最后由 tornadoarps 于
15:44 编辑
我装了windows版本的grads，安装文件名是grads-2.0.a9.oga.1-win32_superpack.exe，装完之后wgrib就可以用了。
你说的100M左右的文件大小，说明你拿到的ec资料是复杂压缩格式，需要转换成简单压缩格式才能被wgrib正确处理。简单压缩格式大小在200M左右。
采用复杂压缩格式下发是为了减小文件大小，节省传输时间和带宽。
密码修改失败请联系qq：
在线时间 小时
楼上正确，就是变成了复杂压缩格式。
密码修改失败请联系qq：
在线时间 小时
非常感谢楼上两位的回答，我的文件是93M，我已经上传到百度网盘http://pan.baidu.com/s/1mgA7jDY，能不能麻烦给看看是不是您说的复杂压缩格式，如果是，那该怎么转换成简单压缩格式呢？不胜感谢！
密码修改失败请联系qq：
在线时间 小时
另外，我还用了openGrads里面的wgrib,就是D:\Program&&Files\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686里的wgrib，版本应该和楼上的grads-2.0.a9.oga.1-win32_superpack.exe一致，用wgrib直接查看版本信息见图，可是解码出来的东西就不是那样了，或许我是拿到的源文件有问题？
我用wgrib shnew.bin -V -d all -bin -nh -o 001.dat &001.txt
结果如图，
001.txt里的内容是：
rec 1:178560:date
W kpds5=135 kpds6=100 kpds7=20 levels=(0,20) grid=255 20 mb 27hr fcst:
&&W=Vertical velocity [Pa s**-1]
&&timerange 0 P1 27 P2 0 TimeU 1&&nx 65535 ny 1 GDS grid 0 num_in_ave 0 missing 0
&&center 98 subcenter 0 process 144 Table 128 scan: WE:NS winds(N/S)
&&latlon: lat&&60.000000 to -10.000000 by 0.250000&&nxny 65535
& && && & long 60.000000 to 150.000000 by 0.250000, (65535 x 1) scan 0 mode 128 bdsgrid 1
内容只有这些，不是内其中的一部分
密码修改失败请联系qq：
在线时间 小时
本帖最后由 tornadoarps 于
21:26 编辑
如果文件只有93M那肯定是复杂压缩格式了。如果你是省级气象部门的工作人员，那你可以从负责资料接收处理的同事那里要到国家气象信息中心下发的“欧洲精细化预报数据格式转换软件安装、使用手册”，里面有详细的复杂压缩格式转简单格式的方法。
如果要不到这个文件，我就简单说一下转换步骤：
1、转换步骤需在linux下实现；
2、安装netcdf
3、安装grib api
4、利用grib api进行转换
netcdf安装：
&&& gunzip netcdf-4.3.0.tar.gz&&& tar xf netcdf-4.3.0.tar&&& mkdir /usr/local/netcdf&&& cd netcdf-4.3.0& &&&& ./configure --prefix=/usr/local/netcdf/ --disable-netcdf-4s& && make&&& make check&&& make install
grib api安装：网址：&&& gunzip grib_api-1.12.3.tar.gz&&& tar xf grib_api-1.12.3.tar&&& mkdir /usr/local/grib_api&&& cd grib_api-1.12.3&&& ./configure --prefix=/usr/local/grib_api --with-netcdf=/usr/local/netcdf --disable-jpeg& && make&&& make check&&& make install
用以下命令进行转换：/usr/local/grib_api/bin/grib_set -r -s packingType=grid_simple $grib_name1&&$grib_name2其中$grib_name1是复杂压缩格式文件名，$grib_name2是简单压缩格式文件名
密码修改失败请联系qq：
在线时间 小时
非常感谢楼上两位的回答，我的文件是93M，我已经上传到百度网盘http://pan.baidu.com/s/1mgA7jDY，能不能麻 ...
我也试了一下你提供的数据，目前的MeteoInfo版本无法正常读取数据。刚下载了netCDF Java最新的库4.5，初步测试能够正常打开该文件并生产数据图层。
Image00108.png (141.46 KB, 下载次数: 0)
23:14 上传
密码修改失败请联系qq：利用python下载ECMWF数据
1.注册ECMWF数据下载账号；2.下载 ECMWF web-API client version in python，下载网址为https://software.ecmwf.int/wiki/display/WEBAPI/Web-API+Downloads，然后pip install 下载的包即可3.登录账号获取下载数据的key，网址https://api.ecmwf.int/v1/key/，然后将获取的信息保存成名为.ecmwfapirc文件，windows放置在 C:\Users\&USERNAME&\.ecmwfapirc，文件内容样例如下：{
: "&a href="https://api.ecmwf.int/v1" "="" style="border-radius: 0 background: 0 border: 0 bottom: float: height: left: line-height: 20 margin: 0 outline: 0 overflow: padding: 0 position: right: top: vertical-align: width: box-sizing: content- min-height: color: rgb(0, 51, 102) !"&https://api.ecmwf.int/v1",
: "XXXXXXXXXXXXXXXXXXXXXX",
"email" : "john."}4.修改需要下载数据的pytho程序样例。重点修改下载的时间，变量（变量编号或变量名），垂直层以及覆盖的区域。ECMWF网站上有关说明：https://softwarepythonchen.ecmwf.int/wiki/display/WEBAPI/Accessing+ECMWF+data+servers+in+batch5.所附程序为下载ECMWF interim气压层的数据#!/usr/bin/python
#-*-coding: UTF-8 -*-
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
# This script downloads ERA-Interim forecasts, on pressure levels.
# Adapt the script to your requirements.
# The data volume for all pressure level data is about 5GB per day, and all pressure level data for more than a single day will exceed the WebAPI limit of 600.000 fields. Thus please restrict the download to what you really need.
# It is likely you need to split your request, this is best done by time periods, ie. first download for month 1, then for month 2, etc.
year=['2012','2013','2014','2015','2016','2017']
month=['05','06','07','08','09']
day=['31','30','31','31','30']
for i in year[0:6]:
for j in month[0:5]:
print i+'-'+j+'-01/to/'+i+'-'+j+'-'+day[k]
server.retrieve({
# Specify the ERA-Interim data archive. Don't change.
"class": "ei",
"dataset": "interim",
"expver": "1",
"stream": "oper",
# pressure levels (levtype:pl), all available levels (levelist)
"levtype": "pl",
"levelist": "100/200/225/250/300/350/400/450/500/550/600/650/700/750/775/800/825/850/875/900/925/950/975/1000",
# forecast (type:fc), from both daily forecast runs (time) with all available forecast steps (step, in hours)
"type": "an",
"time": "00/06/12/18",
"step": "0",
# all available parameters, for codes see http://apps.ecmwf.int/codes/grib/param-db
"param": "60.128/129.128/130.128/131.128/132.128/133.128/135.128/138.128/155.128/157.128",
# two days worth of data
"date": i+'-'+j+'-01/to/'+i+'-'+j+'-'+day[k],
# in 0.75 degrees lat/lon
"grid": "0.75/0.75",
# optionally restrict area to Europe (in N/W/S/E).
#"area": "75/-20/10/60",
# Optionally get output in NetCDF format. However, for NetCDF timestamps (time+step) must not overlap, so use e.g. "time":"00:00:00/12:00:00","step":"12"
"format" : "netcdf",
# set an output file name
"target": "E:\\ERA-Interim\\"+i+"-"+j+"-01to"+i+"-"+j+"-"+day[k]+"-pl.nc"
print "Download finished"
用Python (Anconda)下载欧洲中心 ECMWF 数据
利用Python调用ECMWF欧洲中心API批量下载数据
如何使用python访问ECMWF公共数据集
python调用ECMWF欧洲中心api下载数据
python下载数据集存放位置
使用Python批量下载数据
没有更多推荐了，【图文】ECMWF模式初值场的确定_百度文库
您的浏览器Javascript被禁用，需开启后体验完整功能，
享专业文档下载特权
&赠共享文档下载特权
&10W篇文档免费专享
&每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
ECMWF模式初值场的确定
&&欧洲数值预报初始场的确定
阅读已结束，下载本文到电脑
想免费下载本文？
登录百度文库，专享文档复制特权，积分每天免费拿！
你可能喜欢