基于物方影像匹配的DEM聚类变化检测.pdfVIP

测绘信息与工程 Journal of Geomatics Feb．2OO8；33(1) 35 文章编号：1007—3817(2008)01-0035—02 中图分类号：P231 文献标志码：B 基于物方影像匹配的DEM聚类变化检测 冯甜甜 王 密 潘 俊 (武汉大学测绘遥感信息工程国家重点实验室，武汉市珞喻路129号，430079) 摘 要 提出了DEM的变化检测问题，采用了基于密度的聚类分析方法对匹配所得的初始变化区域进行分析， 实现了地表高程变化区域的检测。 关键词 数字高程模型；变化检测；铅垂线轨迹法；聚类分析；影像匹配 随着自然环境的作用、人类活动以及人与自然的交互作 构影像，即金字塔影像，建立3级左右的金字塔影像层即可 用，地表随着时间的变化而不断变化Ⅲ1]。因而，了解地表的 获得很好的效果。在匹配时，先在最上一层的影像上进行匹 变化情况具有十分重要的意义，只有在确认高程没有发生变 配，再在下一层影像上匹配，将在上一层获得的匹配位置作 化时，所有的二维变化检测方法才能更有效地被使用[_2]。变 为下一层影像匹配搜索范围的初值，从而减少匹配的搜索空 化检测的研究得到了国内外的普遍关注，并逐渐渗透到越来 间，实现从粗到精的处理过程，达到快速、准确的匹配目的。 越多的领域中Ⅲ3]。文献[4]认为，变化检测要充分利用已有 结合铅垂线轨迹法的特点，将已有 DEM的格网点坐标 的4D产品及各种多源数据，努力提高变化检测的准确性和 作为初始值，利用新航空影像及其外方位元素直接解算格网 自动性，并将变化检测由二维向三维扩展。本文将基于密度 点的真实高程值。其具体步骤如下： 的聚类分析方法应用于DEM的变化检测中，以提高变化检 1)已知DEM某一格网点的坐标为(X，Y， )，按精度 测的准确性。 要求确定高程搜索步距 △z。 2)由该点的平面坐标(X，y)以及可能的高程 1 获取初始变化区域 Z 一 Z o+ i…AZ(i Yl，…，0，1，2…) 本文基于铅垂线轨迹法获得 DEM上初始的高程变化 计算左右像点坐标(z：，Y：)和( ， )，以左片像点的坐标 点。铅垂线轨迹法VI L(Vertical Line Locus)是Kern厂解 (z ，Y )为例，其公式如下Ⅲ5]： 析测图仪DSR一11附加CCD相机构成的混合型数字摄影测 I z：一一厂口 (x—x )+b (y—y )+f (z 一z ) 量工作站中的影像相关器所使用的方法 。在中心投影条 口 (X—x )+ 6 (y—y )+f (z 一z ) 件下，空间任意一条铅垂线以中心投影的方式投射到影像上 l Y：：一厂 仍是一条直线，该铅垂线与地面的交点在影像上的构像必定 右片像点( ， )的计算与左片类似。 位于相应的投影差上。以此为基础在立体像对上搜索相应 3)分别以(z：，y：)和( ， )为中心，在左右影像上取 的像点，从而确定该点的高程。利用VLL法搜索同名点， 影像窗口，计算影像的相关系数 。 VLL影像匹配原理图如图1所示。 4)改变 i的值，重复上述步骤2)和3)，得到一系列相关 系数 ，P，，fD2，…，