CN102930091B - 一种面向牙齿牙龈的高效网格修补方法 - Google Patents
一种面向牙齿牙龈的高效网格修补方法 Download PDFInfo
- Publication number
- CN102930091B CN102930091B CN201210409293.XA CN201210409293A CN102930091B CN 102930091 B CN102930091 B CN 102930091B CN 201210409293 A CN201210409293 A CN 201210409293A CN 102930091 B CN102930091 B CN 102930091B
- Authority
- CN
- China
- Prior art keywords
- tooth
- gum
- border
- hole
- hole border
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 210000000515 tooth Anatomy 0.000 claims description 111
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000002485 combustion reaction Methods 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 2
- 238000000513 principal component analysis Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 10
- 238000013461 design Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011960 computer-aided design Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 1
- 208000010641 Tooth disease Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Landscapes
- Image Generation (AREA)
- Prostheses (AREA)
Abstract
本发明公开了一种面向牙齿牙龈的高效网格修补方法,包括:首先,构造牙齿和牙龈信息缺失部分的孔洞边界,将该孔洞边界投影到二维平面上;其次,采用约束Delaunay三角化算法对投影后的孔洞边界内部进行三角化;再次,采用均值坐标原理将三角化后的网格反投影回原三维空间;最后,采用二阶的拉普拉斯算子优化上述修补后的网格曲面。本发明采用PCA算法计算牙齿孔洞边界的投影平面,采用弦长参数化方法将牙龈孔洞边界投影到单位正方形上,减少孔洞边界的变形程度,将三角化的工作降维,大大缩短了三角化的时间,采用均值坐标原理反投影保证了三角网格的质量,最后采用二阶的拉普拉斯算子进行光顺,保证了孔洞边界的一阶连续性,算法简单,高效。
Description
技术领域
本发明涉及几何修补领域,特别涉及一种面向牙齿牙龈的高效网格修补方法。
背景技术
牙齿隐形正畸是一个涉及多学科的交叉应用,首先通过计算机辅助设计完成对数字化数字牙网格模型的编辑,然后通过材料学的相关技术生成医学牙套,用于口腔正畸。数字牙网格模型由几何信息和连接关系两部分构成,几何信息给出了各个顶点在三维空间的位置坐标,连接关系将这些位置信息连接成特殊的多边形曲面。
三维扫描重建后的牙颌模型是一个整体,无法满足对各个牙齿进行重排以及矫治方案设计的需求,因此网格编辑的第一步是将每颗牙齿从整体牙颌网格模型中分割出来。由于牙齿和牙龈存在遮挡,三维扫描仪不能获得相邻的两牙齿之间的邻面数据,当使用网格分割算法分割出牙齿模型以后,牙齿的邻面数据和牙齿之间的牙龈数据是不完整的,因此数字牙网格模型的编辑设计的一个核心问题是如何对分割后的牙齿邻面数据和牙齿之间的牙龈数据进行自动修补。
几何修补的目标是重建孔洞区域的曲面。一种好的孔洞修补方法应该是自动的,运行速度较快(最好满足交互速度),与周围的网格差异(网格密度、形状等)很小,适用于任意网格的任意孔洞。本发明的方法主要是针对于牙齿牙龈网格,该网格的特点是边界不完整或部分边界无邻域。
早期的孔洞修补算法首先检测孔洞的边界,然后对这些边界顶点进行三角化,但这种算法对于较大的孔洞就比较粗糙,不适用。对于大孔洞,现有的算法主要是在三角形内插入顶点,进行细分。PS96(参见PFEIFLER.,SEIDEL H.-P.:Triangular b-splines for blending and filling of polygonalholes.In GI’96:Proceedings of the conference on Graphics interface’96(Toronto,Ontario,Canada,1996),Canadian Information Processing Society,pp.186-193.)通过最小化曲率函数使插入的顶点尽量满足Delaunay原则,Lie03(参见LIEPA P.:Filling holes in meshes.In Proc.Euro-graphicsSymposium on Geometry Processing(2003),pp.200-205.)使新生成的网格满足采样密度以及法向连续性,WWP10(参见WEI M.,WU J.,PANG M.:An integrated approach to filling holes in meshes.In Pro ceedings of the 2010International Conference on Artificial Intelligence and ComputationalIntelligence(Washington,DC,USA,2010),IEEE Computer Society,pp.306-310.)考虑网格内部三角形的角度、二面角和密度。
发明内容
本发明提供了一种面向牙齿牙龈的高效网格修补方法,能够高效的生成高质量的网格,并且与原始网格自然拼接。
一种面向牙齿牙龈的高效网格修补方法,包括步骤:
(1)选取牙齿网格模型边界上两点作为牙齿邻面孔洞边界的起止顶点,构造牙齿邻面孔洞边界;
(2)采用主元分析算法计算出牙齿邻面孔洞边界的第三主元方向;
(3)将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上,然后将平面旋转至XY平面上;
(4)采用受约束的Delaunay三角化算法对XY平面上的投影孔洞边界及其内部进行三角化;
(5)采用均值坐标原理将三角化后的孔洞反投影回原三维空间,得到修补后的网格曲面;
(6)采用二阶的拉普拉斯光顺算法优化上述修补后的网格曲面,得到光顺后的牙齿网格模型;
(7)采用迪杰斯特拉算法寻找相邻牙齿对应顶点之间的最短路径,构造牙龈孔洞边界;
(8)采用弦长参数化方法将牙龈孔洞边界投影到二维平面的单位正方形上;
(9)对投影到二维平面的单位正方形上的牙龈孔洞边界采用步骤(4)~(6)得到修补后的牙龈网格模型。
本发明的面向牙齿牙龈的高效网格修补方法的步骤(3)采用了主元分析算法(该算法具体可参见J.Shlens.A Tutorial on Principal ComponentAnalysis.Institute for Nonlinear Science,UCSD,2005.)计算出牙齿邻面孔洞边界的第三主元方向,进而确定孔洞边界投影平面。
投影到二维平面后,在步骤(4)中还采用了约束Delaunay三角化算法(该算法具体可参见SHEWCHUK,J.R.:Triangle:engineering a 2Dquality mesh generator and Delaunay triangulator,In Proc.Of First Workshopon Applied Computational Geometry,(1996),124-133.)对投影后的孔洞边界进行三角化,大大降低三角化的时间,在该算法的实现上使用了Triangle库(该库具体参见http://www.cs.cmu.edu/~quake/tripaper/triangle0.html,由Jonathan Richard Shewchuk建立,来自Computer Science DivisionUniversity of California at Berkeley,邮箱为jrscs.berkeley.edu)。
得到二维的三角化的结果后,步骤(5)中采用了均值坐标原理(该原理具体可参见FLOATER,M.S.:Mean-value coordinates.Com--puter Aided Geometric Design,20,(2003),19-27.)将三角化后的孔洞反投影回原三维空间,该算法具有保形性,因此能保持二维空间三角形的Delaunay特性。
由于经过上述算法生成的网格模型在边界处会比较生硬,在步骤(6)中本发明还采用了二阶拉普拉斯光顺算法(该算法具体可参见Sorkine O.,Cohen-Or D.Least-squares meshes[C].Shape Modeling Applications,2004:191-199)对原始边界附近以及新生成的网格区域进行优化,得到光滑的经邻面修补的牙齿网格模型。
修补后相邻牙齿邻面数据后,需要对两邻牙之间的牙龈数据进行修补,在步骤(7)中,采用迪杰斯特拉算法计算相邻牙齿对应顶点之间的最短路径,删除两邻牙间最短路径以外的冗余网格,构造牙龈孔洞边界。
由于将牙龈孔洞边界直接投影到最近似平面上可能会产生自交,因此采用弦长参数化方法(该算法具体可参见FLOATER,M.S.,HORMANN,K.:Surface parameterization:a tutorial and survey.In Advances onMulti-resolution in Geometric Modeling,Springer-Verlag,Heidelberg,M.S.F.N.Dodgson and M.Sabin,Eds.,(2004).)将牙龈孔洞边界直接投影到二维平面的单位正方形上。
最后,对牙龈孔洞执行(4)~(6)的步骤得到光滑的经修补的牙龈网格模型。
具体的,步骤(1)中,所述选取牙齿网格模型边界上两点作为牙齿邻面孔洞边界的起止顶点,构造牙齿邻面孔洞边界,包括步骤:
(1.1)用鼠标选择牙网格模型边界上的两点作为牙齿邻面孔洞边界的起止顶点;
(1.2)连接以上选取的起始顶点和截止顶点,离散该连线,离散长度为牙网格模型边界的网格距离平均值,离散后的连线与牙网格模型部分边界构成一个封闭的邻面孔洞边界。
具体地,步骤(2)中,所述的采用主元分析算法计算出牙齿邻面孔洞边界的第三主元方向,包括步骤:
(2.1)将边界上的顶点坐标表示成一个3×n的矩阵,其中n为边界顶点的个数,矩阵列向量为顶点的(x,y,z)坐标;
(2.2)在每个矩阵行向量上减去该行向量的平均值得到矩阵X,XT为X的转置;
(2.3)对XXT进行特征分解,求取特征向量以及所对应的特征值,最小特征值对应的特征向量即为第三主元方向。
具体的,步骤(3)中,所述的将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上,然后将平面旋转至XY平面上,包括步骤:
(3.1)将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上,该平面最近似于孔洞边界的拟合平面;
(3.2)计算以上平面与XY平面之间的旋转矩阵,将孔洞边界顶点的投影坐标乘以旋转矩阵即得到顶点在XY平面上的坐标。
具体的,步骤(4)中,所述的采用约束Delaunay三角化算法对XY平面上的投影孔洞边界及其内部进行三角化,包括步骤:
(4.1)对XY平面上的投影孔洞边界顶点进行Delaunay三角化,形成具有新的连接关系的网格;
(4.2)孔洞边界上的边在新的网格中可能并不存在,在该网格中插入,同时去除相交的边,对该边邻域重新三角化,使其满足Delaunay规则;
(4.3)去除孔洞边界外的三角形;
(4.4)在三角形内插入新的顶点,重新进行Delaunay三角化,并使其满足最小角度和最大面积约束。
具体的,步骤(5)中,所述采用均值坐标原理将三角化后的孔洞反投影回原三维空间,得到修补后的网格曲面,包括步骤:
(5.1)若顶点属于边界顶点,在将该顶点对应投影回原三维坐标顶点;
(5.2)若顶点属于边界内顶点,即新添加的顶点,则利用均值坐标计算出其三维坐标,方法如下:
对于任意顶点v0,其坐标可以用其一环邻域顶点vi坐标表示:
其中,k为v0一环邻域内点的个数,λi为单位化后的顶点vi的权重系数,wi为单位化前的顶点vi的权重系数,αi为边v0vi和v0vi+1的夹角,αi-1为边v0vi-1和v0vi的夹角,顶点vi-1,vi,vi+1呈逆时针顺序。
具体的,步骤(6)中,所述采用二阶的拉普拉斯光顺算法优化上述修补后的网格曲面,得到光顺后的数字牙网格模型,包括步骤:
(6.1)对于任意顶点v0,拉普拉斯算子可表示为
其中,vi为v0的一环邻域顶点,ω(v0,vi)为vi相对于v0的权重,
则高阶的拉普拉斯算子可以表示为这里二阶拉普拉斯算子中k=2;
(6.2)建立如下方程组进行求解,得到经过光顺后的顶点坐标:
其中,P=(v1,....,vp)是孔洞内部的可自由移动的顶点,F=(f1,...,fF)包含边界顶点和边界外围一环邻域内的点,这些点为固定点,IF为单位矩阵。
具体的,步骤(7)中,所述采用迪杰斯特拉算法寻找相邻牙齿对应顶点之间的最短路径,构造牙龈孔洞边界,包括步骤:
(7.1)采用迪杰斯特拉算法寻找相邻牙齿对应顶点(牙齿邻面孔洞边界构造时选取的顶点)之间的最短路径;
(7.2)删除相邻牙齿之间,两条最短路径之上的冗余网格,两条最短路径与两邻牙的邻面边界构成了一条封闭的牙龈孔洞边界;
具体的,步骤(8)中,所述的采用弦长参数化方法将牙龈孔洞边界投影到二维平面的单位正方形上,包括步骤:
(8.1)采用弦长参数化方法将牙龈孔洞边界直接投影到二维平面的单位正方形上。弦长参数化公式如下:
(ti+1-ti)/(ti-ti-1)=||xi+1-xi||/||xi-xi-1||
其中,ti为曲线参数,xi为第i个顶点的坐标,可表示为(x(ti),y(ti),z(ti))。
本发明提出了一种面向牙齿牙龈的高效网格修补算法。该算法大致可分为四个步骤:首先,构造牙齿和牙龈信息缺失部分的孔洞边界,将该孔洞边界投影到二维平面上;其次,采用约束Delaunay三角化算法对投影后的孔洞边界内部进行三角化;再次,采用均值坐标原理将三角化后的网格反投影回原三维空间;最后,采用二阶的拉普拉斯算子优化上述修补后的网格曲面。本发明采用PCA算法计算牙齿孔洞边界的投影平面,采用弦长参数化方法将牙龈孔洞边界投影到单位正方形上,尽可能减少孔洞边界的变形程度,将三角化的工作降维,大大缩短了三角化的时间,采用均值坐标原理反投影保证了三角网格的质量,最后采用二阶的拉普拉斯算子进行光顺保证了孔洞边界的一阶连续性,算法简单,高效。
附图说明
图1为本发明面向牙齿牙龈的高效网格修补方法的流程示意图;
图2为本发明面向牙齿牙龈的高效网格修补方法中牙齿网格模型的邻面示意图;
图3为本发明面向牙齿牙龈的高效网格修补方法中构造的牙齿邻面孔洞边界;
图4为本发明面向牙齿牙龈的高效网格修补方法中CDT算法中一个吉他的边界轮廓示意图;
图5为本发明面向牙齿牙龈的高效网格修补方法中CDT算法中对吉他顶点Delaunay三角化后的结果示意图;
图6为本发明面向牙齿牙龈的高效网格修补方法中CDT算法中插入吉他边界轮廓原本存在的边后的结果示意图;
图7为本发明面向牙齿牙龈的高效网格修补方法中CDT算法中去除吉他边界轮廓外的三角形后的结果示意图;
图8为本发明面向牙齿牙龈的高效网格修补方法中CDT算法中对内部三角形进行Refine后的结果示意图;
图9为本发明面向牙齿牙龈的高效网格修补方法中均值坐标中顶点的一环邻域示意图;
图10为本发明面向牙齿牙龈的高效网格修补方法中相邻牙齿对应顶点之间最短路径示意图;
图11为本发明面向牙齿牙龈的高效网格修补方法中构造的牙龈孔洞边界示意图。
具体实施方式
下面结合附图对本发明邻面修补方法进行详细说明。
一种如图1所示的面向牙齿牙龈的高效网格修补方法,包括步骤:
(1)选取牙齿网格模型边界上两点作为牙齿邻面孔洞边界的起止顶点,构造牙齿邻面孔洞边界。
牙齿邻面孔洞边界的构造进一步包括步骤:
(1.1)图2是牙齿网格模型的邻面示意图,虚线构成了邻面的部分边界,用鼠标选择牙网格模型边界上的两点作为牙齿邻面孔洞边界的起止顶点,在图2中分别为A点、B点;
(1.2)如图3所示,连接以上选取的起始顶点A和截止顶点B,离散该连线,离散长度为牙网格模型边界的网格距离平均值,离散后的连线与牙网格模型部分边界构成一个封闭的邻面孔洞边界。
(2)采用主元分析算法计算出牙齿邻面孔洞边界的第三主元方向,进一步包括步骤:
(2.1)将边界上的顶点坐标表示成一个3×n的矩阵,其中n为边界顶点的个数,矩阵列向量为顶点的(x,y,z)坐标,矩阵表示为:
(2.2)在每个矩阵行向量上减去该行向量的平均值得到矩阵X,XT为X的转置:
(2.3)对XXT进行特征分解,求取特征向量以及所对应的特征值,最小特征值对应的特征向量即为第三主元方向。
(3)将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上,然后将平面旋转至XY平面上,进一步包括步骤:
(3.1)将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上,该平面最近似于孔洞边界的拟合平面;
空间点P0(x,y,z)在空间平面(normal,offset)上的投影坐标为P=P0-normal*((P0-offset)*normal)。
(3.2)计算以上平面与XY平面之间的旋转矩阵,将孔洞边界顶点的投影坐标乘以旋转矩阵即得到顶点在XY平面上的坐标。
以上平面与XY平面之间的夹角等于两平面的法向量夹角θ,然后运用Rodrigues旋转公式,XY平面上的点坐标为其中为两平面的法向量的叉乘。
(4)采用受约束的Delaunay三角化算法对XY平面上的投影孔洞边界及其内部进行三角化,进一步包括如下步骤:
(4.1)对XY平面上的投影孔洞边界顶点进行Delaunay三角化,形成具有新的连接关系的网格;
以图4、图5为例,对该步进行描述,图4是一个吉他的边界轮廓,图5是对边界顶点进行Delaunay三角化后的结果。在本发明中,输入是孔洞边界,然后对孔洞边界顶点进行Delaunay三角化。
(4.2)孔洞边界上的边在新的网格中可能并不存在,在该网格中插入,同时去除相交的边,对该边邻域重新三角化,使其满足Delaunay规则;
以图6为例,对该步进行描述,图6是插入吉他边界轮廓上原本存在的边,去除与这些边相交的边,并对该边邻域重新三角化后的结果。在本发明中,对邻面孔洞边界顶点进行三角化后,也需要插入孔洞边界上原本存在的边。
(4.3)去除孔洞边界外的三角形;
以图7为例,对该步进行描述,图7是去除吉他轮廓以外的三角形后的结果。在本发明中,也需要去除孔洞边界外的三角形。
(4.4)在三角形内插入新的顶点,重新进行Delaunay三角化,并使其满足最小角度和最大面积约束。
以图8为例,对该步进行描述,图8是在三角形内部插入新的顶点并重新进行Delaunay三角化后的结果。在本发明中,给定三角形最大面积参数,使内部三角形面积都不大于以边界平均长度为边的等边三角形的面积。
在该步骤的具体实施中,使用Triangle库实现,输入孔洞边界和给定参数,输出三角化后的结果。
(5)采用均值坐标原理将三角化后的孔洞反投影回原三维空间,得到修补后的网格曲面,进一步包括如下步骤:
(5.1)若顶点属于边界顶点,在将该顶点对应投影回原三维坐标顶点;
(5.2)若顶点属于边界内顶点,即新添加的顶点,则利用均值坐标计算出其三维坐标,方法如下:
结合图9,对于任意顶点v0,其坐标可以用其一环邻域顶点vi坐标表示:
其中,k为v0一环邻域内点的个数,λi为单位化后的顶点vi的权重系数,wi为单位化前的顶点vi的权重系数,αi为边v0vi和v0vi+1的夹角,αi-1为边v0vi-1和v0vi的夹角,顶点vi-1,vi,vi+1呈逆时针顺序。
(6)对于任意顶点v0,拉普拉斯算子可表示为
其中,vi为v0的一环邻域顶点,ω(v0,vi)为vi相对于v0的权重, 这里取ω(v0,vi)=1。
则高阶的拉普拉斯算子可以表示为这里二阶拉普拉斯算子中k=2;
(6.2)建立如下方程组进行求解,得到经过光顺后的顶点坐标:
其中,k=2,P=(v1,....,vp)T是孔洞内部的可自由移动的顶点,p为可自由移动的顶点个数,F=(f1,...,fn)T包含边界顶点和边界外围一环邻域内的点,这些点为固定点,n为固定点的个数,矩阵大小为p*(p+n),IF为n*n的单位矩阵。
(7)采用迪杰斯特拉算法寻找相邻牙齿对应顶点之间的最短路径,构造牙龈孔洞边界,包括步骤:
(7.1)采用迪杰斯特拉算法寻找相邻牙齿对应顶点(牙齿邻面孔洞边界构造时选取的顶点)之间的最短路径,在图10中,A,B和C,D分别为对应顶点。曲线AB,CD分别为对应顶点之间的最短路径。
(7.2)删除相邻牙齿之间,两条最短路径之上的冗余网格,两条最短路径与两邻牙的邻面边界构成了一条封闭的牙龈孔洞边界,在图11中,ABCD为牙龈孔洞边界。
(8)采用弦长参数化方法将牙龈孔洞边界投影到二维平面的单位正方形上,包括步骤:
(8.1)采用弦长参数化方法将牙龈孔洞边界直接投影到二维平面的单位正方形上。弦长参数化公式如下:
(ti+1-ti)/(ti-ti-1)=||xi+1-xi||/||xi-xi-1||
其中,ti为曲线参数,xi为第i个顶点的坐标,可表示为(x(ti),y(ti),z(ti))。
Claims (9)
1.一种面向牙齿牙龈的高效网格修补方法,其特征在于,包括步骤:
(1)选取牙齿网格模型边界上两点作为牙齿邻面孔洞边界的起止顶点,构造牙齿邻面孔洞边界;
(2)采用主元分析算法计算出牙齿邻面孔洞边界的第三主元方向;
(3)将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上,然后将平面旋转至XY平面上;
(4)采用受约束的Delaunay三角化算法对XY平面上的投影孔洞边界及其内部进行三角化;
(5)采用均值坐标原理将三角化后的孔洞反投影回原三维空间,得到修补后的网格曲面;
(6)采用二阶的拉普拉斯光顺算法优化上述修补后的网格曲面,得到光顺后的牙齿网格模型;
(7)采用迪杰斯特拉算法寻找相邻牙齿对应顶点之间的最短路径,构造牙龈孔洞边界;
(8)采用弦长参数化方法将牙龈孔洞边界投影到二维平面的单位正方形上;
(9)对投影到二维平面的单位正方形上的牙龈孔洞边界采用步骤(4)~(6)得到修补后的牙龈网格模型。
2.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(1)中,所述选取牙齿网格模型边界上两点作为牙齿邻面孔洞边界的起止顶点,构造牙齿邻面孔洞边界,包括步骤:
(1.1)用鼠标选择牙齿网格模型边界上的两点作为牙齿邻面孔洞边界的起止顶点;
(1.2)连接以上选取的起始顶点和截止顶点,离散该连线,离散长度为牙齿网格模型边界的网格距离平均值,离散后的连线与牙齿网格模型部分边界构成一个封闭的邻面孔洞边界。
3.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(2)中,所述的采用主元分析算法计算出牙齿邻面孔洞边界的第三主元方向,包括步骤:
(2.1)将边界上的顶点坐标表示成一个3×n的矩阵,其中n为边界顶点的个数,矩阵列向量为顶点的(x,y,z)坐标;
(2.2)在每个矩阵行向量上减去该行向量的平均值得到矩阵X,XT为X的转置;
(2.3)对XXT进行特征分解,求取特征向量以及所对应的特征值,最小特征值对应的特征向量即为第三主元方向。
4.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(3)中,所述的将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上,然后将平面旋转至XY平面上,包括步骤:
(3.1)将牙齿邻面孔洞边界投影到与第三主元方向相垂直的平面上;
(3.2)计算以上平面与XY平面之间的旋转矩阵,将孔洞边界顶点的投影坐标乘以旋转矩阵即得到顶点在XY平面上的坐标。
5.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(4)中,所述的采用受约束的Delaunay三角化算法对XY平面上的投影孔洞边界及其内部进行三角化,包括步骤:
(4.1)对XY平面上的投影孔洞边界顶点进行Delaunay三角化,形成具有新的连接关系的网格;
(4.2)孔洞边界上的边在新的网格中可能并不存在,在该网格中插入,同时去除相交的边,对该边邻域重新三角化;
(4.3)去除孔洞边界外的三角形;
(4.4)在三角形内插入新的顶点,重新进行Delaunay三角化,并使其满足最小角度和最大面积约束。
6.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(5)中,所述采用均值坐标原理将三角化后的孔洞反投影回原三维空间,得到修补后的网格曲面,包括步骤:
(5.1)若顶点属于边界顶点,则将该顶点对应投影回原三维坐标顶点;
(5.2)若顶点属于边界内顶点,则利用均值坐标计算出其三维坐标,方法如下:
对于任意顶点v0,其坐标可以用其一环邻域顶点vi坐标表示:
其中,k为v0一环邻域内点的个数,λi为单位化后的顶点vi的权重系数,wi为单位化前的顶点vi的权重系数,αi为边v0vi和v0vi+1的夹角,αi-1为边v0vi-1和v0vi的夹角,顶点vi-1,vi,vi+1呈逆时针顺序。
7.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(6)中,所述采用二阶的拉普拉斯光顺算法优化上述修补后的网格曲面,得到光顺后的牙齿网格模型,包括步骤:
(6.1)对于任意顶点v0,拉普拉斯算子可表示为
其中,vi为v0的一环邻域顶点,ω(v0,vi)为vi相对于v0的权重,
则高阶的拉普拉斯算子可以表示为这里二阶拉普拉斯算子中k=2;
(6.2)建立如下方程组进行求解,得到经过光顺后的顶点坐标:
其中,P=(v1,....,vp)是孔洞内部的可自由移动的顶点,F=(f1,...,fF)包含边界顶点和边界外围一环邻域内的点,这些点为固定点,IF为单位矩阵。
8.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(7)中,所述采用迪杰斯特拉算法寻找相邻牙齿对应顶点之间的最短路径,构造牙龈孔洞边界,包括步骤:
(7.1)采用迪杰斯特拉算法寻找相邻牙齿对应顶点之间的最短路径;
(7.2)删除相邻牙齿之间,两条最短路径之上的冗余网格,两条最短路径与两邻牙的邻面边界构成了一条封闭的牙龈孔洞边界。
9.如权利要求1所述的面向牙齿牙龈的高效网格修补方法,其特征在于,步骤(8)中,所述的采用弦长参数化方法将牙龈孔洞边界投影到二维平面的单位正方形上,包括步骤:
(8.1)采用弦长参数化方法将牙龈孔洞边界直接投影到二维平面的单位正方形上,弦长参数化公式如下:
(ti+1-ti)/(ti-ti-1)=||xi+1-xi||/||xi-xi-1||
其中,ti为曲线参数,xi为第i个顶点的坐标,可表示为(x(ti),y(ti),z(ti))。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409293.XA CN102930091B (zh) | 2012-10-24 | 2012-10-24 | 一种面向牙齿牙龈的高效网格修补方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210409293.XA CN102930091B (zh) | 2012-10-24 | 2012-10-24 | 一种面向牙齿牙龈的高效网格修补方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102930091A CN102930091A (zh) | 2013-02-13 |
CN102930091B true CN102930091B (zh) | 2015-05-20 |
Family
ID=47644888
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210409293.XA Active CN102930091B (zh) | 2012-10-24 | 2012-10-24 | 一种面向牙齿牙龈的高效网格修补方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102930091B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104504663B (zh) * | 2014-12-29 | 2018-05-01 | 佛山市诺威科技有限公司 | 一种义齿三角网格模型孔洞的迭代修补方法 |
CN104504759B (zh) * | 2014-12-29 | 2018-10-02 | 佛山市诺威科技有限公司 | 一种基于义齿基底冠三角网格的快速过渡缝补的方法 |
CN104599327B (zh) * | 2015-01-30 | 2017-12-26 | 杭州一牙数字口腔有限公司 | 一种牙颌黑三角的自动检测和修复方法 |
CN105551081B (zh) * | 2016-02-05 | 2018-04-10 | 杭州美齐科技有限公司 | 一种虚拟牙龈三角网格构建和随动方法 |
CN106204737B (zh) * | 2016-07-14 | 2019-06-18 | 广东工业大学 | 一种牙齿的咬合表面恢复方法及系统 |
CN106204730B (zh) * | 2016-07-19 | 2019-04-26 | 上海爱圣美科技有限公司 | 一种三维牙齿模型中牙龈随动的方法 |
CN106504318B (zh) * | 2016-10-19 | 2019-08-23 | 青岛达芬奇科技有限公司 | 一种用于无托槽正畸的牙龈生成方法 |
CN108567501B (zh) * | 2017-03-14 | 2021-04-27 | 无锡时代天使医疗器械科技有限公司 | 牙齿邻接面的重建方法 |
CN107067476A (zh) * | 2017-05-11 | 2017-08-18 | 重庆市劢齐医疗科技有限责任公司 | 一种基于高阶拉普拉斯算子的快速网格补洞方法 |
CN109493380B (zh) * | 2017-09-11 | 2022-05-17 | 重庆大学 | 一种岩石节理面抗剪试验中不规则剪切面面积的计算方法 |
CN108992188B (zh) * | 2018-07-12 | 2020-10-20 | 杭州美齐科技有限公司 | 一种单颗牙齿网格模型的修补算法 |
CN109872394B (zh) * | 2019-01-10 | 2021-04-27 | 重庆大学 | 基于最小二乘支持向量机的狭长三角形网格优化方法 |
US10984592B2 (en) | 2019-07-31 | 2021-04-20 | Siemens Industry Software Inc. | Spatial decomposition-based repair of defects in boundary representations |
JP7339846B2 (ja) * | 2019-10-25 | 2023-09-06 | 住友重機械工業株式会社 | モデル化装置及びプログラム |
CN112206064B (zh) * | 2020-09-25 | 2021-12-28 | 雅客智慧(北京)科技有限公司 | 一种数字化牙模生成方法及系统 |
CN112381950B (zh) * | 2020-11-05 | 2024-07-23 | 新拓三维技术(深圳)有限公司 | 一种网格孔洞修补方法、电子设备和计算机可读存储介质 |
CN112998888B (zh) * | 2021-02-01 | 2022-04-22 | 重庆邮电大学 | 一种基于网格投影的义齿模型倒凹去除方法 |
CN117094197B (zh) * | 2023-08-24 | 2024-05-14 | 上海慕灿信息科技有限公司 | 一种基于cdt参数化网格与全局各向同性重新网格化的方法 |
CN118071942B (zh) * | 2024-04-24 | 2024-08-06 | 福州海狸家口腔科技有限公司 | 一种牙冠三维模型的空洞补全方法、介质和设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708591A (zh) * | 2012-04-28 | 2012-10-03 | 浙江大学 | 一种数字牙网格模型的孔洞修补方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060127852A1 (en) * | 2004-12-14 | 2006-06-15 | Huafeng Wen | Image based orthodontic treatment viewing system |
-
2012
- 2012-10-24 CN CN201210409293.XA patent/CN102930091B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708591A (zh) * | 2012-04-28 | 2012-10-03 | 浙江大学 | 一种数字牙网格模型的孔洞修补方法 |
Non-Patent Citations (2)
Title |
---|
三角网格模型的补洞算法研究;田建磊等;《计算机应用》;20090801;第29卷(第8期);第2035-2037,2059页 * |
牙模型编辑系统的研究与实现;杜中义;《中国优秀硕士学位论文全文数据库信息科技辑》;20120715;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN102930091A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102930091B (zh) | 一种面向牙齿牙龈的高效网格修补方法 | |
KR101741402B1 (ko) | 평행 지오데식 곡선들에 의해 모델링되는 부품의 설계 | |
Alliez et al. | Isotropic surface remeshing | |
JP6242497B2 (ja) | ポケット形状を機械加工する工具の経路を決定する方法及びシステム | |
Campen et al. | Quad layout embedding via aligned parameterization | |
CN102708591B (zh) | 一种数字牙网格模型的孔洞修补方法 | |
JP4934789B2 (ja) | 補間処理方法および補間処理装置 | |
Garth et al. | Generation of accurate integral surfaces in time-dependent vector fields | |
KR100707841B1 (ko) | 3차원 목표곡선을 이용한 곡면 변형 장치 및 그 방법 | |
JP2013507679A (ja) | 三次元物体モデルの3dプリントが可能な方法及びシステム | |
Louhichi et al. | CAD/CAE integration: updating the CAD model after a FEM analysis | |
JP4662959B2 (ja) | デザイン最適化のための自由形状変形(FreeFormDeformation)表現の進化的直接操作(EvolutionaryDirectManipulation) | |
US20080225043A1 (en) | Computer-Implemented Process and System for Creating a Parametric Surface | |
US9269189B1 (en) | Conversion of T-spline models to trim-free T-spline models: functional case | |
EP2930691A1 (en) | Fitting sample points with an isovalue surface | |
JP2002501639A (ja) | プログレッシブメッシュの適応細分方法および装置 | |
JP2010277585A (ja) | サンプル点において掃引容積の距離場を再構成する方法 | |
JP2008305372A (ja) | 3次元スキャンデータを用いてロフトサーフィスを計算するためのシステム及び方法 | |
Bronson et al. | Lattice cleaving: A multimaterial tetrahedral meshing algorithm with guarantees | |
EP2930692A1 (en) | Fitting sample points of 3D curves sketched by a user with an isovalue surface | |
US7062353B2 (en) | Method and system for reconstructing a surface | |
JP2010033568A (ja) | 放電加工工程用の電極形状を設計するための装置および方法 | |
Li et al. | On surface reconstruction: A priority driven approach | |
Lee et al. | Fast and intuitive metamorphosis of 3d polyhedral models using smcc mesh merging scheme | |
Zhao et al. | Geometry of axisymmetric 3D origami consisting of triangular facets |
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 |