[发明专利]一种基于惩罚控制竞争学习聚类算法的特征分层匹配方法

摘要

本发明公开了一种基于惩罚控制竞争学习聚类算法的特征分层匹配方法，其采用一种能够自动确定类别数的惩罚控制竞争学习聚类算法得到参考特征集和待匹配特征集的聚类结果，然后用参考特征集中的所有中心单元构建一个父层K‑D树，用每个类中的每个样本作为节点组成一个子层K‑D树，进行匹配时将待匹配特征集的每个类的中心单元作为父查询目标在父层K‑D树中找出最近邻节点，再将待匹配特征集的每个类中的每个样本作为子查询目标，在该类对应的父查询目标的最近邻节点所属类的子层K‑D树中找出该类对应的每个子查询目标的最近邻节点，将每个子查询目标与其最近邻节点作为一个特征匹配对；优点是特征匹配对的数目不多，降低了计算复杂度，提高了图像匹配效率。

主权项

一种基于惩罚控制竞争学习聚类算法的特征分层匹配方法，其特征在于包括以下步骤：①选取不同时间、不同成像条件下获取的包含同一目标对象的尺寸大小一致的两幅图像，随机选择其中一幅图像作为参考图像，将另一幅图像作为待匹配图像；②利用SIFT算子提取参考图像和待匹配图像的所有SIFT特征，将由参考图像的所有SIFT特征构成的集合定义为参考特征集，将由待匹配图像的所有SIFT特征构成的集合定义为待匹配特征集，其中，每个SIFT特征为一个128×1维向量；③采用惩罚控制竞争学习聚类算法，分别对参考特征集和待匹配特征集进行聚类操作，得到参考特征集和待匹配特征集各自的最终的类别数以及每个类的中心单元和每个类各自拥有的样本；所述的步骤③的具体过程为：③‑1、将参考特征集中的每个SIFT特征中的梯度信息作为一个样本，每个样本为一个128×1维的向量；③‑2、令K表示类别数，假设K的初始值为K0，K0∈[10,50]；令m表示迭代次数，m的初始值为1；③‑3、从所有样本中随机选取K0个样本，将选取的每个样本作为一个种子单元，假定每个种子单元代表一个类的中心单元；③‑4、将第c个种子单元yc作为获胜单元，然后从所有样本中随机抽取1个样本，假设该样本为第t个样本，则将该样本记为xt，接着计算除获胜单元yc外的每个种子单元与xt之间的欧氏距离，将最小的欧氏距离值对应的种子单元作为竞争单元，假设该竞争单元为所有种子单元中的第r个种子单元yr，再执行步骤③‑5；其中，1≤c≤K0，c的初始值为1，1≤t≤T，T表示参考特征集对应的所有样本的总个数，1≤r≤K，r≠c；③‑5、根据获胜单元yc在本次迭代前的获胜次数nc，计算获胜单元yc在本次迭代中的学习率，记为αc，其中，nc的初始值为0，const为平衡参数；然后对获胜单元yc进行更新学习以使获胜单元yc朝xt移动靠近，将更新学习后的获胜单元记为ycnew，ycnew＝yc+αc(xt‑yc)；再令yc＝ycnew，并令nc＝nc+1，其中，yc＝ycnew和nc＝nc+1中的“＝”为赋值符号；同样，根据竞争单元yr在本次迭代前的获胜次数nr，计算竞争单元yr在本次迭代中的惩罚率，记为αr，其中，nr的初始值为0，const为平衡参数；然后对竞争单元yr进行惩罚以使竞争单元yr进行反向移动远离xt，将惩罚后的竞争单元记为yrnew，yrnew＝yr‑ρ(xt‑yr)，其中，ρ＝αr×DisRatio，符号“|| ||”为求欧氏距离符号；再令yr＝yrnew，其中，yr＝yrnew中的“＝”为赋值符号；③‑6、令c＝c+1，如果c>K0，则执行步骤③‑7，如果c≤K0，则返回步骤③‑4继续执行；③‑7、从所有样本中随机抽取1个样本，假设该样本为第t个样本，则将该样本记为xt，然后计算每个种子单元与xt之间的距离，接着将最小的距离值对应的种子单元作为获胜单元，假设该获胜单元为所有种子单元中的第c个种子单元yc，将值次小的距离对应的种子单元作为竞争单元，假设该竞争单元为所有种子单元中的第r个种子单元yr，再执行步骤③‑8；其中，1≤c≤K0，1≤r≤K0；③‑8、根据获胜单元yc在本次迭代前的获胜次数nc，计算获胜单元yc在本次迭代中的学习率，记为αc，其中，nc的初始值为0，const为平衡参数；然后对获胜单元yc进行更新学习以使获胜单元yc朝xt移动靠近，将更新学习后的获胜单元记为ycnew，ycnew＝yc+αc(xt‑yc)；再令yc＝ycnew，并令nc＝nc+1，其中，yc＝ycnew和nc＝nc+1中的“＝”为赋值符号；同样，根据竞争单元yr在本次迭代前的获胜次数nr，计算竞争单元yr在本次迭代中的惩罚率，记为αr，其中，nr的初始值为0，const为平衡参数；然后对竞争单元yr进行惩罚以使竞争单元yr进行反向移动远离xt，将惩罚后的竞争单元记为yrnew，yrnew＝yr‑ρ(xt‑yr)，其中，ρ＝αr×DisRatio，符号“|| ||”为求欧氏距离符号；再令yr＝yrnew，其中，yr＝yrnew中的“＝”为赋值符号；③‑9、重复执行步骤③‑7至步骤③‑8共T‑K0次；然后令m＝m+1，如果m>50，则结束迭代，再执行步骤③‑10，如果m≤50，则令nc＝0,nr＝0，再返回步骤③‑4继续执行，其中，m＝m+1中的“＝”为赋值符号；③‑10、将获胜次数小于次的种子单元剔除；然后统计剩余的种子单元的个数，将统计得到的个数作为参考特征集最终的类别数，并将剩余的每个种子单元重新确定为一个类的中心单元；接着计算参考特征集对应的每个样本与重新确定的各个类的中心单元之间的欧氏距离，将每个样本相对应的多个欧氏距离值中最小的欧氏距离值对应的中心单元的类标签作为该样本的类标签，从而得到各个类各自拥有的样本；按照步骤③‑1至步骤③‑10的操作过程，以相同方式获取待匹配特征集的最终的类别数，并确定每个类的中心单元，再确定待匹配特征集对应的每个样本的类标签；④根据参考特征集的每个类的中心单元和每个类中的所有样本，利用K‑D树构建方法构建得到参考特征集对应的父层K‑D树和每个类对应的子层K‑D树，然后将待匹配特征集的每个类的中心单元作为一个父查询目标在父层K‑D树中找出对应的最近邻节点，再将待匹配特征集的每个类中的每个样本作为一个子查询目标，在该类对应的父查询目标的最近邻节点所属类的子层K‑D树中找出该类对应的每个子查询目标的最近邻节点，将每个子查询目标与其最近邻节点作为待匹配图像与参考图像之间的一个特征匹配对，完成特征分层匹配。