发明内容
本发明的目的提供了一种设计中间全颌牙齿数字模型的方法,该方法克服已有的无托槽隐形矫治器只是对需要矫治牙齿从初始状态到目标状态的简单分解的不足,通过加强实际矫治过程中的支抗,把需矫治的牙齿在支抗牙与移动牙之间的转换,并根据该转换得到一系列中间初始全颌牙齿数字模型。同时本发明还基于成熟的数字三维网格处理算法,利用图形硬件的加速功能,提高处理精细牙齿数字模型的速度和稳定性。
为了达到上述目的,本发明采用包括以下5个步骤的技术方案:
1.交互式调整牙齿位置,获取牙齿从初始位置到目标位置的运动量。为了实现此目的,包含3个子步骤:建立牙齿局部运动系、交互式运动操作、添加关键帧。
a)建立牙齿局部坐标系
牙齿的三维数据是和全颌牙齿数字模型整体保持一致的,这使得牙齿只能参照全颌牙齿数字模型整体的坐标系进行移动、旋转,这不符合数字化口腔正畸要求每一颗牙齿可以独立相对自身进行移动、旋转的需求。
建立每一颗牙齿的局部坐标系使得每一颗牙齿可以相对于自身发生移动、旋转运动。
建立牙齿局部坐标系的具体步骤是:
●计算牙齿三维数字网格数据世界坐标(即牙齿原始空间x、y、z坐标值)各维数据各自的平均值X、Y、Z,此平均值即为牙齿三维数字网格的中心点。
●计算协方差矩阵。
公式: 用来计算两个坐标维之间的协方差。
其中:X和Y代表X、Y这两个坐标方向上的平均值,由前面步骤计算。
公式:Cm*n=(ci,j,ci,j=cov(Dimi,Dimj))用于计算整个协方差矩阵。具体展开即为:
●计算协方差矩阵C的特征向量和特征值,根据特征值的大小进行排序,按照特征值大小的递减顺序分别定义特征值对应的特征向量为主轴、次轴、3轴。
通过这3个步骤,就计算得到一个牙齿三维数字网格的局部坐标系:中心点和3个彼此正交的坐标轴。
b)交互式运动操作
通过牙齿周围的运动控件来交互式操作牙齿的运动。具体表现为:在牙齿的周围绘制6个带有双向箭头示意的控件图标,分别表示沿主轴正/负方向平移、沿次轴正/负方向平移、沿3轴正/负方向平移以及绕主轴正方向顺/逆时针旋转、绕次轴正方向顺/逆时针旋转、绕3轴正方向顺/逆时针旋转。
当鼠标发生拖拽动作的时候,通过拾取当前屏幕坐标点并且应用OpenGL矩阵反向计算出此屏幕点当前在三维空间的坐标p(x,y,z)。
拾取点三维空间坐标p用于比对各个运动控件查看是哪一个运动控件发生作用,根据此控件代表的运动意义,以及鼠标在屏幕上运动的距离计算牙齿的运动。
c)添加关键帧
关键帧定义为全颌牙齿数字模型状态的一个快照,其中若干牙齿发生运动。关键帧添加是牙齿位置调整中的关键功能,通过各个关键帧来严格限定牙齿在从初始位置运动到目标位置过程中经历的特殊位置来规划整个牙齿的运动路径,这是本发明可以处理复杂病例的核心技术之一。
为了实现这个目的,本发明提供的牙齿交互式位置调整方法包含了保存当前位置到关键帧以及关键帧的选择、删除等操作。
2.平滑插值分解牙齿运动量为渐进运动步骤。
一个合适的作用周期为2周的无托槽隐形矫治器,其对应牙齿的单步最大运动量应该为0.5mm,该数值可根据患者的年龄,颚骨密度等因素做适当的修正。因此平滑插值分解获得渐进运动步骤为了保证牙齿的单步最大运动量不超过0.5mm。
平滑插值分解牙齿运动量包括计算机自动插值和手动微调两个主要步骤。
a)计算机自动插值
计算机自动插值基于满足三个限制参数:单步最大平移量、单步最大旋转角、单步最大运动量,将运动量自动分解为初始渐进运动步骤的方法。具体的实现方式为:
平移运动插值:直接线性插值
公式 用于分解平移运动,其中,p0为牙齿初始位置的坐标,pd为牙齿目标位置的坐标,n是分解的步骤总数。
●旋转运动插值:
为了克服直接对角度进行插值带来的不平滑以及万向锁等问题,本发明采用四元数的球面线性插值方法,此方法可以在两个四元数之间平滑插值,并且避免万向锁问题。
公式slerp(q0,qd,t)=q0*((q0)-1*qd)t用于计算四元数的球面线性插值,其中q0和qd分别是牙齿在初始位置和目标位置处的旋转角四元数值,t是插值参数,相当于平移运动插值的
●运动数据校验:
运动数据校验是保证单颗牙齿的单步运动量不超过人体承受限度,使用这个模型制备的无托槽隐形矫治器在临床上安全、有效。
在牙齿表面选择4个以上(包括4个)参考点,计算这些参考点在此单步运动中的运动量,以全部参考点的运动量都不超过单步最大运动量作为计算机自动插值生成初始渐进步骤的标准。参考点的选择原则为相邻两步骤间参考点的最大运动量能尽量接近牙齿的最大运动量。例如依据牙齿局部坐标系3个坐标轴,选取接近3个坐标轴两端的顶点,即距离牙齿中心远端的6个顶点作为参考点。
b)手动微调
手动交互式微调为设计人员提供修改初始渐进运动步骤的机会,使运动数据分解能够根据设计人员的需求进行调整生成渐进运动步骤。
为了实现这个目的,需要将前述计算机自动插值的数据进行显示,提供修改接口。
具体的修改操作实现需要经过当前步骤数据修改、当前步骤运动数据校验、前驱步骤数据修改、前驱步骤运动数据校验、后继步骤数据修改、后继步骤运动数据校验等环节。核心思想是在修改当前单步数据导致前驱、后继单步数据发生变化的时候,需要保证当前以及前驱、后继单步运动数据均满足运动数据校验的要求。
3.以需矫治的牙齿在支抗牙与移动牙之间的转换为原则,组合步骤(2)中的牙齿渐进运动步骤得到一系列中间初始全颌牙齿数字模型。
中间全颌牙齿数字模型对应着数字化口腔正畸方法的一个矫治器,矫治方案的制定就是将需矫治牙齿的各个渐进运动步骤结合病例实际情况合理安排到不同的中间初始全颌牙齿数字模型中。
一种基本的矫治方案是所有需矫治的牙齿采用齐头并进式位置调整。该矫治方案需要的中间全颌牙齿数字模型是最少的,设计的矫治周期也是最短的。但是由于力的相互作用,使得真实的矫治过程往往会偏离设计的矫治过程。从而导致在矫治一段时间后,需要对患者重新采集全颌牙齿数据,重新设计矫治方案,最终导致完成整个矫治过程时间更长,花费的人力、物力更多。当牙颌内需矫治的牙齿数量大于无需矫治的牙齿数量时,上述情况出现的频率非常高。所有发生运动的牙齿采用齐头并进式位置调整的矫治方案所以会出现上述情况,主要是因为没有考虑到实际矫治过程中支抗的重要性。
为了加强实际矫治过程中的支抗,减少重复取牙颌模型及重复做矫治方案的次数,本发明提供了组合需矫治牙齿各个渐进运动步骤的接口,包含添加中间初始全颌牙齿数字模型、删除中间初始全颌牙齿数字模型,以及安置需矫治牙齿单个渐进运动步骤到任意中间初始全颌牙齿数字模型的操作。
4.安装附件。
附件主要设计用来增强力的作用,处理一些特殊的牙齿运动,辅助矫治器进行牙齿的运动。前牙的升高、压低,锥形或者原型牙齿的旋转,拔牙病例中间隙的闭合以及牙齿大距离的平移是需要附件进行辅助矫治的主要应用场合。
附件安装涉及矫治器附件容纳和牙齿表面附件安装。一方面,医生根据矫治方案中指明的附件种类、大小等信息在病例需要进行附件安装的牙齿表面安装附件。另一方面,矫治器必须包含可以容纳附件的附件容纳仓,附件容纳仓紧紧包裹着安装在牙齿上的附件,加大矫治器对牙齿施加的力。
为了实现这个目的,需要进行附件三维数字网格和宿主牙齿三维数字网格的融合即布尔合并操作。郭开波,张李超等在《STL模型布尔运算的实现》一文中提出的基于相交环检测的布尔操作方法可以很好的应用在牙齿模型和附件模型的网格合并中。具体的实现步骤如下:
a)读入两个实体,拓扑重构,建立连接关系,建立封闭表面的信息。
b)相交性测试,若有面片相交,则转c,无面片相交则转f。
c)求取交线,跟踪提取交线环。
d)对被相交的三角形进行二次三角划分,利用交线环剖分相交表面。
e)判断剖分得到的子表面相对于另一个实体的位置关系。
f)判断所有非相交表面相对于另一个实体的包含关系。
g)利用布尔算子实现布尔运算。
5.虚拟牙龈自由变形,生成中间全颌牙齿数字模型
牙齿位置的变动会带动牙齿周围的牙龈组织发生变化。牙龈变形计算所得到的更加符合牙齿牵引作用的牙龈数据能提高矫治器在使用过程中的舒适程度,减小矫治器对牙龈组织的挤压。同时,牙龈随着牙齿位置运动而变形的过程能够提高医生与患者在矫治方案设计交流的方便,能够直观地向患者展示整个矫治方案。
虚拟牙龈自由变形是根据矫治方案阶段中的牙齿运动信息对虚拟牙龈网格数据进行变形计算,使中间全颌牙齿数字模型包含变形后的牙龈数据。
为了实现这个目的,本发明采用基于Laplacian算子的模型自由变形算法。将虚拟牙龈分为不同的变形区域,一个变形区域包括一颗牙齿临近牙龈区域与相邻两颗牙齿的边界作为变形区域的固定边界。使用牙齿运动架控制牙龈边界,进而控制牙龈局部区域变形。牙龈整体变形采用局部区域依次变形,最终完成牙龈变形。
实施例1上牙颌间隙关闭单颌矫治案例中的全颌牙齿数字模型生成。
需矫治牙齿为右尖牙(UR3)、右第二切牙(UR2)、右切牙(UR1)、左切牙(UL1)、左第二切牙(UL2)、左尖牙(UL3)。
1.打开软件,导入完成需关闭间隙的牙齿分割修补的全颌牙齿数字模型,如说明书附图:图1上牙颌间隙关闭单颌矫治案例中的完成牙齿分割修补的全颌牙齿数字模型(包括颌面图、唇面图、左颊面图、右颊面图)。
2.通过建立牙齿局部运动系、交互式运动操作、添加关键帧,调整牙齿位置,使牙齿间隙关闭,排列整齐,获取牙齿从初始位置到目标位置的运动量。如说明书附图:图2上牙颌间隙关闭单颌矫治案例中完成位置调整的全颌牙齿数字模型。
3.设定平滑插值分解牙齿单步最大移动量0.5mm,单步最大旋转量5度,单步最大运动量0.5mm,通过平滑插值及交互式调整,分解牙齿运动量为渐进运动步骤。
4.在第三步牙齿渐进运动的基础上,考虑实际矫治中的支抗作用,得到一系列矫治过程中中间初始全颌牙齿数字模型。如说明书附图:图3上牙颌间隙关闭单颌矫治案例中从最初位置到最终位置共7步的中间初始全颌牙齿数字模型(颌面图),图4上牙颌间隙关闭单颌矫治案例中从最初位置到最终位置共7步的中间初始全颌牙齿数字模型(右颊面图),每一步中间初始全颌牙齿数字模型中牙齿位置变化情况见表1-7。
表1第一步中间初始全颌牙齿数字模型中牙齿位置变化情况
第一步 |
UR3 |
UR2 |
UR1 |
UL1 |
UL2 |
UL3 |
主轴移动 |
0.00 |
-0.45 |
-0.46 |
-0.41 |
0.00 |
-0.18 |
次轴移动 |
0.00 |
0.06 |
0.00 |
-0.10 |
0.00 |
-0.25 |
3轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
主轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
次轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
3轴旋转 |
0.00 |
-1.63 |
-1.26 |
0.00 |
0.00 |
0.00 |
表2第二步中间初始全颌牙齿数字模型中牙齿位置变化情况
第二步 |
UR3 |
UR2 |
UR1 |
UL1 |
UL2 |
UL3 |
主轴移动 |
-0.33 |
-0.45 |
0.00 |
-0.42 |
-0.32 |
0.00 |
次轴移动 |
0.15 |
0.07 |
0.00 |
-0.10 |
-0.21 |
0.00 |
3轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
主轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
次轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
3轴旋转 |
2.26 |
-1.63 |
0.00 |
0.00 |
0.32 |
0.00 |
表3第三步中间初始全颌牙齿数字模型中牙齿位置变化情况
第三步 |
UR3 |
UR2 |
UR1 |
UL1 |
UL2 |
UL3 |
主轴移动 |
0.00 |
-0.45 |
-0.46 |
0.00 |
-0.34 |
0.00 |
次轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
-0.21 |
0.00 |
3轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
主轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
次轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
3轴旋转 |
0.00 |
-1.63 |
-1.26 |
0.00 |
0.32 |
0.00 |
表4第四步中间初始全颌牙齿数字模型中牙齿位置变化情况
第四步 |
UR3 |
UR2 |
UR1 |
UL1 |
UL2 |
UL3 |
主轴移动 |
-0.34 |
0.00 |
-0.46 |
-0.42 |
0.00 |
0.00 |
次轴移动 |
0.14 |
0.00 |
0.01 |
-0.10 |
0.00 |
0.00 |
3轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
主轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
次轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
3轴旋转 |
2.26 |
0.00 |
-1.26 |
0.00 |
0.00 |
0.00 |
表5第五步中间初始全颌牙齿数字模型中牙齿位置变化情况
第五步 |
UR3 |
UR2 |
UR1 |
UL1 |
UL2 |
UL3 |
主轴移动 |
0.00 |
-0.45 |
-0.46 |
-0.42 |
0.00 |
0.00 |
次轴移动 |
0.00 |
0.10 |
0.02 |
-0.10 |
0.00 |
0.00 |
3轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
主轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
次轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
3轴旋转 |
0.00 |
-1.63 |
-1.26 |
0.00 |
0.00 |
0.00 |
表6第六步中间初始全颌牙齿数字模型中牙齿位置变化情况
第六步 |
UR3 |
UR2 |
UR1 |
UL1 |
UL2 |
UL3 |
主轴移动 |
-0.34 |
0.00 |
-0.46 |
-0.42 |
-0.34 |
0.00 |
次轴移动 |
0.12 |
0.00 |
0.03 |
-0.10 |
-0.21 |
0.00 |
3轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
主轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
次轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
3轴旋转 |
2.26 |
0.00 |
-1.26 |
0.00 |
0.32 |
0.00 |
表7第七步中间初始全颌牙齿数字模型中牙齿位置变化情况
第七步 |
UR3 |
UR2 |
UR1 |
UL1 |
UL2 |
UL3 |
主轴移动 |
0.00 |
-0.44 |
-0.46 |
-0.42 |
-0.34 |
0.00 |
次轴移动 |
0.00 |
0.11 |
0.04 |
-0.10 |
-0.21 |
0.00 |
3轴移动 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
主轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
次轴旋转 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
0.00 |
3轴旋转 |
0.00 |
-1.63 |
-1.26 |
0.00 |
0.32 |
0.00 |
5.在中间初始全颌牙齿数字模型上左右第二切牙上安装用于增强矫治力的附件模型。图5上牙颌间隙关闭单颌矫治案例中在中间初始全颌牙齿数字模型的左右第二切牙上安装附件模型。
6.设置牙齿位置变动时,虚拟牙龈自由变形的区域。如说明书附图:图6上牙颌间隙关闭单颌矫治案例中设置虚拟牙龈自由变形区域。
7.输出中间全颌牙齿数字模型,保存为三角网格模型。如说明书附图:图7上牙颌间隙关闭单颌矫治案例中第七步的中间全颌牙齿数字模型三角网格模型。