一种抗视角变化的近景影像区域特征匹配方法组成比例

本发明专利技术公开了一种抗视角变化的近景影像区域特征匹配方法，包括如下步骤：提取最大稳定极值区域MSER；确定局部二值模式LBP邻域形状、大小和旋转方向；提取五种LBP特征并计算联合统计直方图；比较各个特征区域联合统计直方图相似度完成粗匹配；基于RANSAC与基本矩阵估计剔除粗差点；对匹配区域进行质心偏移修正完成精匹配。本发明专利技术利用一种改进的局部二值模式LBP计算方法，直接对不规则MSER区域中的每一个像素进行编码，避免了传统方法中拟合SIFT匹配区域所带来的拟合误差和边界误差，其本身具有计算复杂度低、无需训练学习等优点，其对于线性的亮度变化保持稳定，对非线性的光照变化存在容错能力，提高了区域特征匹配的可靠性，降低了计算复杂度。

【技术实现步骤摘要】
一种抗视角变化的近景影像区域特征匹配方法
本专利技术涉及一种抗视角变换的近景影像区域特征匹配方法，本专利技术所属
为属于摄影测量领域。
技术介绍
影像匹配即找寻影像之间的同名点是影像定向和模型连接的前提和关键。在一些狭窄的隧洞中，沿着主光轴采集序列影像，相邻影像是沿摄影机光轴移动所构成的平行序列影像，得到的是一种尺度渐进变化的影像序列，且所成图像间会因为相机所处位置的调整而呈现较大尺度和视角变换的情况，其特点包括：基线不是大致平行于影像而是大致垂直于影像，影像的重叠范围是沿着影像的辐射线方向，影像间尺度变化较大、影像重叠区域的分辨率差别较大、同名点视角差异大。其对于一般立体像对而言具有很大的差别，不能采用针对传统立体影像的匹配方法。最大稳定极值区域（MSER）是一种十分具有影响力的区域检测算子，图像中的最大稳定极值区域指的是一个由彼此相连的一个像元集合，该集合中的每个像元的像素值全部大于或者全部小于包围这个区域的边界的每个像元的像素值，该区域具有严格意义上的仿射不变性。对于视角变换，MSER在处理结构型图像和纹理型图像都有十分明显的优越性。传统的最大稳定极值区域匹配方法是将特征区域归一化后用SIFT算子描述，但是SIFT描述MSER特征区域有其不足，第一，在构造描述算法前需要把不规则区域进行椭圆形状拟合，该过程会产生拟合误差和降低匹配性能，特别是对于大视角变换图像更为明显；第二，其基于几何固定位置（圆形邻域）划分子区域，这样会产生边界误差；第三，其计算复杂度高，导致大大提高了图像特征提取的时间复杂度。
技术实现思路
为了解决上述现有技术中减少拟合匹配区域所带来的拟合误差和边界误差的影响，本专利技术提供给了一种对于线性的亮度变化保持稳定，对非线性的光照变化存在一定的容错能力，提高了区域特征匹配的可靠性，并降低计算复杂度的抗视角变化的近景影像区域特征匹配方法。本专利技术的目的是这样实现的：一种抗视角变换的近景影像区域特征匹配方法，包括如下步骤：第一步，提取匹配影像内的最大稳定极值区域MSER；第二步，采用直接最小二乘法对提取出的形状不规则的最大稳定极值区域拟合椭圆，并以此确定局部二值模式LBP邻域大小和编码起始方向；第三步，提取五种LBP特征；第四步，计算各特征区域的五种LBP特征的联合统计直方图；第五步，利用相关系数比较各个特征区域的联合统计直方图相似度，从而进行特征区域的初步匹配；第六步，利用RANSAC与八点法相结合的策略估计基本矩阵，剔除粗差点；第七步，利用最小二乘匹配法对已配对特征区域的椭圆中心进行微调，得到更为精确的匹配点对；上述方法工作原理包括：最大稳定极值区域原理、局部二值模式原理和计算机视觉原理；其工作过程如下：①先分别提取要进行特征匹配的两张影像上的最大稳定极值区域；②对两张影像上的所有最大稳定极值区域拟合椭圆，并确定所有最大稳定极值区域局部二值模式特征的椭圆邻域大小和起始编码方向；③计算两张影像上所有最大稳定极值区的五种局部二值模式特征；④计算两张影像上所有最大稳定极值区的局部二值模式特征直方图；计算两张影像上特征区域的局部二值模式特征直方图相似度，完成粗匹配估计基本矩阵，剃除粗差点；利用最小二乘匹配法对粗匹配结果进行微调，得到更为精确的匹配结果；第二步中所述LBP邻域为椭圆形，其大小为：若某一特征区域的拟合椭圆长半轴为A、短半轴为B，则LBP椭圆形邻域的长半轴为a、b，他们之间的关系为：a=A/10,b=B/10；第二步中所述主方向：采用椭圆右侧长半轴方向为主方向，即以此方向为LBP初始编码方向，顺时针进行编码；第三步中所述五种LBP特征，计算公式为别为：其中，gc为邻域中心点的像素灰度值，c1为图像中所有像素灰度的均值；其中，P为邻域采样点个数，L为不同半径邻域的个数，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值，gc为邻域中心点的像素灰度值；其中，P为邻域采样点个数，L为不同半径邻域的个数，，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值，gc为邻域中心点的像素灰度值，而是在第l个半径邻域中的的均值；其中，P为邻域采样点个数，L为不同半径邻域的个数，gp+1,Rl为第l重半径邻域中第p+1个采样点的像素灰度值，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值；其中，P为邻域采样点个数，L为不同半径邻域的个数，gp+1,Rl为第l重半径邻域中第p+1个采样点的像素灰度值，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值，为第l个半径领域中所有的均值；第五步中所述相关系数的阈值可根据实际情况设置，推荐阈值为0.95，高于这一阈值即认为是一对相匹配特征区域；积极有益效果：与现有技术相比，本专利技术所提供的抗视角变换的近景影像区域特征匹配方法具有两方面的优点：采取了不同半径LBP邻域参与计算，提高了区域特征匹配的可靠性；对椭圆拟合区域不做归一化处理，降低了计算复杂度，节省了计算时间。附图说明图1为本专利技术实施例的流程图。具体实施方式下面结合附图及具体实施例，对本专利技术进一步的说明：实施例根据图1所示的方法流程，以“地质勘探洞近景摄影测量影像匹配”为应用实例，对本专利技术进一步阐明：沿着洞轴方向一前一后拍了两张影像，属于平行摄影测量方式，即摄影基线与相机主光轴一致，对这两张影像进行特征区域的匹配；过程如下：（1）分别提取两张影像上的最大稳定极值区域（MSER）；（2）假设一张影像上有n个最大稳定极值区域，对所有提取出的形状不规则的最大稳定极值区域进行椭圆拟合，得到n个椭圆，其长短半轴分别为A1,B1,A2,B2,...,An,Bn,其长半轴与屏幕坐标系y轴正方向的夹角分别为α1,α2,...,αn；则各最大稳定极值区域的椭圆形LBP邻域也可以确定：其长半轴分别为a1=A1/10,a2=A2/10,...,an=An/10;短半轴分别为b1=B1/10,b2=B2/10,...,bn=Bn/10；椭圆形LBP邻域的旋转角度也可以确定，分别为α1,α2,...,αn；若选择多重LBP邻域来计算，则第2重邻域长短半轴分别为第1重邻域的2倍，第3重邻域长短半轴分别为第一重邻域的3倍，以此类推；（3）按以下公式计算n个最大稳定极值区域各自的5种局部二值模式（LBP）特征：其中，gc为邻域中心点的像素灰度值，c1为图像中所有像素灰度的均值；其中，P为邻域采样点个数，L为不同半径邻域的个数，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值，gc为邻域中心点的像素灰度值；其中，P为邻域采样点个数，L为不同半径邻域的个数，，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值，gc为邻域中心点的像素灰度值，而是在第l个半径邻域中的的均值；其中，P为邻域采样点个数，L为不同半径邻域的个数，gp+1,Rl为第l重半径邻域中第p+1个采样点的像素灰度值，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值；其中，P为邻域采样点个数，L为不同半径邻域的个数，gp+1,Rl为第l重半径邻域中第p+1个采样点的像素灰度值，gp,Rl为第l重半径邻域中第p个采样点的像素灰度值，为第l个半径领域中所有的均值；（4）计算两张影像上所有最大稳定极值区的5种局部二值模式特征直方图，将5种直方图归一化并连接在一起；本文档来自技高网...

