Skip to main content

GMT 地球を描く

Posted in

とりあえず、丸い地球を描く方法です。

正射方位図法

無限遠に光源を置き、平面に投影します。遥か彼方から眺めた地球ですね。指定方法は以下の通りです。

 -Jg中心経度/中心緯度/[範囲/]縮尺
 -JG中心経度/中心緯度/[範囲/]幅

日本を中心に描いてみると以下のようになります。1行目はスケールを指定して、2行目は幅を指定して描いていますが、どちらでも同じような図を描くはずです。地球を1/200000000にすれば6.4cm位になりますから。

pscoast -Jg135/35/1:2e+8 -Rg -Bg30 -Dc -A1000 -G102/255/102 -P > globe1.eps
pscoast -JG135/35/6.4 -Rg -Bg30 -Dc -A1000 -G102/255/102 -P > globe2.eps

globe1.png

-Jgあるいは-JGの2番目のパラメータとして、描く範囲を指定できます。これは中心からの角度(度)で指定します。規定値の90度だと半球すべてですね。

1行目だとスケール指定なので、上の図と同じスケールとなりますが、2行目だと幅指定なので、切り出した部分が幅に合わせられます。

pscoast -Jg135/35/45/1:2e+8 -Rg -Bg30 -Dc -A1000 -G102/255/102 -P > globe3.eps
pscoast -JG135/35/45/6.4 -Rg -Bg30 -Dc -A1000 -G102/255/102 -P > globe4.eps

globe3.png