文章快速检索  
  高级检索
顾及拓扑关系的立体影像直线特征可靠匹配算法
王竞雪1,2, 朱庆1, 王伟玺3     
1. 西南交通大学地球科学与环境工程学院, 四川 成都 611756;
2. 辽宁工程技术大学测绘与地理科学学院, 辽宁 阜新 123000;
3. 深圳大学智慧城市研究院, 广东 深圳 518060
摘要:针对单直线匹配过程中缺乏考虑邻近直线特征之间关系,纹理断裂处单一直线描述符的弱可靠性,提出了一种顾及拓扑关系的立体影像直线特征可靠匹配算法。该算法首先根据直线间距离、角度等基本拓扑关系对参考影像、搜索影像上提取的直线进行编组;然后将编组得到的直线组作为匹配基元,充分利用直线特征组内的拓扑关系,依次采用核线约束、单应矩阵约束、象限约束、不规则三角形区域灰度相关约束对其进行匹配;最后将同名直线组分裂为两对同名单直线、并对分裂后的结果进行整合、拟合、检核等后处理,得到“一对一”的同名直线。选取典型纹理特征的航空影像和近景影像进行参数分析及直线匹配试验,结果表明,本文算法能获取可靠的直线匹配结果。
关键词:直线匹配    直线组匹配    拓扑关系    核线约束    单应矩阵    
Reliable Line Matching Algorithm for Stereo Images with Topological Relationship
WANG Jingxue1,2, ZHU Qing1, WANG Weixi3     
1. Faculty of Geosciences and Environmental Engineering, Southwest Jiaotong University, Chengdu 611756, China;
2. School of Geomatics, Liaoning Technical University, Fuxin 123000, China;
3. Research Institute for Smart Cities, Shenzhen University, Shenzhen 518060, China
Foundation support: The National Natural Science Foundation of China(Nos.41631174;41101452); General Scientific Research Project Supported by Liaoning Provincial Education Department (No.LJYL010)
First author: WANG Jingxue(1981—),female,PhD,associate professor,majors in digital image matching.E-mail:xiaoxue1861@163.com
Corresponding author: ZHU Qing, E-mail:zhuq66@263.net
Abstract: Because of the lack of relationships between matching line and adjacent lines in the process of individual line matching, and the weak reliability of the individual line descriptor facing on discontinue texture, this paper presents a reliable line matching algorithm for stereo images with topological relationship. The algorithm firstly generates grouped line pairs from lines extracted from the reference image and searching image according to the basic topological relationships such as distance and angle between the lines. Then it takes the grouped line pairs as matching primitives, and matches these grouped line pairs by using epipolar constraint, homography constraint, quadrant constraint and gray correlation constraint of irregular triangle in order. And finally, it resolves the corresponding line pairs into two pairs of corresponding individual lines, and obtains one to one matching results after the post-processing of integrating, fitting, and checking. This paper adopts digital aerial images and close-range images with typical texture features to deal with the parameter analysis and line matching, and the experiment results demonstrate that the proposed algorithm in this paper can obtain reliable line matching results.
Key words: line matching     pair-wise line matching     topological relations     epipolar constraint     homography    

直线特征是人类视觉信息中最显著、最直观、最具有代表性的特征。城区建筑物轮廓边缘、屋顶和立面都蕴含着大量的直线特征,包含丰富的语言信息和几何意义,能够较好地表达建筑物结构特征,因此,将直线特征作为匹配基元并实现其可靠匹配是基于影像精细三维建模的基础[1-4]。目前,摄影测量和计算机视觉领域对直线匹配进行了大量的卓有成效的研究。现有的直线匹配方法可以分为单直线匹配和组直线匹配两类。

单直线匹配方法主要利用单一直线的几何属性信息(如斜率、长度、位置、重叠度、距离等)或邻域灰度信息的特征相似性,结合一定的几何约束对直线逐一进行匹配。文献[5]采用分层次直线符号法,利用直线段的长度、方向、位置和强度、对比度等特征进行直线的相似性匹配;结合核线约束,文献[6]利用直线邻域灰度相关实现直线匹配;结合同名三角网约束,文献[7-8]提出移动窗口的自适应直线相关方法,有效解决位于纹理断裂处的直线匹配问题;结合同名点约束,文献[9]通过窗口行向量平均化计算相似性;文献[10]利用物方平面约束,通过计算参考直线和候选直线间重叠距离和垂直距离,确定同名直线;文献[11]将待匹配线及与其共面的同名点作为整体进行描述及匹配;在无任何约束条件下,文献[12-13]提出了MSLD均值-标准差直线描述子(mean-standard deviation line descriptor, MSLD),该描述子通过统计像素支撑域内每个子区域4个方向的梯度向量构建描述子矩阵,提高描述子的可区分性及稳健性。目前单直线匹配研究取得了一定的进展[14-18],但受限于单一直线描述符,匹配过程中没有考虑邻近直线特征之间关系,匹配结果模糊性一直存在。

