grdcut

官方文档

grdcut

简介

从一个网格文件中裁剪出一个子区域的网格文件

grdcut 模块可以从一个网格文件中根据条件裁剪出一个新的网格文件。子区域的选择有如下几种方式：

  • -R 选项：直接指定子区域的范围

  • -Z 选项：通过检查网格点的值间接限定子区域的范围

  • -S 选项：通过指定与特定点的距离间接限定子区域的范围

  • -J 和 -R 选项：对于倾斜投影，可以使用这两个选项获得覆盖整个倾斜区域的子区域

语法

gmt grdcut ingrid -Goutgrid -Rregion [ -Jparameters ] [ -N[nodata] ] [ -Slon/lat/radius[+n] ] [ -V[level] ] [ -Z[min/max][+n|N|r] ] [ -fflags ] [ --PAR=value ]

必选选项

ingrid

输入网格文件名

-Goutgrid

输出网格文件名

可选选项

-Jprojection (more …)

设置地图投影方式

-N[nodata]

允许新网格的区域范围超过原网格的区域范围。

一般情况下，若指定的区域范围大于输入网格的区域范围，超出的部分会被自动忽略，实际的输出网格的区域范围会自动适应输入网格的区域范围。使用 -N 选项，则超出的区域范围内，网格节点会被赋以指定的值，默认赋值为NaN，也可以自己指定其值为 nodata。

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

指定数据范围

指定要截取的网格子区域的范围。

-Slon/lat/radius[+n]

指定圆心位置（lon/lat）及其半径（radius，可指定 距离单位），程序会自动计算出一个包含了圆上及圆内所有网格点的矩形区域。 +n 表示将所有矩形区域内但不在圆内的节点的值设为NaN。

-V[level] (more …)

设置 verbose 等级 [w]

-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

距离单位

GMT支持多种不同的距离单位，以及三种不同的球面距离计算方式。详情见 单位 和 -j 选项。

示例

使用 -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_06m -Goutput.nc -S45/30/500k+n

对于一个倾斜投影的地图，由于 grdcut 获得的区域是由经线、纬度划分出来的矩形区域，而倾斜投影地图的边界不是由经纬线构成的，所以需要使用 -J 和 -R 选项指定一个更大的矩形区域以包含倾斜投影地图的研究区域:

gmt grdcut @earth_relief_02m -R160/20/220/30+r -Joc190/25.5/292/69/1 -Gdata.nc