Pada postingan sebelumnya, saya telah membagikan script untuk menampilkan peta slab subduksi. Kali ini saya akan berbagi cara menampilkan model slab subduksi tersebut secara 3 dimensi. Sebelumnya downloadlah data model slab subduksi dari USGS disini serta data topografi disini. Saya menggunakan GMT 5.1.0, untuk GMT versi sebelumnya mungkin akan ada command yang tidak akan berjalan.
Selanjutnya buatlah script berikut di notepad dan simpan dalam format ".bat".
set data_slab=sum_slab1.0_clip.grd
set data_topo=indo.nc
set output=slab3d.ps
grd2cpt %data_slab% -Cjet -Z> slab.cpt
set output=slab3d.ps
grd2cpt %data_slab% -Cjet -Z> slab.cpt
makecpt -Cglobe -Z > elev.cpt
grdview %data_slab% -Jx1 -R92/120/-12/6/-700/1 -Cslab.cpt -JZ5 -Qi720 -p100/15 -K -Bx2+l"Longitude" -By2+l"Latitude" -B+b -Bz200+l"Depth (km)" -Y5> %output%
grdview %data_topo% -Jx -R92/120/-12/6/-700000/1000 -Celev.cpt -JZ5 -Qi720 -p100/15 -K -O -t50>> %output%
psscale -D12/-1/17/0.4h -Cslab.cpt -B50:"Depth to Slab (km)": -O -Y-1>> %output%
ps2raster %output% -Tf -A -P
grdview %data_slab% -Jx1 -R92/120/-12/6/-700/1 -Cslab.cpt -JZ5 -Qi720 -p100/15 -K -Bx2+l"Longitude" -By2+l"Latitude" -B+b -Bz200+l"Depth (km)" -Y5> %output%
grdview %data_topo% -Jx -R92/120/-12/6/-700000/1000 -Celev.cpt -JZ5 -Qi720 -p100/15 -K -O -t50>> %output%
psscale -D12/-1/17/0.4h -Cslab.cpt -B50:"Depth to Slab (km)": -O -Y-1>> %output%
ps2raster %output% -Tf -A -P
Script diatas terdiri dari 6 perintah:
- grdcpt berfungsi untuk membuat file warna berdasarkan data grid dalam hal ini data slab subduksi.
- makecpt berfungsi untuk membuat file warna berdasarkan default warna dari GMT.
- grdview pertama berfungsi untuk menampilkan data slab subduksi. Bagian -R merupakan batas min dan max longitude, latitude, serta elevasi. Bagian -Qi720 artinya output peta nanti resolusinya 720 dpi. Bagian -p100/15 berfungsi untuk mengubah sudut pandang peta, dalam hal ini azimuth peta adalah 100 derajat dan sudut elevasinya adalah 15 derajat.
- grdview kedua hampir sama seperti grdview pertama, hanya karena data topografi satuannya adalah meter sedangkan data slab satuannya adalah kilometer, maka pada bagian -R elevasinya dikali 1000. Bagian -t50 artinya transparansi peta sebesar 50 %. Transparansi ini hanya akan berfungsi jika output peta berformat PDF, sehingga anda harus mengkonversinya menjadi format gambar secara manual.
- psscale berfungsi untuk menampilkan skala kedalaman slab.
-ps2raster berfungsi untuk menkonversi file menjadi format raster. Bagian -Tf artinya output file berformat PDF.
Sekian untuk kali ini, semoga bermanfaat.
Selamat Malam Pak Yosi Setiawan,
ReplyDeletePerkenalkan Nama saya I Dewa Made, salah satu mahasiswa yang mengambil major Geomatika khususnya di bidang kegempaan. Saya ingin bertanya terkait plotting 3D slab berserta aftershock distributionnya. untuk pengeplotan slab 3D saya sudah berhasil melakukannya namun ketika ditambahkan data aftershock distribution dalam bentuk format (long, lat dan depth) tidak muncul gambar aftershocknya Pak. pertanyaannya adalah:
1. Haruskan data saya dirubah ke grd file dengan command xyz2grd?
2. Apakah dengan format long lat dan depth juga bisa di plot ke 3D berserta slab dengannya?
Terima Kasih sebelumnya.