Re: [分享] 望山圖集

[evanzxcv (麥粒小螺絲)]

看了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
推 withdream:我覺得光是能做出來就超級強的了
推 OscarJeff:推新竹跟十八尖山!                (原po的網域好威orz)
→ OscarJeff:看了一下 果然是金牌學弟@@ 興趣很廣泛阿 真不錯
→ evanzxcv:其實只要有DEM資料,程式的部分不難寫
推 amjads:Astgtm 那個高程資料位置本身就有偏差, 除非轉換前已經有設
→ amjads:控制點進行較正, 否則山頭一定會有很多地方會對不到的,
→ amjads:我的作法就是不較正, 但觀測點直接放在高程資料區域的最高
→ amjads:點, 然後在將此座標與方位角計算同步就好了, 偏差會比較小
→ amjads:另外我的 40m 本身就是用 Astgtm , 缺三貂角的部份用 80m
→ amjads:下去補, 30m 解析度取 40m 空間樣本, 視域範圍就如同下面這
→ amjads:些圖一樣: http://ppt.cc/,uug
→ amjads:十八尖山因為可以展望大霸, 就算用 30m 解析下去改出山形
→ amjads:應該還是不夠美觀, 應該要用 10~20m 解析度才比較有看頭~
推 piebear:麻省理工IP @@ 拜一下~~~
→ evanzxcv:三貂角那一塊我是用目測經建版等高線圖手動加上去的...
→ evanzxcv:花了我將近一個小時..
→ evanzxcv:大霸山形出不來也可能是我極座標解析度的關係,我弄成50m
→ evanzxcv:(原始資料是經緯兩方30m*30m,我轉極座標0.02度*50m)
→ evanzxcv:也就是半徑100km內切線方向35m*徑向50m來畫圖
推 jijung:天啊! 麻省理工IP @@ 拜一下~~~ 看見強者!
※ 編輯: evanzxcv        來自: 18.95.6.17           (03/02 09:58)
推 xhung:厲害...