gmtinfo

官方文件: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|f|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

對 -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[p|f|s]<dx>[/<dy>[/<dz>...]]

調整最大/小值使得其是給定的增量的最近倍數，返回的格式爲 -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[xy]

爲繪製誤差棒增加額外的列，與 -I 選項一起使用生成數據並直接用於接下來 plot -E 繪製誤差棒。

  • -Sx 將數據的第三列(2)作爲水平誤差棒的數據
  • -Sy 將數據的第四列(3)作爲垂直誤差棒的數據
  • -Sxy 將數據的第三(2)和第四(3)列作爲水平和垂直誤差棒的數據
-T<dz>[+c<col>]
以 -Tzmin/zmax/dz 的形式報告第 <col> 列的最值， <col> 的默認值爲0。該選項不能與 -I 選項連用。