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

缥渺哈哈

mrsoberli 发表于 2020-9-17 17:33
修改下面两个文件即可，在里面增加雷达站的经纬度和站号信息；

pickle文件修改代码示例

非常感谢！！！

jack_ln

求教个问题。将pycinrad升级到1.6.1之后。

1. liner = fig.geoax.gridlines(draw_labels=True,linewidth=2, color='gray', alpha=0.5, linestyle='--')
   
2. liner.top_labels = False
   
3. liner.right_labels = False
   
4. liner.xformatter = LONGITUDE_FORMATTER
   
5. liner.yformatter = LATITUDE_FORMATTER

liner = fig.geoax.gridlines(draw_labels=True,linewidth=2, color='gray', alpha=0.5, linestyle='--')<br /> liner.top_labels = False<br /> liner.right_labels = False<br /> liner.xformatter = LONGITUDE_FORMATTER<br /> liner.yformatter = LATITUDE_FORMATTER

在图上绘制网格，网格绘制成功，但是标注经纬度，左侧和底部，都没有标注。不知道是什么原因，谢谢了。
cartopy 0.18.0
matplotlib 2.2.5

上面的问题解决了，修改安装目录visualize/ppi.py中的代码。

1. self.geoax.set_position([0.05, 0.2, 0.75, 0.8])
   
2. ax2 = self.fig.add_axes([0.05, 0.01, 0.75, 0.17])

复制代码

张小晚

eeeee 发表于 2019-1-7 19:42
斜体是模块底层部分
文件读取部分(cinrad.io)
目前cinrad模块支持SA/SB/SC/CA/CB/CC雷达以及新版本标准数 ...

楼主大大，我用cinrad.io.PUP来读取pup的37号组合反射率数据，但报错了。
脚本是这样的
from cinrad.io import PUP
f = PUP(r'D:\radar\data\20140521\pup_data\CR\37\20140521.090000.00.37.200')
data = f.get_data()
print(data)
报错内容：
File "E:\anaconda\lib\site-packages\requests\adapters.py", line 516, in send
    raise ConnectionError(e, request=request)

ConnectionError: HTTPSConnectionPool(host='raw.githubusercontent.com', port=443): Max retries exceeded with url: /Unidata/MetPy/v0.12.1/staticdata/sfstns.tbl (Caused by NewConnectionError('<urllib3.connection.VerifiedHTTPSConnection object at 0x000001F8EAFE8F48>: Failed to establish a new connection: [Errno 11004] getaddrinfo failed'))
我想请问一下，我哪一步开始错了呀！！！

gongxi12315

eeeee 发表于 2019-1-7 19:42
斜体是模块底层部分
文件读取部分(cinrad.io)
目前cinrad模块支持SA/SB/SC/CA/CB/CC雷达以及新版本标准数 ...

好厉害，请问大神你这里任意两点的剖面图生成要多少秒啊？

jack_ln

楼主好，
用plot_cross_section绘制剖面图如下图：
用section绘制剖面图如下图:
section的图像在高层缺失了，plot_cross_section的图像没有这个问题，想问问怎么解决section高层缺失的问题呢？谢谢

淡水河

标准格式的37号组合反射率产品是不是不支持呢？

伊萱

花容天下重莲 发表于 2020-9-9 21:23
老师想问下为什么双击main_ui.pyw没有反应，进不到GUI界面？

先运行一个用来编写py的软件，再open main_ui.pyw

伊萱

求助：运行报错“Only VPPI scan mode is supported”，要怎么解决？非常感谢

相关程序：
header = np.frombuffer(f.read(1024), CC_header)
scan_mode = header["ucScanMode"][0]
        if scan_mode < 100:
            raise NotImplementedError("Only VPPI scan mode is supported")

