一种基于弯制点角距比和的正畸弓丝变角度划分方法与流程

本发明专利涉及一种基于弯制点角距比和的正畸弓丝变角度划分方法，属于正畸弓丝弯制技术领域。

背景技术：

错颌畸形是危机人体健康的第三大口腔疾病，在人群中呈现较高的发病率，现代口腔医学中，固定矫治是一种常用且有效的正畸治疗手段，而正畸弓丝的弯制是固定矫治技术的关键，近年来，深受数字化制造技术的影响，传统的口腔制造加工工艺正发生革命性变化，口腔正畸领域也受益于数字化技术，正畸矫治器中弓丝的加工正在向数字化发展，然而弯制点顺序的高效规划是实现弓丝数字化加工，提高机器人弯制正畸弓丝效率以及保证弯制精度的前提。

在个性化弓丝弯制点顺序规划时，个性化正畸弓丝上存在弯制点分布相对分散，即正畸弓丝曲线上各弯制点的单位弯制点密度都很小的情况，当弓丝曲线上各弯制点的单位弯制点密度值均小于根据机器人弯制正畸弓丝而设定的单位弯制点密度上限值时，由于其弯制点分布的特殊性，在对这类个性化正畸弓丝进行变角度域划分时，现有的正畸弓丝成形控制点弯制顺序规划计算量相对较大，规划期间耗时相对较大，目前正畸弓丝弯制技术领域缺少对此类个性化正畸弓丝弯制点顺序规划的方法，无法对此类个性化正畸弓丝实现高效数字化弯制。

技术实现要素：

针对上述问题，本发明提出一种基于弯制点角距比和的正畸弓丝变角度划分方法，解决目前正畸弓丝弯制技术领域缺少针对此类每个弯制点单位弯制点密度都相对较小的正畸弓丝的高效弯制顺序规划方法，以避免对具有特殊属性的此类正畸弓丝的弯制点进行划分时出现划分区域内弯制难度较大的情况，有利于弯制顺序的合理规划，避免弯制过程中出现干涉的情况，进而实现正畸弓丝高效的数字化弯制。

本发明为解决上述问题所采取的方案为：一种基于弯制点角距比和的正畸弓丝变角度划分方法，所述方法的具体实现过程为：

步骤一、变角度域划分数据导入：

根据患者有i个成形控制点的个性化正畸弓丝曲线，输入个性化正畸弓丝曲线成形控制点信息集m＝{m1，m2，m3，...，mi}，mi＝(ui，vi，wi)'，为每个个性化正畸弓丝曲线成形控制点的坐标，每一个正畸弓丝曲线成形控制点mi均对应一个成形控制点机器人运动信息单元ni，输入的成形控制点的机器人运动信息集为n＝{n1，n2，n3，...，ni}，ni表示机器人在弯制该点时的成形控制点坐标及弯制角度，ni＝(ui，vi，wi，αi)'，ui、vi、wi为该成形控制点mi的坐标，αi为机器人作用在成形控制点mi上的弯制角度，将正畸弓丝曲线成形控制点信息集m、成形控制点的机器人运动信息集n，单位弯制点密度上限值ρmax输入到正畸弓丝弯制系统中，并且根据机器人弯制正畸弓丝的特点设定变角度域弯制点角距比和阀值(∑e)max；

个性化正畸弓丝曲线两端点坐标为ms(us，vs，ws)，me(ue，ve，we)，将弓丝曲线左端点ms设定为成形控制点m0，将弓丝曲线右端点me设定为成形控制点mi+1，以正畸弓丝曲线左右端点的中点定义为圆心o，将成形控制点mk作为划分变角度域的起点，k的初始值为k＝0，0≤k≤i，成形控制点mk与圆心o之间的连线定义为划分半径rk，表示划分变角度域bk内划分的弯制点个数，的初始值为

步骤二、正畸弓丝曲线坐标转换及划分变量定义：

将个性化正畸弓丝成形控制点信息集中各成形控制点的坐标mi＝(ui,vi,wi)'中的wi赋值为0，即令wi＝0，获得正畸弓丝转换曲线m′；

针对变角度划分方法而言，定义弯制点mk的角距比为弯制点mk的弯制角度αk与其弯制距离的比值，即其中表示弯制点mk与弯制点mk+1之间的弧长，表示变角度域bk内所有弯制点角距比的累加和，个性化正畸弓丝上任意一弯制点mx的单位弯制点密度其中βx表示弯制点mx的划分半径rx与弯制点mx+1的划分半径rx+1之间形成的角度，1≤x≤i，计算个性化正畸弓丝上各弯制点的单位弯制点密度当各弯制点的单位弯制点密度均小于单位弯制点密度上限值ρmax时，跳转至步骤三；

