grdcut¶
| 官方文档: | grdcut |
|---|---|
| 简介: | 从一个网格文件中裁剪出一个子区域的网格文件 |
grdcut 模块可以从一个网格文件中根据条件裁剪出一个新的网格文件。子区域的选择有如下几种方式：
- -R 选项：直接指定子区域的范围
- -Z 选项：通过检查网格点的值间接限定子区域的范围
- -S 选项：通过指定与特定点的距离间接限定子区域的范围
- -J 和 -R 选项：对于倾斜投影，可以使用这两个选项获得覆盖整个倾斜区域的子区域
语法¶
gmt grdcut ingrid -Goutgrid -Rregion [ -Jparameters ] [ -N[nodata] ] [ -Slon/lat/radius[unit][+n] ] [ -V[level] ] [ -Z[min/max][+n|N|r] ] [ -fflags ] [ --PAR=value ]
可选选项¶
- -Jprojection (more …)
- 设置地图投影方式。
- -N[nodata]
允许新网格的区域范围超过原网格的区域范围。
一般情况下，若指定的区域范围大于输入网格的区域范围，超出的部分会被自动忽略，实际的输出网格的区域范围会自动适应输入网格的区域范围。使用 -N 选项，则超出的区域范围内，网格节点会被赋以指定的值，默认赋值为NaN，也可以自己指定其值为 nodata。
- -Rxmin/xmax/ymin/ymax[+r][+uunit] (more …)
- 指定数据范围
- -Slon/lat/radius[unit][+n]
- 指定圆心位置（lon/lat）及其半径（radius，单位为 unit），程序会自动计算一个矩形区域，该矩形区域包含了圆上及圆内所有网格点。 +n 表示将所有矩形区域内但不在圆内的节点的值设为NaN。
- -V[level] (more …)
- 设置 verbose 等级 [c]
- -Z[min/max][+n|N|r]
确定一个新的矩形区域，使得所有该区域外的节点值都在指定的Z值范围 min/max 外 min 和 max 默认值为正负无穷，可以用 - 表示无穷。
默认情况下，值为NaN的节点会被自动忽略。
- +n 表示将 NaN 节点视作在给定的Z值范围之外，则保证了新生成的子区域内是无NaN的
- +r 表示将 NaN 节点视作在给定的Z值范围之内，因而子区域的外部是无NaN的
- +N 若新生成的子区域的某条边界上的所有节点都是NaN，则去除该条边界上的这一行或列
- -f[i|o]colinfo (more …)
- 指定输入或输出列的数据类型
- -^ 或 -
- 显示简短的帮助信息，包括模块简介和基本语法信息（Windows下只能使用 -）
- -+ 或 +
- 显示帮助信息，包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
- 显示完整的帮助信息，包括模块简介、基本语法以及所有选项的说明
- --PAR=value
- 临时修改GMT参数的值，可重复多次使用。参数列表见 gmt.conf
示例¶
使用 -R 选项直接指定子区域的范围:
gmt grdcut @earth_relief_02m -Gregion_0_30_-30_30.nc -R0/30/-30/30
使用 -S 选项生成一个包含了原点(45,30)周围 500 km 以内的所有点的矩形区域，并设置矩形区域内圆外的节点值为NaN:
gmt grdcut @earth_relief_02m -Goutput.nc -S45/30/500k+n
使用 -Z 选项使得子区域外被除去的部分大于0:
gmt grdcut @earth_relief_02m -Goutput.nc -Z-/0
对于一个倾斜投影的地图，由于 grdcut 获得的区域是由经线、纬度划分出来的矩形区域，而倾斜投影地图的边界不是由经纬线构成的，所以需要使用 -J 和 -R 选项指定一个更大的矩形区域以包含倾斜投影地图的研究区域:
gmt grdcut @earth_relief_02m -R160/20/220/30+r -Joc190/25.5/292/69/1 -Gdata.nc