CN100553587C - 在三维网格牙颌模型上快速精确探测牙弓线的方法 - Google Patents
在三维网格牙颌模型上快速精确探测牙弓线的方法 Download PDFInfo
- Publication number
- CN100553587C CN100553587C CNB2007100254281A CN200710025428A CN100553587C CN 100553587 C CN100553587 C CN 100553587C CN B2007100254281 A CNB2007100254281 A CN B2007100254281A CN 200710025428 A CN200710025428 A CN 200710025428A CN 100553587 C CN100553587 C CN 100553587C
- Authority
- CN
- China
- Prior art keywords
- point
- arch wire
- biteplane
- initial
- final
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000003993 interaction Effects 0.000 claims abstract description 6
- 238000013507 mapping Methods 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 5
- 238000005211 surface analysis Methods 0.000 claims description 4
- 241000282341 Mustela putorius furo Species 0.000 claims description 2
- 238000011960 computer-aided design Methods 0.000 abstract description 2
- 230000002452 interceptive effect Effects 0.000 abstract 1
- 230000002159 abnormal effect Effects 0.000 description 5
- 210000002455 dental arch Anatomy 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 210000001847 jaw Anatomy 0.000 description 2
- 238000007620 mathematical function Methods 0.000 description 2
- 239000011505 plaster Substances 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 208000003227 Jaw Abnormalities Diseases 0.000 description 1
- DKMVJQCQTCLYIF-UHFFFAOYSA-M Methylbenactyzium bromide Chemical compound [Br-].C=1C=CC=CC=1C(O)(C(=O)OCC[N+](C)(CC)CC)C1=CC=CC=C1 DKMVJQCQTCLYIF-UHFFFAOYSA-M 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001373 regressive effect Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Landscapes
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
一种在三维网格牙颌模型上快速精确探测牙弓线的方法，属计算机辅助设计领域。本方法包括简单人机交互确定咬合平面、筛选初始参照点、拟合初始牙弓线、提取最终参照点、拟合最终牙弓线五个步骤。本发明的方法只包含有简单的人机交互操作，自动化程度高，同时避免了使用常用的复杂计算方法，以线性计算为主，能够快速、稳定、精确地在三维牙颌模型上探测牙弓线，在计算机辅助口腔正畸领域有着重要的应用价值。
Description
所属技术领域
本发明涉及计算机辅助设计领域，特别是计算机辅助口腔正畸治疗中，三维牙颌模型上的牙弓线探测方法。
背景技术
牙弓形状在正畸诊断和治疗中具有重要的意义，正畸医生需要绘制出上下颌的牙弓线以判断错颌畸形的程度，测量牙弓进行拥挤度分析是正畸诊断重要的定量评价依据。修正异常的牙弓形态，建立与患者个体匹配的个性化牙弓形状，是口腔正畸的重要目标。正畸医生最早是根据Bonwill-Hawley氏原理来绘制个别弓形图，但是绘制的过程极为复杂、繁琐。后来各种数学模型被用来模拟牙弓形态，比如椭圆线函数、抛物线函数等复杂的数学函数，绘制由数学函数表示的牙弓线更为繁琐、难度更大。
随着计算机技术在口腔正畸治疗中的广泛应用，人们开始使用计算机来绘制牙弓线。Begole首先采集牙颌1∶1的二维图像，在二维图像上手动拾取参照点，使用三次样条曲线拟合参照点得到牙弓线。三维测量技术应用于口腔正畸治疗以后，可以很容易地获得数字化的三维牙颌模型，于是出现了在三维数字模型上绘制牙弓线的方法，Joon Huang等在三维数字牙颌模型上手动选取参照点来拟合牙弓线，Schicho等在三维数字模型上人机交互定义牙弓线，这些手动的或者交互式的算法需要正畸医生的频繁操作、费时费力。Kondo在论文Robust archdetection and tooth segmentation in 3D images of dental plaster models(In：Proceedings International Workshop on Medical Imaging and AugmentedReality[C]，Hong Kong，China，2001：241-246)中介绍了一种自动探测牙弓线的方法，首先确定咬合平面，依据咬合平面调整牙颌模型的坐标，然后把三维模型转换为二维的灰度图像，灰度图像的亮度值代表坐标z值(即高度值)，在灰度图像上使用梯度方向分析(gradient orientation analysis)方法寻找梯度方向间断的点作为初始参照点，使用初始参照点拟合一条初始牙弓线，依据初始牙弓线在灰度图像上提取最终特征点，拟合提取到的最终特征点得到最终牙弓线，这种自动探测牙弓线的方法有良好的适应性，但是由于要把三维模型转换为二维的灰度图像，在转换的过程中部分网格顶点的信息会丢失，如果保证网格顶点信息不丢失则会引入大量的冗余数据，在探测初始特征点时使用的梯度方向分析方法涉及复杂的计算，计算数据量很大的精细模型会耗费大量的时间。
发明内容
针对现有方法的不足，提出了一种直接在三维牙颌模型上快速探测牙弓线的方法，不需要从三维到二维的转换，避免精度损失或引入冗余数据；避免使用复杂的计算，可以实现稳定、快速、精确的探测。
本方法的技术方案如下：
步骤1：简单人机交互确定咬合平面：在三维模型上交互选取四个参照点：左右两颗第一恒磨牙的近中颊侧尖点和左右两颗第一尖磨牙的颊侧尖点，使用趋势面分析方法确定咬合平面；
步骤2：筛选网格顶点得到初始参照点：首先将网格顶点中其法矢与咬合平面法矢夹角在0°～15°范围内的顶点保留，然后，保留下来的顶点中到咬合平面的距离在0～5mm范围内的顶点作为初始参照点集；
步骤3：确定初始牙弓线：把步骤2中在三维模型上选中的初始参照点集映射到二维平面，在二维平面上用四次多项式拟合映射点，得到初始牙弓线；
步骤4：依据初始牙弓线提取最终参照点：将初始牙弓线自动等分为40～60段，在每一段中选取一个初始参照点作为最终参照点，其具体方式为：当该段内由步骤2选取的所有初始参照点均在咬合平面牙根侧时，选取离咬合平面最近的一个初始参照点为该段的最终参照点，其它情况，则以咬合平面牙尖侧离咬合平面最远的一个初始参照点为最终参照点；
(5)、确定最终牙弓线：采用与第(3)步同样的方法，得到最终牙弓线。
本发明的方法直接在三维模型上探测牙弓线，不需要从三维到二维的转换，能够最大限度的利用已有三维牙颌模型的数据或者避免引入冗余数据，使用线性计算筛选初始参照点集，避免了使用复杂的计算，能够快速、自动、精确的提取患者个性化牙颌的牙弓线，算法稳定、可靠，在计算机正畸矫治有重要的临床和工程应用价值。
附图说明
图1牙弓线探测流程示意图。
图2确定咬合平面参照点示意图。
图2中符号名称：1、4-第一恒磨牙近中颊侧尖点；2、3-第一尖磨牙的颊侧尖点
图3初始牙弓线示意图。
图4等分初始牙弓线提取最终参照点示意图。
图5最终结果示意图。
具体实施方案
本发明的核心思想是利用咬合平面法矢与网格顶点法矢的夹角以及网格顶点到咬合平面的距离来筛选初始参照点集、利用拟合的初始牙弓线对初始参照点集进行二次筛选得到最终参照点，并最终快速、自动、精确拟合最终的牙弓线。
在利用本发明的方法探测牙弓线之前，首先利用国产高精度三维激光测量系统扫描得到石膏模型表面的点云数据，使用逆向工程软件处理扫描得到的点云数据，经过网格化后得到三维网格模型，由于在测量中不可避免的遇到噪声、扫描盲区等缺陷，需要对得到的模型进行光顺、修补等处理，最终得到比较理想的三维数字模型。该步过程已经是非常成熟的技术。
下面详细描述本发明的牙弓线探测算法。作为一种具体的实现方案，本算法主要包括五个步骤：简单人机交互确定咬合平面、筛选网格顶点得到初始参照点集、拟合初始牙弓线、依据初始牙弓线提取最终参照点、拟合最终牙弓线。下面结合附图对其逐一介绍。
(1)、简单人机交互确定咬合平面：如图2所示，在三维模型上交互选取四个参照点：左第一恒磨牙的近中颊侧尖点为第1点，左第一尖磨牙的颊侧尖点为第2点，右第一尖磨牙的颊侧尖点为第3点，右第一恒磨牙的近中颊侧尖点为第4点，使用趋势面分析方法确定咬合平面。
具体为：按照上述1、2、3、4点的顺序拾取参照点，由于直接在三维模型上选中所需要的参照点比较困难，所以程序不要求直接选中目标点，算法会自动在拾取点的附近搜索目标点，具体的搜索过程为：用拾取的四个点确定一个参照平面，远离牙根的方向定义为参照平面的正法向，计算拾取点周围一定范围内的网格顶点到参照面的距离，一般取距离拾取点2mm以内的网格顶点，把在参照面负法向侧的点到平面的距离定义为负值，选择距离值最大的点作为目标点。使用搜索到的目标点利用趋势面分析的方法来确定咬合平面，即把观测值Zj(xj，yj)的变化分解成两个部分：Zj(xj，yj)＝f(xj，yj)+σj，式中：f(xj，yj)为趋势值，σj为剩余值，采用回归的方法求得趋势值和剩余值，即根据已知数据Z的回归方程f(x，y)＝ax+by+c，使得 达到极小，求出系数得到咬合平面方程。其中xj，yj代表目标点的x，y坐标值，Zj代表目标点的z坐标值。
(2)、筛选网格顶点得到初始参照点集：首先将网格顶点中其法矢与咬合平面法矢夹角在0°～15°范围内的顶点保留，然后，保留下来的顶点中到咬合平面的距离在0～5mm范围内的顶点作为初始参照点集。
具体为：网格顶点表示为xi、yi和zi为网格顶点的坐标值，
为网格顶点的法矢，网格所有的顶点集合SV定义为： 集合SNOR表示顶点法矢
与咬合平面法矢
的夹角小于给定角度值
的顶点的集合，其定义为：
其中
表示顶点法矢
与咬合平面法矢
的夹角。集合SNOR *是集合SNOR的一个真子集，表示集合SNOR中到咬合面的距离小于给定值τ的顶点，SNOR *的定义为 dis(Vi，FOCC)表示顶点Vi与咬合平面FOCC的距离。集合SNOR *中的点即为初始参照点集。
(3)、确定初始牙弓线：把三维模型上选中的初始参照点集映射到XOY二维平面上，用四次多项式拟合映射点，得到初始牙弓线。
具体为：把初始参照点集沿咬合面法矢方向映射到咬合面上，然后再把咬合面上的点映射到XOY坐标平面上，以XOY坐标平面的坐标原点为中心，在XOY坐标平面内旋转投影点，使左右两颗第一恒磨牙的近中颊侧尖点在XOY坐标面上映射点的连线与XOY坐标平面的X坐标轴平行，带权误差ε的定义为：ε＝∑[w(x)(ax4+bx3+cx2+dx+e-y)]2，式中的w(x)是对误差使用的一个权因子，目的是为了使牙弓线比较精确的通过前牙的参照点，权因子的定义为： xm为XOY面上所有映射点的X坐标的中值，x为投影点的x坐标值。图3显示了得到的初始牙弓线。
(4)、依据初始牙弓线提取最终参照点：将初始牙弓线自动等分为40～60段，在每一段中选取一个初始参照点作为最终参照点，其具体方式为：当该段内所有初始参照点均在咬合平面牙根侧时，选取离咬合平面最近的一个初始参照点为该段的最终参照点，其它情况，则以咬合平面牙尖侧离咬合平面最远的一个初始参照点为最终参照点。
具体为：如图4所示，首先等分初始牙弓线，在等分点处插入法平面，相邻两个法平面之间的初始参照点集是初始参照点集合SNOR *的一个真子集，这样得到集合SNOR *的一个划分SNOR *1，SNOR *2，……，SNOR *n，对于任意i≠j≤n，
在集合SNOR *i中寻找一个在咬合面的负法向侧离咬合面最近(该集合中的点全部在咬合面的负法向侧)或在咬合面正法向侧离咬合面最远(该集合中的点部分或全部在咬合面的正法向侧)的点，得到一个顶点集合Sfinal，即最终参照点集合。
(5)、确定最终牙弓线：把三维模型上的最终参照点映射到二维平面上，用四次多项式拟合映射点，得到最终牙弓线。
具体为：把集合Sfinal按照前面所述的拟合初始牙弓线的方法，首先把集合Sfinal中的点映射到XOY坐标平面上，再用四次多项式拟合映射点，使用同样的带权误差确定多项式系数，得到最终的牙弓线，如图5所示。
Claims (2)
1、一种在三维网格牙颌模型上快速精确探测牙弓线的方法，其特征在于：
步骤1：简单人机交互确定咬合平面：在三维模型上交互选取四个参照点：左右两颗第一恒磨牙的近中颊侧尖点和左右两颗第一尖磨牙的颊侧尖点，使用趋势面分析方法确定咬合平面；
步骤2：筛选网格顶点得到初始参照点集：首先将网格顶点中其法矢与咬合平面法矢夹角在0°～15°范围内的顶点保留，然后，保留下来的顶点中到咬合平面的距离在0～5mm范围内的顶点作为初始参照点集；
步骤3：确定初始牙弓线：把步骤2中在三维模型上选中的初始参照点集映射到二维平面，在二维平面上用四次多项式拟合映射点，得到初始牙弓线；
步骤4：依据初始牙弓线提取最终参照点：将初始牙弓线自动等分为40～60段，在每一段中选取一个初始参照点作为最终参照点，其具体方式为：当该段内由步骤2选取的所有初始参照点均在咬合平面牙根侧时，选取离咬合平面最近的一个初始参照点为该段的最终参照点，其它情况，则以咬合平面牙尖侧离咬合平面最远的一个初始参照点为最终参照点；
步骤5：确定最终牙弓线：采用与步骤3同样的方法，得到最终牙弓线。
2、根据取利要求1所述的在三维网格牙颌模型上快速精确探测牙弓线的方法，其特征在于：在第一步确定咬合平面中，计算机会自动在交互选取的四个参照点附近搜索目标点，具体搜索方法为，用拾取的四个点确定一个参照平面，远离牙根的方向定义为参照平面的正法向，计算拾取点周围2mm以内的网格顶点到参照面的距离，把在参照面负法向侧的点到平面的距离定义为负值，选择距离值最大的点作为目标点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100254281A CN100553587C (zh) | 2007-07-31 | 2007-07-31 | 在三维网格牙颌模型上快速精确探测牙弓线的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2007100254281A CN100553587C (zh) | 2007-07-31 | 2007-07-31 | 在三维网格牙颌模型上快速精确探测牙弓线的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101103940A CN101103940A (zh) | 2008-01-16 |
| CN100553587C true CN100553587C (zh) | 2009-10-28 |
Family
ID=38998157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2007100254281A Expired - Fee Related CN100553587C (zh) | 2007-07-31 | 2007-07-31 | 在三维网格牙颌模型上快速精确探测牙弓线的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100553587C (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101543431B (zh) * | 2008-03-28 | 2012-02-29 | 上海交通大学医学院附属第九人民医院 | 一种正牙个体化牙弓弓形的绘制方法 |
| CN102609965B (zh) * | 2012-02-13 | 2014-12-03 | 北京朗视仪器有限公司 | 一种牙科cbct生成全景图的数据处理方法 |
| CN103494653B (zh) * | 2013-10-17 | 2015-10-14 | 福州大学 | 一种牙颌生长发育规律的数字化测量方法 |
| CN107658129B (zh) * | 2017-08-09 | 2020-11-27 | 佛山市众拓科技有限公司 | 一种基于可逻辑编程控制的五瓣花型电磁线圈的绕线方法 |
| CN107564090A (zh) * | 2017-09-07 | 2018-01-09 | 杭州美齐科技有限公司 | 一种三维牙颌模型渲染方法 |
| CN108814742B (zh) * | 2018-04-04 | 2021-12-07 | 正雅齿科科技(上海)有限公司 | 牙颌托盘的设计方法、装置、用户终端及存储介质 |
| CN110269703B (zh) * | 2019-06-18 | 2021-02-09 | 成都玻尔兹曼智贝科技有限公司 | 下颌牙齿模型的咬合曲线的生成方法 |
| CN112022408A (zh) * | 2020-07-09 | 2020-12-04 | 常州飞诺医疗技术有限公司 | 一种基于双目视觉获得人体牙弓曲线的方法 |
| CN115588006B (zh) * | 2022-11-11 | 2023-11-21 | 四川大学 | 一种标准化牙弓形态的提取方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1260698A (zh) * | 1997-06-20 | 2000-07-19 | 矫正技术公司 | 逐渐移位牙齿的方法和装置 |
| CN1686058A (zh) * | 2005-04-28 | 2005-10-26 | 上海隐齿丽医学技术有限公司 | 计算机辅助牙齿隐形正畸系统 |
-
2007
- 2007-07-31 CN CNB2007100254281A patent/CN100553587C/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1260698A (zh) * | 1997-06-20 | 2000-07-19 | 矫正技术公司 | 逐渐移位牙齿的方法和装置 |
| CN1686058A (zh) * | 2005-04-28 | 2005-10-26 | 上海隐齿丽医学技术有限公司 | 计算机辅助牙齿隐形正畸系统 |
Non-Patent Citations (2)
| Title |
|---|
| 口腔修复体咬合面虚拟调整. 张翔等.中国机械工程，第17卷第20期. 2006 |
| 口腔修复体咬合面虚拟调整. 张翔等.中国机械工程，第17卷第20期. 2006 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101103940A (zh) | 2008-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100553587C (zh) | 在三维网格牙颌模型上快速精确探测牙弓线的方法 | |
| CN108470365B (zh) | 一种基于上下牙颌的牙弓线绘制方法 | |
| US7474932B2 (en) | Dental computer-aided design (CAD) methods and systems | |
| CN103494653B (zh) | 一种牙颌生长发育规律的数字化测量方法 | |
| ES2627810T3 (es) | Método y sistema de análisis para el análisis geométrico de datos de exploración de estructuras orales | |
| CN103886145B (zh) | 牙齿预备体数字化模型设计方法 | |
| CN108986111B (zh) | 一种用于计算机辅助隐形正畸的三维牙颌模型分割方法 | |
| US8029277B2 (en) | Method and system for measuring tooth displacements on a virtual three-dimensional model | |
| CN108665533A (zh) | 一种通过牙齿ct图像和三维扫描数据重建牙列的方法 | |
| Vág et al. | A novel method for complex three-dimensional evaluation of intraoral scanner accuracy | |
| CN103871097B (zh) | 基于牙齿预备体的数据柔性融合方法 | |
| CN105030364B (zh) | 牙齿横纵向倾斜角度的测量方法 | |
| CN104715099B (zh) | 用于牙齿冠桥功能性咬合面的虚拟牙合架设计方法 | |
| CN101828974A (zh) | 种植义齿个性化定位导板的制造方法 | |
| Zhao et al. | Interactive tooth segmentation of dental models | |
| CN101706971A (zh) | 一种牙齿模型中牙冠的自动分割方法 | |
| CN113052902B (zh) | 牙齿治疗监测方法 | |
| CN112790879B (zh) | 一种牙齿模型的牙轴坐标系构建方法及系统 | |
| CN104867131A (zh) | 基于数字化模型的牙冠数据提取方法 | |
| Ding et al. | Morphology and mechanical performance of dental crown designed by 3D-DCGAN | |
| CN105389444A (zh) | 一种个性化种植牙的龈缘曲线设计方法 | |
| CN110811653A (zh) | 一种测量正畸过程中牙齿三维移动的方法 | |
| CN115619773B (zh) | 一种三维牙齿多模态数据配准方法及系统 | |
| EP4144324A1 (en) | Intelligent design method for digital model for oral digital impression instrument | |
| CN101950430A (zh) | 基于曲面断层片的三维牙齿重建方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091028 Termination date: 20140731 |
|
| EXPY | Termination of patent right or utility model |