主要考虑了太阳方位角对阴影的影响.docVIP

利用卫星影像量算建筑物高度的方法 石家庄市二中 田园诗 杨博洋 指导老师 董双发 汪冰2 摘要：利用卫星影像上建筑物与其阴影之间的几何关系，提出了高度纹理法和三角截面阴影法量算建筑物高度的方法。并以石家庄市电视塔等建筑物为例，运用多种方法分别进行量算，通过计算误差再加以比较：高度纹理法和三角截面阴影法比其它量算方法过程更为简单，计算精度更高。 关键词：卫星影像，建筑物高度，高度纹理法，三角截面阴影法 选题背景和研究目的 测量是进行科学研究、搞好工程技术的一项基础性工作，提高测量精确度、简化测量方法是测量界的共同追求。测量建筑物高度是测量工作中的一项重要内容。为寻找较好的测量方法，长期以来，科学家们进行了大量的实验研究。近年来越来越多的人利用高分辨率卫星影像进行建筑物高度估算的研究。Cheng和Thiel等利用SPOT全色图像中建筑物阴影直接估算建筑物高度，并在42座建筑物高度测量中获得3.69m的均方根误差[1]。何国金等也进行了同样估算，在以北京市为例的试验中，建筑物高度分级抽样验证准确率达80％以上，但由于SPOT全色图像分辨率相对较低，其介绍的方法已进行了简化，仅对太阳方位角和卫星方位角相差180°的情况进行了计算[2]。在此之后，谢军飞等研究太阳高度角、卫星高度角、太阳方位角、卫星方位角等的几何关系后又对上述方法进行了一些改进[3]。但上述方法阴影起点较难确定，存在着误差较大、计算方法和判断过程复杂的缺点。能不能找到一种简便易行、快速有效、误差较小的建筑物高度测量方法呢？ 本人姑父从事遥感航测工作，因此对测量和遥感航测新技术有较多的接触了解。通过初、高中阶段数学物理等课程学习，了解了测量、误差以及影响测量精确度因素等方面的知识。测量建筑物高度，有很多种方法，但都存在着误差，有的误差较大。因而产生了利用高分辨率卫星影像，运用建筑物高度与其阴影之间的几何关系，进行简易快速精确测量的想法。通过与核工业部遥感航测中心的董双发、 汪滨两位老师探讨和研究，在董老师、（韩）老师的启发和指导下，设计了利用比其它量算方法过程更为简单的高度纹理法和三角截面阴影法，来有效减小误差，提高计算精度。 二、原理与方法 地球上某个地点和某一时间，我们可以根据某个人的身高计算出他影子的长度。因为我们可以知道他所处地球纬度的地方当时太阳的高度角，再依据“勾股定理”中知道两角一边而求另一边，即可完成计算。卫星影像是卫星在太空中对地球表面拍摄的照片，目前高分辨率的卫星影像可以十分清楚地识别城市建筑物的大小和形状，对其高度也可以利用阴影的长度依据“勾股定理”来加以推算。 在前人研究的基础上，作者提出了直接利用影像上建筑物投影高度计算建筑物实际高度和利用建筑物阴影长度计算建筑物高度的方法。此方法很好的解决了卫星影像上阴影的长度不易确定而加大计算误差的问题，可以更直接、简便、准确地利用高分辨率卫星影像进行建筑物实际高度的计算。 （一）建筑物高度量算的一般方法 通常利用卫星影像上阴影计算建筑物高度的研究方法主要根据影像获取参数的不同，太阳高度角的不同等，可分为以下几种情况： 设H为建筑物实际高度，S为建筑物阴影长，M为图像上成像阴影长，ω为太阳高度角，θ为卫星高度角，β为太阳方位角，α为卫星方位角。 1.何国金等基于SPOT数据提出的阴影法[2] 其研究以北京市为例，以SPOT数据为基础，提出了卫星传感器（照相仪器）扫描方向与太阳照射方向不同的两种方法，如下： ⑴卫星传感器扫描方向与太阳照射方向相同 此种情况下，影像上的成像影长M可测，其与建筑物实际高度H的关系式为： M=H/tanω-H/tanθ ⑵卫星传感器扫描方向与太阳照射方向相反 此种情况下，建筑物实际阴影长S与M相同可测，其与建筑物实际高度H的关系式为： H＝M×tanω 2.谢军飞等基于IKONOS数据提出的阴影法[3] 谢军飞等在研究太阳高度角、卫星高度角、太阳方位角、卫星方位角等的几何关系后提出以下三种情况的计算方法。 ⑴太阳方位角和卫星方位角相等 H=M×tanω×tanθ/（tanθ- tanω） ⑵太阳方位角和卫星方位角相差180° H= S×tanω ⑶太阳方位角和卫星方位角差在0—180°之间 假设太阳和卫星在建筑物的同一侧，在平坦的场地上，设EF为墙面法线；作辅助线 AD平行于CO；EC为建筑物的实际阴影长度；ED为卫星成像的阴影部分；BC为建筑物的计算高度 H；∠SEC为太阳方位角β；∠SAC为卫星方位角α；∠CEB为太阳高度角ω；∠BAC为卫星高度角θ。 因此，在ΔAEC中， ∠ECA＝α-β 在ΔABC中， AC＝BC／tanθ＝H／tanθ EC＝BC／tanω＝H／tanω 根据上式和余弦公式 EA2= EC 2+AC2—2EC·AC·cos∠ECA 再根据正弦