【python】CMPAS原始数据处理

Python pygrib库读取GRB2数据
以CMPAS数据为例

一、GRIB格式说明：

GRIB是世界气象组织（WMO）开发的一种用于交换和存储规则分布数据的二进制文件格式。最初GRIB表示“二进制格点”（GRIdded Binary），后来扩展为“二进制的通用规则分布信息”（General Regularly-distributed Information in Binary form）。

GRIB 是与计算机无关的压缩的二进制编码,主要用来表示数值天气预报的产品资料。现行的GRIB 码版本有GRIB1 和GRIB2 两种格式。 GRIB2较之GRIB1具有更多优点而被广泛使用。如：表示多维数据、模块性结构、支持多种压缩方式、IEEE标准浮点表示法等。

以cmpas降水资料为例子，该数据格式为grb2格式。

二、利用python进行读取

需要用到的库：pygrib
1.安装pygrib库，pip install pygrib安装不成功的话，可以使用
conda install pygrib
pip install pygrib -i 清华镜像源地址

pip install pygrib

2.开始读取

import pygrib
import os
import numpy as np
filen = r"E:\Z_SURF_C_BABJ_20230901001213_P_CMPA_RT_CHN_0P01_HOR-PRE-2023090100.GRB2"
grbs = pygrib.open(filen)
for grb in grbs:
    print(grb.name)
    lats,lons = grb.latlons()
    ## 对数据进行裁切，裁切范围：lat1~lat2，lon1~lon2
    data1, lats1, lons1 = grb.data(lat1=28, lat2=32, lon1=109, lon2=112)

3.读取后的结果，可以写入nc

import netCDF4 as nc
import numpy as np
import os

def create_nc(lat_value,lon_value,data,outname):
    # 创建一个新的数据集
    newdataset = nc.Dataset(outname,'w')
    # 创建维度 维度名，维度长度(默认None)
    nn = len(lat_value)
    mm = len(lon_value)
    newdataset.createDimension('lat',nn)
    newdataset.createDimension('lon',mm)
    # 创建变量 变量名，数据类型(numpy中的数据类型),维度(如果不给维度，则变量是标量)
    lat = newdataset.createVariable('lat',np.float32,('lat'))
    lon = newdataset.createVariable('lon',np.float32,('lon'))
    tp_temp = newdataset.createVariable('tp',np.float32,('lat','lon'))
    # 给变量赋值
    lat[:] = lat_value
    lon[:] = lon_value
    tp_temp[:] = data
    newdataset.close()