Instalasi GMT pada Ubuntu 12.04

Kali ini saya akan berbagi cara menginstall GMT di Ubuntu 12.04. Untuk mendapatkan GMT kita bisa langsung dari website GMT. Berikut adalah langkah-langkahnya :

Continue

Konversi Shapefile Menjadi Format GMT

Pada kesempatan kali ini, saya akan berbagi cara mengkonversi data shapefile (SHP) menjadi format GMT agar bisa diplot menggunakan GMT. Pertama-tama, downloadlah software ogr2gui disini dan disini. Setelah itu ekstrak kedua file tersebut dan jadikan dalam satu folder. Selanjutnya jalankan software tersebut. Pada bagian source masukkan file SHP yang akan anda konversi dan pada bagian target pilih format GMT dan beri nama filenya. Selanjutnya klik tombol Execute. Selesailah sudah konversinya.

Continue

Plotting Sinyal Gempa Dengan GMT

Kali ini saya akan berbagi cara menampilkan sinyal gempa dengan GMT. Sebelumnya carilah event gempa yang anda inginkan, selanjutnya carilah sinyal dari event tersebut disini atau disini. Data yang anda download nantinya akan berformat SEED. Rubahlah menjadi format ASCII dan pisahkan menurut komponen-komponennya. Berikut adalah contoh data yang sudah saya rubah menjadi ASCII dan dipisah menurut komponenya. Data tersebut adalah data event gempa tanggal 1 Juli 2011 dari stasiun Lembang (LEM). Panjang data adalah 73 sekon. Episenter gempa terletak pada koordinat 108.82 BT - 8.19 LS. Selain itu dapat anda tambahkan pula data mekanisme fokal dari Global CMT. Berikut adalah data mekanisme fokal yang saya dapatkan :

108.82 -8.19 103 -3.72 5.55 -1.83 5.82 -1.70 0.29 24 104 -11 010177A          

Simpan dengan nama "jawa.gmt".

Episenter gempa terletak pada kolom 1 dan 2. Tetapi saya akan memplotnya pada koordinat di kolom 11 dan 12.

Sedangkan data episenter gempa yang saya dapatkan adalah :

2011    7    1    14    33    22.8  -8.19    108.82   103    4.6    4.6

Simpan dengan nama "maineq.dat".

Continue

Membuat Peta Anomali Gravitasi dan Magnetik Dengan GMT

Pada kesempatan kali ini, saya akan berbagi cara membuat peta anomali gravitasi dan magnetik. Sebelumnya downloadlah data gravitasi dan magnet dahulu. Ada banyak model data gravitasi dan magnetik yang bisa anda download. Contohnya seperti di TOPEX, EarthByte dan NOAA untuk data gravitasi, dan GEOMAG untuk data magnet. Pada contoh kali ini saya akan menggunakan data gravitasi dari Earthbyte dan data model magnet EMAG3 dari GEOMAG.

Sebelum anda buat petanya, pertama harus anda pahami format data yang akan anda download. Pada data gravitasi dari EarthByte, format datanya adalah netCDF dengan satuan mGal. Data tersebut dikali 10 oleh EarthByte, sehingga kita harus membaginya dengan 10 kembali agar didapat nilai aslinya. Untuk membagi data tersebut kita gunakan perintah "grdmath".

Sedangkan pada data model magnetik EMAG3 dari GEOMAG format datanya adalah ASCII. Sehingga anda harus melakukan gridding untuk menampilkannya di peta. Sebelum melakukan gridding, data tersebut harus difilter dengan perintah "blockmean", "blockmedian", atau "blockmode". Setelah itu baru dilakukan gridding dengan perintah "surface".

Continue

Membuat Peta Hyposenter Gempa Dengan GMT

Kali ini saya akan berbagi cara membuat peta hyposenter dengan GMT. Sebenarnya cara berikut tidak terbatas pada data gempa saja, tetapi bisa dilakukan pada semua data yang memiliki data elevasi . Sebelumnya buatlah data gempa seperti berikut :

 PDE    1973  01 02 005320.30  -9.85  117.43  66  5.5 mbGS    ... .......     
 PDE    1973  01 02 022709.20   1.03  126.21  61  5.4 mbGS    ... .......     
 PDE    1973  01 07 221702.60   5.68  127.30  79  5.1 mbGS    ... .......     
 PDE    1973  01 09 061425.30   6.98  126.15  51  5.3 mbGS    ... .......     
 PDE    1973  01 13 111020.40  -2.70  101.27 105  5.2 mbGS    ... .......      

.............................................................

Simpan dengan nama "gempa.xyz".

Kolom yang kita pakai adalah kolom ke 7 dan 6 untuk koordinat, kolom ke 8 untuk data elevasi, dan kolom 9 untuk klasifikasi warna berdasarkan magnitude. Anda dapat merubah format data diatas, tetapi anda juga harus merubah script "gawk" agar sesuai dengan kolom yang anda buat.

Continue

Seleksi Data Spasial Dengan GMT

Pada kesempatan kali ini, saya akan menunjukkan cara melakukan seleksi data berdasarkan koordinatnya. Saya menggunakan data gempa di Indonesia dan akan saya beri warna merah untuk gempa di laut dan kuning untuk gempa di darat. Serta warna hijau untuk gempa dengan radius 750 km dari kota Palu, dan warna putih untuk gempa dengan batas poligon di Sumatera. Untuk melakukan seleksi data, perintah yang digunakan adalah "gmtselect".

Continue