AX

ホーム>資料>パソコン

パソコンのプチ情報

前の記事   次の記事

b15.GIF GMT講座 その9.5:自作データでコンター図(2)

今度はデータ値そのものだけを用いてブロック状のコンター図を描いてみましょう。

1)まずデータを用意します。
データその1:観測値データ(X, Z, 値)
例:ファイル名:ex2_1.dat
----
-20 3 57.8158648
-20 6 67.53058295
-10 3 57.8158648
-10 6 58.78505576
-10 9 68.35012886
-10 13 89.77596304
0 3 56.54866776
0 6 57.84863714
0 9 66.87494623
0 13 89.63269659
0 18 124.3051796
10 3 57.1769863
10 6 56.87366011
10 9 68.84185641
10 13 87.51579535
20 3 60.35025887
20 6 75.06752972

データその2:多角形データ
データ値の前に「> -Z」をつけて、その下にデータ取得箇所を取り巻くように多角形の頂点の座標をいれます。
例:ファイル名:ex2_2.dat
----
> -Z  57.8158648
-25  1.5 
-25  4.5 
-15  4.5 
-15  1.5 
-25  1.5 
> -Z  67.53058295
-25  4.5 
-25  7.5 
-15  7.5 
-15  4.5 
-25  4.5 
> -Z  57.8158648
-15  1.5 
-15  4.5 
-5  4.5 
-5  1.5 
-15  1.5 
> -Z  58.78505576
-15  4.5 
-15  7.5 
-5  7.5 
-5  4.5 
-15  4.5 
> -Z  68.35012886
-15  7.5 
-15  11 
-5  11 
-5  7.5 
-15  7.5 
> -Z  89.77596304
-15  11 
-15  16 
-5  16 
-5  11 
-15  11 
> -Z  56.54866776
-5  1.5 
-5  4.5 
5  4.5 
5  1.5 
-5  1.5 
> -Z  57.84863714
-5  4.5 
-5  7.5 
5  7.5 
5  4.5 
-5  4.5 
> -Z  66.87494623
-5  7.5 
-5  11 
5  11 
5  7.5 
-5  7.5 
> -Z  89.63269659
-5  11 
-5  16 
5  16 
5  11 
-5  11 
> -Z  124.3051796
-5  16 
-5  22 
5  22 
5  16 
-5  16 
> -Z  57.1769863
5  1.5 
5  4.5 
15  4.5 
15  1.5 
5  1.5 
> -Z  56.87366011
5  4.5 
5  7.5 
15  7.5 
15  4.5 
5  4.5 
> -Z  68.84185641
5  7.5 
5  11 
15  11 
15  7.5 
5  7.5 
> -Z  87.51579535
5  11 
5  16 
15  16 
15  11 
5  11 
> -Z  60.35025887
15  1.5 
15  4.5 
25  4.5 
25  1.5 
15  1.5 
> -Z  75.06752972
15  4.5 
15  7.5 
25  7.5 
25  4.5 
15  4.5 

2)これをコンター図にするGMTスクリプトを下記に示します。
(スクリプト名:contour2.gmt)
------------------------------------------------------------
#!/bin/tcsh -f
# GMT script for making bathymetry map (Japan)

gmtdefaults -D >! .gmtdefaults
gmtset LABEL_FONT_SIZE 12 HEADER_FONT_SIZE 18

makecpt -Crainbow -Z -I -T40/150/5 >! tmp.cpt
psxy ex2_2.dat -JX6.6/-3 -R-22/22/0/20 -Ctmp.cpt \
-M -Ba5/a5WsNe:."ex2": -Xa1.0i -Ya2.0i -P -K >! tmp.ps

psscale -Ctmp.cpt -D3.5/1.3/3.0/0.2h \
-B50/:"App.Resistivity(Ohm-m)": -P -K -O >> tmp.ps

# Plot observed location
awk '{print $1, $2}' < ex2_1.dat | \
psxy -JX -R -Sc0.1i -G0 -Xa1.0i -Ya2.0i -P -K -O >> tmp.ps

echo 0 0 18 0 1 1 VES_Survey | \
pstext -JX4.5/1 -R0/16/0/1 -Xa3.5i -Ya6.0i -P -O >> tmp.ps

gv tmp.ps
mv tmp.ps ex2.ps
rm tmp*.dat tmp*.grd tmp*.cpt *grd.shade
rm .gmt*
------------------------------------------------------------

psxy の-Mオプションで多角形をプロットできるようになります。多角形内の色もファイルから読み取っています。

[関連画像]

 
図1

前の記事

目次に戻る

次の記事