解析grib格式数据库eccodes的安装

解析grib格式数据库eccodes的安装

• 一、说明
• 二、安装eccodes库
• 三、测试代码

一、说明

1、环境：在rethat7.4下使用eccodes-2.8.0解析grib数据。
2、eccodes下载地址：https://confluence.ecmwf.int//display/ECC/Releases
3、需要cmake编译安装，cmake下载地址：https://cmake.org/download/

二、安装eccodes库

1、解压cmake
gunzip cmake-3.12.0.tar.gz
tar -xvf cmake-3.12.0.tar

2、安装cmake
cd cmake-3.12.0/
./bootstrap
make
make install
如果不安装cmake的话，会提示下面的错误：

3、解压eccodes库
gunzip eccodes-2.8.0-Source.tar.gz
tar -xvf eccodes-2.8.0-Source.tar

4、安装eccodes库
cd eccodes-2.8.0-Source/
mkdir build
cd build
cmake …/…/eccodes-2.8.0-Source -DCMAKE_INSTALL_PREFIX=/usr/local/lib
如果提示找不到libgfortran.so.3链接库，可以使用下面的命令：
cmake -DCMAKE_INSTALL_PREFIX=/usr/local/lib …/…/eccodes-2.8.0-Source -DENABLE_FORTRAN=OFF

编译后正确的截图如下：

然后 make make install 即可

三、测试代码

测试eccodes库的使用，使用的是官方的例子

/*
 * Copyright 2005-2018 ECMWF.
 *
 * This software is licensed under the terms of the Apache Licence Version 2.0
 * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
 *
 * In applying this licence, ECMWF does not waive the privileges and immunities granted to it by
 * virtue of its status as an intergovernmental organisation nor does it submit to any jurisdiction.
 */

/*
 * C Implementation: grib_get_data
 *
 * Description: how to get lat/lon/values from a GRIB message
 *
 */
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include "eccodes.h"

using namespace std;

int main (int argc, char **argv)
{
    int err = 0;
    size_t i = 0;
    FILE *in = NULL;
    const char *filename = "/mnt/hgfs/winShareUbuntu/reduced_latlon_surface.grib1";
    codes_handle *h = NULL;
    long numberOfPoints = 0;
    const double missing = 9999.0;
    double *lats, *lons, *values;       /* arrays */

    in = fopen (filename, "r");
    if (!in) {
        printf ("ERROR: unable to open input file %s\n", filename);
        return 1;
    }

    /* create new handle from a message in a file */
    h = codes_handle_new_from_file (0, in, PRODUCT_GRIB, &err);
    if (h == NULL) {
        printf ("Error: unable to create handle from file %s\n", filename);
        return 1;
    }

    CODES_CHECK (codes_get_long (h, "numberOfPoints", &numberOfPoints), 0);
    CODES_CHECK (codes_set_double (h, "missingValue", missing), 0);

    lats = (double *) malloc (numberOfPoints * sizeof (double));
    if (!lats) {
        printf ("unable to allocate %ld bytes\n", (long) (numberOfPoints * sizeof (double)));
        return 1;
    }
    lons = (double *) malloc (numberOfPoints * sizeof (double));
    if (!lons) {
        printf ("unable to allocate %ld bytes\n", (long) (numberOfPoints * sizeof (double)));
        free (lats);
        return 1;
    }
    values = (double *) malloc (numberOfPoints * sizeof (double));
    if (!values) {
        printf ("unable to allocate %ld bytes\n", (long) (numberOfPoints * sizeof (double)));
        free (lats);
        free (lons);
        return 1;
    }

    CODES_CHECK (codes_grib_get_data (h, lats, lons, values), 0);

    for (i = 0; i < numberOfPoints; ++i) {
        if (values[i] != missing) {
            printf ("%f %f %f\n", lats[i], lons[i], values[i]);
        }
  }

    free (lats);
    free (lons);
    free (values);
    codes_handle_delete (h);

    fclose (in);
    return 0;
}

运行的结果截图：

说明：最近才开始学习ecodes库，也在继续学习中，待出成果再继续分享。。。