一种基于SURF特征点匹配和RANSAC模型估计的牙齿X-ray影像匹配方法组成比例

本发明专利技术涉及一种基于SURF特征点匹配和RANSAC模型估计的牙齿X-ray影像匹配方法，包括：（1）图像采集；（2）图像增强；（3）特征提取；（4）特征的匹配；（5）计算相关匹配系数；（6）匹配结果；（7）若图库中还有没被匹配的参考图像，则返回(1)中继续在图库中选取一张新的参考图像，否则退出；（8）匹配成功，记录下参考图像对应的个人信息，退出。本发明专利技术利用SURF算法来提取输入牙齿影像和参考牙齿影像的特征点，然后利用RANSAC算法对各特征点进行匹配，最后通过匹配后的两幅图像之间的相关系数来测定两幅图像的匹配程度。实验表明，该算法实现了高查准率和高实时性。

【技术实现步骤摘要】
一种牙齿X-ray影像匹配方法
本专利技术涉及一种牙齿X-ray影像匹配方法，属于刑侦信息的

技术介绍
无名尸骨的身份鉴定，是司法领域的关键课题。随着恶性刑事案件以及爆炸、空难、海啸、火灾、车祸等灾难性事故的不断出现，诸多基于有机物质的身份鉴定手段，如唇印、指纹、虹膜等，常常无法起到理想的身份鉴定作用，甚至完全丧失价值。对于此类难以鉴别的无名尸骨，牙齿形态将成为独当一面的身份特征。牙齿是人体最坚硬的组织，具有较强的耐热性和抗腐蚀性。对于不同的个体，牙齿特征完全相同的概率仅有25亿分之一，因此被法齿学家称作人的另一张“身份证”。鉴于这种稳定性和差异性，基于牙齿特征的身份鉴定将成为死亡者身份鉴定最有潜力和价值的个体识别标志群。并且多数国家的司法局确认了法齿学在法庭科学中的重要地位。随着口腔医学的发展，人们自我保健、爱牙、护牙的意识不断提升，大众牙科医疗已趋于普及。牙齿影像已成为口腔诊疗中的必备资料，特别是全口曲面断层片投照具有辐射小成本低的特点，目前牙科诊疗中运用最广泛的标准初始检验手段。大量积累和保存的牙病病历档案，为基于牙齿特征的身份鉴定，提供了必要条件。近十年来，人类牙齿X-ray影像的匹配算法，一直是法齿学中的热点课题。世界各地的科学家都相继发表了相关领域的文献。2003年12月，AnilK.Jain和HongChen在文献[1]中提出下列算法：首先利用灰度值累和法来确定牙齿图像的分割边界，并利用该分割边界把牙齿图像分割成单颗牙齿；然后利用贝叶斯概率方法和概率方法，对牙冠部分提取边缘，并利用局部灰度值优化算法，对牙根部分提取边缘；再把提取到的边缘点作为特征点，将输入图像的边缘进行适当的相似变换后，计算输入图像边缘和参考图像边缘的欧式距离；最后选择对应欧式距离最小的若干幅参考图像作为最终的结果图像。该算法虽然比较经典，但提取边缘的过程仍然是个半自动化的算法。并且对于光照不均匀的图像，牙齿分割算法比较繁琐。对于边缘模糊的图像，牙齿的边缘很难提取。2004年8月，MohammadH.Mahoor和MohamedAbdel-Mottaleb把医学中的牙齿分类引入牙齿影像处理中[2]，并按照医学上的标准牙齿序列对牙齿进行标号。在匹配环节，只把对应标号的牙齿进行匹配，从而提高了系统的实时性。2004年12月，OmaimaNomir和MohamedAbdel-Mottaleb提出了新的牙齿分割算法和牙齿图像匹配算法[3]。在牙齿分割过程中，该算法把原图像Canny边缘处灰度值的均值作为阈值，将图像分割成牙齿部分和背景部分。在匹配过程中，该算法选取了边缘上若干个曲率极大点作为特征点，并把这些特征点之间的几何距离作为特征向量，从而减小了匹配的复杂度。2005年4月，GamalFahmy等人在文献[4]中系统地提出了自动化牙齿影像识别系统（ADIS）的架构。该架构主要包含数字图像存储器（DIR，包括存储功能和查询功能），潜在检索匹配模块（包括咬翼片图像的边缘提取功能和牙根尖图像的牙冠牙根边缘提取功能），图像匹配模块（包括增强、分割、匹配三个环节）三部分。模块化提高了系统的易维护性。2005年8月，AnilK.Jain和HongChen在文献[5]中优化了两个关键算法，即特征提取和相似度匹配。特征提取不仅提取了牙齿的边缘，还提取了牙齿的修补区域。相似度匹配算法通过后验概率将边缘距离和修补区面积融合，先计算出单颗牙齿层面的距离，再计算出图像层面的距离，最后根据图像层面的距离来检索最匹配的图像。多特征的融合提高了系统的有效性。2006年4月，DiaaEldinM.Nassar把人工神经网络算法引入到医学牙齿影像识别[6]。该神经网络分为特征提取和判断匹配两个层次。特征提取层利用非线性滤波器对输入图像和参考图像进行特征过滤，然后将过滤得到的特征传递到判断匹配层。最后判断匹配层判断两幅图像是否匹配。该算法可以通过自学习过程自动修正网络权系数，使得系统具有更明显的灵活性。2006年6月，EyadHajSaid等人利用形态学算法对牙齿图像进行去噪、分割和标记[7]，简化了传统的灰度值累和法，从而进一步简化了整个系统的设计。2007年5月，OmaimaNomir和MohamedAbdel-Mottaleb提出了多分辨率配准的算法[8]，认为图像的边缘特征可以先从低分辨率进行粗匹配，粗匹配后滤除一部分参考图像，然后将输入图像与剩下的参考图像进行高分辨率的精细匹配，最后筛选出若干幅最佳匹配图像。该算法通过逐层筛选的方法，提高了系统的实时性。2007年6月，OmaimaNomir和MohamedAbdel-Mottaleb将牙齿影像的表观特征引入鉴定系统[9]，认为牙齿特征不止包括边缘特征，也包括灰度特征。本文利用引力场算法求解出牙齿图像中的灰度特征点，并将灰度特征点与边缘特征点融合，使得最终的匹配效果更佳可靠。2008年6月，OmaimaNomir和MohamedAbdel-Mottaleb对牙齿的三种特征进行了融合[10]。他们整合了牙齿的边缘特征、表观特征、多分辨率特征，并分别在特征层、匹配层、决策层进行融合。实验结果证明了融合后的特征可以提高系统的整体性能。同时，DiaaEldinNassar等人提出了一种建立牙齿数据库的方法[11]。该方法分为两个步骤。首先，根据牙齿的视觉特征，将图库中的每一颗牙齿进行分类存储（门齿、犬齿、前臼齿、臼齿）。其次，根据牙齿的类别和位置进行标记，同时验证牙齿是否被正确分类。数据库建立完成后，可以将输入的牙齿图像与相应位置的牙齿图像进行对比，从而降低了计算量，提高了系统的实时性。2009年11月，Febriliyan等人提出了一种在全齿图像中自动提取臼齿的算法[12]。作者认为，臼齿的形态包含很多身份特征信息，可以完全通过臼齿的特征来完成身份鉴定。该算法首先选择全齿图像的关键区域，然后分割上颌牙和下颌牙，最后通过边缘信息来选择所需要提取的臼齿。但该算法仍然是一个半自动化的过程，需要用户在全齿图像中选择三个关键点。2013年4月，VijayakumariPushparaj等人将牙齿的相片也列入特征提取范畴[13]。本文指出，牙齿影像和牙齿相片都包含身份特征，因此引入牙齿相片的分析，可以提高鉴定系统的可靠性。并且统计数据显示，在牙齿相片中，上门齿的匹配就能达到70%的命中率。此外，除了利用牙齿的边缘特征外，本文还提出了牙齿的骨架特征，从而进一步提高了系统的可靠性。但以上的算法都存在着缺陷：首先，牙齿分割会破坏牙齿的特征。因为牙齿的特征不仅仅集中在边缘上，牙缝的形状、齿间的重合区域、牙齿的生长角度等也都可以作为特征，而牙齿分割就会割裂牙缝的形状。其次，所有文献都把牙齿的边缘作为特征，但是边缘未必是最有效的特征，因此特征提取的过程都缺少自适应性。再次，医学影像往往都是质量较差、噪声较多、光照不均匀、边缘较模糊的图像，图像分割和边缘提取都很难达到理想的效果，这样就使得系统的灵活性不高。最后，很多文献都把咬翼片图像作为输入图像，但咬翼片图像的拍摄需要被拍摄人咬住纸板翼片来拍摄，而对于实际应用中的无名尸骨，这种拍摄显然是不可行的。
技术实现思路
针对现有技术的不足，本专利技术提供一种基于SURF和RAN本文档来自技高网...

