PyGMTの基本

ここでは，簡単な地図の作成を通じてPyGMTの使い方の基本を学びます．

準備¶

Jupyter Notebookの立ち上げ¶

本節以降では，Jupyter NotebookでインタラクティブにPyGMTを実行していくことにします．そこで，（minicondaによるインストールの場合は仮想環境をactivateしたのちに）端末から

jupyter notebook

によりJupyter Notebookを起動します．すると自動的に既定のブラウザが立ち上がります．別のブラウザウィンドウからも，http://localhost:8888 にアクセスするとJupyterにアクセスができるはずです．この記事もそうしてJupyter Notebook上で作成しています．

モジュールのimport¶

PyGMTを使うためには，まず import します．

import pygmt

これで，pygmt.XXX という形でPyGMT内のオブジェクトを呼び出せるようになりました．

また，画像の保存場所として fig ディレクトリを作成しておきましょう．Jupyter Notebookでは行の先頭に ! をつけるか，あるいはセルの1行目に %%bash と書くことでシェルコマンドを実行できます．

%%bash
# fig ディレクトリが存在しなければ作成する
if [ ! -d fig ]; then
  mkdir fig
fi

PyGMTは活発な更新がなされているため，このページに含まれているサンプルコードも，バージョンによっては修正が必要であったり，警告がでたりします．自分の使っているPyGMTのバージョンを確認するには，以下のようにします．

print(pygmt.__version__)

そのほか，pygmt.show_versions() という関数では，関連ライブラリのバージョンやその格納状況まで詳しい情報が得られます．

PyGMTによる作図の基本¶

pygmtでは，まず

fig = pygmt.Figure()

で図のオブジェクトを作成し，GMTのコマンドに相当する

fig.XXXXX()

を重ねていくことで作図していきます．最後に

fig.show()

とすればJupyter Notebook上に表示されますし，

fig.savefig(filename)

とするとfilenameの拡張子に応じて画像ファイルとして保存できます．もうすこし融通の効く（その代わりにやや複雑な）画像保存の方法はこちらで説明します．

基本の地図¶

まずは単純な，従来のGMTでの gmt pscoast コマンドに相当する地図を作図してみましょう．各オプションの意味はコード内コメントを参照してください．

fig = pygmt.Figure()

fig.coast(
    projection  = 'N138/12c',                         # 投影法とサイズ．メルカトルなら 'M12c' など．
    region      = (128, 146, 30, 46),                 # 領域は西端・東端・南端・北端の順にリストかタプルで与える．
    shorelines  = 'default,black',                    # 海岸線のペンの設定（後述）
    area_thresh = 100,                                # 細かい領域の描画下限 (km^2) 値を大きくすると小さな島が消える
    resolution  = 'f',                                # 'c', 'l', 'i', 'h', 'f' の順に高くなる
    land        = '249/243/202',                      # 陸地の色を 0-255の値の R/G/B で指定
    water       = '190/220/250',                      # 水の色を 0-255の値の R/G/B で指定
    map_scale   = '142.5/32/32/400',                  # 距離スケールの水平・縦位置と距離サイズ
    frame       = ['WSen+tFigure', 'xafg', 'yafg']    # 南西に軸名・タイトル指定（+t) およびXY方向それぞれの軸情報
)

fig.show()
fig.savefig('fig/pygmt-coastline.png')

Jupyter Notebookにおける画像の表示¶

結果をJupyter Notebook内で見るだけならば fig.show()を指定します．画像として保存するなら fig.savefig ですが，pngのほか pdfなども指定可能です．また，show=Trueオプションをつけると，その環境で標準のビューワ（たとえばmacOSならPreview.app）が起動します．

また，fig.show() に width=オプションで数字を指定すると，表示される図の大きさを指定できます．

投影法¶

fig.coastのprojectionオプションで地図投影法を指定しています．ここではロビソン図法（-JN）を用いています．この図法では，中心線の経度（138）と図の横幅の大きさ（12c）を指定する必要があります．投影法の特徴により，図の縦の大きさは描画範囲と横幅から自動的に決まります．

よく使われるメルカトル図法の場合は中心の経度は必要なく，たとえば -JM12c のように指定します．この場合も縦方向の大きさは指定できません．

