GMT6.5のサンプルスクリプトを試す1 basemap~linear~

GMT（Generic mapping tools）6.5の公式HPに記載されている例・サンプルスクリプトが示されているが、図が示されていないので描いてみる。

アルファベット順にbasemapから順に試してみよう。

引用元
GMT6.5.0 Documentation basemap

最初は地図ではないグラフなどを描くためのもののようだ。

gmt begin linear
  gmt basemap -R0/9/0/5 -Jx1c -Bf1a2 -Bx+lDistance -By+l"No of samples" -BWeSn
gmt end show

データが入っていないため枠だけが形成された。
ごくありふれた線形のグラフ枠。
-Rで範囲を指定。X軸：0~9、Y軸：0~5
-Jで表示方法を指定。x1c：線形スケール1cm単位。小文字xではプロットするスケールと同じになる。
-Bで境界、軸を指定。複数個入力してもよいらしい。
　　-Bf1a2　f1：細かい目盛りの間隔(１ずつ)
　　　　　　a2：注釈の間隔（２ずつ）
　　-Bx+lDistance　x+l：X軸側のラベルを設定（“Distance”という文字を入れる）
　　-By+l”No of samples”　 y+l：Y軸のラベルを設定（“No of samples”という文字を入れる）””で囲うとスペースを含めたラベル、タイトル名を入れられる。スペースがないのであれば””なしでもよい模様。
　　-BWeSn　WeSn：軸・ラベルを記載する場所を示す
　　　　　　　　　　West左,East右,South下,North上の頭文字をそれぞれ大文字で示すと、軸とラベルを表示、小文字にすると軸のみ表示、記載しないと無くなる。

枠だけだとつまらないので、適当に数字を入力してプロットしてみる。
せっかくなので、適当にデータを用意してプロットしてみることにした。
sample1.txtで$y=\frac{1}{2} x^2$を以下のように準備。

0 0
1 0.5
2 1

plotはとりあえず適当に赤丸をプロットするようにした。
また、beginで出力ファイル形式を指定できるが、指定しないとPDFになり図として扱いにくいのでpngとした。

gmt begin linear2 png
  gmt basemap -R0/9/0/5 -Jx1c -Bf1a2 -Bx+lDistance -By+l"No of samples" -BWeSn
  gmt plot sample1.txt -Sc0.2 -Gred -W0.5
gmt end