python2 解析 Grib/Grib2数据

python2 解析 Grib/Grib2数据

sudo apt-get install gfortran cmake  zlib1g-dev  libpng-dev jasper libjpeg-dev libgrib-api-dev python-tk libgeos++-dev
pip install -i  numpy matplotlib  pyproj geos pygrib   eccodes

install jasper

mkdir build
cd build
cmake ../ -DCMAKE_INSTALL_PREFIX=/usr/local/jasper
make -j 64
sudo make install

install ecCodes

  tar -xzf eccodes-2.6.0-Source.tar.gz
  mkdir build
  cd build
  cmake  ../ -DCMAKE_INSTALL_PREFIX=/usr/local/eccodes/
  make -j 64
  make install

  1. 配置环境变量
  nano ~/.bashrc

  export ECCODES_HOME=/usr/local/eccodes
  export PATH


  source ~/.bashrc


  2. GRIB编解码的Python开发环境配置

    vi /home/ubuntu/.local/lib/python2.7/site-packages/gribapi.pth



  3. 此时导入gribapi测试,如未报错,则说明配置成功

    python -c "from gribapi import *"

测试能否读取 grib文件

#!/usr/bin/env python
# -*- coding:utf-
from eccodes import *

ifile = open('example.grib')
while 1:
    igrib = codes_grib_new_from_file(ifile)
    if igrib is None: break

    date = codes_get(igrib,"dataDate")
    levtype = codes_get(igrib,"typeOfLevel")
    level = codes_get(igrib,"level")
    values = codes_get_values(igrib)
    print (date,levtype,level,values[0],values[len(values)-1])



tar -xzf grib_api-1.28.0-Source.tar.gz
cd grib_api-1.28.0-Source
mkdir build
cd build
cmake  ../ -DCMAKE_INSTALL_PREFIX=/usr/local/grib_api/
make -j 64
make install



### 0611


pygrib-2.0.3]$ python3 install


## install basemap-1.0.7

<!-- tar zxvf basemap-1.0.7.tar.gz

cd basemap-1.0.7/geos-3.3.3/
export GEOS_DIR=/usr/local
./configure --prefix=$GEOS_DIR 
make -j 16
make install -->

pip install --user git+
<!-- Collecting git+ -->

python -c "from mpl_toolkits.basemap import Basemap"# 测试安装成功

cd examples/



  • read ECTHIN OK

    g++ main.cpp -L/usr/local/eccodes/lib/ -leccodes

    export LD_LIBRARY_PATH=/usr/local/eccodes/lib/:$LD_LIBRARY_PATH

    export PYTHONPATH=/usr/local/eccodes/lib/:$PYTHONPATH

 * Copyright 2005-2018 ECMWF.
 * This software is licensed under the terms of the Apache Licence Version 2.0
 * which can be obtained at
 * 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 = "/home/ubuntu/chenwei/C1D06031200060600001";
    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;