一、水璉礫岩10號橋剖面成份進行分析

地球科學資料處理
本次所要決定研究區域為台灣東部海岸山脈北段的水璉礫岩
本研究主要針對的沉積礫石的成份和顆粒大小進行分析
本次將使用GMT，EXCEL， MATLAB 作為工具
因水璉礫岩的露頭甚多，本研就選用鹽寮10號橋剖面進行分析
首先十號橋剖面的大約所在位置如圖下

紅色框內的橫線為10號橋剖面大約所在的位置位在北緯23.8到24度，東經121.4到121.5之間所在的最高的高度為400多公尺。

以下為程式碼

project -C121.5/23.87 -E121.6/23.85 -G0.1 -N -Q > temp1.d
gawk "{print $1, $2}" temp1.d > temp2.d
grdtrack temp2.d -G"/Numerical process/GMT_DTM_data/etopo2.grd" -V > temp3.d
gawk "{print $1, $3}" temp3.d > temp4.d
makecpt -Crelief -T-8000/3000/1000 -Z -V > taiwan.cpt
gmtset BASEMAP_TYPE PLAIN
psbasemap -JM5i -R119/123/21/26 -B1nwSE -V -K > taiwan.ps
grdimage "/Numerical process/GMT_DTM_data/etopo2.grd" -JM5i -R119/123/21/26 -I"/Numerical process/GMT_DTM_data/etopo2_temp_shade.grd" -Ctaiwan.cpt -V -K -O >> taiwan.ps
psxy -JM5i -R119/123/21/26 eq_nu_phi1.txt -Sc -Cdepth.cpt -: -V -K -O >> taiwan.ps
psxy temp2.d -JM5i -R119/123/21/26 -W10/red -V -K -O >> taiwan.ps
psxy poly.txt -JM5i -R119/123/21/26 -Ss1 -W10/red -V -K -O >> taiwan.ps
pscoast -JM5i -R119/123/21/26 -Df -W3 -Na -Ia -T119.5/22.5/1.5 -L119.5/21.75/22/50k -V -K -O >> taiwan.ps
psxy temp4.d -JX3/2i -R121.5/121.6/-100/500 -Sc0.1 -B1/100nwSE -X6.5i -Y4i -V -K -O >> taiwan.ps
pause

以放大的圖形來看10號橋剖面的話，紅色框框將鹽寮10號橋剖面框出來，利用橫剖面可以發現10號橋剖面為一個峽谷地形，因為被河流所切穿，因此高差可接近100~200公尺之間。因此本研究得以溯溪進行剖面的測量和樣本的採集。另一方面則是etopo.grd 和twd_small.grd兩者的解析度不一樣，以水璉鹽寮10號橋地區這種小範圍的研究，利用twd_small.grd做出來的效果較佳。

程式碼

project -C121.588/23.86 -E121.588/23.88 -G0.01 -N -Q > temp1.dgawk "{print $1, $2}" temp1.d > temp2.d

grdtrack temp2.d -G"/Numerical process/GMT_DTM_data/twdtm_small.grd" -V > temp3.dgawk "{print $2, $3}" temp3.d > temp4.d

gmtset BASEMAP_TYPE PLAINmakecpt -Crelief -T-100/1000/100 -Z -V > taiwan.cpt

psbasemap -JM5i -R121.50/121.65/23.85/23.95 -B0.1nwSE -V -K > shuilien.ps

grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -JM5i -R121.50/121.65/23.85/23.95 -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -Ctaiwan.cpt -V -K -O >> shuilien.ps

psxy temp2.d -JM5i -R121.50/121.65/23.85/23.95 -W10/red -V -K -O >> shuilien.

pspsxy poly.txt -JM5i -R121.50/121.65/23.85/23.95 -Ss1 -W10/red -V -K -O >> shuilien.ps

pscoast -JM5i -R121.50/121.65/23.85/23.95 -Df -W3 -Na -Ia -T119.5/22.5/1.5 -L119.5/21.75/22/50k -V -K -O >> shuilien.ps

psxy temp4.d -JX3/2i -R23.86/23.88/0/400 -Sc0.1 -B1/100nwSE -X6.5i -Yi -V -K -O >> shuilien.ps

pause

最後再做3-d立體圖形來看10 號橋剖面
被切穿出現的的峽谷地形，紅色框框所框之位置為鹽寮十號橋剖面之位置。

程式碼
makecpt -Crainbow -T0/1000/200 -Z -V > taiwan.cpt

grdview "/Numerical process/GMT_DTM_data/twdtm_small.grd" -Jm50 -Jz0.0015 -R121.50/121.65/23.85/23.95/0/500 -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -Ctaiwan.cpt -Qi -E80/50 -B1/1/400 -V -K > taiwan.ps

psxy poly1.txt -JM5i -R121.50/121.65/23.85/23.95/0/500 -Ss1.5 -W10/red -V -K -O >> taiwan.ps

psscale -D20c/3/5/0.2 -Ctaiwan.cpt -B500 -V -O >> taiwan.ps

pause