组直线匹配多利用满足一定条件的直线组作为匹配基元进行匹配。文献[19]将定义的干线对作为匹配基元,采用全局最优化思想完成异源图像的匹配及配准;文献[20]根据几何特性、灰度特性为直线分配可信度,以组直线为匹配单元,利用马氏距离计算组直线之间的特征相似性进行匹配;文献[21-23]提出了一种组直线匹配方法,首先参考影像上的直线组需满足距离、角度、半邻域灰度相似性3个条件,综合利用核线约束、几何属性约束、空间直方图、灰度相关等约束条件确定搜索影像上的同名直线组。研究结合Daisy描述符、匹配结果冗余及成对质量检核3个方面对匹配结果进行后处理,减少匹配结果不确定性;文献[24]将三维空间共面的直线组作为匹配基元,根据拓扑关系确定直线组,对直线组交点构建描述符,根据标准互相关函数确定同名直线组,该方法对纹理匮乏区域均能取得较好的匹配结果,但对于拍摄倾角变化较大的立体影像还需进一步完善。在现有组直线匹配方法的基础上,本文提出一种顾及拓扑关系的立体影像直线特征可靠匹配算法。该算法将满足一定拓扑关系的直线组作为匹配基元,充分利用直线组间几何信息、拓扑关系等多重约束对直线组进行匹配,并通过对同名直线组分裂、匹配结果后处理获取可靠的直线匹配结果。

1 原理

在已有LSD直线提取结果和渐进三角网约束点密集匹配结果的基础上,本文提出的顾及拓扑关系的立体影像直线特征可靠匹配算法分为特征直线组提取、直线组匹配和匹配结果后处理3个阶段。第1阶段,根据直线间距离、角度等拓扑关系对参考影像、搜索影像提取的直线分别进行特征编组,提取满足预设条件的特征直线组;第2阶段,充分利用直线组间拓扑关系,依次采用核线约束、单应矩阵约束、象限约束和三角形区域灰度相关约束对参考影像、搜索影像上提取的直线组进行匹配,得到两张影像上的同名直线组;第3阶段,将匹配得到的同名直线组分解为同名单直线,并对结果进行后处理,主要包括对结果中“一对多”“多对一”“多对多”匹配结果进行整合,剔除错误的“多”直线匹配结果,并对正确的“多”直线进行拟合,最终得到“一对一”的同名单直线。

1.1 拓扑关系约束的直线编组

采用LSD算法[25]对参考影像、搜索影像分别进行直线提取,结果表示为L={l1, l2, …, lI},R={r1, r2, …, rJ},IJ分别为参考影像、搜索影像上提取直线的数目。对参考影像、搜索影像提取的直线分别进行编组,得到特征直线组。对影像上任一直线段li,两端点分别为ab,长度为len,将其作为目标直线,直线{li+1, li+2, …, lI}为待编组的直线,具体编组过程如下:

(1) 以目标直线li为中心建立矩形窗口作为目标直线邻域窗口。如图 1所示,沿目标直线两端及垂直直线方向的两侧,分别扩展d个像素长度,得到大小为(2d+1)×(len+2d)的矩形窗口为目标直线的邻域窗口。

图 1 直线编组 Fig. 1 Finding intersection line pairs

(2) 依次判断待编组的直线与目标直线邻域窗口关系,如某一待编组直线li在目标直线邻域窗口内或与邻域窗口相交,且与目标直线li的夹角θ在一定的角度范围内,则该直线li与目标直线li构成一特征直线组,记为PL(li, li, CL),CL为特征直线组的交点坐标。

(3) 分别选择lili端点中与交点CL距离较近的点作为lili的起点,与点CL距离较远的点作为直线lili的终点。如图 1所示,be分别为两直线的起点,ac分别为两直线的终点。

对每张影像上提取的直线依次进行编组,分别得到参考影像和搜索影像上的特征直线组集合,记为PL={P1L(li1, li1, C1L), P2L(li2, li2, C2L), …, PML(liM, liM, CML)},PR={P1R(rj1, rj1, C1R), P2R(rj2, rj2, C2R), …, PNR(rjN, rjN, CNR)},其中,MN分别为参考影像和搜索影像上直线组的数目。图 2为参考影像、搜索影像上编组得到的直线组,其中圆点为直线组中两条直线的交点。

