Interpolasi Surface dan Clip Peta Dengan GMT

Untuk artikel kali ini, saya akan berbagi cara melakukan interpolasi dengan GMT. Di dalam GMT sendiri ada beberapa perintah untuk melakukan interpolasi, tapi kali ini yang akan saya bahas adalah dengan perintah surface. Beberapa perintah yang digunakan kali ini adalah blockmean dan surface. Blockmean digunakan untuk melakukan filter data. Selain blockmean bisa pula digunakan perintah blockmedian dan blockmode, semua ini dianjurkan sebelum melakukan perintah surface agar nantinya tidak terjadi kesalahan data. Keterangan lebih lengkap bisa dibaca di manual GMT bagian blockmedian, blockmean, dan blockmode.

Sebelumnya, buatlah file berikut :

106.845    -6.211    155.4676242    LT    Jakarta
106.743    -6.302    132.999889    LT    Ciputat
107.609    -6.914    172.0670772    LT    Bandung
106.793    -6.589    126.4239398    LT    Bogor
106.927    -6.919    112.5143378    LT    Sukabumi
106.15    -6.12    62.0182684    LT    Serang
106.106    -6.308    60.90117673    LT    Pandeglang
106.252    -6.564    67.98173778    LT    Lebak
105.991    -6.308    69.06197946    LT    Malingping
106.564    -6.983    83.71995025    LT    Pelabuhan Ratu
107    -6.233    194.3086206    LT    Bekasi
107.286    -6.306    291.4413973    LT    Krawang
108.287    -6.393    193.5966186    LT    Indramayu
108.567    -6.717    110.8217199    LT    Cirebon
108.478    -6.98    97.9745866    LT    Kuningan
108.357    -7.326    93.45444372    LT    Ciamis
108.2    -7.333    109.4457299    LT    Tasikmalaya
107.909    -7.248    127.2282693    LT    Garut
107.145    -6.818    152.8897652    LT    Cianjur
107.759    -6.571    284.604707    LT    Subang
107.953    -6.833    177.3904801    LT    Sumedang
107.443    -6.561    269.7958866    LT    Purwakarta

Simpan dengan nama "pga.xyz". Kolom pertama adalah bujur, kolom kedua adalah lintang, kolom ketiga adalah nilai percepatan tanah dalam gal, kolom keempat adalah indentasi nama tempat yang akan diplot nantinya yang merupakan kombinasi dari L, C, R (left, center, right) dan T, M, B (top, middle, bottom), kolom keempat adalah nama lokasi.

Data diatas hanya digunakan untuk contoh dan bukan data sebenarnya. Anda dapat merubahnya sesuai data yang anda miiliki.

Berikut adalah scriptnya :

#Setting parameter peta
set F=clip.ps
set R=105/109/-8/-5.5 

#Membuat file palette warna untuk laut
makecpt -Cglobe -Z > elev.cpt

#Menampilkan topografi dan bathymetri
grdimage indo.nc -JM9.5i -Celev.cpt -K -R%R% -Sn > %F%

#Membuat arah mata angin dan skala peta
psbasemap -R -JM -O -K -Lf105.6/-7.4/17/100+l+jt -Tf108.6/-5.8/0.9i/2 --LABEL_FONT_SIZE=14p --HEADER_FONT_SIZE=14p --HEADER_OFFSET=0.05i>> %F%

#Melakukan filter
blockmean pga.xyz -R%R% -I0.001 > mean.xyz

#Melakukan gridding
surface mean.xyz -R -I0.001 -Gdata.grd

#Melakukan clip (buka)
pscoast -R -JM -O -K -Dh -Gc >> %F%

#Membuat file palette warna untuk percepatan
grd2cpt data.grd -Cseis > accelero.cpt

#Menampilkan percepatan
grdimage data.grd -JM9.5i -O -K -Caccelero.cpt -R -Sn >> %F%

#Melakukan clip (tutup)
pscoast -R -JM9.5i -O -K -Q >> %F%

#Plot titik data PGA dan judul peta
psxy pga.xyz -R%R% -JM9.5i -B1f0.5:."Peta Percepatan Tanah Provinsi Jawa Barat": -Sc0.05i -Gblack -K -O --HEADER_FONT_SIZE=24p >> %F%

#Menampilkan nama lokasi data PGA
gawk "{print $1, $2, 10, 0, 1, $4, $5}" pga.xyz | pstext -R%R% -JM -K -O -N -Dj-0.5/-0.5 >> %F%

#Menampilkan skala percepatan
psscale -Caccelero.cpt -D3.5/1/8.5/0.4 -B80:"Percepatan (gal)": -O -X1.3 -Y1 -I -Dh -Al --LABEL_FONT_SIZE=14p>> %F%

#Mencetak peta
ps2raster %F% -Gc:\programs\gs\gs8.53\bin\gswin32c -Tj -P -FPGA

Beberapa hal yang perlu diperhatikan antara lain :

• Jarak spasi grid pada blockmean dan surface harus sama, disini saya gunakan jarak 0.001 derajat, jika anda mau merubahnya kedalam menit maka tambahkan m di akhir opsi -I , tambahkan c jika detik, e jika meter, k jika kilometer, i jika mil, dan n jika mil laut.
• Perintah pscoast pertama digunakan untuk melakukan clip di wilayah daratan. Jika anda ingin melakukan clip wilayah lautan, ubahlah opsi -Gc menjadi -Sc.
• Perintah pscoast kedua digunakan untuk menutup perintah clip, ditandai dengan opsi -Q.

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