步骤三、计算变角度域弯制点角距比和：

a)判断是否成立；

如果成立，说明已划分至弓丝曲线的右端点，则用公式计算变角度域bk的弯制点角距比累加和，将弯制点mk和弯制点之间的变角度域定义为变角度域区间bk，变角度域区间bk的角距比累加和为跳转至步骤四，统计已划分的弯制点数目q；

如果不成立，说明未划分至弓丝曲线的右端点，则用公式计算变角度域bk的弯制点角距比累加和，并且判断是否成立；

如果成立，则可以继续增加该变角度域内弯制点的个数，即令跳转至步骤三a；

b)如果不成立，说明该变角度域内不能再增加弯制点，跳转至步骤四；

步骤四、定义变角度域区间：

将弯制点mk和弯制点之间的变角度域定义为变角度域区间bk，并且变角度域区间bk内从弯制点mk到弯制点包括弯制点但是不包括弯制点mk的所有弯制点已划分，统计已划分的弯制点数目q，变角度域区间bk的弯制点角距比累加为跳转至步骤五；

步骤五、判断是否继续进行变角度域划分：

判断q＝i是否成立；

如果q＝i不成立，则继续进行变角度域划分，由于当前弯制点不符合弯制要求，故从当前弯制点的上一个弯制点继续开始变角度划分，即以划分变角度域的起点，令跳转至步骤一；

如果q＝i成立，说明所有的弯制点均已被划分，跳转至步骤六；

步骤六、获得最终弯制点弯制顺序：

输出变角度域区间信息集b＝{b1，b2，...，bg}，每个变角度域区间(b1，b2，...，bg)所对应的角度值分别为(θ1，θ2，...，θg)，并且θ1+θ2+...+θg＝π，比较每个变角度域区间(b1，b2，...，bg)的弯制点角距比和得到则以弯制点角距比和为指标将g个变角度域区间降序排列，得到降序变角度域区间信息集c＝{c2，c4，...，cg}，在任何一个变角度域区间内，以各弯制点的单位弯制点密度ρmk为指标对弯制点进行降序排列，将单位弯制点密度降序排列顺序定义为该变角度域区间内的弯制点弯制顺序，进而得到个性化正畸弓丝曲线成形控制点坐标降序矩阵m1{m4，m6，...，ms}和机器人运动降序信息集n1{n4，n6，...，ns}，其中ms表示变角度域区间bg中的弯制点，输出最终弯制点弯制顺序m1{m2，m4，...，mg}，n1{n2，n4，...，ng}，程序结束。

本发明的有益效果为：

1.本发明在对正畸弓丝进行变角度划分时，提出了单位弯制点密度的概念，从而对弯制点在正畸弓丝曲线上的分布情况进行量化描述，并且根据机器人弯制正畸弓丝的特点设定了单位弯制点密度上限值ρmax，从而限定了此时正畸弓丝曲线上弯制点分布的集中程度，方便下一步利用弯制点角距比和对正畸弓丝曲线进行变角度域划分。

2.同本发明人同日申报的发明专利《一种用于正畸弓丝弯制规划的平面变角度划分方法》相比，本专利先对正畸弓丝曲线上弯制点分布的集中程度进行判断，即判断正畸弓丝曲线上各弯制点的单位弯制点密度均满足划分要求后，再对其进行变角度划分，此时划分时只需计算变角度域弯制点角距比累加和，从而减少了循环次数，降低计算规模，提高了划分效率。

3.本发明在根据变角度域弯制点角距比累加和对各变角度区间进行降序排列之后，在各个变角度区间内部，使用每个弯制点的单位弯制点密度的降序排列作为该区间内弯制点的排列顺序，最后可以输出合理的弯制顺序，保证了正畸弓丝弯制规划的完整性。

4.同本发明人已经授权的授权公告号为cn107714203b的发明专利《一种等角度划分的正畸弓丝弯制顺序规划方法》相比，本方法是在平面变角度划分的基础上，针对具有特定属性的一类弓丝，充分考虑此类正畸弓丝曲线上弯制点分布信息个性化特点，即患者个性化正畸弓丝上弯制点的单位弯制点密度都相对较小，各弯制点的单位弯制点密度均小于规定的上限值，提出了基于变角度域弯制点角距比和，用于划分此类正畸弓丝曲线的划分方法，使得划分过程并非以一种无依据的均匀化标准进行划分，而是为了使划分的变角度域符合弯制复杂程度的上限要求，不断改变划分变角度域的角度值，生成了一系列满足变角度域弯制点角距比和要求的变角度域区间，提高了正畸弓丝成形控制点弯制顺序规划方法的合理性，有效避免了弯制机器人产生空程无效动作、弯制过程相互干涉动作和弯制运动复杂动作，充分发挥弯制机器人优势的最大化，明显提高了弯制效率。

