Membuat Peta Arah Vektor dengan GMT

Kali ini, saya akan membagikan cara membuat peta arah vektor. Peta ini bisa dipakai untuk menggambarkan arah dan besarnya pergeseran suatu titik. Disini saya akan menggunakan perintah psxy untuk menggambarkan arah vektor tersebut.

Pertama-tama, buatlah file yang berisi koordinat (bujur, lintang), arah azimuth dalam derajat, dan jarak pergeseran. Berikut ini contoh file yang saya gunakan dengan nama file data.dat :

101.7     -0.617   42    3.877
139.06    -2.69    33    0.433
139.34    -2.8     19    0.223
134       -2.24    11    2.268
137.96    -1.91    21    0.012
106.77    -7.22    7     2.294
137.4     -4.01    30    0.266
106.9     -7.06    31    2.52
139.05    -4.07    50    0.54
96.62      4.5     163   0.052
136.16    -2.93    178   0.154
..............................

 

Kemudian ketiklah perintah berikut :

set R=92/108/-8/8
set F=vektor.ps

psbasemap -JM13c -R%R%  -Ba3g3WSne -K > %F%
pscoast -JM -R -K -O -Dh -Glightgray -Slightblue -Wthinnest >> %F %
gawk "{print $1, $2, $3, $4*50}" data.dat | psxy -JM -R -S=0.15+e -Gred -W0.3 -O -K >> %F%

echo H 8 1 Keterangan > legenda
pslegend -Dx0i/0.5i+w1i/0.5i+jTL -JM -R -O -K -F+glightyellow legenda >> %F%

echo 92.3 -7.4 90 100 | psxy -JM -R -S=0.15+e -Gred -W0.3 -O -K >> %F%
echo 94.5 -7.4 2 mm | pstext -JM -R -O -F+f8 >> %F%

psconvert %F% -A -TG -P

Kemudian simpan dengan ekstensi .bat dan jalankan.

Pada bagian gawk "{print $1, $2, $3, $4*50}" data.dat | psxy -JM -R -S=0.15+e -Gred -W0.3 -O -K >> %F%, saya menggunakan nilai 50 sebagai pengali jarak. Data yang saya gunakan (data.dat) berisi data pergeseran dalam satuan mm. Sedangkan psxy membacanya dalam satuan km. Sehingga disini saya menggambarkan jarak 1 mm akan dibaca sebagai jarak 50 km oleh GMT.

Kemudian pada bagian echo 92.3 -7.4 90 100 | psxy -JM -R -S=0.15+e -Gred -W0.3 -O -K >> %F%, saya menggunakan nilai azimuth 90 derajat dan jarak 100 km yang akan saya jadikan jarak acuan pada legenda. Pastikan bagian -S=0.15+e pada saat kita menampilkan legenda sama dengan baris saat mengeplotnya pada peta.

Terakhir, pada bagian echo 94.5 -7.4 2 mm | pstext -JM -R -O -F+f8 >> %F%, saya menggunakan nilai 2 mm yang didapat dari jarak 100 km (legenda) dibagi dengan 50 km (pengali pada peta).

Selamat mencoba.