python绘制CPEFS人影区域模式预报产品图-CSDN博客

本文链接：https://blog.csdn.net/weixin_45863084/article/details/147311122

模式说明

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)] = 0
        flipped_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)