作者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