[0006] (1)计算输入图像中每个特征点iptYq描述子向量与注册图像中所有特征点fpctX描述子向量的汉明距离，取出最小距离mindis的一对点(iptYq,fpctXp)，1≤p≤n,1≤q≤m和次小距离secdis的一对点(iptYq,fpctXt)，1≤t≤n,1≤q≤m。若mindis
[0007] (2)通过RANSAC(Random Sample Consensus的缩写，它是根据一组包含异常数据的样本数据集，计算出数据的数学模型参数，得到有效样本数据的算法)方法对匹配点对集合kpsList中的错误匹配点对进行剔除。通过随机选点
迭代计算单应矩阵，剔除不满足该单应矩阵变换关系的点对，直到找到一个最佳单应矩阵。
[0008] (3)通过单应矩阵将注册图像提取的所有特征点坐标映射到输入图像
坐标系下，计算在局部区域内余弦相似度cosine最大且点对欧式距离和特征点主方向偏差在给定
阈值范围内的点对，记为一对匹配点对。最终得到所有的匹配点对。
[0009] 在实际操作过程中，上述方法存在如下问题：
[0010] (1)由于指纹图像局部区域相似性较大且图片会有畸变，往往存在很多假匹配点对比真匹配点对的汉明距离更小，因此通过描述子向量汉明距离的方式计算匹配点对的准确率较低。
[0011] (2)RANSAC方法是一种随机抽样，每次的选点的结果不一致。首先迭代次数无法控制，其次因为是随机选择，每次选取计算单应矩阵的点可能不同，从而导致每次匹配结果可能不同。对相同的一对图像进行匹配，多次实验可能得到不同的匹配结果，甚至原本可以匹配的图像在某次匹配中判断为图像不匹配。
发明内容
[0012] 针对上述
现有技术存在的问题，本发明提供一种基于有序抽样的特征匹配算法，能够计算出更加精准稳定的单应矩阵，从而在匹配点对时效果更好。
[0013] 为了实现上述目的，本发明采用如下技术方案：一种基于有序抽样的特征匹配算法，记待匹配的两幅图像分别为注册图像和输入图像，fpctX1，fpctX2，...，fpctXn为注册图像提取的n个特征点；iptY1，iptY2，...，iptYm为输入图像提取的m个特征点；其特征在于，整个特征匹配算法具体包括如下步骤：
[0014] S1、计算全局单向最优点对：计算注册图像中的每个特征点描述子向量与输入图像中所有特征点描述子向量的余弦相似度cosine，并记录余弦相似度cosine值最大的一对点(fpctXp,iptYq)，1≤p≤n,1≤q≤m为单向最优点对，最终得到注册图像指向输入图像的全局单向最优点对List1，记点对个数Count1；同理，计算输入图像指向注册图像的全局单向最优点对List2，记点对个数Count2；
[0015] S2、计算全局双向最优点对：基于以上得到的两个全局单向最优点对List1、List2计算全局双向最优点对List3，记点对个数为Count3；
[0016] S3、双向最优点对排序：根据特征点之间的相关性对全局双向最优点对List3进行排序；
[0017] S4、单应矩阵：取出t个最相关的双向最优点对计算出一个单应矩阵，并计算此单应矩阵的畸变值confid，4≤t≤Count3；通过此单应矩阵对List3中所有的注册图像点坐标进行映射，统计满足该单应矩阵变换关系的双向最优点对个数sum；遍历List3集合，多次选点计算不同的单应矩阵；
[0018] 1)若当前计算的单应矩阵下满足变换关系的点对个数sum比之前记录的sum大，则保留当前计算的单应矩阵，并更新sum与confid的值；
[0019] 2)若当前计算的单应矩阵下满足变换关系的点对个数sum与之前记录的sum相等，并且畸变值confid比之前的畸变值confid小，则保留当前计算的单应矩阵，并更新sum与confid；
[0020] 最终，计算得到了一个最优的单应矩阵；
[0021] S5、特征点匹配：通过最优的单应矩阵将注册图像提取的所有特征点坐标映射到输入图像坐标系下，计算在局部区域内余弦相似度cosine最大且点对欧式距离和特征点主方向偏差在给定阈值范围内的点对，记为一对匹配点对；最终得到所有的匹配点对；
[0022] S6、计算得分：根据匹配的点对数进行打分，保留匹配得分作为此注册图像与输入图像的匹配得分。
[0023] 所述步骤S2中双向最优点对是指若全局单向最优点对List1中存在点对(fpctXi,iptYj)，且全局单向最优点对List2存在点对(iptYj,fpctXi)，则把点对(fpctXi,iptYj)记录为一对双向最优点对，其中，1≤i≤n,1≤j≤m；全局双向最优点对List3是指所有双向最优点对的集合。
[0024] 与现有的技术相比，本发明先通过特征点描述子向量余弦相似度计算出双向最优匹配点对，再根据特征点的相关性把粗匹配下得到的双向最优点对进行排序，然后选择相关性相近的t个点对计算一个变换矩阵，然后对双向最优点对进行变换得到匹配的点对数；迭代计算变换矩阵，直到找到一个最佳变换矩阵，然后对所有的点对进行变换得到最终匹配结果。本发明能够计算出更加精准稳定的单应矩阵，从而在匹配点对时效果更好，可确保对一对图像进行多次匹配时得到的匹配结果相同，避免出现原本可以匹配的图像在某次匹配中判断为图像不匹配的情况。
附图说明
具体实施方式
[0026] 下面结合附图对本发明实施中的技术方案进行清楚，完整的描述，所描述的
实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0027] 如图1所示，本发明实施例提供的一种基于有序抽样的特征匹配算法，记待匹配的两幅图像分别为注册图像和输入图像，fpctX1，fpctX2，...，fpctXn为注册图像提取的n个特征点；iptY1，iptY2，...，iptYm为输入图像提取的m个特征点；整个特征匹配算法包括如下步骤：
[0028] S1、计算全局单向最优点对：计算注册图像中的每个特征点描述子向量与输入图像中所有特征点描述子向量的余弦相似度cosine，并记录余弦相似度cosine值最大的一对点(fpctXp,iptYq)，1≤p≤n,1≤q≤m为单向最优点对，最终得到注册图像指向输入图像的全局单向最优点对List1，记点对个数Count1；同理，计算输入图像指向注册图像的全局单向最优点对List2，记点对个数Count2；其中，全局单向最优点对List1是指所有注册图像指向输入图像的单向最优点对集合，Count1＝n；全局单向最优点对List2是指所有输入图像指向注册图像的单向最优点对集合，Count2＝m；
[0029] S2、计算全局双向最优点对：基于以上得到的两个全局单向最优点对List1、List2计算全局双向最优点对List3，记点对个数为Count3；其中，双向最优点对是指若全局单向最优点对List1中存在点对(fpctXi,iptYj)，且全局单向最优点对List2存在点对(iptYj,fpctXi)，则把点对(fpctXi,iptYj)记录为一对双向最优点对，其中，1≤i≤n,1≤j≤m；全局双向最优点对List3是指所有双向最优点对的集合；
[0030] S3、双向最优点对排序：根据特征点(指双向最优点对集合里的点)之间的相关性对全局双向最优点对List3进行排序；
[0031] S4、单应矩阵：取出t个最相关的双向最优点对计算出一个单应矩阵，并计算此单应矩阵的畸变值confid，4≤t≤Count3；通过此单应矩阵对List3中所有的注册图像点坐标进行映射，统计满足该单应矩阵变换关系的双向最优点对个数sum；遍历List3集合，多次选点计算不同的单应矩阵，例如有六个点对，排序后为1,2,3,4,5,6，t取4，则有以下几个选点方式：1,2,3,4；1,2,3,5；1,2,3,6；2,3,4,5；2,3,4,6；3,4,5,6；
[0032] 1)若当前计算的单应矩阵下满足变换关系的点对个数sum比之前记录的sum大，则保留当前计算的单应矩阵，并更新sum与confid的值；
[0033] 2)若当前计算的单应矩阵下满足变换关系的点对个数sum与之前记录的sum相等，并且畸变值confid比之前的畸变值confid小，则保留当前计算的单应矩阵，并更新sum与confid；
[0034] 最终，计算得到了一个最优的单应矩阵；
[0035] S5、特征点匹配：通过最优的单应矩阵将注册图像提取的所有特征点坐标映射到输入图像坐标系下，计算在局部区域内余弦相似度cosine最大且点对欧式距离和特征点主方向偏差在给定阈值范围内的点对，记为一对匹配点对；最终得到所有的匹配点对；
[0036] S6、计算得分：根据匹配的点对数进行打分，保留匹配得分作为此注册图像与输入图像的匹配得分。
[0037] 综上所述，本发明通过特征点描述子向量余弦相似度计算出双向最优匹配点对，再根据特征点的相关性把粗匹配下得到的双向最优点对进行排序，然后选择相关性相近的t个点对计算一个变换矩阵，然后对双向最优点对进行变换得到匹配的点对数；迭代计算变换矩阵，直到找到一个最佳变换矩阵，然后对所有的点对进行变换得到最终匹配结果。本发明能够计算出更加精准稳定的单应矩阵，从而在匹配点对时效果更好，可确保对一对图像进行多次匹配时得到的匹配结果相同，避免出现原本可以匹配的图像在某次匹配中判断为图像不匹配的情况。
[0038] 对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神和基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附
权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
[0039] 此外，应当理解，虽然本
说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。
