[GMT] PSCOUPE 迷惑的 -Aa…/Dmin/Dmax[f]

今天去中研院當本年度地物暑期生的GMT小助教,雖然大多數時間是 De 同學們的 bug ,但是也有學到新東西,今天分享的是 PSCOUPE 指令中 -Aa 剖面最後兩個參數:「Dmin/Dmax[f]」的差異。

前言:

PSCOUPE 是可以幫我們化很多不同種類剖面的好用指令,舉凡地震震央深度分布、震源機制解的側投影都可以幫忙完成,非常之強大。選擇剖面種類主要是用「-A」的參數來指定,今天我們用 -Aa 來繪製地震深度的分佈圖,主要切了一個臺灣東北部的剖面如下圖:

2017 GMT 中研院暑期生 ex08

2017 GMT 中研院暑期生 ex08

官方說明:

-Aalon1/lat1/lon2/lat2/dip/p_width/dmin/dmax[f]

lon and lat are the longitude and latitude of points 1 and 2 limiting the length of the cross-section. dip is the dip of the plane on which the cross-section is made. p_width is the width of the cross-section on each side of a vertical plane or above and under an oblique plane. dmin and dmax are the distances min and max from horizontal plane, along steepest descent direction. Add f to get the frame from the cross-section parameters.

指令

主要的變數說明如下:

  • $cpt :使用的色盤。
  • $beglo: 起點的經度。
  • $begla:起點的緯度。
  • $endlo:終點的經度。
  • $endla:終點的緯度。

這邊我們 -Aa 切了一個垂直地表的剖面(dip=90),並設定寬度範圍為30公里(p_width),最後我就來嘗試四種不同的結果:

  1. (左上)原圖:-R0/200/0/250  -Aa$beglo/$begla/$endlo/$endla/90/30/0/190
  2. (右上)更改Dmin/Dmax範圍嘗試指定50-100公里深的點:-R0/200/0/250  -Aa$beglo/$begla/$endlo/$endla/90/30/100/150
  3. (右下)在 Dmin/Dmax 加上 f 的選項:-R0/200/0/250  -Aa$beglo/$begla/$endlo/$endla/90/30/100/150f
  4. (左下)改變 -R的範圍與右下圖比較:-R0/200/100/150  -Aa$beglo/$begla/$endlo/$endla/90/30/0/190

實驗結果

2017 GMT 中研院暑期生 ex08
  • 從 1、2  圖比較我們可以發現,在沒有指定 f 的情形下,Dmin/Dmax 對於剖面的內容是沒有影響的。
  • 2、3 圖比較的話,加上 f 後剖面僅繪製在「Dmin/Dmax」深度範圍內的資料,無視 -R 選項中設定的 ymin/ymax (倒數兩個選項),並且將資料延伸到 -JX 設定的繪圖範圍內。這邊可以看到有點奇妙的是在 X 大於 200 公里的位置還有點,或許在程式調整投影深度時,會影響到投影範圍設定,使有些資料點落在 200 公里以上。
  • 比較 1、3、4:4 圖改變了 -JX 的Y寬度(比例隨 -R 範圍調整保持與圖 1 相同),非常像直接切下圖一 100 至 150 公里範圍內部分,也沒有像圖 3 位於在 200 公里外有奇妙的資料點殘留。

結論

  • 在沒有設定「f」的狀況下,「Dmin/Dmax」對於剖面繪圖設定是沒有影響的。
  • 指定「f」以後,圖形會變形,繪圖範圍被調整為-JX設定的大小。
  • 指定「f」以後,有些原本並未出現的點會投影出來。

如果有對指令比較熟的朋友也請不吝指教!


請多多指教!

這個網站採用 Akismet 服務減少垃圾留言。進一步瞭解 Akismet 如何處理網站訪客的留言資料