![](/img/trans.png)
[英]Using XArray.isel to access data in GRIB2 file from a specific location?
[英]NDFD GRIB2 how to fix mirrored data when using xarray
我在新英格兰提取风速 grib 文件的代码代码:
import pandas as pd
import numpy as np
import requests
import cfgrib
import xarray as xr
resp = requests.get('https://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/DF.gr2/DC.ndfd/AR.neast/VP.001-003/ds.wspd.bin', stream=True)
f = open('..\\001_003wspd.grib2', 'wb')
f.write(resp1.content)
f.close()
xr_set = xr.load_dataset('..\\001_003wspd.grib2', engine='cfgrib')
xr_set.si10[0].plot(cmap=matplotlib.pyplot.cm.coolwarm)
如您所见,它从东到西每隔一条线镜像一次。 缅因州是最明显的。
我相信这不是代码问题,而是未正确编写的文件。 如果你每两行只取一行,你会得到一张正确的地图:
import numpy as np
import requests
import xarray as xr
from fs.tempfs import TempFS
resp = requests.get('https://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/DF.gr2/DC.ndfd/AR.neast/VP.001-003/ds.wspd.bin', stream=True)
with TempFS() as tempfs:
path = tempfs.getsyspath("001_003wspd.grib2")
f = open(path, 'wb')
f.write(resp.content)
f.close()
ds = xr.load_dataset(path, engine='cfgrib')
ds = ds.isel(y=np.arange(len(ds.y))[1::2])
ds.si10.isel(step=15).plot(cmap="coolwarm", x='longitude', y='latitude')
问题与 grib 文件的排序方式有关。 如果执行以下命令,您可以看到顺序。
$wgrib2 -grid in.grib2
您可能会看到类似“(2345 x 1597) 输入 WE|EW:SN output WE:SN”的内容
WE|EW:SN 需要更改为 WE:SN。 要更改顺序,请执行以下命令
wgrib2 in.grib2 -ijsmall_grib 1:2345 1:1597 out.grib2
现在您将看到 out.grib2 具有正确的 WE:SN 顺序并且可以与 xarray 一起使用
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.