5、同本发明人同日申报的发明专利《一种基于弯制点密度的正畸弓丝变角度划分方法》相比，虽然两方法都适用于具有特殊属性的一类个性化正畸弓丝曲线，但对于《一种基于弯制点密度的正畸弓丝变角度划分方法》中提及的方法所侧重的前提是每个弯制点的单位角距比都满足设定要求，进而仅以变角度域弯制点密度作为划分变角度域的依据，本方法所侧重的前提是弯制点的单位弯制点密度满足设定要求，进而仅以变角度域弯制点角距比和作为划分变角度域的依据，两种方法在进行正畸弓丝弯制顺序规划时应用情况不同，因此本方法的提出与另种方法互相补偿，进而完善了正畸弓丝弯制顺序规划的系列方法。

附图说明

为了易于说明，本发明由下述的具体实施及附图作以详细描述。

图1为一种基于弯制点角距比和的正畸弓丝变角度划分方法流程图；

图2为单位弯制点密度计算示意图；

图3为基于弯制点角距比和变角度划分个性化正畸弓丝曲线示意图；

图4为个性化正畸弓丝成形控制点分布示意图；

图5为变角度划分个性化正畸弓丝曲线完成后规划弯制点弯制顺序示意图；

具体实施方式

为使本发明专利的目的、技术方案和优点更加清楚明了，下面通过附图中示出的具体实施例来描述本发明专利，但是应该理解，这些描述只是示例性的，而并非要限制本发明专利的范围，此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明专利的概念。

实施例1：如图1、图2、图3、图4、图5所示，本具体实施方式采用以下技术方案：一种基于弯制点角距比和的正畸弓丝变角度划分方法，所述方法的具体实现过程为：

步骤一、变角度域划分数据导入：

根据患者有i个成形控制点的个性化正畸弓丝曲线，输入个性化正畸弓丝曲线成形控制点信息集m＝{m1，m2，m3，...，mi}，mi＝(ui，vi，wi)'，为每个个性化正畸弓丝曲线成形控制点的坐标，每一个正畸弓丝曲线成形控制点mi均对应一个成形控制点机器人运动信息单元ni，输入的成形控制点的机器人运动信息集为n＝{n1，n2，n3，...，ni}，ni表示机器人在弯制该点时的成形控制点坐标及弯制角度，ni＝(ui，vi，wi，αi)'，ui、vi、wi为该成形控制点mi的坐标，αi为机器人作用在成形控制点mi上的弯制角度，将正畸弓丝曲线成形控制点信息集m、成形控制点的机器人运动信息集n，单位弯制点密度上限值ρmax输入到正畸弓丝弯制系统中，并且根据机器人弯制正畸弓丝的特点设定变角度域弯制点角距比和阀值(∑e)max；

个性化正畸弓丝曲线两端点坐标为ms(us，vs，ws)，me(ue，ve，we)，将弓丝曲线左端点ms设定为成形控制点m0，将弓丝曲线右端点me设定为成形控制点mi+1，以弓丝曲线左右端点的中点定义为圆心o，将成形控制点mk作为划分变角度域的起点，k的初始值为k＝0，0≤k≤i，成形控制点mk与圆心o之间的连线定义为划分半径rk，表示划分变角度域bk内划分的弯制点个数，的初始值为

步骤二、正畸弓丝曲线坐标转换及划分变量定义：

将个性化正畸弓丝成形控制点信息集中各成形控制点的坐标mi＝(ui,vi,wi)'中的wi赋值为0，即令wi＝0，获得正畸弓丝转换曲线m′；

针对变角度划分方法而言，定义弯制点mk的角距比为弯制点mk的弯制角度αk与其弯制距离的比值，即其中表示弯制点mk与弯制点mk+1之间的弧长，表示变角度域bk内所有弯制点角距比的累加和，弯个性化正畸弓丝上任意一弯制点mx的单位弯制点密度其中βx表示弯制点mx的划分半径rx与弯制点mx+1的划分半径rx+1之间形成的角度，1≤x≤i，计算个性化正畸弓丝上各弯制点的单位弯制点密度当各弯制点的单位弯制点密度均小于单位弯制点密度上限值ρmax时，跳转至步骤三；

步骤三、计算变角度域弯制点角距比和：

a)判断是否成立；

