使用python读取SMAP grib数据

Question

在阅读SMAP数据(grib格式)使用xarray以cfgrib作为引擎，变量显示为未知的数据的经度值约为539，这是不可能的。我假设变量被标识为“未知”可能是由于缺少SMAP数据的grib表。使用xarray读取数据的总结如下:

< xarray。数据set> Dimensions: (latitude: 1920, longitude: 2560) Coordinates: time datetime64[ns] ... step timedelta64[ns] ... level float64 ... * latitude (latitude) float64 -89.95 -89.86 -89.77 ... 89.77 89.86 89.95 * longitude (longitude) float64 -179.9 -179.8 -179.6 ... 179.6 179.8 539.9 valid_time datetime64[ns] ... Data variables: unknown (latitude, longitude) float32 ... Attributes: GRIB_edition: 2 GRIB_centre: consensus GRIB_centreDescription: Consensus GRIB_subCentre: 65535 Conventions: CF-1.7 institution: Consensus history: 2023-01-17T14:45 GRIB to CDM+CF via cfgrib-0.9.1...

wgrib2的输出如下:

表4.3 = 255(未定义)对于20150402_20150404_anom1。Grb2 ** 1:0:d=1970010100:var discipline=1 master_table=2 parcat =255 parm=255:reserved:anl:

Answer 1

似乎GRIB格式被用于分发数据(可能是为了存储，还是为了简化处理?)，而卫星SMAP不符合GRIB标准(请参阅GRIB标准支持的可用参数)在这里)．这就是为什么输入代码是255(即。“失踪”)。

我下载了文件20150402 _20150404_anom1.grb2并使用xarray 0.18.2和cfgrib 0.9.10.1在Python中打开。我看不到异常的经度。

< xarray。数据set> Dimensions: (latitude: 1920, longitude: 2560) Coordinates: time datetime64[ns] 1970-01-01 step timedelta64[ns] 00:00:00 level float64 0.0 * latitude (latitude) float64 -89.95 -89.86 -89.77 ... 89.77 89.86 89.95 * longitude (longitude) float64 -179.9 -179.8 -179.6 ... 179.6 179.8 179.9 valid_time datetime64[ns] 1970-01-01 Data variables: unknown (latitude, longitude) float32 nan nan nan nan ... nan nan nan Attributes: GRIB_edition: 2 GRIB_centre: consensus GRIB_centreDescription: Consensus GRIB_subCentre: 65535 Conventions: CF-1.7 institution: Consensus history: 2023-01-17T15:53 GRIB to CDM+CF via cfgrib-0.9.1...

您可以在任何情况下使用Python读取和处理数据“未知”作为变量名:

smap = Ds .open_dataset(…)变量['unknown'] lon = ds。变量['经度']lat = ds.variables['经度']