【技术保护点】
1.一种抗视角变换的近景影像区域特征匹配方法，其特征在于，包括如下步骤：第一步，提取匹配影像内的最大稳定极值区域MSER；第二步，采用直接最小二乘法对提取出的形状不规则的最大稳定极值区域拟合椭圆，并以此确定局部二值模式LBP邻域大小和编码起始方向；第三步，提取五种LBP特征；第四步，计算各特征区域的五种LBP特征的联合统计直方图；第五步，利用相关系数比较各个特征区域的联合统计直方图相似度，从而进行特征区域的初步匹配；第六步，利用RANSAC与八点法相结合的策略估计基本矩阵，剔除粗差点；第七步，利用最小二乘匹配法对已配对特征区域的椭圆中心进行微调，得到更为精确的匹配点对。

【技术特征摘要】
1.一种抗视角变换的近景影像区域特征匹配方法，其特征在于，包括如下步骤：第一步，提取匹配影像内的最大稳定极值区域MSER；第二步，采用直接最小二乘法对提取出的形状不规则的最大稳定极值区域拟合椭圆，并以此确定局部二值模式LBP邻域大小和编码起始方向；第三步，提取五种LBP特征；第四步，计算各特征区域的五种LBP特征的联合统计直方图；第五步，利用相关系数比较各个特征区域的联合统计直方图相似度，从而进行特征区域的初步匹配；第六步，利用RANSAC与八点法相结合的策略估计基本矩阵，剔除粗差点；第七步，利用最小二乘匹配法对已配对特征区域的椭圆中心进行微调，得到更为精确的匹配点对。2.根据权利要求1所述的一种抗视角变换的近景影像区域特征匹配方法，其特征在于：第二步中所述LBP邻域为椭圆形，其大小为：若某一特征区域的拟合椭圆长半轴为A、短半轴为B，则LBP椭圆形邻域的长半轴为a、b，他们之间的关系为：a=A/10,b=B/10。3.根据权利要求1所述的一种抗视角变换的近景影像区域特征匹配方法，其特征在于：第二步中所述主方向：采用椭圆右侧长半轴方向为主方向，即以此方向为LBP初始编码方向，顺时针进行编码。4...

【专利技术属性】
技术研发人员：朱邦彦，储征伟，毛羽丰，刘文伍，段伟，
申请(专利权)人：南京市测绘勘察研究院股份有限公司，
类型：发明
国别省市：江苏,32