如果成立，说明已划分至弓丝曲线的右端点，则用公式计算变角度域bk的弯制点角距比累加和，将弯制点mk和弯制点之间的变角度域定义为变角度域区间bk，变角度域区间bk的角距比累加和为跳转至步骤四，统计已划分的弯制点数目q；

如果不成立，说明未划分至弓丝曲线的右端点，则用公式计算变角度域bk的弯制点角距比累加和，并且判断是否成立；

如果成立，则可以继续增加该变角度域内弯制点的个数，即令跳转至步骤三a；

b)如果不成立，说明该变角度域内不能再增加弯制点，跳转至步骤四；

步骤四、定义变角度域区间：

将弯制点mk和弯制点之间的变角度域定义为变角度域区间bk，并且变角度域区间bk内从弯制点mk到弯制点包括弯制点但是不包括弯制点mk的所有弯制点已划分，统计已划分的弯制点数目q，变角度域区间bk的弯制点角距比累加为跳转至步骤五；

步骤五、判断是否继续进行变角度域划分：

判断q＝i是否成立；

如果q＝i不成立，则继续进行变角度域划分，由于当前弯制点不符合弯制要求，故从当前弯制点的上一个弯制点继续开始变角度划分，即以为划分变角度域的起点，令跳转至步骤一；

如果q＝i成立，说明所有的弯制点均已被划分，跳转至步骤六；

步骤六、获得最终弯制点弯制顺序：

输出变角度域区间信息集b＝{b1，b2，...，bg}，每个变角度域区间(b1，b2，...，bg)所对应的角度值分别为(θ1，θ2，...，θg)，并且θ1+θ2+...+θg＝π，比较每个变角度域区间(b1，b2，...，bg)的弯制点角距比和得到则以弯制点角距比和为指标将g个变角度域区间降序排列，得到降序变角度域区间信息集c＝{c2，c4，...，cg}，在任何一个变角度域区间内，以各弯制点的单位弯制点密度为指标对弯制点进行降序排列，将单位弯制点密度降序排列顺序定义为该变角度域区间内的弯制点弯制顺序，进而得到个性化正畸弓丝曲线成形控制点坐标降序矩阵m1{m4，m6，...，ms}和机器人运动降序信息集n1{n4，n6，...，ns}，其中ms表示变角度域区间bg中的弯制点，输出最终弯制点弯制顺序m1{m2，m4，...，mg}，n1{n2，n4，...，ng}，程序结束。

实施例2：如图2、图3、图4、图5所示，在一条包含i＝22个弯制点(正畸弓丝的左右端点不参与弓丝弯制)的个性化正畸弓丝曲线进行基于弯制点角距比和变角度划分的弯制顺序规划过程中，假设最终得到的弯制变角度域区间的个数为n＝8，根据在步骤六中变角度域区间信息集b＝{b1，b2，...，bg}，得到各变角度域区间的弯制点个数为每个弯制变角度域区间(b1，b2，...，bg)的弯制点角距比和为比较各变角度域区间的弯制点角距比和，存在以各变角度域区间的弯制点角距比和为指标将8个变角度域区间降序排列，从而得到降序变角度域区间信息集为c1＝{c7，c5，c3，c4，c1，c8，c2，c6}，规定在任何一个变角度域区间内，以各弯制点的单位弯制点密度为指标对弯制点进行降序排列，将单位弯制点密度降序排列顺序定义为该变角度域区间内的弯制点弯制顺序，则各变角度域区间的弯制点顺序为b7＝(m17，m16，m18)，b5＝(m11，m12)，b3＝(m6，m7)，b4＝(m9，m8，m10)，b1＝(m1，m3，m4，m2)，b8＝(m20，m19，m22，m21)，b2＝(m5)，b6＝(m14，m13，m15)，结合变角度域区间排列顺序，进而得到正畸弓丝曲线成形控制点坐标矩阵m1＝{m17，m16，m18，m11，m12，m6，m7，m9，m8，m10，m1，m3，m4，m2，m20，m19，m22，m21，m5，m14，m13，m15}和机器人弯制信息集n1＝{n17，n16，n18，n11，n12，n6，n7，n9，n8，n10，n1，n3，n4，n2，n20，n19，n22，n21，n5，n14，n13，n15}，输出最终弯制点弯制顺序m1＝{m17，m16，m18，m11，m12，m6，m7，m9，m8，m10，m1，m3，m4，m2，m20，m19，m22，m21，m5，m14，m13，m15}、n1＝{n17，n16，n18，n11，n12，n6，n7，n9，n8，n10，n1，n3，n4，n2，n20，n19，n22，n21，n5，n14，n13，n15}程序结束。