首页 / 专利库 / 医疗设备 / 出牙环 / 基于灰度直方图的单牙及其牙髓分割方法

基于灰度直方图的单牙及其牙髓分割方法

阅读:434发布:2020-05-11

专利汇可以提供基于灰度直方图的单牙及其牙髓分割方法专利检索,专利查询,专利分析的服务。并且本 发明 提出一种基于灰度直方图的单牙及其 牙髓 分割方法。 牙齿 的CBCT图像中,不同图层灰度值差异较大,难以通过单一 阈值 进行三维牙齿、牙髓分割,但在同一图层中,牙齿灰度较周围组织具有显著特征。本发明对每一图层绘制灰度直方图,通过直方图信息,找到其最适合的灰度进行二值分割,获得牙齿边界所形成的闭环区域,并填充闭环区域内的孔洞,即可快速分割牙齿与周围组织。在得到分割完成的单牙后,运用相似方法即可分割出该牙齿的牙髓。以上单颗牙齿CBCT三维图像牙齿及其牙髓分割及建模方法,实现了牙齿与周围组织、牙齿与牙髓的快速分割。,下面是基于灰度直方图的单牙及其牙髓分割方法专利的具体信息内容。

1.一种单颗牙齿CBCT图像的牙齿牙髓分割方法,包含以下步骤:
第一部分:单个CBCT牙齿图像的分割算法
步骤1.采集单颗牙齿CBCT三维图像序列,作为原始牙齿CT图像,对待处理的单颗牙齿X光片图像进行预处理;
步骤1.1.采集单颗牙齿CBCT三维图像序列,对于该单颗牙齿CBCT三维图像,以牙冠向上、牙根向下为正方向,读入图像序列,该序列内的每一张图层作为分割算法应用截面;
步骤1.2.将步骤1.1中选择的特定方向的二维图像序列灰度值归一化,调整图像灰度到(0,255)区间,得到的二维图像序列灰度图,用于提供给步骤2;
步骤2.遍历步骤1.2中得到的预处理后的二维图像序列中的每一个图层;
步骤2.1.将步骤1.2中得到的预处理后的每一层图像序列绘制对应图片的灰度直方图,得到相应的灰度直方图;
步骤2.2.对步骤2.1的灰度直方图,选择20次多项式对图像进行函数拟合,拟合得到的函数图像;
所述拟合函数为
x:直方图中每个像素的灰度值;
步骤2.3,在步骤2.2拟合得到的函数图像中,选取靠近灰度值最大值一侧(x→200)的第一个波谷对应的灰度值,作为计算得到的波谷值(min_global)。
步骤2.4,设置一个阈值偏移量(de_weight),调整最终分割阈值(final_global),所述最终分割阈值等于由步骤2.3计算得到的波谷值减去阈值偏移量;
Snumber:表示图层在图像序列的位置
a:该图像序列的图层总数min_global;
所述最终分割阈值计算:
final_global=min_global-de_weight
步骤2.5,用步骤2.4得到的最终分割阈值,执行阈值分割算法,将图像二值化,获得牙齿封闭轮廓的初步分割二值图像结果;
步骤2.6,对步骤2.5得到的二值图像,以(0,0)为起始种子点,自动选取背景点进行泛洪填充,得到泛洪填充后的图像(flood_image);
步骤2.7,对步骤2.6得到的泛洪填充后的图像取反,得到牙齿内部为填充图像(reverse_image)。
步骤2.8,将步骤2.5得到的二值图片与步骤2.7得到的牙齿内部为填充图像(reverse_image)做“并”运算得到完整的牙齿分割截面;
步骤2.9,对步骤2.8得到的完整的牙齿分割截面,统计图片内所有连通区域的数量(N),并计算各个面积(Si,i=1,2,3,……,N),计算连通区域的平均面积(Savg),将小于平均面积Savg的区域视为噪点,全部清除,并通过开、闭运算形态学操作;
步骤3,将步骤2.9得到的一层层图像排列起来,构成分割结束后的单牙二值图像序列;
步骤4,对步骤3得到的二值图像序列与原始牙齿CT图像序列求并集,得到最终单颗牙齿CT图像分割结果;
第二部分:牙齿的牙髓图像的分割算法
步骤5,读入步骤4得到的图像序列,作为分割算法应用截面;
步骤6,对步骤5读入的图像序列中每张图,以(0,0)为起始种子点,自动选取背景点进行泛洪填充,将背景灰度值全部填充为255,得到填充后的图像;
步骤7,遍历操作步骤6,得到所有填充后的图像序列;
步骤7.1,将步骤6中得到的预处理后的图像序列绘制每一层图片的灰度直方图,得到灰度直方图;
步骤7.2,对步骤7.1的灰度直方图,选择20次多项式对图像进行函数拟合x:直方图中每个像素的灰度值;
步骤7.3在函数图像中,选取靠近灰度值最大值一侧(x→200)的第一个波谷对应的灰度值,作为计算得到的波谷值(Endo_min_global);
步骤7.4设置一个牙髓阈值偏移量(Endo_de_weight),计算最终牙髓分割阈值(Endo_final),最终牙髓分割阈值等于由步骤7.3计算得到的波谷值减去牙髓阈值偏移量Endo_de_weight;
Snumber:图层在截面的位置,
a:该截面的图层总数Endo_min_global;
计算最终牙髓分割阈值:Endo_final=Endo_min_global–Endo_de_weight;
步骤7.5.采用步骤7.4得到的最终牙髓分割阈值,执行阈值分割算法,将图像二值化,获得牙髓闭轮廓的初步分割二值图像分割结果;
步骤7.6对步骤7.5得到的牙髓分割截面,统计图片内所有连通区域的数量(N),并计算各个面积(Si,i=1,2,3,……,N),然后计算连通区域的平均面积(Savg),将小于平均面积的区域视为噪点,全部清除,再通过开、闭运算形态学操作,
步骤8,将步骤7.6得到的一层层图像排列起来,构成分割结束后的单牙牙髓二值图像序列;
步骤9,对步骤8得到的二值图像序列与步骤4所得的单颗牙牙齿CT图像序列求并集,得到最终单颗牙齿牙髓CT图像分割结果。

