WANtaroHP (GMT: 各種事例)

ページ概要

このページでは，以下の事項について紹介しています．

日本地図

GMTによる作図のためのスクリプトなど

世界地図

GMTによる作図のためのスクリプトなど

世界地図

大圏航路

シュミットネット (ランベルト正積方位図法)

Script and figure by GMT

東アジア震央マップ

準備

地震データはUSGSのウェブページよりダウンドードしました．

データ書式は以下のとおりです．

      Link to USGS Home Page 


  NEIC: Earthquake Search Results


                    U.  S.  G E O L O G I C A L  S U R V E Y

                     E A R T H Q U A K E  D A T A  B A S E



 FILE CREATED:  Sun Oct 30 05:53:38 2011
 Global Search   Earthquakes=    383175
 Catalog Used: PDE
 Date Range:   1900/01/01   to    2011/10/31
 Magnitude Range:   3.0  -  10.0
 Depth Range:    0  -  100
 Data Selection: Historical & Preliminary Data


 CAT    YEAR  MO DA  ORIG TIME   LAT    LONG  DEP  MAGNITUDE  IEM DTSVNWG DIST
                                                              NFO          km
                                                              TF

 PDE    1973  01 01 034609.80  -9.21  150.63  41  5.3 mbGS    ... .......
 PDE    1973  01 01 052229.80 -15.01 -173.96  33  5.0 mbGS    ... .......
 PDE    1973  01 01 114237.50 -35.51  -16.21  33  6.0 MsGS    ... .......
.......... .......... ..........
 PDE-Q  2011  10 28 084809.40  38.49   43.77  10  4.6 mbGS    ... .......
 PDE-Q  2011  10 28 090837.95 -40.63  126.42  10  4.7 mbGS    ... .......
 PDE-Q  2011  10 28 091846.27  35.53  -97.37   5  3.5 LgGS    ... .......

------------------------------------------------------------------------

*USGS National Earthquake Information Center*
USGS Privacy Statement  | Disclaimer
 	FirstGov 

GMTによる作図のためのスクリプトなど

(注) 凡例描画用入力データはGMT作図用バッチファイルの中で作成しています．

各種説明図のサンプル(説明図作成事例)

GMT による作図のためのスクリプトなど (埋設パイプモデル)

GMT による作図のためのスクリプトなど (厚肉円筒モデル)

気候比較図(気温と降水量)

概要

• 日本国内 4 都市の，月平均気温と月間雨量を示したグラフを作成しました．
• GMT の -X および -Y オプションにより，グラフの左下座標をずらしながら，4 枚のグラフを 1 ファイルに配置しています．
• GMT 実行用バッチファイルは，コントロール部と主要作図部に分割し，主要作図部を繰り返し呼び出すことにより4 都市分のグラフを作成しています．
• 4 都市の入力データは，gnuplot による気候比較図に用いているものと同一です．このため，GMT で読み込みやすいように，コントロール用バッチファイル内で，gawk によりファイルの分割・整形を行っています．

P オプションと画像の回転

• この事例では，作図用バッチファイル内の GMT コマンドでは，-P オプションを指定していないため，A4 横タイプの eps ファイルが作成されています．
• -P オプションを指定しない場合，gsview で確認する場合は問題ありませんが，png 画像にする場合，90 度回転させる必要があります．
• このため，作成された eps 画像を，ImageMagick コマンド(注*)により，90 度回転させて png 画像に変換しています．

GMT による作図のためのスクリプトなど

時系列図(気温と雨量の経時間変化)

概要

• 月平均気温と月間雨量の経時変化を示すグラフです．
• 気温は赤い実線で，雨量は水色の棒グラフの外枠をマゼンタでなぞって作成しています．
• GMTでは色名シアンは使えますが色名アクアは使えませんので注意！

GMT による作図のためのスクリプトなど

(注) 凡例描画用入力データはGMT作図用バッチファイルの中で作成しています．

時系列図(河川流量の経時間変化)

概要

• 河川流量の日平均値の経時変化を示すグラフです．
• ついでに流況図を重ねて，簡易な流況表を右上につけてみました．
• 流況表のデータは，流況図を作る際ソートしたデータより自動的に拾っています．

GMT による作図のためのスクリプトなど

(注) 凡例描画用入力データはGMT作図用バッチファイルの中で作成しています．

普通軸グラフ(地震動加速度時刻歴)

