- 积分
- 60
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-3-19
- 最后登录
- 1970-1-1
|
发表于 2022-10-24 15:10:03
|
显示全部楼层
本帖最后由 laoping 于 2022-10-24 15:11 编辑
老师,我也是用你类似的代码“file2 = '20190516.000133.01.19.552' #旧pupf2 = PUP(file2)
print(f2)
data2 = f2.get_data()
file3 = '20190516.000318.01.19.558' #旧pup
f3 = PUP(file3)
print(f3)
data3 = f3.get_data()
allData.append(data0)
allData.append(data2)
#allData.append(data3)
gm=cinrad.calc.GridMapper(allData)
grid = gm(0.001)
print(type((grid)))
print(grid)
#data5 = data4(0.1)
fig = cinrad.visualize.PPI(grid,add_city_names=False) # extent参数可以控制绘制的经纬度范围
fig(r'./')”
为什么会报错“/home/laoping/anaconda3/envs/cinrad/lib/python3.8/site-packages/cinrad/calc.py:425: RuntimeWarning: divide by zero encountered in true_divide return np.ma.average(inp, weights=1 / wgt, axis=2)
/home/laoping/anaconda3/envs/cinrad/lib/python3.8/site-packages/numpy/ma/extras.py:620: RuntimeWarning: invalid value encountered in multiply
wgt = wgt*(~a.mask)
Traceback (most recent call last):
File "readpupold.py", line 32, in <module>
grid = gm(0.01)
File "/home/laoping/anaconda3/envs/cinrad/lib/python3.8/site-packages/cinrad/calc.py", line 455, in __call__
r_attr["nyquist_vel"],
KeyError: 'nyquist_vel'”呢?
|
|