grdselect

贡献者:

周茂

最近更新时间:

2022-12-11


官方文档:

grdselect

简介:

使用一定的空间标准对 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 子选项,将经纬度方向的间隔进行一定程度的舍入以应用到最终的区域。 lh 分别采用最小的网格间隔或者最大的网格间隔,或者可以指定间隔 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]

反转选择结果:

  • C 列出不包含 -C 设置的任何点的网格

  • D 列出和 -D 指定的分辨率不匹配的网格

  • F 列出不与 -F 指定多边形存在交集的网格

  • N 列出不满足 -N 选项的网格

  • R 列出不与 -R 相交的网格

  • W 列出不在 -W 范围内的网格

  • Z 列出不在 -Z 范围内的网格

  • r 列出不与 -rreg 指定的配准方式相同的网格

-Llinefile

指定一个多段线文件,列出至少被其中的一条线穿过的网格。如果输入网格为 3-D, 则忽略 z 方向

-Mmargins

扩展通过 -A 输出的网格范围,扩展大小为 margins ,其可以设置为单个值, 即在所有方向均使用相同大小的扩展;也可设置为斜杠分隔的两个值,即分别设置 x 和 y 方向的扩展;也可设置为斜杠分隔的 4 个值,分别设置西东南北四个方向 的扩展。默认不使用扩展。对于地理坐标网格,可使用 d , ms 分别表示度分秒。

-Nl|h[n]
  • h 只列出至少含有 n 个 NaN 值的网格

  • l 只列出最多含有 n 个 NaN 值的网格

n 默认为 0。该选项不能用于图片

-W[wmin]/[wmax]

只列出数据范围在 wminwmax 的网格,如果没给出 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

相关模块

gmtselect, grd2xyz, grdedit, grdinfo