grdfill
- 官方文档
- 简介
grdfill 模块读入一个文件数据，并向数据中的“洞”填充数据。“洞”通常指值为NaN的节点，但用户也可以使用其它准则指定“洞”。
语法
gmt grdfill ingrid -Amode[arg] -Goutgrid [ -Rregion ] [ -L[p] ] [ -V[level] ] [ -fflags ] [ --PAR=value ]
必须选项
- ingrid
输入网格文件
- -Amode[arg]
填充“洞”所使用的算法
目前支持两种算法：
cvalue 使用某个常数填充“洞”
n[radius] 使用最近的非NaN值填充“洞”，默认的搜索半径为半径默认取 \(r^2 = \sqrt{X^2 + Y^2}\)，其中 X 和 Y 分别是X方向和Y方向的网格数目。也可以指定 radius 为搜索半径（单位是节点数）。
- -Goutgrid
输出网格文件
可选选项
- -N[nodata]
所有值等于 nodata 的节点都被认为为“hole”，默认值为NaN
- -Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)
指定数据范围
该选项定义了要处理了子区域范围。
- -L[p]
不填充“洞”，仅列出每个“洞”所处的子区域的范围
-G 选项会被忽略。-Lp 表示输出每个子区域对应的闭合多边形。
- -V[level] (more …)
设置 verbose 等级 [w]
- -f[i|o]colinfo (more …)
指定输入或输出列的数据类型
- -^ 或 -
显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
- -+ 或 +
显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
- --PAR=value
临时修改GMT参数的值，可重复多次使用。参数列表见 gmt.conf
示例
检测网格文件中所有包含NaN的区域，并列出这些矩形区域的边界坐标:
gmt grdfill data.grd -L > wesn_listing.txt
检测网格文件中所有包含NaN的区域，并以多段文件的形式输出这些矩形区域对应的闭合多边形:
gmt grdfill data.grd -Lp > NaN_regions.txt
将网格文件中所有NaN值替换为999.0:
gmt grdfill data.grd -Ac999 -Gno_NaNs_data.grd
将网格文件中所有NaN值用最近的非NaN值替代:
gmt grdfill data.grd -An -Gno_NaNs_NN_data.grd