本发明专利技术提供一种基于平面透视投影的倾斜影像匹配方法，属于倾斜影像匹配技术领域。本发明专利技术包括以下步骤：首先利用倾斜影像的外方位元素和内方位元素，将倾斜影像通过透视投影到测区内的平均高程面上，获得近似的正射影像；然后利用SIFT算法对近似正射影像进行匹配，并对匹配结果进行基于基础矩阵和单应矩阵双重约束的粗差剔除；最后将获得的匹配点通过逆透视变换至原始影像上，获得原始影像的匹配结果。本发明专利技术相对于其他的匹配方法，能够有效地提高倾斜影像匹配速度、准确率以及匹配点对数。
An oblique image matching method based on plane perspective projection
【技术实现步骤摘要】
一种基于平面透视投影的倾斜影像匹配方法
本专利技术涉及倾斜影像匹配
，具体地说是一种基于平面透视投影的倾斜影像匹配方法。
技术介绍
相对于传统的航空摄影测量，倾斜摄影测量能够实现建筑物的正直影像和其立面纹理的倾斜影像相结合，大大减少了建筑物侧面纹理信息的采集工作，提高了三维建模效率。然而由于倾斜摄影独特的拍摄方式，使得倾斜影像存在几何变形大，色差大，遮挡严重，数据海量等特点，这给倾斜影像匹配带来了巨大的困难。目前，针对倾斜影像匹配并不存在一种普遍适用的匹配方法。Harris-Affine、Hessian–Affine、MSER、IBR、EBR、Salient等是常用的仿射不变局部特征检测子，广泛用于影像匹配算法中，其中MSER、Harris-Affine&Hessian–Affine效果最好，但它们都不具备完全仿射不变性，对倾斜影像匹配效果不佳；ASIFT是一种具有完全仿射不变性的匹配算法，对于倾斜影像匹配效果较好，但它采取近似穷举的方法获取影像的最佳经度和纬度，效率很低；另外一些算法沿用了ASIFT思想进行基于影像仿射变换来进行匹配，但各有优缺点。
技术实现思路
本专利技术的技术任务是针对现有技术的不足，提供一种基于平面透视投影的倾斜影像匹配方法。该方法在某种程度上对影像进行了透视变换，较之ASIFT算法对影像进行仿射变换，经过透视变换的影像之间差异会更小，更容易匹配。本专利技术解决其技术问题所采用的技术方案是：本专利技术的一种基于平面透视投影的倾斜影像匹配方法，包括以下步骤：步骤1，利用原始倾斜影像S的外方位元素和内方位元素，将原始倾斜影像S通过透视投影到测区内的平均高程面Hm上，获得近似的正射影像，作为平面透视投影之后的目标影像D；步骤2，获取目标影像对之间的重叠区域，作为待匹配区域：对于两张待匹配的原始倾斜影像，通过步骤1，分别获取它们的目标影像对及其范围，根据多边形求交运算获取目标影像对之间的重叠区域四边形abcd，将重叠区域abcd作为待匹配区域；步骤3，采用SIFT算法对目标影像对进行匹配，使用左右一致性检验、唯一性约束等对粗差点进行剔除，获得粗匹配结果；通过单应矩阵和基础矩阵双重约束对粗匹配结果进行筛选，获得最终的精匹配结果；步骤4，将步骤3获得的目标影像对上的匹配点对，利用共线方程反投影至原始倾斜影像对上，作为最终匹配结果。作为本专利技术进一步改进的，所述步骤1具体实现步骤如下：步骤1.1，将原始倾斜影像S的四个角点(xi,yi)(其中i＝0,1,2,3)，通过透视投影投影至目标影像D上，对应目标影像D上的四个角点(x′i,y′i)(其中i＝0,1,2,3)；步骤1.2，利用目标影像四个角点(x′i,y′i)，确定目标影像D的边界、宽W′和高H′，建立W′×H′的空白目标影像D；步骤1.3，灰度填充法填充目标影像D。作为本专利技术进一步改进的，步骤1.1中所述目标影像D所在的高程面固定，为整个测区内的平均高程Hm；原始倾斜影像S上任意一像点(x,y)，通过公式(1)经过透视投影至目标影像D上的物方点坐标为(X，Y，Hm)，(X，Y)为目标影像D的像点坐标(x′i,y′i)；其中，f为焦距，a1、a2、a3、b1、b2、b3、c1、c2、c3为旋转矩阵的元素，Xs、Ys、Zs为摄影中心坐标。作为本专利技术进一步改进的，步骤1.2中，目标影像D的边界(Xmax,Xmin,Ymax,Ymin)、宽W′和高H′为：作为本专利技术进一步改进的，步骤1.2中，使用比例因子s约束目标影像D的大小，使目标影像D和原始影像S之间大小一致；其中W为原始倾斜影像S的宽度。作为本专利技术进一步改进的，步骤1.3中，每一个空白目标影像D的像素点(x′,y)′，由公式(5)获得，它对应着物方点坐标(X,Y,Z)，利用共线方程(6)将物方点(X,Y,Z)反投影至原始倾斜影像S上像素点(x,y)；作为本专利技术进一步改进的，步骤1.3中，采用双线性内插法计算目标影像D像素点(x′,y′)的灰度值：g′(x′,y′)＝g(x,y)＝(1-i)*(1-j)*g(x-i,y-j)+(1-i)*j*g(x-i,y-j+1)+i*(1-j)*g(x-i+1,y-j)+i*j*g(x-i+1,y-j+1)(7)其中g′为目标影像D对应的灰度值，g为原始倾斜影像S对应的灰度值，i，j分别为像点横纵坐标的小数部分。作为本专利技术进一步改进的，步骤1.3中灰度填充法具体操作如下：从左至右，从上至下逐像素遍历目标影像D，对于空白目标影像D的每个像素点(x′,y′)，通过反投影获得它在原始倾斜影像S上的像平面坐标(x,y),并通过灰度内插获得(x,y)点的灰度值，作为目标影像D上像素点(x′,y′)的灰度值。作为本专利技术进一步改进的，步骤2和步骤3中，用四边形abcd和透视投影影像D的面积比作为影像对之间的重叠度，如果重叠度大于阈值，则对影像对的重叠区域执行匹配，否则认为影像对间不存在匹配点，不对影像对进行匹配。本专利技术的一种基于平面透视投影的倾斜影像匹配方法与现有技术相比，所产生的有益效果是，本专利技术的技术方案是将原始的倾斜影像进行透视投影到一个平面上，对经过透视投影的影像进行匹配，获取匹配结果。相对于ASIFT算法对影像先进行仿射变换再进行匹配的方式，本专利技术经过透视变换，使得影像对之间的几何差异更小，相似特征更加明显，更容易匹配。附图说明附图1为本专利技术采用的流程图。附图2为影像平面透视投影示意图。附图3为利用F矩阵剔除影像匹配点对粗差流程图。附图4为平面透视投影前后对比图。具体实施方式为了更好的理解本专利技术的技术方案，下面结合附图1-4对本专利技术做进一步的详细说明。本专利技术的实施例对三对倾斜航空影像进行匹配。参照图1，本专利技术实施例流程的步骤如下：步骤1，对待匹配影像对中的每一张影像进行平面透视投影。利用原始倾斜影像S的外方位元素和内方位元素，将原始倾斜影像S通过透视投影到测区内的平均高程面Hm上，获得近似的正射影像，作为平面透视投影之后的目标影像D。实施例的步骤1中平面透视投影具体实现步骤如下：步骤1.1，按照公式(1)将原始倾斜影像S的四个角点(xi,yi)(其中i＝0,1,2,3)，通过透视投影投影至目标影像D上的物方点(Xi，Yi，Hm)(其中i＝0,1,2,3)，(xi,yi)对应目标影像D上的四个角点(x′i,y′i)(其中i＝0,1,2,3)。其中，f为焦距，a1、a2、a3、b1、b2、b3、c1、c2、c3为旋转矩阵的元素，Xs、Ys、Zs为摄影中心坐标。步骤1.2，利用目标影像四个角点(x′i,y′i)，根据公式(2)和公式(3)分别确定目标影像D的边界，以及目标影像本文档来自技高网...
【技术保护点】
1.一种基于平面透视投影的倾斜影像匹配方法，其特征在于，包括以下步骤：/n步骤1，利用原始倾斜影像S的外方位元素和内方位元素，将原始倾斜影像S通过透视投影到测区内的平均高程面H
【技术特征摘要】
1.一种基于平面透视投影的倾斜影像匹配方法，其特征在于，包括以下步骤：
步骤1，利用原始倾斜影像S的外方位元素和内方位元素，将原始倾斜影像S通过透视投影到测区内的平均高程面Hm上，获得近似的正射影像，作为平面透视投影之后的目标影像D；
步骤2，获取目标影像对之间的重叠区域，作为待匹配区域：对于两张待匹配的原始倾斜影像，通过步骤1，分别获取它们的目标影像对及其范围，根据多边形求交运算获取目标影像对之间的重叠区域四边形abcd，将重叠区域abcd作为待匹配区域；
步骤3，采用SIFT算法对目标影像对进行匹配，使用左右一致性检验、唯一性约束等对粗差点进行剔除，获得粗匹配结果；通过单应矩阵和基础矩阵双重约束对粗匹配结果进行筛选，获得最终的精匹配结果；
步骤4，将步骤3获得的目标影像对上的匹配点对，利用共线方程反投影至原始倾斜影像对上，作为最终匹配结果。
2.根据权利要求1所述的基于平面透视投影的倾斜影像匹配方法，其特征在于，所述步骤1具体实现步骤如下：
步骤1.1，将原始倾斜影像S的四个角点(xi,yi)(其中i＝0,1,2,3)，通过透视投影投影至目标影像D上，对应目标影像D上的四个角点(x′i,y′i)(其中i＝0,1,2,3)；
步骤1.2，利用目标影像四个角点(x′i,y′i)，确定目标影像D的边界、宽W′和高H′，建立W′×H′的空白目标影像D；
步骤1.3，灰度填充法填充目标影像D。
3.根据权利要求2所述的基于平面透视投影的倾斜影像匹配方法，其特征在于，步骤1.1中所述目标影像D所在的高程面固定，为整个测区内的平均高程Hm；原始倾斜影像S上任意一像点(x,y)，通过公式(1)经过透视投影至目标影像D上的物方点坐标为(X，Y，Hm)，(X，Y)为目标影像D的像点坐标(x′i,y′i)；
其中，f为焦距，a1、a2、a3、b1、b2、b3、c1、c2、c3为旋转矩阵的元素，Xs、Ys、Zs为摄影中心坐标。
4....
【专利技术属性】
技术研发人员：李治远，吴子敬，胡俊，豆虎林，
申请(专利权)人：国家海洋局第一海洋研究所，
类型：发明
国别省市：山东;37
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。