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

eeeee

atmosphereic 发表于 2021-3-28 09:19
老师您好，请问运行程序画图画到一半的时候（能够画一半的图出来）报错
    ValueError: buffer size must ...

可能是模块没有正确识别文件类型，报错的是什么数据啊

atmosphereic

eeeee 发表于 2021-3-28 22:55
可能是模块没有正确识别文件类型，报错的是什么数据啊

老师，我的数据是下面这个类型的
Z_RADR_I_Z9010_20060706000200_O_DOR_SA_CAP.bin
我进行的操作是指定一个仰角，循环画一天的，最终应该有240幅图，在出了115幅图后就报了这个错

老师，我还想请教您一个问题，我在保存图片时只指定了保存目录，然后图片的名称是自定义的，请问老师这个自定义名称是按什么格式呀
Z9010_20060706000240_0.5_150_REF.png     原文件 Z_RADR_I_Z9010_20060706000200_O_DOR_SA_CAP.bin
Z9010_20060706000843_0.6_150_REF.png     原文件 Z_RADR_I_Z9010_20060706000800_O_DOR_SA_CAP.bin
这是其中一些图片的名称，但是这个时间我在原文件里并未找到（秒级数字不同）

单纯xiao可爱

eeeee 发表于 2021-3-28 22:53
grid = gm(0.1)
然后把grid传入PPI里

多谢老师！！我再试试

eeeee

atmosphereic 发表于 2021-3-29 11:19
老师，我的数据是下面这个类型的
Z_RADR_I_Z9010_20060706000200_O_DOR_SA_CAP.bin
我进行的操作是指定 ...

文件里面是存有扫描时间精确到秒的，只是文件名没有反映出来。
报错的数据可以发一份到我的邮箱吗，274555447@qq.com

atmosphereic

eeeee 发表于 2021-4-4 22:56
文件里面是存有扫描时间精确到秒的，只是文件名没有反映出来。
报错的数据可以发一份到我的邮箱吗，

已发送数据，请老师查收，谢谢老师呀

eeeee

atmosphereic 发表于 2021-4-5 11:40
已发送数据，请老师查收，谢谢老师呀

我感觉是里面有几个数据损坏了，没法读出来，之后可以考虑支持一下读取这种不完整的文件

atmosphereic

eeeee 发表于 2021-4-5 15:48
我感觉是里面有几个数据损坏了，没法读出来，之后可以考虑支持一下读取这种不完整的文件

好呢好呢，谢谢老师！非常期待

蔡小亦

大神你好，我在去年12月23日下载的CyanideCN-PyCINRAD-master，成功安装并可以画出REF、VEL、VIL、ET等图。有两个问题请教您：
1、现在想画CC、ZDR、KDP等数据，在您上一个帖子里面看到您有更新一个PyCINRAD的压缩包，我下载后不知道该如何安装使用？还是我之前下载的CyanideCN-PyCINRAD-master里面已经包含了这个压缩包了呢？
2、我在github上和两个帖子（包括评论）找了很久都没有看到如何画CC、ZDR、KDP代码块，请问是否有示例呢？
以上，非常感谢！

更新一下程序。
新增功能：
1.CC1.0雷达的多普勒速度和基本反射率绘制
2.CA雷达多普勒速度和基本反射率绘制
3.双偏振雷达的基本反射率，差分反射率，协相关系数，差分传播相移率的绘制
4.优化雷达类型的识别
5.一些其他Bug的修复
6.组合反射率绘制

单纯xiao可爱

老师我又来打扰了。。我想做组合反射率的雷达拼图，拼图之前用的基本发射率成功了。但是计算了组合反射率后没有了'distance'，想问一下这时候咋整呀

1. cr3 = cinrad.calc.quick_cr([f3.get_data(i, 200, 'REF') for i in f3.angleindex_r])
   
2. cr4 = cinrad.calc.quick_cr([f4.get_data(i, 200, 'REF') for i in f4.angleindex_r])
   
3. print(cr3)
   
4. print(cr4)
   
5. gm = GridMapper([cr3,cr4])
   
6. grid = gm(0.01)
   
7. fig = cinrad.visualize.PPI(grid, dpi=200, style='white',add_city_names=True)[code]KeyError            
   
8. [code]<xarray.Dataset>
   
9. Dimensions:    (latitude: 1000, longitude: 1000)
   
