ＧＭＴたまごくらぶ

ＧＭＴ(The Generic Mapping Tools)を使っていろんな図を書きましょう。
ＧＭＴは地図を描く能力が高いので一気にメジャーになりましたが、それ以外の 一般的なＸＹプロットとかを描くのにも便利です。UNIXでシェルスクリプトとかを 書けるようになるともっと便利に使えます。

ところが、ドキュメントが英語であること、ドキュメントに含まれる例がかなり 高度であること、かなり柔軟なプログラム群であること、などの理由で最初に 「ＧＭＴ文化」に慣れるまで、少し時間がかかります。

「たまごくらぶ」の発想とは、簡単なスクリプトを紹介しつつ、はじめの一歩を お手伝いしようというものです。

似たようなことを考えた人はいっぱいいて、ＧＭＴ関連のリンクもいろいろと あります。

• 京大防災研堀川晴央氏のＧＭＴよいよいガイド 。ver3.3用もあります。
  マニュアル全部読むのもめんどくさい人向け。

  神戸大筧氏ののＧＭＴガイド

  ここよりもはるかに親切。

  SOEST U Hawaii （ＧＭＴのご本尊）

  ＧＭＴで何ができるのか？、

  SOEST U HawaiiでのＧＭＴの授業（！）

  （これさえ読めばたまごくらぶは要らないか？）

1. ＸＹプロットを描く１（psbasemap,psxy; awk, shell variables）
   
2. 地図を描く１（pscoast）
3. ＸＹプロットを描く２（gmtset,pstext）
4. 地図を描く２＋コンターマップを描く１（grdcontour,grdimage）
5. ＧＭＴ最前線：ピンクの縞々が描きたい

クックブックをちゃんと読む（日本語版を作る？）

２次元データのプロット：psxyを用いてiasp91の Ｓ波速度をプロットしよう。 （完成品と、簡単な解説） 。

もう少し線を太くするには、どうすればいいか。
線はなく、記号でプロットするのには？
Ｙ軸の数値が、上へいくほど大きくなるようにするには ？

複数のデータをプロットしよう

psxyを３回使って、iasp91モデル全部 をプロットしよう （完成品と、簡単な解説）。

-K と -Oに注意！（GMTのポストスクリプトの構造とは）
awkは便利

一度の手間で大きな効果を上げる

カラーと白黒の図を一つのスクリプトで描きわけたりすることができるのが、 GMTのスクリプトの良いところです。

シェル変数を使おう
色や線種をまとめて指定（結果はおんなじ）、簡単な 解説）。
もうちょっと、複雑なケース （同じもの）。

エラーバーを書く

結構、ある状況なので、例を一つ

カラーの図での、エラーバーの指定
データは表面波の位相速度です （結果）。

簡単な地図の書き方。簡単なメルカトル図法

pscoastの大事なオプションは-D (high, intermediate, low, crude)
川は -I、国境は -N。縮尺を書くときは-Lなのだ。

図中に文字を書こう
pstextを使う。フォントの種類。
地図中に字を書こう。簡単な場合。
pstextの入力のjustifyって指定した座標と書く字の 相対位置を決めるんです。
困ったときのpstext -L
.gmtcommandsの役割。
gmtsetと.gmtdefaults
タイトルの大きさを変えたい／緯度経度の表記を 変える
センチメートルの場合。
スクリプトはまとめてここ

おまけ：

色に溺れる
あなたのカラープリンターの性能をテストする
白黒プリンター用もできました 【もっと賢い方法を考え付いた人は教えてね】