Friday, 10 May 2013

Membuat Peta Seismisitas Berdasarkan Kedalaman dan Magnitude

By May 10, 2013 18 comments
Mumpung masih bersemangat nulis di blog, kali ini saya akan berbagi cara membuat peta seismisitas berdasarkan kedalaman dan magnitude sehingga peta yang dihasilkan nantinya akan lebih menarik. Kali ini saya menggunakan data dari USGS. Langkah pertama, jelas anda harus mendownload data gempanya dulu. Berikut contoh data yang sudah didownload :
2013-05-10 01:06:05.07     3.750     99.224     10.0    4.2    mb    us
2013-05-09 22:32:24.32     1.684     93.226     23.1    4.7    mb    us
2013-05-09 21:26:35.81     2.526     91.399     99.3    4.4    mb    us
2013-05-09 19:30:11.51     1.948     97.970     21.7    4.5    mb    us
..................................................
dimana kolom 1 adalah tanggal, kolom 2 jam, kolom 3 lintang, kolom 4 bujur, kolom 5 kedalaman, kolom 6 magnitude, kolom 7 jenis magnitude dan kolom 8 sumbernya. Selanjutnya simpanlah data tersebut, disini saya simpan dengan nama "gempa.dat".


Kali ini kita akan menampilkan gempa menurut kedalaman dengan warna, serta gempa menurut magnitude dengan ukuran yang berbeda. Buatlah file warna seperti berikut :
0        255 0 0        100    255 0 0
100    255 255 0    300    255 255 0
300    0 255 0        800    0 255 0
lalu simpanlah dengan nama "tabel.cpt".
kolom 1 dan 5 adalah batas kedalaman yang akan kita bedakan warnanya, kolom 2,3,4 dan 6,7,8 adalah warna dalam format RGB.

Selanjutnya buatlah baris-baris perintah berikut dan simpan dengan ekstensi ".bat" :

pscoast -R90/100/0/7 -JM20c -Ba10g5 -Ggreen3 -S140/235/255 -Dh -N1-W1 -K -X3 -Y6> map.ps
gawk "{ print $4, $3, $5, $6*0.03}" gempa.dat | psxy -R -JM -O -K -Ctabel.cpt -Sci -Wthin -H >> map.ps
echo H 14 1 SEISMICITY OF SUMATERA > legenda
echo D 0 1p  >> legenda
echo N 3   >> legenda
echo V 0 1p   >> legenda
echo S 0.1i c 0.1i red 0.25p 0.2i Shallow depth (0-100 km)   >> legenda
echo S 0.1i c 0.1i yellow 0.25p 0.2i Intermediate (100-300 km)  >> legenda
echo S 0.1i c 0.1i green 0.25p 0.2i Very deep (\076 300 km)   >> legenda
echo V 0 1p   >> legenda
echo D 0 1p   >> legenda
echo N 7   >> legenda
echo V 0 1p   >> legenda
echo S 0.1i c 0.09i - 0.25p 0.3i M 3   >> legenda
echo S 0.1i c 0.12i - 0.25p 0.3i M 4   >> legenda
echo S 0.1i c 0.15i - 0.25p 0.3i M 5   >> legenda
echo S 0.1i c 0.18i - 0.25p 0.3i M 6   >> legenda
echo S 0.1i c 0.21i - 0.25p 0.3i M 7   >> legenda
echo S 0.1i c 0.24i - 0.25p 0.3i M 8   >> legenda
echo S 0.1i c 0.27i - 0.25p 0.3i M 9   >> legenda
echo V 0 1p   >> legenda
echo D 0 1p  >> legenda
pslegend -Dx4.75i/0.3i/8i/0.97i/TC -J -R -O -K -F legenda -Gwhite -X-2 -Y-2>> map.ps

Jangan bingung melihat baris-baris diatas, baris-baris dengan akhiran "legenda" adalah perintah untuk membuat legenda pada peta. Anda dapat merubahnya sesuka anda. Setelah dieksekusi, hasilnya akan seperti berikut :


Selamat mencoba dan semoga bermanfaat. Salam olahraga...

