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）：-Aa：报告所有文件的总最值-Af：分别报告各个文件自己的最值-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 列最值所在的记录行。
l|h表示返回最小、最大值所在的记录行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 选项连用。