GRIB文件是分段编码的自描述二进制格式,不可用fread直接读结构体;必须使用eccodes等专业库逐消息解析,注意handle管理、类型匹配与跨平台路径编码问题。

GRIB文件不是标准二进制,不能用fread直接读结构体
很多刚接触气象数据的人会尝试把GRIB当普通二进制文件处理:开个FILE*、fread进一个自定义struct——立刻失败。GRIB是分段(message)编码格式,每段含元数据+压缩的格点/谱数据,且支持多种编码(PNG、JPEG2000、complex packing),还带校验和。硬解析不仅工作量大,还极易因版本(GRIB1/GRIB2)、中心编号、模板编号差异出错。
实操建议:
- 必须用专业库,首选
eccodes(ECMWF官方维护,C/C++接口稳定,支持GRIB1/2、NetCDF混合读取) - 避免用过时的
grib_api(已归档,无新特性支持)或纯Python方案(如pygrib)在C++项目中桥接调用(性能损耗+异常传递复杂) - Windows下优先用vcpkg安装:
vcpkg install eccodes;Linux/macOS用包管理器或源码编译,注意启用ENABLE_FORTRAN=OFF减少依赖
用eccodes读单个GRIB消息:从codes_handle*开始
GRIB文件本质是多个独立消息拼接,eccodes不提供“读整个文件到内存”接口,而是逐消息迭代。核心流程是打开文件 → 获取下一个handle → 提取值 → 释放handle。
常见错误现象:codes_handle_new_from_file返回NULL但没检查;读完没调codes_handle_delete导致句柄泄漏(尤其循环读多消息时);用codes_get_double读整数型key(如numberOfPoints)不报错但值截断。
立即学习“C++免费学习笔记(深入)”;
关键步骤:
- 用
codes_open_file打开文件,模式必须是"rb"(文本模式在Windows下会误判换行符) - 循环调用
codes_handle_new_from_file,返回NULL表示无更多消息(不是错误) - 读数值优先用
codes_get_long(整数)或codes_get_double(浮点),不要依赖codes_get自动类型转换 - 必须为每个handle配对调用
codes_handle_delete,否则内存和内部缓存持续增长
示例片段(获取经纬度格点数):
long ni = 0, nj = 0;
int err = codes_get_long(handle, "Ni", &ni);
if (err == CODES_SUCCESS) {
codes_get_long(handle, "Nj", &nj);
}
读取格点数据数组:用codes_get_double_array而非逐点get
实际业务中,99%的需求是拿到double*形式的二维场(如温度、风速)。有人试图循环调用codes_get_double读values(i),结果极慢且可能触发内部重复解码(GRIB2常对values做二级压缩)。
正确做法是预分配足够内存,一次性提取:
- 先用
codes_get_size查"values"长度(返回元素个数,不是字节数) - 分配
double* data = new double[size],再调codes_get_double_array(handle, "values", data, &size) - 注意
size是输入输出参数:调用前填期望长度,返回时写入实际读取数;若传入长度小于实际值,函数返回CODES_BUFFER_TOO_SMALL - GRIB默认按“扫描线顺序”(即先行后列),对应数组索引
data[j * ni + i]为第j行第i列(0起始),与常见图像坐标系一致
跨平台路径与编码问题:文件名含中文或空格时fopen失败
Windows下用std::string构造路径传给codes_open_file,若路径含中文(如"D:/气象数据/202405.grib"),在UTF-8编译环境下fopen底层调用会失败(Windows API要求UTF-16)。这不是eccodes的bug,而是C标准库限制。
解决方案取决于构建环境:
- MSVC项目:改用
codes_open_file_utf8(需eccodes ≥ 2.27.0),它内部转UTF-16 - MinGW或跨平台代码:路径先用
std::filesystem::u8path转std::string(C++17),再确保编译时定义_CRT_SECURE_NO_WARNINGS并链接legacy_stdio_definitions.lib - 最稳方式:服务端部署统一用英文路径,客户端做映射(避免在生产环境碰编码边界)
容易被忽略的是,GRIB文件本身可能带BOM(尤其Windows生成工具导出),虽不影响eccodes解析,但会导致codes_open_file返回NULL(日志显示Cannot read from file)。用file -i filename.grib确认是否为application/octet-stream,而非text/plain。


















