CN104751516B - 开髓位置的参照点的获取装置 - Google Patents
开髓位置的参照点的获取装置 Download PDFInfo
- Publication number
- CN104751516B CN104751516B CN201510167190.0A CN201510167190A CN104751516B CN 104751516 B CN104751516 B CN 104751516B CN 201510167190 A CN201510167190 A CN 201510167190A CN 104751516 B CN104751516 B CN 104751516B
- Authority
- CN
- China
- Prior art keywords
- mrow
- munderover
- msub
- sigma
- target tooth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种开髓位置的参照点的获取装置，属于口腔治疗技术领域。所述装置包括：三维模型建立模块，用于从包含目标牙齿的计算机断层扫描CT切片中截取目标牙齿处中间位置的图像，通过对所述图像切片进行分割计算，得到目标牙齿的三维模型；根管分割模块，用于从包含目标牙齿的CT切片中分割出目标牙齿的根管，并分离出预设长度的上段根管；交点计算模块，用于提取所述上段根管的中心线，计算所述中心线与所述目标牙齿三维模型的交点；开髓位置确定模块，用于以所述中心线与所述目标牙齿三维模型的交点为开髓位置的参照点。本发明不需要去除整个髓室顶且能更精确找到根管的位置，显著减小了开髓口的直径并最大限度保存髓室侧壁的牙本质。
Description
技术领域
本发明涉及口腔治疗技术领域，特别涉及一种开髓位置的参照点的获取装置。
背景技术
根管治疗是现代口腔治疗的常规手段，主要用于治疗牙髓病和根尖周病。在治疗过程中，根管治疗包括开髓、根管预备、根管消毒和根管填充等步骤。其中，开髓是根管治疗的关键步骤，现有的开髓方法包括：①去除髓室顶；②去除部分髓室侧壁的牙本质以消除器械进入根管的阻力。
在实现本发明的过程中，发明人发现现有技术至少存在以下问题：
由于医生不能基于根管上段中轴在患牙牙冠方向的延长线确定切削路径，在开髓过程中很容易导致开髓口的设计过大，使得很多健康牙本质被不必要的去除，削弱了牙体组织的抗力，患牙易于发生纵折。此外，在开髓过程中也很容易出现切削部位、方向和深度不准确，从而导致开髓口过大、髓室侧壁和底壁台阶和穿孔形成等问题。
发明内容
为了解决现有技术的问题，本发明实施例提供了一种开髓位置的参照点的获取装置，所述装置包括：
三维模型建立模块，用于从包含目标牙齿的计算机断层扫描CT切片中截取患牙处中间位置的图像，通过对图像进行分割计算，得到目标牙齿的三维模型；
根管分割模块，用于从包含目标牙齿的CT切片中分割出目标牙齿的根管，并分离出预设长度的上段根管；
交点计算模块，用于提取所分割根管的中心线，计算中心线与目标牙齿三维模型的交点；
开髓位置确定模块，用于以中心线与目标牙齿三维模型的交点为开髓位置的参照点。
可选地，所述三维模型建立模块还用于
S1、在所述切片中勾画出目标牙齿的初始轮廓，轮廓内部区域为前景，属于目标牙齿区域，轮廓外部区域为后景；
S2、计算轮廓中的代表前景的像素IF和代表后景的像素IB；
其中，IF和IB分别通过式(1)和式(2)计算；
其中p为像素点，F为前景区域，B为后景区域，I(p)为p点的灰度值，为前景中灰度值为I(p)的像素点的百分比，为后景中灰度值为I(p)的像素点的百分比；
S3、通过式(3)和式(4)设计数据项，利用式(5)设计平滑项，从而建立能量函数
Cs＝|Ip-IF| (3)，
Ct＝|Ip-IB| (4)，
Cp＝0.5 (5)；
其中，Cs代表前景数据项，Ct代表后景数据项，Cp代表平滑项，Ip为像素点p的灰度值；
S4、利用现有的最大流算法对能量函数进行求解从而得到图像的分割结果；
S5、对所述分割结果进行腐蚀操作，将上一步的分割结果减去腐蚀操作的结果得到目标区域的轮廓；
S6、将所述目标区域的轮廓分别作为位于当前CT切片上一片和下一片CT切片的输入重复S2至S6得到最终的牙齿分割结果；
S7、对所得的分割结果进行表面建模，得到目标牙齿的三维模型，所述模型由三角网格组成，只包含目标牙齿的表面轮廓几何位置信息。
可选地，所述交点计算模块还用于
S1、从所述目标牙齿的三维模型提取根管表面模型；
S2、设中心线的数学表达形式为目标是求的a，b，c和d的值；
S3、利用式(7)计算a，b，c和d的值从而得到中心线；
其中，m为根管表面模型所包含的点的数量，xi,yi,zi为根管表面模型第i个点的空间坐标。
可选地，所述交点计算模块还用于
S1、确定所述目标牙齿三维模型上表面所包含的多个三角面片，所述多个三角面片中只有一个三角形面片所在的平面与中心线的交点在该三角面片内部；
S2、计算所述目标牙齿三维模型上表面每一个三角面片所在平面与中心的交点，如果交点在三角面片内部，则该交点即为所求交点。
本发明实施例提供的技术方案带来的有益效果是：
本发明利用图像分割算法，首先从牙术前CT切片中分割出目标牙齿，建立目标牙齿的三维模型，然后，计算该目标牙齿指定长度的上段根管的中心线与目标牙齿的三维模型的交点，以该交点为参照点，利用导航***，不需要去除整个髓室顶且更精确找到根管的位置，显著减小了开髓口的直径并最大限度保存髓室侧壁的牙本质，与传统开髓技术相比，由于保留了更多的健康牙体组织，患牙折裂的机会减少；并且更准确的切削部位、方向和深度的控制可以有效避免开髓口过大、髓室侧壁和底壁台阶和穿孔形成等并发症。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的开髓位置的参照点的获取方法流程图；
图2是本发明实施例提供的开髓位置的参照点的获取装置结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的开髓位置的参照点的获取方法流程图。参见图1，该实施例包括：
101、从包含目标牙齿的计算机断层扫描CT切片中截取所述目标牙齿处中间位置的图像，通过对所述切片进行分割计算，得到所述目标牙齿的三维模型；
可选地，在本发明实施例中，可以利用图像分割算法(如Graph-Cut、LevelSet、Snake等)，从包含目标牙齿的计算机断层扫描CT切片中截取患牙处于中间位置的图像。
在本发明实施例中，从包含目标牙齿的CT切片中截取所述患牙处中间位置的图像，通过对切片进行分割计算，得到目标牙齿的三维模型包括：
S1、在所述切片中勾画出目标牙齿的初始轮廓，轮廓内部区域为前景，属于目标牙齿区域，轮廓外部区域为后景；
S2、计算轮廓中的代表前景的像素IF和代表后景的像素IB；
其中，IF和IB分别通过式(1)和式(2)计算；
其中p为像素点，F为前景区域，B为后景区域，I(p)为p点的灰度值，为前景中灰度值为I(p)的像素点的百分比，为后景中灰度值为I(p)的像素点的百分比；
S3、通过式(3)和式(4)设计数据项，利用式(5)设计平滑项，从而建立能量函数
Cs＝|Ip-IF| (3)，
Ct＝|Ip-IB| (4)，
Cp＝0.5 (5)；
其中，Cs代表前景数据项，Ct代表后景数据项，Cp代表平滑项，Ip为像素点p的灰度值；
S4、利用现有的最大流算法对能量函数进行求解从而得到图像的分割结果；
S5、对所述分割结果进行腐蚀操作，将上一步的分割结果减去腐蚀操作的结果得到目标区域的轮廓；
S6、将所述目标区域的轮廓分别作为位于当前CT切片上一片和下一片CT切片的输入重复S2至S6得到最终的牙齿分割结果；
S7、对所得的分割结果进行表面建模，得到目标牙齿的三维模型，所述模型由三角网格组成，只包含目标牙齿的表面轮廓几何位置信息。
102、从包含目标牙齿的CT切片中分割出所述目标牙齿的根管，并分离出预设长度的上段根管；
可选地，在本发明实施例中，可以利用图像分割算法(如Graph-Cut、LevelSet、Snake等)，从包含目标牙齿的CT切片中分割出所述目标牙齿的根管。
在本发明实施例中，分离出牙髓下预设长度的根管，优选为2-5mm。
103、提取所述上段根管的中心线，计算所述中心线与所述目标牙齿三维模型的交点；
可选地，在本发明实施例中可以利用中心线提取算法(如最小二乘法、Skeleton等)提取所分割根管的中心线。
在本发明实施例中，可选地，提取所分割根管的中心线包括：
S1、从所述目标牙齿的三维模型提取根管表面模型；
S2、设中心线的数学表达形式为目标是求的a，b，c和d的值；
S3、利用式(7)计算a，b，c和d的值从而得到中心线；
其中，m为根管表面模型所包含的点的数量，xi,yi,zi为根管表面模型第i个点的空间坐标。
可选地，在本发明实施例中，计算所述中心与所述目标牙齿三维模型的交点包括：
S1、确定所述目标牙齿三维模型上表面所包含的多个三角面片，所述多个三角面片中只有一个三角形面片所在的平面与中心线的交点在该三角面片内部；
S2、计算所述目标牙齿三维模型上表面每一个三角面片所在平面与中心的交点，如果交点在三角面片内部，则该交点即为所求交点。
104、以所述中心线与所述目标牙齿三维模型的交点为开髓位置的参照点，根据具体图像分析并调整切入点位置，利用导航***最终确定所述目标牙齿的开髓位置。
在本发明实施例中，导航方式可以是光学球定位或导板定位，可以将开髓位置精确到0.5mm。
本发明实施例提供的方法，利用图像分割算法，首先从牙术前CT切片中分割出目标牙齿，建立目标牙齿的三维模型，然后，计算该目标牙齿指定长度的上段根管的中心线与目标牙齿的三维模型的交点，以该交点为参照点，利用导航***，不需要去除整个髓室顶且更精确找到根管的位置，显著减小了开髓口的直径并最大限度保存髓室侧壁的牙本质，与传统开髓技术相比，由于保留了更多的健康牙体组织，患牙折裂的机会减少；并且更准确的切削部位、方向和深度的控制可以有效避免开髓口过大、髓室侧壁和底壁台阶和穿孔形成等并发症。
图2是本发明实施例提供的一种开髓位置的参照点的获取装置结构示意图。参见图2，该装置包括：
三维模型建立模块21，用于从包含目标牙齿的计算机断层扫描CT切片中截取所述患牙处于中间位置的图像，通过对所述切片进行分割计算，得到所述目标牙齿的三维模型；
根管分割模块22，用于从包含目标牙齿的CT切片中分割出所述目标牙齿的根管，并分离出预设长度的上段根管；
交点计算模块23，用于提取所述上段根管的中心线，计算所述中心线与所述目标牙齿三维模型的交点；
开髓位置确定模块24，用于以所述中心线与所述目标牙齿三维模型的交点为参照点，根据具体图像分析并调整切入点位置，利用导航***最终确定所述目标牙齿的开髓位置。
可选地，所述三维模型建立模块21还用于
S1、在所述切片中勾画出目标牙齿的初始轮廓，轮廓内部区域为前景，属于目标牙齿区域，轮廓外部区域为后景；
S2、计算轮廓中的代表前景的像素IF和代表后景的像素IB；
其中，IF和IB分别通过式(1)和式(2)计算；
其中p为像素点，F为前景区域，B为后景区域，I(p)为p点的灰度值，为前景中灰度值为I(p)的像素点的百分比，为后景中灰度值为I(p)的像素点的百分比；
S3、通过式(3)和式(4)设计数据项，利用式(5)设计平滑项，从而建立能量函数
Cs＝|Ip-IF| (3)，
Ct＝|Ip-IB| (4)，
Cp＝0.5 (5)；
其中，Cs代表前景数据项，Ct代表后景数据项，Cp代表平滑项，Ip为像素点p的灰度值；
S4、利用现有的最大流算法对能量函数进行求解从而得到图像的分割结果；
S5、对所述分割结果进行腐蚀操作，将上一步的分割结果减去腐蚀操作的结果得到目标区域的轮廓；
S6、将所述目标区域的轮廓分别作为位于当前CT切片上一片和下一片CT切片的输入重复S2至S6得到最终的牙齿分割结果；
S7、对所得的分割结果进行表面建模，得到目标牙齿的三维模型，所述模型由三角网格组成，只包含目标牙齿的表面轮廓几何位置信息。
可选地，所述交点计算模块23还用于
S1、从所述目标牙齿的三维模型提取根管表面模型；
S2、设中心线的数学表达形式为目标是求的a，b，c和d的值；
S3、利用式(7)计算a，b，c和d的值从而得到中心线；
其中，m为根管表面模型所包含的点的数量，xi,yi,zi为根管表面模型第i个点的空间坐标。
可选地，所述交点计算模块23还用于
S1、确定所述目标牙齿三维模型上表面所包含的多个三角面片，所述多个三角面片中只有一个三角形面片所在的平面与中心线的交点在该三角面片内部；
S2、计算所述目标牙齿三维模型上表面每一个三角面片所在平面与中心的交点，如果交点在三角面片内部，则该交点即为所求交点。
本发明实施例提供的装置，利用图像分割算法，首先从牙术前CT切片中分割出目标牙齿，建立目标牙齿的三维模型，然后，计算该目标牙齿指定长度的上段根管的中心线与目标牙齿的三维模型的交点，以该交点为参照点，利用导航***，不需要去除整个髓室顶且更精确找到根管的位置，显著减小了开髓口的直径并最大限度保存髓室侧壁的牙本质，与传统开髓技术相比，由于保留了更多的健康牙体组织，患牙折裂的机会减少；并且更准确的切削部位、方向和深度的控制可以有效避免开髓口过大、髓室侧壁和底壁台阶和穿孔形成等并发症。
需要说明的是：上述实施例提供的装置，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种开髓位置的参照点的获取装置，其特征在于，所述装置包括：
三维模型建立模块，用于从包含目标牙齿的计算机断层扫描CT切片中截取所述目标牙齿处中间位置的图像，通过对所述切片进行分割计算，得到所述目标牙齿的三维模型；
根管分割模块，用于从包含目标牙齿的CT切片中分割出所述目标牙齿的根管，并分离出预设长度的上段根管；
交点计算模块，用于提取所上段根管的中心线，计算所述中心线与所述目标牙齿三维模型的交点；
开髓位置确定模块，用于以所述中心线与所述目标牙齿三维模型的交点为开髓位置的参照点。
2.如权利要求1所述的装置，其特征在于，所述三维模型建立模块还用于
S1、在所述切片中勾画出目标牙齿的初始轮廓，轮廓内部区域为前景，属于目标牙齿区域，轮廓外部区域为后景；
S2、计算轮廓中的代表前景的像素IF和代表后景的像素IB；
其中，IF和IB分别通过式(1)和式(2)计算；
<mrow>
<msub>
<mi>I</mi>
<mi>F</mi>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>&Element;</mo>
<mi>F</mi>
</mrow>
</munder>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>F</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mi>B</mi>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>&Element;</mo>
<mi>B</mi>
</mrow>
</munder>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>p</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>B</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中p为像素点，F为前景区域，B为后景区域，I(p)为p点的灰度值，为前景中灰度值为I(p)的像素点的百分比，为后景中灰度值为I(p)的像素点的百分比；
S3、通过式(3)和式(4)设计数据项，利用式(5)设计平滑项，从而建立能量函数
Cs＝|Ip-IF| (3)，
Ct＝|Ip-IB| (4)，
Cp＝0.5 (5)；
其中，Cs代表前景数据项，Ct代表后景数据项，Cp代表平滑项，Ip为像素点p的灰度值；
S4、利用现有的最大流算法对能量函数进行求解从而得到图像的分割结果；
S5、对所述分割结果进行腐蚀操作，将上一步的分割结果减去腐蚀操作的结果得到目标区域的轮廓；
S6、将所述目标区域的轮廓分别作为位于当前CT切片上一片和下一片CT切片的输入重复S2至S6得到最终的牙齿分割结果；
S7、对所得的分割结果进行表面建模，得到目标牙齿的三维模型，所述模型由三角网格组成，只包含目标牙齿的表面轮廓几何位置信息。
3.如权利要求1所述的装置，其特征在于，所述交点计算模块还用于
S1、从所述目标牙齿的三维模型提取根管表面模型；
S2、设中心线的数学表达形式为目标是求的a，b，c和d的值；
S3、利用式(7)计算a，b，c和d的值从而得到中心线；
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>b</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>m</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mi>d</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>d</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<mi>m</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，m为根管表面模型所包含的点的数量，xi,yi,zi为根管表面模型第i个点的空间坐标。
4.如权利要求1所述的装置，其特征在于，所述交点计算模块还用于
S1、确定所述目标牙齿三维模型上表面所包含的多个三角面片，所述多个三角面片中只有一个三角形面片所在的平面与中心线的交点在该三角面片内部；
S2、计算所述目标牙齿三维模型上表面每一个三角面片所在平面与中心的交点，如果交点在三角面片内部，则该交点即为所求交点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510167190.0A CN104751516B (zh) | 2015-04-09 | 2015-04-09 | 开髓位置的参照点的获取装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510167190.0A CN104751516B (zh) | 2015-04-09 | 2015-04-09 | 开髓位置的参照点的获取装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104751516A CN104751516A (zh) | 2015-07-01 |
| CN104751516B true CN104751516B (zh) | 2018-01-09 |
Family
ID=53591134
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510167190.0A Active CN104751516B (zh) | 2015-04-09 | 2015-04-09 | 开髓位置的参照点的获取装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104751516B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104958118B (zh) * | 2015-07-20 | 2016-08-17 | 武汉大学 | 一种用于牙齿开髓牙体预备的智能机器人 |
| CN105596091B (zh) * | 2016-03-14 | 2019-09-27 | 四川大学 | 用于根管治疗的导航装置及其制备方法 |
| CN106308952A (zh) * | 2016-10-27 | 2017-01-11 | 苏州迪凯尔医疗科技有限公司 | 基于导航技术的开髓器械定位方法及装置 |
| CN106510865B (zh) * | 2016-11-11 | 2019-08-09 | 北京大学口腔医学院 | 牙齿开髓洞形预备的引导装置和制作该装置的方法 |
| CN108932716B (zh) * | 2017-05-26 | 2020-09-22 | 无锡时代天使医疗器械科技有限公司 | 用于牙齿图像的图像分割方法 |
| US11416999B2 (en) * | 2018-05-29 | 2022-08-16 | Medicim Nv | Methods, systems, and computer programs for segmenting a tooth's pulp region from an image |
| CN109846561A (zh) * | 2018-12-26 | 2019-06-07 | 苏州迪凯尔医疗科技有限公司 | 数字化根管预备的方法、装置及*** |
| CN117814934B (zh) * | 2024-01-31 | 2024-06-14 | 武汉欧若博机器人有限公司 | 一种根管预备过程可视化方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103826566A (zh) * | 2011-09-05 | 2014-05-28 | 美特瑞牙科有限公司 | 用于3d根管治疗计划的方法和*** |
| CN104367387A (zh) * | 2014-10-13 | 2015-02-25 | 浙江工业大学 | 模拟自然牙力学性能的个性化金属牙桩及其制作方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10130450B2 (en) * | 2013-05-14 | 2018-11-20 | Ipg Photonics Corporation | Method and apparatus for laser induced thermo-acoustical streaming of liquid |
| US9579166B2 (en) * | 2013-08-07 | 2017-02-28 | Jeffrey I Lasner | Root canal probe tool and method of removing a broken instrument fragment from a root canal |
-
2015
- 2015-04-09 CN CN201510167190.0A patent/CN104751516B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103826566A (zh) * | 2011-09-05 | 2014-05-28 | 美特瑞牙科有限公司 | 用于3d根管治疗计划的方法和*** |
| CN104367387A (zh) * | 2014-10-13 | 2015-02-25 | 浙江工业大学 | 模拟自然牙力学性能的个性化金属牙桩及其制作方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104751516A (zh) | 2015-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104751516B (zh) | 开髓位置的参照点的获取装置 | |
| CN110189352B (zh) | 一种基于口腔cbct图像的牙根提取方法 | |
| CN104809723B (zh) | 基于超体素和图割算法的三维肝脏ct图像自动分割方法 | |
| CN106327487B (zh) | 基于x射线冠脉造影图像的冠状动脉血流储备分数计算方法 | |
| CN104143184B (zh) | 一种肺部分割的方法 | |
| CN105574933B (zh) | 物体全方位轮廓精确绘制方法 | |
| CN106355582B (zh) | 一种基于形状先验的肝脏医学图像分割的方法 | |
| CN105118056A (zh) | 基于三维形态学的冠脉血管自动提取方法 | |
| CN102048550A (zh) | 一种自动生成肝脏3d图像并准确定位肝脏血管支配区域的方法 | |
| CN104077806B (zh) | 基于城市建筑三维模型的自动分体提取方法 | |
| CN108038862A (zh) | 一种交互式医学图像智能分割建模方法 | |
| CN103606190A (zh) | 单张人脸正面照片自动转换为三维人脸模型的方法 | |
| CN102254317A (zh) | 牙种植导航中牙弓曲面的自动提取方法 | |
| CN108230342A (zh) | 一种基于心脏解剖结构的左右心室水平集分割方法 | |
| CN107895364A (zh) | 一种用于虚拟手术术前规划的三维重建*** | |
| CN106806021A (zh) | 一种基于人体器官3d模型的vr手术模拟***及方法 | |
| CN106504277A (zh) | 一种改进的icp点云自动配准方法 | |
| CN114757960B (zh) | 一种基于cbct图像的牙齿分割与重建方法及存储介质 | |
| CN107123163A (zh) | 一种基于多视角立体视觉的植株三维重建*** | |
| CN106408576A (zh) | 基于三维超声图像的感兴趣区域的自动分割方法及*** | |
| CN106780497A (zh) | 一种基于统计信息的器官血管树自动提取方法 | |
| CN108537750A (zh) | 一种快速的脊柱ct图像三维重建方法 | |
| CN107610095A (zh) | 基于图像融合的心脏ct冠脉全自动分割方法 | |
| CN106327536A (zh) | 一种基于切面点云的领围测量方法 | |
| CN105574882B (zh) | 基于胸部横断面ct图像的肺部分割提取方法及*** |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220516 Address after: 430079 Luojia mountain, Wuchang, Wuhan City, Hubei Province Patentee after: WUHAN University Address before: 430079 No. 237-051-2 Luoyu Road, Hongshan District, Wuhan City, Hubei Province Patentee before: Fan Bing |
|
| TR01 | Transfer of patent right |