使用データ:世界の陸地と海洋地形1分(解像度約1km)グリッドデータは下記URLからダウンロード可能です。全データは450MBくらいあるので、緯度経度で利用範囲を指定し、部分的にダウンロードすることも可能です。set文で出力するポストスクリプトファイル名の変数psfileを設定。このポストスクリプトファイルにイメージが追加されていく。set文で描く範囲を示す変数wakuを設定。左端/右端/下端/上端の順。set文で緯度経度枠のデザインを指定。a:何度毎に数値を入れるか。f:何度刻みで枠に目盛を入れるか。g:何度刻みで格子を描くか(ここでは0なので図内に格子が入らない)。psbasemap:画像の下地の設定。-JbのbはAlbers図法で投影されることを示します。grdimage:grddata変数で定義されたグリッドデータを指定された情報で投影します。pscoast:指定された情報で海岸線を投影します。変数名の前後に%を付けるのはWindowsでシェルスクリプトを動かすときの約束事らしいのですが、詳しいことはよく知りません。grddata=..\Gridone.grdとなります。以上のスクリプトを実行すると、Output_image.psというファイルが出来上がります。このファイルのイメージを見るには、ポストスクリプトファイルを閲覧可能なGsviewやAdobe社製イラストレーターなどのソフトが必要です。世界地図を描く場合、上記の1km解像度のデータだと処理が長引く恐れがあるので、5km解像度くらいのものが適当だと思います。それでも、下記スクリプトで得られる結果は350MBくらいのサイズになります。5km解像度のデータはNOAAの以下のサイトからDLできます。世界の地形と海洋のデータはこのサイト中、カバレージがGLOBAL_LAND_and_OCEANSのもので、Both(両方)となっているリンクところをクリックするとDLできます。%ofn%以上のスクリプトを実行すると、glb_ekr.epsというポストスクリプトファイルができます。PCのスペックにもよりますが、gsviewでイメージを開くのに数分〜10分程度かかります。jpgなどの閲覧しやすい画像に変換しておくことで、その後の処理が楽です。都市名をプロットするスクリプトは上記のスクリプトに以下の行を付加し、下記の都市データを同一フォルダ内に置けばOKです。コマンドの詳しい説明は省略します。gawk文のあたりはちょっと複雑そうに見えますが、パイプコマンドを示す"|"のところで左と右に分かれています。pstext用のデータをパイプで右に受け流しているのです。左からなにかがやってきて、右にそのまま受け流すワケですね。!?どこかで聞いたような、、、。※主要都市名をプロットする場合のデータ形式サンプル。経度)で都市のデータがある場合、地図上に都市名をプロットできます。し、以下の並びでデータが入っているものとします。アーカイブフォトアルバムリンク集サイトマップお問い合せ。