利用NOAA(美國國家海洋及大氣觀測局)的資料畫SST(海水表面溫度)圖/Peta suhu lapisan air laut teratas pakai GMT,data dari NOAA AMERIKA SYARIKAT

網路上有很多開放的性的資料庫，我們可以善加利用，如果要畫一張海水表面溫度圖，可以選擇美國國家海洋大氣觀測局的資料庫 http://ghrsst.nodc.noaa.gov/acessdata.html 的資料，但是在讓GMT畫圖的時候，我們必須對資料進行轉換.

Dalam internet, kita boleh dapat banyak data-data dari katalog atau FTP terbuka, Kita boleh guna katalog terbuka buat peta. Kalau nak buat sesuatu tajuk peta, seperti  suhu permukaan air laut dunia.
Di sini ialah satu contoh NOAA(National Oceanographic Data Center Amerika syarikat) katalog pautan http://ghrsst.nodc.noaa.gov/acessdata.html,Sebelum pakai data yang didapat dari katalog terbuka ini, kita harus buat konversi file format.

1.下載的介面


http://ghrsst.nodc.noaa.gov/acessdata.html > HTTP: http://data.nodc.noaa.gov/ghrsst/
>L4 > Glob > Navo > K10_SST > 可以選擇年度有資料天數的SST溫度資料，點選天數之後下載檔案格式為 ******SST.nc.bz2的資料

2.將資料下載完畢後解壓縮。然後開啟Global Mapper，選擇所解壓縮的.nc檔案
會出現一個視窗，選擇ANALYSED_SST


開啟之後選擇 >Tools > control centre 就會出現以下的畫面

點選圖層，滑鼠左鍵按兩下會出現 elevation options 的視窗 選擇 alter elevation values > Valid elvation range的Minimum data將無資料的點轉換成0。

最後則將.nc檔案轉成XYZ檔
流程為File> export Raster and elevation data > export XYZ grid > options > coordinate seperator > space, 就可以轉成XYZ檔了。

3.其餘的部份則是利用GMT來畫圖
GMT的程式碼如下

minmax analysised20110414.xyz %找出XYZ的最大和最小值%
gawk "{print $1, $2 , $3*0.1}" analysised20110414.xyz > sst01.xyz %Z溫度值為給予*0.1%
makecpt -Crainbow -T0/35/5 -V -Z > temp.cpt%製作調色盤%
xyz2grd sst01.xyz -Gsst.grd -I0.1 -R-180/180/-90/90 %將xyz值轉成grd檔案格式%
grdimage sst.grd -R-180/180/-80/80 -Ctemp.cpt -Jj150/0.025i -V -K > sst01.ps %將轉好的grd檔案讀入%
pscoast -R-180/180/-80/80 -Jj150/0.025i -B30 -G0 -O -V >> sst01.ps %畫海岸線，海岸以內-G填為黑色%
pause

不同投影法的台灣,加部分河川流域/Jenis projection yang lain dan tambah basin sungai

程式碼/ Skript GMT
set R= -R120.45/22.23/122.2/24.2r
set J=-Joc121.1/23.5/21/108/2.5i
makecpt -Cglobe -T-4000/4000/100 -Z -V > usa.cpt
makecpt -Cglobe -T-8000/8000/100 -Z -V > sea.cpt
psscale -D9c/5c/5c/0.5c -Cusa1.cpt -I0.5 -B1000 -V -K > qixing_coast1.ps
psbasemap %R% %J% -B1f1nwES -G241/252/255 -V -K -O >> qixing_coast1.ps
pscoast %R% %J% -Na --HEADER_FONT=33 -T121.8/23.8/1 -Lf121.15/22.2/22.2/40k -V -O -K >> qixing_coast1.ps
pscoast %R% %J% -Df -Gc -V -K -O >> qixing_coast1.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" %R% %J% -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -Cusa.cpt -V -K -O >>
qixing_coast1.ps
pscoast %R% %J% -Q -V -O -K >> qixing_coast1.ps
pscoast %R% %J% -Df -W2/red -Na -V -O -K >> qixing_coast1.ps
psxy liwu.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy lele.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy jinlun.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy luye.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy mukua.asc %R% %J% -W5/red -V -O -K >> qixing_coast1.ps
psxy xuhai.asc %R% %J% -W5/red -X-0.02 -Y-0.02 -V -O -K >> qixing_coast1.ps
psxy shuilien.asc %R% %J% -W5/blue -V -O -K >> qixing_coast1.ps
psxy xiukuluan.asc %R% %J% -W5/blue -V -O -K >> qixing_coast1.ps
psxy madagida.asc %R% %J% -W5/blue -V -O -K >> qixing_coast1.ps
pause

