读SMAP格丽数据使用python

Question

而阅读SMAP数据(格丽格式)使用与cfgrib xarray引擎,作为显示的变量未知的的经度数据有一个价值约539,不应该是可能的。我假设变量被确定为“未知的可能是由于缺少SMAP格丽表数据。总结的数据读取使用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纪律的master_table = 1 = 2 parmcat = 255改= 255:保留:退火:

Answer 1

似乎格丽格式是用于分发的数据(可能存储,便于处理?),而卫星SMAP格丽标准(见下不适合可用参数支持一种文件格式标准在这里)。这就是为什么是255(即输入代码。“失踪”)。

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

< 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“未知”变量名:

ds = xr.open_dataset (…) smap = ds。变量(未知的)朗= ds。变量(经度的)lat = ds.variables(“纬度”)