PyCINRAD读取探测中心天气雷达拼图系统v3产品

PyCINRAD  

PyCINRAD号称国内最强大的雷达数据处理工具，支持多种CINRAD雷达数据格式，兼容性强，横扫CMA雷达基数据和产品，国内首个支持“探测中心拼图3.0”。快速提取、计算雷达衍生产品，如组合反射率CR、回波顶高ET、降水粒子分类HCL等。作为一个开源项目，PyCINRAD对所有用户免费开放，鼓励社区贡献和改进。据悉已遍布CMA的各个系统角落。    

PyCINRAD官网 https://pycinrad.cn

本文介绍如何使用PyCINRAD读取：探测中心天气雷达拼图系统v3产品

import warnings
import cinrad
import cartopy
from cinrad.visualize import PPI
import numpy as np
import matplotlib.pyplot as plt
import datetime
%matplotlib inline
warnings.filterwarnings("ignore")
 
 
basePath = "d:/temp/"
print(cinrad.__version__)
print(np.__version__)
print(cartopy.__version__)
1.9.1
2.0.2
0.23.0

1、打开数据

探测中心天气雷达拼图系统v3产品

# MocMosaic
nFiles = (basePath+"cinrad/bz2/Z_RADA_C_BABJ_20240925080617_P_DOR_ACHN_CREF_20240925_080000.bin")
f = cinrad.io.MocMosaic(nFiles)
dt = f.get_data()
dt

 Size: 208MB
Dimensions:    (latitude: 4200, longitude: 6200)
Coordinates:
  * latitude   (latitude) float64 34kB 12.2 12.21 12.22 ... 54.18 54.19 54.2
  * longitude  (longitude) float64 50kB 73.0 73.01 73.02 ... 135.0 135.0 135.0
Data variables:
    CR         (latitude, longitude) float64 208MB nan nan nan ... nan nan nan
Attributes:
    scan_time:        2024-09-25 08:00:00
    site_code:        MOC
    site_name:        MOC
    tangential_reso:  nan
    range:            nan
    elevation:        0

2、画图 

# dt = dt.rename({"CREF": "CR"})  # 重命名，因为PPI并不支持CAP.v1.9.1以下需要这句
fig = PPI(dt, style="black", extent=[105, 110, 30, 35], add_city_names=True)

 

 PyCINRAD有丰富的接口，看接口名称应该都能看懂吧

cinrad.io.StandardData    
cinrad.io.StandardPUP    
cinrad.io.MocMosaic    
cinrad.io.SWAN(filename,product="CR")    
cinrad.io.CinradReader(filename,radar_type="CC")        
cinrad.io.PhasedArrayData    
cinrad.io.PUP    
cinrad.io.read_auto

#1. 标准格式基数据StandardData

标准格式基数据: 标准格式文件名中有“FMT”、基数据是O_DOR Z_RADR_I_Z9999_20231025220414_O_DOR_SA_CAP_FMT.bin

#2. 标准格式产品StandardPUP

标准格式产品: 标准格式文件名中有“FMT”、产品数据是P_DOR 这里文件名中有HCL，代表是HCL水凝物分类产品 Z_RADR_I_Z9999_20231231160439_P_DOR_SAD_HCL_250_230_5_FMT.bin 下面这种老格式的也是标准格式产品 Z9737_20231126150649Z_CR_00_37

#3. 探测中心天气雷达拼图系统v3产品.MocMosaic

ACHN表示是全网拼图，CREF表示是产品类型是CREF组合反射率 Z_RADA_C_BABJ_20231212010615_P_DOR_ACHN_CREF_20231212_010000

#4. SWAN

以前的老SWAN3产品.MCR表示拼图组合反射率 Z_OTHE_RADAMCR_20220525085400.bin

#5. 老格式的国产基数据.CinradReader

和标准格式的唯一区别是文件名中没有“FMT”（这里指的是大部分新情况） Z_RADR_I_Z9731_20240511070524_O_DOR_SAD_CAP.bin

#6. 相控阵PhasedArrayData

标准格式的相控阵雷达基数据(2023以后). 文件名一般有AXPT，就是纳睿雷达的产品。 Z_RADR_I_ZBJ02_20210815155836_O_DOR_DXK_CAR
Z_RADR_I_ZA601_20240415183600_O_DOR_AXPT0364_CRA_FMT
Z_RADR_I_ZS999_20231220000000_O_DOR_AXPT0364_CRA_.bz2

#7: PUP

NEXRAD Level 3 (NIDS) product files. 已经被淘汰的格式，一言难尽，格式太乱，能用则用吧。

#万能读取read_auto

如果你搞不清楚，那请用cinrad.io.read_auto(your_radar_file)来读取。