图 2 参考影像、搜索影像上编组得到的直线组 Fig. 2 Grouped line pairs from two images

1.2 顾及拓扑关系的直线组匹配

以参考参考、搜索影像上编组得到的直线组作为匹配基元进行匹配,匹配过程中依次采用核线约束、单应矩阵约束、象限约束、不规则三角形区域灰度相关约束。对参考影像上任一直线组,确定搜索影像上满足上述四种约束条件的特征直线组为其同名直线组。

1.2.1 核线约束

图 3(a)所示,参考影像上任一特征直线组PmL(lim, lim, CmL),计算其交点CmL在搜索影像上的核线。计算搜索影像上各特征直线组交点{C1R, C2R, …, CNR}到核线的距离,距离小于一定阈值Td的特征直线组为候选直线组。如图 3(b)所示,其中交点C1RC2RCnR对应的直线组为候选直线组。

图 3 多重约束示意图 Fig. 3 Diagram of multiple constraints

1.2.2 单应矩阵约束

单应矩阵表示两个平面之间的可逆齐次变换,同一平面上的点,在两张影像上成像的同名点坐标记为u =[xL, yL]Tv =[xR, yR]T,则u点可通过式(1)变换到v

(1)

式中,H为单应矩阵,通过两张影像上至少4对同名点进行求解。对于实际影像而言,一般都存在地形起伏或者景深变换,因此全局单应矩阵不适用于整张影像。本文对任一待匹配直线组分别选取其邻域范围内同名点计算局部单应矩阵,用于近似表示该直线组所在局部影像区域之间的对应关系,用于限定核线上的搜索范围。如图 3(a)所示,以参考直线组交点CmL为中心,以两直线limlim终点方向所构成扇形区域内包含的同名点计算局部单应矩阵H图 3中实心三角点为用于计算局部单应矩阵的同名点。为了确保用于单应矩阵计算所需同名点的数目,迭代扩大扇形半径直到同名点数目满足为止。

根据局部单应矩阵H和式(1),将参考影像上参考直线组交点CmL映射到搜索影像上,得到点,如图 3(b)所示;计算搜索影像上各候选直线组交点到点的距离,剔除距离大于阈值Th的候选直线组。

1.2.3 象限约束

根据局部单应矩阵H,将参考影像上参考直线组PmL(lim, lim, CmL)映射到搜索影像上,得到;如图 3(b)所示,在搜索影像上,对建立以为原点的局部坐标系,x轴和y轴分别与像素坐标系的x轴和y轴平行,记录中直线和直线终点所在象限分别为Q1Q2;对搜索影像上各候选直线组分别以其交点为原点建立局部坐标系,判断每个候选直线组中两直线终点所在象限是否为Q1Q2,是,保留该候选直线组,否则,舍弃该候选直线组。

1.2.4 不规则三角形区域灰度相关约束

对满上述3个约束的候选直线组,采用不规则三角形区域灰度相关约束确定最终同名直线组。

1.2.4.1 建立单直线对应关系

首先根据直线组中两直线所在象限建立参考直线组和候选直线组中单直线对应关系。固定参考影像上参考直线组PmL(lim, lim, CmL)中lim为起始边,lim为终止边。根据直线和直线所在象限Q1Q2确定候选直线组中与其对应的起始边和终止边,这里分两种情况:

Q1Q2时,即参考直线组中两条直线在两个不同的象限,候选直线组PnR(rin, rin, CnR)中,如果rin终点所在象限为Q1,则rin为起始边,rin为终止边,反之,如果rin终点所在象限为Q1,则rin为起始边,rin为终止边,更新候选直线组为PnR(rin, rin, CnR)。

Q1=Q2时,即参考直线组PmL(lim, lim, CmL)中两条直线在同一象限,根据向量叉积判断从limlim是顺时针还是逆时方向旋转,对应的,判断候选直线组PnR(rin, rin, CnR)中从rinrin是否与从limlim的旋转方向一致,如果不一致,更新候选直线组为PnR(rin, rin, CnR),否则,候选直线组不变。

1.2.4.2 对应三角形区域提取

图 3(a)所示,以参考直线组交点CmL、起始边终点a、终止边终点c构成三角形ΔaCmLc所包含的区域作为灰度相关区域,并采用扫描线算法确定三角形区域内所包含的像素点坐标及其灰度值,灰度值记为G={g1, …, gi, …, gNg},Ng为区域内包含的像素点数目。

