影像匹配基本算法(摄影测量)

这篇具有很好参考价值的文章主要介绍了影像匹配基本算法(摄影测量)。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

影像匹配实质上是在两幅(或多幅)影像之间识别同名点,它是计算机视觉及数字摄影测量的核心问题。由于早期的研究一般使用相关技术解决影像匹配问题,所以影像匹配常常被称为影像相关。

根据影像匹配过程中匹配基元选择的不同,影像匹配方法可分为像方基元匹配和地面元匹配两类。像方基元匹配是以影像上的像元或特征(点特征、线特征、面特征)为待匹配基元,利用像方信息(灰度或特征向量等)进行匹配准则的定义与计算,在像方侧完成同名像点的搜索确定。像方基元匹配方法总体上可分为局部影像匹配方法和全局影像匹配方法两类。另外,对于局部或全局匹配方法得到的视差,通常还要进行视差优化处理,以滤除错误视差、填补视差空洞等。而地面元影像匹配是以已知物方平面坐标(X,Y)的地面点作为待匹配基元,影像匹配的目的是确定地面基元的物方高程Z。由于地面元影像匹配方法能直接确定地物的空间三维坐标, 其也被称作基于物方的影像匹配方法,在数字表面模型生成、数字线划图绘制等地理信息产品生产过程中具有重要的应用。

像方匹配基元中局部匹配方法:SAD、SSD、NCC\ZNCC、Census-Transform、Mutual Information、PatchMatch
像方匹配基元中全局匹配方法:Graph Cut、Belief Propagation、Dynamic Programming
像方匹配基元中半全局匹配方法:SGM

地面元匹配方法:PMVS、VLL/MVLL

当然啦,本文的重点是摄影测量里的影像匹配算法,虽然有些算法与前面局部匹配方法中的相似,但叫法不同。摄影测量中基于像方的匹配算法(常见的五种算法):1.相关函数法;2.协方差函数法;3.差平方和法;4.差绝对值和法;5.相关系数法。实验表明,相关系数法影像匹配算法的成功率最高,所以这里就简单介绍一下基于像方的匹配算法——相关系数法以及基于物方的匹配算法——铅垂线轨迹法(VLL-Vertical Line Locus),还有摄影测量中最常用的最小二乘影像匹配算法。

1.基于像方的匹配算法——相关系数法

相关系数法是标准化的协方差函数,协方差函数除以两信号的方法即得相关系数。g(x,y)与g’(x’,y’)的相关系数为:
影像匹配基本算法(摄影测量)
其中,
影像匹配基本算法(摄影测量)

由离散灰度数据对相关系数的估计为:
影像匹配基本算法(摄影测量)

考虑到计算工作量,相关系数的实用公式为
影像匹配基本算法(摄影测量)

或化为
影像匹配基本算法(摄影测量)

优点:相关函数是衡量左右影像两个灰度矢量X与Y相似性的一个数值指标。由于相关系数是标准化的协方差函数,因此当目标影像的灰度与搜索影像的灰度之间存在线性畸变时,仍能较好地评价他们之间的相似性程度。(相关系数时灰度线性变换的不变量)。

2.基于物方的匹配算法——VLL影像匹配算法

基于物方的影像匹配可以理解为高程直接解求的影像匹配,已知待定点平面坐标(X,Y),只需确定其高程Z,可以沿过点(X,Y)的铅垂线到左右影像上的投影直线进行匹配,这就是铅垂线轨迹法(Vertical Line Locus),当匹配完成时,点(X,Y)的高程也同时获得了。

