スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

gmt覚え書き

gmtの覚え書きです
国土地理院の数値データをGMTで描けるように
してみました。

国土地理院n数値標高データからぴっくあっぷしました。

以下描いたGMTのコマンドデータ


surface out.txt -Gout.grd -R133.625/133.75/34.5/34.583333333 -I0.0001 -V
makecpt -Cglobe -T-300/300/10 -V -Z > color.cpt

grdgradient out.grd -Gout.gradient -A45 -Ne0.4

R133.625/133.75/34.5/34.583333333 -P -K >a.ps

pscoast -JM5 -R133.625/133.75/34.5/34.583333333 -Ba0.05/a0.01 -Df -Gc -K -Xc -Yc -P >b.ps
grdimage out.grd -JM5 -Iout.gradient -Ccolor.cpt -R133.625/133.75/34.5/34.583333333 -K -O -P >> b.ps
pscoast -Q -O -K -P >>b.ps
pscoast -J -R -Df -Wthick/black -O -P >>b.ps
スポンサーサイト

コメントの投稿

非公開コメント

プロフィール

かみさまみならい

Author:かみさまみならい
FC2ブログへようこそ!

最近の記事
最近のコメント
最近のトラックバック
月別アーカイブ
カテゴリー
ブロとも申請フォーム

この人とブロともになる

ブログ内検索
RSSフィード
リンク
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。