概要

• 3つの地震動(EW方向，NS方向，UD方向)の加速度時刻歴をプロットした事例です．
• 縦軸の値はそれぞれの地震動の加速度の最大値または最小値と一致するように設定しています．これは，awkスクリプトにより加速度の最大値(最小値)を検索するとともに軸設定を行うバッチファイルを作成・読み込むことで GMT に渡しています．管理人の知る限り，gnuplot ではできない技だと思います．

作図用入力データ書式

2011/03/11 14:46:00 TCGH16 EW2 Max.acc=1196.697  | Comment
dt,0.010                                         | Sampling period (sec)
ndata,30000                                      | Number of data
-1.890                                           | Acceleration (gal)
-1.888                                           | ...
-1.888                                           | ...
-1.891                                           |
.....                                            |
.....                                            |

GMT による作図のためのスクリプトなど

両対数グラフ(加速度応答スペクトル)

概要

• 3つの地震動(EW方向，NS方向，UD方向)による加速度応答スペクトルを1枚のグラフに重ねがきした事例です．

作図用入力データ書式

2011/03/11 14:46:00 TCGH16 EW2 Max.acc=1196.697             | Comment
nt,200                                                      | Number of data
period(sec)  , acc(gal)     , vel(kine)    , dis(cm)        |
2.000000e-002, 1.196764e+003, 2.924900e-001, 1.213241e-002  | Period and
2.063122e-002, 1.198466e+003, 3.113090e-001, 1.292889e-002  |    Response
2.128236e-002, 1.200244e+003, 3.314467e-001, 1.377853e-002  | .....
.....                                                       |
.....                                                       |

GMT による作図のためのスクリプトなど

(注) 凡例描画用入力データはGMT作図用バッチファイルの中で作成しています．

両対数グラフ(3重応答スペクトル)

概要

• 3つの地震動(EW方向，NS方向，UD方向)による応答スペクトルを3重応答スペクトル図として示したものです．
• このグラフでは，まず速度応答スペクトルを両対数グラフとして描画し，そこに加速度と変位の軸を重ね書きしています．

作図用入力データ書式

2011/03/11 14:46:00 TCGH16 EW2 Max.acc=1196.697             | Comment
nt,200                                                      | Number of data
period(sec)  , acc(gal)     , vel(kine)    , dis(cm)        |
2.000000e-002, 1.196764e+003, 2.924900e-001, 1.213241e-002  | Period and
2.063122e-002, 1.198466e+003, 3.113090e-001, 1.292889e-002  |    Response
2.128236e-002, 1.200244e+003, 3.314467e-001, 1.377853e-002  | .....
.....                                                       |
.....                                                       |

GMT による作図のためのスクリプトなど

(注) 凡例描画用入力データはGMT作図用バッチファイルの中で作成しています．

片対数グラフ(フーリエスペクトル比)

概要

• 2つの地震動によるフーリエスペクトルとその比を片対数グラフに表した事例です．
• スペクトル比のピークはawkスクリプトで検索し，GMTで表示させています．
• この事例では，W軸をスペクトル比，E軸をスペクトルとして描画しています．
• 出力事例1では，3本のグラフの線が重なって見にくくなっています．
• 出力事例2ではE軸の範囲をずらしています．また軸数値描画のためのコマンドラインオプション -B に p を加えることにより，べき乗表示としてみました．

作図用入力データ書式

i, frq(Hz)                 , sp1               , sp2               , spr
1, 0.0000000000000000      , 68.644469660094956, 93.489593441899203, 1.3619391905106006
2, 3.05175781250000000E-003, 68.646266058682215, 93.493209037733777, 1.3619562199903368
3, 6.10351562500000000E-003, 68.651654422441766, 93.504054788801866, 1.3620073044916456
.....
.....

GMT による作図のためのスクリプトなど

(注) 凡例描画用入力データはGMT作図用バッチファイルの中で作成しています．

ヒストグラム(最大日雨量・年間雨量)

概要

• GMTのもつ機能pshistogramでヒストグラムを作成した事例です．
• この事例では軸範囲は指定せず(-Rを指定しない)，縦軸･横軸の範囲はGMTが自動的に定めています．
• またGMTのもつ機能minmaxで入力データのデータ数・最小値・最大値を調べ，グラフ内に描画しています．

GMT による作図のためのスクリプトなど