GridMapper做雷达拼图，无法显示
tooleek opened this issue · 10 comments
`for binFile in binFileList:
f = StandardData(binFile)
refData = f.get_data(0, drange, 'REF')
allData.append(refData)
gm = GridMapper(allData)
grid = gm(0.1)
fig = PPI(grid, dpi=75, add_city_names=False)
fig._save(outPath)`
请问，使用该方法做雷达拼图，但是最终显示的图缺少雷达数据，是什么原因呢
代码看起来没问题，确定基数据里是有回波的吗
嗯，基数据是有回波的
具体拼图画出来是啥样的
三个站的单站图如下：
合成后的图如下：
代码如下：
执行过程中，抛了如下异常：
D:\Software\Anaconda3\lib\site-packages\cinrad\calc.py:425: RuntimeWarning: divide by zero encountered in true_divide
return np.ma.average(inp, weights=1 / wgt, axis=2)
D:\Software\Anaconda3\lib\site-packages\numpy\ma\extras.py:621: RuntimeWarning: invalid value encountered in multiply
wgt = wgt*(~a.mask)
D:\Software\Anaconda3\lib\site-packages\cinrad\visualize\ppi.py:233: RuntimeWarning: All-NaN slice encountered
"Max: {:.1f}{}".format(np.nanmax(var), unit[self.dtype]),
我把参数调成 grid = gm(0.001)，就有数据了，看来是0.1太小了这个原因导致的
现在有一个新的问题，就是组合反射率的雷达拼图，要抛异常，找不到 distance
组合反射率目前不支持拼图
我也是类似的问题，“file2 = '20190516.000133.01.19.552' #旧pup
f2 = PUP(file2)
data2 = f2.get_data()
gm=cinrad.calc.GridMapper([data0,data2])
grid = gm(0.1)”前面一步没有报错，最后一步“grid = gm(0.1)”却会报错“ 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'”请问这是为什么？
我也是类似的问题，“file2 = '20190516.000133.01.19.552' #旧pup f2 = PUP(file2) data2 = f2.get_data() gm=cinrad.calc.GridMapper([data0,data2]) grid = gm(0.1)”前面一步没有报错，最后一步“grid = gm(0.1)”却会报错“ 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'”请问这是为什么？
@muyi151 这个算是一个bug，后续会修复
我在对多站雷达基本反射率做拼图时结果如图，画单站图时并不会如此，请问这是什么原因呢？不知道该怎么修改，请求帮助
