多点最小二乘匹配的可变权阵阵列松弛算法.pdfVIP

第23卷 第 3期 绘 学 报 Vo1．23，No．3 ACTA GEODAETICA etCARTOGRAPH lCA SINICA AuE ． i994 !7 多点最小二乘匹配的可变权阵 阵列松弛算法 张祖勋 吴晓良 武汉删’绘科技大学’ 。oo ¨ 提【要l 本文详细介绍 了多点最小二采匹配的常规算法和 由阵列代数表示的阵 列松弛算法，并把 Rauhala限于等权的情形推广到不等权 (变权)的情形，给 出7实 现阵列松弛的步骤，最后从计算效率、匹配精度对常规算法，等权阵列松弛和不等 权阵列松 弛进 行 了比较、分析 ，揭 示 了变权 阵列松 弛法 巨大的实用潜 力。 【关键词】 多点 (整体)影像匹配 阵列松弛法 权阵 金字塔影像结构 1 前 言 多点最小二乘匹配 “ ·“、“’ ，”是在单点最小二乘匹配基础上发展而来。它是 一种 整体匹配方法 ，与其它整体匹配技术，如随机最优化途径 。 ，正则化影像匹 配 “”，神 经 元引导的随机立体匹配一样，都具有较高的可靠性。但它们共 同缺 点 是 计 算 效 率 极 低。 Rosenholm用低通滤波以加快收敛速度，但并不总是有效 ”，Li用多级格网法结合多级金 字塔影像匹配技术有效地提高了匹配速度，但与实际要求仍相遥远 ”‘，因此，欲提高最小二 乘匹配的效率，其关键在于对大型矩阵求逆 。 阵列代数是对矩阵、张量代数的扩展，是解决多维数据快速处理的一个有力数学工具。 自1968年 Rauhala提 出阵列代数 以来，它被用于许多领域 ” ”。它的基本理论载于 ISP 1976年 Helsinki大会论文集中 “A Revlew ofAarrayAlgebra 一文 “ 。 应用阵列代数处理数据一般应具备两个条件：一是系统大型矩阵线性可分I二是所处理的 数据位于规则格网的节点上 l990年 Rauhala揭示了非线性最小二乘和强健估计的阵列代数 理论 ”，进一步拓宽 了阵列代数应用范围。阵列代数应用于影像匹配可追溯到l977年 ““， 把阵列代数与有限元方法结台 “ “，用于最小二乘匹配，这与其他学者，如Barnard Rosenholm “，Helava… ，worbe1 ”，Ebnsr和 Helpke“提出的多点(或称整体) 匹配，无不有异曲同工之效，而 Raaha18由于采用阵列松弛，在计算效率上独 胜 一筹。本 文将对多点最小二乘匹配的常规算法 ，阵列松弛算法，以及如何用阵列松弛实现 作 详 细 讨 论，把 Rauhala只能针对等权的庸形扩展到不等权情形。最后给出实验算例，对常观算法、 等权阵列松弛，不等权阵列松弛的结果进行比较、分析。 本戈 1993年9月 9日收到．截稿 日期为1994年 1月 31日． 本文为国家自然科学基金资助项 目． 测 绘 学 报 23卷 2 多点最小二乘匹配的常规算法 设g·(， )， t(， )依次表示左右影像：对于左影像中任一 像 元 g。 ，Y·)及 它在右像影中的同名点go2，Y：)，我们有下列观测方程 (假设无辐射畸变存在)： I ( l，Y1)一 ￡(2，Yi)=押(xl，Y1) (1) n0 ，y )是影像噪声，如果在核线影像对上匹配，即只存在左右视差 ，则有； I ( ，)一g2( + 。， )=”( ， ) (2) 式 (2)中视差 。可由该像元邻近四个视差格网点的视差双线性 内插表示，见图 l。设节点 间的距离为1，像元P到节点 (f，，)的距离为 ，d ，0d，也 ≤l： 。=(1一d1)(1一dt) +dI(1一d2)