Barnes滤波

镜灵17

雨落森林 发表于 2024-5-26 10:22
可能是输入没输入对，你可以详细介绍一下你是怎么做的吗？比如pybarnes是什么版本？以及，你是如何输入的 ...

pybarnes版本是0.2.0；步骤是这样：先读取wrfout文件的位势高度，然后插值到500hPa上（没有缺测，这一步print出来都正常，也能画出来），再用Barnes滤波得到大尺度的位势高度，这步做完之后print出来的数据很多都是nan，也画不出来了

dst = ncdataset('文件路径')

geopt = getvar(dst, "geopt", 60)
p       = getvar(dst, "pressure", 60)
hgt    = interplevel(geopt, p, 500)

f = BarnesFilter(hgt)
hgt_g = f.lowpass(g=0.3, c=30000) #低通滤波,保留大尺度
hgt_g = hgt_g*units.m**2*units.s**-2

雨落森林

镜灵17 发表于 2024-5-27 10:52
pybarnes版本是0.2.0；步骤是这样：先读取wrfout文件的位势高度，然后插值到500hPa上（没有缺测，这一步p ...

理论上这个应该报错的才对，0.2版本我已经不想支持非常规的xr数组了，这个没有报错，应该是把不正常的东西识别成了经纬度，建议手动给它赋予经纬度，像这样：

f = BarnesFilter(hgt.data, lon=hgt.XLONG.data, lat=hgt.XLAT.data)
hgt_g = f.lowpass(g=0.3, c=30000) #低通滤波,保留大尺度

镜灵17

雨落森林 发表于 2024-5-28 23:33
理论上这个应该报错的才对，0.2版本我已经不想支持非常规的xr数组了，这个没有报错，应该是把不正常的东 ...

非常感谢！现在输出的数据看起来正常了！

hgsndg

请问可以直接输入三维数据对每一层进行滤波吗？

雨落森林

hgsndg 发表于 2024-7-8 15:02
请问可以直接输入三维数据对每一层进行滤波吗？

可以的，你可以试试

侯祎林

楼主，请问代码运行时可以关闭
Caculating the first revision...
Caculating the second revision...
的显示么

雨落森林

侯祎林 发表于 2024-8-2 09:33
楼主，请问代码运行时可以关闭
Caculating the first revision...
Caculating the second revision...

可以，lowpass或者bandpass的时候加个print_progress=False就可以取消了，这个输出主要是因为大的数组滤波的时候很慢，输出进度来看的，小数据的话很快就好了，反倒是累赘。

侯祎林

雨落森林 发表于 2024-8-2 15:00
可以，lowpass或者bandpass的时候加个print_progress=False就可以取消了，这个输出主要是因为大的数组滤 ...

我的代码是这样的
f = BarnesFilter(data_, lon=lon_1, lat=lat_1)
g = 0.3
c = 150000
data_b = f.lowpass(g=g, c=c, print_progress=False)
但是程序仍然生成
Caculating the first revision...
Caculating the second revision...
请问怎么解决呢

雨落森林

侯祎林 发表于 2024-8-2 15:13
我的代码是这样的
f = BarnesFilter(data_, lon=lon_1, lat=lat_1)
g = 0.3

已经修复了这个bug并重新发布了，可以重新安装一下试试看，pip uninstall pybarnes卸载，然后pip install pybarnes==0.2.0 -i https://pypi.python.org/simple用pypi源重新安装一下

侯祎林

雨落森林 发表于 2024-8-3 12:20
已经修复了这个bug并重新发布了，可以重新安装一下试试看，pip uninstall pybarnes卸载，然后pip install ...

现在可以了，谢谢大佬