说明书全文

基于灰度直方图的单牙及其牙髓分割方法

技术领域

[0001] 本发明属于数字口腔图像分割领域,公开了针对单颗牙齿CBCT三维图像的牙齿与牙髓分割方法。

背景技术

[0002] 牙齿CBCT的分割在临床的研究具有重要意义:已有文献表明在牙髓治疗前了解患牙的解剖形态有利于治疗效果;牙齿的根管形态错综复杂,通过牙齿牙髓形态的分割重建可用于形态学研究。
[0003] 目前牙齿CBCT分割的难点
[0004] 由于牙科CBCT放射剂量小,时间短,导致图像噪声干扰多,牙齿本身的成像不均匀,图像灰度不统一,且牙槽骨牙根密度接近使得图像中牙齿界限模糊,较难分辨。

发明内容

[0005] 为了解决背景技术中存在的问题,本发明提出一种单颗牙齿CBCT三维图像牙齿及其牙髓分割及建模方法,实现了牙齿与周围组织、牙齿与牙髓的快速分割。
[0006] 本发明需要保护的技术方案为:
[0007] 一种单颗牙齿CBCT图像的牙齿牙髓分割方法,包含以下步骤:
[0008] 第一部分:单个CBCT牙齿图像的分割算法
[0009] 步骤1.采集单颗牙齿CBCT三维图像序列,作为原始牙齿CT图像,对待处理的单颗牙齿X光片图像进行预处理;
[0010] 步骤1.1.采集单颗牙齿CBCT三维图像序列,对于该单颗牙齿CBCT三维图像,以牙冠向上、牙根向下为正方向,读入图像序列,该序列内的每一张图层作为分割算法应用截面;
[0011] 步骤1.2.将步骤1.1中选择的特定方向的二维图像序列灰度值归一化,调整图像灰度到(0,255)区间,得到的二维图像序列灰度图,用于提供给步骤2,便于后续灰度分析;
[0012] 步骤2.遍历步骤1.2中得到的预处理后的二维图像序列中的每一个图层。
[0013] 步骤2.1.将步骤1.2中得到的预处理后的每一层图像序列绘制对应图片的灰度直方图,得到相应的灰度直方图;
[0014] 步骤2.2.对步骤2.1的灰度直方图,选择20次多项式对图像进行函数拟合,拟合得到的函数图像;
[0015] 所述拟合函数为
[0016] x:直方图中每个像素的灰度值
[0017] 由于牙齿灰度大多集中在(100,200)区间,为提高函数拟合效果,选取该区间进行函数拟合。
[0018] 步骤2.3,在步骤2.2拟合得到的函数图像中,选取靠近灰度值最大值一侧(x→200)的第一个波谷对应的灰度值,作为计算得到的波谷值(min_global)。
[0019] 步骤2.4,设置一个阈值偏移量(de_weight),调整最终分割阈值(final_global),所述最终分割阈值等于由步骤2.3计算得到的波谷值减去阈值偏移量。
[0020]
[0021] Snumber:表示图层在图像序列的位置
[0022] a:该图像序列的图层总数min_global;
[0023] 所述最终分割阈值计算:
[0024] final_global=min_global-de_weight
[0025] 步骤2.5,用步骤2.4得到的最终分割阈值,执行阈值分割算法(该算法过程已为现有技术),将图像二值化,获得牙齿封闭轮廓的初步分割二值图像结果(th2_image)。
[0026] 步骤2.6,对步骤2.5得到的二值图像,以(0,0)为起始种子点,自动选取背景点进行泛洪填充,得到泛洪填充后的图像(flood_image),由于牙齿外轮廓为封闭图形,此时未填充的区域仅剩牙齿内部。
[0027] 步骤2.7,对步骤2.6得到的泛洪填充后的图像取反,得到牙齿内部为填充图像(reverse_image)。
[0028] 步骤2.8,将步骤2.5得到的二值图片(th2_image)与步骤2.7得到的牙齿内部为填充图像(reverse_image)做“并”运算(th2_image or reverse_image)得到完整的牙齿分割截面。
[0029] 步骤2.9,对步骤2.8得到的完整的牙齿分割截面,统计图片内所有连通区域的数量(N),并计算各个面积(Si,i=1,2,3,……,N),计算连通区域的平均面积(Savg),将小于平均面积Savg的区域视为噪点,全部清除,并通过开、闭运算等形态学操作,进一步优化分割结果。
[0030]
[0031] 步骤3,将步骤2.9得到的一层层图像排列起来,构成分割结束后的单牙二值图像序列。
[0032] 步骤4,对步骤3得到的二值图像序列与原始牙齿CT图像序列求并集,得到最终单颗牙齿CT图像分割结果。
[0033] 第二部分:牙齿的牙髓图像的分割算法
[0034] 步骤5,读入步骤4得到的图像序列,作为分割算法应用截面;
[0035] 步骤6,对步骤5读入的图像序列中每张图,以(0,0)为起始种子点,自动选取背景点进行泛洪填充,将背景灰度值全部填充为255,得到填充后的图像。
[0036] 步骤7,遍历操作步骤6,得到所有填充后的图像序列。
[0037] 步骤7.1,将步骤6中得到的预处理后的图像序列绘制每一层图片的灰度直方图,得到灰度直方图;
[0038] 步骤7.2,对步骤7.1的灰度直方图,选择20次多项式对图像进行函数拟合[0039] x:直方图中每个像素的灰度值
[0040] 由于牙齿灰度大多集中在(100,200)区间,为提高函数拟合效果,选取该区间进行函数拟合。
[0041] 步骤7.3在函数图像中,选取靠近灰度值最大值一侧(x→200)的第一个波谷对应的灰度值,作为计算得到的波谷值(Endo_min_global)。
[0042] 步骤7.4设置一个牙髓阈值偏移量(Endo_de_weight),计算最终牙髓分割阈值(Endo_final),最终牙髓分割阈值等于由步骤7.3计算得到的波谷值减去牙髓阈值偏移量Endo_de_weight;
[0043] Snumber:图层在截面的位置,
[0044] a:该截面的图层总数Endo_min_global;
[0045] 计算最终牙髓分割阈值:Endo_final=Endo_min_global–Endo_de_weight。
[0046] 步骤7.5.采用步骤7.4得到的最终牙髓分割阈值,执行阈值分割算法(为现有算法),将图像二值化,获得牙髓闭轮廓的初步分割二值图像分割结果。
[0047] 步骤7.6对步骤7.5得到的牙髓分割截面,统计图片内所有连通区域的数量(N),并计算各个面积(Si,i=1,2,3,……,N),然后计算连通区域的平均面积(Savg),将小于平均面积的区域视为噪点,全部清除,再通过开、闭运算等形态学操作,进一步优化分割结果。
[0048]
[0049] 步骤8,将步骤7.6得到的一层层图像排列起来,构成分割结束后的单牙牙髓二值图像序列。
[0050] 步骤9,对步骤8得到的二值图像序列与步骤4所得的单颗牙牙齿CT图像序列求并集,得到最终单颗牙齿牙髓CT图像分割结果。
[0051] 牙齿的CBCT图像中,不同图层灰度值差异较大,难以通过单一阈值进行三维牙齿、牙髓分割,但在同一图层中,牙齿灰度较周围组织具有显著特征。本发明对每一图层绘制灰度直方图,通过直方图信息,找到其最适合的灰度进行二值分割,获得牙齿边界所形成的闭环区域,并填充闭环区域内的孔洞,即可快速分割牙齿与周围组织。在得到分割完成的单牙后,运用相似方法即可分割出该牙齿的牙髓。本发明通过一系列操作,处理由于分割产生的细小噪点,从而改善分割效果。
[0052] 本发明的有益效果是:
[0053] 本发明基于灰度进行图像分割,计算量小,能快速准确的从CBCT数据中分割出完整的整颗牙齿,在此基础上运用相似方法能快速分离重建此颗牙齿的牙髓,有利于牙髓治疗前了解患牙的解剖形态从而提高治疗效果。由于本发明在图像中的分割是自动的,故可大大提高效率。附图说明
[0054] 图1为单个CBCT牙齿图像的分割算法流程图
[0055] 图2为在图1完成单颗牙齿分割后,该牙齿的牙髓图像的分割算法流程图。
[0056] 图3a为本发明算法应用图像截面(待分割图像序列的其中一张截图),图3b、图3c是单颗牙齿三维图像序列中的另外两个截面,可以看出图3a所示截面牙齿灰度具有统一而明显的特征,而其他两个截面牙齿灰度不统一。
[0057] 图4a、图4b、图4c显示了绘制灰度直方图拟合函数的过程,其中图4a是待分割的牙齿截面图,图4b是是该图像的灰度直方图,图4c是在灰度直方图中选取灰度值在(100,200)区间,用多项式拟合出来的函数图像,所需的特定波谷在图中已标出。
[0058] 图5对比了是否设置阈值偏移量对牙冠部分的阈值分割效果的影响。图5a为待分割的牙齿截面图;图5b为该图的灰度直方图在(100,200)区间拟合函数图像,所需的特定波谷已标出;图5c为阈值设置为图5b中波谷对应灰度值的分割情况;图5d为阈值设置为图5b中波谷对应灰度值加上阈值偏移量后的分割情况。可以看出对于牙冠部分设置阈值偏移量后,能得到封闭牙齿轮廓,保证下一步泛洪填充达到理想效果。
[0059] 图6是本发明中对初步分割后牙齿截面图进行牙齿轮廓内的孔洞填充的过程,目的是获得完整的牙齿分割区域。图6a为阈值分割后的图像,图6b为泛洪填充后的图像,图6c为泛洪填充后取反的图像,图6d为图6a与图6c求并集的图像,即完成牙齿内部孔洞填充的图像。
[0060] 图7对上述完整的牙齿分割截面进行去噪处理,图7a为去噪前图像,图7b为去噪后图像,
[0061] 图8为二值图像序列与原始牙齿图像序列求并集,获得的单颗牙齿图像最终分割结果(以序列中两个截面为例)。
[0062] 图9为在分割出单颗牙齿的基础上,运用类似方法分割其牙髓的结果。图9a为分割后牙齿截面图像,图9b为该牙齿的牙髓截面分割结果图像。具体实施方案
[0063] 下面结合附图,进一步说明本发明的技术方案。
[0064] 按照本发明方法的实施例及其具体实施过程如下:
[0065] CBCT图像不同层牙齿的灰度值差异较大,难以通过单一阈值进行整个三维图像的牙齿、牙髓分割。但以牙冠向上、牙根向下为正方向的这一图像序列中,每一层牙齿灰度特征较为明显且统一,故设计此算法针对每一层图像进行自动分割,从而获取三维分割结果。
[0066] 图3为CBCT三维图像的三个视,本实施例选择图3a所示的视角作为算法实施截面,该截面中,牙齿区域较其他部分有统一且明显的灰度特征,其灰度值更大。
[0067] 第一部分
[0068] 选定视角之后,在CT图像序列中,规定图层序列号小值对应牙根,图层序列号大值对应牙冠,对该序列每张图片进行以下操作:
[0069] (11)将图像灰度值归一化到(0,255)区间便于后续处理
[0070] (12)以图4a为例,绘制其灰度直方图(图4b),选择灰度值在(100,200)区间,用20次多项式f(x)进行拟合,得到图4c。图4c中x轴大值一侧(x→200)的第一个波谷对应的灰度值(如箭头所示),即为计算所得的分割阈值。
[0071]
[0072] (13)设置一个阈值偏移量de_weight,调整最终分割阈值(final_global)。最终分割阈值等于由(12)计算得到的波谷值减去阈值偏移量。
[0073]
[0074] Snumber:图层在截面的位置a:该截面的图层总数min_global:计算得到的波谷值[0075] final_global=min_global-de_weight
[0076] 设置阈值偏移量的目的为:若用计算所得的阈值分割图像,容易产生不封闭区域(如图5c),不利接下来的填充孔洞操作,因此加入阈值偏移量提高分割效果(如图5d)。
[0077] (14)使用最终阈值对图像进行二值分割,得到图6a,接下来需要对牙齿内部进行孔洞填充。在待填充的图层中,以(0,0)为起始种子点,自动选取背景点进行泛洪填充,得到图6b,对泛洪填充后的图像取反,得到牙齿内部为填充图像图6c。将外轮廓图像与牙齿内部为填充图像做“并(or)”运算,得到完整的二值牙齿分割截面图6d。
[0078] (15)对上述完整的二值牙齿分割截面进行去噪处理。统计图片内所有连通区域的数量,并计算各个面积。计算连通区域的平均面积,将小于平均面积的区域视为噪点,全部清除,通过开、闭运算等形态学操作,进一步优化分割结果(对比图7a,图7b)。
[0079] 完成上述操作之后,得到一个完整的牙齿三维二值数组,将此数组与原始三维图像数组求交集,即可获得最终分割结果。(其中选取两个截面作为示例,如图8)。
[0080] 第二部分接下来用相似思路完成牙髓分割。对于已分割的单颗CBCT牙齿图像,选择上述序列作为分割算法应用截面,对序列中每一张图片进行如下操作:
[0081] (21)对二维图像序列中每张图,以(0,0)为起始种子点,自动选取背景点进行泛洪填充,将背景灰度值全部填充为255,得到填充后的图像。
[0082] (22)重复(12)中所有步骤,得到自动计算的牙髓特定波谷值(Endo_min_global)。
[0083] (23)设置一个阈值偏移量(Endo_de_weight),计算最终牙髓分割阈值(Endo_final),最终牙髓分割阈值等于由(22)计算得到的波谷值减去牙髓阈值偏移量。
[0084]
[0085] Snumber:图层在截面的位置a:该截面的图层总数Endo_min_global:
[0086] 计算得到的波谷值
[0087] Endo_final=Endo_min_global–Endo_de_weigh(24)执行阈值分割算法,将图像二值化,获得牙髓分割结果。
[0088] (25)对上述完整的牙髓分割截面序列,用(15)中所阐述的方法进行去噪处理,通过开、闭运算等形态学操作优化分割结果,得到最终该颗牙齿的牙髓图像(图9)。
[0089] 本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举,本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式,本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