3D台灣加構造線/ dimensi 3 dan garisan struktur Taiwan

GMT程式碼/Skript GMT/Script GMT
makecpt -Cglobe -T-11000/8000/100 -Z -V > sea.cpt
psscale -D20c/12.75c/8c/0.5c --ANNOT_FONT_SIZE=8 -B5/:meter: -E -Csea1.cpt -I0.5 -B2000 -V -K > qixing_coast2.ps
#grdsample "/Numerical process/GMT_DTM_data/Taiwan_topo01.grd" -L -I0.02 -T -R119/124/21/26 -Ql -fg -GTaiwan_topo01_1by1.grd -V
grdgradient Taiwan_topo01_1by1.grd -G"/Numerical process/GMT_DTM_data/Taiwan_topo01_1by1_temp_shade.grd" -A60 -Ne0.9 -V
grdview Taiwan_topo01_1by1.grd -R119/124/21/26/-20000/4000 -Jm1 -Jz0.000075 -N-20000/gray -I"/Numerical process/GMT_DTM_data/Taiwan_topo01_1by1_temp_shade.grd" -E140/40 -B1/1/10000:"meter":SEwnZ -Qi -Csea1.cpt -V -K -O >> qixing_coast2.ps
pscoast -R -Jm1 -Jz0.000075 -W3 -E140/40 -Df -Tx123/25/5 -Z10 -O -V -K >>qixing_coast2.ps
psxyz deformationfront01.xyz -R -Jm1 -Jz0.000075 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz deformationfront02.xyz -R -Jm1 -Jz0.000075 -W4,red,- -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz deformationfront03.xyz -R -Jm1 -Jz0.000075 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz lishanfault02.xyz -R -Jm1 -Jz0.000075 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -Z10 -O -V -K >> qixing_coast2.ps
psxyz Luzontrench.xyz -R -Jm1 -Jz0.000055 -W4,red,- -E140/40 -Z6500 -O -V -K >> qixing_coast2.ps
psxyz Ryukyutrench01.xyz -R -Jm1 -Jz0.000095 -W4,red,- -Sf0.6/0.1irt -E140/40 -Gred -Z-2750 -O -V -K >> qixing_coast2.ps
psxyz Ryukyutrench02.xyz -R -Jm1 -Jz0.000065 -W4,red -Sf0.6/0.1irt -E140/40 -Gred -Z-5600 -O -V -K >> qixing_coast2.ps
psxyz chuchifault.xyz -R -Jm1 -Jz0.000070 -W4/red -Sf0.5/0.1ilt -E140/40 -Gred -O -V -Z2100 -K >> qixing_coast2.ps
psxyz lishanfault.xyz -R -Jm1 -Jz0.000068 -W4/red -Sf0.5/0.1irt -E140/40 -Gred -Z3500 -O -V >> qixing_coast2.ps
pause

利用地震資料來畫地震分布,深度以及規模/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

將KML檔轉成其他格式的網站/Laman web untuk file peta konversi

如果要將某個特定地區(POI)利用GOOGLE EARTH繪製出來的polygon 或是point等,經過存檔之後為KML檔,如果需要轉檔,推介的網站為

Laman web ini untuk file pernukaran format lain, kalau ada polygon, detik atau garisan (POI), boleh cuba gunakan.


http://www.poieditor.com/poi_convert

2010年1月23日環境學院體驗營

圖1 體驗營活動位置圖

這次的野外我參與了上午的活動 , 這主要的路線主要是在奇萊鼻海崖和七星潭兩個部份,
這次參與體驗營的高中生都來自花蓮縣內的各高中職 , 主要是介紹奇萊鼻海崖上的珊瑚礁 , 米崙斷層 , 民意斷層 , 以及七星潭的海階和米崙礫岩.



圖2 奇萊鼻燈塔以及海階上的珊瑚礁

因日本人Yamaguchi. M.,以及Ota, Y, 在2004 年發表的一篇文章中提到有六個海階 , 而根據他們的研究認為花蓮到台東有許多海階 , 而這些海階的形成是因為地震造成地表抬升.再加上在震間的的抬升 , 故可推算現今距離海水面20-30公尺的位置 , 再加上珊瑚定年得到的資料為2000-4000年之間 , 故利用先在測量的震間地表抬升量 , 故可以推估地表抬升的速度.12345
而誘發此地區地震主要為米崙斷層造成了此海階的抬升.

