自动采样厚度检测的非线性口腔CT全景图像合成方法技术

自动采样厚度检测的非线性口腔CT全景图像合成方法,包括15个步骤。本发明专利技术的有益效果如下：第一、前期不需要人工设定断层图像范围。第二、能够自动计算和识别包含上下颌的断层图像范围，实现真正的全自动牙弓曲线提取。第三、在提取牙弓曲线的同时自动检测虚拟全景的采样厚度，克服人工设置过大或过小导致的虚拟全景合成存在的模糊或不全的问题，实现真正的全自动提取过程。本发明专利技术和现有方法相比，极大的提高的虚拟全景图像的清晰度和对比度，能够获得和真实X射线全景图像相近似的效果。

【技术实现步骤摘要】
自动采样厚度检测的非线性口腔CT全景图像合成方法
本专利技术属于口腔CT重建数据的图像后处理
，特别涉及自动采样厚度检测的非线性口腔CT全景图像合成方法。
技术介绍
头部CT重建获取的断层图像所生成的虚拟全景图像是口腔CT图像后处理的一项重要功能，虚拟全景图像在口腔手术及口腔种植有广泛应用，虚拟全景图像在牙齿正畸、压槽骨及颞颌关节等疾病的诊断和治疗方面提供了上下颌的全局信息，精确牙弓线的生成是牙科全景图像生成必不可少的一个重要步骤，直接决定了虚拟全景图像生成的正确性和最终效果。现有技术中，牙弓曲线检测方法可以分为直接方法和间接方法两类。直接方法是依次对输入的断层图像进行牙齿分割和形态学操作，获取所有断层的牙弓曲线，并从所有的牙弓曲线中选择最优的结果为最后的牙弓曲线。间接方法首先生成体数据不同角度的投影图像，在投影图像上经过分割和形态学操作获取牙弓曲线。直接方法可实现全自动牙弓曲线提取，但由于需要对所有输入断层图像进行操作，其计算量较大，运行时间较长，无法达到临床实际的需求。另外由于直接方法是从多个断层提取的牙弓曲线中获取最优结果，其本质是用局部来代表全局牙弓曲线，操作较大的偏差。间接方法具有速度快和直接表达全局牙弓曲线的优势，但是间接方法无法实现全自动的提取，需要人工设置断层图像的范围。直接方法和间接方法都存在的一个主要问题是虚拟全景图像的合成都采用ray-sum方法，如图10所示，合成的虚拟全景图像非常模糊，对比度差，直接影响医生对虚拟全景图像的观察。现有技术中的牙弓曲线检测方法还有一个缺点，是需要人工设置虚拟全景的采样厚度，设置过大则虚拟全景图像会变的模糊；反之虚拟全景图像无法包含所有的上下颌解决结构信息。因此针对现有技术不足，提供一种自动采样厚度检测的非线性口腔CT全景图像合成方法以解决现有技术不足甚为必要。
技术实现思路
本专利技术的目的在于避免现有技术的不足之处而提供一种自动采样厚度检测的非线性口腔CT全景图像合成方法。该自动采样厚度检测的非线性口腔CT全景图像合成方法能够极大的提高的虚拟全景图像的清晰度和对比度，能够获得和真实X射线全景图像相近似的效果。本专利技术的上述目的通过以下技术措施实现：提供一种自动采样厚度检测的非线性口腔CT全景图像合成方法,步骤包括有：步骤一、输入口腔CT重建后体数据，得到体数据最大密度投影图像；步骤二、对步骤一得到的体数据最大密度投影图像进行二维高斯滤波，得到滤波体数据最大密度投影图像；步骤三、计算步骤二得到的滤波体数据最大密度投影图像的直方图，得到滤波体数据最大密度投影图像分割阈值；步骤四、用步骤三得到的滤波体数据最大密度投影图像分割阈值对步骤二得到的滤波体数据最大密度投影图像进行分割，得到上下颌中牙齿在序列图像中范围；步骤五、由步骤四得到的上下颌中牙齿在序列图像中范围生成得到横断面最大密度投影图像；步骤六、统计步骤五得到的横断面最大密度投影图像的直方图，得到横断面最大密度投影图像分割阈值；步骤七、根据步骤六得到的横断面最大密度投影图像分割阈值对步骤五得到的横断面最大密度投影图像进行分割，得到横断面最大密度投影图像分割结果；步骤八、对步骤七得到的横断面最大密度投影图像分割结果进行轮廓统计，然后滤波，得到牙弓轮廓；步骤九、对步骤八得到的牙弓轮廓进行细化操作，得到初步牙弓曲线；步骤十、统计步骤九得到的初步牙弓曲线的端点，得到两端点初步牙弓曲线；步骤十一、对步骤十得到的两端点初步牙弓曲线进行控制点采样操作，得到初步牙弓曲线进控制点；步骤十二、根据步骤十一得到的初步牙弓曲线进控制点构建B-Cubic样条曲线，得到正式牙弓曲线；步骤十三、根据步骤八得到的牙弓轮廓和步骤十二得到的正式牙弓曲线计算，得到虚拟全景图像采样厚度；步骤十四、对步骤十二得到的正式牙弓曲线进行等像素间隔采样，并计算采样所在点的法向量，并根据步骤十三得到的虚拟全景图像采样厚度对步骤一的口腔CT重建后体数据进行切割，得到全景容积图像和虚拟全景容积数据；步骤十五、沿步骤十四的正式牙弓曲线的法向量的方向对步骤十四的虚拟全景容积数据和全景容积图像进行合成，得到虚拟全景图像。优选的，上述体数据最大密度投影图像为体数据冠状面最大密度投影图像或者体数据矢状面的最大密度投影图像。优选的，上述步骤三具体为，步骤3.1、计算步骤二得到的滤波体数据最大密度投影图像的直方图，进入步骤3.2；步骤3.2、计算直方图的峰值，进入步骤3.3；步骤3.3、从大至小排列直方图的灰度值得到第一个直方图峰值，进入步骤3.4；步骤3.4、采用高斯函数逼近第一个直方图峰值，取右侧置信区间98％的值，得到滤波体数据最大密度投影图像分割阈值；优选的，上述步骤六具体为；步骤6.1、统计步骤五得到的横断面最大密度投影图像的直方图，进入步骤6.2；步骤6.2、计算的直方图的峰值，进入步骤6.3；步骤6.3、从大至小排列直方图的灰度值，得到第一个直方图峰值并定义为横断面最大密度投影图像分割阈值；优选的，上述步骤八具体为，步骤8.1、对步骤七得到的横断面最大密度投影图像分割结果进行轮廓统计，得到轮廓面积，进入步骤8.2；步骤8.2、提取面积最大的轮廓进行横坐标和纵坐标的一维高斯滤波，得到滤波后轮廓，进入步骤8.3；步骤8.3、对滤波后轮廓进行绘制，得到牙弓轮廓。优选的，上述步骤十具体为，步骤10.1、统计步骤九得到的初步牙弓曲线的端点的个数，进入步骤10.3，步骤10.2、统计去除分叉后初步牙弓曲线的端点的个数，进入步骤10.3，步骤10.3、如果端点的个数多于两个，则进入步骤10.4，如何端点个数等于两个，则进入步骤10.5；步骤10.4、对步骤九得到的初步牙弓曲线进行去除分叉的形态学操作，得到去除分叉后初步牙弓曲线，进入步骤10.2；步骤10.5、定义为两端点初步牙弓曲线，进入步骤十一。优选的，上述步骤十一具体为，对步骤十得到的两端点初步牙弓曲线进行控制点采样操作，等间隔或选取M个控制点，得到初步牙弓曲线进控制点，且M≥3，M为正整数。优选的，上述M＝11。优选的，上述步骤十三具体为，步骤13.1、对步骤十二得到的正式牙弓曲线进行等间距采集，得到等间距采样点，进入步骤13.2；步骤13.2、在步骤13.1的等间距采样点处，计算正式牙弓曲线的法向量，进入步骤13.3；步骤13.3、沿步骤13.3的法向量方向，计算步骤八得到的牙弓轮廓边缘的距离，得到牙弓轮廓边缘距离，进入步骤13.4；步骤13.4、统计步骤13.3得到牙弓轮廓边缘距离的直方图，得到距离直方图，进入步骤13.5；步骤13.5、计算步骤13.4得到的距离直方图的最大峰值，并取最大峰值的±20范围的平均值，得到虚拟全景图像采样厚度。优选的，上述步骤十四具体为，步骤14.1、对步骤十二得到的正式牙弓曲线进行等像素间隔采样，得到等像素间隔采样点，进入步骤14.2；步骤14.2、并计算等像素间隔采样点所在位置的法向量，沿法向量方向，并根据步骤十三得到的虚拟全景图像采样厚度对步骤一的口腔CT重建后体数据进行切割，进入步骤14.3；步骤14.3、根据步骤14.1得到的等像素间隔采样点和步骤14.2得到的切割图，依次得到所有等像素间隔采样点的进入切割，每个等像素间隔采样点各得到全景容积图像，进入步骤14本文档来自技高网...