图 3(b)所示,分别计算参考直线组起始边终点a、终止边终点c在搜索影像上的核线HaHc,分别计算HaHc与候选直线组中起始边和终止边的交点,记为a′和c′。搜索影像上以候选直线组交点、a′点和c′点构成三角形所包含的区域作为灰度相关区域。

利用参考影像、搜索影像灰度相关区域的3个顶点坐标,根据式(2)计算两灰度相关区域仿射变换参数

(2)

式中,a0a1a2b0b1b2表示仿射变换的6参数;xy表示参考影像上灰度相关区域顶点坐标,x′和y′表示搜索影像上对应的灰度相关区域对应顶点坐标。

根据仿射变换参数和参考影像上灰度相关区域ΔaCkLc内包含的像素点坐标计算搜索影像上灰度相关区域内对应像素点的坐标,并根据双线性插值方法计算得到对应像素点的灰度值,记为F={f1, …, fi, …, fNg}。

1.2.4.3 灰度相关计算

根据相关系数式(3)计算参考影像、搜索影像灰度相关区域的灰度相关系数ρ

(3)

式中,gifi分别为参考影像、搜索影像灰度相关区域内对应像素的灰度值;gf分别为对应的灰度相关区域内所有像素灰度值的平均值。选取相关系数最大且大于一定阈值Tρ的灰度相关区域对应的候选直线组为PmL(lim, lim, CmL)的同名特征直线组。

1.3 匹配结果后处理

根据直线组中单直线应关系将匹配得到的每对同名直线组分裂为两对同名单直线,由于直线编组过程中同一直线可能出现在多个直线组中,参与多次匹配,会得到多次匹配结果。同时由于直线提取结果断裂和错误匹配,匹配结果中存在“一对多”“多对一”“多对多”匹配情况,因此需要对匹配结果进行后处理。在确保匹配结果中对应关系唯一性的基础上,对结果中非“一对一”的“多”直线进行检核,剔除错误的“多”直线匹配结果,并对正确的“多”直线进行拟合,最终得到“一对一”的同名单直线。

首先,根据文献[10]中匹配结果整合算法对分裂得到的同名单直线依次实现从左到右,再从右到左的“一对多”匹配结果整合。然后,对整合后同名直线对应关系中“一对多”“多对一”“多对多”结果中的“多”直线进行拟合,如果拟合直线斜率与原有的多条直线斜率均一致且“多”直线间重叠度为0,则用拟合直线代替原有的“多”直线,否则将该组匹配结果删除。最终得到参考影像上与搜索影像上“一对一”的同名直线。

2 试验与分析

为了验证本文算法的有效性,共采用一组航空影像、两组近景影像用于直线匹配试验。第1组为图 4(a)所示的DMC数字航空影像的子影像对,影像集中在人工建筑物密集区域,提取的直线特征多为建筑物边缘,影像大小分别为499× 463像素、487×450像素;第2组为图 4(b)所示存在尺度变化的近景影像对,两张影像之间尺度比约为1:1.5,影像大小均为600×800像素;第3组为图 4(c)所示存在旋转的近景影像对,两张影像旋转角度约为40°,影像上地物细节信息较为丰富。试验过程中,首先对图 4所示3组影像实现渐进三角网约束点密集匹配及LSD直线提取,然后在此基础上进行直线匹配。匹配过程中参数阈值设置如下:Td=5,Th=20,Tρ=0.7。

图 4 影像数据 Fig. 4 Test images

2.1 直线编组参数选取

本文采用直线组作为匹配基元进行直线匹配,而直线编组结果对直线匹配结果产生直接的影响。本文编组过程中共包含直线邻域扩展长度d和直线夹角θ角度范围两个参数,因此本文采用第1组航空影像数据,在dθ取不同值的情况下进行直线编组及匹配试验,分析编组过程中不同参数值对编组的影响,进而研究其对匹配结果的影响,通过对比分析确定编组参数用于后续匹配。

