13.29 grdsample

官方文档:grdsample 简介:对网格文件做重采样

该命令读取一个网格文件，并对其做插值以生成一个新的网格文件。新旧网格文件的区别在于：

1. 不同的配准方式（ -r 或 -T ）
2. 不同的网格间隔或网格节点数（ -I ）
3. 不同的网格范围（ -R ）

默认使用 bicubic 插值方式，可以使用 -n 选项设置其他插值方式。该命令可以很安全地将粗网格插值为细网格；反之，将细网格插值为粗网格时，则可能存在混淆效应，因而需要在插值前使用 grdfft 或 grdfilter 对网格文件做滤波。

13.29.1 必选选项

<in_grdfile>

要重采样的2D网格文件

-G<out_grdfile>

重采样生成的网格文件

13.29.2 可选选项

-I<xinc>[<unit>][+e|+n]/<yinc>[<unit>][+e|+n]

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

• <xinc> X方向的网格间隔
• <yinc> Y方向的网格间隔
• <unit> 网格间隔的单位。对于地理坐标系而言默认值为度，可以设置单位为 m|s|e|f|k|M|n|u 。
• +e 微调X和Y方向范围的最大值，使得其是网格间隔的整数倍（默认会微调网格间隔以适应给定的数据范围）
• +n 表明 <xinc> 和 <yinc> 不是网格间隔，而是X和Y方向的节点数。此时会根据节点数、网格区域范围以及网格配准方式重新计算网格间隔。

注意：

1. 若 <yinc> 设置为0，则表示其与 <xinc> 相同
2. 若使用 -R<grdfile> 选项，则网格网格已经根据网格文件自动初始化，此时依然可以使用 -I 覆盖网格文件中的间隔值

-R<w>/<e>/<s>/<n>

指定新网格的数据范围。

若只使用 -R 选项，则等效于使用 grdcut 或 grdedit -S 。

-T

交换网格文件的配准方式。即若输入是网格线配准，则输出为像素点配准；若输入是像素点配准，则输出为网格线配准。

-n[b|c|l|n][+a][+b<BC>][+c][+t<threshold>]

重采样时使用的插值算法，见 -n 选项 一节。

13.29.3 示例

将5x5弧分的数据采样成1x1弧分:

gmt grdsample hawaii_5by5_topo.nc -I1m -Ghawaii_1by1_topo.nc

将网格线配准的网格文件修改为像素配准的网格文件:

gmt grdsample surface.nc -T -Gpixel.nc