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

eeeee

蓝色随风 发表于 2021-1-22 21:52
您好，我想请问下绘制剖面时为什么底层缺失

因为雷达有仰角

蓝色随风

eeeee 发表于 2021-1-23 14:04
因为雷达有仰角

麻烦再请教一下，剖面图如何增加colorbar

Jaxon

老师你好，我想请问下，我画出来的PPI图片为什么没有右侧的色标以及站点信息，显示如图的警告？

eeeee

Jaxon 发表于 2021-3-8 15:28
老师你好，我想请问下，我画出来的PPI图片为什么没有右侧的色标以及站点信息，显示如图的警告？

我这边无法复现你的错误，那个警告不用管

zytcerulean

请问，PUP STI产品数据一般从哪里获取呢？

zytcerulean

1. from cinrad.io import PUP
   
2. f = PUP('E:/radar/Z_RADR_I_Z9411_20190904011619_P_DOR_SA_STI_NUL_345_NUL.bin')
   
3. data = f.get_data()

复制代码

请问 读取pup产品报错 是什么原因呢

Traceback (most recent call last):
  File "C:\Users\jhg\Desktop\python\readpup.py", line 21, in <module>
    f = PUP('E:/radar/Z_RADR_I_Z9411_20190904011619_P_DOR_SA_STI_NUL_345_NUL.bin')
  File "E:\software\Anaconda3\lib\site-packages\cinrad\io\level3.py", line 31, in __init__
    self.dtype = self._det_product_type(f.prod_desc.prod_code)
  File "E:\software\Anaconda3\lib\site-packages\cinrad\io\level3.py", line 94, in _det_product_type
    raise RadarDecodeError('Unsupported product type {}'.format(spec))
cinrad.error.RadarDecodeError: Unsupported product type 58

单纯xiao可爱

老师不好意思打扰了，我想问下visualize.PPI能直接画出雷达拼图吗，我是用了cinrad.calc.GridMapper后，再用的PPI

1. import cinrad
   
2. from cinrad.io import CinradReader, StandardData
   
3. from cinrad.calc import GridMapper
   
4. from cinrad.visualize import PPI
   
5. import matplotlib.ticker as mticker
   
6. from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
   
7. import matplotlib.pyplot as plt
   
8. 
   
9. plt.rcParams['font.sans-serif'] = ['SimHei'] # 用来正常显示中文标签
   
10. plt.rcParams['axes.unicode_minus'] = False # 用来正常显示负号
    
11. plt.rcParams['font.family']=['SimHei']
    
12. 
    
13. ele = 0 #第1个仰角
    
14. radius = 200 #绘制图像的范围大小
    
15. 
    
16. f = CinradReader('E:\\421\\Z_RADR_I_Z9421_20201117120209_O_DOR_CD_CAP.bin')
    
17. h = CinradReader('E:\\240\\Z_RADR_I_Z9240_20201117120402_O_DOR_SC_CAP.bin')
    
18. #j = CinradReader('E:\\411\\Z_RADR_I_Z9411_20201117115900_O_DOR_SA_CAP.bin')
    
19. 
    
20. r1 = f.get_data(ele,radius,'REF')
    
21. r2 = h.get_data(ele,radius,'REF')
    
22. #r3 = j.get_data(ele,radius,'REF')
    
23. #gm = GridMapper([r1, r2])
    
24. #grid = gm(0.1)
    
25. gm = cinrad.calc.GridMapper([r1, r2])
    
26. fig = cinrad.visualize.PPI(gm, dpi=200, add_city_names=True)
    

复制代码

1. AttributeError                            Traceback (most recent call last)
   
2. <ipython-input-63-bb4965aa25cd> in <module>
   
3.      24 #grid = gm(0.1)
   
4.      25 gm = cinrad.calc.GridMapper([r1, r2])
   
5. ---> 26 fig = cinrad.visualize.PPI(gm, dpi=200, add_city_names=True)
   
6. 
   
7. E:\anaconda\lib\site-packages\cinrad-1.6.2-py3.8-win-amd64.egg\cinrad\visualize\ppi.py in __init__(self, data, fig, norm, cmap, nlabel, label, dpi, highlight, coastline, extent, section, style, add_city_names, plot_labels, **kwargs)
   
8.      80     ):
   
9.      81         self.data = data
   
10. ---> 82         self.dtype = get_dtype(data)
    
11.      83         self.settings = {
    
12.      84             "cmap": cmap,
    
13. 
    
14. E:\anaconda\lib\site-packages\cinrad-1.6.2-py3.8-win-amd64.egg\cinrad\common.py in get_dtype(data)
    
15.       6
    
16.       7 def get_dtype(data: Dataset) -> str:
    
17. ----> 8     all_data = list(data.keys())
    
18.       9     geo_var_name = ["longitude", "latitude", "height", "x_cor", "y_cor", "RF"]
    
19.      10     for i in geo_var_name:
    
20. 
    
21. AttributeError: 'GridMapper' object has no attribute 'keys'

复制代码

atmosphereic

老师您好，请问运行程序画图画到一半的时候（能够画一半的图出来）报错
    ValueError: buffer size must be a multiple of element size
这个应该怎么处理呀，报错的程序是下面这行代码
    f = CinradReader(sa_radar_file[nFiles]) #老版本数据
谢谢老师的解答

atmosphereic

grug 发表于 2019-1-8 12:54
楼主给力，我用python3.6版本pip安装cinrad==1.2， from cinrad.io import CinradReader总会显示ModuleNotF ...

您好，请问您添加经纬度成功了吗，想向您请教一下

eeeee

单纯xiao可爱 发表于 2021-3-25 16:50
老师不好意思打扰了，我想问下visualize.PPI能直接画出雷达拼图吗，我是用了cinrad.calc.GridMapper后，再 ...

grid = gm(0.1)
然后把grid传入PPI里