【技术保护点】
一种基于SURF特征点匹配和RANSAC模型估计的牙齿X‑ray影像匹配方法，其特征在于，该匹配方法包括步骤如下：(1)采集输入图像；在图库中选取一张新的参考图像；(2)将输入图像与参考图像进行灰度线性增强；(3)分别在输入图像和参考图像中检测SURF特征点，并对各个SURF特征点进行SURF特征描述；(4)将输入图像中的SURF特征描述向量与参考图像中的SURF描述向量进行初步匹配；根据初步匹配所产生的SURF匹配样本，估算输入图像与参考图像之间的仿射变换模型；(5)根据所估算出的最佳仿射变换模型，将输入图像变换为与参考图像大小相同的新图像；(6)计算新图像的有效区域与参考图像重合部分的灰度相关系数：若相关系数大于90%，则输入图像与参考图像是匹配的，跳至(8)；否则就是不匹配的，跳至(7)；(7)若图库中还有没被匹配的参考图像，则返回(1)，否则退出；(8)匹配成功，记录下参考图像对应的个人信息，退出。

【技术特征摘要】
1.一种基于SURF特征点匹配和RANSAC模型估计的牙齿X-ray影像匹配方法，其特征在于，该匹配方法包括步骤如下：(1)采集输入图像；在图库中选取一张新的参考图像；(2)将输入图像与参考图像进行灰度线性增强；(3)分别在输入图像和参考图像中检测SURF特征点，并对各个SURF特征点进行SURF特征描述；(4)将输入图像中的SURF特征描述向量与参考图像中的SURF特征描述向量进行初步匹配；根据初步匹配所产生的SURF匹配样本，估算输入图像与参考图像之间的仿射变换模型；(5)根据所估算出的最佳仿射变换模型，将输入图像变换为与参考图像大小相同的新图像；(6)计算新图像的有效区域与参考图像重合部分的灰度相关系数：若灰度相关系数大于90％，则输入图像与参考图像是匹配的，跳至(8)；否则就是不匹配的，跳至(7)；(7)若图库中还有没被匹配的参考图像，则返回(1)，否则退出；(8)匹配成功，记录下参考图像对应的个人信息，退出；所述匹配方法，包括具体步骤如下：所述步骤(1)具体为：(11)图像采集图像的采集包括输入图像的采集及参考图像的采集，所述的输入图像为待核实牙齿身份的牙齿X-ray影像，所述参考图像为现有牙齿X-ray影像数据库中的任一牙齿X-ray影像；所述步骤(2)具体为：(21)图像增强对输入图像及参考图像的增强方法如下：所述的增强方法是灰度线性增强，即按照式(i)所示的通式，对图像中每一像素点的灰度进行修改：在式(i)中，(m，n)是图像的像素坐标，原始灰度值为I(m,n)，变换后的灰度值为I'(m,n)，min(I)为图像I中的最小灰度值，max(I)为图像I中的最大灰度值；所述步骤(3)具体为：(31)特征提取SURF特征点匹配包括SURF特征点检测、SURF描述子生成和SURF描述子匹配；利用所述SURF特征点检测、SURF描述子生成对经步骤(2)处理后的输入图像和参考图像进行特征提取：a.SURF特征点检测首先对输入图像和参考图像进行不同尺度的高斯滤波，如式(ii)所示ImgG(x,y,σ)＝Img(x,y)*G(x,y,σ)，(ii)在式(ii)中，(x，y)为图像的坐标，ImgG(x，y，σ)为高斯滤波后的图像，增加了尺度变量σ；Img(x，y)为原始图像；G(x，y，σ)的定义如式(iii)所示；然后，将滤波后的图像ImgG(x，y，σ)代入Hessian行列式，如式(iv)所示：解得使式(iv)取得局部极值的点(xk，yk，σk)，其中(xk，yk)即为输入图像或参考图像中SURF特征点的坐标，σk为SURF特征点(xk，yk)所对应的特征尺度，这样就完成了SURF特征点检测的步骤；b.SURF描述子生成对每一个特征点取一定大小的邻域：对于SURF特征点(xk，yk)，其邻域设为以SURF特征点(xk，yk)为中心，以15σk为半径的圆域Ck；为了保证所生成SURF描述子的尺度不变性，对圆域Ck的半径进行统一化，即以σk为横纵间隔，对圆域Ck进行下采样；所得到的下采样点组成了一个归一化的圆域Ckn，中心为SURF特征点(xk，yk)，半径为15个像素；然后继续在所述归一化的圆域Ckn内生成SURF特征点(xk，yk)的SURF描述子；为了保证所生成SURF描述子的旋转不变性，对归一化圆域Ckn的梯度主方向进行统一化：首先解出圆域Ckn内每一点的梯度；对圆域Ckn内每一点的两个梯度分量进行高斯加权；高斯函数的中心点为SURF特征点(xk，yk)，尺度σ为3.3；然后根据圆域Ckn内全部加权后的梯度进行主方向计算；计算规则如式(v)所示：在式(v)中，θmain(Ckn)为圆域Ckn内的主方向角度，即SURF特征点(xk，yk)的主方向；(xw，yw)为圆域Ckn内高斯加权后的梯度向量；解得每个SURF特征点的主方向后，将所有的SURF特征点旋转到统一的方向；设旋转后的圆域Ckn为C′kn，仍然以SURF特征点(xk，yk)为中心，以15个像素为半径；为了提取SURF特征点(xk，yk)的SURF描述子向量，首先对圆域C′kn进行高斯加权，以SURF特征点(xk，yk)为中心，以3.3为尺度σ；然后以SURF特征点(xk，yk)为中心，以圆域c′kn的坐标轴为对称轴，作一个边长为20个像素的正方形Rkn，并求取正方形区域Rkn内各点的梯度向量(xw',yw')；再将正方形区域Rkn平均分割成4×4的子区域，每个子区域的大小均为5×5；最后对每个5×5的子区域，统计如式(vi)所示的4维特征向量；式(vi)中，Rkn(i)为正方形区域Rkn中的第i个子区域；由于式(vi)中的特征向量中不包含特征点的位置信息，而仅仅包含特征点邻域内的梯度信息，因此是平移不变的；由于对正方形区域Rkn中的每个子区域都能够提取到如式(vi)所示的4维特征向量，那么对于正方形区域Rkn的4×4个子区域，总共提取到4×4×4＝64维平移不变的特征向量，这64维的向量，就是SURF特征点(xk，yk)的SURF描述子；所述步骤(4)具体为：(41)特征的匹配该步骤包括SURF描述子匹配、仿射变换模型估计、对输入图像进行仿射变换：c.SURF描述子是用来描述SURF特征点的特...

【专利技术属性】
技术研发人员：孙国霞，刘伟，孙惠强，李惠，周舟，
申请(专利权)人：山东大学，
类型：发明
国别省市：山东;37