基于K近邻图的误匹配点迭代检测方法技术

本发明专利技术针对现有误匹配点检测方法的不足，提出了一种基于K近邻图的误匹配点迭代检测方法，可以提高误匹配点检测精度。处理过程为：(1)将匹配点分为若干组，对每一组匹配点分别按照步骤(2)-(5)进行误匹配点检测；(2)在原始图像上构建K近邻图；(3)对于每一对匹配点，通过K近邻图获得距离最近的K对匹配点，然后使用K对匹配点建立变换模型，将该对匹配点和变换模型的误差作为该对匹配点的误差；(4)删除误差最大的S匹配点；(5)迭代检测终止判断。如果所有剩余匹配点的误差都小于给定阈值，则停止迭代；否则N＝N-S，返回步骤(2)。

【技术实现步骤摘要】
基于K近邻图的误匹配点迭代检测方法
本专利技术涉及图像处理技术，具体地说，是一种基于K近邻图的误匹配点迭代检测方法，用于检测并删除遥感图像匹配中产生的误匹配点，提高遥感图像匹配精度。
技术介绍
遥感图像匹配是遥感图像处理中的一个重要步骤，在遥感图像配准、目标识别、目标跟踪、变化检测、图像拼接等领域得到了广泛的应用。遥感图像匹配是通过一定的匹配算法在两幅图像之间寻找的高精度匹配点。在遥感图像匹配中，待匹配的遥感图像称为原始图像，作为匹配基准的图像称为基准图像。近几年来，随着图像匹配技术的不断发展，匹配精度和自动化水平不断提高，尤其是以SIFT为代表的尺度、旋转不变特征匹配算法，使得遥感图像匹配的精度得到了极大的提升。但是，由于光照、地物变化、不同成像角度和成像时间、几何变形等因素的影响，无论采用何种图像匹配算法，总会产生误匹配点。因此，对于遥感图像匹配，误匹配点检测是一个必须步骤，其精度直接影响图像匹配的精度。在遥感图像匹配中，随机采样一致性方法(RandomSampleConsensus，RANSAC)和多项式拟合方法通常用于误匹配点检测。RANSAC方法是从一组包含异常数据的样本数据集中，估算最优模型参数(给定的两幅图像之间的变换模型)的迭代方法。在误匹配点检测中，通常使用基本矩阵或者单应矩阵作为变换模型，和变换模型误差较大的匹配点就是误匹配点。多项式拟合方法是假设所有正确的匹配点满足某一个多项式模型，然后使用最小二乘拟合等方法求解多项式的系数，最后计算每个匹配点和多项式模型的误差，误差大于设定阈值的匹配点就是误匹配点。虽然RNASAC方法和多项式拟合方法的原理、计算方法不同，但都是假设所有正确的匹配点满足一个统一的模型。但是对于部分覆盖范围较大的遥感图像，受复杂地形、几何变形等的影响，所有正确匹配点无法满足同一个模型，因此，常用的RANSAC方法和多项式拟合方法无法直接用于该类遥感图像的误匹配检测。
技术实现思路
本专利技术针对现有误匹配点检测技术的不足，提出了一种基于K近邻图的误匹配点迭代检测方法，准确的完成误匹配点检测，进一步提高图像匹配的精度。本专利技术实施的前提是使用图像匹配方法获得了匹配点。本专利技术的技术方案如下：一种基于K近邻图的误匹配点迭代检测方法，其特征在于包括以下步骤：(1)匹配点分组。对于待检测的N对匹配点，按照原始图像上匹配点的坐标，将匹配点分为若干组，然后对每一组匹配点分别按照步骤(2)-(5)进行误匹配点检测；(2)在原始图像上构建K近邻图。将原始图上的每一个匹配点作为K近邻图的一个顶点，对于每个顶点，如果其它顶点是该顶点的K近邻点，则这两个顶点存在一条无向边；顶点和顶点之间的无向边构成了K近邻图。(3)计算每一对匹配点的误差。对于每一对匹配点，通过K近邻图获得距离最近的K对匹配点，然后使用K对匹配点建立变换模型，计算每对匹配点和变换模型的误差作为该对匹配点的误差。(4)删除误差最大的匹配点。所有匹配点的误差计算完成后，删除误差最大的S对匹配点；(5)迭代检测终止判断。如果所有剩余匹配点的误差都小于给定阈值，则停止迭代；否则N＝N-S，返回步骤(2)。进一步，步骤(1)中的匹配点分组，采用网格划分的方法，首先对原始图像进行网格划分，网格大小为M*N，然后根据原始图像上匹配点的坐标，将所有匹配点分配到不同的网格；M和N的取值根据原始图像大小来确定，保证每一个网格都有上百个匹配点。进一步，步骤(2)中K的取值根据步骤(3)选择的变换模型来确定，K的取值要大于计算变换模型参数所需要的最少匹配点个数。进一步，步骤(2)中的K近邻图为无向图，使用邻接表存储，记为G＝(V，E)，其中V为顶点的集合，E为边的集合；构建准则如下：把原始图像中的每一个匹配点作为K近邻图的一个顶点，构成顶点集合V＝v1，...vN，N为原始图像中匹配点的总个数；对于顶点vi，如果顶点vj是vi的K最近邻域中的点，则vi和vj之间存在一条无向边；顶点vi的K最近邻域中的点是距离顶点vi最近的K个顶点；两个顶点之间的距离是顶点对应的原始图像上匹配点的欧式距离。进一步，步骤(3)中的变换模型为n阶多项式、基本矩阵、单应矩阵。进一步，步骤(3)中的匹配点误差计算公式为：上式中，(xi，yi)为匹配点对Mi在原始图像上的坐标，(x′i，y′i)为匹配点对Mi在基准图像上的坐标，按照构建的变换模型换算到原始图像上的坐标。本专利技术与现有技术相比具有下列优点：通过构建K近邻图，实现了只在图像局部建立变换模型。由于图像局部区域范围较小，因此，所有正确的匹配点都能满足同一个变换模型，从而很好的克服了传统RANSAC方法和多项式拟合方法的缺点，提高了误匹配点检测的精度。附图说明图1基于K近邻图的误匹配点迭代检测流程图具体实施方式现在结合附图，描述本专利技术的一种具体实施方式。图1基于K近邻图的误匹配点迭代检测流程图，包括五个步骤：(1)匹配点分组。对于大数据量的遥感影像，由于获得的匹配点较多，如果直接进行误匹配点检测，会导致速度比较慢，因此，首先对图像匹配获得的匹配点进行分组。匹配点分组采用网格划分的方法，首先对原始图像进行网格划分，网格大小为M*N，然后根据原始图像上匹配点的坐标，将所有匹配点分配到不同的网格。M和N的取值根据原始图像大小来确定，保证每一个网格都有上百个匹配点。不同网格的匹配点形成若干组匹配点。对于每一组匹配点，分别按照步骤(2)-(5)进行误匹配点检测。(2)在原始图像上构建K近邻图。图是数据之间关系的一种表示方式，由顶点的集合和边的集合(顶点之间的关系)构成。如果边没有方向，则称为无向图，如果边有方向，则称为有向图。本专利技术专利中，构建的K近邻图为无向图，记为G＝(V，E)，其中V为顶点的集合，E为边的集合。其中K的取值根据步骤(3)选择的变换模型来确定，K的值要大于计算变换模型参数所需要的最少匹配点个数。具体构建准则如下：把原始图像中的每一个匹配点作为K近邻图的一个顶点，构成顶点集合V＝v1，...vN，N为原始图像中匹配点的总个数。N的初始值为自动匹配方法获得的匹配点总个数，随着不断迭代检测和删除，N的值会不断减小。对于顶点vi，如果顶点vi是vi的K最近邻域中的点，则vi和vi之间存在一条无向边。顶点vi的K最近邻域中的点是距离顶点vi最近的K个顶点。两个顶点之间的距离是顶点对应的原始图像上匹配点的欧式距离。图的存储方法包括邻接矩阵和邻接表。邻接矩阵是一个n*n的矩阵，矩阵中存储边的信息。而邻接表是另一种存储结构，每个顶点vi的所有邻接点构成一个线性表，其优点是可以方便的获取顶点的K个邻接点。由于本专利技术专利中，需要通过K近邻图获得距离某一匹配点最近的K对匹配点，因此采用邻接表的存储方法。(3)计算每一对匹配点的误差。构建K近邻图后，对于每一个匹配点，使用其K个最近邻的匹配点计算其相对于给定模型的误差。对于每一对匹配点Mi，i＝1，...，N，首先根据原始图像上的匹配点，从K近邻图中获取K对匹配点，然后使用K对匹配点建立给定的变换模型。遥感图像中通常使用n阶多项式、基本矩阵、单应矩阵等变换模型。可以选择任意一种变换模型进行误差计算。在本专利技术的实施例中，使用一阶多项式模型，多项式的系数采用最小二乘拟合方法计算获得；最本文档来自技高网...

