Re: [分享] 望山圖集

看板Hiking (登山)作者 (麥粒小螺絲)時間16年前 (2010/02/24 08:48), 編輯推噓6(6015)
留言21則, 7人參與, 最新討論串2/2 (看更多)
看了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

02/24 09:09, , 1F
我覺得光是能做出來就超級強的了
02/24 09:09, 1F

02/24 09:10, , 2F
推新竹跟十八尖山! (原po的網域好威orz)
02/24 09:10, 2F

02/24 09:19, , 3F
看了一下 果然是金牌學弟@@ 興趣很廣泛阿 真不錯
02/24 09:19, 3F

02/24 09:28, , 4F
其實只要有DEM資料,程式的部分不難寫
02/24 09:28, 4F

02/24 21:30, , 5F
Astgtm 那個高程資料位置本身就有偏差, 除非轉換前已經有設
02/24 21:30, 5F

02/24 21:31, , 6F
控制點進行較正, 否則山頭一定會有很多地方會對不到的,
02/24 21:31, 6F

02/24 21:32, , 7F
我的作法就是不較正, 但觀測點直接放在高程資料區域的最高
02/24 21:32, 7F

02/24 21:34, , 8F
點, 然後在將此座標與方位角計算同步就好了, 偏差會比較小
02/24 21:34, 8F

02/24 21:36, , 9F
另外我的 40m 本身就是用 Astgtm , 缺三貂角的部份用 80m
02/24 21:36, 9F

02/24 21:41, , 10F
下去補, 30m 解析度取 40m 空間樣本, 視域範圍就如同下面這
02/24 21:41, 10F

02/24 21:42, , 11F
些圖一樣: http://ppt.cc/,uug
02/24 21:42, 11F

02/24 21:48, , 12F
十八尖山因為可以展望大霸, 就算用 30m 解析下去改出山形
02/24 21:48, 12F

02/24 21:50, , 13F
應該還是不夠美觀, 應該要用 10~20m 解析度才比較有看頭~
02/24 21:50, 13F

02/24 23:05, , 14F
麻省理工IP @@ 拜一下~~~
02/24 23:05, 14F

02/25 00:06, , 15F
三貂角那一塊我是用目測經建版等高線圖手動加上去的...
02/25 00:06, 15F

02/25 00:07, , 16F
花了我將近一個小時..
02/25 00:07, 16F

02/25 00:10, , 17F
大霸山形出不來也可能是我極座標解析度的關係,我弄成50m
02/25 00:10, 17F

02/25 00:12, , 18F
(原始資料是經緯兩方30m*30m,我轉極座標0.02度*50m)
02/25 00:12, 18F

02/25 00:14, , 19F
也就是半徑100km內切線方向35m*徑向50m來畫圖
02/25 00:14, 19F

02/25 18:21, , 20F
天啊! 麻省理工IP @@ 拜一下~~~ 看見強者!
02/25 18:21, 20F
※ 編輯: evanzxcv 來自: 18.95.6.17 (03/02 09:58)

03/05 23:34, , 21F
厲害...
03/05 23:34, 21F
文章代碼(AID): #1BX7Tb_r (Hiking)
討論串 (同標題文章)
文章代碼(AID): #1BX7Tb_r (Hiking)