基于双方向图的快速局部不变特征描述子的构造方法技术

本发明专利技术提供了一种基于双方向图的快速局部不变特征描述子的构造方法，应用于图像处理技术领域；包括以下步骤：预处理、划分子区域、选取采样点、选取采样点的邻域点、预插值、获取邻域点的像素值、计算采样点的主方向区间和次方向区间、构造FRDOH子描述子、构造FRDOH描述子。该构造方法可使FRDOH描述子具有旋转不变性，提高FRDOH描述子的运算速度，同时可提高FRDOH描述子的稳定性和可区分性。

【技术实现步骤摘要】
基于双方向图的快速局部不变特征描述子的构造方法
本专利技术涉及图像处理
，更具体地说，涉及一种图像处理描述子的构造方法。
技术介绍
随着智能化时代的发展，基于机器视觉的自动识别技术得到了广泛的应用和关注。然而计算机远不及人脑的高度智能化，让计算机从大量视觉信息中提取出有用信息并准确识别目标便成了一个值得研究的课题，而各种局部不变性特征提取算法则提供了解决这一难题的有效手段。因此，局部不变特征在许多方面得到了应用，包括目标识别、3D重建、场景分类和全景拼接等。通常对局部不变特征的研究主要包括三个方面，分别是特征检测、特征描述以及特征匹配。在特征检测方面目前已有多种检测算子，例如Harris检测算子和Hessian检测算子分别基于Harris矩阵和Hessian矩阵检测角点，在此基础上应用Lindeberg提出的自动尺度选择理论，可以得到带有仿射不变性的Harris-Affine检测子和Hessian-Affine检测子。TimorKadir和MichaelBrady于2001年提出了Kadir-Brady显著性检测子，并以此来提取图像的显著性特征。在Kadir-Brady显著性检测子的基础上，通过后续不断的研究获得了具有仿射不变性和鲁棒性的两个改进版本。FAST(FeaturesfromtheAcceleratedSegmentTest)检测子则是通过比较邻域信息来找到兴趣点，而Matas在文献中则提出了一种检测区域的方法，称为最大稳定极值区域(MSERs)，该方法是通过检测图像中灰度最稳定的局部区域来获得具有仿射不变性的局部特征。在特征描述方面，特征描述子是一种图像局部结构特征的定量化数据描述，一个理想的特征描述子应该具备高的鲁棒性、可区分性和速度，这里的速度包括描述子构建速度和匹配速度。因此设计一个描述子的核心在于构造一个在鲁棒性、可区分性和速度方面都有优异表现的局部不变特征描述子。获得特征的描述子后则需要适当的特征匹配来完成物体识别。一般是先求取两个描述子之间的距离(如欧氏距离)，然后用最短距离或最短/次短距离比来计算相应的匹配分数，最后将匹配分数与预设的阈值相比来决定识别结果。局部不变特征描述子中影响最广的是SIFT描述子。SIFT描述子是在确定的主方向上进行梯度直方图的统计，因而具有良好的稳定性和准确性，，尤其是在一些具有尺度变换、旋转变换、视角变换和仿射变换的场合，SIFT描述子的准确性要远优于当时的其他描述子。但是，由于其主方向的选取比较费时，易受到统计误差的影响，且梯度向量的计算过程比较繁杂，因此SIFT描述子在实际应用中受到一定影响。但由于SIFT描述子应对一般图像变换的能力比较全面，许多研究人员都愿意在SIFT的基础上进行各种优化和改进。比如PCA-SIFT描述子是将主成分分析法跟SIFT描述子进行了综合，大大降低了SIFT的维数，也提高了描述子的准确率；GLOH描述子则是用极坐标表示代替栅格来进行梯度方向直方图的统计，实验证明这种表示方法能有效提高描述子的鲁棒性；SURF描述子引入了哈尔小波变换，在保证SIFT的部分优势的基础上实现降维并且极大提高了速度；DAISY描述子则是利用了梯度方向和像素位置这两个因素，在仿射变化和光强线性变换时具有良好的效果，但是由于其复杂的计算，DAISY描述子在速度方面则表现一般。由于SIFT描述子是通过寻找主方向的方法实现旋转不变性，而这种方法容易受到统计误差的影响，并非真正意义上的旋转不变，因此SIFT描述子的稳定性容易受到影响。SIFT的其他衍生的描述子沿用了SIFT描述子的主方向寻找方法，因此也存在这种缺陷。此外还有另外两个比较新的描述子，MROGH和RATMIC。MROGH描述子在计算采样点的梯度时运用了实现旋转不变性的方法，使得描述过程中不需要特意计算特征区域的主方向，这样描述子在多种情况下能更加稳定；同时MROGH描述子增加了描述子的可区分性和对非线性光照变化的适应能力。但是MROGH描述子引入了多支持区域的方法来提高整体的可区分能力，MROGH描述子的整体耗时增加了几倍，所以速度方面的表现难以令人满意。考虑到速度方面的需求，RATMIC描述子在MROGH描述子的基础上进行了修改。由于上面提到的多支持区域比较耗时，RATMIC描述子放弃了这种提高描述子可区分能力的方法，而是用了另外一种表示采样点值的方法来代替MROGH描述子中所用的计算梯度值的方法。这种方法是将采样点周边的4个邻域点进行比较并编码成为二进制串，即对每个点的表示维数从MROGH描述子中的8维变成了15维。从而提高描述子的可区分性。整体上看，RATMIC描述子在速度方面要远快于MROGH等高耗时的描述子，但在可区分性方面的表现略差于MROGH描述子。鉴于上述情况，需要设计出一种速度快、可区分性能好的描述子。
技术实现思路
本专利技术的目的在于克服现有技术中的缺点与不足，提供一种基于双方向图的快速局部不变特征描述子的构造方法；该构造方法可使FRDOH描述子具有旋转不变性，提高FRDOH描述子的运算速度，同时可提高FRDOH描述子的稳定性和可区分性。为了达到上述目的，本专利技术通过下述技术方案予以实现：一种基于双方向图的快速局部不变特征描述子的构造方法，其特征在于，包括以下步骤：第一步，进行预处理，预处理包括获取特征点Pj，j∈{1,...,KP}，其中KP为特征点Pj的总数；获取特征点Pj是通过带仿射不变性的检测子进行检测实现；第二步，通过预插值方法建立特征区域中亚像素像素值的查找表：对特征区域进行亚像素细分形成具有亚像素点的预插值区域，并采用双线性插值方式获得亚像素点的像素值，形成亚像素像素值查找表；第三步，将特征区域内所有像素点按照灰度值大小进行排序形成像素点序列，然后将像素点序列平均分成k个区间，不同区间分别对应着不同子区域，每个区间的像素点映射到相应的子区域，从而把特征区域划分成k个子区域；并设定各个子区域的子区域标记S,S∈{0,...,k-1}；第四步，采样点C是指特征区域内除去特征点Pj外的任意一个点，Ci(i＝1,...,12)则是以采样点C为中心选取的邻域点，对采样点C进行FRDOH子描述子函数构造；采样点C的FRDOH子描述子函数构造包括以下步骤：A步，选取采样点C的邻域点Ci(i＝1,...,12)：以采样点C为中心、射线为对角线方向、2L为对角线长度确定正方形一，将正方形一的四个顶点和四边中点设定为邻域点C1～邻域点C8；然后以正方形一的四个顶点作为正方形二的四边中点来确定正方形二，将正方形二的四个顶点设定为邻域点C9～邻域点C12；B步，将特征区域上的邻域点C1～邻域点C12映射到预插值区域上；依次获取与邻域点C1～邻域点C12最近邻的亚像素点，邻域点C1～邻域点C12的像素值分别等于最近邻的亚像素点的像素值；最近邻的亚像素点的像素值通过亚像素像素值查找表获得；C步，沿正方形一的对角线建立x-y坐标系一，将邻域点C1～邻域点C8的像素值投影到x-y坐标系一中，将主方向区间划分为m级，计算采样点C主方向区间O1，O1∈{0,...,m-1}；D步，沿正方形二的对角线建立x-y坐标系二，将邻域点C9～邻域点C12的像素值投影到x-y坐标系二中，将次方向本文档来自技高网...

