一种用于机器人辅助牙体预备的邻面备牙曲线生成方法技术

一种用于机器人辅助牙体预备的邻面备牙曲线生成方法，它涉及机器人辅助牙体预备技术领域，本发明专利技术将邻面备牙阶段规划为机器人的三维曲线轨迹规划，基于医生的标准化后牙手工预备体模型，获取邻面备牙的曲线数据，根据提取的邻面备牙的曲线数据完成了邻面曲线的反算与插补，进而实现了机器人的三维曲线轨迹规划，建立一种用于机器人辅助牙体预备的邻面备牙曲线生成方法。技术要点为：预备体模型预处理；提取标准化邻面备牙曲线及离散点；计算节点矢量U和控制点P；曲线插补。本发明专利技术根据提取的邻面备牙的曲线数据基于NURBS曲线完成了邻面曲线的反算与插补，进而实现了机器人辅助牙体预备的三维曲线轨迹规划，为实现机器人辅助牙体预备奠定了基础。

【技术实现步骤摘要】
一种用于机器人辅助牙体预备的邻面备牙曲线生成方法
本专利技术专利涉及一种用于机器人辅助牙体预备的邻面备牙曲线生成方法，属于机器人辅助牙体预备

技术介绍
龋齿是导致牙体缺损的重要原因，严重的影响了人们的口腔健康。口腔修复是治疗牙体缺损的重要手段，牙体预备是修复过程中的必要治疗环节，是指医生对患龋齿处的牙齿进行定量去除并形成预期三维形状的操作过程。在传统的牙体预备过程中，需要依赖于医生手工的操作并结合丰富的临床经验对患龋牙齿进行大量的重复的微细调整。然而目前我国医患比失衡的现状无法满足当前极大的牙体预备需求量，因此需要引入机器人和辅助软件来代替或辅助医生高效的完成牙体预备工作，不仅能够缓解医患失衡的难题，而且还能有效的提高备牙质量和口腔治疗的效果。在机器辅助完成牙体预备的过程中，机器人末端的运动轨迹决定了牙体预备的精度，进而影响治疗的效果和成功率。因此，获取合理的邻面备牙轨迹，着重考虑预备体与后牙修复体的贴合程度，从而减少因牙体预备不当所形成口腔继发疾病的情况，是目前机器人辅助牙体预备领域的研究热点。
技术实现思路
针对上述问题，本专利技术提出一种用于机器人辅助牙体预备的邻面备牙曲线生成方法，解决目前机器人辅助牙体预备
缺少获得邻面备牙曲线的方法，以提高牙体预备的精度和成功率，进而实现机器人辅助牙体预备。本专利技术为解决上述问题所采取的方案为：一种用于机器人辅助牙体预备的邻面备牙曲线生成方法，所述方法的具体实现过程为：步骤一、预备体模型预处理根据医生临床的备牙模型进行扫描，得到obj格式的标准化预备体模型，选用GeomagicWrap逆向工程软件对obj格式的标准化备牙三维模型进行处理，使用“网格医生”命令对模型中的非流体边、自相交、高度折射边或钉状物等三角形片问题进行修复，使用“松弛多边形”命令通过增加三角形的数量使其表面更加光滑方便进一步的曲面化，使用“自动曲面化”设置“几何图形类型”为“机械”，设置“曲面片计数”为“自动评估”，“曲面细节”适中，“曲面拟合”的“适应性公差”为0.0097，将生成的曲面片通过“构造栅格”命令，在每个曲面片内构建细分网格，完成从三角形片组成的obj格式三维模型输出具有曲面特征和边界特征的step格式的三维模型；步骤二、提取标准化邻面备牙曲线及离散点将获得的step格式的预备体模型导入到Creo三维软件中，设置车针直径D的2/3部分在即将预备的肩台内部，剩余的D/3位于牙体外部，车针轴线距牙体边缘曲线CL的距离为D/6，复制标准化预备体的牙体边缘曲线CL作为参考曲线，再使用“偏移”命令输入偏移量D/6对CL完成偏移得到邻面备牙曲线CG，在邻面备牙曲线CG上按离散点数量NQ和离散点间距LQ两个参数获得多个离散点Q；步骤三、邻面备牙曲线的反算将预处理阶段中的离散点Q通过基于NURBS的反求计算得到邻面备牙曲线，反算目标为计算节点矢量U和控制顶点P，引入NURBS曲线数学模型，ωi(i＝0,1,...,n)为权重因子；Pi(i＝0,1,...,n)为控制顶点，其数量为n+1；p为NURBS曲线次数；Ni,p(u)为基函数；公式(1)中的Ni,p(u)为：式中：ui是非均匀节点向量U的元素，如公式(3)所示；m+1为节点向量U的长度；m、p、n三者的关系为m＝n+p+1；a和b通常为0和1；由于邻面备牙曲线是一条封闭的三维曲线，因此首末数据点重合Q0＝Qm，该邻面备牙曲线为4自由度3次NURBS曲线，首末3个控制点依次相等，得到等式：P0＝Pn-2，P1＝Pn-1，P2＝Pn，将NURBS曲线的有理分式用矩阵的形式表示：式中：t∈[0,1],令特殊的Δi＝ui+1-ui,m12＝1-m11-m13，根据确定的首末两端条件，得到封闭曲线情况下的控制点P的计算公式：式中各个参数如公式(7)所示：式中：i＝1,2,...,n-2；将预处理阶段得到的离散点Q通过均匀参数化法、积累弦长参数化法和向心参数化法分别计算邻面曲线的节点矢量U，积累弦长参数化法所得曲线反应了离散点Q按照弦长的分布情况，备牙曲线更加接近预备体模型所提取的标准化邻面备牙曲线CG，提高了备牙曲线的预备精度。步骤四、邻面备牙曲线插补邻面备牙曲线插补原理是利用时间序列{t1,t2,...,tk,...,tn-1,tn}分割参数序列{u1,u2,...,uk,uk+1,...,un-1,un}，进而得到插补点序列{C(u1),C(u2),...,C(uk),...,C(un-1),C(un)}，插补的核心计算就是利用插补周期T求得uk和uk+1之间的关系，进一步通过C(uk)得到C(uk+1)；首先，对NURBS曲线的导数进行推导，得到基于NURBS的邻面备牙曲线对u的一阶导数C(1)(u)；式中：C(1)(u)为曲线对u的一阶导数；继续对曲线求二阶导数得到C(2)(u)：式中：C(2)(u)为曲线对u的二阶导数，基函数Ni,p(m)(u)的表达式为：将NURBS曲线上的进给速率表示为：整理变换公式(11)得到公式(12)：式中：Cx(1)(uk)为曲线x方向的一阶导；Cy(1)(uk)为曲线y方向的一阶导；Cz(1)(uk)为曲线z方向的一阶导。继续对公式(12)求二阶导数，得到：式中：Cx(2)(uk)为曲线x方向的二阶导；Cy(2)(uk)为曲线y方向的二阶导；Cz(2)(uk)为曲线z方向的二阶导；整理得到了基于NURBS的邻面备牙曲线的几何特性和运动特性之间的关系:本专利技术的有益效果为：1、本专利技术提出了一种用于机器人辅助牙体预备的邻面备牙曲线生成方法，在标准化预备体模型的预处理阶段对模型进行离散化曲面造型，对原有的三角形片进行连接和重建，进一步进行曲面拟合，使处理后的预备体模型具备曲面特征和边界特征，提高所提取备牙曲线的精度和拟合度。2、本专利技术在邻面备牙曲线的反算过程中，选取NURBS来表达邻面备牙曲线，利用节点间距不同，让控制点对于邻面备牙曲线获得更加灵活的控制范围，邻面备牙曲线为4自由度3次NURBS曲线，使首末3个控制点依次相等，提高了邻面备牙曲线的运算速率，并使得插补后的曲线的拐点数量。3、本专利技术在邻面备牙曲线的反算过程中，考虑到预备过程中对目标牙齿邻牙的无效接触，利用累弦长参数化法所得曲线的特性，反应了离散点Q按照弦长的分布情况，有效提高了备牙曲线的预备精度，邻面备牙曲线保持一致且光滑连续。4、本专利技术以标准化后牙预备体模型作为参照，将邻面备牙阶段规划为机器人的三维曲线轨迹规划，通过对标准后牙全冠预备体模型的预处理，提取了邻面备牙的曲线数据，根据提取的邻面备牙的曲线数据基于NURBS曲线完成了邻面曲线的反算与插补，进而实现了机器人辅助牙体预备本文档来自技高网...

