2014美國地球物理聯盟年會
2014AGU(American Geophysical Union Fall meetings), Grand Canyon and Canyonlands National Park
2014 Konvensyen Pakatan Geofizikal Amerika, Taman Negara Grand Canyon dan Canyonlands
2015年1月3日 星期六
2013 EGU VIENNA, 歐洲地球科學聯合會(維也納)
2013年EGU歐洲地球科學聯合會的一些照片
Pictures of 2013European Geoscience Union meetings
Gambar yang menyertai 2013 European Geoscience Union meetings
2012年12月7日 星期五
利用GMT的點(X,Y)如果有Z值,如何畫成contour/Use coordinate to plot the contour map by coordinate point and z value./Pakai detik koordinat melukis suatu peta contour sama data nilai Z.
Melukis Peta contour, sama dengan fungsi blockmean, blockmode, triangulate, nearneigbor dan lain-lain
分別由X,Y,Z的資料,Z資料可以為深度,高度等等
XYZ data represented the coordinate (X,Y,Z)
Data XYZ bermakna detik koordinat (X,Y,Z)
沒有加入色調的等深線圖
no color bathymetry map
Peta dalam dasar laut tidak berwarna
no color bathymetry map
Peta dalam dasar laut tidak berwarna
加入等深線色調的圖
with color pallett
Peta yang ber warna
with color pallett
Peta yang ber warna
Triangulate
 | ||||||||||
