使用pygrib写入插入的grib2数据会导致grib文件不可用

Question

我试图使用pygrib从grib2文件读取数据，使用python插入它，并将其写入另一个文件。我尝试过pygrib和eccodes，两者都产生同样的问题。输出文件的大小增加了3倍，但是当我试图在Weather and Climate Toolkit等应用程序中查看数据时，它列出了所有变量，但是"没有数据"。当绘制。如果我使用相同的脚本，不插入数据，而只是将其写入新文件，它在WCT中工作得很好。如果我使用wgrib2，它列出了所有的grib消息，但如果我使用wgrib2 -V，它工作在未改变的数据上，但产生错误";***致命错误:不支持:代码表5.6=0 ***&";对于插值数据。我在我的python脚本做错了什么?下面是我编写该文件的示例(使用pygrib 2.05和2.1.3的结果相同)。我使用了一个基本的hrrr文件作为示例。

import pygrib
import numpy as np
import sys
def writeNoChange():
# This produces a useable grib file.
filename = 'hrrr.t00z.wrfprsf06.grib2'
outfile = 'test.grib2'
grbs = pygrib.open(filename)
with open(outfile, 'wb') as outgrb:
for grb in grbs:
msg = grb.tostring()
outgrb.write(msg)
outgrb.close()
grbs.close()
def writeChange():
# This method produces a grib file that isn't recognized by WCT

filename = 'hrrr.t00z.wrfprsf06.grib2'
outfile = 'testChange.grib2'
grbs = pygrib.open(filename)
with open(outfile, 'wb') as outgrb:
for grb in grbs:
vals = grb.values * 1
grb['values'] = vals
msg = grb.tostring()
outgrb.write(msg)
outgrb.close()
grbs.close()
#-------------------------------
if __name__ == "__main__":
writeNoChange()
writeChange()

Answer 1

表5.6 for GRIB2 (https://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc/)与"ORDER OF SPATIAL differentiing"有关。

由于某些原因，当你修改grb['values']时，它设置grb[' orderofspatialdifference '] = 0，这"wgrib2 -V"不喜欢。因此，在改变'values'之后，将' orderofspatialdifference '更改为初始值:

orderOfSpatialDifferencing = grb['orderOfSpatialDifferencing']
grb['values']= [new values]
grb['orderOfSpatialDifferencing'] = orderOfSpatialDifferencing

这在让wgrib2 -V工作方面对我有效，但弄乱了数据。可能第5节中的其他一些变量也需要修改。