在匹配参数一致的条件下,本文对θ∈[45°, 135°]、θ∈[30°, 150°]、θ∈[15°, 165°]3个取值范围及d分别取5、10、15、20、25值的情况下进行匹配试验,结果分别如表 1表 2表 3所示。其中,第1列为d的取值;第2、3列为参考影像、搜索影像上编组得到的直线组的数目;第4列为匹配得到的同名直线组的数目;第5列为同名直线组分裂后得到同名直线的数目。但是由于匹配结果中存在“一对多”“多对一”“多对多”对应关系,因此结果中包含参考影像、搜索影像上直线数目不相同,第6、7列为匹配结果中分别包含的参考影像、搜索影像上提取直线的数目;第8列为经过对“一对多”“多对一”“多对多”结果整合、拟合、检核后得到的最终的同名直线数;第9列为正确匹配的数目,第10列为匹配的正确率。

表 1 θ∈[45°, 135°]的直线匹配结果 Tab. 1 Line matching results by different parameter d and θ∈[45°, 135°]
d 上编组直线组数目 同名直线组数目 分裂后同名直线数目 结果中包含直线的数目 整合后同名直线数目 正确匹配数目 匹配正确率/(%)
参考影像 搜索影像 参考影像 搜索影像
5 129 97 52 91 89 85 80 79 98.7
10 258 253 132 171 165 152 131 128 97.7
15 420 399 233 253 237 216 189 184 97.3
20 584 561 333 306 277 250 204 201 98.5
25 790 767 475 383 313 288 175 171 97.7

表 2 θ∈[30°, 150°]的直线匹配结 Tab. 2 Line matching results by different parameter d and θ∈[30°, 150°]
d 上编组直线组数目 同名直线组数目 分裂后同名直线数目 结果中包含直线的数目 整合后同名直线数目 正确匹配数目 匹配正确率/(%)
参考影像 搜索影像 参考影像 搜索影像
5 135 100 54 95 93 89 84 83 98.8
10 271 260 136 175 169 156 135 132 97.8
15 440 410 239 258 241 220 191 185 96.9
20 614 576 339 310 280 254 206 203 98.5
25 825 792 488 389 316 292 174 169 97.1

表 3 θ∈[15°, 165°]的直线匹配结果 Tab. 3 Line matching results by different parameter d and θ∈[15°, 165°]
d 上编组直线组数目 同名直线组数目 分裂后同名直线数目 结果中包含直线的数目 整合后同名直线数目 正确匹配数目 匹配正确率/(%)
参考影像 搜索影像 参考影像 搜索影像
5 151 114 55 97 95 91 86 85 98.8
10 308 294 144 185 178 162 137 131 95.6
15 489 460 252 269 249 224 190 183 96.3
20 684 649 351 321 285 256 204 200 98.0
25 916 885 510 401 321 294 167 162 97.0

图 4(a)两张影像上提取直线的数目分别为398条、406条。一方面,由于一条直线可能存在多个组中,因此编组得到的直线组数会多于直线提取的数目;另一方面,由于夹角和邻域范围的设定,会存在一些直线未被编组,因此编组得到直线组数也会少于直线提取的数目。

表 1表 2表 3中可以看出随着θ角度范围增大、d值增大,参考影像、搜索影像上编组得到的直线组数目越多,匹配得到的同名直线组数目越多,对应分裂后同名单直线数目越多,且包含参考影像、搜索影像上直线数目越多。但是随着结果中“一对多”“多对一”“多对多”对应关系增多,匹配后处理后得到的同名直线的数目并没有因此而增多。这是由于本文在对“多”结果的拟合和检核的过程中采用了“一刀切”的方法,导致部分正确的匹配结果也被删除掉,影响最终同名直线的数目。尤其在误匹配较多的情况下,同名直线数目会大幅度的减小,例如,当d=25时,参考影像、搜索影像上编组后得到直线组的重复率较小,直线组匹配正确率不高,错误匹配较多,最终同名直线数目较小。后续将针对后处理中合并和检核过程进行针对性的研究。此外,针对图 4(a)所示影像,匹配结果中误匹配多发生在人字形的屋顶处,主要是由于不同视角产生的直线投影不同导致直线提取结果的不同而引起的错误匹配。

2.2 本文算法匹配结果及不同方法对比分析

