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

eeeee

pixi 发表于 2019-7-8 20:01
对雷达数据不甚至了解，只知道个大概。请教一下，比如我想取某个射线路径方向的反射率因子，按照帖子的示例 ...

通过方位角索引出每个仰角所对应的那条径向，然后拼在一起

eeeee

dongfangxing 发表于 2019-7-5 12:43
@楼主，谢谢共享，推广python处理雷达数据，我按照您的方法全部安装好了，应用你写的代码（PPI.py)运行时， ...

cartopy的这个功能还没有在最新的release当中包含，可以去GitHub更新一下，或者等cartopy发布下一个版本

jack_ln

eeeee 发表于 2019-6-15 23:52
换个中文字体就可以了

谢谢指点，
具体操作
Lib\site-packages\cinrad-1.3.3-py3.7-win32.egg\cinrad\visualize目录下，修改__init__.py里面的字体设置就可以。

jack_ln

@eeeee 你好，
画图的过程中遇到下面的问题，
1,画REF，ele = 0时画图正常，ele = 1时提示错误如下图，生成的图片没有数据。
2、画VEL，ele = 0提示错误如下图，生成的图片没有数据；ele = 1时画图正常
ele = 0时，Elevation angle: 0.615234375
ele = 1时，Elevation angle: 0.52734375
是基数据的问题造成的吗？忘大神给指导指导

jack_ln

@eeeee 老师好，
用您的这个程序处理SC雷达数据时，提示下面的错误。
  File "C:\Program Files\Python37-32\lib\site-packages\cinrad-1.3.3-py3.7-win32.egg\cinrad\io\io.py"
, line 144, in __init__
    self._CD_handler(f)
  File "C:\Program Files\Python37-32\lib\site-packages\cinrad-1.3.3-py3.7-win32.egg\cinrad\io\io.py"
, line 237, in _CD_handler
    header['obs']['ssecond'][0]) - utc_offset
ValueError: year 13886 is out of range
提示日期处理错误。
io.py中这段提示出错
    def _CD_handler(self, f:Any):
        header = np.frombuffer(f.read(CD_dtype.itemsize), CD_dtype)
        el_num = header['obs']['stype'][0] - 100 # VOL
        if el_num == 5:
            self.task_name = 'VCP31'
        elif el_num == 9:
            self.task_name = 'VCP21'
        elif el_num == 14:
            self.task_name = 'VCP11'
        self.scantime = datetime.datetime(header['obs']['syear'][0], header['obs']['smonth'][0], header['obs']['sday'][0],
                                          header['obs']['shour'][0], header['obs']['sminute'][0],
                                          header['obs']['ssecond'][0]) - utc_offset






尝试修改了一下cinrad-1.3.3-py3.7-win32.egg\cinrad\io\_radar_struct\CD.py中
observation_param = [("stype", "B"),
                     ("syear", "u2"),
                     ("smonth", "B"),
                     ("sday", "B"),
                     ("shour", "B"),
                     ("sminute", "B"),
                     ("ssecond", "B"),
                     ("Timep", "B"),
                     ("smillisecond", "u4"),
                     ("calibration", "B"),
                     ("intensityI", "B"),
                     ("VelocityP_1", "B"),
                     ("layerparam", layer_param_dtype, 30),
                     ("RHIA", "u2"),
                     ("RHIL", "i2"),
                     ("RHIH", "i2"),
                     ("Eyear", "u2"),
                     ("Emonth", "B"),
                     ("Eday", "B"),
                     ("Ehour", "B"),
                     ("Eminute", "B"),
                     ("Esecond", "B"),
                     ("Etenth", "B")]
u2修改为B之后，又出现month 必须在1-12之间。不知道是不是这里的错误。


数据名称Z_RADR_I_Z****_20190703091841_O_DOR_SC_CAP.bin，望指点下这个错误怎么修改。谢谢

eeeee

jack_ln 发表于 2019-7-14 11:20
@eeeee 你好，
画图的过程中遇到下面的问题，
1,画REF，ele = 0时画图正常，ele = 1时提示错误如下图，生 ...

最低的两个仰角分别扫两次，0和2是反射率，1和3是速度

eeeee

jack_ln 发表于 2019-7-14 15:59
@eeeee 老师好，
用您的这个程序处理SC雷达数据时，提示下面的错误。
  File "C:\Program Files\Python37 ...

SC雷达理论上不会出现这个错误的诶，我不太清楚

jack_ln

eeeee 发表于 2019-7-15 10:30
最低的两个仰角分别扫两次，0和2是反射率，1和3是速度

明白了，谢谢指点

jack_ln

eeeee 发表于 2019-7-15 10:34
SC雷达理论上不会出现这个错误的诶，我不太清楚

谢谢，在你的指导下，明白了是我数据的问题

skyinglll

老师好，请问一下，我安装了这个GUI界面，为什么只能显示REF雷达反射率，而不能显示VEL呢？我想看看那个径向移动速度，请问要怎么做呢？