2011年7月3日 星期日

利用地震資料來畫地震分布,深度以及規模/Pakai data gempa bumi dari katalog terbuka di internet


主要是利用地震資料(可以利用美國地震資料中心 http://earthquake.usgs.gov/earthquakes/eqarchives/epic/epic_rect.php,台灣地震科學中心TEC http://tecdc.earth.sinica.edu.tw/ 下載資料),選取所需要的經緯度,所下載的資料。資料將會有不同的屬性,其中會有經緯度,規模以及深度。再將這些資料利用GMT畫出來。

Pakai data Gempa bumi dari katalog terbuka(Pautan USGS amerika syarikat http://earthquake.usgs.gov/earthquakes/eqarchives/epic/epic_rect.php,Taiwan TEC http://tecdc.earth.sinica.edu.tw/ ) ,pilih keluasan longtitude dan latitude, download data. Dalam data yang dapat itu ada properti longtitude, latititude, kekuatan gempa bumi(magnitude), dalam(depth).Pakai GMT skript boleh baca file ini.

以下是GMT程式碼/ Skript GMT

makecpt -Crainbow -T0/300/5 -V -Z > depth.cpt %先建立調色盤/ tetapkan color pallett%
psbasemap -JM6i -R114/128/3/19 -B2f2NWes -V -K >philipines.ps
%設定圖框大小,經緯度範圍,投影方法等參數/ setting keluasan peta dan jenis projection%
psxy -JM6i -R114/128/3/19 eq_nu_phi2.txt -Sc -Cdepth.cpt -: -V -K -O >> philipines.ps
%讀入所下載的地震資料檔案,-Sc 為simbol的意思,c為circle/ baca data download dari katalog, -Sc bermakna simbol, c bererti bulat%
pscoast -JM6i -R114/128/3/19 -Dh -W3 -Na -Ia -T115.5/10.25/1.25 -L115.5/8.5/8.5/100k -V -K -O >> philipines.ps
%畫出海岸線, -Dh為線條的解析度高 ,另外還有Di,Dc,Dh,Df可選擇,-W海岸線的線條粗細;後面加/─類似-W3/100灰階色,也可以改成彩色的RGB的色階可寫成-W3/0/0/255;-T為指北,前兩個參數為經緯度,第三個為箭頭大小;-L則為比例尺前面三個為經緯度,第三個為所要設定的間隔為多少/Fungsi ini melukis garisan pantai, -Dh, Di, Dc atau Df boleh dipilih, ini ialah resolusi garisan yang hendak ditunjuk, -W bermakna garisan pantai tebal atau nipis, contoh nya -W3, belakangnya boleh tambah warna dalam jenis kelabu gred(gray grade), atau jenis RGB. -T bermakna penunjuk arah utara, nilai di petak depan dan kedua bermakna longtitude dan latitude, ketiga iaitu saiz tanda penunjuk; -L ialah skela yang ingin dipakai, petak depan dan kedua longtitude dan latitude juga, ketiga selang skela.%
psscale -D1.25c/14c/5c/0.25 -Cdepth.cpt -B50 -V -O >> philipines.ps
%表示深度的scale,將所建好的色盤讀入/skela dalam gempa bumi, baca color pallett yang setting didepan.%
pause

1 則留言:

  1. 您好,想請問您是如何將點畫作透明的呢?我參照您的指令提示調整了psxy,加了-:(不確定是不是這個),卻無法展點...
    謝謝您!

    回覆刪除