假设在物方有一条铅垂线轨迹,则它的像片上的投影也是一条直线。这就是说VLL与地面交点A在像片上的构像必定位于相应的投影差位于的直线上。利用VLL法搜索其相应的像点a1与a2,从而确定A点的高程。这个过程与人工在解析测图仪或立体测图仪上的过程十分相似。
影像匹配基本算法(摄影测量)
其步骤为:(1)给定地面点的平面坐标(X,Y)与近似最低高程Zmin,高程搜索步距 Δ \Delta ΔZ可由所要求的高程精度确定。
(2)由地面点平面坐标(X,Y)与可能的高程Zi = Zmin + i * Δ \Delta ΔZ (i = 0,1,2,…),计算左右像坐标(xi’,yi’)与(xin,xin)
影像匹配基本算法(摄影测量)
(3)分别以(xi’,yi’)与(xin,xin)为中心在左右影像上取影像窗口,计算器匹配调度,如相关系数 ρ \rho ρ(也可以利用其他测度)。
(4)将i的值增加1,重复(2)(3)两步,得到 ρ \rho ρ0, ρ \rho ρ1,…, ρ \rho ρn,取其最大者 ρ \rho ρk,其对应高程为Zk = Zmin + k * Δ \Delta ΔZ,则认为地面点A的高程Z = Zk
(5)还可以利用 ρ \rho ρk及其相邻的几个相关系数拟合一条抛物线,以其极值对应的高程作为A的高程,以进一步提高精度,或者以更小的步距 Δ \Delta ΔZ在一个小范围内重复以上过程。

3.最小二乘影像匹配算法

最小二乘法在影像匹配中的应用是在20世纪八十年代发展起来的。德国Ackermann教授提出了一种新的影像匹配方法——最小二乘影像匹配。由于该方法充分利用了影像窗口内的信息进行平差计算,使影像匹配可以达到1/10甚至1/100像素的高精度,即影像匹配精度可达到子像素等级。

首先来介绍一下影响灰度的系统变形有两大类:辐射畸变、几何畸变。
辐射畸变如何产生:1.照明及被摄物体辐射面的方向;2.大气与摄像机物镜所产生的衰减;3.摄影处理条件的差异以及影响数字化过程中所产生的误差等等。
影像匹配基本算法(摄影测量)
几何畸变如何产生:1.摄影机方位不同所产生的影像的透视畸变;2.影像的各种畸变;3.由于地形坡度所产生的影响即便;4.竖直航空摄影的情况下,地形高差是几何畸变的主要因素。
影像匹配基本算法(摄影测量)

最小二乘的匹配原理:
影像匹配中判断影像匹配的度量很多,其中有一种是“灰度差的平方和最小”。若将灰度差记为余差v,则上述判断可以写为: Σ \Sigma Σvv = min。因此,它与最小二乘的原则是一致的。但是,在一般情况下,它没有考虑影像灰度中存在着系统误差,仅仅认为影像灰度只存在偶然误差(随机噪声n),即
影像匹配基本算法(摄影测量)
这就是一般的按 Σ \Sigma Σvv = min的原则进行影像匹配的数字模型。若在此系统中引入系统变形的参数,同时按最小二乘 Σ \Sigma Σvv = min的原则,解求这些参数,就构成了最小二乘匹配系统。

最小二乘匹配优点:
  1.最小二乘影像匹配中可以非常灵活地引入各种已知参数和条件,从而可以进行整体平差;
  2.最小二乘影像匹配即可以解决“单点”的影响匹配问题,以求得“视差”,也可以直接解求其空间坐标;
  3.最小二乘影像匹配同时解决“多点”影像匹配或“多片”影像匹配;
  4.最小二乘影像匹配可引入“粗差检测”,从而大大地提高影像匹配的可靠性;
  5.甚至还可以用于解决影像遮蔽问题。

3.2基于辐射改正的LSM(仅考虑辐射的线性畸变的最小二乘匹配——相关系数)
影像匹配基本算法(摄影测量)

有n个像素点时误差方程为:
影像匹配基本算法(摄影测量)

按照 Σ \Sigma Σvv = min的原理,可得法方程式
影像匹配基本算法(摄影测量)
最后解得:
影像匹配基本算法(摄影测量)
为什么说仅考虑辐射的线性畸变的最小二乘匹配与相关系数等价呢,接下来解释一下:
将解得的h0,h1带入下面公式,得到 Σ \Sigma Σvv
影像匹配基本算法(摄影测量)

已知相关系数的定义如下:
影像匹配基本算法(摄影测量)

所以将 Σ \Sigma Σvv用相关系数表示得到:
影像匹配基本算法(摄影测量)

以下式相关系数与信噪比之间的关系可见,以相关系数最大作为影像匹配搜索同名点的准则,其实质是搜索“信噪比为最大”的灰度序列。
影像匹配基本算法(摄影测量)

3.3仅考虑影像相对位移的一维最小二乘匹配

