一個 GMT 命令由“gmt + 模塊 + 選項(xiàng) + 參數(shù)”構(gòu)成,寫成如下形式: gmt module -Axx+bxxxx -Bxx+axxxx? gmt 是 GMT 中“唯一”的一個二進(jìn)制文件,所有 GMT 命令必須以 gmt 開頭 ? module 是用于完成某個特定操作的模塊的名字 ? -A 是模塊 module 提供的選項(xiàng),xx 為選項(xiàng) -A 的參數(shù) ? +b 是選項(xiàng) -A 的子選項(xiàng),xxxx 為該子選項(xiàng)的參數(shù) 一個完整的示例: gmt pscoast -R0/20/0/20 -JM6i -Ggray -Wthin -B5 -B+t"Title with spaces" -V -P > map.ps? 不以 - 開頭的參數(shù),都會被當(dāng)做文件,GMT 會嘗試去讀取 ? 子選項(xiàng)以 + 開頭,后接單個字符以及子選項(xiàng)的參數(shù) 官方文檔中子選項(xiàng)稱為 modifier,可以譯為“修飾符”,本文檔中統(tǒng)一使用“子選項(xiàng)” 子選項(xiàng)以 + 開頭,后接單個字符以及子選項(xiàng)的參數(shù) 各選項(xiàng)間以空格分隔,選項(xiàng)內(nèi)部不能有空格 ? 選項(xiàng)內(nèi)部的字符串,若存在空格,應(yīng)用單引號或雙引號括起來 STDOUT)用于顯示輸出數(shù)據(jù),標(biāo)準(zhǔn)錯誤流(STDERR)用于顯示錯誤消息 1. 將標(biāo)準(zhǔn)輸出流重定向到數(shù)據(jù)文件中 2. 將標(biāo)準(zhǔn)錯誤流重定向到日志文件中 3. 將標(biāo)準(zhǔn)輸出流通過管道傳遞給下一個命令 對于標(biāo)準(zhǔn)輸出流: ? > :將標(biāo)準(zhǔn)輸出流重定向到新文件中。若該文件已存在,則覆蓋文件中原內(nèi)容;若該文件不存在,則 創(chuàng)建該文件 ? >> :將標(biāo)準(zhǔn)輸出流追加到文件中。若文件已存在,則將標(biāo)準(zhǔn)輸出流追加到已有文件后面;若文件不 存在,則創(chuàng)建該文件 對于標(biāo)準(zhǔn)錯誤流,重定向符號是 2> 和 2>> 。這里的 2 表示標(biāo)準(zhǔn)錯誤流,大于號的含義與標(biāo)準(zhǔn)輸出 流相同。

?-K? used in the starting plot;-K -O used in processing plots;-O used in the final plot;

管道 舉例如下,假設(shè)文件 input.dat 中包含了一系列地震的經(jīng)度、緯度和震級共三列數(shù)據(jù),想要在圖上 畫很多圓表示地震的位置,圓的大小表示震級的大小。可以用類似如下命令: gmt psxy input.dat -Rxxx ... > test.ps 此時 psxy 模塊會讀取 input.dat 文件的內(nèi)容作為其輸入。 也可以使用管道: cat input.dat | gmt psxy ... > test.ps cat 命令會讀取 input.dat 的內(nèi)容并將其發(fā)送到標(biāo)準(zhǔn)輸出流,由于使用了管道,標(biāo)準(zhǔn)輸出流中的內(nèi) 容被 gmt psxy 接收作為自己的標(biāo)準(zhǔn)輸入流。 當(dāng)然還可以使用 gawk gawk '{print $1, $2, $3/10}' input.dat | gmt psxy ... > test.ps gawk 會讀取 input.dat 的內(nèi)容,并對數(shù)據(jù)做簡單處理并輸出。

?

4. 通過管道輸入 管道可以將前一個命令的標(biāo)準(zhǔn)輸出作為后一個命令的標(biāo)準(zhǔn)輸入: cat event.loc | gmt psxy -R70/140/20/60 -JM6i -B5 -Sc0.2c -Gred -P > test.ps 5. Here Documents 示例如下,兩個 EOF 之間的所有數(shù)據(jù)都會被傳遞給 GMT: gmt psxy -R70/140/20/60 -JM6i -B5 -Sc0.2c -Gred -P > map.ps << EOF 100.0 40.0 110.0 45.0 EOF

?Quoting:

想要用 psxy 繪制某數(shù)據(jù)時,需要提供數(shù)據(jù)的范圍 -R ,而 gmtinfo 模塊可以用于計算并輸出 數(shù)據(jù)的范圍,即需要將 gmtinfo 的輸出作為 psxy 的一個選項(xiàng)。 比如: $ gmt info in.dat -I1/1 -R0/10/0/10 $ gmt psxy in.dat -JX10c -R0/10/0/10 > map.ps 上面的做法需要人工干預(yù),不適合腳本自動化,可以利用倒引號將 gmtinfo 的輸出保存到變量中: #!/bin/bash R=`gmt info input -I1/1` gmt psxy in.dat -JX10c $R > map.ps 上面的例子還可以進(jìn)一步簡化。此處變量 $R 只需要用一次,因而沒有必要把 gmtinfo 的輸出信息 保存到變量中,可以直接在 psxy 命令中使用倒引號: $ gmt psxy in.dat -JX10c `gmt info in.dat -I1/1` > map.ps

?

?

EOF 可以被替換成其他任意字符(比如 END ),只要保證開始和結(jié)束的符 號一致即可

?

轉(zhuǎn)載于:https://www.cnblogs.com/seisjun/p/6201954.html

總結(jié)

以上是生活随笔為你收集整理的CMT learning的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。