CN109431632B - 一种模拟牙龈变形的方法 - Google Patents
一种模拟牙龈变形的方法 Download PDFInfo
- Publication number
- CN109431632B CN109431632B CN201811158336.5A CN201811158336A CN109431632B CN 109431632 B CN109431632 B CN 109431632B CN 201811158336 A CN201811158336 A CN 201811158336A CN 109431632 B CN109431632 B CN 109431632B
- Authority
- CN
- China
- Prior art keywords
- points
- point
- deformable
- deformation
- gum
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61C—DENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
- A61C7/00—Orthodontics, i.e. obtaining or maintaining the desired position of teeth, e.g. by straightening, evening, regulating, separating, or by correcting malocclusions
- A61C7/002—Orthodontic computer assisted systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Software Systems (AREA)
- Dentistry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Geometry (AREA)
- Epidemiology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Cosmetics (AREA)
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
本发明公开了一种模拟牙龈变形的方法。该方法的主要步骤包括:在牙龈模型上确定控制点、变形点和固定点;求取这三类点的拉普拉斯表达式;利用拉普拉斯变形公式和每个点的旋转矩阵,进而求出每个点拉普拉斯表达式旋转后的值,迭代多次,求取变形后点的坐标;将变形后点的坐标更新到牙龈模型中。本方法可以在牙龈模型变形过程中,保持牙龈模型本身的微小细节,并且变形后的牙龈模型更加真实。
Description
技术领域
本发明属于牙科隐形矫治领域,涉及一种模拟牙龈变形的方法。
背景技术
随着人们物质生活水平的不断提高,牙齿的整洁被看成是仪容仪表的重要部分,备受大众关注。在正畸市场逐渐从幼龄化转向成人化发展的趋势下,隐形矫治器满足了人们在正畸过程中美观的需求。由于每一个患者的牙龈形态各异,在计算机虚拟设计隐形牙套的过程中,需要对每个个体的牙龈变形进行模拟,从而实现个性化定制。在这个变形过程中,保留患者本身牙龈的细节信息,具有很大意义。
对于模型的变形研究中,主要包括三种方式:有限元模型、质点弹簧模型和基于拉普拉斯公式的变形。
其中,有限元模型具有精度高、能解决各种复杂问题和详细描述变形过程等优点,但是受计算机计算能力的限制,无法达到实时变形的效果,因此无法应用到牙龈变形中。
质点弹簧模型能够满足实时性的要求,并且在牙龈变形中具有广泛的应用。例如,中国专利,专利号201610231481.6提供了一种牙龈变形仿真的方法和装置,该专利就是基于质点弹簧模型所设计的,但是该专利的问题在于其稳定性差、不能保留牙龈模型上的细节。
基于拉普拉斯的变形应用最为广泛,它不改变模型的拓扑关系,只改变点的位置,性能稳定并且能够达到实时性要求。但是普通拉普拉斯变形后,细小结构部位出现一些毛刺,对精度和美观度有很大的影响。
还有将质点弹簧模型和拉普拉斯的变形结合到一起进行的,例如,中国专利,专利号201610878247.X提出的一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法,这种方法不仅仅需要计算每个点的拉普拉斯表示,还需要建立一个质点运动系统,计算量很大,并且这个方法并未将拉普拉斯表示的旋转考虑在内,变形后模型的细小部位精度可能会降低。
发明内容
为了解决背景技术中的问题,本发明基于拉普拉斯公式来模拟牙龈变形,在迭代过程中,增加了每个点对应的拉普拉斯表达的旋转矩阵,不仅计算量较小,并且使变形后的曲面更加光滑,并且保留了牙龈的细小形状,大大提高了模拟的效率和精度。
本发明的具体技术方案是:
本发明提出了一种模拟牙龈变形的方法,基本步骤如下:
【1】获取牙齿和牙龈的模型,并确定两者之间的位置关系;
【2】牙龈模型中特征点的选取,并记录特征点的初始坐标值;所述特征点包括控制点、可变形点以及固定点;
控制点为牙齿模型和牙龈模型相交的轮廓上的所有点,所述控制点随着牙齿移动做刚性移动;
根据所有控制点形成的轮廓,向外扩展N圈,将扩展N圈后区域内的所有点作为可变形点,所述可变形点随着牙齿移动而变形;
根据最外圈可变形点形成的轮廓,再次向外扩展一圈,扩展之后轮廓上的所有点为固定点;
【3】求取每个特征点初始坐标值所对应的拉普拉斯表达式;
【4】通过牙龈模型上点的拓扑结构,求取所有特征点所对应拉普拉斯矩阵;
【5】通过控制点与可变形点的连接关系,得到可变形点随着控制点运动的线性方程;
【6】获取变形后可变形点和控制点的坐标值;
【6.1】移动牙齿模型,牙龈模型上的所有控制点随之移动,从而得到控制点变形后坐标值;
【6.2】将控制点变形后坐标值更新到步骤【5】的线性方程中,得到变形后可变形点坐标值;
【6.3】利用步骤【6.2】得到的变形后可变形点的坐标值以及步骤【2】中可变形点的初始坐标值,通过坐标变换,得到可变形点所对应的拉普拉斯表达的旋转矩阵;
【6.4】将步骤【3】中的拉普拉斯表达式和步骤【6.3】中的拉普拉斯表达的旋转矩阵执行第一次旋转变换,得到可变形点对应的新的拉普拉斯表达式;
【6.5】通过步骤【6.4】得到的新的拉普拉斯表达式,重新求取可变形点的坐标值;
【6.6】重复执行步骤【6.3】至【6.5】至少一次,得到最终可变形点的坐标值。
具体地说,所述步骤【6.3】中任意可变形点所对应的拉普拉斯表达的旋转矩阵的具体计算过程是:
设任意一个可变形点的初始坐标值vi;利用步骤【6.2】得到可变形点对应的坐标v″i,进而得到v″i的转置v″i T,根据坐标变换,得到一个3×3的矩阵ti:
ti=vi·v″i T
对ti进行奇异值分解,得到
ti=uiεimi
其中,ui和mi是3×3阶酉矩阵,εi是半正定3×3阶对角矩阵;
那么所述任意可变形点对应的旋转矩阵为
Ti=mi*ui T
其中,ui T为ui的转置。
具体地说,所述步骤【6.6】中重复执行步骤【6.3】至【6.5】1-30次。
本发明的有益效果是:
本发明基于拉普拉斯表达式以及每个特征点对应的拉普拉斯表达的旋转矩阵所设计的方法,较之现有的模拟仿真方法,计算机只需进行拉普拉斯及拉普拉斯旋转变换,计算量减小,可以在200ms内执行完成,提高了模拟的效率。同时使变形后的牙龈曲面更加光滑,并且保留了牙龈的细小形状,提升了模拟精度。
附图说明
图1为本发明的流程框图;
图2为去掉了一颗牙齿的牙齿模型图;
图3为在图2的基础上标记了控制点、可变形点和固定点的牙齿模型图;
图4为在图3的基础上变形后所展示出的牙齿模型图。
图5为本发明方法模拟牙龈变形后的效果图。
具体实施方式
下面结合一个具体实例和附图对本发明的技术方案进行详细的介绍:
在正畸治疗过程中,由于每个正畸患者的牙齿移动方向和角度等信息不确定,需要根据每个患者的实际情况,实时的模拟牙龈变形,从而让正畸医师能够精准有效的进行治疗工作。
本发明设计一种模拟牙龈变形方法,参见图1,其具体步骤是:
步骤【1】通过计算机获取患者牙齿和牙龈的模型,并确定两者之间的位置关系;
步骤【2】牙龈模型中特征点的选取,并记录特征点的初始坐标值;所述特征点包括控制点、可变形点以及固定点;
控制点为牙齿模型和牙龈模型相交的轮廓上的所有点,本例中记录控制点的个数为k;控制点与牙齿具有相同的移动量和旋转量,随着牙齿做刚性运动。
选取控制点的边缘,向外扩展5-20圈,作为可变形点,本例中个数记录为m。
将可变形点向外扩一圈,获得固定点个数记为l。这三类点,构成了所有特征点,其中所有特征点个数记为t;t=k+m+l。
步骤【3】求取每个特征点初始坐标值所对应的拉普拉斯表达式;
设,牙龈模型上任意一个特征点vi的拉普拉斯表达式Δ(vi),其可以用一个三维向量来表示:
步骤【4】通过牙龈模型上点的拓扑结构,求取所有特征点所对应拉普拉斯矩阵;
LV=Δ(2)
其中:V为所有特征点初始坐标组成的t×3的行列式,Δ为所有特征点对应的拉普拉斯表达式,L为拉普拉斯矩阵。
步骤【5】通过控制点与可变形点的连接关系,得到可变形点随着控制点运动的线性方程;该线性方程具体如下:
其中:
Lf与步骤【4】中拉普拉斯矩阵相对应,取控制点对应的行,满足Lf∈L;
Ic是k×k的单位矩阵,k为控制点个数;
Inv是1×1的单位矩阵;1为固定点个数
V′表示所有特征点坐标值的行列式,第一次迭代与步骤【4】中值相同,每次迭代后更新;
Δf是与公式(2)中对应的可变形点的拉普拉斯表达式;
Vc′是k个控制点随着牙齿变形后对应的坐标值组成的k×3的行列式;
Vnv是固定点的坐标所组成的1×3的行列式;
步骤【6】获取变形后可变形点和控制点的坐标值;
步骤【6.1】移动牙齿模型,牙龈模型上的所有控制点随之移动,,得到控制点变形后坐标值Vc′;
步骤【6.2】将控制点变形后坐标值Vc′代入到步骤【5】的线性方程中,得到变形后可变形点坐标值V″;
步骤【6.3】利用步骤【6.2】得到的变形后可变形点的坐标值V″以及步骤【2】中可变形点的初始坐标值,通过坐标变换,得到可变形点所对应的拉普拉斯表达的旋转矩阵;
具体实施过程如下:设可变形点的初始坐标值vi,其中vi∈V;利用步骤【6.2】,得到与它对应的点坐标v″i,符合v″i∈V″;进而得到v″i的转置v″i,由坐标变换的知识,得到一个3×3的矩阵ti:
ti=vi·v″i T
对ti进行奇异值分解,得到
ti=uiεimi
其中,ui和mi是3×3阶酉矩阵,εi是半正定3×3阶对角矩阵。
那么这个点对应的旋转矩阵为
Ti=mi*ui T
其中ui T为ui的转置;
步骤【6.4】将步骤【3】中的拉普拉斯表达式和步骤【6.3】中的拉普拉斯表达的旋转矩阵执行第一次旋转变换,得到可变形点对应的新的拉普拉斯表达式;
变换后的拉普拉斯表达式为Δ′(vi)=Δ(vi)·Ti,
步骤【6.5】通过步骤【6.4】得到的新的拉普拉斯表达式,重新求取可变形点的坐标值;具体做法是:将公式(3)中的Δ(vi)替换为Δ′(vi),重新求取变形后点的坐标。
【6.6】重复执行步骤【6.3】至【6.5】多次(本实施中采用25次),得到最终可变形点的坐标值。
如图5所示,利用本发明进行拉普拉斯表达式的旋转后变形的更加光滑,可变形点与固定点的衔接更加自然,不会出现现有方法中会出现的明显“台阶”和“凹槽”结构,使得模拟效果更加贴近真实牙龈变形的情况,大大提高了模拟的精度和准确度。
Claims (3)
1.一种模拟牙龈变形的方法,其特征在于,包括以下步骤:
【1】获取牙齿和牙龈的模型,并确定两者之间的位置关系;
【2】选取牙龈模型中的特征点,并记录所有特征点的初始坐标值;所述特征点包括控制点、可变形点以及固定点;
控制点为牙齿模型和牙龈模型相交的轮廓上的所有点,所述控制点随着牙齿移动做刚性移动;
根据所有控制点形成的轮廓,向外扩展N圈,将扩展N圈后区域内的所有点作为可变形点,所述可变形点随着牙齿移动而变形;N取值为5-20;
根据最外圈可变形点形成的轮廓,再次向外扩展一圈,扩展之后轮廓上的所有点为固定点;
【3】求取每个特征点初始坐标值所对应的拉普拉斯表达式;
【4】通过牙龈模型上点的拓扑结构,求取所有特征点所对应拉普拉斯矩阵;
【5】通过控制点与可变形点的连接关系,得到可变形点随着控制点运动的线性方程;
【6】获取变形后可变形点和控制点的坐标值;
【6.1】移动牙齿模型,牙龈模型上的所有控制点随之移动,从而得到控制点变形后坐标值;
【6.2】将控制点变形后坐标值代入到步骤【5】的线性方程中,得到变形后可变形点坐标值;
【6.3】利用步骤【6.2】得到的变形后可变形点的坐标值以及步骤【2】中可变形点的初始坐标值,通过坐标变换,得到可变形点所对应的拉普拉斯表达的旋转矩阵;
【6.4】将步骤【3】中的拉普拉斯表达式和步骤【6.3】中的拉普拉斯表达的旋转矩阵执行第一次旋转变换,得到可变形点对应的新的拉普拉斯表达式;
【6.5】通过步骤【6.4】得到的新的拉普拉斯表达式,重新求取可变形点的坐标值;
【6.6】重复执行步骤【6.3】至【6.5】至少一次,得到最终可变形点的坐标值。
3.根据权利要求1所述的模拟牙龈变形的方法,其特征在于:所述步骤【6.6】中重复执行步骤【6.3】至【6.5】1-30次。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811158336.5A CN109431632B (zh) | 2018-09-30 | 2018-09-30 | 一种模拟牙龈变形的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811158336.5A CN109431632B (zh) | 2018-09-30 | 2018-09-30 | 一种模拟牙龈变形的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109431632A CN109431632A (zh) | 2019-03-08 |
CN109431632B true CN109431632B (zh) | 2020-09-11 |
Family
ID=65544526
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811158336.5A Active CN109431632B (zh) | 2018-09-30 | 2018-09-30 | 一种模拟牙龈变形的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109431632B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111973298A (zh) * | 2020-07-30 | 2020-11-24 | 崔军 | 一种用于口腔种植修复附着龈增宽的装置及方法 |
KR102520630B1 (ko) * | 2020-10-28 | 2023-04-12 | 주식회사 메디트 | 삼차원 구강 모델을 처리하는 방법 및 장치 |
US20230397972A1 (en) * | 2020-10-28 | 2023-12-14 | Medit Corp. | Method and device for processing three-dimensional oral cavity model and computer-readable recording medium |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101276483B (zh) * | 2008-04-14 | 2012-04-18 | 中山大学 | 一种实现平移敏感的Laplacian网格编辑方法 |
KR101295610B1 (ko) * | 2012-11-02 | 2013-08-12 | 주식회사 리얼오쏘 | 치아 데이터 컬러링 장치 |
CN104546151A (zh) * | 2013-10-09 | 2015-04-29 | 北京大学口腔医学院 | 数控激光自动化牙体预备方法及装备 |
CN105551081B (zh) * | 2016-02-05 | 2018-04-10 | 杭州美齐科技有限公司 | 一种虚拟牙龈三角网格构建和随动方法 |
CN105930566B (zh) * | 2016-04-10 | 2019-04-26 | 西安科技大学 | 牙龈变形仿真的方法和装置 |
FR3050374A1 (fr) * | 2016-04-22 | 2017-10-27 | H43 Dev | Procede de prediction d'une situation dentaire |
CN106204739A (zh) * | 2016-07-19 | 2016-12-07 | 上海爱圣美科技有限公司 | 一种三维牙齿模型中牙龈重建的方法 |
CN106204730B (zh) * | 2016-07-19 | 2019-04-26 | 上海爱圣美科技有限公司 | 一种三维牙齿模型中牙龈随动的方法 |
CN106504318B (zh) * | 2016-10-19 | 2019-08-23 | 青岛达芬奇科技有限公司 | 一种用于无托槽正畸的牙龈生成方法 |
CN108053484A (zh) * | 2017-12-12 | 2018-05-18 | 武汉大学 | 一种基于法向平均的常平均曲率曲面构建方法 |
-
2018
- 2018-09-30 CN CN201811158336.5A patent/CN109431632B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109431632A (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103208133B (zh) | 一种图像中人脸胖瘦的调整方法 | |
US11666416B2 (en) | Methods for simulating orthodontic treatment | |
CN106504329B (zh) | 一种基于牙体长轴的质点弹簧模型的牙龈变形仿真方法 | |
US11642195B2 (en) | Visual presentation of gingival line generated based on 3D tooth model | |
CN109431632B (zh) | 一种模拟牙龈变形的方法 | |
CN106023288B (zh) | 一种基于图像的动态替身构造方法 | |
US8830269B2 (en) | Method and apparatus for deforming shape of three dimensional human body model | |
US20210186659A1 (en) | 2d-to-3d tooth reconstruction, optimization, and positioning frameworks using a differentiable renderer | |
CN105030347A (zh) | 基于牙根信息的数字化正畸排牙方法及基于该方法的托槽转移装置设计方法 | |
CN110619681B (zh) | 一种基于欧拉场形变约束的人体几何重建方法 | |
Zheng et al. | A novel 3D morphing approach for tooth occlusal surface reconstruction | |
JP2003044873A (ja) | 顔の3次元モデルの作成方法及びその変形方法 | |
CN105389444B (zh) | 一种个性化种植牙的龈缘曲线设计方法 | |
CN106504318B (zh) | 一种用于无托槽正畸的牙龈生成方法 | |
CN102236913A (zh) | 一种体网格变形方法 | |
Zhang et al. | Computer-aided design of dental inlay restoration based on dual-factor constrained deformation | |
Badin et al. | Three-dimensional linear modeling of tongue: Articulatory data and models | |
WO2011054177A1 (zh) | 设计中间全颌牙齿数字模型的方法 | |
Ma et al. | Animating visible speech and facial expressions | |
CN110163957A (zh) | 一种基于唯美人脸程序的表情生成系统 | |
CN113379920A (zh) | 一种鼻部整形区域重构方法 | |
CN111839728A (zh) | 一种鼻整形仿真方法 | |
CN116824032A (zh) | 产生牙龈三维数字模型的方法 | |
CN116543096A (zh) | 产生牙颌三维数字模型的方法 | |
JP2005074136A (ja) | 顔形状のモデリング方法及び顔形状モデリング用プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |