编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

读SMAP格丽数据使用python

而阅读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:保留:退火:

回答

似乎格丽格式是用于分发的数据(可能存储,便于处理?),而卫星SMAP不适合在格丽标准(见可用参数支持格丽标准[这][1])。这就是为什么是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... ``` You can read and process the data in Python in any case using `'unknown'` as the variable name: ``` ds = xr.open_dataset(...) smap = ds.variables['unknown'] lon = ds.variables['longitude'] lat = ds.variables['latitude'] ``` [![2]][2] [1]: https://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc/grib2_table4-1.shtml [2]: https://i.stack.imgur.com/OYe7f.jpg

编辑总结

取消

\ begingroup美元我xarray版本是2022.12.0出现异常经度值。这个问题消失了降级的时候出现0.18.2。 \ endgroup美元
- - - - - -irok94
1月18日在9:04
\ begingroup美元请注意,时间戳在xarray和wgrib输出不正确。解决时间问题呢? \ endgroup美元
- - - - - -irok94
1月19日在41
\ begingroup美元一种文件格式的数据文件包含三天复合材料与文件名中提供的开始和结束日期。这可能是为什么时间戳在格丽是“错误的”,只有一个二维数组。我会用cdo如果你真的想要改变该文件的时间戳:cdo settaxis, 2015-04-02, 20150402 _20150404_anom1 00:00, 1天。grb2 out.grb2 \ endgroup美元
- - - - - -英格瓦•卢卡斯
在1月19日15:43

添加一个评论 |

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复

创建代码用引号”或式子~篱笆
' ' '
像这样
' ' '
添加语言标识符来突出显示的代码
”“python
def函数(foo):
打印(foo)
' ' '
把段落之间的回报
在结束linebreak添加2个空间
_italic_或* *的* *
报价将>开始的
做链接(尽量使用https)
< https://example.com >
(例子)(https://example.com)
< a href = " https://example.com " > < / >示例
MathJax方程罪\ \ ^ 2θ美元

格式帮助»
回答帮助»

MathJax帮助»