【技术保护点】
一种基于K近邻图的误匹配点迭代检测方法，用于检测并删除遥感图像匹配中产生的误匹配点，提高遥感图像匹配精度，其特征在于包括如下步骤：(1)匹配点分组，对于待检测的N对匹配点，按照原始图像上匹配点的坐标，将匹配点分为若干组，然后对每一组匹配点分别按照步骤(2)‑(5)进行误匹配点检测；(2)构建K近邻图，使用原始图上的匹配点构建K近邻图；(3)计算每一对匹配点的误差，对于每一对匹配点，通过K近邻图获得距离最近的K对匹配点，然后使用K对匹配点建立变换模型，计算每对匹配点和变换模型的误差作为该对匹配点的误差；(4)删除误差最大的匹配点，所有匹配点的误差计算完成后，删除误差最大的S对匹配点；(5)迭代检测终止判断，如果所有剩余匹配点的误差都小于给定阈值，则停止迭代；否则N＝N‑S，返回步骤(2)。

【技术特征摘要】
1.一种基于K近邻图的误匹配点迭代检测方法，用于检测并删除遥感图像匹配中产生的误匹配点，提高遥感图像匹配精度，其特征在于包括如下步骤：(1)匹配点分组，对于待检测的N对匹配点，按照原始图像上匹配点的坐标，将匹配点分为若干组，然后对每一组匹配点分别按照步骤(2)-(5)进行误匹配点检测；(2)构建K近邻图，使用原始图上的匹配点构建K近邻图，K的取值根据步骤(3)选择的变换模型来确定，K的取值要大于计算变换模型参数所需要的最少匹配点个数；(3)计算每一对匹配点的误差，对于每一对匹配点，通过K近邻图获得距离最近的K对匹配点，然后使用K对匹配点建立变换模型，计算每对匹配点和变换模型的误差作为该对匹配点的误差；(4)删除误差最大的匹配点，所有匹配点的误差计算完成后，删除误差最大的S对匹配点；(5)迭代检测终止判断，如果所有剩余匹配点的误差都小于给定阈值，则停止迭代；否则N＝N-S，返回步骤(2)。2.根据权利要求1中所述的方法，其特征在于：步骤(1)中的匹配点分组，采用网格划分的方法，首先对原始图像进行网格划分，网格大小为M*N，然后根据原始图像上匹配点的坐标，将所有匹配点分配...

【专利技术属性】
技术研发人员：单小军，胡昌苗，李宏益，郑柯，唐亮，唐娉，
申请(专利权)人：中国科学院遥感与数字地球研究所，
类型：发明
国别省市：北京;11