【技术保护点】
1.一种用于机器人辅助牙体预备的邻面备牙曲线生成方法，其特征在于：所述方法的具体实现过程为：/n步骤一、预备体模型预处理/n根据医生临床的备牙模型进行扫描，得到obj格式的标准化预备体模型，选用GeomagicWrap逆向工程软件对obj格式的标准化备牙三维模型进行处理，使用“网格医生”命令对模型中的非流体边、自相交、高度折射边或钉状物等三角形片问题进行修复，使用“松弛多边形”命令通过增加三角形的数量使其表面更加光滑方便进一步的曲面化，使用“自动曲面化”设置“几何图形类型”为“机械”，设置“曲面片计数”为“自动评估”，“曲面细节”适中，“曲面拟合”的“适应性公差”为0.0097，将生成的曲面片通过“构造栅格”命令，在每个曲面片内构建细分网格，完成从三角形片组成的obj格式三维模型输出具有曲面特征和边界特征的step格式的三维模型；/n步骤二、提取标准化邻面备牙曲线及离散点/n将获得的step格式的预备体模型导入到Creo三维软件中，设置车针直径D的2/3部分在即将预备的肩台内部，剩余的D/3位于牙体外部，车针轴线距牙体边缘曲线C

【技术特征摘要】
1.一种用于机器人辅助牙体预备的邻面备牙曲线生成方法，其特征在于：所述方法的具体实现过程为：
步骤一、预备体模型预处理
根据医生临床的备牙模型进行扫描，得到obj格式的标准化预备体模型，选用GeomagicWrap逆向工程软件对obj格式的标准化备牙三维模型进行处理，使用“网格医生”命令对模型中的非流体边、自相交、高度折射边或钉状物等三角形片问题进行修复，使用“松弛多边形”命令通过增加三角形的数量使其表面更加光滑方便进一步的曲面化，使用“自动曲面化”设置“几何图形类型”为“机械”，设置“曲面片计数”为“自动评估”，“曲面细节”适中，“曲面拟合”的“适应性公差”为0.0097，将生成的曲面片通过“构造栅格”命令，在每个曲面片内构建细分网格，完成从三角形片组成的obj格式三维模型输出具有曲面特征和边界特征的step格式的三维模型；
步骤二、提取标准化邻面备牙曲线及离散点
将获得的step格式的预备体模型导入到Creo三维软件中，设置车针直径D的2/3部分在即将预备的肩台内部，剩余的D/3位于牙体外部，车针轴线距牙体边缘曲线CL的距离为D/6，复制标准化预备体的牙体边缘曲线CL作为参考曲线，再使用“偏移”命令输入偏移量D/6对CL完成偏移得到邻面备牙曲线CG，在邻面备牙曲线CG上按离散点数量NQ和离散点间距LQ两个参数获得多个离散点Q；
步骤三、计算节点矢量U和控制点P
将预处理阶段中的离散点Q通过基于NURBS的反求计算得到邻面备牙曲线，反算目标为计算节点矢量U和控制顶点P，引入NURBS曲线数学模型，



ωi(i＝0,1,...,n)为权重因子；Pi(i＝0,1,...,n)为控制顶点，其数量为n+1；p为NURBS曲线次数；Ni,p(u)为基函数；公式(1)中的Ni,p(u)为：



式中：ui是非均匀节点向量U的元素，如公式(3)所示；



m+1为节点向量U的长度；m、p、n三者的关系为m＝n+p+1；a和b通常为0和1；
由于邻面备牙曲线是一条封闭的三维曲线，因此首末数据点重合Q0＝Qm，该邻面备牙曲线为4自由度3次NURBS曲线，首末3个控制点依次相等，得到等式：P0＝Pn-2，P1＝Pn-1，P2＝Pn，将NURBS曲线的有理分式用矩阵的形式表示：





...

【专利技术属性】
技术研发人员：姜金刚，钱伟，黄致远，张永德，唐德栋，曾阳，
申请(专利权)人：哈尔滨理工大学，
类型：发明
国别省市：黑龙江;23