18 comments:

  1. Unknown10 May 2013 at 18:00

    jiah salam olah raga

    ReplyDelete
  2. Om BAMBS10 May 2013 at 18:07

    coba ukuran lingkaran pada skala 5, 6, 7 lebih identik lagi...soalnya 6 & 7 hampir sama ukuran fontnya yuh...kurang greget....

    ReplyDelete
  3. Yosi10 May 2013 at 20:45

    ketularan bung brama bang..haha
    oke pak, itu di baris ke 15 sampai 21 tinggal dirubah aja yg 0.09i sampai 0.27 i nya

    ReplyDelete
  4. Unknown11 May 2013 at 00:57

    wahh kren bray...

    gw pnh cb pke pallete'a, kl klasifikasi dangkal (0-100), menengah (100300km) kyk'a cptnya mst :
    0 255 0 0 100.1 255 0 0
    100.1 255 255 0 300.1 255 255 0
    300.1 0 255 0 800 0 255 0

    kl ga,, kdalaman 100 msk di menengah, 300 masuk di kategori gempa dalam.. cb di cek..

    ReplyDelete
  5. Yosi11 May 2013 at 13:06

    oh iya bener juga om eki, :D

    ReplyDelete
  6. Unknown13 May 2013 at 18:51

    mantab om... sukses teross..!!

    ReplyDelete
  7. Unknown13 May 2013 at 19:23

    bray.. mo nanya ni.. yang baris echo S dst.... itu cmd bwt legendnya aj ato bwt simbol yg di petanya jg?? mo cb bwt legend magnitude'a dlm klmpk nih..

    ReplyDelete
  8. Yosi13 May 2013 at 22:45

    itu cuman buat legendanya doang blay, yg di peta simbolnya didapat dari gawk | psxy di baris 2 itu, ukuran lingkarannya didapat dari $6*0.03, $6 itu kolom magnitudenya di kolom 6,
    itulah sebabnya kenapa di echo S ukuran mag 3 kubuat 0.09i soalnya mag 3 dikali 0.03 jadinya 0.09, begitu jg yg lain

    ReplyDelete
  9. Unknown14 May 2013 at 14:01

    oo... brarti cmd -Sci itu simbol circle dlm inch,, tz size'a diambil dr data $6 gt bray..

    brarti bsar simbol mag 3.2 = 3.2*0.03inch? ga sma persis dgn mag 3 y?

    ReplyDelete
    Replies
    1. Yosi15 May 2013 at 15:55

      betul sekali mas bro, mag 3.2=0.096 inch
      makanya klo misalnya mau dibikin kelas2 kedalamannya, saranku dipisah pas di file .dat nya

      Delete
  10. Unknown15 May 2013 at 10:21

    Ka, saya coba buat tapi gempa.dat nya kayak gak kebaca. Udah dibuat desimalnya sama juga msih ga kebaca. Kenapa ya ka?

    ReplyDelete
    Replies
    1. Yosi15 May 2013 at 15:56

      format data gempanya sama kyk di contoh?

      Delete
  11. Periwinkle_ina30 May 2013 at 09:08

    Udah bisa ka :)

    ReplyDelete
  12. Unknown18 February 2015 at 16:18

    This comment has been removed by the author.

    ReplyDelete
  13. Afifah Nur28 February 2016 at 23:04

    Mas Yosi saya sudah coba scrip di atas tapi hasilnya hanya peta sumatra jawa tanpa seismositas, dan legenda.
    mohon bantuaanya.
    salam

    ReplyDelete
  14. Yosi4 June 2016 at 22:53

    apakah program gawk sudah bisa jalan di komputernya sebelumnya?

    ReplyDelete
  15. Anonymous28 May 2020 at 13:42

    mas yosi, saya coba ketik script sama persis seperti atas tapi kok legendanya tidak muncul ya?

    ReplyDelete
  16. Madha22 September 2020 at 15:30

    Saya sudah mencoba, tetapi tidak muncul bagian legendanya. Gawknya dapat dijalankan. Tolong bantuannya?. Terima kasih.

    ReplyDelete

Search