GMTには多種多様な投影法が用意されていますが，ここではその網羅的な説明はしません．必要に応じて他解説や公式マニュアルを参照してください．

ペンと色¶

GMTにおける線の描画は

太さ (太さの名前あるいは数値）
色（色名あるいは RGB）
パターン（パターン名あるいは数値）

の3つの情報を , で組み合わせます．上記のスクリプトの場合は太さが default，色が black を組み合わせています．パターンは実線・破線・点線などですが，省略したことでデフォルトの実線が選ばれています．

太さの選択肢には faint, default, thinnest, thinner, thin, thick, thicker, thickest, fat, fatter, fattest があります．左ほど細いです．色名は主な名前は使えるようです．パターンはdotted, dashed, solidのほか，たとえば..だと点線，--だと破線，-.-だと一点鎖線のような指定もできます．

色は名前のほか，光の三原色 Red(R), Green(G), Blue(B)の大きさを0-255までの数値で表して，R/G/B で表現することもできます．

地図描画 Tips: `resolution` と `area_thresh`¶

海岸線の解像度は5種類ありますが，粗すぎても細かすぎても綺麗な図になりません．大は小を兼ねると思うかもしれませんが，細かすぎると微細な地形の特徴まで描画されるため，細かな島やリアス式海岸などでごちゃごちゃした印象になりがちです．解像度だけを変えた図を以下に示します．

fig = pygmt.Figure()

fig.shift_origin(3, 25)

fig.coast(
    projection  = 'N138/8c',                          
    region      = (128, 146, 30, 46),                 
    shorelines  = 'default,black',                    
    resolution  = 'f',                                
    land        = '249/243/202',                      
    water       = '190/220/250',                      
    map_scale   = '142.5/32/32/400',                  
    frame       = ['WSen+tresolution = f', 'xafg', 'yafg'] 
)

fig.shift_origin(10, 0)
fig.coast(
    projection  = 'N138/8c',                          
    region      = (128, 146, 30, 46),                 
    shorelines  = 'default,black',                    
    resolution  = 'h',                                
    land        = '249/243/202',                      
    water       = '190/220/250',                      
    map_scale   = '142.5/32/32/400',                  
    frame       = ['WSen+tresolution = h', 'xafg', 'yafg']
)

fig.shift_origin(-10, -11)
fig.coast(
    projection  = 'N138/8c',                          
    region      = (128, 146, 30, 46),                 
    shorelines  = 'default,black',                    
    resolution  = 'i',                                
    land        = '249/243/202',                      
    water       = '190/220/250',                      
    map_scale   = '142.5/32/32/400',                  
    frame       = ['WSen+tresolution = i', 'xafg', 'yafg']
)

fig.shift_origin(10, 0)
fig.coast(
    projection  = 'N138/8c',                          
    region      = (128, 146, 30, 46),                 
    shorelines  = 'default,black',                    
    resolution  = 'l',                                
    land        = '249/243/202',                      
    water       = '190/220/250',                      
    map_scale   = '142.5/32/32/400',                  
    frame       = ['WSen+tresolution = l', 'xafg', 'yafg']
)

fig.show()

高解像度ほどよいかというと，特に瀬戸内海などは小さな島がたくさん描画されて，真っ黒になってしまっています．そこで活躍するのが area_thresh オプションです．このオプションを指定すると，数値で与えられた面積（km $^2$ ）以下のサイズのものが描画されなくなります．そこで，海岸線の解像度は高めにしつつ，area_threshを適宜設定すると綺麗な海岸線が描画できます．いくつかarea_threshの数値を変えた例を示しましょう．

fig = pygmt.Figure()

fig.shift_origin(3, 25)

fig.coast(
    projection  = 'N138/8c',                         
    region      = (128, 146, 30, 46),                 
    shorelines  = 'default,black',                    
    resolution  = 'f',  
    land        = '249/243/202',                      
    water       = '190/220/250',                      
    map_scale   = '142.5/32/32/400',                  
    area_thresh = 0, 
    frame       = ['WSen+tarea_thresh=0', 'xafg', 'yafg']
)

fig.shift_origin(10, 0)
fig.coast(
    projection  = 'N138/8c',                          # 投影法とサイズ．メルカトルなら 'M12c' など．
    region      = (128, 146, 30, 46),                 # 領域はリストかタプルで与える．
    shorelines  = 'default,black',                    # 海岸線のペンの設定．
    resolution  = 'h',                                # 'c', 'l', 'i', 'h', 'f' の順に高くなる
    land        = '249/243/202',                      # 陸地の色を 0-255の値の R/G/B で指定
    water       = '190/220/250',                      # 水の色を 0-255の値の R/G/B で指定
    map_scale   = '142.5/32/32/400',                  # 距離スケールの水平・縦位置と距離サイズ
    area_thresh = 100, 
    frame       = ['WSen+tarea_thresh=100', 'xafg', 'yafg']
)

fig.shift_origin(-10, -11)
fig.coast(
    projection  = 'N138/8c',                          # 投影法とサイズ．メルカトルなら 'M12c' など．
    region      = (128, 146, 30, 46),                 # 領域はリストかタプルで与える．
    shorelines  = 'default,black',                    # 海岸線のペンの設定．
    resolution  = 'i',                                # 'c', 'l', 'i', 'h', 'f' の順に高くなる
    land        = '249/243/202',                      # 陸地の色を 0-255の値の R/G/B で指定
    water       = '190/220/250',                      # 水の色を 0-255の値の R/G/B で指定
    map_scale   = '142.5/32/32/400',                  # 距離スケールの水平・縦位置と距離サイズ
    area_thresh = 10000, 
    frame       = ['WSen+tarea_thresh=10000', 'xafg', 'yafg']

)

fig.shift_origin(10, 0)
fig.coast(
    projection  = 'N138/8c',                          # 投影法とサイズ．メルカトルなら 'M12c' など．
    region      = (128, 146, 30, 46),                 # 領域はリストかタプルで与える．
    shorelines  = 'default,black',                    # 海岸線のペンの設定．
    resolution  = 'l',                                # 'c', 'l', 'i', 'h', 'f' の順に高くなる
    land        = '249/243/202',                      # 陸地の色を 0-255の値の R/G/B で指定
    water       = '190/220/250',                      # 水の色を 0-255の値の R/G/B で指定
    map_scale   = '142.5/32/32/400',                  # 距離スケールの水平・縦位置と距離サイズ
    area_thresh = 100000, 
    frame       = ['WSen+tarea_thresh=100000', 'xafg', 'yafg']

)

fig.show()

後半は明らかにやりすぎですね．日本列島全体を表示するなら，area_thresh=100程度がほどよいように思います．もっと細かな地域を表示する場合は，適宜area_threshの値を小さくするとよいでしょう．

`basemap` とパラメタの省略¶

上記の例では fig.coast ひとつですべてを描画しましたが，軸の描画を fig.basemap に任せることもできます：

fig = pygmt.Figure()


fig.coast(
    projection  = 'N138/12c',          
    region      = (128, 146, 30, 46),  
    shorelines  = 'default,black',     
    area_thresh = 100,                 
    resolution  = 'f',                 
    land        = '249/243/202', 
    water       = '190/220/250', 
    map_scale   = '142.5/32/32/400',   
)

fig.basemap(
    frame       = ['WSen+tFigure', 'xafg', 'yafg'] 
)

fig.show()

上記の例では，fig.basemap のオプションは frame しかないことに注目してください．

必須と思える projection や region といったオプションは，その前の fig.coast で用いられた値が自動的に引き継がれています．

エイリアス¶

PyGMTでは，従来のGMTと比べてオプションの意味が多くの人にとってわかりやすい名前がつけられてる一方，従来オプション名のエイリアスも設定されています．たとえば先程の地図は，以下のようなエイリアスオプション名を用いて作成することもできます．初めて見る人には全く意味不明のオプションでしょうが，GMTに馴染みのある人にはむしろ使いやすいかもしれません．

fig = pygmt.Figure()

fig.coast(
    J = 'N138/12c',                        
    R = '128/146/30/46',                
    W = 'default,black',                   
    A = 100,                               
    D = 'f',                               
    G = '249/243/202',                     
    S = '190/220/250',                     
    L = '142.5/32/32/400',                 
    B = ['WSen+tFigure', 'xafg', 'yafg'] 
)

オプション名のエイリアスなので当然ですが，結果として得られる画像は全く同じです．ただし，

Short-form parameter (A) is not recommended. 
Use long-form parameter 'area_thresh' instead.

のような警告がたくさん出てきます．

従来のGMTユーザーは，この警告を逆手にとって，まず自分の知っている short-form オプションでpygmtを動かしてみて，そこで表示される SyntaxWarning を見て現代的な long-form parameter 名を調べるようにすると，ドキュメントを引く手間が省けてよいかもしれません．

さらに進んだ応用としては，公式には実装，ドキュメント化されていない機能の利用 があります．PyGMTの公式マニュアルを読んでいると，GMTの各モジュールに対して，必ずしもすべてのオプションが実装，ドキュメント化されているわけではないことがありますしかし，そのような場合でもshort-formオプションをつかうと「動いてしまう」ことがしばしばあるよう．もちろんこれは自己責任での利用になりますし，今後PyGMTの開発が進むことで減っていくことでしょう．

従来型のGMTスクリプトとの比較¶

これまでに紹介した PyGMTのスクリプトと，それらと（ほとんど）同じ図を出力するbashによる従来型のGMTスクリプトを比較してみます．以下のタブをクリックするとコードが表示されます．なお，比較しやすいようできるだけ改行を多めに書いています．

pygmt original

pygmt alias

bash classic

bash modern

fig = pygmt.Figure()

fig.coast(
    projection  = 'N138/12c',
    region      = (128, 146, 30, 46),
    shorelines  = 'default,black',
    area_thresh = 100,
    resolution  = 'f',
    land        = '230/240/220',
    water       = '220/235/250',
    map_scale   = '142.5/32/32/400',
    frame       = ['WSen+tFigure', 'xafg', 'yafg']
)

fig = pygmt.Figure()

fig.coast(
    J = 'N138/12c',                        
    R = (128, 146, 30, 46),                
    W = 'default,black',                   
    A = 100,                               
    D = 'f',                               
    G = '230/240/220',                     
    S = '220/235/250',                     
    L = '142.5/32/32/400',                 
    B = ['WSen+tFigure', 'xafg', 'yafg'] 
)

#!/bin/bash

# 海岸線と枠を描画
gmt pscoast -JN138/12c \
            -R128/146/30/46 \
            -Wdefault,black \
            -A100 \
            -Df \
            -G230/240/220 \
            -S220/235/250 \
            -L142.5/32/32/400 \
            -BWSen+t"Figure" -Bxafg -Byafg \
> fig/classic-coastline.ps

# Postscript画像をPNGに変換
gmt psconvert -Tg -A -P -Qg2 -Qt2 fig/classic-coastline.ps

# 不要な中間生成物であるPostscriptを削除
rm -f fig/classic-coastline.ps

#!/bin/bash

# modern mode では gmt begin ~ end の間でコマンド実行．
gmt begin fig/modern-coastline png

    gmt coast \
        -JN138/12c \
        -R128/146/30/46 \
        -Wdefault,black \
        -A100 \
        -Df \
        -G230/240/220 \
        -S220/235/250 \
        -L142.5/32/32/400 \
        -BWSen+t"Figure" -Bxafg -Byafg 
        
gmt end

このような単一のコマンドからなるスクリプトではあまり差が大きくありませんが，コードの可読性は一番左のpygmt original がもっとも高いでしょう．また，PyGMTとGMT6 modern modeでは複数のGMTコマンドを用いる図においてミスを犯しやすい -O と -K オプションが不要であったり，現代では直接用いることはあまりない（と思われる）Postscriptファイルから他形式に変換する手間が省けています．

準備¶

Jupyter Notebookの立ち上げ¶

モジュールのimport¶

PyGMTによる作図の基本¶

基本の地図¶

Jupyter Notebookにおける画像の表示¶

投影法¶

ペンと色¶

地図描画 Tips: resolution と area_thresh¶

basemap とパラメタの省略¶

エイリアス¶

従来型のGMTスクリプトとの比較¶

地図描画 Tips: `resolution` と `area_thresh`¶

`basemap` とパラメタの省略¶