其中header 值为：
[(b'RSTM\x01 \x01 \x01       ', b'                Z9691   XiShua', b'
ngBanNa_Z9691       ', b"      \x02+\xacA'\xf1\xcaB?  s\x05   P\xafE{\x14n?fff?\x06+  #  ", b'          ', b'                    ', b'                ', b'       VCP31    ', 538976288, 538976288, 538976288, 8224, 8224, 32, 32, 32, 32, 32, 32, 32, 86, 67, 80, 51, 49, 32, 32, 32, 32, 538976288, 8224, 8224, 8224, 8224, 8224, 32, [32, 32, 32], [ 32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,   3,  32,  32,  32,  32,  32,  32,  32,  63,  32,  32,  97, 185, 123,  95,   7,  32,  32,  32, 236,  81, 162,  65,  32,  32,  32,  32, 225, 122, 166,  65,  32,  32,  32,  32, 219, 112,   2,  67,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,   1,  32,  32,  32,  32,  32,  32,  32,  32,  32, 175,  67,  32,  32,  32,  32,   1,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  63,  29,  59, 179,  67, 121,  51, 222,  64,  32,  32, 128,  63,  32,  32,  64,  64,  63,  32,  32,  32,  63,  32,  32,  32,  72,  27,   6,  32,  32,  32,  32,  32,  32,  32,  32,  32, 120,  32,  32,  32,   8,  32,  32,  32,   1,  32,  32,  32, 111,  18,  63, 141, 187, 149,  64,  30,  32,  32,  32,  32,  32,  32,  32,  30,  32,  32,  32,  32,  32,  32,  32,  14,  32,  32,  32, 154, 153,  63,  32,  32,  32,  63,  32,  32,  32,  66,  32,  32,  63,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,   8,  32,  32,  32,  12,  32,  32,  32,  14,  32,  32,  32,  14,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,   1,  32,  32,  32,   2,  32,   5,  32,   1,  32,  32,  32,  50,  32,  50,  32,  63,  32,  63,  32,  63,  32,  94,   1,  63,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32, 124,   1, 120,   1, 118,   1, 119,   1, 123,   1, 123,   1, 124,   1,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,   1,  32,  32,  32,   1,  32,  32,  32,  32,  32,  97,  68,  32,  32,  32,  32,   1,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  63, 121,  54, 158,  65, 198, 140, 177,  65,  32,  32, 128,  63,  52,  51, 211,  64,  63,  32,  32,  32,  63,  32,  32,  32,  72,  27,   6,  32,  32,  32,  32,  32,  32,  32,  32,  32, 128,  32,  32,  32,   8,  32,  32,  32,   1,  32,  32,  32, 111,  18,  63, 108, 131,  64,  65,  30,  32,  32,  32,  32,  32,  32,  32,  30,  32,  32,  32,  32,  32,  32,  32,  14,  32,  32,  32, 154, 153,  63,  32,  32,  32,  63,  32,  32,  32,  66,  32,  32,  63,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,   8,  32,  32,  32,  12,  32,  32,  32,  14,  32,  32,  32,  14,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,   1,  32,  32,  32,   2,  32,   5,  32,   1,  32,  32,  32,  50,  32,  50,  32,  63,  32,  63,  32,  63,  32,  94,   1,  63,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32,  32], [32, 32], 24911904, 24510840, 24838519, 379, 380, 8224, 8224, 8224, 8224, 8224, 8224, 8224, 8224, 8224, 32, b' ', 32, 32, 1, 32, 32, 32, 32, 32, 32, 32, 32, [ 32, 175,  67,  32,  32,  32,  32,   1,  32,  32,  32,  32,  32,  32,  32,  32,  32,  63,  74, 153,  82,  66,  63,  91,  66,  32,  32, 128,  63,  32,  32,  64,  64,  63,  32,  32,  32,  63,  32,  32,  32,  72,  27,   6,  32,  32,  32,  32,  32,  32,  32,  32,  32, 120,  32,  32,  32,   8,  32,  32,  32,   1,  32,  32,  32, 111,  18,  63, 141, 187, 149,  64,  30,  32,  32,  32,  32,  32,  32,  32,  30,  32,  32,  32,  32,  32,  32,  32,  14,  32,  32,  32], 154, 538984345)]

伊萱

找到原因了，转换日期这里写漏了几个[0]：
    def _CC_handler(self, f: Any):
        header = np.frombuffer(f.read(1024), CC_header)
        scan_mode = header["ucScanMode"][0]
        print('运行  _CC_handler')
        print(header["ucScanMode"][0])
        if scan_mode < 100:
            raise NotImplementedError("Only VPPI scan mode is supported")
        stop_angle = scan_mode - 100
        self.scantime = (
            datetime.datetime(
                header["ucEYear1"][0] * 100 + header["ucEYear2"][0],
                header["ucEMonth"][0],
                header["ucEDay"][0],
                header["ucEHour"][0],
                header["ucEMinute"][0],
                header["ucESecond"][0],
            )
            - utc_offset
        )

伊萱

楼主好，用你的程序成功画出了SA和CC雷达图，非常感谢分享。现在想画我们新建成的一部CC雷达，是C波段双偏振雷达，请问适用吗？