【技术保护点】
1.自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于，步骤包括有：步骤一、输入口腔CT重建后体数据，得到体数据最大密度投影图像；步骤二、对步骤一得到的体数据最大密度投影图像进行二维高斯滤波，得到滤波体数据最大密度投影图像；步骤三、计算步骤二得到的滤波体数据最大密度投影图像的直方图，得到滤波体数据最大密度投影图像分割阈值；步骤四、用步骤三得到的滤波体数据最大密度投影图像分割阈值对步骤二得到的滤波体数据最大密度投影图像进行分割，得到上下颌中牙齿在序列图像中范围；步骤五、由步骤四得到的上下颌中牙齿在序列图像中范围生成得到横断面最大密度投影图像；步骤六、统计步骤五得到的横断面最大密度投影图像的直方图，得到横断面最大密度投影图像分割阈值；步骤七、根据步骤六得到的横断面最大密度投影图像分割阈值对步骤五得到的横断面最大密度投影图像进行分割，得到横断面最大密度投影图像分割结果；步骤八、对步骤七得到的横断面最大密度投影图像分割结果进行轮廓统计，然后滤波，得到牙弓轮廓；步骤九、对步骤八得到的牙弓轮廓进行细化操作，得到初步牙弓曲线；步骤十、统计步骤九得到的初步牙弓曲线的端点，得到两端点初步牙弓曲线；步骤十一、对步骤十得到的两端点初步牙弓曲线进行控制点采样操作，得到初步牙弓曲线进控制点；步骤十二、根据步骤十一得到的初步牙弓曲线进控制点构建B‑Cubic样条曲线，得到正式牙弓曲线；步骤十三、根据步骤八得到的牙弓轮廓和步骤十二得到的正式牙弓曲线计算，得到虚拟全景图像采样厚度；步骤十四、对步骤十二得到的正式牙弓曲线进行等像素间隔采样，并计算采样所在点的法向量，并根据步骤十三得到的虚拟全景图像采样厚度对步骤一的口腔CT重建后体数据进行切割，得到全景容积图像和虚拟全景容积数据；步骤十五、沿步骤十四的正式牙弓曲线的法向量的方向对步骤十四的虚拟全景容积数据和全景容积图像进行合成，得到虚拟全景图像。...

