GrADSのStation Dataは、作成するまで非常に手間がかかります。一方、GMTではタブ区切りor空白区切りのテキストデータを使用することができます。 さすがにGrADSでパイプ(|)入力させることは難しいですが、テキストデータから読み込むスクリプトを作成しましたのでご紹介します。 なお、このページで使用したデータは2007年4月10日00UTCにおける、QC後のSSMISの22GHz垂直偏波の輝度温度(K)などです。 |
GrADSスクリプトを使用して描いた図 |
環境変数GASCRPに設定されているディレクトリに GrADS_Scripts.tar.gz を展開する。(右クリックして保存。うまく展開できない場合)
環境変数GASCRPを設定していない場合 >>
テキストで [ x方向の位置 y方向の位置 値 ] をタブ区切りor空白区切りで1行ずつ並べます。カンマ区切り(CSV形式)のデータは扱えませんので要注意!
実際の地点データなら [ 経度 緯度 値 ] の順になります。経度は-180〜180でも0〜360でもOKです。
改行コードは「LF
」にしてください。(「CR+LF
」の場合は行末にスペースを入れると上手く動くかも?)
smcol
と一致している部分も多いので詳細は省略。異なる部分を中心に以下で例を示します。
使用したいファイル名が tbb22.txt だとします。
と入力すると、単にforground colorでプロットされます。 以下のような標準出力も出てきます。1, 2行目は色未設定時に出力されます。4行目以降にプロットデータの統計量が表示されます。 別の単色で描くには、
等とします。 |
|
3列目の値を使用して、値に応じて色が変化するように描画するには、
このとき、先ほど示した標準出力 |
|
描画領域を指定するには
としますが、
とします。(なぜかというと、コントロールファイルを開いていないので、 |
|
プロットマークの形と大きさを指定するには
と設定する代わりに、
とします。 |
|
出力時に機械間引きを行うには
|
|
描画領域の背景を塗りつぶすには
なお、陸域でマスクしているように見えるかもしれませんが、そのようなことはありません(元データで陸域を除いている)。 |
|
プロットする点の間を直線でつなぐには、
とします。 なお、分かりやすいように点を大きくしています。領域内から外、あるいは外から中へデータが連続している場合には、枠の縁まで線を引きます。ただし、領域外の点同士の場合には直線を引くことができません。 また、データはてきとーに作ったただのサンプルです。 |
|
点ではなく、文字をプロットすることもできます。その場合には入力データの4列目以降に文字を並べ、
とすれば、色は2(赤)、中心位置(0)に、やや大きめ(1.5倍)の文字をプロットできます。 |
|
以下のように色々設定すると、右の図のようになります。
|
|
2つの変数の比較など、地図を必要としない散布図を作成することもできます。まず、
として、投影法を変えてから、
とします。 |