当前位置：首页 > news >正文

python绘制CPEFS人影区域模式预报产品图

news 来源：原创 2025/6/1 23:12:47

模式说明

CMA-CPEFS 自然云预报产品。
1.自然云预报产品是指未施加任何云催化影响下的数值预报产品，用于每日定时预报自然状态下的云和降水发展演变
2.产品数据采用GrADS数据格式存储方式，由一个GrADS文本格式的数据描述文件和多个GrADS二进制格式的数据存储文件组成
3.数据分发时使用bzip2软件进行压缩，压缩后数据文件名以“.dat.xz”结尾，具体命名规则如下：
Z_NAFP_C_BABJ_YYYYMMDDHH_P_WMC-CPEFS-AREA_FF.dat.xz
4.产品包含三维预报和二维预报产品：
温、压、湿、风，位势高度、水汽混合比、云水、雨水混合比、雨滴数浓度、冰晶比含水量、雪比含水量、霰比含水量、雷达反射率、云量等多种。
5.用途包含，基础产品图、人影作业建模、作业潜力区识别、云雾仿真、气象算法，产品融合

绘制产品

以下案例代码使用云顶温度绘制，其他产品绘制代码过大不展示可以绘制

fileName = os.path.basename(file_path)time = fileName.split('_')[4]moment = fileName.split('_')[-1].split('.')[0]with open(file_path, 'r') as f:data1 = np.fromfile(f, dtype=np.dtype('>f')) data = data1.reshape((int(891), int(595), 356), order='F') # data[np.where(data == 1.e+30)] = 0 vardata = data[:, :, 352].T  # np.place(rainnc, rainnc == 1.e+30, np.nan)vardata[np.where(vardata == 1.e+30)] = 0flipped_data = np.flipud(vardata)custom_colors1 = ['#FFFFFF', '#D588FD', '#9600B4', '#C8649B', '#FFB4B4', '#FF6464', '#FF2828','#FCCE32', '#E3C048', '#EEF74A', '#69EB38', '#06C80C','#71E8A6', '#5E9CFB', '#D0CBF3', '#E9EEF5']cmap = colors.ListedColormap(custom_colors1)plt.plot(flipped_data, alpha=0)# 绘制填色图plt.imshow(flipped_data, cmap=cmap, interpolation='nearest')# 保存图形plt.savefig('a.png', bbox_inches='tight', pad_inches=0, transparent=True)