假设两个一维灰度函数g1(x),g2(x),除随机噪声外,g2(x)相对于g1(x)只存在零次几何变形——位移量 Δ \Delta Δx。
影像匹配基本算法(摄影测量)

最小二乘影像匹配是非线性系统,必须进行迭代,迭代过程收敛的速度取决于初值,因此必须已知初始匹配的结果。(先用相关系数匹配结果作为初始值)

3.4单点最小二乘匹配

两个二位影像之间的几何变形,不仅存在着相对位移,而且还存在着图形变化。只有充分考虑影像间的几何变形,才能获得最佳的影像匹配。
影像匹配基本算法(摄影测量)

但是由于影像匹配窗口的尺寸均很小,所以一般只考虑一次畸变。几何变形公式如下:
影像匹配基本算法(摄影测量)

若同时再考虑右影像相对于左影像的线性灰度畸变,可得
影像匹配基本算法(摄影测量)

线性化误差方程式:
影像匹配基本算法(摄影测量)

初始值分别为:
影像匹配基本算法(摄影测量)

接下来的误差方程及法方程为:
影像匹配基本算法(摄影测量)

其中矩阵C的值为:
影像匹配基本算法(摄影测量)

由于在数字影像匹配中,灰度均是按规则格网排列的离散阵列,且采样间隔为常数 Δ \Delta Δ,可被视为单位长度,偏导数均用差分代替。
影像匹配基本算法(摄影测量)

最小二乘匹配流程图如下所示:
影像匹配基本算法(摄影测量)

最小二乘影像匹配的迭代过程:
(1)几何变形改正;
(2)重采样;
(3)辐射畸变改正;
(4)计算左影像窗口与经过几何、辐射改正后的右方影像窗口的灰度阵列g1与h0+h1·g2(x2,y2)之间的相关系数 ρ \rho ρ,判断是否需要继续迭代;
(5)采用最小二乘影像匹配,解求变形参数的改正值dh0,dh1…。
(6)计算变形参数;
(7)计算最佳匹配的点位。

在计算最佳匹配点位过程中,匹配精度却决于一个影像灰度的梯度,用梯度的平方为权,在左影像窗口内对坐标作加权平均:
影像匹配基本算法(摄影测量)
以它作为目标点坐标,它的同名点坐标可由最小二乘影像匹配所求得得几何变换参数求得:
影像匹配基本算法(摄影测量)

为了进一步提高其可靠性与精度,例如,附带共线条件的最小二乘相关以及与VLL法结合的最小二乘影像匹配方法都得到了广泛的研究。

3.5最小二乘匹配精度

最小二乘匹配算法,可以根据 σ \sigma σ0以及法方程式系数矩阵的逆矩阵,同时求得其精度指标。
以最简单的一维最小二乘匹配为例:
影像匹配基本算法(摄影测量)

3.5多点最小二乘影像匹配

一般影像匹配包括人工观测在内,要想在缺乏纹理的区域完成匹配和观测是十分困难的,这是因为一般匹配方法每次只能匹配一个点。
如果采用双线性有限元内插与最小二乘匹配相结合的方法,同时答解各网格点的视差,再通过周围点视差平滑约束,可以求得这些地区(小范围纹理缺乏地区)的视差值,从而完成局部整体影像匹配,这就是多点最小二乘影像匹配。
文章来源地址https://www.toymoban.com/news/detail-445137.html

