PyCINRAD模块处理雷达基数据指南（发布1.7.0版本）(修复经纬度标注问题)

不不106

eeeee 发表于 2019-9-4 16:46
这个不是我的模块的问题，可能是有模块的冲突

好的，谢谢答复

不不106

eeeee 发表于 2019-9-4 16:46
这个不是我的模块的问题，可能是有模块的冲突

你好，我在运行模板时，第一次是正常出图的，但第二遍时就会出现重新下载模板的错误。原程序如下，错误如图：
import cinrad
from cinrad.io import CinradReader, StandardData
from cinrad.visualize import PPI
import matplotlib.ticker as mticker
f = cinrad.io.CinradReader("D:/leida/nb/Z_RADR_I_Z9574_20180304110000_O_DOR_SA_CAP.bin.bz2")
r= f.get_data(0,230,'REF')
fig = cinrad.visualize.PPI(r)
fig('D:\\leida\\nb')

不不106

不不106 发表于 2019-9-5 15:06
你好，我在运行模板时，第一次是正常出图的，但第二遍时就会出现重新下载模板的错误。原程序如下，错误如 ...

问题已解决，谢谢

一二宿

求助！
使用StandarData读取南京的S波段双偏振业务雷达，显示：

1. 1.   File "C:\Users\Marcus\Anaconda3\envs\cartopy\lib\site-packages\cinrad-1.4.2-py3.7.egg\cinrad\io\level2.py", line 437, in __init__
      
   2.     self._parse()
      
   3.   File "C:\Users\Marcus\Anaconda3\envs\cartopy\lib\site-packages\cinrad-1.4.2-py3.7.egg\cinrad\io\level2.py", line 462, in _parse
      
   4.     self.task_name = merge_bytes(task['task_name'][0]).decode()
      
   5. UnicodeDecodeError: 'utf-8' codec can't decode byte 0xd5 in position 7: invalid continuation byte

   复制代码

<div class=\"blockcode\"><div id=\"code_Tx6\"><ol><li>  File "C:\Users\Marcus\Anaconda3\envs\cartopy\lib\site-packages\cinrad-1.4.2-py3.7.egg\cinrad\io\level2.py", line 437, in __init__<br /> <li>    self._parse()<br /> <li>  File "C:\Users\Marcus\Anaconda3\envs\cartopy\lib\site-packages\cinrad-1.4.2-py3.7.egg\cinrad\io\level2.py", line 462, in _parse<br /> <li>    self.task_name = merge_bytes(task['task_name'][0]).decode()<br /> <li>UnicodeDecodeError: 'utf-8' codec can't decode byte 0xd5 in position 7: invalid continuation byte</ol></div><em onclick=\"copycode($('code_Tx6'));\">复制代码</em></div>

尝试用了解压缩后的数据，显示：

1. UnicodeDecodeError: 'utf-8' codec can't decode byte 0xd0 in position 7: invalid continuation byte
   

UnicodeDecodeError: 'utf-8' codec can't decode byte 0xd0 in position 7: invalid continuation byte<br />

读取单偏振的业务雷达，以及广州S波段双偏振业务雷达都没有问题。
这是为什么呢

sskkyy2008

请教老师，能否根据给定的某个经度、纬度，把该点的反射率 或者速度 值提取出来？？？

sskkyy2008

sskkyy2008 发表于 2019-9-24 16:28
请教老师，能否根据给定的某个经度、纬度，把该点的反射率 或者速度 值提取出来？？？

import cinrad
import numpy as np
from cinrad.io import CinradReader, StandardData
from scipy import spatial
f = CinradReader("C:/Users/gsy/Desktop/RADA_CHN_DOR_L2_O-Z9837-CA-CAP-20190820003000.bin/2017082108000.05V") #老版本数据
f.set_code('Z9816')
f._update_radar_info()
r = f.get_data(0, 230, 'REF')
tree = spatial.KDTree(np.dstack((r.lon.ravel(),r.lat.ravel()))[0])
t2=tree.query(np.array([104.74,31.78]))
#print(t2)
t3=np.unravel_index(t2[1],r.data.shape)
t3=r.data[t3]
print(t3)

qxjzkx@163.com

Traceback (most recent call last):
  File "D:\pys\read_rad.py", line 14, in <module>
    fig= cinrad.visualize.PPI(rl,dpi=75)
  File "D:\Python37\lib\site-packages\cinrad\visualize\ppi.py", line 94, in __init__
    self._plot(**kwargs)
  File "D:\Python37\lib\site-packages\cinrad\visualize\ppi.py", line 134, in _plot
    dtype = self.data.dtype
AttributeError: 'list' object has no attribute 'dtype'

eeeee

qxjzkx@163.com 发表于 2019-10-29 08:43
Traceback (most recent call last):
  File "D:\pys\read_rad.py", line 14, in
    fig= cinrad.visua ...

qxjzkx@163.com

底图省界还是有点问题，其他很好，牛B。

eeeee

qxjzkx@163.com 发表于 2019-10-30 10:59
绘图问题：
1、0到5dBZ白灰色区域偏多
2、强度偏弱，这是一次强对流天气过程

麻烦你看清楚这是VIL不是反射率