作者evanzxcv (麥粒小螺絲)
看板Hiking
標題Re: [分享] 望山圖集
時間Wed Feb 24 08:48:34 2010
看了amjads大的模擬展望圖,我就在想怎麼沒有新竹地區的,
畢竟我從小在新竹長大的,三不五時就會去爬十八尖山,
山頂展望其實還不錯,就是樹木多了點.
求人不如求己,與其請amjads大幫忙做,不如自己弄一張!
新竹十八尖山東南180度展望:
http://web.mit.edu/ypwang/maps/shibajianshan.png
後續更新:
2/28 玉山主峰360度展望:
http://web.mit.edu/ypwang/maps/yushan.png
3/01 龜山島401高地190度展望:
http://web.mit.edu/ypwang/maps/gueishan.png
==以下是心得/插曲/技術性閒聊==
我使用的高程資料是NASA ASTER GDEM,解析度為一角分(約30m),
程式全部用MATLAB寫,跑下去才發現要把這東西弄到在PDA上即時繪製是不太可能的:
首先把高程資料以雙線性內插法轉成以觀測點為原點的極座標投影,
(極座標網格解析度設定:方位角0.02度,距離50m)
原本打算把全台灣的高程資料整合成一個大矩陣然後一次轉好,
但出現了記憶體不足的問題,只好一塊一塊慢慢弄,這個就至少要半小時左右.
然後把極座標投影的高程轉成極座標投影的仰角矩陣,這個只要幾秒鐘,
接著就是把仰角矩陣轉成展望圖,這個大概要五分鐘左右.
最後就是把圖片手動加上山頭名稱,使用的是amjads大寫的試算表版山頭方位角計算器.
(理論上是可以把這個整合進程式裡,但如何過濾看不見的山頭、
如何避免文字標籤重疊,就不是我短時間內可以想到的了...)
在做玉山展望圖的時候,因為牽涉到近距離的山頭,
而我早有所聞SRTM的玉山比實際矮很多,不敢斷然確定ASTER的玉山高度是準確的,
保險起見就把高度找了出來,才3846m而已,整整矮了106m!
然後又找了附近的幾座三千公尺以上山頭,都比實際高度低50m上下,
於是就把玉山高度代入3900m進去,跑出來的圖還算滿準確的.
但數位資料畢竟還是有本身的限制--
我把十八尖山展望圖做出來以後,第一眼竟然找不到五指山! orz
後來順著鵝公髻山邊找下去才找到"勉強"呈五指狀的山頭....
(幸好鵝公髻山的形狀是準確的...)
然後大霸尖山的那一塊大石頭也不見了,圖上看起來就像一座很普通的山.
最後我想說的是手動加註解超費時費眼力的:
1.假山頭超級超級多!!
2.要確認方位角、仰角、距離三個參數,才能準確標上山名.
3.方位角計算器算出來的山頭位置,跟圖上的山頭位置往往有所出入,
也就是我們所看到的山頭並不一定是那座山真正的最高點.這通常跟稜線走向有關.
所以在加註解時不能光把線拉到準確的角度,還要手動適度調整.
4.如果稜線有很多條在同一個角度範圍,往往會漏掉很多山頭沒有標註,
例如:由西往東標註的過程中,有時候是標到"XX山東峰"才發現"XX山"漏掉了,如此之類.
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 18.95.6.42
1F：推 withdream:我覺得光是能做出來就超級強的了 02/24 09:09
2F：推 OscarJeff:推新竹跟十八尖山! (原po的網域好威orz) 02/24 09:10
3F：→ OscarJeff:看了一下 果然是金牌學弟@@ 興趣很廣泛阿 真不錯 02/24 09:19
4F：→ evanzxcv:其實只要有DEM資料,程式的部分不難寫 02/24 09:28
5F：推 amjads:Astgtm 那個高程資料位置本身就有偏差, 除非轉換前已經有設 02/24 21:30
6F：→ amjads:控制點進行較正, 否則山頭一定會有很多地方會對不到的, 02/24 21:31
7F：→ amjads:我的作法就是不較正, 但觀測點直接放在高程資料區域的最高 02/24 21:32
8F：→ amjads:點, 然後在將此座標與方位角計算同步就好了, 偏差會比較小 02/24 21:34
9F：→ amjads:另外我的 40m 本身就是用 Astgtm , 缺三貂角的部份用 80m 02/24 21:36
10F：→ amjads:下去補, 30m 解析度取 40m 空間樣本, 視域範圍就如同下面這 02/24 21:41
12F：→ amjads:十八尖山因為可以展望大霸, 就算用 30m 解析下去改出山形 02/24 21:48
13F：→ amjads:應該還是不夠美觀, 應該要用 10~20m 解析度才比較有看頭~ 02/24 21:50
14F：推 piebear:麻省理工IP @@ 拜一下~~~ 02/24 23:05
15F：→ evanzxcv:三貂角那一塊我是用目測經建版等高線圖手動加上去的... 02/25 00:06
16F：→ evanzxcv:花了我將近一個小時.. 02/25 00:07
17F：→ evanzxcv:大霸山形出不來也可能是我極座標解析度的關係,我弄成50m 02/25 00:10
18F：→ evanzxcv:(原始資料是經緯兩方30m*30m,我轉極座標0.02度*50m) 02/25 00:12
19F：→ evanzxcv:也就是半徑100km內切線方向35m*徑向50m來畫圖 02/25 00:14
20F：推 jijung:天啊! 麻省理工IP @@ 拜一下~~~ 看見強者! 02/25 18:21
※ 編輯: evanzxcv 來自: 18.95.6.17 (03/02 09:58)
21F：推 xhung:厲害... 03/05 23:34