Menghitung Jarak Dengan GMT

Akhirnya bisa posting lagi setelah beberapa hari nggak ada postingan. Maklum, nemu mainan baru di laptop.. :D . Oke, pada postingan kali ini, saya akan berbagi cara menghitung jarak dengan GMT. Jika pada ArcGIS kita tinggal tarik-tarik garis antar titik yang mau diukur jaraknya, tidak seperti itu pada GMT. Kita perlu menggunakan perintah grdmath untuk menghitung jarak lingkaran dari node-node pada grid terhadap suatu koordinat (lintang, bujur) dalam derajat, sehingga perlu kita kali dengan 111.13 agar menjadi kilometer (1 derajat = 111.13 km). Selanjutnya menggunakan perintah grdtrack untuk melakukan sampling file grid 2D sepanjang suatu garis (dalam hal ini jarak antar titik).

Berikut adalah scriptnya :

#Parameter gempa
set lon=97
set lat=-2
set depth=12 Km
set mag=5.0 mB
set OT=00:01:40 UTC
set Tanggal=12 Mei 2013

#Set parameter peta
set F=gempa.ps
set R=85/107/-6/6

#Peta dasar
psbasemap -JM25 -R%R% -Ba5g5WSne:: -K -Y5.5> %F%
makecpt -Cglobe -Z > elev.cpt
grdimage indo.nc -R -JM -Celev.cpt -K -O -Na -P >> %F%
pscoast -JM -R -K -O -P -Dh -Ir/0.1p,blue -Lf88/-5.2/17/500+l+jr --LABEL_FONT_SIZE=13p -Tf105.3/4.5/0.75i/2 --HEADER_FONT_SIZE=11p

--HEADER_OFFSET=0.05i >>%F%
psxy -R -JM -W10.0 -Sf0.8i/0.1ilt -Dh -Gblack -m -O -K trench.gmt>> %F%
psxy -R -JM -Wthicker -O -K transform.gmt -Dh -m >> %F%
psxy -R -JM -Wthinner,- -m -O -K GSFML_SF_FZ_KM.gmt >> %F%
psxy -R -JM -Wthinnest,- -m -O -K GSFML_SF_ER_KM.gmt >> %F%

#Plot jarak episenter
echo 95.3 5.55 B.ACEH LT > banda.dat
echo %lon% %lat% >> banda.dat
psxy banda.dat -R -JM -O -K -Wthick,red,- >> %F%
echo 98.7 3.59 MEDAN LT > medan.dat
echo %lon% %lat% >> medan.dat
psxy medan.dat -R -JM -O -K -Wthick,red,- >> %F%
echo 100.36 -0.95 PADANG LT > padang.dat
echo %lon% %lat% >> padang.dat
psxy padang.dat -R -JM -O -K -Wthick,red,- >> %F%
echo 102.3 -3.79 BENGKULU LT > bengkulu.dat
echo %lon% %lat% >> bengkulu.dat
psxy bengkulu.dat -R -JM -O -K -Wthick,red,- >> %F%

#Membuat kotak merah pada kota
echo 95.3 5.55 98 -7.7 B.ACEH LT Km kota dari> city.d
echo 98.7 3.59 98 -8.2 MEDAN LT Km kota dari>> city.d
echo 100.36 -0.95 98 -8.7 PADANG LT Km kota dari>> city.d
echo 102.3 -3.79 98 -9.2 BENGKULU LT Km kota dari>> city.d
psxy -R -J -O -K -Ss0.2 -Gred -Wthinnest city.d >>%F%
gawk "{print $1, $2, 9, 0, 1, $6, $5}" city.d | pstext -R%R% -J -O -K -Dj0.15/0 -N >> %F%

#Sampling dan menghitung jarak episenter
grdmath -Rg -I1 %lon% %lat% SDIST 111.13 MUL = dist.grd
grdtrack city.d -Gdist.grd > track.d
gawk "{print  $3, $4, 12, 0, 1, $6, int($NF+0.5), $7, $9, $8, $5}" track.d | pstext -R%R% -JM -K -O -N>> %F%

#Parameter gempa
echo 94 7 18 0 1 LT Recent Earthquake | pstext -R%R% -JM -K -O -N>> %F%
echo 85 -7 18 0 1 LT Parameter : | pstext -R%R% -JM -K -O -N>> %F%
echo 85 -7.7 12 0 1 LT Date = %Tanggal% | pstext -R%R% -JM -K -O -N>> %F%
echo 85 -8.2 12 0 1 LT Origin Time = %OT% | pstext -R%R% -JM -K -O -N>> %F%
echo 85 -8.7 12 0 1 LT Longitude = %lon% | pstext -R%R% -JM -K -O -N>> %F%
echo 85 -9.2 12 0 1 LT Latitude = %lat% | pstext -R%R% -JM -K -O -N>> %F%
echo 85 -9.7 12 0 1 LT Magnitude = %mag% | pstext -R%R% -JM -K -O -N>> %F%
echo 85 -10.2 12 0 1 LT Depth = %depth% | pstext -R%R% -JM -K -O -N>> %F%

#Plot data gempa berbentuk bintang
echo %lon% %lat% 10 0 1 > datagempa.dat
psxy -R -JM -Sa0.7c -W1 -Gyellow -O datagempa.dat >> %F%

#Print Peta
ps2raster %F% -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -Fjarak

Sekian untuk postingan kali ini. Selamat mencoba dan semoga bermanfaat. Salam Orang Indonesia...