GMT(General Mapping Tool)

2015年12月15日

GMTは，いろいろな投影法を使って地図を描くことができます．また，地図の中に，観測点や観測点の変位，あるいは地震のメカニズムなどを描くことができます．地図だけでなく，２次元のグラフや横軸が時間のグラフも描くことができます．


１．まずはエディタを使って以下の文字（GMTのコマンド）を打ち込んでみましょう．ファイル名はたとえば，ex01.shとします．

pscoast -R0/360/-70/70 -Jm4.e-2 -Ba60f30/a30f15 -Dc -W1 -Gblack -P > GMT_mercator.ps

このできたファイルを以下のように実行するとメルカトル図法で書かれた世界地図がかけます．
[nakao@kokufu2 nakao] sh ex01.sh リターン
[nakao@kokufu2 nakao] gs GMT_mercator.ps リターン
と入力すると世界地図が画面に現れるはずです．修了するときはquitリターンと入力します．

以下のように”-Gblack”を書かないと陸地は白いままです．
pscoast -R0/360/-70/70 -Jm4.e-2 -Ba60f30/a30f15 -Dc -W1 -P > GMT_mercator.ps

２．ex01.shに以下の３行を加えてhawaiiの位置に赤丸を書きます．このとき先ほどのpscoastの行に−Kを追加します．

pscoast -R0/360/-70/70 -Jm4.e-2 -Ba60f30/a30f15 -Dc -W1 -P -K > GMT_mercator.ps
psxy -R -Jm -Sc0.5 -Gred -W1/red  -O << END >> GMT_mercator.ps
-155.5 19.6
END

再度ex01.shを実行します．gsを使ってできた絵を見ます．

３．さらに以下の３行を加えてHawaiiという文字を書きます．

pscoast -R0/360/-70/70 -Jm4.e-2 -Ba60f30/a30f15 -Dc -W1 -Gblack -P -K > GMT_mercator.ps
psxy -R -Jm -Sc0.5 -Gred -W1/red -O -K << END >> GMT_mercator.ps
-155.5 19.6
END
pstext -R -Jm -Gorange -O <<END >> GMT_mercator.ps
-155.5 19.6 24 0.0 0 1 Hawaii
END

４．ハワイ島の移動方向と年率を矢印で書き込む．最後の３行を加えます．

pscoast -R0/360/-70/70 -Jm4.e-2 -Ba60f30/a30f15 -Dc -W1 -Gblack -P -K > GMT_mercator.ps
psxy -R -Jm -Sc0.5 -Gred -W1/red -O -K << END >> GMT_mercator.ps
-155.5 19.6
END
pstext -R -Jm -Gorange -O -K <<END >> GMT_mercator.ps
-155.5 19.6 24 0.0 0 1 Hawaii
END
psvelo -R -Jm -W1/blue -Gblue -O -Sr0.05/0.05 << END >> GMT_mercator.ps
-157.86 21.30 -61.4 32.4
END

５．ハワイ周辺で発生した地震のメカニズム解を描きます．以下のスクリプトのうち最後の３行を加えます．

pscoast -R0/360/-70/70 -Jm4.e-2 -Ba60f30/a30f15 -Dc -W1 -Gblack -P -K > GMT_mercator.ps
psxy -R -Jm -Sc0.5 -Gred -W1/red -O -K << END >> GMT_mercator.ps
-155.5 19.6
END
pstext -R -Jm -Gorange -O -K <<END >> GMT_mercator.ps
-155.5 19.6 24 0.0 0 1 Hawaii
END
psvelo -R -Jm -W1 -O -Sr0.05/0.05 -K << END >> GMT_mercator.ps
-157.86 21.30 -61.4 32.4
END
psmeca -R -Jm -Sa0.5 -O -W1 -C1 << END >> GMT_mercator.ps
-155.454 19.430 7.0 137.0 35.0 8.0 6.7 -145.0 0.0
END