|
程式碼/script/skript
blockmean petroleum_game.csv -I0.1 -R0/20/0/30 -V > nf_filt.dat
#triangulate petroleum_game.csv -Gpetroleum_game1.grd -I0.2 -R0/20/0/30 -V
#nearneighbor petroleum_game.csv -Gpetroleum_game1.grd -I0.2 -N2 -R0/20/0/30 -S50 -V
surface nf_filt.dat -Gpetroleum_game1.grd -I0.1 -R0/20/0/30 -S25 -T0.2 -V
makecpt -Crainbow -T-6900/-4500/100 -Z -V > nf.cpt
psbasemap -JX4/6i -R0/20/0/30 -Y-3i -B10f1g1nwSE -V -K > Zuolan_formation.ps
grdimage petroleum_game1.grd.grd -JX4i/6i -R -Cnf.cpt -Inf_shade.grd -V -K -O >> Zuolan_formation.ps
psxy petroleum_game.csv -Sc0.05 -Gblack -R -J -V -K -O >> Zuolan_formation.ps
grdcontour petroleum_game1.grd -JX4/6i -B10f1g1nwSE -R -C50 -A200 -V -K -O >> Zuolan_formation.ps
pause
#triangulate petroleum_game.csv -Gpetroleum_game1.grd -I0.2 -R0/20/0/30 -V
#nearneighbor petroleum_game.csv -Gpetroleum_game1.grd -I0.2 -N2 -R0/20/0/30 -S50 -V
surface nf_filt.dat -Gpetroleum_game1.grd -I0.1 -R0/20/0/30 -S25 -T0.2 -V
makecpt -Crainbow -T-6900/-4500/100 -Z -V > nf.cpt
psbasemap -JX4/6i -R0/20/0/30 -Y-3i -B10f1g1nwSE -V -K > Zuolan_formation.ps
grdimage petroleum_game1.grd.grd -JX4i/6i -R -Cnf.cpt -Inf_shade.grd -V -K -O >> Zuolan_formation.ps
psxy petroleum_game.csv -Sc0.05 -Gblack -R -J -V -K -O >> Zuolan_formation.ps
grdcontour petroleum_game1.grd -JX4/6i -B10f1g1nwSE -R -C50 -A200 -V -K -O >> Zuolan_formation.ps
pause
也可嘗試利用triangulate 或nearneighbor, blockmean 也可改成 blockmode, blockmedian等.
You can try another function like triangulate, nearneighbor, blockmean, blockmode, blockmedian
Anda boleh menggunakan fungsi yang lain bagi memuaskan peta yang lukis. Seperti triangulate, nearneighbor, blockmean, blockmode, blockmedian
You can try another function like triangulate, nearneighbor, blockmean, blockmode, blockmedian
Anda boleh menggunakan fungsi yang lain bagi memuaskan peta yang lukis. Seperti triangulate, nearneighbor, blockmean, blockmode, blockmedian
2012年11月23日 星期五
地形圖鋸齒狀的修正和精化/smooth the topographic by low resolusion DEM/ Menghaluskan peta topographic dari resolusi rendah DEM
地形圖原始資料的解析度不好,或者是有鋸齒狀,修改的的方法則是在grdimage 中加入-Evalue,例如說-E10, -E50, -E100
Low resolution of DEM raw data, we can use grdimage -Evalue, eg. -E10, -E50, -E100 to smooth the rugged boundary.
Data DEM resolusi rendah boleh menggunakan fungsi grdimage, tambah nilai -Evalue, misalnya. -E10, -E50, -E100 bisa menghaluskan garisan yang bengkang-bengkok.
圖1 為原始讀入的DEM資料,圖二和圖三則是在grdimage中分別加入-E50和-E100的地形圖
Figure 1 represent the DEM raw data, Figure2 and Figure 3 is using the function grdimage add the -E50 and -E100 value respectively.
Gambar 1 menunjuk Data DEM asli, Gambar 2 dan Gambar 3 pakai fungsi grdimage, tambah nilai -E50 dan -E100.
程式碼/ script/ skript
makecpt -Cglobe -T-199/200/50 -Z -V > milun.cpt
psbasemap -R121.58/121.65/23.95/24.06 -JM3i -B0.05f0.05NWes -G200/200/225 -V -K > magnifymilun.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -E100 -Cmilun.cpt -V -K -O >> magnifymilun.ps
psxy hualian01_1.txt -R121.58/121.65/23.95/24.06 -JM3i -Wwhite -G0/200/200 -I1 -V -O -K >> magnifymilun.ps
grdcontour "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -A20+w60 -C10 -G20c -L1/100 -Wared -V -O -K >> magnifymilun.ps
psxy meilun_river.asc -R121.58/121.65/23.95/24.06 -JM3i -Y0.1 -W10/blue -V -K -O >> magnifymilun.ps
psxy moea_fault.gmt -R121.57/121.65/23.95/24.06 -JM3i -MX -W5/red -Sf0.5/0.1irt -X0.1 -Y0.4 -Gred -V -K -O >> magnifymilun.ps
psscale -Cmilun.cpt -D5.1i/1.35i/2.88i/0.4i -Ac -B200:Topographic:/:m: -V -O >> magnifymilun.ps
Low resolution of DEM raw data, we can use grdimage -Evalue, eg. -E10, -E50, -E100 to smooth the rugged boundary.
Data DEM resolusi rendah boleh menggunakan fungsi grdimage, tambah nilai -Evalue, misalnya. -E10, -E50, -E100 bisa menghaluskan garisan yang bengkang-bengkok.
 | |||||||||