【技术保护点】
一种基于双方向图的快速局部不变特征描述子的构造方法，其特征在于，包括以下步骤：第一步，进行预处理，预处理包括获取特征点Pj，j∈{1,...,KP}，其中KP为特征点Pj的总数；获取特征点Pj是通过带仿射不变性的检测子进行检测实现；第二步，通过预插值方法建立特征区域中亚像素像素值的查找表：对特征区域进行亚像素细分形成具有亚像素点的预插值区域，并采用双线性插值方式获得亚像素点的像素值，形成亚像素像素值的查找表；第三步，将特征区域内所有像素点按照灰度值大小进行排序形成像素点序列，然后将像素点序列平均分成k个区间，不同区间分别对应着不同子区域，每个区间的像素点映射到相应的子区域，从而把特征区域划分成k个子区域；并设定各个子区域的子区域标记S,S∈{0,...,k‑1}；第四步，采样点C是指特征区域内除去特征点Pj外的任意一个点，Ci(i＝1,...,12)则是以采样点C为中心选取的邻域点，对采样点C进行FRDPH子描述子函数构造；采样点C的FRDPH子描述子函数构造包括以下步骤：A步，选取采样点C的邻域点Ci(i＝1,...,12)：以采样点C为中心、射线为对角线方向、2L为对角线长度确定正方形一，将正方形一的四个顶点和四边中点设定为邻域点C1～邻域点C8；然后以正方形一的四个顶点作为正方形二的四边中点来确定正方形二，将正方形二的四个顶点设定为邻域点C9～邻域点C12；B步，将特征区域上的邻域点C1～邻域点C12映射到预插值区域上；依次获取与邻域点C1～邻域点C12最近邻的亚像素点，邻域点C1～邻域点C12的像素值分别等于最近邻的亚像素点的像素值；最近邻的亚像素点的像素值通过亚像素像素值查找表获得；C步，沿正方形一的对角线建立x‑y坐标系一，将邻域点C1～邻域点C8的像素值投影到x‑y坐标系一中，将主方向区间划分为m级，计算采样点C主方向区间O1，O1∈{0,...,m‑1}；D步，沿正方形二的对角线建立x‑y坐标系二，将邻域点C9～邻域点C12的像素值投影到x‑y坐标系二中，将次方向区间划分为n级，计算采样点C次方向区间O2，O2∈{0,...,n‑1}；E步，获取采样点C的子区域标记S；F步，建立采样点C的FRDPH子描述子：采样点C的FRDPH子描述子函数为γ(C)=O1+O2×m+S×m×nFRDOH(C)=φ(γ(C))=(0,...,0,1γ(C),0,...,0)]]>其中，γ(C)是采样点C映射到FRDPH子描述子中的位置；FRDOH(C)的长度＝m×n×s，且γ(C)位置处的值为1，其他位置的值为0；第五步，判断采样点C是否遍历特征区域Region(j)中的所有点：若是，则跳至第六步；否则，重复执行第四步；第六步，构造FRDOH描述子：FRDOH描述子函数为subdes(j)=ΣC∈Region(j)FRDOH(C),j∈{1,...,KP}]]>其中，KP为特征点Pj的总数，Region(j)为任一特征点Pj的特征区域。...