【技术特征摘要】
1.自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于，步骤包括有：步骤一、输入口腔CT重建后体数据，得到体数据最大密度投影图像；步骤二、对步骤一得到的体数据最大密度投影图像进行二维高斯滤波，得到滤波体数据最大密度投影图像；步骤三、计算步骤二得到的滤波体数据最大密度投影图像的直方图，得到滤波体数据最大密度投影图像分割阈值；步骤四、用步骤三得到的滤波体数据最大密度投影图像分割阈值对步骤二得到的滤波体数据最大密度投影图像进行分割，得到上下颌中牙齿在序列图像中范围；步骤五、由步骤四得到的上下颌中牙齿在序列图像中范围生成得到横断面最大密度投影图像；步骤六、统计步骤五得到的横断面最大密度投影图像的直方图，得到横断面最大密度投影图像分割阈值；步骤七、根据步骤六得到的横断面最大密度投影图像分割阈值对步骤五得到的横断面最大密度投影图像进行分割，得到横断面最大密度投影图像分割结果；步骤八、对步骤七得到的横断面最大密度投影图像分割结果进行轮廓统计，然后滤波，得到牙弓轮廓；步骤九、对步骤八得到的牙弓轮廓进行细化操作，得到初步牙弓曲线；步骤十、统计步骤九得到的初步牙弓曲线的端点，得到两端点初步牙弓曲线；步骤十一、对步骤十得到的两端点初步牙弓曲线进行控制点采样操作，得到初步牙弓曲线进控制点；步骤十二、根据步骤十一得到的初步牙弓曲线进控制点构建B-Cubic样条曲线，得到正式牙弓曲线；步骤十三、根据步骤八得到的牙弓轮廓和步骤十二得到的正式牙弓曲线计算，得到虚拟全景图像采样厚度；步骤十四、对步骤十二得到的正式牙弓曲线进行等像素间隔采样，并计算采样所在点的法向量，并根据步骤十三得到的虚拟全景图像采样厚度对步骤一的口腔CT重建后体数据进行切割，得到全景容积图像和虚拟全景容积数据；步骤十五、沿步骤十四的正式牙弓曲线的法向量的方向对步骤十四的虚拟全景容积数据和全景容积图像进行合成，得到虚拟全景图像。2.根据权利要求1所述的自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于：所述体数据最大密度投影图像为体数据冠状面最大密度投影图像或者体数据矢状面的最大密度投影图像。3.根据权利要求2所述的自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于：所述步骤三具体为，步骤3.1、计算步骤二得到的滤波体数据最大密度投影图像的直方图，进入步骤3.2；步骤3.2、计算直方图的峰值，进入步骤3.3；步骤3.3、从大至小排列直方图的灰度值得到第一个直方图峰值，进入步骤3.4；步骤3.4、采用高斯函数逼近第一个直方图峰值，取右侧置信区间98％的值，得到滤波体数据最大密度投影图像分割阈值。4.根据权利要求3所述的自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于：所述步骤六具体为；步骤6.1、统计步骤五得到的横断面最大密度投影图像的直方图，进入步骤6.2；步骤6.2、计算的直方图的峰值，进入步骤6.3；步骤6.3、从大至小排列直方图的灰度值，得到第一个直方图峰值并定义为横断面最大密度投影图像分割阈值。5.根据权利要求4所述的自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于：所述步骤八具体为，步骤8.1、对步骤七得到的横断面最大密度投影图像分割结果进行轮廓统计，得到轮廓面积，进入步骤8.2；步骤8.2、提取面积最大的轮廓进行横坐标和纵坐标的一维高斯滤波，得到...

【专利技术属性】
技术研发人员：贠照强，赵磊，杨烁，阳维，冯前进，
申请(专利权)人：南方医科大学，
类型：发明
国别省市：广东,44