到了这里,关于影像匹配基本算法(摄影测量)的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 摄影测量(tip1):空间前方交会与空间后方交会

    1.1 定义 已知 A, B 平面坐标,在已知点 A, B 分别瞄准未知点 P,从而求出点 P 平面坐标 1.2 分类 -角度前方交会- 通过观测角 A 和 B 求出点 P 坐标的定位方法 -方位角前方交会- 通过观测方位角  和  求出点 P 坐标的定位方法  1.3 定位精度 点P的精度为 可见,为了减少测量误差,

    2024年02月10日
    浏览(28)
  • 目前最流行的无人机摄影测量软件有哪些?各有什么特点?

    推荐阅读 : 《无人机航空摄影测量精品教程》 包括:无人机航测外业作业流程(像控点布设、航线规划、仿地飞行、航拍)和内业数据处理软件(Pix4d、CC、EPS、PhotoScan、Globalmapper)像控点权重调配、空三加密、DOM、DSM、DEM,DLG生成、等高线生成、高程点提取,点云分类过滤

    2024年02月03日
    浏览(51)
  • 大疆精灵4A无人机航空摄影测量外业数据采集完整操作流程 - 点云处理

    无人机在航空摄影测量领域中发挥着重要的作用,能够高效地获取大范围地理信息数据。本文将介绍大疆精灵4A无人机的航空摄影测量外业数据采集的完整操作流程,并重点讨论点云处理的相关内容。以下是详细的操作步骤和源代码示例。 准备工作 在开始操作之前,需要完成

    2024年02月06日
    浏览(219)
  • 道路匹配MapMatching:HMM模型、维特比算法Viterbi、道路匹配基本算法ST、STD、IVVM算法介绍

    我曾经做过有关 道路匹配(MapMatching) 的相关研究,学习过几个重要的道路匹配算法,我将先对 重要的匹配模型:隐马尔科夫模型(HMM) 进行介绍,再介绍 维特比算法Viterbi ,最后对 ST、STD、IVMM 三种算法做一个简单的介绍,供大家参考。 隐马尔科夫模型(Hidden Markov Model,简

    2024年02月06日
    浏览(35)
  • C++医院影像科PACS源码:三维重建、检查预约、胶片打印、图像处理、测量分析等

    PACS连接DICOM接口的医疗器械(如CT、MRI、CR、DR、DSA、各种窥镜成像系统设备等),实现图像无损传输,实现DICOM胶片打印机回传打印功能,支持各种图像处理,可以进行窗技术调节,与登记台管理系统共享数据,提供智能化专家诊断知识库,提供图像相互转换功能(DICOM转换为

    2024年02月07日
    浏览(41)
  • 摄影入门 | 相机的基本原理

    小孔成像实验中,点燃蜡烛,会在小孔另一面的白纸上看到一个倒立的烛焰。 此现象可以用来解释物理学原理: 光在同种均匀介质中,在不受引力作用干扰的情况下沿直线传播 。 这样,我们就用一种最简单的方法在白纸上获得了蜡烛烛焰的图像。 如果将小孔成像实验中的

    2024年02月07日
    浏览(36)
  • 【视觉SLAM入门】5.1. 特征提取和匹配--FAST,ORB(关键点描述子),2D-2D对极几何,本质矩阵,单应矩阵,三角测量,三角化矛盾

    为什么重要?我们是在做什么事? 特征提取和匹配: 首先是两幅图像的特征提取,然后是对应特征点的匹配。接下来的工作是根据得到的匹配点对,估计相机的运动,具体根据相机分为三种方法: 单目相机:2D-2D: 对极几何 方法 双目或者RGBD相机: 3D-3D: ICP 方法 一个3D点

    2024年02月13日
    浏览(45)
  • Multisim软件常用仪表的使用与一些基本测量方法

    本文作者使用multisim14.0中文破解版 (图1) (图2) 函数信号发生器用法讲解: 函数信号发生器 可产生三类五种波形,分别是 正弦信号、三角波、锯齿波、方波和脉冲波 (图3) 右框可调单位,如振幅: (图4) 示例:脉冲波 1KHz 20% 0-5V的脉冲信号 步骤: (图5) 示波器用

    2024年02月07日
    浏览(42)
  • 医学影像PACS系统源码: 三维重建基本后处理方法的介绍和说明

    多层面重建(MPR) 最大密度投影(MIP) 最小密度投影(MinIP) 表面阴影遮盖(SSD) 容积漫游技术(VRT) 曲面重建(CPR) 虚拟内镜技术(VE)   1、MPR MPR(Multi-Planar Reformatting),多平面重建,是将扫描范围内所有的轴位图像叠加起来再对某些标线标定的重组线所指定的组织进行

    2024年02月09日
    浏览(57)
  • C#处理医学影像(四):基于Stitcher算法拼接人体全景脊柱骨骼影像

    在拍摄脊柱或胸片时,经常会遇到因设备高度不够需要分段拍摄的情况, 对于影像科诊断查阅影像时希望将分段影像合并成一张影像,有助于更直观的观察病灶, 以下图为例的两个分段影像:       我们使用OpenCVSharp中的Stitcher类的Stitch方法,导入两张图像并拼接:  但结果

    2024年02月02日
    浏览(50)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包