爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: eeeee

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

  [复制链接]

新浪微博达人勋

 楼主| 发表于 2021-3-28 22:55:24 | 显示全部楼层
atmosphereic 发表于 2021-3-28 09:19
老师您好,请问运行程序画图画到一半的时候(能够画一半的图出来)报错
    ValueError: buffer size must ...

可能是模块没有正确识别文件类型,报错的是什么数据啊
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-29 11:19:58 | 显示全部楼层
本帖最后由 atmosphereic 于 2021-3-29 11:25 编辑
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
这是其中一些图片的名称,但是这个时间我在原文件里并未找到(秒级数字不同)
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-31 13:17:37 | 显示全部楼层
eeeee 发表于 2021-3-28 22:53
grid = gm(0.1)
然后把grid传入PPI里

多谢老师!!我再试试
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-4-4 22:56:59 | 显示全部楼层
atmosphereic 发表于 2021-3-29 11:19
老师,我的数据是下面这个类型的
Z_RADR_I_Z9010_20060706000200_O_DOR_SA_CAP.bin
我进行的操作是指定 ...

文件里面是存有扫描时间精确到秒的,只是文件名没有反映出来。
报错的数据可以发一份到我的邮箱吗,274555447@qq.com
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-4-5 11:40:45 | 显示全部楼层
eeeee 发表于 2021-4-4 22:56
文件里面是存有扫描时间精确到秒的,只是文件名没有反映出来。
报错的数据可以发一份到我的邮箱吗,

已发送数据,请老师查收,谢谢老师呀
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-4-5 15:48:58 | 显示全部楼层
本帖最后由 eeeee 于 2021-4-5 15:52 编辑
atmosphereic 发表于 2021-4-5 11:40
已发送数据,请老师查收,谢谢老师呀

我感觉是里面有几个数据损坏了,没法读出来,之后可以考虑支持一下读取这种不完整的文件
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-4-6 10:25:02 | 显示全部楼层
eeeee 发表于 2021-4-5 15:48
我感觉是里面有几个数据损坏了,没法读出来,之后可以考虑支持一下读取这种不完整的文件

好呢好呢,谢谢老师!非常期待
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-4-7 15:00:03 | 显示全部楼层
大神你好,我在去年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.组合反射率绘制

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-4-9 17:14:24 | 显示全部楼层
老师我又来打扰了。。我想做组合反射率的雷达拼图,拼图之前用的基本发射率成功了。但是计算了组合反射率后没有了'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]
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-4-10 14:49:39 | 显示全部楼层
单纯xiao可爱 发表于 2021-4-9 17:14
老师我又来打扰了。。我想做组合反射率的雷达拼图,拼图之前用的基本发射率成功了。但是计算了组合反射率后 ...

组合反射率暂时不支持拼图
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表