接著之後則是沿著七星潭腳踏車步道沿線可以觀察到與柏油路面呈現垂直的破裂面 , 此破裂面可以推測為民意斷層活動所造成的.



圖3 此圖在七星潭所拍攝的露頭(上覆為米崙層 , 下覆為米崙礫岩)

然後則到達米崙礫岩的部份, 米崙礫岩和上部所疊覆的米崙層呈現一種交角不整合的情況 , 故兩者在沉積的時間上應該是不同的.

謝孟龍和鄧屬予在1994年的文章提到 , 在米崙礫岩內部所採到的漂木進行碳14定年發現米崙礫岩的年代老於5萬5千年.在米崙礫岩上部的米崙層大部分受到上部的珊瑚礁經過雨水的溶解的碳酸鈣而造成灘岩的地形.


圖4 灘岩

最後則是介紹米崙斷層以及斷層造成的斷層崖和斷層湖這兩種相對的地形. 以東華大學美崙校區內的湖和空軍基地內的湖可以發現米崙層沿線出現的這兩種地形,這種地形出現主要是平移斷層在小範圍呈現不連續的現象, 受壓的部份部份會形成壓力脊, 拉張的部份則會形成凹陷盆地, 而且地表會出現斷層崖或線性構造.


圖5 平移斷層在小範圍上會有不連續的現象，在受壓的區域可能形成壓力脊，在拉張處常形成凹陷盆地，而地表則會出現斷層崖或線性構造(例如斷裂)。

利用GMT -pscoast clip的方法/Satu fungsi dalam GMT- pscoast clip

主要是利用clip 的方法 , 將兩張圖疊在一起 , 以上為例台灣的底圖較世界的底圖解析度較高
故利用此方法做疊合 , 如果遇到的問題是解析度較好的台灣底圖沒有海底的資料 , 這樣海洋的部份則以調色盤的-1 - 1 之間變化範圍內的-1為起始點 ,如果在相同色盤的情況下,所出現的情況是解析度較高的台灣島較低處和海底地形較低處的顏色相似 , 這似乎在顯示地形高低的情況是不可能發生的 , 故做兩個色表再加以套疊.

Fungsi ini ialah melakukan cara clip, terutamanya pakai dua lapisan yang berkaitan atas bawah, diatas itu resolusi tinggi(resolution),dan dibawah itu lebih rendah tetapi resolusi tinggi lapisan itu data yang meliputi itu kecil daripada resolusi yang rendah itu. Kita mungkin hanya dapat permukaan dasar laut pixel saiz besar.

以下是程式碼/Skript GMT

set R=-R118/127/17/26 %設定範圍/ masukan keluasan%
set J=-JM4i %設定投影方法/ masukan projection yang mahu di pakai%
makecpt -Cpanoply -T-5500/4000/100 -Z -V > usa.cpt %設定色盤/ pilih color pallett %
makecpt -Crainbow -T-8000/2000/100 -Z -V > sea.cpt %設定色盤/ pilih color pallett %
psbasemap %R% %J% -B2f2nwES -V -K > qixing_coast1.ps
grdimage "/Numerical process/GMT_DTM_data/etopo2.grd" %R% %J% -I"/Numerical process/GMT_DTM_data/etopo2_temp_shade.grd" -Csea.cpt -V -K -O >> qixing_coast1.ps
grdcontour "/Numerical process/GMT_DTM_data/etopo2.grd" %R% %J% -A1000+w60+s4+ap -G10c -L-11000/-1000 -Wared -V -K -O >> qixing_coast1.ps %此為加入等高線的指令/ Fungsi ini tambah garisan kontor(contour)%
pscoast %R% %J% -Df -Gc -V -K -O >> qixing_coast1.ps %加入-Gc 之後 , 就可以進行剪圖的動作/ Fungsi ini ialah fungsi Clip, tambah simbol bendera(flag) -Gc dapat melakukan fungsi klip %
grdimage "/Numerical process/GMT_DTM_data/etopo2.grd" %R% %J% -I"/Numerical process/GMT_DTM_data/etopo2_temp_shade.grd" -Cusa.cpt -V -K -O >> qixing_coast1.ps
pscoast %R% %J% -Df -Gc -V -K -O >> qixing_coast1.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" %R% %J% -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -Cusa.cpt -V -K -O >> qixing_coast1.ps
pscoast %R% %J% -Df -W2 -Na -V -O >> qixing_coast1.ps
pause