10. Coordinates:
    
11.   * longitude  (longitude) float64 119.3 119.3 119.3 119.3 ... 124.0 124.0 124.0
    
12.   * latitude   (latitude) float64 37.11 37.11 37.12 37.12 ... 40.7 40.71 40.71
    
13. Data variables:
    
14.     CR         (longitude, latitude) float64 nan nan nan nan ... nan nan nan nan
    
15. Attributes:
    
16.     elevation:        0
    
17.     range:            200
    
18.     scan_time:        2020-11-17 13:05:03.526000
    
19.     site_code:        Z9411
    
20.     site_name:        大连
    
21.     site_longitude:   121.64888888888889
    
22.     site_latitude:    38.90972222222222
    
23.     tangential_reso:  1.0
    
24.     nyquist_vel:      30.73
    
25.     task:             VCP21
    
26. <xarray.Dataset>
    
27. Dimensions:    (latitude: 1000, longitude: 1000)
    
28. Coordinates:
    
29.   * longitude  (longitude) float64 119.9 119.9 119.9 119.9 ... 124.6 124.6 124.6
    
30.   * latitude   (latitude) float64 38.87 38.87 38.88 38.88 ... 42.47 42.47 42.47
    
31. Data variables:
    
32.     CR         (longitude, latitude) float64 nan nan nan nan ... nan nan nan nan
    
33. Attributes:
    
34.     elevation:        0
    
35.     range:            200
    
36.     scan_time:        2020-11-17 13:05:03.346000
    
37.     site_code:        Z9417
    
38.     site_name:        营口
    
39.     site_longitude:   122.27194444444444
    
40.     site_latitude:    40.672777777777775
    
41.     tangential_reso:  1.0
    
42.     nyquist_vel:      30.72
    
43.     task:             VCP21

复制代码

                     Traceback (most recent call last)
F:\anacondo\lib\site-packages\xarray\core\dataset.py in _construct_dataarray(self, name)
   1258         try:
-> 1259             variable = self._variables[name]
   1260         except KeyError:

KeyError: 'distance'

During handling of the above exception, another exception occurred:

KeyError                                  Traceback (most recent call last)
<ipython-input-20-b91e71348677> in <module>
      7 print(cr3)
      8 print(cr4)
----> 9 gm = GridMapper([cr3,cr4])
     10 grid = gm(0.01)
     11 fig = cinrad.visualize.PPI(grid, dpi=200, style='white',add_city_names=True)

F:\anacondo\lib\site-packages\cinrad\calc.py in __init__(self, fields, max_dist)
    380         self.data_ravel = np.ma.hstack([i[self.dtype].values.ravel() for i in fields])
    381         self.dist_ravel = np.hstack(
--> 382             [
    383                 np.broadcast_to(i["distance"], i["longitude"].shape).ravel()
    384                 for i in fields

F:\anacondo\lib\site-packages\cinrad\calc.py in <listcomp>(.0)
    381         self.dist_ravel = np.hstack(
    382             [
--> 383                 np.broadcast_to(i["distance"], i["longitude"].shape).ravel()
    384                 for i in fields
    385             ]

F:\anacondo\lib\site-packages\xarray\core\dataset.py in __getitem__(self, key)
   1357
   1358         if hashable(key):
-> 1359             return self._construct_dataarray(key)
   1360         else:
   1361             return self._copy_listed(np.asarray(key))

F:\anacondo\lib\site-packages\xarray\core\dataset.py in _construct_dataarray(self, name)
   1259             variable = self._variables[name]
   1260         except KeyError:
-> 1261             _, name, variable = _get_virtual_variable(
   1262                 self._variables, name, self._level_coords, self.dims
   1263             )

F:\anacondo\lib\site-packages\xarray\core\dataset.py in _get_virtual_variable(variables, key, level_vars, dim_sizes)
    169         ref_var = dim_var.to_index_variable().get_level_variable(ref_name)
    170     else:
--> 171         ref_var = variables[ref_name]
    172
    173     if var_name is None:

KeyError: 'distance'[/code][/code]

eeeee

单纯xiao可爱 发表于 2021-4-9 17:14
老师我又来打扰了。。我想做组合反射率的雷达拼图，拼图之前用的基本发射率成功了。但是计算了组合反射率后 ...

组合反射率暂时不支持拼图