MeteoinfoLab网格到网格插值

lovechang1314

王老师好，最近使用EC预报数据做网格插值，但总是出错，得到的结果全部都是NAN。。

无论我尝试一个小的范围还是就直接插值到数据本身的网格都是输出nan，我原以为是版本问题，但试过[linint2](http://meteothink.org/docs/meteo ... ctions/linint2.html)之后，发现运行正常，实在找不到问题，请王老师帮看看...

下面是我的代码
```
fn = addfile(r'D:\ECMWF\20211226\00\W_NAFP_C_ECMF_20211226064448_P_C1D12260000010300001.bz2')
gh = fn['Geopotential_Height_isobaric'][:]
ghdata =gh[0]
gx = ghdata.dimvalue(1)
gy = ghdata.dimvalue(0)
ngy = interpolate.linint2(ghdata, gx, gy)
print(ngy)

```
图片为上面代码输出结果，真的尝试各种验证，缺失找不到问题了。。。
关于数据，我想王老师那边应该也能拿到EC的C1D数据，我就不上传了

lovechang1314

又看了一下linint2的网页，感觉和示例代码不太一样，文档中params的形参顺序是x/y/z，但是示例中是z/x/y，不知道是这个原因不，晚点再试试

MeteoInfo

你代码获取的 gx, gy 并非数据X, Y 坐标的值。

lovechang1314

不是应该这么获取吗？

lovechang1314

王老师是这样的，我按照linint2里的lon和lat生成了X/Y，依然是nan

lovechang1314

我再试试王老师，看了一些别的示例，XY似乎应该从文件获取，而不是变量。谢谢您，有结果我及时在贴下回复，如果还有问题也请您指导

lovechang1314

看了些别的代码示例，似乎应该直接从文件对象获取X/Y，而不是变量对象，但这个var.dimvalue获取的是什么呢？回头我再尝试一下，有问题还请您指导

lovechang1314

王老师，还是不行额。。。

我感觉dimvalue获取的方式是可以的，因为我在其他地方也看到dimvalue获取X和Y。
即便是我自己设定X和Y也不可以，例如下代码：
```
lon = arange(70,140,1.5)
lat = arange(50,5,-1.5)
ngy = interpolate.linint2(ghdata, gx, gy)
```
这样获取的插值网格仍然是nan，求指导。。。

MeteoInfo

lovechang1314 发表于 2022-3-27 14:06
王老师，还是不行额。。。

我感觉dimvalue获取的方式是可以的，因为我在其他地方也看到dimvalue获取X和Y ...

一个变量可能有多维，比如一个四维变量 v（时间、高度、纬度、经度），那么 v.dimvalue(0) 获取第一个维的数据，也就是时间，v.dimvalue(2) 和 v.dimvalue(3) 分别是纬度和经度值。

还有一个问题需要注意，就是数据的投影。

lovechang1314

MeteoInfo 发表于 2022-3-28 09:00
一个变量可能有多维，比如一个四维变量 v（时间、高度、纬度、经度），那么 v.dimvalue(0) 获取第一个维 ...

王老师，如您所说和顺序有关，我检查了一下，获取的gx和gy没有问题，是这个顺序。

投影的话EC的预报数据都是等经纬度投影，应该不涉及这个转换的问题，很奇怪。

如果是带投影的网格到网格应该如何做呢？