通过2.1节不同参数对比分析试验,本文选取参数θ∈[30°, 150°],d=20对图 4所示影像进行直线匹配试验。为了与本文算法进行更加客观的对比分析,在初始同名点及直线提取结果相同的条件下,本文实现了文献[13]中基于MSLD描述子的直线匹配方法,实现过程中在原文献描述子基础上,增加了利用同名三角约束匹配候选直线,对其结果也采用本文的后处理算法进行处理。文献[13]根据最小欧氏距离与次最小欧氏距离的比值小于阈值TN确定同名直线,TN=0.8。两种方法直线匹配结果显示如图 5-图 7所示,人工目视判读直线匹配结果是否正确,其中红色线为正确匹配直线,蓝色线为错误匹配直线。两种算法直线匹配结果统计如表 4所示,其中第1列为试验影像,第2、3列为参考影像、搜索影像上提取直线的数目,第4列为采用的算法,第5列为匹配得到同名直线的数目,后5列表示的含义与表 1所述含义一致。从中可以看出对于图 4(a)航空影像和图 4(b)存在尺度变化影像,本文算法取得较好的结果。航空影像上特征线多位于建筑物边缘,而不同影像上同一建筑物两侧纹理信息的一致性较弱,MSLD构建描述子过程中同时考虑直线两侧信息,这对不同影像上同一建筑物边缘相似性的计算会产生影响,而本文算法利用直线组构建相关区域,对于位于同一屋顶面上的直线组不会受纹理断裂的影响;对于尺度变化影像,主要是因为MSLD描述子本身没有考虑尺度变化,而本文算法在确定三角形灰度相关区域时用到仿射变换,因此对于小尺度变化的影像具有一定的适应性;但对于图 4(c)所示纹理丰富,细节较多的影像,MSLD的结果要好于本文算法。这是由于MSLD描述子利用梯度特征,其稳健性强于单纯基于区域的灰度相关,尤其对于细节信息丰富的影像,MSLD能产生较好的匹配效果,产生较少的错误匹配。而本文单纯基于区域灰度相关的匹配结果可靠性较弱,初始匹配会产生较多的错误匹配,后处理后删除掉大部分错误匹配,但是同名直线的数目也大幅度减小。后续研究将会考虑用基于梯度的描述子代替传统的灰度相关。

图 5 图 4(a)所示影像不同算法的直线匹配结果 Fig. 5 Line matching results by different algorithms for figure 4(a)

图 6 图 4(b)所示影像不同算法的直线匹配结果 Fig. 6 Line matching results by different algorithms for figure 4(b)

图 7 图 4(c)所示影像不同算法的直线匹配结果 Fig. 7 Line matching results by different algorithms for figure 4(c)

表 4 不同算法直线匹配结果 Tab. 4 Line matching results by different algorithms
立体影像对 上编组直线组数目 算法 同名直线数目 结果中包含直线的数目 整合后同名直线数目 正确匹配数目 匹配正确率/(%)
参考影像 搜索影像 参考影像 搜索影像
图 4(a) 398 406 MSLD 216 216 198 188 183 97.3
本文算法 310 280 254 206 203 98.5
图 4(b) 689 884 MSLD 123 123 122 122 119 97.5
本文算法 347 287 308 200 197 98.5
图 4(c) 1170 844 MSLD 423 423 322 277 274 98.9
本文算法 408 353 283 170 166 97.6

此外,匹配过程中涉及3个参数阈值TdThTρ的设定。其中,TdTh主要用于核线约束和单应矩阵约束匹配候选,在本文算法匹配过程中这两种约束较为宽泛,在一定范围内受阈值影响较小,后续象限约束较为严格。相关系数阈值采用经验化阈值,取为0.7。

3 结论

本文探索了一种顾及拓扑关系的立体影像直线特征可靠匹配算法。该算法主要特点表现在以下几个方面:①在已有单应矩阵的基础上,提出了象限约束,使得特征编组具有任意性,不用考虑任何对应关系,简化了编组过程;②对位于同一面片上的直线组,能有效避免纹理断裂对匹配产生的影响,提高匹配的可靠性,同时,三角形灰度区域建立过程中引入的仿射变换使得该方法适用于存在一定尺度变化的影像;③对“一对多”“多对一”“多对多”的匹配结果进行整合,并对整合结果中的“多”直线进行检核,对正确的“多”直线进行拟合,剔除错误的“多”直线匹配结果,得到“一对一”的同名直线。此算法一方面解决了直线提取过程中产生的断裂,另一方面提高了匹配的正确率。

下一步研究工作:①目前算法没有充分利用匹配结果中的冗余信息,同时匹配结果后处理过程中,采用“一刀切”政策,部分正确的匹配结果被删除,使得匹配直线数目大大减少,后续将重点研究充分利用冗余信息提高后处理的高效性及有效性;②现有研究是基础性和探索性的,后续将现有算法中基于灰度的描述符替换成基于梯度信息的描述符,提高描述符的稳定性和可靠性;③进一步将考虑未被编组的直线,使其加入匹配。


