20051114北海道の白図ができたら、等雨量線図にもチャレンジしてみましょう。以下1〜4の準備をしてください。作業フォルダを設定しよう。ここでは\GMTの下に、\gmtsvript、\topodata、\raindatフォルダを作ってください。別に名前はなんでもいいのですが、あとで分かりやすくするためこうしています。形でファイルしてください。(contour.batはここからダウンロード可能)また、観測地点毎の総雨量を面的に補完するため、surfaceコマンドで作ったグリッドデータtempr19920808-09.grd(約20MBあります)もここのフォルダに置いておきます。(hokkaido.rivはここからダウンロード可能)※hokkaido.rivの作り方はここをクリックしてください。フォルダに全道アメダスデータo920808.csvを置きます。(ここからダウンロード可能)。このデータには1992年8月8日〜9日の全道アメダス観測所の観測所コード、緯度、経度と総雨量の値が入っています。ここで用いている全道アメダスのデータには、欠測のある観測所のデータは除いています。アメダス観測所の緯度経度などの諸量は以下サイトからも入手できます。いうポストスクリプト形式の北海道海岸線+等雨量線図が出来上がっています(いるはず!)。これを見るには、イラストレータかgsviewというビューワが必要です。上の様な絵ができているのを見ることができるはずです。(結果のポストスクリプトファイルはここからダウンロードできます)タイトルや河川名などを日本語で修正する場合、GMTで作ったポストスクリプトファイルをイラストレータで開き、編集したほうが便利かもしれません。基本的な文法はWINDOWSのコマンドスクリプトです。set文は、ある変数になにかの値を代入する命令。数字でも文字列でもいいようです。set文で定義された変数は、別のコマンド中で%記号に挟まれて使用できます。そもそもUNIX系のコマンドawk(オークと読むらしい)をwindowsで使えるようにしたコマンド。ここではファイル中のデータ列の何番目と何番目を抽出する、といった機能しか使ってませんが、他にも機能があるようです。現段階ではこんないい加減な説明しかできず、心苦しい限りです。あとは変数をいろいろいじってみてどうなるか、各種サイトで公開されている用例を見ながら、実践あるのみです!上図で使用したデータは環境研究室で平成16年度に整理したデータをベースにしています。アメダスのデータであればWEB上で公開されていますので、自分でもデータを入手して等雨量線図を作ることができますよね?北海道海岸線+等雨量線図用スクリプト本スクリプトは\gmtscript\フォルダに入っているものとする。↓読み込みファイル名は8文字以内じゃないとうまく読み込んでくれないようです。of:最終的な出力イメージ(ポストスクリプト)ファイル名を格納する変数。このpsファイルにイメージが重ねられていく。0mmから250?まで50?ピッチで色を割り当てている。-zで連続的な色の変化を指定。↑wakuは描画の範囲(緯度経度)を示す。↑a:軸目盛数字の間隔(度)、b:白黒帯の間隔(度)、g:線の間隔(度)。r:1度あたり何cmで書くかという意味。↓基本の地図の枠を作る。この枠のイメージをベースに図が追加されていく。−Oは追加のとき必須。↓surfaceコマンドで面的に補間し、グリッドデータ化。このグリッドデータが等雨量線図に必要。グリッドデータを作ってしまえば、以下のコマンド2行をREM文にして実行しない。等雨量線図を線でなく面的に描く場合はこのコマンドを実行。↑ここでは%raindat%の3列目$3と2列目のデータ$2をパイプ"|"コマンドでpsxyに渡している。3列目と2列目かというと、このデータは3列目が経度(X軸)、2列目が緯度(Y軸)であるから。↓全道河川(%hokriv%)を描く。Surfaceコマンドでグリッド化した雨データを使っている。Sコマンドで指定しているのは水域のベタ塗りの色で、ここでは白(255)を指定。これを指定しなければ海の上にも等雨量線図がうねうね出てしまい、シマラない絵になってしまう。%of%REM以上、ランダムに配置されている雨量観測所のデータをグリッドデータにするあたりが面倒ですが、それさえクリアすれば、それほど面倒ではないと思います。各コマンドで定義している変数の意味はほかにもっと丁寧な説明がされているサイトがありますので、調べて見てください。するともっとエレガントな使い方があるかも。アーカイブフォトアルバムリンク集サイトマップお問い合せ。