grdselect
- 贡献者:
- 最近更新时间:
2022-12-11
- 官方文档:
- 简介:
使用一定的空间标准对 2—D 网格、图像或 3-D 网格进行范围筛选或者确定并集
grdselect 读取一个或者多个图像或者网格,可对数据数据xxx
语法
gmt grdselect source1 source2 … [ -Ai|u[+il|h|inc] ] [ -Cpointfile ] [ -Dinc ] [ -E[b] ] [ -Fpolygonfile[+i|o] ] [ -G ] [ -I[C][D][F][L][N][R][W][Z][r] ] [ -Llinefile ] [ -Mmargins ] [ -Nl|h[n] ] [ -Rregion ] [ -V[level] ] [ -W[wmin/wmax] ] [ -Z[zmin/zmax]] ] [ -fflags ] [ -hheaders ] [ -oflags ] [ -rreg ] [ --PAR=value ]
必须选项
ingrid[=ID|?varname][+bband][+ddivisor][+ninvalid][+ooffset][+sscale]
输入网格名。通过追加 =ID 可指定 网格格式 [默认为 =nf];追加 ?varname 可指定 NetCDF 变量 [默认为 GMT 找到的第一个 2-D 网格]。对网格追加子选项可进行额外设置:
+b 选取一个波段 band(仅用于图片)[默认为 0]
+d 将网格除以一个数 divisor [默认为 1]
+n 将网格中的 invalid 替换为 NaN
+o 将网格中的值进行一定的偏移,即加一个数 offset [默认为 0]
+s 将网格中的值缩放,即乘以 scale [默认为 1]
注 :offset 操作位于 scale 之前。
备注
不可将 2—D 和 3—D 文件融合
可选选项
- -Ai|u[+il|h|inc]
i 求交集
u 求并集
可追加 +i 子选项,将经纬度方向的间隔进行一定程度的舍入以应用到最终的区域。 l 和 h 分别采用最小的网格间隔或者最大的网格间隔,或者可以指定间隔 inc 。如果不使用 -A ,则只列出所有通过以下选项 -C -D -F -L -N -R -W -Z 设置的测试的文件名,与 -I 同时使用,则列出未通过测试的文件名
- -Cpointfile
指定一个多段点文件 pointfile ,列出至少包含文件中一个点的网格。如果输入网格为 3-D,则会检查 pointfile 中的点是否包含 z 坐标,然后检查其是否位于网格中
- -Dinc
只列出网格间隔为 inc 的网格
- -E[b]
格式化输出数据,所有字段以制表符分隔置于同一行。其格式为: w e s n {b t} v0 v1 。只有输入数据为 3-D 网格时,才输出括号中的字段。默认输出 -Rw/e/s/n{/b/t} 字符串。追加 b 选项可输出区域范围的闭合多边形。v0 v1 为区域内所有数据的数值范围。如果输入网格的分辨率不同,则需要使用 -A+i 指定输出分辨率,若使用 Eb ,则无需 -A+i 选项
- -Fpolygonfile[+i|o]
指定一个多段闭合多边形文件。列出所有至少与其中一个多边形有交集的网格。使用 +i 和 +o 分别列出完全位于多边形内和外的网格。如果输入数据为 3-D,则忽略 z 方向。
- -G
为了列出所有信息,强制下载所有可能的的远程文件的区块
- -I[C][D][F][L][N][R][W][Z][r]
反转选择结果:
- -Llinefile
指定一个多段线文件,列出至少被其中的一条线穿过的网格。如果输入网格为 3-D,则忽略 z 方向
- -Mmargins
扩展通过 -A 输出的网格范围,扩展大小为 margins ,其可以设置为单个值,即在所有方向均使用相同大小的扩展;也可设置为斜杠分隔的两个值,即分别设置 x 和 y 方向的扩展;也可设置为斜杠分隔的 4 个值,分别设置西东南北四个方向的扩展。默认不使用扩展。对于地理坐标网格,可使用 d , m 或 s 分别表示度分秒。
- -Nl|h[n]
h 只列出至少含有 n 个 NaN 值的网格
l 只列出最多含有 n 个 NaN 值的网格
n 默认为 0。该选项不能用于图片
- -W[wmin]/[wmax]
只列出数据范围在 wmin 到 wmax 的网格,如果没给出 wmin, 则默认为 -infinity,没给出 wmax ,则默认为 +infinity
- -Z[zmin]/[zmax]
与 -W 相同,但用于 z 值
- -V[level] (more …)
设置 verbose 等级 [w]
- -f[i|o]colinfo (more …)
指定输入或输出列的数据类型
- -h[i|o][n][+c][+d][+msegheader][+rremark][+ttitle] (more …)
跳过或生成指定数目的头段记录
- -ocols[,…][,t[word]] (more …)
设置输出数据列(0表示第一列,t 表示文本列)
- -r[g|p] (more …)
设置网格配置方式 [默认为网格线配准]
- -^ 或 -
显示简短的帮助信息,包括模块简介和基本语法信息(Windows下只能使用 -)
- -+ 或 +
显示帮助信息,包括模块简介、基本语法以及模块特有选项的说明
- -? 或无参数
显示完整的帮助信息,包括模块简介、基本语法以及所有选项的说明
- --PAR=value
临时修改GMT参数的值,可重复多次使用。参数列表见 配置参数
示例
寻找所有网格的并集,并在最终结果中扩展 2 度:
gmt grdcut @earth_relief_01d -R1/3/1/3 -Gsubset1.nc
gmt grdcut @earth_relief_01d -R2/5/2/5 -Gsubset2.nc
gmt grdselect *.nc -Au -M2
寻找所有网格的交集:
gmt grdselect *.nc -Ai
寻找所有网格的交集,并将结果扩展 2 度,然后输出边界多边形:
gmt grdselect *.nc -Ai -M2 -Eb > wesn_polygon.txt
列出所有至少包含 10 个 NaN 值的像素配准方式的网格:
gmt grdselect *.nc -Nh10 -rp
列出完全位于 wesn_polygon.txt 中包含的多边形内的网格:
gmt grdselect *.nc -Fwesn_polygon.txt+i