参考文献
[1] ZHANG Li.Automatic Digital Surface Model(DSM)Generation from Linear Array Images[D].Zurich:Institute of Geodesy and Photogrammetry, 2005. http://dx.doi.org/10.3929/ethz-a-005055636
[2] 江万寿. 航空影像多视匹配与规则建筑物自动提取方法研究[D]. 武汉: 武汉大学, 2004.
JIANG Wanshou.Multiple Aerial Image Matching and Automatic Building Detection[D].Wuhan:Wuhan University, 2004. http://cdmd.cnki.com.cn/Article/CDMD-10486-2006031556.htm
[3] 文贡坚. 一种基于特征编组的直线立体匹配全局算法[J]. 软件学报, 2006, 17(12): 2471–2484.
WEN Gongjian. A Global Algorithm for Straight Line Stereo Matching Based on Feature Grouping[J]. Journal of Software, 2006, 17(12): 2471–2484.
[4] ZHANG Lilian, KOCH R. An Efficient and Robust Line Segment Matching Approach Based on LBD Descriptor and Pairwise Geometric Consistency[J]. Journal of Visual Communication and Image Representation, 2013, 24(7): 794–805. DOI:10.1016/j.jvcir.2013.05.006
[5] 席学强, 王润生. 基于直线特征的图像-模型匹配算法[J]. 国防科技大学学报, 2000, 22(6): 70–74.
XI Xueqiang, WANG Runsheng. An Algorithm of Image-model Matching Based on Straight Line Features[J]. Journal of National University of Defense Technology, 2000, 22(6): 70–74.
[6] SCHMID C, ZISSERMAN A.Automatic Line Matching Across Views[C]//Proceedings of the 2013 IEEE Conference on Computer Vision and Pattern Recognition.San Juan, Puerto Rico:IEEE, 1997:666-671.
[7] 张云生, 朱庆, 吴波, 等. 一种基于三角网约束的立体影像线特征多级匹配方法[J]. 武汉大学学报(信息科学版), 2013, 38(5): 522–527.
ZHANG Yunsheng, ZHU Qing, WU Bo, et al. A Hierarchical Stereo Line Matching Method Based on a Triangle Constraint[J]. Geomatics and Information Science of Wuhan University, 2013, 38(5): 522–527.
[8] WU Bo, ZHANG Yunsheng, ZHU Qing. Integrated Point and Edge Matching on Poor Textural Images Constrained by Self-adaptive Triangulations[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2012, 68(1): 40–55.
[9] 梁艳, 盛业华, 张卡, 等. 利用局部仿射不变及核线约束的近景影像直线特征匹配[J]. 武汉大学学报(信息科学版), 2014, 39(2): 229–233.
LIANG Yan, SHENG Yehua, ZHANG Ka, et al. Linear Feature Matching Method Based on Local Affine Invariant and Epipolar Constraint for Close-range Images[J]. Geomatics and Information Science of Wuhan University, 2014, 39(2): 229–233.
[10] 王竞雪, 宋伟东, 王伟玺. 同名点及高程平面约束的航空影像直线匹配算法[J]. 测绘学报, 2016, 45(1): 87–95.
WANG Jingxue, SONG Weidong, WANG Weixi. Line Matching Algorithm for Aerial Image Based on Corresponding Points and Z-plane Constraints[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(1): 87–95. DOI:10.11947/j.AGCS.2016.20140527
[11] FAN Bin, WU Fuchao, HU Zhanyi.Line Matching Leveraged By Point Correspondences[C]//Proceedings of 2010 IEEE Conference on Computer Vision and Pattern Recognition.San Francisco:IEEE, 2010:390-397. http://ieeexplore.ieee.org/xpls/icp.jsp?arnumber=5540186
[12] 王志衡, 吴福朝. 均值-标准差描述子与直线匹配[J]. 模式识别与人工智能, 2009, 22(1): 32–39.
WANG Zhiheng, WU Fuchao. Mean-standard Deviation Descriptor and Line Matching[J]. Pattern Recognition and Artificial Intelligence, 2009, 22(1): 32–39.
[13] WANG Zhiheng, WU Fuchao, HU Zhanyi. MSLD:A Robust Descriptor for Line Matching[J]. Pattern Recognition, 2009, 42(5): 941–953. DOI:10.1016/j.patcog.2008.08.035
[14] AL-SHAHRI M, YILMAZ A. Line Matching in Wide-Baseline Stereo:A Top-down Approach[J]. IEEE Transactions on Image Processing, 2014, 23(9): 4199–4210.
[15] WANG Jingxue, WANG Weixi, LI Xiaoming, et al.Line Matching Algorithm for Aerial Image Combining Image and Object Space Similarity Constraints[C]//International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences:XLI-B3.Prague, Czech Republic:ISPRS, 2016:783-788.
[16] 缪君, 储珺, 张桂梅. 一种仿射不变的直线描述子与直线匹配[J]. 电子学报, 2015, 43(12): 2505–2512.
MIAO Jun, CHU Jun, ZHANG Guimei. An Affine Invariant Line Descriptor and Line Matching[J]. Acta Electronica Sinica, 2015, 43(12): 2505–2512. DOI:10.3969/j.issn.0372-2112.2015.12.024
[17] 宋伟东, 朱红, 王竞雪, 等. 多重约束下的近景影像线特征匹配方法[J]. 中国图象图形学报, 2016, 21(6): 764–770.
SONG Weidong, ZHU Hong, WANG Jingxue, et al. Line Feature Matching Method based on Multiple Constraints for Close-range Images[J]. Journal of Image and Graphics, 2016, 21(6): 764–770. DOI:10.11834/jig.20160609
[18] 张跃强, 苏昂, 朱遵尚, 等. 结合局部描述和拓扑约束的直线段匹配算法[J]. 国防科技大学学报, 2014, 36(6): 25–30.
ZHANG Yueqiang, SU Ang, ZHU Zunshang, et al. A Line Matching Algorithm Combining Local Description with Topological Constraint[J]. Journal of National University of Defense Technology, 2014, 36(6): 25–30. DOI:10.11887/j.cn.201406005
[19] 高峰, 文贡坚, 吕金建. 基于干线对的红外与可见光最优图像配准算法[J]. 计算机学报, 2007, 30(6): 1014–1021.
GAO Feng, WEN Gongjian, LV Jinjian. An Optimal Algorithm for IR/Visual Image Registration Based on Main-line-pairs[J]. Chinese Journal of Computers, 2007, 30(6): 1014–1021.
[20] 张浩, 才辉, 张光新, 等. 一种新的基于边缘拟合的图像配准方法[J]. 光电子·激光, 2009, 20(1): 103–107.
ZHANG Hao, CAI Hui, ZHANG Guangxin, et al. A Novel Image Registration Method Based on Edge Fitting[J]. Journal of Optoelectronics·Laser, 2009, 20(1): 103–107.
[21] OK A O, WEGNER J D, HEIPKE C, et al. Matching of Straight Line Segments from Aerial Stereo Images of Urban Areas[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2012, 74(1): 133–152.
[22] OK A Ö, WEGNER J D, HEIPKE C, et al.A Stereo Line Matching Technique for Aerial Images Based on a Pair-wise Relation Approach[C]//Proceedings of the Istanbul Workshop 2010 on Modeling of Optical Airborne and Space borne Sensors:XXXVⅢ-1/W17.Istanbul, Turkey:ISPRS, 2010.
[23] OK A Ö, WEGNER J D, HEIPKE C, et al.A New Straight Line Reconstruction Methodology from Multi-spectral Stereo Aerial Images[C]//PAPARODITIS N, PIERROT-DESEILLIGNY M, MALLET C, et al.IAPRS, Vol.XXXVⅢ.Saint-Mandé, France:ISPRS, 2010:25-30.
[24] KIM H, LEE S. Simultaneous Line Matching And Epipolar Geometry Estimation Based on the Intersection Context of Coplanar Line Pairs[J]. Pattern Recognition Letters, 2012, 33(10): 1349–1363. DOI:10.1016/j.patrec.2012.03.014
[25] VON DIOI R G, JAKUBOWICZ J, MOREL J M, et al. LSD:A Fast Line Segment Detector with a False Detection Control[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2010, 32(4): 722–732. DOI:10.1109/TPAMI.2008.300
http://dx.doi.org/10.11947/j.AGCS.2017.20170162
中国科学技术协会主管、中国测绘地理信息学会主办。
0

文章信息

王竞雪,朱庆,王伟玺
WANG Jingxue, ZHU Qing, WANG Weixi
顾及拓扑关系的立体影像直线特征可靠匹配算法
Reliable Line Matching Algorithm for Stereo Images with Topological Relationship
测绘学报,2017,46(11):1850-1858
Acta Geodaetica et Cartographica Sinica, 2017, 46(11): 1850-1858
http://dx.doi.org/10.11947/j.AGCS.2017.20170162

文章历史

收稿日期:2017-04-06
修回日期:2017-09-19

相关文章

工作空间