本发明属于图像处理和口腔医学技术领域，涉及一种基于牙齿舌(腭)侧的全景图的生成方法。
背景技术：
传统的牙齿检查需要用户前往诊所。牙科医生需要用口腔内窥镜观察每一刻牙齿。另一种方式为，在护士的协助下完成对患者口腔全景图的拍摄。这种拍摄方法成本高(时间成本和费用成本)，效率较低，患者体验不好。对于牙齿的拍摄通常也分为咬合面，唇(颊)面和舌(腭)面三个部分。舌(腭)侧在人们日常刷牙时难以刷到，也是烟斑或牙结石存在的区域。因此，对舌(腭)侧全景图像的拍摄对于牙科医生具有重要的实用价值。
技术实现要素：
本发明的目的在于提供基于局部优化算法的牙齿舌腭侧全景图像拼接方法，本发明的有益效果是方法简单，易于操作，稳定性高。牙齿舌(腭)侧是烟斑和牙结石存在区域。采用这个方法获取的全景图有良好色彩和形态的保真度，利于牙医的口腔检查。
本发明所采用的技术方案是按照以下步骤进行：
1.牙齿舌(腭)侧的全景图像是由微距摄像头进行采集。摄像头固定在预设于牙齿舌(腭)侧的附近的轨道上，拍摄时摄像头沿轨道滑动。在拍摄过程中，摄像头光轴与牙齿舌(腭)侧所在曲面的法线方向所成角度保持近似稳定。
2.所获取的牙齿舌(腭)侧视频按照其帧率保存为离散图片。在平面直角坐标系和平面极坐标系中对离散图片进行筛选。筛选原则是，设立关联性阈值，删除关联度高于所设阈值的图像。
3.根据拍摄视频和模型形状对于牙齿舌(腭)侧所形成柱面的半径进行估计。
4.根据第3点中，估计得到的半径对每一张图片进行柱面投影。
5.利用自适应sift算法，计算相邻两张舌(腭)侧投影图像的特征点坐标，及其对应描述子。并根据这些信息进行特征点匹配。匹配过程中，利用特征点坐标信息剔除靠的过近的匹配特征点对。
6.在两张舌(腭)侧投影图像之间，根据其特征点的匹配情况，利用局部(两张图间)光束平差法进行图像变换参数的计算(透视变换参数)。
7.将得到得到两张舌(腭)侧投影图像按照计算得到的参数拼接为一张图像。
8.以得到的这张局部舌(腭)侧局部全景图为基础，重复步骤5－7，依次将后续图像拼接到局部图上，直到所有图像都完成拼接。最后形成牙齿舌(腭)侧的全景图。
进一步，阈值选取通过以下过程实现:a.给定一组连续拍摄舌侧图像，任意选择某一图片作为基准图片，计算相邻图像与它的关联值大小；b.从按照关联值的大小，由高到低将对应的图像与基准图像两两拼接，直到拼接图像发生明显畸变或错误为止，这时对应的关联值标记为待选阈值；c.重复以上过程若干次，得到若干个待选阈值，选择最大值作为算法的阈值。
附图说明
图1是本发明方法流程示意图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
本发明如图1所示，第一步：舌(腭)侧的全景图像是通过微距摄像头采集。摄像头固定在预设于舌(腭)侧附近的轨道上，拍摄时摄像头沿轨道滑动。在滑动过程中，应保证摄像头稳定，不出现剧烈的视场变化。稳定拍摄是指：摄像头光轴与牙齿舌(腭)侧所在平面所成角度近似不变。如出现不稳定抖动(例如：因为患者牙齿形态畸形而造成的不稳定抖动)，其抖动数据应被惯性传感器捕获，并反馈至算法中。
第二步：所获取的牙齿舌(腭)侧视频按照其帧率保存为离散图片。在平面直角坐标系和平面极坐标系中，对离散图片进行筛选。具体做法是：事先设定好在这两个坐标系下设定阈值。利用相位关联算法计算这两张图在平面直角坐标系和平面极坐标系中相位关联的最大值。在相同坐标系中，将得到的相位关联最大值与其阈值相比。当有任意一个值大于阈值时，则删除序号更大的图片。
第三步：根据轨道本身设计参数，拍摄到的图像及其当前位置对应的柱面曲率半径。
第四步：根据第三步获得曲率半径，对拍摄到的图像进行柱面投影。
第五步：利用自适应sift算法，计算相邻两舌(腭)侧投影图像(分别标记为第i张和第i+1张)的特征点坐标及其对应描述子。并根据这些信息进行特征点匹配。匹配的原则是，计算每一个特征点与其它特征点之间在描述子对应向量的欧式距离，选择距离最小的一个作为匹配的特征点。此外，当匹配完毕，还需要筛除同一张舌(腭)侧图像中空间位置靠的比较近的特征点匹配对。这一过程通过设定距离阈值来完成,距离阈值通常为8个像素。
第六步：在两张舌(腭)侧投影图像之间，根据其特征点的匹配情况，利用局部(两张图间)光束平差法进行图像变换参数的计算(透视变换参数)。具体可以分成两个步骤：第一、利用ransac算法获取一个粗匹配参数，与此同时获得这个参数对应的具有相似几何一致性的特征点集；第二、根据这些筛选出来的点，利用光束平差法完成两张图像间的精匹配过程。
第七步：利用第四步得到的匹配参数对两张舌(腭)侧投影图像进行拼接。拼接的流程是：首先，生成与图像尺寸等大的拼接模版，模版中每个像素的大小由这个像素的位置距离中心位置的距离决定；其次，将第i+1张图及其对应的模版按照第四步中得到的变换参数映射到将i张图及其对应的坐标空间中；比较两幅模版中每个像素的数值大小，保留大的那一个。最后，按照模版的数值对舌(腭)侧投影图像中像素进行加权求和，得到到将i张和到将i+1张的全景图mi，i+1。
第八步：对第i+1和i+2张图重复以上第5－7步。将i+2张图映射到全景图mi，i+1中，并进行粘合，获得到全景图mi，i+2.重复这一过程获得所有图像的全景图。
本发明提出的基于口腔全景图像的牙检报告生成方法能有效简化舌(腭)侧拍摄的过程，让患者在家中可以完成舌(腭)侧全景图像的拍摄和拼接。这简化了获得舌(腭)侧全景图的获取成本，同时也提高了牙医的诊断效率。
以上所述仅是对本发明的较佳实施方式而已，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施方式所做的任何简单修改，等同变化与修饰，均属于本发明技术方案的范围内。