一种Micro-CT图像分析量化评估磨牙牙槽骨的方法及装置制造方法
【专利摘要】本发明公开一种Micro-CT图像分析量化评估磨牙牙槽骨的方法，包括以下步骤：获取包含待分析磨牙的牙槽骨与牙齿的Micro-CT图像；将牙齿图像与牙槽骨图像分离，并将牙齿图像进行3D运算，得到牙齿的3D图像；用球形包围的方法确定感兴趣区域：以待分析磨牙的根分叉中心点到近中根根尖孔中心点距离的平均值，作为球形包围的球形半径r，调整球形包围的中心至待分析磨牙的根分叉中心处，将2D的球形包围区域3D化，最终将感兴趣区域及感兴趣区域内的牙槽骨图像进行3D运算，得到感兴趣区域的总体积、感兴趣区域内牙槽骨的骨小梁体积、骨表面积，并据此计算出感兴趣区域内牙槽骨的骨小梁数目、骨小梁厚度、骨小梁分离度和结构模型指数。
【专利说明】-种Micro-CT图像分析量化评估磨牙牙槽骨的方法及装 CP3
【技术领域】
[0001] 本发明涉及计算机辅助评估牙槽骨领域，具体而言，涉及一种Micro-CT图像分析 量化评估磨牙牙周炎余留牙槽骨的方法及装置。
【背景技术】
[0002] 牙槽骨（alveolar bone)也称为牙槽突（alveolar process)，是上下颌骨包围和 支持牙根的部分。牙槽骨的吸收是牙周炎的主要表现和诊断标准，牙槽骨吸收主要表现在 牙槽骨高度的降低和牙槽骨密度减少。由于牙槽骨内、外骨板均由骨密质构成，中间夹以骨 松质，牙槽骨唇、颊侧骨板较薄，且有多数小孔通向其内的骨松质。因此牙槽骨是非均匀的， 是含气的多孔结构，而且其骨密度很接近牙齿结构，所以在评价牙槽骨吸收时，如何量化牙 槽骨是非常困难的。
[0003] 很多学者用Micro-CT (micro computed tomography,微计算机断层扫描技术）评 估实验性牙周炎的牙槽骨吸收情况。目前应用Micro-CT图像分析牙槽骨时，是按照密度进 行评估的。但是牙根与牙槽骨的密度又非常接近，常规Micro-CT图像操作对此自动识别能 力又比较差；如果手动方法进行人为识别，又会费时费力，并且存在较大人为误差。如：
[0004] 采用2D线性分析的方法评估牙槽骨吸收，与组织切片类似，需要选取很多张结 构清晰、标志点全面的2D切片进行测量，操作复杂繁琐，相对精确度差（Lu Sh eng-Hua， Evidence-Based Complementary and Alternative Medicine, 2013 ；Suda,Infection and immunity, 2013等）。还有学者用Micro-CT，扫描实验性牙周炎的颌骨标本，进行牙槽骨吸 收的3D体积分数评估。该方法的3D感兴趣区域（ROI)边界是在2D切片上先进行定点， 然后3D重建之后进行体积参数分析；但是以上述及的解剖标志点，不能在同一张2D切片 上全部显示出来，这样定点就会存在困难和误差，而且工作量大，实际操作比较困难，重复 性和准确性都不高（（iuimarSes，Journal of periodontal research，2011;Trombettae， Journal of dental research, 2011 ；Bak, Journal of periodontology, 2010 等）。而密 歇根大学Park等在2007年提出评估实验性牙周炎2D及3D牙槽骨吸收的手动方法，相对 准确，有一定的可重复性，但是程序繁琐，自动化程度低，实际可操作性较差（Park，Journal of per iodontology,2007)〇
[0005] 所以，研究探讨一种准确具有高度重复性且自动化程度高的方法，对评估牙槽骨 吸收情况具有指导性意义。
【发明内容】
[0006] 本发明提供一种Micro-CT图像分析量化评估磨牙牙槽骨的方法及装置，用以评 估牙周炎牙槽骨余留量，克服现有技术的缺点，提供一种可靠、准确、可重复的Micro-CT图 像分析方法。
[0007] 为达到上述目的，本发明提供了一种Micro-CT图像分析量化评估磨牙牙槽骨的 方法，包括以下步骤：
[0008] Sl :获取包含待分析磨牙的牙槽骨与牙齿的Micro-CT图像；
[0009] S2 :以待分析磨牙的根分叉、根尖孔作为标志点选取感兴趣区域，具体步骤为：
[0010] S21 :调整所述牙槽骨与牙齿的Micro-CT图像的阈值（Threshold)后，进行牙齿 图像与牙槽骨图像的分离，具体为：通过对所述牙槽骨与牙齿的Micro-CT图像边缘进行 消蚀（Erode)调整，利用区域增长（Region Growing)方法标记牙齿图像，通过布尔运算 (Boolean Operations)提取牙槽骨图像，将消蚀过程中调整过的牙槽骨图像膨胀（Dilate) 到消蚀前的相应像素水平，从而分离牙齿图像与牙槽骨图像；
[0011] S22 :将分离出来的牙齿图像进行3D运算，得到牙齿的3D图像；
[0012] S23 :用球形包围（Sphere)的方法确定感兴趣区域的位置，以待分析磨牙的根分 叉中心点到近中根根尖孔中心点距离的平均值，作为球形包围的球形半径r，调整球形包围 的中心至待分析磨牙的根分叉中心处，将2D的球形包围区域3D化，该球形包围区域即为感 兴趣区域；
[0013] S3 :将感兴趣区域进行3D运算，得到3D的球形包围区域的总体积TV ;将感兴趣区 域内的牙槽骨图像进行3D运算，得到感兴趣区域内牙槽骨的骨小梁体积BV、骨表面积BS ; 根据球形包围区域的总体积TV、感兴趣区域内牙槽骨的骨小梁体积BV及骨表面积BS，计算 出感兴趣区域内的牙槽骨骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分离度Tb. Sp和结构 模型指数SMI。
[0014] 其中，计算骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分离度Tb. Sp和结构模型 指数SMI的公式是：
[0015] Tb. Th = 2BV/BS ；
[0016] Tb. N = BS/2TV ；
[0017] Tb. Sp = 2 (TV-BV) /BS ；
[0018] SMI = [BV(dBS/dr)]/BS2。
[0019] 本发明还提供一种Micro-CT图像分析量化评估磨牙牙槽骨的装置，包括：
[0020] 图像获取模块，用于获取包含待分析磨牙的牙槽骨与牙齿的Micro-CT图像；
[0021] 感兴趣区域划定模块，用于以待分析磨牙的根分叉、根尖孔作为标志点选取感兴 趣区域；所述感兴趣区域划定模块包括分离模块、牙齿3D运算模块、感兴趣区域3D化模块， 其中：
[0022] 所述分离模块，用于调整所述牙槽骨与牙齿的Micro-CT图像的阈值 (Threshold)，进行牙齿图像与牙槽骨图像的分离，具体为：通过对所述牙槽骨与牙齿的 Micro-CT图像边缘进行消蚀（Erode)调整，利用区域增长（Region Growing)方法标记牙齿 图像，通过布尔运算（Boolean Operations)提取牙槽骨图像，将消蚀过程中调整过的牙槽 骨图像膨胀（Dilate)到消蚀前的相应像素水平，从而分离牙齿图像与牙槽骨图像；
[0023] 所述牙齿3D运算模块，用于将分离出来的牙齿图像进行3D运算，得到牙齿的3D 图像；
[0024] 所述感兴趣区域3D化模块，用球形包围（Sphere)的方法确定感兴趣区域的位置， 以待分析磨牙的根分叉中心点到近中根根尖孔中心点距离的平均值，作为球形包围的球形 半径r，调整球形包围的中心至待分析磨牙的根分叉中心处，将2D的球形包围区域3D化，该 球形包围区域即为感兴趣区域；
[0025] 感兴趣区域3D运算模块：将感兴趣区域进行3D运算，得到球形包围区域的总体 积TV;将感兴趣区域内的牙槽骨图像进行3D运算，得到感兴趣区域内牙槽骨的骨小梁体积 BV、骨表面积BS ;根据球形包围区域的总体积TV、感兴趣区域内牙槽骨的骨小梁体积BV及 骨表面积BS，计算出感兴趣区域内的牙槽骨骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分 离度Tb. Sp和结构模型指数SMI。
[0026] 与现有技术相比，本发明的有益效果体现在：
[0027] 本发明提供一种Micro-CT图像分析量化评估磨牙牙槽骨的方法及装置，应 用区域增长的标记功能，依据像素连续性的特点，能准确分离感兴趣区域（Region of Interest，R0I)内的牙齿图像和牙槽骨图像；并且选取客观的解剖学标志点，例如根分叉、 根尖孔等进行定点，选取R0I，可以进行余留牙槽骨的标准化数据分析，操作简便，准确性、 可重复性好，3D分析方法高效省时（预期分析一个样本用时15?20min)。
【专利附图】
【附图说明】
[0028] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。
[0029] 图1为本发明一实施例的Micro-CT图像分析量化评估磨牙牙槽骨的方法实现的 流程图。
【具体实施方式】
[0030] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0031] 如图1所示，为本发明一实施例的Micro-CT图像分析量化评估磨牙（包括第一磨 牙、第二磨牙、第三磨牙）牙槽骨的方法实现的流程图。本发明提供的一种Micro-CT图像 分析量化评估磨牙牙槽骨的方法，包括以下步骤：
[0032] 在本发明的一个实施例中，利用C57小鼠牙槽骨与牙齿组织的Micro-CT图像，修 剪（Crop)掉多余组织的图像，获取包含待分析磨牙的牙槽骨与牙齿的Micro-CT图像。
[0033] 调整所述牙槽骨与牙齿的Micro-CT图像的阈值（Threshold)后，进行牙齿图像与 牙槽骨图像的分离，具体为：对所述牙槽骨与牙齿的Micro-CT图像边缘进行消蚀（Erode) 调整，利用区域增长（Region Growing)方法标记牙齿图像，通过布尔运算（Boolean Operations)提取牙槽骨图像，将消蚀过程中调整过的牙槽骨图像膨胀（Dilate)到消蚀 (Erode)前的相应像素水平，从而分离牙齿图像与牙槽骨图像；
[0034] 将分离出来的牙齿图像进行3D运算，得到牙齿的3D图像；
[0035] 用球形包围（Sphere)的方法确定感兴趣区域的位置：以待分析磨牙的根分叉中 心点到近中根根尖孔中心点距离的平均值，作为球形包围的球形半径r，调整球形包围的中 心至待分析磨牙的根分叉中心处，将2D的球形包围区域3D化，该球形包围区域即为感兴趣 区域；
[0036] 将感兴趣区域进行3D运算，得到3D的球形包围区域的总体积TV ;将感兴趣区域 内的牙槽骨图像进行3D运算，得到感兴趣区域内牙槽骨的骨小梁体积BV、骨表面积BS ;根 据球形包围区域的总体积TV、感兴趣区域内牙槽骨的骨小梁体积BV及骨表面积BS，计算出 感兴趣区域内的牙槽骨骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分离度Tb. Sp和结构模 型指数SMI。
[0037] 在本发明的一个实施例中，计算骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分离 度Tb. Sp和结构模型指数SMI的公式是：
[0038] Tb. Th = 2BV/BS ；
[0039] Tb. N = BS/2TV ；
[0040] Tb. Sp = 2 (TV-BV) /BS ；
[0041] SMI = [BV(dBS/dr)]/BS2 ;其中dBS/dr表示骨表面积对球形包围的半径r求微分。
[0042] 本发明还提供一种Micro-CT图像分析量化评估磨牙牙槽骨的装置，包括：
[0043] 图像获取模块，用于获取包含待分析磨牙的牙槽骨与牙齿的Micro-CT图像；
[0044] 感兴趣区域划定模块，用于以待分析磨牙的根分叉、根尖孔作为标志点选取感兴 趣区域；所述感兴趣区域划定模块包括分离模块、牙齿3D运算模块、感兴趣区域3D化模块， 其中：
[0045] 所述分离模块，用于调整所述牙槽骨与牙齿的Micro-CT图像的阈值 (Threshold)，对牙齿图像与牙槽骨图像进行分离，具体为：通过对所述牙槽骨与牙齿的 Micro-CT图像边缘进行消蚀（Erode)调整，利用区域增长（Region Growing)标记牙齿图 像，通过布尔运算（Boolean Operations)提取牙槽骨图像，将消蚀过程中调整过的牙槽骨 图像膨胀（Dilate)到消蚀前的相应像素水平，从而分离牙齿图像与牙槽骨图像；
[0046] 所述牙齿3D运算模块，用于将分离出来的牙齿图像进行3D运算，得到牙齿的3D 图像；
[0047] 所述感兴趣区域3D化模块，用球形包围的方法确定感兴趣区域的位置，选择待分 析磨牙的根分叉中心点到近中根根尖孔中心点距离的平均值，作为球形包围的球形半径r， 调整球形包围的中心至待分析磨牙的根分叉中心处，将2D的球形包围区域3D化，该球形包 围区域即为感兴趣区域；
[0048] 感兴趣区域3D运算模块：将感兴趣区域进行3D运算，得到球形包围区域的总体 积TV;将感兴趣区域内的牙槽骨图像进行3D运算，得到感兴趣区域内牙槽骨的骨小梁体积 BV、骨表面积BS ;根据球形包围区域的总体积TV、感兴趣区域内牙槽骨的骨小梁体积BV及 骨表面积BS，计算出感兴趣区域内的牙槽骨骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分 离度Tb. Sp和结构模型指数SMI。
[0049] 将本发明一个实施例的3D分析牙槽骨的量化结果用2D分析的数据进行线性回归 验证，得到的的结果是R 2X). 996,两者结果趋于一致性，说明了本方法得到的结果准确性。
[0050] 本领域普通技术人员可以理解：附图只是一个实施例的示意图，附图中的模块或 流程并不一定是实施本发明所必须的。
[0051] 本领域普通技术人员可以理解：实施例中的装置中的模块可以按照实施例描述分 布于实施例的装置中，也可以进行相应变化位于不同于本实施例的一个或多个装置中。上 述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块。
[0052] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽 管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然 可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换； 而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范 围。
【权利要求】
1. 一种Micro-CT图像分析量化评估磨牙牙槽骨的方法，其特征在于，包括以下步骤： 51 :获取包含待分析磨牙的牙槽骨与牙齿的Micro-CT图像； 52 :以待分析磨牙的根分叉、根尖孔作为标志点选取感兴趣区域，具体步骤为： 521 :调整所述牙槽骨与牙齿的Micro-CT图像的阈值后，进行牙齿图像与牙槽骨图像 的分离，具体为：通过对所述牙槽骨与牙齿的Micro-CT图像的边缘进行消蚀调整，利用区 域增长方法标记牙齿图像，通过布尔运算提取牙槽骨图像，将消蚀过程中调整过的牙槽骨 图像膨胀到消蚀前的相应像素水平，从而分离牙齿图像与牙槽骨图像； 522 :将分离出来的牙齿图像进行3D运算，得到牙齿的3D图像； 523 :用球形包围的方法确定感兴趣区域的位置：以待分析磨牙的根分叉中心点到近 中根根尖孔中心点距离的平均值，作为球形包围的球形半径r，调整球形包围的中心至待分 析磨牙的根分叉中心处，将2D的球形包围区域3D化，该球形包围区域即为感兴趣区域； 53 :将感兴趣区域进行3D运算，得到3D的球形包围区域的总体积TV ;将感兴趣区域内 的牙槽骨图像进行3D运算，得到感兴趣区域内牙槽骨的骨小梁体积BV、骨表面积BS ;根据 球形包围区域的总体积TV、感兴趣区域内牙槽骨的骨小梁体积BV及骨表面积BS，计算出感 兴趣区域内牙槽骨的骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分离度Tb. Sp和结构模型 指数SMI。
2. 如权利要求1所述的Micro-CT图像分析量化评估磨牙牙槽骨的方法，其特征在于， 其中计算骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分离度Tb. Sp和结构模型指数SMI的 公式是： Tb.Th = 2BV/BS ； Tb. N = BS/2TV ； Tb.Sp = 2 (TV-BV)/BS ； SMI = [BV(dBS/dr)]/BS2。
3. -种Micro-CT图像分析量化评估磨牙牙槽骨的装置，其特征在于，包括： 图像获取模块，用于获取包含待分析磨牙的牙槽骨与牙齿的Micro-CT图像； 感兴趣区域划定模块，用于以待分析磨牙的根分叉、根尖孔作为标志点选取感兴趣区 域；所述感兴趣区域划定模块包括分离模块、牙齿3D运算模块、感兴趣区域3D化模块，其 中： 所述分离模块，用于调整牙槽骨与牙齿的Micro-CT图像的阈值后，进行牙齿图像与牙 槽骨图像的分离，具体为：通过对所述牙槽骨与牙齿的Micro-CT图像边缘进行消蚀调整， 利用区域增长方法标记牙齿图像，通过布尔运算提取牙槽骨图像，将消蚀过程中调整过的 牙槽骨图像膨胀到消蚀前的相应像素水平，从而分离牙齿图像与牙槽骨图像； 所述牙齿3D运算模块，用于将分离出来的牙齿图像进行3D运算，得到牙齿的3D图像； 所述感兴趣区域3D化模块，用球形包围的方法确定感兴趣区域，以待分析磨牙的根分 叉中心点到近中根根尖孔中心点距离的平均值，作为球形包围的球形半径r，调整球形包围 的中心至待分析磨牙的根分叉中心处，将2D的球形包围区域3D化，该球形包围区域即为感 兴趣区域； 感兴趣区域3D运算模块：将感兴趣区域进行3D运算，得到球形包围区域的总体积TV ; 将感兴趣区域内的牙槽骨图像进行3D运算，得到感兴趣区域内牙槽骨的骨小梁体积BV、骨 表面积BS ;根据球形包围区域的总体积TV、感兴趣区域内牙槽骨的骨小梁体积BV及骨表 面积BS，计算出感兴趣区域内牙槽骨的骨小梁数目Tb. N、骨小梁厚度Tb. Th、骨小梁分离度 Tb. Sp和结构模型指数SMI。
【文档编号】G06T17/00GK104240241SQ201410452800
【公开日】2014年12月24日 申请日期:2014年9月5日 优先权日:2014年9月5日
【发明者】林江, 张旭, 刘雪飞 申请人:哈尔滨医科大学