Question

官方文档

gmtinfo

简介

从表数据中提取极值信息

最小示例

提取文件中各列的极值:

$ gmt info input.dat
input.dat: N = 6    <1/2>   <1/3>   <0/3.2>

从输出中可以看到，输入文件中有6行3列，三列数据的范围分别是1到2、1到3和0到3.2。

可选选项

-Aa|t|s

输入数据为多段数据或多文件时的处理方式（默认为 -Aa ）：

1. -Aa ：报告所有文件的总最值

2. -Af ：分别报告各个文件自己的最值

3. -As ：分别报告各段数据自己的最值

-C

输出时，每列的最大最小值分别占一列

在输出每一列的最大最小值时，默认会用 <min/max> 的格式输出，不易于脚本处理。该选项会将每列的最大最小值分别输出为一列:

$ gmt info input.dat -C
1   2   1   3   0   3.2

该选项经常与 -i 及 -o 选项一起使用，以提取部分列的最值信息。

-D[dx[/dy]]

对 -I 选项的结果做整体偏移，使得区域范围的中心与原数据的中心对齐。

比如数据的真实范围是 0/3 ，使用 -I2 选项会得到 0/4 ，在此基础上使用 -D 选项，则变成 -0.5/3.5 。

-EL|l|H|h[col]

返回第 col 列最值所在的记录行。

1. l|h 表示返回最小、最大值所在的记录行

2. L|H 表示返回绝对值最小、最大值所在的记录行

若有多行均满足要求，则只返回第一行。若 col 未指定，则默认是最后一列。

-F[i|d|t]

返回记录统计数据

• -Fi 返回总文件数、总数据段数、总数据记录数、总头段记录数、总记录数

• -Fd 返回每个数据段的信息：表号、段号、行数、开始记录号、结束记录号

• -Ft 与 -Fd 类似，但会对每个输入文件重置段号、开始记录号和结束记录号

-I[b|e|f|p|s]dx[/dy[/dz…][+e|r|R]

调整最大/小值使得其是给定的增量的最近倍数，返回的格式为 -Rw/e/s/n

$ gmt info input.dat -I2/2
-R0/2/0/4

如果设置 <dx> 的值为 - ，则表示以 -Rw/e/s/n 的格式输出真实的 -R 范围:

$ gmt info input.dat -I-
-R1/2/1/3

当 -I 选项和 -C 选项一起使用时，输出不再是 -Rw/e/s/n 格式:

$ gmt info input.dat -C -I2
0   2   0   4   0   3.2

上面的例子中， -I 选项后只接了一个增量，此时第二列数据也会使用同样的增量（出于向后兼容性考虑）。若只想要为第一列指定增量，第二列不指定，可以使用 -Ip 选项:

$ gmt info input.dat -C -Ip2
0   2   1   3   0   3.2

对于某些算法，比如FFT或者surface命令， -R 选项的不同取值会影响到运算的速度， -If 会输出优化后的 -R 以得到最快FFT结果， -Is 会输出优化后的 -R 以得到最快的surface结果。

-L

确定多个表数据（ -Af ）或多段数据（ -As ）的共同范围。若与 -I 选项连用，则最终的结果会向内近似以使得极值位于真实数据范围内。

-S[x][y]

为绘制误差棒增加额外的列，与 -I 选项一起使用生成数据并直接用于接下来 plot -E 绘制误差棒。

• -Sx 将数据的第三列(2)作为水平误差棒的数据

• -Sy 将数据的第四列(3)作为垂直误差棒的数据

• -Sxy 将数据的第三(2)和第四(3)列作为水平和垂直误差棒的数据

-Tdz[w|d|h|m|s][+ccol]

以 -Tzmin/zmax/dz 的形式报告第 col 列的最值， col 的默认值为0。该选项不能与 -I 选项连用。