CN111588504A - 一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法 - Google Patents
一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法 Download PDFInfo
- Publication number
- CN111588504A CN111588504A CN202010486033.7A CN202010486033A CN111588504A CN 111588504 A CN111588504 A CN 111588504A CN 202010486033 A CN202010486033 A CN 202010486033A CN 111588504 A CN111588504 A CN 111588504A
- Authority
- CN
- China
- Prior art keywords
- bending
- ball
- sphere
- points
- arch wire
- 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.)
- Granted
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/12—Brackets; Arch wires; Combinations thereof; Accessories therefor
- A61C7/20—Arch wires
-
- 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
-
- 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/02—Tools for manipulating or working with an orthodontic appliance
- A61C7/026—Tools for manipulating or working with an orthodontic appliance for twisting orthodontic ligature wires
Landscapes
- Health & Medical Sciences (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Epidemiology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
本发明公开一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法,它涉及正畸弓丝弯制技术领域,本发明根据患者的个性化正畸弓丝曲线,基于正畸弓丝曲线弯制点信息集、弯制点的机器人弯制信息集,结合机器人弯制正畸弓丝的运动特点,建立一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法。技术要点为:变半径球域划分数据导入;设定球域限制参数的上限值;确定划分球域的半径和球心;定义合理弯制球域;判断是否继续进行球域划分;以球域弯制点密度为指标,将各球域降序排列,规定球域的弯制点顺序,输出最终弯制点弯制顺序T1和R1。本发明以变半径球域对空间弓丝曲线划分区域,可以提高正畸弓丝弯制规划的准确性,避免了机器人弯制正畸弓丝过程中出现干涉的问题。
Description
技术领域
本发明专利涉及一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法,属于正畸弓丝弯制技术领域。
背景技术
错颌畸形是危及人体健康的第三大口腔疾病,呈现较高的发病率,现代口腔医学中,固定矫治是一种常用且有效的正畸治疗手段,而正畸弓丝的弯制是固定矫治技术的关键,近年来,深受数字化制造技术的影响,传统的口腔制造加工工艺正发生革命性变化,口腔正畸领域也受益于数字化技术,正畸矫治器中弓丝的加工正在向数字化发展。
在机器人弯制个性化正畸弓丝的过程中,个性化正畸弓丝与机器人弯制手爪间可能发生干涉,干涉即个性化正畸弓丝与机器人弯制手爪发生碰撞,干涉发生后会大大影响个性化正畸弓丝的弯制精度,进而影响矫治效果,致使所弯制的个性化正畸弓丝无法应用于临床治疗;研究发现,在正向弯制个性化正畸弓丝的过程中,正向弯制即由未弯制的正畸弓丝弯制成复杂的成形弓丝,干涉往往是由于不合理的成形控制点弯制顺序引起的,而合理的正向弯制成形控制点弯制顺序可以有效地避免干涉的发生,获得合理的正向弯制成形控制点弯制顺序是实现正畸弓丝数字化弯制的必要前提,但目前正畸弓丝弯制技术领域缺少规划正向弯制成形控制点弯制顺序的方法,难以实现正畸弓丝数字化弯制。
对于正畸弓丝弯制规划的划分方法的研究,在本发明人已经授权的授权公告号为CN107647925B的发明专利《一种用于正畸弓丝弯制规划的圆域划分方法》属于一种平面弯制规划方法,采用一种平面等半径圆域划分方法,在正畸弓丝曲线上进行区域划分,最终对每个区域进行排序,从而得出最终弯制点的弯制顺序,虽然这种方法对正畸弓丝弯制规划具有一定应用价值,但是仍然需要做进一步的优化,首先该方法需要在空间正畸弓丝曲线转化为平面曲线的前提下才能进行弯制规划的圆域划分,空间曲线转化为平面曲线后难免会存在个别弯制点投影偏差、错位或相互叠加的情况,从而使得弯制规划后的实际弯制顺序并不能达到理想效果,其次由于该方法仅以一种无依据的均匀化标准对正畸弓丝曲线进行划分,所划分的圆域区间内通常会存在弯制点密集程度过大或过小的情况,即划分后的区间未充分考虑到正畸弓丝曲线上弯制点分布信息的个性化特点,例如当正畸弓丝上存在空间分布的特殊功能曲时,普通的平面弯制规划方法就不能完全发挥其有效性,从而不能有效避免弯制机器人产生空程无效动作或弯制过程相互干涉动作,不利于发挥弯制机器人优势的最大化,不能明显提高弯制效率。
发明内容
针对上述问题,本发明提出一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法,解决目前正畸弓丝弯制技术领域缺少对空间正畸弓丝曲线直接进行弯制规划,进而获得合理正畸弓丝弯制顺序的方法,在弯制顺序规划过程中提出合理的划分依据,对所划分区域的弯制点密集程度及弯制复杂程度定量约束,最终获得符合正畸弓丝曲线上弯制点分布信息个性化特点的弯制顺序,从而发挥弯制机器人优势的最大化,避免了曲线转换导致的弯制顺序规划偏差,保证正畸弓丝弯制过程正常运作,提高正畸弓丝弯制规划的效率,避免了机器人弯制正畸弓丝过程中出现干涉的问题。
上述目的主要通过以下方案实现:
本发明的一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法,所述方法的具体实现过程为:
步骤一、变半径球域划分数据导入:
根据患者有i个弯制点的正畸弓丝曲线,计算并输入正畸弓丝曲线弯制点信息集T={t1,t2,t3,...,ti},ti=(xi,yi,zi)′为每个正畸弓丝曲线弯制点的坐标,在每个弯制点ti上机器人执行不同的弯制运动,每一个正畸弓丝曲线弯制点ti均对应一个弯制点机器人弯制信息单元ri,输入弯制点的机器人弯制信息集为R={r1,r2,r3,...,ri},ri=(xi,yi,zi,αi)′表示机器人在弯制该点时的弯制点坐标及弯制角度,αi为机器人作用在弯制点ti上的弯制角度;
步骤二、球域限制参数的设定
定义球域弯制点个数,用符号表示,球域弯制点个数为半径为的球域an内的弯制点个数;定义球域弯制点密度,用符号表示,球域弯制点密度是球域an内个弯制点在半径为的球域空间内紧密程度的量化描述,规定球域弯制点密度的单位为个/mm3,为正畸弓丝曲线上第n个变半径划分球域an的半径值;定义弯制点角距比,用符号E表示,弯制点角距比是对单个弯制点的弯制复杂程度的量化描述,第j个弯制点的弯制点角距比规定特别地,由于第一个弯制点t1无需弯制,规定弯制点t1的弯制点角距比E1=0;定义球域弯制点角距比和,用符号表示,是对划分球域an所划分的弯制点整体弯制复杂程度的量化描述,表示正畸弓丝曲线上第n个变半径划分球域an内的所有弯制点的弯制点角距比之和,当变半径划分球域an内的弯制点分别为tq+1、tq+2、tq+3、…、时,规定αj为作用在弯制点tj处的弯制角度,表示作用在弯制点tj处弯制距离,即弯制点tj-1与tj之间曲线段的长度,j的取值范围为q表示正畸弓丝曲线上已完成变半径球域划分的球域内的所有弯制点个数,即未进行球域划分时q=0,即q的初始值为q=0,q的取值范围为0≤q≤i,将球域弯制点密度球域弯制点角距比和和球域弯制点个数统称为球域限制参数,分别对球域限制参数的上限值进行限定,设定的上限值Qmax,设定的上限值为ρmax,设定的上限值(∑E)max,在整个球域划分过程中Qmax、ρmax和(∑E)max恒为常数,特别地Qmax=5,跳转至步骤三;
步骤三、确定划分球域的半径和球心:
划分球域以弯制点tq+1为起始点进行,分别取弯制点tq+1与tq+2、tq+3、…、之间的直线段,依次记为线段tq+1tq+2、tq+1tq+3、…、将线段tq+1tq+2、tq+1tq+3、…、中长度最大的线段记为即分别表示线段tq+ 1tq+2、tq+1tq+3、…、的长度,则在空间正畸弓丝曲线上即将生成的第n个划分球域an的球心为线段的中点,半径为线段长度的一半此时恰好有两个弯制点落到球域边界面上,且新生成的球域an刚好能划分步骤三中预先规定的弯制点tq+1、tq+2、tq+3、…、规定球域边界面所截的空间正畸弓丝曲线段上的所有弯制点被该球域所划分,当生成的球域边界面通过弯制点时,则该弯制点也被球域所划分,已被划分的弯制点所在的正畸弓丝曲线段将不会再被其他球域划分;的初始值为n的初始值为n=1,即首次划分第1个球域a1时预先规定该球域划分到的弯制点刚好达到上限值,此时所能划分到的弯制点分别为t1、t2、t3、t4、t5,且t1为划分球域a1的起始点;
步骤四、定义合理弯制球域:
具体为:
如果成立且成立,说明以直线段的中点为球心,以为半径的划分球域的球域弯制点密度没有超过所设定的球域弯制点密度上限值ρmax且球域弯制点角距比和没有超过所设定的球域弯制点角距比和上限值(∑E)max,则将以直线段的中点为球心,以为半径的包含正畸弓丝曲线段的划分球域定义为合理弯制球域an,计算正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q,跳转至步骤五;
如果不成立或不成立,且存在说明此时划分球域的球域弯制点个数仅为1个,将以tq+1为球心,以弯制点tq+1与相邻弯制点tq+2之间直线距离的一半为半径生成的包含空间正畸弓丝曲线段的划分球域定义为合理弯制球域an,则该合理弯制球域an上仅包含一个弯制点tq+1,计算空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q,跳转至步骤五;
步骤五、判断是否继续进行球域划分:
判断空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q是否与弯制点个数i相等,
具体为:
如果空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q与弯制点个数i不相等,则继续进行球域划分,令n=n+1,即表示对下一个球域进行划分,此时,
若i-q<5且i-q≠1,说明此时正畸弓丝曲线上剩余未被划分的弯制点少于5个但超过1个,则令即进行下一个球域的首次划分时令球域能划分到的弯制点个数为正畸弓丝曲线上剩余未被划分的弯制点个数,跳转至步骤三;
若i-q<5且i-q=1,说明此时正畸弓丝曲线上剩余未被划分的弯制点只有最后1个弯制点ti,将以ti为球心,以ti-1与ti之间直线距离的一半为半径生成的包含空间正畸弓丝曲线段的划分球域定义为合理弯制球域an,则该合理弯制球域an内仅包含一个弯制点ti,跳转至步骤六;
如果正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q与弯制点个数i相等,说明所有的弯制点均已被合理弯制球域划分,输出合理弯制球域信息集合A1={a1,a2,...,an},跳转至步骤六;
步骤六、输出最终弯制顺序:
计算每个合理弯制球域(a1,a2,...,an)的球域弯制点密度获得球域弯制点密度信息集比较各合理弯制球域的球域弯制点密度,假设得到则以球域弯制点密度为指标将n个球域降序排列,从而得到降序合理弯制球域信息集为A1={a3,a1,...,as},规定在任何一个弯制球域,以每个弯制点角距比E为指标对每个球域所划分的弯制点进行降序排列,将降序弯制点角距比所对应弯制点的顺序定义为该球域弯制点的弯制顺序,进而得到正畸弓丝曲线成形控制点坐标矩阵T1={t7,t8,...,tm}和机器人弯制信息集R1={r7,r8,...,rm},其中tm表示第s个合理弯制球域中弯制点角距比最小的弯制点,输出最终弯制点弯制顺序T1={t7,t8,...,tm}、R1={r7,r8,...,rm},程序结束。
本发明的有益效果为:
1、本发明针对空间正畸弓丝弯制规划方法,对球域限制参数进行定义,提出了球域弯制点个数和球域弯制点密度的概念,在球域空间内对弯制点的紧密程度量化描述,提出了球域弯制点角距比和的概念,对划分球域内弯制点弯制复杂程度进行量化描述,设定的上限值Qmax,可保证每个划分球域弯制点个数不大于上限值,设定的上限值ρmax和的上限值(∑E)max,从而可在一条正畸弓丝曲线上形成若干个符合设定要求的合理弯制球域,避免了所划分的球域出现弯制点密度过大、弯制复杂程度过高的现象,从而最大程度避免了机器人在弯制过程发生干涉的问题,通过对等球域限制参数的定义,便于在空间正畸弓丝弯制规划中,利用球域限制参数作为一种规划指标进行划分弯制点,为空间正畸弓丝弯制规划提供了理论依据。
2、本发明采用变半径球域划分的方法,划分区域的球心和半径是由空间区域内的弯制点的分布状况所决定的,空间区域内弯制点的个数的变化引起了球域的球心以及半径的改变,能够更加充分的反映出区域弯制点的信息,便于系统对划分球域数据的计算及选取,提高了本方法的划分效率。
3、本发明在划分完所有弯制点后,以针对球域划分定义的球域弯制点密度为指标对n个球域进行降序排列,得到了降序合理弯制球域信息集,使得弯制球域的排列为由密到疏,从最大程度上避免干涉的产生,并规定在任何一个弯制球域内,以降序弯制点角距比所对应弯制点的顺序定义为该球域弯制点的弯制顺序,使得在每一个球域内的弯制点的弯制难度由复杂到简单过渡,从而使得机器人弯制的优势最大化,保证了每个弯制点有唯一确定的弯制顺序,提高了正畸弓丝弯制规划的可操作性和准确性。
4、同本发明人已经授权的授权公告号为CN107647925B的发明专利《一种用于正畸弓丝弯制规划的圆域划分方法》相比,本专利所提出的方法无需转化空间曲线就能进行弯制规划,避免了个别弯制点投影偏差、错位或相互叠加的情况,使得弯制顺序规划后的规划顺序更加准确,实际弯制效果更加理想,另外依据所定义的球域限制参数对所划分区域的弯制点进行弯制复杂程度及密集程度的定量约束,从而导致划分球域的半径为符合球域限制参数的规定而发生变化,最终得到符合正畸弓丝曲线上弯制点分布信息个性化特点的不同半径的合理弯制球域,有效避免完成划分后的各球域内弯制点的密集程度、弯制复杂程度相差悬殊的情况,从而可以发挥弯制机器人优势的最大化,保证正畸弓丝弯制过程正常运作,提高正畸弓丝弯制规划的效率,避免了机器人弯制正畸弓丝过程中出现干涉的问题。
5、同本发明人同日申请的发明专利《一种用于正畸弓丝弯制顺序规划的平面变半径圆域划分方法》相比,两方法都对区域划分的限制参数进行了定义,且划分思路大致相同,但本方法是基于空间正畸弓丝曲线进行弯制顺序规划,无需对正畸弓丝曲线进行平面转化,因此避免了个别弯制点投影偏差、错位或相互叠加的情况,保证了规划后的弯制顺序精准度更高,实际弯制效果更加理想,从而进一步发挥弯制机器人优势的最大化,避免了机器人弯制正畸弓丝过程中出现干涉的问题。
附图说明
为了易于说明,本发明由下述的具体实施及附图作以详细描述。
图1为一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法流程图;
图2为个性化正畸弓丝弯制点分布示意图;
图3为空间变半径球域划分个性化正畸弓丝曲线示意图;
具体实施方式
为使本发明专利的目的、技术方案和优点更加清楚明了,下面通过附图中示出的具体实施例来描述本发明专利,但是应该理解,这些描述只是示例性的,而并非要限制本发明专利的范围,此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明专利的概念。
实施例1:如图1、图2、图3所示,本具体实施方式采用以下技术方案:一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法,所述方法的具体实现过程为:
步骤一、变半径球域划分数据导入:
根据患者有i个弯制点的正畸弓丝曲线,计算并输入正畸弓丝曲线弯制点信息集T={t1,t2,t3,...,ti},ti=(xi,yi,zi)′为每个正畸弓丝曲线弯制点的坐标,在每个弯制点ti上机器人执行不同的弯制运动,每一个正畸弓丝曲线弯制点ti均对应一个弯制点机器人弯制信息单元ri,输入弯制点的机器人弯制信息集为R={r1,r2,r3,...,ri},ri=(xi,yi,zi,αi)′表示机器人在弯制该点时的弯制点坐标及弯制角度,αi为机器人作用在弯制点ti上的弯制角度;
步骤二、球域限制参数的设定
定义球域弯制点个数,用符号表示,球域弯制点个数为半径为的球域an内的弯制点个数;定义球域弯制点密度,用符号表示,球域弯制点密度是球域an内个弯制点在半径为的球域空间内紧密程度的量化描述,规定球域弯制点密度的单位为个/mm3,为正畸弓丝曲线上第n个变半径划分球域an的半径值;定义弯制点角距比,用符号E表示,弯制点角距比是对单个弯制点的弯制复杂程度的量化描述,第j个弯制点的弯制点角距比规定特别地,由于第一个弯制点t1无需弯制,规定弯制点t1的弯制点角距比E1=0;定义球域弯制点角距比和,用符号表示,是对划分球域an所划分的弯制点整体弯制复杂程度的量化描述,表示正畸弓丝曲线上第n个变半径划分球域an内的所有弯制点的弯制点角距比之和,当变半径划分球域an内的弯制点分别为tq+1、tq+2、tq+3、…、时,规定αj为作用在弯制点tj处的弯制角度,表示作用在弯制点tj处弯制距离,即弯制点tj-1与tj之间曲线段的长度,j的取值范围为q表示正畸弓丝曲线上已完成变半径球域划分的球域内的所有弯制点个数,即未进行球域划分时q=0,即q的初始值为q=0,q的取值范围为0≤q≤i,将球域弯制点密度球域弯制点角距比和和球域弯制点个数统称为球域限制参数,分别对球域限制参数的上限值进行限定,设定的上限值Qmax,设定的上限值为ρmax,设定的上限值(∑E)max,在整个球域划分过程中Qmax、ρmax和(∑E)max恒为常数,特别地Qmax=5,跳转至步骤三;
步骤三、确定划分球域的半径和球心:
划分球域以弯制点tq+1为起始点进行,分别取弯制点tq+1与tq+2、tq+3、…、之间的直线段,依次记为线段tq+1tq+2、tq+1tq+3、…、将线段tq+1tq+2、tq+1tq+3、…、中长度最大的线段记为即分别表示线段tq+ 1tq+2、tq+1tq+3、…、的长度,则在空间正畸弓丝曲线上即将生成的第n个划分球域an的球心为线段的中点,半径为线段长度的一半此时恰好有两个弯制点落到球域边界面上,且新生成的球域an刚好能划分步骤三中预先规定的弯制点tq+1、tq+2、tq+3、…、规定球域边界面所截的空间正畸弓丝曲线段上的所有弯制点被该球域所划分,当生成的球域边界面通过弯制点时,则该弯制点也被球域所划分,已被划分的弯制点所在的正畸弓丝曲线段将不会再被其他球域划分;的初始值为n的初始值为n=1,即首次划分第1个球域a1时预先规定该球域划分到的弯制点刚好达到上限值,此时所能划分到的弯制点分别为t1、t2、t3、t4、t5,且t1为划分球域a1的起始点;
步骤四、定义合理弯制球域:
具体为:
如果成立且成立,说明以直线段的中点为球心,以为半径的划分球域的球域弯制点密度没有超过所设定的球域弯制点密度上限值ρmax且球域弯制点角距比和没有超过所设定的球域弯制点角距比和上限值(∑E)max,则将以直线段的中点为球心,以为半径的包含正畸弓丝曲线段的划分球域定义为合理弯制球域an,计算正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q,跳转至步骤五;
如果不成立或不成立,且存在说明此时划分球域的球域弯制点个数仅为1个,将以tq+1为球心,以弯制点tq+1与相邻弯制点tq+2之间直线距离的一半为半径生成的包含空间正畸弓丝曲线段的划分球域定义为合理弯制球域an,则该合理弯制球域an上仅包含一个弯制点tq+1,计算空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q,跳转至步骤五;
步骤五、判断是否继续进行球域划分:
判断空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q是否与弯制点个数i相等,
具体为:
如果空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q与弯制点个数i不相等,则继续进行球域划分,令n=n+1,即表示对下一个球域进行划分,此时,
若i-q<5且i-q≠1,说明此时正畸弓丝曲线上剩余未被划分的弯制点少于5个但超过1个,则令即进行下一个球域的首次划分时令球域能划分到的弯制点个数为正畸弓丝曲线上剩余未被划分的弯制点个数,跳转至步骤三;
若i-q<5且i-q=1,说明此时正畸弓丝曲线上剩余未被划分的弯制点只有最后1个弯制点ti,将以ti为球心,以ti-1与ti之间直线距离的一半为半径生成的包含空间正畸弓丝曲线段的划分球域定义为合理弯制球域an,则该合理弯制球域an内仅包含一个弯制点ti,跳转至步骤六;
如果正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q与弯制点个数i相等,说明所有的弯制点均已被合理弯制球域划分,输出合理弯制球域信息集合A1={a1,a2,…,an},跳转至步骤六;
步骤六、输出最终弯制顺序:
计算每个合理弯制球域(a1,a2,…,an)的球域弯制点密度获得球域弯制点密度信息集比较各合理弯制球域的球域弯制点密度,假设得到则以球域弯制点密度为指标将n个球域降序排列,从而得到降序合理弯制球域信息集为A1={a3,a1,…,as},规定在任何一个弯制球域,以每个弯制点角距比E为指标对每个球域所划分的弯制点进行降序排列,将降序弯制点角距比所对应弯制点的顺序定义为该球域弯制点的弯制顺序,进而得到正畸弓丝曲线成形控制点坐标矩阵T1={t7,t8,…,tm}和机器人弯制信息集R1={r7,r8,…,rm},其中tm表示第s个合理弯制球域中弯制点角距比最小的弯制点,输出最终弯制点弯制顺序T1={t7,t8,…,tm}、R1={r7,r8,…,rm},程序结束。
实施例2:如图2、图3所示,在一条包含i=18个弯制点的个性化正畸弓丝曲线进行平面变半径球域划分的正畸弓丝弯制顺序规划过程中,假设最终得到的合理弯制球域的个数为n=6,各球域弯制点个数分别为在步骤六中计算每个合理弯制球域(a1,a2,…,an)的球域弯制点密度获得球域弯制点密度信息集比较各合理弯制球域的球域弯制点密度,存在以球域弯制点密度为指标将n个球域降序排列,从而得到降序合理弯制球域信息集为A1={a5,a2,a3,a4,a1,a6},规定在任何一个弯制球域,以降序弯制点角距比所对应弯制点的顺序定义为该球域弯制点的弯制顺序,此时各合理弯制球域的弯制点顺序为a5=(t14,t15,t13),a2=(t5,t4,t6),a3=(t10,t8,t7,t9),a4=(t11,t12),a1=(t3,t2,t1),a6=(t13,t14,t15),根据球域内弯制点的排序,进而可得到正畸弓丝曲线成形控制点坐标矩阵T1={t14,t15,t13,t5,t4,t6,t10,t8,t7,t9,t11,t12,t3,t2,t1,t13,t14,t15}和机器人弯制信息集R1={r14,r15,r13,r5,r4,r6,r10,r8,r7,r9,r11,r12,r3,r2,r1,r13,r14,r15},输出最终弯制点弯制顺序T1={t14,t15,t13,t5,t4,t6,t10,t8,t7,t9,t11,t12,t3,t2,t1,t13,t14,t15}、R1={r14,r15,r13,r5,r4,r6,r10,r8,r7,r9,r11,r12,r3,r2,r1,r13,r14,r15},程序结束。
Claims (1)
1.一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法,其特征在于:所述方法的具体实现过程为:
步骤一、变半径球域划分数据导入:
根据患者有i个弯制点的正畸弓丝曲线,计算并输入正畸弓丝曲线弯制点信息集T={t1,t2,t3,…,ti},ti=(xi,yi,zi)′为每个正畸弓丝曲线弯制点的坐标,在每个弯制点ti上机器人执行不同的弯制运动,每一个正畸弓丝曲线弯制点ti均对应一个弯制点机器人弯制信息单元ri,输入弯制点的机器人弯制信息集为R={r1,r2,r3,…,ri},ri=(xi,yi,zi,αi)′表示机器人在弯制该点时的弯制点坐标及弯制角度,αi为机器人作用在弯制点ti上的弯制角度;
步骤二、球域限制参数的设定
定义球域弯制点个数,用符号表示,球域弯制点个数为半径为的球域an内的弯制点个数;定义球域弯制点密度,用符号表示,球域弯制点密度是球域an内个弯制点在半径为的球域空间内紧密程度的量化描述,规定球域弯制点密度的单位为个/mm3,为正畸弓丝曲线上第n个变半径划分球域an的半径值;定义弯制点角距比,用符号E表示,弯制点角距比是对单个弯制点的弯制复杂程度的量化描述,第j个弯制点的弯制点角距比规定特别地,由于第一个弯制点t1无需弯制,规定弯制点t1的弯制点角距比E1=0;定义球域弯制点角距比和,用符号表示,是对划分球域an所划分的弯制点整体弯制复杂程度的量化描述,表示正畸弓丝曲线上第n个变半径划分球域an内的所有弯制点的弯制点角距比之和,当变半径划分球域an内的弯制点分别为tq+1、时,规定αj为作用在弯制点tj处的弯制角度,表示作用在弯制点tj处弯制距离,即弯制点tj-1与tj之间曲线段的长度,j的取值范围为q表示正畸弓丝曲线上已完成变半径球域划分的球域内的所有弯制点个数,即未进行球域划分时q=0,即q的初始值为q=0,q的取值范围为0≤q≤i,将球域弯制点密度球域弯制点角距比和和球域弯制点个数统称为球域限制参数,分别对球域限制参数的上限值进行限定,设定的上限值Qmax,设定的上限值为ρmax,设定的上限值(∑E)max,在整个球域划分过程中Qmax、ρmax和(∑E)max恒为常数,特别地Qmax=5,跳转至步骤三;
步骤三、确定划分球域的半径和球心:
划分球域以弯制点tq+1为起始点进行,分别取弯制点tq+1与之间的直线段,依次记为线段将线段中长度最大的线段记为即 分别表示线段 的长度,则在空间正畸弓丝曲线上即将生成的第n个划分球域an的球心为线段的中点,半径为线段长度的一半此时恰好有两个弯制点落到球域边界面上,且新生成的球域an刚好能划分步骤三中预先规定的弯制点规定球域边界面所截的空间正畸弓丝曲线段上的所有弯制点被该球域所划分,当生成的球域边界面通过弯制点时,则该弯制点也被球域所划分,已被划分的弯制点所在的正畸弓丝曲线段将不会再被其他球域划分;的初始值为n的初始值为n=1,即首次划分第1个球域a1时预先规定该球域划分到的弯制点刚好达到上限值,此时所能划分到的弯制点分别为t1、t2、t3、t4、t5,且t1为划分球域a1的起始点;
步骤四、定义合理弯制球域:
具体为:
如果成立且成立,说明以直线段的中点为球心,以为半径的划分球域的球域弯制点密度没有超过所设定的球域弯制点密度上限值ρmax且球域弯制点角距比和没有超过所设定的球域弯制点角距比和上限值(∑E)max,则将以直线段的中点为球心,以为半径的包含正畸弓丝曲线段的划分球域定义为合理弯制球域an,计算正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q,跳转至步骤五;
如果不成立或不成立,且存在说明此时划分球域的球域弯制点个数仅为1个,将以tq+1为球心,以弯制点tq+1与相邻弯制点tq+2之间直线距离的一半为半径生成的包含空间正畸弓丝曲线段的划分球域定义为合理弯制球域an,则该合理弯制球域an上仅包含一个弯制点tq+1,计算空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q,跳转至步骤五;
步骤五、判断是否继续进行球域划分:
判断空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q是否与弯制点个数i相等,
具体为:
如果空间正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q与弯制点个数i不相等,则继续进行球域划分,令n=n+1,即表示对下一个球域进行划分,此时,
若i-q<5且i-q≠1,说明此时正畸弓丝曲线上剩余未被划分的弯制点少于5个但超过1个,则令即进行下一个球域的首次划分时令球域能划分到的弯制点个数为正畸弓丝曲线上剩余未被划分的弯制点个数,跳转至步骤三;
若i-q<5且i-q=1,说明此时正畸弓丝曲线上剩余未被划分的弯制点只有最后1个弯制点ti,将以ti为球心,以ti-1与ti之间直线距离的一半为半径生成的包含空间正畸弓丝曲线段的划分球域定义为合理弯制球域an,则该合理弯制球域an内仅包含一个弯制点ti,跳转至步骤六;
如果正畸弓丝曲线上所有已被合理弯制球域划分的弯制点个数q与弯制点个数i相等,说明所有的弯制点均已被合理弯制球域划分,输出合理弯制球域信息集合A1={a1,a2,…,an},跳转至步骤六;
步骤六、输出最终弯制顺序
计算每个合理弯制球域(a1,a2,…,an)的球域弯制点密度获得球域弯制点密度信息集比较各合理弯制球域的球域弯制点密度,假设得到则以球域弯制点密度为指标将n个球域降序排列,从而得到降序合理弯制球域信息集为A1={a3,a1,…,as},规定在任何一个弯制球域,以每个弯制点角距比E为指标对每个球域所划分的弯制点进行降序排列,将降序弯制点角距比所对应弯制点的顺序定义为该球域弯制点的弯制顺序,进而得到正畸弓丝曲线成形控制点坐标矩阵T1={t7,t8,…,tm}和机器人弯制信息集R1={r7,r8,...,rm},其中tm表示第s个合理弯制球域中弯制点角距比最小的弯制点,输出最终弯制点弯制顺序T1={t7,t8,...,tm}、R1={r7,r8,...,rm},程序结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010486033.7A CN111588504B (zh) | 2020-06-01 | 2020-06-01 | 一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010486033.7A CN111588504B (zh) | 2020-06-01 | 2020-06-01 | 一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111588504A true CN111588504A (zh) | 2020-08-28 |
CN111588504B CN111588504B (zh) | 2021-05-28 |
Family
ID=72188201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010486033.7A Active CN111588504B (zh) | 2020-06-01 | 2020-06-01 | 一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111588504B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114983593A (zh) * | 2022-04-18 | 2022-09-02 | 哈尔滨理工大学 | 一种基于向量共线的正畸弓丝弯制点误差评价方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040029067A1 (en) * | 2002-08-09 | 2004-02-12 | Wool Arthur L. | Three segment orthodontic arch wire having uniform flexural rigidity |
EP1744698A1 (de) * | 2004-04-28 | 2007-01-24 | Christoph Von Mandach | Orthodontisches bracket |
CN106803276A (zh) * | 2017-01-16 | 2017-06-06 | 哈尔滨理工大学 | 一种个性化正畸弓丝曲线交互调整方法 |
CN106803018A (zh) * | 2017-01-16 | 2017-06-06 | 哈尔滨理工大学 | 一种个性化正畸弓丝参数化表达方法 |
CN107647925A (zh) * | 2017-11-14 | 2018-02-02 | 哈尔滨理工大学 | 一种用于正畸弓丝弯制规划的圆域划分方法 |
CN107714203A (zh) * | 2017-11-14 | 2018-02-23 | 哈尔滨理工大学 | 一种等角度划分的正畸弓丝弯制顺序规划方法 |
CN108746415A (zh) * | 2018-06-16 | 2018-11-06 | 哈尔滨理工大学 | 一种利用机器人弯制带圈泪滴曲的弯制规划方法 |
-
2020
- 2020-06-01 CN CN202010486033.7A patent/CN111588504B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040029067A1 (en) * | 2002-08-09 | 2004-02-12 | Wool Arthur L. | Three segment orthodontic arch wire having uniform flexural rigidity |
EP1744698A1 (de) * | 2004-04-28 | 2007-01-24 | Christoph Von Mandach | Orthodontisches bracket |
CN106803276A (zh) * | 2017-01-16 | 2017-06-06 | 哈尔滨理工大学 | 一种个性化正畸弓丝曲线交互调整方法 |
CN106803018A (zh) * | 2017-01-16 | 2017-06-06 | 哈尔滨理工大学 | 一种个性化正畸弓丝参数化表达方法 |
CN107647925A (zh) * | 2017-11-14 | 2018-02-02 | 哈尔滨理工大学 | 一种用于正畸弓丝弯制规划的圆域划分方法 |
CN107714203A (zh) * | 2017-11-14 | 2018-02-23 | 哈尔滨理工大学 | 一种等角度划分的正畸弓丝弯制顺序规划方法 |
CN108746415A (zh) * | 2018-06-16 | 2018-11-06 | 哈尔滨理工大学 | 一种利用机器人弯制带圈泪滴曲的弯制规划方法 |
CN109894545A (zh) * | 2018-06-16 | 2019-06-18 | 哈尔滨理工大学 | 一种利用机器人弯制带圈泪滴曲的弯制规划方法 |
Non-Patent Citations (1)
Title |
---|
姜金刚等: "《基于有限点展成法的正畸弓丝成形控制点规划》", 《仪器仪表学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114983593A (zh) * | 2022-04-18 | 2022-09-02 | 哈尔滨理工大学 | 一种基于向量共线的正畸弓丝弯制点误差评价方法 |
CN114983593B (zh) * | 2022-04-18 | 2023-11-17 | 哈尔滨理工大学 | 一种基于向量共线的正畸弓丝弯制点误差评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111588504B (zh) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111588502B (zh) | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 | |
CN107647925B (zh) | 一种用于正畸弓丝弯制规划的圆域划分方法 | |
CN107714203B (zh) | 一种等角度划分的正畸弓丝弯制顺序规划方法 | |
CN111588499B (zh) | 一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法 | |
CN107822722B (zh) | 一种用于正畸弓丝弯制运动规划的有限点寻回展成方法 | |
CN111588504B (zh) | 一种用于正畸弓丝弯制顺序规划的空间变半径球域划分方法 | |
CN116644558A (zh) | 一种基于误差评定域的正畸弓丝误差评价方法 | |
CN111588497B (zh) | 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法 | |
CN114943685A (zh) | 一种基于轮廓降维法的正畸弓丝误差评价方法 | |
CN111588496B (zh) | 一种用于正畸弓丝弯制规划的平面变角度划分方法 | |
CN111588494B (zh) | 一种基于弯制点密度的正畸弓丝变角度划分方法 | |
CN111588505B (zh) | 一种用于正畸弓丝弯制顺序规划的平面变半径圆域划分方法 | |
CN111588503B (zh) | 一种基于弯制点密度的正畸弓丝变半径圆域划分方法 | |
CN111588491B (zh) | 一种基于正畸弓丝弯制点密度的空间等半径球域划分半径确定方法 | |
CN111588501B (zh) | 一种正畸弓丝弯制规划等半径圆域划分半径确定方法 | |
CN115035196A (zh) | 一种基于弯制点复杂度判断的正畸弓丝误差率评价方法 | |
CN114972184B (zh) | 一种基于权值比例法的正畸弓丝误差评价方法 | |
CN115040275A (zh) | 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法 | |
CN115024840A (zh) | 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 | |
CN111588493B (zh) | 一种基于弯制点角距比和的正畸弓丝变角度划分方法 | |
CN115024841A (zh) | 一种基于复杂度区分法的正畸弓丝误差评价方法 | |
CN111588500B (zh) | 一种正畸弓丝弯制顺序规划等角度划分角度确定方法 | |
CN111588495B (zh) | 一种基于正畸弓丝弯制点单位角距比和的等角度划分角度确定方法 | |
CN111588498B (zh) | 一种基于正畸弓丝弯制点密度的等角度划分角度确定方法 | |
CN114943058B (zh) | 一种基于位置误差判断的正畸弓丝误差波动度评价方法 |
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 |