Question

官方文档

grdblend

简介

将多个部分重叠的网格文件合并成一个网格文件

grdblend 模块会读取多个网格文件以及一个合并参数文件，并将多个网格文件合并成一个网格文件。

合并过程中会遇到如下几种情况：

• 某个节点只有一个值：使用该值填充该节点

• 某个节点有多个值：对多个值计算其权重均值并赋值给该节点

• 某个节点没有值：该节点的值由

  gmt grdblend [ blendfile | grid1 grid2 … ]

  grid1 grid2

  要合并的网格文件

  -Goutgrid

  输出的网格文件名

  -Ixinc[+e|n][/yinc[+e|n]]

  指定X和Y方向的网格间隔

  • xinc 和 yinc 为 X 和 Y 方向的网格间隔。对于地理坐标，可以指定网格间隔单位 [默认单位为度]

  • +e 微调X和Y方向范围的最大值，使得其是网格间隔的整数倍（默认会微调网格间隔以适应给定的数据范围）

  • +n 表明 xinc 和 yinc 不是网格间隔，而是X和Y方向的节点数。此时会根据节点数、网格区域范围以及网格配准方式重新计算网格间隔。

  注意：

  • 若 yinc 设置为0，则表示其与 xinc 相同

  • 若使用 -Rgrdfile 选项，则网格间隔和配准方式已经根据网格文件自动初始化，此时依然可以使用 -I 和 -r 覆盖相应的值

  -Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)

  指定数据范围

  可选选项

  blendfile

  合并参数文件。

  文件中包含了要合并的网格文件的文件名列表，其格式为:

  网格文件名（必须）  -R选项（可选）  相对权重（可选）
  

  • 在考虑权重的情况下，-R 范围外的部分会给零权重，-R 范围内的部分则会使用指定的相对权重，在边界处会加上2D余弦taper权重。

  • 若相对权重为负值，则整个反过来，即范围内的会给零权重，范围外的会给权重的绝对值。

  • 假设要合并几个网格文件，则可以设置合并参数文件 blend.job 的内容为:

    piece_1.nc -R0/30/-90/90 1
    piece_2.nc -R25/50/-90/90 1.5
    piece_3.nc -R45/80/-90/90 0.9
    piece_4.nc -R80/160/-90/90 1
    

    执行如下命令即可实现数据合并:

    gmt grdblend blend.job -Gblend.nc -R0/160/-90/90 -I1m/1m -V
    

    将所有网格文件 MB_*.nc 以相同权重合并:

    gmt grdblend MB_*.nc -Gblend.nc -R0/360/-90/90 -I1m/1m -V
    

    相关模块

    grd2xyz, grdconvert, grdedit, grdsample