ecCode库根据具体气象参数获取经纬度以及参数值

GIS开发小白

已于 2022-04-24 20:28:54 修改

阅读量819

点赞数 2

文章标签： c语言 c++ ubuntu

于 2022-04-24 20:27:53 首次发布

本文链接：https://blog.csdn.net/m0_58207883/article/details/124390776

版权

本文介绍了如何利用ecCodes库编译安装后的实例代码，通过指定气象参数ID，来获取格网点的经纬度及对应参数值。在阅读中，读者将学习到如何查找气象参数ID，并检查GRIB文件包含的参数，以便有效地提取所需气象数据。

摘要由CSDN通过智能技术生成

从github下载的ecCodes源码编译安装后，根据其提供的实例代码，调用接口，根据某一具体的气象参数来获取所有格网点的经纬度以及参数值。

具体代码如下:

#include <stdio.h>
#include <stdlib.h>

#include <eccodes.h>

#define MAX_VAL_LEN 1024 

int main(int argc, char** argv)
{
    int err              = 0;
    size_t i             = 0;
    FILE* in             = NULL;

    const char* filename = "/home/he/grib_data/gfs_data/gfs.t06z.sfluxgrbf103.grib2";

    codes_handle* h      = NULL;
    long numberOfPoints  = 0;
    const double missing = 9999.0;
    double *lats, *lons, *values; /* arrays */


    long paramId;
    double value;
    size_t vlen = MAX_VAL_LEN;
    codes_index* index;

    in = fopen(filename, "rb");
    if (!in) {
        fprintf(stderr, "Error: unable to open input file %s\n", filename);
        return 1;
    }

    /* create new handle from a message in a file */
    index = codes_index_new_from_file(0, filena