爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: eeeee

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

  [复制链接]

新浪微博达人勋

发表于 2019-9-5 08:24:37 | 显示全部楼层
eeeee 发表于 2019-9-4 16:46
这个不是我的模块的问题,可能是有模块的冲突

好的,谢谢答复
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-9-5 15:06:48 | 显示全部楼层
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')
QQ截图20190905150055.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-9-6 14:44:58 | 显示全部楼层
不不106 发表于 2019-9-5 15:06
你好,我在运行模板时,第一次是正常出图的,但第二遍时就会出现重新下载模板的错误。原程序如下,错误如 ...

问题已解决,谢谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-9-12 16:15:36 | 显示全部楼层
求助!
使用StandarData读取南京的S波段双偏振业务雷达,显示:
    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
    复制代码

尝试用了解压缩后的数据,显示:

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


读取单偏振的业务雷达,以及广州S波段双偏振业务雷达都没有问题。
这是为什么呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-9-24 16:28:13 | 显示全部楼层
请教老师,能否根据给定的某个经度、纬度,把该点的反射率 或者速度 值提取出来???
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2019-9-25 16:40:10 | 显示全部楼层
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)
密码修改失败请联系微信:mofangbao
回复 支持 3 反对 0

使用道具 举报

新浪微博达人勋

发表于 2019-10-29 08:43:12 | 显示全部楼层
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'
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-10-29 11:25:17 | 显示全部楼层
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 ...

QQ截图20191028202355.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-10-30 10:59:11 | 显示全部楼层
本帖最后由 qxjzkx@163.com 于 2019-11-14 17:21 编辑

底图省界还是有点问题,其他很好,牛B。
Z9531_20160614120050_0.0_230_VIL.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-10-30 11:31:47 | 显示全部楼层
qxjzkx@163.com 发表于 2019-10-30 10:59
绘图问题:
1、0到5dBZ白灰色区域偏多
2、强度偏弱,这是一次强对流天气过程

麻烦你看清楚这是VIL不是反射率
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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