NDFD GRIB2 如何在使用 xarray 时修复镜像数据

Question

我在新英格兰提取风速 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)

这给出：

如您所见，它从东到西每隔一条线镜像一次。缅因州最明显。

Answer 1

我相信这不是代码问题，而是文件写入不正确。如果你每两行只取一行，你会得到一个正确的地图：

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')

Answer 2

问题与 grib 文件的排序方式有关。如果执行以下命令，您可以看到顺序。

$wgrib2 -grid in.grib2

您可能会看到类似“(2345 x 1597) 输入 WE|EW:SN 输出 WE:SN”之类的内容

WE|EW:SN 需要改为WE:SN。要更改顺序，请执行以下命令

wgrib2 in.grib2 -ijsmall_grib 1:2345 1:1597 out.grib2

现在您将看到 out.grib2 具有正确的 WE:SN 顺序并且可以与 xarray

一起使用

NDFD GRIB2 如何在使用 xarray 时修复镜像数据

NDFD GRIB2 how to fix mirrored data when using xarray

python

python-xarray

grib

cfgrib

ndfd