| 圖1/Figure1/Gambar1 |
.JPG) |
| 圖2/Figure2/ Gambar 2 |
.JPG) |
| 圖3/Figure3/Gambar3 |
Figure 1 represent the DEM raw data, Figure2 and Figure 3 is using the function grdimage add the -E50 and -E100 value respectively.
Gambar 1 menunjuk Data DEM asli, Gambar 2 dan Gambar 3 pakai fungsi grdimage, tambah nilai -E50 dan -E100.
程式碼/ script/ skript
makecpt -Cglobe -T-199/200/50 -Z -V > milun.cpt
psbasemap -R121.58/121.65/23.95/24.06 -JM3i -B0.05f0.05NWes -G200/200/225 -V -K > magnifymilun.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -E100 -Cmilun.cpt -V -K -O >> magnifymilun.ps
psxy hualian01_1.txt -R121.58/121.65/23.95/24.06 -JM3i -Wwhite -G0/200/200 -I1 -V -O -K >> magnifymilun.ps
grdcontour "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -A20+w60 -C10 -G20c -L1/100 -Wared -V -O -K >> magnifymilun.ps
psxy meilun_river.asc -R121.58/121.65/23.95/24.06 -JM3i -Y0.1 -W10/blue -V -K -O >> magnifymilun.ps
psxy moea_fault.gmt -R121.57/121.65/23.95/24.06 -JM3i -MX -W5/red -Sf0.5/0.1irt -X0.1 -Y0.4 -Gred -V -K -O >> magnifymilun.ps
psscale -Cmilun.cpt -D5.1i/1.35i/2.88i/0.4i -Ac -B200:Topographic:/:m: -V -O >> magnifymilun.ps
2012年7月18日 星期三
台灣東部北花蓮地形/Peta permukaan darat Hualien, timur Taiwan/Contour map nothern Hualien, eastern Taiwan
修改色盤的方法-利用word pad 開啟.cpt的色盤後修改.
menukarakan warna color pallett yang di pakai, buka kan file .cpt dengan word pad.
Change the color pallett by word pad,
根據以上的的設定,例如1公尺到5公尺範圍,設定RGB為相同的值,
mengikut contoh setting di atas, 1meter sampai 5meter ketinggian, pakai nilai RGB sama,
According the example, the RGB value are setting from 1 meter to 5 meter high with same value.
GTM 程式碼/Skript GMT/ Script GMT
makecpt -Cglobe -T-490/500/50 -Z -V > milun.cpt %移除-Z也可以變成非連續的色標, nilai -Z boleh dibuangkan kalau tidak mahu warna berterusan yang pakai untuk peta%
psbasemap -R121.5/121.65/23.84/24.1 -JM3i -B0.05f0.05NWes -G200/200/225 -V -K > shuilien.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.5/121.65/23.84/24.1 -JM3i -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -Cmilun4.cpt -V -K -O >> shuilien.ps
#grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.5/121.65/23.84/24.1 -JM3i -Csea1.cpt -V -K -O >> shuilien.ps
psxy hualian01_1.txt R121.5/121.65/23.84/24.1 -JM3i -Wwhite -G0/200/200 -I1 -V -O -K >> shuilien.ps
grdcontour "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.5/121.65/23.84/24.1 -JM3i -A20+w60 -C10 -G20c -L1/100 -V -O -K >> shuilien.ps
psxy moea_fault.gmt -R121.5/121.65/23.84/24.1 -JM3i -MX -W5/red -X0.1 -Y-0.1 -Sf0.5/0.1irt -Gred -V -O >> shuilien.ps
#pscoast -R121.5/121.65/23.85/24.06 -JM3i -Df -W2/0 -Na -V -O >> shuilien.ps
pause
%移除-Z也可以變成非連續的色表, nilai -Z boleh dibuangkan kalau tidak mahu warna berterusan yang pakai untuk peta%
makecpt -CCrainbow -T-5/55/5
psbasemap -R121.58/121.65/23.95/24.06 -JM3i -B0.05f0.05NWes -G200/200/225 -V -K > magnifymilun.ps
grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -I"/Numerical process/GMT_DTM_data/twdtm_small_temp_shade.grd" -C milun.cpt -V -K -O >> magnifymilun.ps
#grdimage "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -Csea1.cpt -V -K -O >> magnifymilun.ps
psxy hualian01_1.txt -R121.58/121.65/23.95/24.06 -JM3i -Wwhite -G0/200/200 -I1 -V -O -K >> magnifymilun.ps
grdcontour "/Numerical process/GMT_DTM_data/twdtm_small.grd" -R121.58/121.65/23.95/24.06 -JM3i -A20+w60 -C10 -G20c -L1/100 -Wared -V -O -K >> magnifymilun.ps
psxy meilun_river.asc -R121.58/121.65/23.95/24.06 -JM3i -Y0.1 -W10/blue -V -K -O >> magnifymilun.ps
psxy moea_fault.gmt -R121.57/121.65/23.95/24.06 -JM3i -MX -W5/red -Sf0.5/0.1irt -X0.1 -Y0.4 -Gred -V -O >> magnifymilun.ps
pause
2011年8月2日 星期二
海平面上升,台灣地區被海水淹蓋的地區/Paras laut meningkat, kawasan bajir Taiwan./Where flood in Taiwan by Climate change?
一般我們常常聽到排放溫室氣體造成溫度升高，南北極等冰川融化，海平面上升。其實地球公轉以及自轉也都會帶來氣候變化-米蘭科維奇的循環-
假設如果真的海平面上升，利用DEM來看上升10m與上升20m台灣哪些地區被淹沒?
Kita selalu dapat maklumat melepas gas-gas rumah hijau mengakibatkan cuaca semakin panas dan suhu meningkat lebih tinggi. Ais di Artik dan Antartik cair, paras laut meningkat.
Rotasi bumi atau revolusi bumi juga menyebabkan iklim dan cuaca berubah- dari Milankovitch kitaran- Kalau paras laut meningkat, kita boleh cuba menggunakan DEM buat simulasi dalam ketinggian 10 meter, 20 meter, dan dapat mengetahui di mana tempat yang dibanjiri, contoh yang pakai dibawah ialah DEM Negeri Taiwan.
海平面還未上升的情況
Gambar pertama- paras laut belum meningkat
Fig. 1 normal sea level
海平面上升10m台灣西部沿海鄉鎮(再加上地下水超抽),沉積物的壓密作用造成的沉降速率(沉積物隨著時間,重力而壓密,孔隙率變小, 例如美國常受到颶風侵襲的紐澳良地區是非常典型的例子)，台北盆地，高雄地區被淹沒
Gambar kedua- Paras laut meningkat 10meter, tempat yang berwarna biru gelap bermakna kawasan yang dibanjir, Kawasan yang berkenaan itu ialah, Basin Taipei, kawasan Kaoshiung dan plain Taiwan barat.
海平面上升20m台灣西部沿海，台北盆地，高雄地區，嘉南平原大部分人口稠密區淹沒
不過台灣還有很多地方還沒被海水淹沒。
Gambar ketiga- paras laut yang meningkat 20 meter, kecuali kawasan yang sebut didepan, Plain Chiayi-Tainan juga dibanjiri. Kawasan berikut mempunyayi populasi tinggi.Dalam gambar ini juga memberi maklumat yang penting kerana yang banyak lagi tempat lebih tinggi daripada paras laut selepas meningkat 20 meter.
從這三張圖海平面上升，人類必須面對的問題
1.遷徙
2.居住環境改變
3.改變生活型態
4.政府未來的政策
Dalam tiga gambar yang sebut diatas, masyarakat akan menghadapi masalah paras laut meningkat, jadi cara apa atau idea boleh digunakan.
1.meninggal kawasan yang rendah,
2.lingkungan dan kawasan hidup berubah.
3.Gaya hidup berubah
4.Cadangan atau dasar-dasar negeri masa depan.
GMT程式碼/Skript GMT
set R=-R118.5/123.5/21/26
set J=-JM4i
#makecpt -Cglobe -T-4000/4000/100 -Z -V > taiwan.cpt
psbasemap %J% %R% -B1nWSe -V -K > qixing_coast.ps
#grdgradient Taiwan_topo02.grd -GTaiwan_topo02_shape.grd -A60 -Ne0.6 -V
grdimage Taiwan_topo02.grd %J% %R% -Ctaiwan1_20m.cpt -V -K -O >> qixing_coast.ps %taiwan1_20m.cpt在這邊分別畫出不同海平面的調色盤，調色盤必需開啟檔案.cpt修改/Di sini pakai pallett warna paras laut yang tidak sama,file pallett warna itu format .cpt, jadi buka .cpt dari wordpad dan menukar warna yang anda pilih. %
psscale -D12c/5c/5c/0.5c -Ctaiwan2.cpt -I0.5 -B500 -V -O >> qixing_coast.ps
pause
修改.cpt如下，開啟.cpt檔案，如果要做10m海平面，之後多加一列1m到11m間，顏色則選擇同樣色調的，修改後儲存就可以/gambar di bawah menunakar warna di color pallett, buka file .cpt, kalau nak pilih 1meter sampai 10meter ketinggian, angka dalam RGB dalam 1meter sampai 10meter sila pilih angka yang sama seperti di barisan kedua di gambar bawah.
假設如果真的海平面上升，利用DEM來看上升10m與上升20m台灣哪些地區被淹沒?
Kita selalu dapat maklumat melepas gas-gas rumah hijau mengakibatkan cuaca semakin panas dan suhu meningkat lebih tinggi. Ais di Artik dan Antartik cair, paras laut meningkat.
Rotasi bumi atau revolusi bumi juga menyebabkan iklim dan cuaca berubah- dari Milankovitch kitaran- Kalau paras laut meningkat, kita boleh cuba menggunakan DEM buat simulasi dalam ketinggian 10 meter, 20 meter, dan dapat mengetahui di mana tempat yang dibanjiri, contoh yang pakai dibawah ialah DEM Negeri Taiwan.
海平面還未上升的情況
Gambar pertama- paras laut belum meningkat
Fig. 1 normal sea level
海平面上升10m台灣西部沿海鄉鎮(再加上地下水超抽),沉積物的壓密作用造成的沉降速率(沉積物隨著時間,重力而壓密,孔隙率變小, 例如美國常受到颶風侵襲的紐澳良地區是非常典型的例子)，台北盆地，高雄地區被淹沒
Gambar kedua- Paras laut meningkat 10meter, tempat yang berwarna biru gelap bermakna kawasan yang dibanjir, Kawasan yang berkenaan itu ialah, Basin Taipei, kawasan Kaoshiung dan plain Taiwan barat.
海平面上升20m台灣西部沿海，台北盆地，高雄地區，嘉南平原大部分人口稠密區淹沒
不過台灣還有很多地方還沒被海水淹沒。
Gambar ketiga- paras laut yang meningkat 20 meter, kecuali kawasan yang sebut didepan, Plain Chiayi-Tainan juga dibanjiri. Kawasan berikut mempunyayi populasi tinggi.Dalam gambar ini juga memberi maklumat yang penting kerana yang banyak lagi tempat lebih tinggi daripada paras laut selepas meningkat 20 meter.
從這三張圖海平面上升，人類必須面對的問題
1.遷徙
2.居住環境改變
3.改變生活型態
4.政府未來的政策
Dalam tiga gambar yang sebut diatas, masyarakat akan menghadapi masalah paras laut meningkat, jadi cara apa atau idea boleh digunakan.
1.meninggal kawasan yang rendah,
2.lingkungan dan kawasan hidup berubah.
3.Gaya hidup berubah
4.Cadangan atau dasar-dasar negeri masa depan.
GMT程式碼/Skript GMT
set R=-R118.5/123.5/21/26
set J=-JM4i
#makecpt -Cglobe -T-4000/4000/100 -Z -V > taiwan.cpt
psbasemap %J% %R% -B1nWSe -V -K > qixing_coast.ps
#grdgradient Taiwan_topo02.grd -GTaiwan_topo02_shape.grd -A60 -Ne0.6 -V
grdimage Taiwan_topo02.grd %J% %R% -Ctaiwan1_20m.cpt -V -K -O >> qixing_coast.ps %taiwan1_20m.cpt在這邊分別畫出不同海平面的調色盤，調色盤必需開啟檔案.cpt修改/Di sini pakai pallett warna paras laut yang tidak sama,file pallett warna itu format .cpt, jadi buka .cpt dari wordpad dan menukar warna yang anda pilih. %
psscale -D12c/5c/5c/0.5c -Ctaiwan2.cpt -I0.5 -B500 -V -O >> qixing_coast.ps
pause
修改.cpt如下，開啟.cpt檔案，如果要做10m海平面，之後多加一列1m到11m間，顏色則選擇同樣色調的，修改後儲存就可以/gambar di bawah menunakar warna di color pallett, buka file .cpt, kalau nak pilih 1meter sampai 10meter ketinggian, angka dalam RGB dalam 1meter sampai 10meter sila pilih angka yang sama seperti di barisan kedua di gambar bawah.
2011年7月21日 星期四
xyz2grd的功能以及GMT在windows的安裝/ Fungsi xyz2grd dan intall GMT dalam windows/Function xyz2grd and install GMT for windows
DEM具備有經緯度以及高層三個資料，如果我們從別處獲得資料後要將不同的資料格式給GMT讀必須要做轉檔的方法，例如 .XYZ 的資料格式
DEM ialah salah satu file yang memiliki 3 data, iaitu data longtitude dan latitude dan ketinggian. Kalau kita dapat data peta bukan GMT boleh dibaca, jadi kita mesti buat konversi file itu kepada format yang boleh dibaca, contohnya data format .XYZ.
利用Global mapper開啟的XYZ檔的DEM
然後選擇overlay control center進入conntrol center之後則選擇metadata
可看到以下的視窗
pakai software Global mapper, buka file format XYZ,file ini ialah DEM
selepas itu pilih overlay control center, masuk ke dalam control center dan pilih metadata,
sterusnya kita dapat melihat windows dibawah.
將Data 的Upper left ，Upper right，Lower left ，Lower right PIXEL WIDTH, PIXEL HEIGHT
抄下 最後將資訊寫入GMT程式碼內
Sila salinkan data Upper left, Upper right, Lower left dan Lower right PIXEL WIDTH, PIXEL HEIGHT, selepas itu maklumat yang salin itu tulis ke dalam skript GMT seperti dibawah.
GMT程式碼/Skript GMT
xyz2grd Taiwan_topo.xyz -GTaiwan_topo01.grd -R117/125/18/27 -I/0.0166666666666667/0.0166666666666667 -V %-I為PIXEL的大小/ -I bermakna saiz PIXEL%
pause
這樣就轉成了GMT可讀的grd檔了
Jadi sudah berjaya pernukaran format kepada file format grd, grd boleh dibaca oleh GMT.
另外則是轉貼/
Ini pautan yang macam mana install GMT dalam Windows.
GMT在window 下的安裝
DEM ialah salah satu file yang memiliki 3 data, iaitu data longtitude dan latitude dan ketinggian. Kalau kita dapat data peta bukan GMT boleh dibaca, jadi kita mesti buat konversi file itu kepada format yang boleh dibaca, contohnya data format .XYZ.
利用Global mapper開啟的XYZ檔的DEM
然後選擇overlay control center進入conntrol center之後則選擇metadata
可看到以下的視窗
pakai software Global mapper, buka file format XYZ,file ini ialah DEM
selepas itu pilih overlay control center, masuk ke dalam control center dan pilih metadata,
sterusnya kita dapat melihat windows dibawah.
將Data 的Upper left ，Upper right，Lower left ，Lower right PIXEL WIDTH, PIXEL HEIGHT
抄下 最後將資訊寫入GMT程式碼內
Sila salinkan data Upper left, Upper right, Lower left dan Lower right PIXEL WIDTH, PIXEL HEIGHT, selepas itu maklumat yang salin itu tulis ke dalam skript GMT seperti dibawah.
GMT程式碼/Skript GMT
xyz2grd Taiwan_topo.xyz -GTaiwan_topo01.grd -R117/125/18/27 -I/0.0166666666666667/0.0166666666666667 -V %-I為PIXEL的大小/ -I bermakna saiz PIXEL%
pause
這樣就轉成了GMT可讀的grd檔了
Jadi sudah berjaya pernukaran format kepada file format grd, grd boleh dibaca oleh GMT.
另外則是轉貼/
Ini pautan yang macam mana install GMT dalam Windows.
GMT在window 下的安裝
訂閱：
文章 (Atom)