【技术特征摘要】
1.一种基于双方向图的快速局部不变特征描述子的构造方法，其特征在于，包括以下步骤：第一步，进行预处理，预处理包括获取特征点Pj，j∈{1,...,KP}，其中KP为特征点Pj的总数；获取特征点Pj是通过带仿射不变性的检测子进行检测实现；第二步，通过预插值方法建立特征区域中亚像素像素值的查找表：对特征区域进行亚像素细分形成具有亚像素点的预插值区域，并采用双线性插值方式获得亚像素点的像素值，形成亚像素像素值的查找表；第三步，将特征区域内所有像素点按照灰度值大小进行排序形成像素点序列，然后将像素点序列平均分成k个区间，不同区间分别对应着不同子区域，每个区间的像素点映射到相应的子区域，从而把特征区域划分成k个子区域；并设定各个子区域的子区域标记S,S∈{0,...,k-1}；第四步，采样点C是指特征区域内除去特征点Pj外的任意一个点，Ci(i＝1,...,12)则是以采样点C为中心选取的邻域点，对采样点C进行FRDOH子描述子函数构造；采样点C的FRDOH子描述子函数构造包括以下步骤：A步，选取采样点C的邻域点Ci(i＝1,...,12)：以采样点C为中心、射线PjC为对角线方向、2L为对角线长度确定正方形一，将正方形一的四个顶点和四边中点设定为邻域点C1～邻域点C8；然后以正方形一的四个顶点作为正方形二的四边中点来确定正方形二，将正方形二的四个顶点设定为邻域点C9～邻域点C12；B步，将特征区域上的邻域点C1～邻域点C12映射到预插值区域上；依次获取与邻域点C1～邻域点C12最近邻的亚像素点，邻域点C1～邻域点C12的像素值分别等于最近邻的亚像素点的像素值；最近邻的亚像素点的像素值通过亚像素像素值查找表获得；C步，沿正方形一的对角线建立x-y坐标系一，将邻域点C1～邻域点C8的像素值投影到x-y坐标系一中，将主方向区间划分为m级，计算采样点C主方向区间O1，O1∈{0,...,m-1}；D步，沿正方形二的对角线建立x-y坐标系二，将邻域点C9～邻域点C12的像素值投影到x-y坐标系二中，将次方向区间划分为n级，计算采样点C次方向区间O2，O2∈{0,...,n-1}；E步，获取采样点C的子区域标记S；F步，建立采样点C的FRDOH子描述子：采样点C的FRDOH子描述子函数为其中，γ(C)是采样点C映射到FRDOH子描述子中的位置；FRDOH(C)的长度＝m×n×s，且γ(C)位置处的值为1，其他位置的值为0；第五步，判断采样点C是否遍历特征区域Region(j)中的所有点：若是，则跳至第六步；否则，重复执行第四步；第六步，构造FRDOH描述子：FRDOH描述子函数为其中，KP为特征点Pj的总数，Region(j)为任一特征点Pj的特征区域。2.根据权利...

【专利技术属性】
技术研发人员：康文雄，陈晓鹏，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：广东;44