title: 常数场（Constant fields）

weight: 2

Constant fields

如果要素场的所有数值都相等，那么使用常规的 GRIB 消息会将同一个数值重复存储 N 次，效率非常低。

ecCodes 使用常数场保存这种类型的数据，该常数应该是唯一保存的数值，数据段应该为空，因此常数场非常小，并且被精确地编码。

使用如下命令可以创建一个常数场：

  1. grib_set -d 1 in.grib out.grib

常数场中未定义打包参数（packing parameters），即 bitsPerValue=0。

常数场的问题

加载一个常数场：

  1. codes_grib_new_from_file(infile, igrib)

设置一些非常数的数据值：

  1. codes_set(igrib, 'values', values)

注意：此时还不知道打包参数

将消息写入到文件中：

  1. codes_write(igrib, outfile)

我们会得到什么样的 packingError？

因为常数场未设置打包参数，因此 ecCodes 不知道我们需要什么样的小数精度。 一个安全的选择是设置 bitsPerValue = 24。

更好的实践是在设置数值前设置 decimalPrecision 或 bitsPerValue。

  1. codes_grib_new_from_file(infile,igrib)
  2. codes_set(igrib,’decimalPrecision’,4)
  3. codes_set(igrib,’values’,values)
  4. codes_write(igrib,outfile)
  1. codes_grib_new_from_file(infile,igrib)
  2. codes_set(igrib,’bitsPerValue’,16)
  3. codes_set(igrib,’values’,values)
  4. codes_write(igrib,outfile)

练习

  1. 创建一个常数 GRIB 文件。

    1.  $ grib_set -d 1 t.grib2 c.grib2
    2.  $ grib_set -s packingType=grid_simple c.grib2 c1.grib2
    3.  $ grib_dump c1.grib2 
    4.  ***** FILE: c1.grib2 
    5.  #==============   MESSAGE 1 ( length=179 )                 ==============
    6.  GRIB {
    7.    # Meteorological products (grib2/tables/4/0.0.table)  
    8.    discipline = 0;
    9.    editionNumber = 2;
    10.    ...
    11.    values(1036800) =  {
    12.    1, 1, 1, 1, 1, 
    13.    1, 1, 1, 1, 1, 
    14.    1, 1, 1, 1, 1, 
    15.    ...
    16.    } 
    17.    #-READ ONLY- maximum = 1;
    18.    #-READ ONLY- minimum = 1;
    19.    #-READ ONLY- average = 1;
    20.    #-READ ONLY- numberOfMissing = 0;
    21.    #-READ ONLY- standardDeviation = 0;
    22.    #-READ ONLY- skewness = 0;
    23.    #-READ ONLY- kurtosis = 0;
    24.    #-READ ONLY- isConstant = 1;
    25.    #-READ ONLY- getNumberOfValues = 1036800;
    26.  }

  2. 设置 values={23.26, 42.51, 61.22, 45.95}，打印 packingError 和 bitsPerValue。

    创建一个 filter 规则文件 step02.filter。

    1.  $ cat step02.filter 
    2.  set values={23.26, 42.51, 61.22, 45.95};
    3.  write;
    4.  print "packingError = [packingError]";
    5.  print "bitsPerValue = [bitsPerValue]";

    执行命令：

    1.  $ grib_filter -o step02.grib2 step02.filter c1.grib2
    2.  packingError = 0.00500122070312
    3.  bitsPerValue = 12

译注：此处应该无法读取 packingError 值，但实际运行结果可以显示，可能我用的方式有问题。

  1. 设置 decimalPrecision=1，再设置相同的数据值，再一次打印 packingError 和 bitsPerValue。

    创建一个规则文件 step03.filter。

    1.  $ cat step03.filter 
    2.  set decimalPrecision=1;
    3.  set values={23.26, 42.51, 61.22, 45.95};
    4.  write;
    5.  print "packingError = [packingError]";
    6.  print "bitsPerValue = [bitsPerValue]";

    执行命令：

    1.  $ grib_filter -o step03.grib2 step03.filter c1.grib2
    2.  packingError = 0.0500007629395
    3.  bitsPerValue = 9

  2. 比较两个文件的大小和 packingError。

    两个文件大小相似，但 packingError 不同，见步骤 2 和 3 的输出。

    1.  185 step02.grib2
    2.  184 step03.grib2