CN115024838B - 一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 - Google Patents
一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 Download PDFInfo
- Publication number
- CN115024838B CN115024838B CN202210402124.7A CN202210402124A CN115024838B CN 115024838 B CN115024838 B CN 115024838B CN 202210402124 A CN202210402124 A CN 202210402124A CN 115024838 B CN115024838 B CN 115024838B
- Authority
- CN
- China
- Prior art keywords
- orthodontic archwire
- bending point
- curve
- theoretical
- actual
- 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
- 238000005452 bending Methods 0.000 title claims abstract description 575
- 238000011156 evaluation Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 claims description 50
- 230000009191 jumping Effects 0.000 claims description 26
- 230000009466 transformation Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 6
- 210000004513 dentition Anatomy 0.000 claims description 3
- 230000036346 tooth eruption Effects 0.000 claims description 3
- 238000011158 quantitative evaluation Methods 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000012854 evaluation process Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 208000025157 Oral disease Diseases 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 208000030194 mouth disease Diseases 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/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/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
- A61C2007/004—Automatic construction of a set of axes for a tooth or a plurality of teeth
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
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
本发明公开了一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法,它涉及正畸弓丝弯制评价领域,本发明基于空间变换后的理论正畸弓丝曲线弯制点信息集、空间变换后的实际正畸弓丝曲线信息集,首前先进行理论正畸弓丝曲线弯制点复杂度的计算方法进而判断弯制点特征,结合实际正畸弓丝曲线弯制点的曲率误差波动度、线误差波动度和角误差波动度,建立一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法。本发明针具有不同特征的正畸弓丝弯制点采用误差波动度评价方案,提高了该正畸弓丝评价方法的适用范围,并通过计算实际正畸弓丝曲线弯制点相应应评价的曲率误差波动度、线误差波动度和角误差波动度,实现了对正畸弓丝弯制稳定性的量化评估。
Description
技术领域
本发明专利涉及一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法,属于正畸弓丝弯制评价技术领域。
背景技术
错颌畸形是危及人体健康的第三大口腔疾病,呈现较高的发病率,现代口腔医学中,固定矫治是一种常用且有效的正畸治疗手段,而正畸弓丝的弯制是固定矫治技术的关键,在传统临床应用中,正畸弓丝基本依赖于专业医师手工弯制,难以保证精度;近年来,深受数字化制造技术的影响,传统的口腔制造加工工艺正发生革命性变化,口腔正畸领域也受益于数字化技术,正畸矫治器中弓丝的加工正在向数字化发展。但是现阶段在正畸弓丝弯制完成后,仍然需要医师依据经验对正畸弓丝进行评价,判断是否满足使用要求,该方法严重依赖医师的临床经验,难以实现正畸弓丝的量化评价,且仅能实现对单个正畸弓丝曲线弯制的准确性评价,而无法通过分析大量依据同一理论正畸弓丝曲线弯制出的实际正畸弓丝曲线的误差率数据进而实现对正畸弓丝弯制稳定性的量化评价;此类方法的缺失导致医师和正畸弓丝机器人无法获得具有针对性的指导和改进方向,阻碍了医师弯制技术的进步,严重影响了临床矫治效果,阻碍了正畸弓丝机器人机械结构的优化限制了正畸弓丝弯制算法的迭代升级,制约了正畸弓丝弯制机器人的发展;综上,目前正畸弓丝弯制评价技术领域亟待一种能够精确量化评价正畸弓丝曲线弯制稳定性的方法。
发明内容
针对上述问题,本发明提出了一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法,解决目前正畸弓丝评价领域中缺少针对正畸弓丝弯制稳定性的量化评价方法,在正畸弓丝评价过程中对正畸弓丝弯制点提出合理的判断依据,并针对具有不同特征的正畸弓丝弯制点进行不同的误差波动度评价,提高了该正畸弓丝评价方法的适用范围,并通过正畸弓丝弯制点相应误差波动度的计算,完成对正畸弓丝弯制稳定性的量化评估。
一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法,所述方法的具体实现过程为:
步骤一、理论正畸弓丝曲线数据及实际正畸弓丝曲线数据导入:
以右手定则建立o-xyz三维正畸弓丝误差标定坐标系w,以正畸医师根据患者牙列形态设计的具有n个弯制点的理论正畸弓丝曲线,计算并输入理论正畸弓丝曲线弯制点信息集P'T={Tp'1,Tp'2,Tp'3,...,Tp'i,...,Tp'n},Tp'i=(Tα'i,Tβ'i',Tγ'i,Td'i)为理论正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,i的取值范围为1≤i≤n,Tα'i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Tβ'i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Tγ'i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Td'i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;理论正畸弓丝曲线左端点为ps,理论正畸弓丝曲线右端点为pf,ps和pf之间连线的中点为To',对理论正畸弓丝曲线进行空间变换:令点To'与三维正畸弓丝误差标定坐标系w的原点o重合,理论正畸弓丝曲线左端点ps位于y轴负半轴,理论正畸弓丝曲线右端点pf位于y轴正半轴,且理论正畸弓丝曲线与x轴正半轴无交点,再令理论正畸弓丝曲线沿y轴正方向顺时针旋转,直至理论正畸弓丝曲线与x轴正半轴出现交点,将理论正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系w中的最终位姿,计算并输入处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpi,...,Tpn},Tpi=(Tαi,Tβi,Tγi,Tdi)为处于最终位姿下的理论正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Tαi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Tβi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Tγi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Tdi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
以根据理论正畸弓丝曲线弯制出的m条具有n个弯制点的实际正畸弓丝曲线,计算并输入实际正畸弓丝曲线信息集为第j条实际正畸弓丝曲线弯制点信息集,j的取值范围为1≤j≤m,为第j条实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;第j条实际正畸弓丝曲线左端点为jp's,第j条实际正畸弓丝曲线右端点为jp'f,jp's和jp'f之间连线的中点为对第j条实际正畸弓丝曲线进行空间变换:令点与三维正畸弓丝误差标定坐标系w的原点o重合,第j条实际弓丝曲线左端点jp's位于y轴负半轴,第j条实际正畸弓丝曲线右端点jp'f位于y轴正半轴,且第j条实际正畸弓丝曲线与x轴正半轴无交点,再令第j条实际正畸弓丝曲线沿y轴正方向顺时针旋转,直至第j条实际正畸弓丝曲线与x轴正半轴出现交点,将第j条实际正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系中的最终位姿,计算并输入m条处于最终位姿下的实际正畸弓丝曲线信息集为第j条处于最终位姿下的实际正畸弓丝曲线弯制点信息集,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
步骤二、实际正畸弓丝曲线弯制点复杂度的设定:
定义理论正畸弓丝曲线弯制点复杂度,用符号Cr表示,Cr是对理论正畸弓丝曲线弯制点的弯制难易程度的综合量化描述,理论正畸弓丝曲线弯制点的Cr值越高,即该弯制点在弯制时越困难,规定理论正畸弓丝曲线第i个弯制点的复杂度表示为 表示理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比,规定 TEi表示理论正畸弓丝曲线第i个弯制点的弯制点角距比,弯制点角距比是对正畸弓丝曲线上单个弯制点的弯制复杂程度的量化描述,规定 Tθi为作用在理论正畸弓丝曲线弯制点Tpi处的弯制角度,表示作用在理论正畸弓丝曲线第i个弯制点处的弯制距离,即理论正畸弓丝曲线弯制点Tpi-1与Tpi之间的曲线段的长度,对于理论正畸弓丝曲线第1个弯制点Tp1,表示弯制点Tp1到理论正畸弓丝曲线左端点ps之间的曲线段长度,TEmin为理论正畸弓丝曲线弯制点角距比的最小值,TEmax为理论正畸弓丝曲线弯制点角距比的最大值;规定的分界值为(TE*)b;表示理论正畸弓丝曲线第i个弯制点的归一化弯制点密度,规定 Tρi表示理论正畸弓丝曲线第i个弯制点的弯制点密度,弯制点密度是对正畸弓丝曲线上单个弯制点与相邻弯制点间紧密程度的量化描述,规定公式中的数值1表示为1个弯制点,Tli表示理论正畸弓丝曲线第i个弯制点与其距离最近的弯制点之间的直线距离,即 表示理论正畸弓丝曲线第i-1个弯制点与理论正畸弓丝曲线第i个弯制点之间的距离,表示理论正畸弓丝曲线第i个弯制点与理论正畸弓丝曲线第i+1个弯制点之间的距离,当i=1时,规定 表示理论正畸弓丝曲线第1个弯制点与理论正畸弓丝曲线左端点ps之间的直线距离,表示理论正畸弓丝曲线第1个弯制点与理论正畸弓丝曲线第2个弯制点之间的直线距离,当i=n时,规定 表示理论正畸弓丝曲线第n-1个弯制点与理论正畸弓丝曲线第n个弯制点之间的直线距离,表示理论正畸弓丝曲线第n个弯制点与理论正畸弓丝曲线右端点pf之间的直线距离,Tρmin为理论正畸弓丝曲线弯制点角距比的最小值,Tρmax为理论正畸弓丝曲线弯制点角距比的最大值;复杂度公式中的数值2表示在计算复杂度时考虑了归一化弯制点角距比和归一化弯制单密度两个参数;规定理论正畸弓丝曲线弯制点的复杂度Cr的分界值为Crb;
步骤三、实际正畸弓丝曲线曲率误差波动度、线误差波动度、角误差波动度的设定:
定义实际正畸弓丝曲线的曲率误差波动度,用符号σc表示,σc是对实际正畸弓丝曲线弯制曲率稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的曲率误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的曲率误差率,规定 TKi表示理论正畸弓丝曲线第i个弯制点的曲率,表示第j条实际正畸弓丝曲线第i个弯制点的曲率,表示m条实际正畸弓丝曲线的第i个弯制点的曲率误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度的上限值为(σc)max;定义实际正畸弓丝曲线的线误差波动度,用符号σd表示,σd是对实际正畸弓丝曲线弯制距离稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的线误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的线误差率,规定 表示m条实际正畸弓丝曲线的第i个弯制点的线误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的线误差波动度的上限值为(σd)max;定义实际正畸弓丝曲线的角误差波动度,用符号σa表示,σa是对实际正畸弓丝曲线弯制角度稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的角误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的角误差率,规定其中为理论正畸弓丝曲线第i个弯制点的角Tαi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 为理论正畸弓丝曲线第i个弯制点的角Tβi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 为理论正畸弓丝曲线第i个弯制点的角Tγi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 表示m条实际正畸弓丝曲线的第i个弯制点的角误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的角误差波动度的上限值为(σa)max;
步骤四、理论正畸弓丝曲线弯制点复杂度及归一化弯制点角距比判断:
按照公式计算理论正畸弓丝曲线上第i个弯制点的复杂度,即1Cr表示实际正畸弓丝曲线上第1个弯制点Tp1的复杂度,按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比,即表示实际正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,i的初始值为i=1,判断iCr>Crb是否成立,判断是否成立,
具体为:
a)理论正畸弓丝曲线弯制点复杂度验证;
如果iCr>Crb成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的复杂度大于复杂度分界值Crb,则该理论正畸弓丝曲线第i个弯制点的复杂度较大,即弯制难度大,因此在对实际正畸弓丝曲线第i个弯制点的波动度评价过程中,除了需要评价实际正畸弓丝曲线第i个弯制点的线误差波动度和角误差波动度,还需要评价实际正畸弓丝曲线第i个弯制点的曲率误差波动度,跳转至步骤五a);
如果iCr>Crb不成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的复杂度不大于复杂度分界值Crb,则该理论正畸弓丝曲线第i个弯制点的复杂度较小,但仍存在理论正畸弓丝曲线弯制点归一化弯制点角距比较大,而归一化弯制点密度较小这一情况;归一化弯制点角距比包含实际正畸弓丝曲线的形状特征,因此需要对理论正畸弓丝曲线归一化弯制点角距比进行验证,跳转至步骤五b);
b)理论正畸弓丝曲线归一化弯制点角距比验证;
如果成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比大于归一化弯制点角距比分界值(TE*)b,则该理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比较大,该点处的正畸弓丝曲线形状较为复杂,因此在对实际正畸弓丝曲线第i个弯制点的波动度评价过程中,同样需要对实际正畸弓丝曲线第i个弯制点的曲率误差波动度进行评价,跳转至步骤五a);
如果不成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比不大于归一化弯制点角距比分界值(TE*)b,则理论正畸弓丝曲线第i个弯制点的复杂度和归一化弯制点角距比均较小,即该弯制点的弯制难度较小且在该点处的正畸弓丝曲线形状较为简单,因此仅需要评价该弯制点的线误差波动度和角误差波动度,跳转至步骤五b);
步骤五、实际正畸弓丝曲线误差波动度评价:
根据iσd、iσa和iσc的计算公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度、角误差波动度和曲率误差波动度,i的初始值为i=1;
a)实际正畸弓丝曲线曲率误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度,判断iσc≤(σc)max是否成立,
具体为:
如果iσc≤(σc)max成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度在允许范围内,则跳转至步骤五b);
如果iσc≤(σc)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出许可范围,正畸弓丝评价结束;
b)实际正畸弓丝曲线线误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度,判断iσd≤(σd)max是否成立,
具体为:
如果iσd≤(σd)max成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度在允许范围内,则跳转至步骤五c);
如果iσd≤(σd)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的线误差波动度超出许可范围,正畸弓丝评价结束;
c)实际正畸弓丝曲线角误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的角误差波动度,判断iσa≤(σa)max是否成立,
具体为:
如果iσa≤(σa)max成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度在允许范围内,则跳转至步骤六;
如果iσa≤(σa)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的角误差波动度超出许可范围,正畸弓丝评价结束;
步骤六、判断m条实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对m条实际正畸弓丝曲线上所有弯制点进行评价,则令i=i+1,即表示对下一组实际正畸弓丝曲线弯制点进行评价,跳转至步骤四;
如果i=n成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行评价,且所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,则输出所有实际正畸弓丝曲线弯制点其相应应评价的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,正畸弓丝评价结束。
本发明的有益效果为:
1、本发明针对复杂度Cr的计算方法,提出了以归一化弯制点角距比和归一化弯制点密度作为参数的计算方法;归一化弯制点角距比值的大小能够反映正畸弓丝弯制点自身的弯制难易程度,而归一化弯制点密度值则能够反映因与相邻弯制点距离的远近而增加的额外的弯制难度;通过对弯制点角距比值以及弯制点密度值进行归一化处理,解决了因弯制点角距比和弯制点密度的量纲及数值范围的不同而不适于进行算数运算的问题,从而能够通过对归一化角距比值和归一化密度值进行算数运算,实现了对正畸弓丝弯制点的弯制难度的综合和全面的量化评价,为后续正畸弓丝复杂度判断提供了依据。
2、本发明针对正畸弓丝评价方法,提出了以弯制点复杂度为理论正畸弓丝曲线弯制点特征的判断参数,设定理论正畸弓丝曲线弯制点复杂度的分界值为Crb,当理论正畸弓丝曲线弯制点的复杂度小于所设定的分界值Crb时,根据复杂度计算可能出现的归一化弯制点角距比较大而归一化弯制点密度较小这种情况,又提出以弯制点归一化弯制点角距比作为辅助判断参数,设定正畸弓丝曲线弯制点归一化弯制点角距比的分界值为(TE*)b,通过对理论正畸弓丝曲线弯制点的复杂度和归一化弯制点角距比的计算,完成对理论正畸弓丝曲线弯制点处的曲线复杂程度及弯制困难程度的量化表示,进而为是否需要对实际正畸弓丝曲线弯制点进行曲率误差波动度的评价提供了判断依据。
3、本发明通过计算m条实际正畸弓丝曲线上的第i个弯制点的线误差的方差,实现了对实际正畸弓丝曲线第i个弯制点的线误差波动度iσd的计算,设定实际正畸弓丝曲线弯制点线误差波动度的上限值为(σd)max;通过计算m条实际正畸弓丝曲线上的第i个弯制点的平均偏置误差的方差,实现了对实际正畸弓丝曲线第i个弯制点的角误差波动度iσa的计算,设定实际正畸弓丝曲线弯制点角误差波动度的上限值为(σa)max;通过对实际正畸弓丝曲线弯制点的线误差波动度和角误差波动度的计算,实现了对实际正畸弓丝曲线位置弯制稳定性的量化评价。
4、本发明针对需要进行形状弯制稳定性评价的实际正畸弓丝曲线弯制点,提出了实际正畸弓丝曲线弯制点曲率误差波动度σc的概念,并设定加权曲率误差的上限值为(σc)max;实际正畸弓丝曲线弯制点处的曲率能够反映实际正畸弓丝弯制点处的弯制形状,针对需要进行形状弯制稳定性评价的实际正畸弓丝曲线弯制点,通过计算m条实际正畸弓丝曲线上的第i个弯制点的曲率误差的方差,实现了对实际正畸弓丝曲线形状弯制稳定性的量化评价。
5、同本发明人同日申报的发明专利《一种基于曲率误差预先判断的正畸弓丝弯制点误差波动度评价方法》相比,本发明所提供的方法不要求实际正畸弓丝曲线的弯制点预先满足复杂度,归一化弯制点角距比的约束,另外同本发明人同日申报的发明专利《一种基于位置误差判断的正畸弓丝弯制点误差波动度评价方法》相比,本发明所提供的方法不要求理论正畸弓丝曲线弯制点的复杂度和归一化弯制点角距比均小于其相应的分界值这一约束条件;与两种方法相比,本方法对实际正畸弓丝曲线弯制点进行判断分类,对不同类别的实际正畸弓丝曲线弯制点应用不同的实际正畸弓丝曲线弯制点误差波动度评价,从而使得该方法对所有能够进行稳定性评价的正畸弓丝曲线都适用,因此在正畸弓丝误差波动度评价中具有普遍适用性和全面性。
6、本发明针对理论正畸弓丝曲线弯制点复杂度Cr的计算方法,提出了以归一化弯制点角距比和归一化弯制点密度作为参数的计算方法;归一化弯制点角距比值的大小能够反映正畸弓丝曲线弯制点自身的弯制难易程度,而归一化弯制点密度值则能够反映因与相邻弯制点距离的远近而增加的额外的弯制难度;通过对弯制点角距比值以及弯制点密度值进行归一化处理,解决了因弯制点角距比和弯制点密度的量纲及数值范围的不同而不适于进行算数运算的问题,从而能够通过对归一化角距比值和归一化密度值进行算数运算,通过通过判断正畸弓丝曲线弯制点的复杂度及归一化弯制点角距比即可完成对正畸弓丝弯制点在形状及位置的弯制难度的判断,进而实现对正畸弓丝弯制点的弯制难度的综合和全面的量化评价,为后续对弯制点采用什么类型的误差波动度评价提供依据。
附图说明
为了易于说明,本发明由下述的具体实施及附图作以详细描述。
图1为一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法流程图;
图2为理论正畸弓丝曲线在进行空间变换前的位置示意图;
图3为为m条实际正畸弓丝曲线在空间变换前的位姿示意图;
图4为第12个实际正畸弓丝曲线弯制点线误差波动度大于上限值的m条处于最终位姿下的实际正畸弓丝曲线示意图;
图5为m条处于最终位姿下的弯制点误差波动度均在允许范围内的实际正畸弓丝曲线示意图;
具体实施方式
为使本发明专利的目的、技术方案和优点更加清楚明了,下面通过附图中示出的具体实施例来描述本发明专利,但是应该理解,这些描述只是示例性的,而并非要限制本发明专利的范围,此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明专利的概念。
实施实例1:如图1、图2、图3、图4、图5所示,本具体实施方式采用以下技术方案:一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法,所述方法的具体实现过程为:
步骤一、理论正畸弓丝曲线数据及实际正畸弓丝曲线数据导入:
以右手定则建立o-xyz三维正畸弓丝误差标定坐标系w,以正畸医师根据患者牙列形态设计的具有n个弯制点的理论正畸弓丝曲线,计算并输入理论正畸弓丝曲线弯制点信息集P’T={Tp’1,Tp'2,Tp'3,...,Tp’i,...,Tp'n},Tp’i=(Tα’i,Tβ’i,Tγ’i,Td’i)为理论正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,i的取值范围为1≤i≤n,Tα’i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Tβ’i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Tγ’i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Td’i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;理论正畸弓丝曲线左端点为ps,理论正畸弓丝曲线右端点为pf,ps和pf之间连线的中点为To',对理论正畸弓丝曲线进行空间变换:令点To'与三维正畸弓丝误差标定坐标系w的原点o重合,理论正畸弓丝曲线左端点ps位于y轴负半轴,理论正畸弓丝曲线右端点pf位于y轴正半轴,且理论正畸弓丝曲线与x轴正半轴无交点,再令理论正畸弓丝曲线沿y轴正方向顺时针旋转,直至理论正畸弓丝曲线与x轴正半轴出现交点,将理论正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系w中的最终位姿,计算并输入处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpi,...,Tpn},Tpi=(Tαi,Tβi,Tγi,Tdi)为处于最终位姿下的理论正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Tαi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Tβi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Tγi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Tdi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
以根据理论正畸弓丝曲线弯制出的m条具有n个弯制点的实际正畸弓丝曲线,计算并输入实际正畸弓丝曲线信息集为第j条实际正畸弓丝曲线弯制点信息集,j的取值范围为1≤j≤m,为第j条实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;第j条实际正畸弓丝曲线左端点为jp's,第j条实际正畸弓丝曲线右端点为jp'f,jp's和jp'f之间连线的中点为对第j条实际正畸弓丝曲线进行空间变换:令点与三维正畸弓丝误差标定坐标系w的原点o重合,第j条实际弓丝曲线左端点jp's位于y轴负半轴,第j条实际正畸弓丝曲线右端点jp'f位于y轴正半轴,且第j条实际正畸弓丝曲线与x轴正半轴无交点,再令第j条实际正畸弓丝曲线沿y轴正方向顺时针旋转,直至第j条实际正畸弓丝曲线与x轴正半轴出现交点,将第j条实际正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系中的最终位姿,计算并输入m条处于最终位姿下的实际正畸弓丝曲线信息集为第j条处于最终位姿下的实际正畸弓丝曲线弯制点信息集,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
步骤二、实际正畸弓丝曲线弯制点复杂度的设定:
定义理论正畸弓丝曲线弯制点复杂度,用符号Cr表示,Cr是对理论正畸弓丝曲线弯制点的弯制难易程度的综合量化描述,规定理论正畸弓丝曲线第i个弯制点的复杂度表示为 表示理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比,规定 TEi表示理论正畸弓丝曲线第i个弯制点的弯制点角距比,弯制点角距比是对正畸弓丝曲线上单个弯制点的弯制复杂程度的量化描述,规定 Tθi为作用在理论正畸弓丝曲线弯制点Tpi处的弯制角度,表示作用在理论正畸弓丝曲线第i个弯制点处的弯制距离,即理论正畸弓丝曲线弯制点Tpi-1与Tpi之间的曲线段的长度,对于理论正畸弓丝曲线第1个弯制点Tp1,表示弯制点Tp1到理论正畸弓丝曲线左端点ps之间的曲线段长度,TEmin为理论正畸弓丝曲线弯制点角距比的最小值,TEmax为理论正畸弓丝曲线弯制点角距比的最大值;规定的分界值为(TE*)b;表示理论正畸弓丝曲线第i个弯制点的归一化弯制点密度,规定 Tρi表示理论正畸弓丝曲线第i个弯制点的弯制点密度,弯制点密度是对正畸弓丝曲线上单个弯制点与相邻弯制点间紧密程度的量化描述,规定公式中的数值1表示为1个弯制点,Tli表示理论正畸弓丝曲线第i个弯制点与其距离最近的弯制点之间的直线距离,即 表示理论正畸弓丝曲线第i-1个弯制点与理论正畸弓丝曲线第i个弯制点之间的直线距离,表示理论正畸弓丝曲线第i个弯制点与理论正畸弓丝曲线第i+1个弯制点之间的直线距离,当i=1时,规定 表示理论正畸弓丝曲线第1个弯制点与理论正畸弓丝曲线左端点ps之间的直线距离,表示理论正畸弓丝曲线第1个弯制点与理论正畸弓丝曲线第2个弯制点之间的直线距离,当i=n时,规定 表示理论正畸弓丝曲线第n-1个弯制点与理论正畸弓丝曲线第n个弯制点之间的直线距离,表示理论正畸弓丝曲线第n个弯制点与理论正畸弓丝曲线右端点pf之间的直线距离,Tρmin为理论正畸弓丝曲线弯制点密度的最小值,Tρmax为理论正畸弓丝曲线弯制点密度的最大值;复杂度公式中的数值2表示在计算复杂度时考虑了归一化弯制点角距比和归一化弯制单密度两个参数;规定理论正畸弓丝曲线弯制点的复杂度Cr的分界值为Crb;
步骤三、实际正畸弓丝曲线弯制点曲率误差波动度、线误差波动度、角误差波动度的设定:
定义实际正畸弓丝曲线弯制点的曲率误差波动度,用符号σc表示,σc是对实际正畸弓丝曲线弯制点弯制曲率稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的曲率误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的曲率误差率,规定 TKi表示理论正畸弓丝曲线第i个弯制点的曲率,表示第j条实际正畸弓丝曲线第i个弯制点的曲率,表示m条实际正畸弓丝曲线的第i个弯制点的曲率误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度的上限值为(σc)max;定义实际正畸弓丝曲线弯制点的线误差波动度,用符号σd表示,σd是对实际正畸弓丝曲线弯制点弯制距离稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的线误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的线误差率,规定 表示m条实际正畸弓丝曲线的第i个弯制点的线误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的线误差波动度的上限值为(σd)max;定义实际正畸弓丝曲线弯制点的角误差波动度,用符号σa表示,σa是对实际正畸弓丝曲线弯制角度稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的角误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的角误差率,规定其中为理论正畸弓丝曲线第i个弯制点的角Tαi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 为理论正畸弓丝曲线第i个弯制点的角Tβi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 为理论正畸弓丝曲线第i个弯制点的角Tγi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 表示m条实际正畸弓丝曲线的第i个弯制点的角误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的角误差波动度的上限值为(σa)max;
步骤四、理论正畸弓丝曲线弯制点复杂度及归一化弯制点角距比判断:
按照公式计算理论正畸弓丝曲线上第i个弯制点的复杂度,即1Cr表示实际正畸弓丝曲线上第1个弯制点Tp1的复杂度,按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比,即表示实际正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,i的初始值为i=1,判断iCr>Crb是否成立,判断是否成立,
具体为:
a)理论正畸弓丝曲线弯制点复杂度验证;
如果iCr>Crb成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的复杂度大于复杂度分界值Crb,则跳转至步骤五a);
如果iCr>Crb不成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的复杂度不大于复杂度分界值Crb,则对理论正畸弓丝曲线归一化弯制点角距比进行验证,跳转至步骤四b);
b)理论正畸弓丝曲线归一化弯制点角距比验证;
如果成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比大于归一化弯制点角距比分界值(TE*)b,跳转至步骤五a);
如果不成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比不大于归一化弯制点角距比分界值(TE*)b,则跳转至步骤五b);
步骤五、实际正畸弓丝曲线弯制点误差波动度评价:
根据iσd、iσa和iσc的计算公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度、角误差波动度和曲率误差波动度,i的初始值为i=1;
a)实际正畸弓丝曲线弯制点曲率误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度,判断iσc≤(σc)max是否成立,
具体为:
如果iσc≤(σc)max成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度在允许范围内,则跳转至步骤五b);
如果iσc≤(σc)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,正畸弓丝评价结束;
b)实际正畸弓丝曲线弯制点线误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度,判断iσd≤(σd)max是否成立,
具体为:
如果iσd≤(σd)max成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度在允许范围内,则跳转至步骤五c);
如果iσd≤(σd)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,正畸弓丝评价结束;
c)实际正畸弓丝曲线弯制点角误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的角误差波动度,判断iσa≤(σa)max是否成立,
具体为:
如果iσa≤(σa)max成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度在允许范围内,则跳转至步骤六;
如果iσa≤(σa)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的角误差波动度超出许可范围,正畸弓丝评价结束;
步骤六、判断m条实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对m条实际正畸弓丝曲线上所有弯制点进行评价,则令i=i+1,即表示对下一组实际正畸弓丝曲线弯制点进行评价,跳转至步骤四;
如果i=n成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行评价,且所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,则输出所有实际正畸弓丝曲线弯制点其相应应评价的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,正畸弓丝评价结束。
实施实例2:如图4所示,在m=20条包含n=16弯制点的实际正畸弓丝曲线上进行基于弯制点复杂度判断的正畸弓丝误差波动度评价过程中,i的初始值为i=1,设定弯制点复杂度的分界值为Crb=0.5,归一化弯制点角距比的分界值(TE*)b=0.4,曲率误差波动度的上限值为(σc)max=0.07,线误差波动度的上限值为(σd)max=0.05,角误差波动度的上限值为(σa)max=0.06;根据步骤四中计算理论正畸弓丝曲线第1个弯制点复杂度1Cr=0.53,通过比较得出1Cr>Crb,则跳转至步骤五a)进行实际正畸弓丝曲线曲率误差波动度、线误差波动度和角误差波动度评价,根据iσc、iσd、iσa的计算公式,得出20条实际正畸弓丝曲线的第1个弯制点的曲率误差波动度1σc=0.05<(σc)max,再判断20条实际正畸弓丝曲线的第1个弯制点的线误差波动度1σd=0.043<(σd)max,最后判断20条实际正畸弓丝曲线的第1个弯制点的角误差波动度1σa=0.04<(σa)max,跳转至步骤六,根据判断此时1=16不成立,则令i=i+1,跳转至步骤四进行理论正畸弓丝曲线第2个弯制点的复杂度以及弯制点角距比的判断,根据相应的计算公式得出2Cr=0.37,通过比较得出2Cr<Crb,则跳转至步骤四b),判断得出则跳转至步骤五b)进行实际正畸弓丝曲线线误差波动度和角误差波动度评价,根据iσd、iσa的计算公式,得出20条实际正畸弓丝曲线的第2个弯制点的线误差波动度2σd=0.044<(σd)max,再判断20条实际正畸弓丝曲线的第2个弯制点的角误差波动度2σa=0.047<(σa)max,跳转至步骤六,根据判断此时2=16不成立,则令i=i+1,跳转至步骤四进行理论正畸弓丝曲线第3个弯制点的复杂度以及弯制点角距比的判断,根据此步骤,计算得出20条实际正畸弓丝曲线的第3、第4、第5、第6、第7、第8、第9、第10、第11个弯制点其相应应评价的曲率误差波动度、线误差波动度和角误差波动度均小于相应的上限值,在计算20条实际正畸弓丝曲线的第12个弯制点时,得出12Cr=0.48<Crb,则跳转至步骤五b)进行实际正畸弓丝曲线线误差波动度和角误差波动度评价,根据iσd、iσa的计算公式,得出20条实际正畸弓丝曲线的第12个弯制点的线误差波动度12σd=0.063>(σd)max,则输出实际正畸弓丝曲线第12个弯制点的线误差波动度超出允许范围,不再对后续实际正畸弓丝曲线弯制点进行评价,正畸弓丝评价结束。
实施实例3:如图5所示,在m=20条包含n=16弯制点的实际正畸弓丝曲线上进行基于弯制点复杂度判断的正畸弓丝误差波动度评价过程中,i的初始值为i=1,设定弯制点复杂度的分界值为Crb=0.5,归一化弯制点角距比的分界值(TE*)b=0.4,曲率误差波动度的上限值为(σc)max=0.07,线误差波动度的上限值为(σd)max=0.05,角误差波动度的上限值为(σa)max=0.06;根据步骤四中计算理论正畸弓丝曲线第1个弯制点复杂度1Cr=0.61,通过比较得出1Cr>Crb,则跳转至步骤五a)进行实际正畸弓丝曲线曲率误差波动度、线误差波动度和角误差波动度评价,根据iσc、iσd、iσa的计算公式,得出20条实际正畸弓丝曲线的第1个弯制点的曲率误差波动度1σc=0.042<(σc)max,再判断20条实际正畸弓丝曲线的第1个弯制点的线误差波动度1σd=0.031<(σd)max,最后判断20条实际正畸弓丝曲线的第1个弯制点的角误差波动度1σa=0.029<(σa)max,跳转至步骤六,根据判断此时1=16不成立,则令i=i+1,跳转至步骤四进行理论正畸弓丝曲线第2个弯制点的复杂度以及弯制点角距比的判断,根据相应的计算公式得出2Cr=0.45,通过比较得出2Cr<Crb,则跳转至步骤四b),判断得出则跳转至步骤五b)进行实际正畸弓丝曲线线误差波动度和角误差波动度评价,根据iσd、iσa的计算公式,得出20条实际正畸弓丝曲线的第2个弯制点的线误差波动度2σd=0.028<(σd)max,再判断20条实际正畸弓丝曲线的第2个弯制点的角误差波动度2σa=0.011<(σa)max,跳转至步骤六,根据判断此时2=16不成立,则令i=i+1,跳转至步骤四进行理论正畸弓丝曲线第3个弯制点的复杂度以及弯制点角距比的判断,根据此步骤,计算得出20条实际正畸弓丝曲线的第3、第4、第5、第6、第7、第8、第9、第10、第11、第12、第13、第14、第15、第16个弯制点其相应应评价的曲率误差波动度、线误差波动度和角误差波动度均小于相应的上限值,跳转至步骤六,判断此时16=16成立,则输出所有实际正畸弓丝曲线弯制点其相应应评价的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,正畸弓丝评价结束。
Claims (1)
1.一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法,其特征在于:所述方法的具体实现过程为:
步骤一、理论正畸弓丝曲线数据及实际正畸弓丝曲线数据导入:
以右手定则建立o-xyz三维正畸弓丝误差标定坐标系w,以正畸医师根据患者牙列形态设计的具有n个弯制点的理论正畸弓丝曲线,计算并输入理论正畸弓丝曲线弯制点信息集P′T={Tp′1,Tp′2,Tp′3,...,Tp′i,...,Tp′n},Tp′i=(Tα′i,Tβ′i,Tγ′i,Td′i)为理论正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,i的取值范围为1≤i≤n,Tα′i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Tβ′i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Tγ′i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Td′i为理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;理论正畸弓丝曲线左端点为ps,理论正畸弓丝曲线右端点为pf,ps和pf之间连线的中点为To',对理论正畸弓丝曲线进行空间变换:令点To'与三维正畸弓丝误差标定坐标系w的原点o重合,理论正畸弓丝曲线左端点ps位于y轴负半轴,理论正畸弓丝曲线右端点pf位于y轴正半轴,且理论正畸弓丝曲线与x轴正半轴无交点,再令理论正畸弓丝曲线沿y轴正方向顺时针旋转,直至理论正畸弓丝曲线与x轴正半轴出现交点,将理论正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系w中的最终位姿,计算并输入处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpi,...,Tpn},Tpi=(Tαi,Tβi,Tγi,Tdi)为处于最终位姿下的理论正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Tαi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Tβi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Tγi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Tdi为处于最终位姿下的理论正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
以根据理论正畸弓丝曲线弯制出的m条具有n个弯制点的实际正畸弓丝曲线,计算并输入实际正畸弓丝曲线信息集为第j条实际正畸弓丝曲线弯制点信息集,j的取值范围为1≤j≤m,为第j条实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,为第j条实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;第j条实际正畸弓丝曲线左端点为jp′s,第j条实际正畸弓丝曲线右端点为jp′f,jp′s和jp′f之间连线的中点为对第j条实际正畸弓丝曲线进行空间变换:令点与三维正畸弓丝误差标定坐标系w的原点o重合,第j条实际弓丝曲线左端点jp′s位于y轴负半轴,第j条实际正畸弓丝曲线右端点jp′f位于y轴正半轴,且第j条实际正畸弓丝曲线与x轴正半轴无交点,再令第j条实际正畸弓丝曲线沿y轴正方向顺时针旋转,直至第j条实际正畸弓丝曲线与x轴正半轴出现交点,将第j条实际正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系中的最终位姿,计算并输入m条处于最终位姿下的实际正畸弓丝曲线信息集为第j条处于最终位姿下的实际正畸弓丝曲线弯制点信息集,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,为第j条处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
步骤二、实际正畸弓丝曲线弯制点复杂度的设定:
定义理论正畸弓丝曲线弯制点复杂度,用符号Cr表示,Cr是对理论正畸弓丝曲线弯制点的弯制难易程度的综合量化描述,规定理论正畸弓丝曲线第i个弯制点的复杂度表示为 表示理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比,规定 TEi表示理论正畸弓丝曲线第i个弯制点的弯制点角距比,弯制点角距比是对正畸弓丝曲线上单个弯制点的弯制复杂程度的量化描述,规定 Tθi为作用在理论正畸弓丝曲线弯制点Tpi处的弯制角度,表示作用在理论正畸弓丝曲线第i个弯制点处的弯制距离,即理论正畸弓丝曲线弯制点Tpi-1与Tpi之间的曲线段的长度,对于理论正畸弓丝曲线第1个弯制点Tp1,表示弯制点Tp1到理论正畸弓丝曲线左端点ps之间的曲线段长度,TEmin为理论正畸弓丝曲线弯制点角距比的最小值,TEmax为理论正畸弓丝曲线弯制点角距比的最大值;规定的分界值为(TE*)b;表示理论正畸弓丝曲线第i个弯制点的归一化弯制点密度,规定 Tρi表示理论正畸弓丝曲线第i个弯制点的弯制点密度,弯制点密度是对正畸弓丝曲线上单个弯制点与相邻弯制点间紧密程度的量化描述,规定公式中的数值1表示为1个弯制点,Tli表示理论正畸弓丝曲线第i个弯制点与其距离最近的弯制点之间的直线距离,即 表示理论正畸弓丝曲线第i-1个弯制点与理论正畸弓丝曲线第i个弯制点之间的直线距离,表示理论正畸弓丝曲线第i个弯制点与理论正畸弓丝曲线第i+1个弯制点之间的直线距离,当i=1时,规定 表示理论正畸弓丝曲线第1个弯制点与理论正畸弓丝曲线左端点ps之间的直线距离,表示理论正畸弓丝曲线第1个弯制点与理论正畸弓丝曲线第2个弯制点之间的直线距离,当i=n时,规定 表示理论正畸弓丝曲线第n-1个弯制点与理论正畸弓丝曲线第n个弯制点之间的直线距离,表示理论正畸弓丝曲线第n个弯制点与理论正畸弓丝曲线右端点pf之间的直线距离,Tρmin为理论正畸弓丝曲线弯制点密度的最小值,Tρmax为理论正畸弓丝曲线弯制点密度的最大值;复杂度公式中的数值2表示在计算复杂度时考虑了归一化弯制点角距比和归一化弯制单密度两个参数;规定理论正畸弓丝曲线弯制点的复杂度Cr的分界值为Crb;
步骤三、实际正畸弓丝曲线弯制点曲率误差波动度、线误差波动度、角误差波动度的设定:
定义实际正畸弓丝曲线弯制点的曲率误差波动度,用符号σc表示,σc是对实际正畸弓丝曲线弯制点弯制曲率稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的曲率误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的曲率误差率,规定 TKi表示理论正畸弓丝曲线第i个弯制点的曲率,表示第j条实际正畸弓丝曲线第i个弯制点的曲率,表示m条实际正畸弓丝曲线的第i个弯制点的曲率误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度的上限值为(σc)max;定义实际正畸弓丝曲线弯制点的线误差波动度,用符号σd表示,σd是对实际正畸弓丝曲线弯制点弯制距离稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的线误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的线误差率,规定 表示m条实际正畸弓丝曲线的第i个弯制点的线误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的线误差波动度的上限值为(σd)max;定义实际正畸弓丝曲线弯制点的角误差波动度,用符号σa表示,σa是对实际正畸弓丝曲线弯制角度稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的角误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的角误差率,规定其中为理论正畸弓丝曲线第i个弯制点的角Tαi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 为理论正畸弓丝曲线第i个弯制点的角Tβi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 为理论正畸弓丝曲线第i个弯制点的角Tγi与第j条实际正畸弓丝曲线第i个弯制点的角之间的误差率,规定 表示m条实际正畸弓丝曲线的第i个弯制点的角误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的角误差波动度的上限值为(σa)max;
步骤四、理论正畸弓丝曲线弯制点复杂度及归一化弯制点角距比判断:
按照公式计算理论正畸弓丝曲线上第i个弯制点的复杂度,即1Cr表示实际正畸弓丝曲线上第1个弯制点Tp1的复杂度,按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比,即表示实际正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,i的初始值为i=1,判断iCr>Crb是否成立,判断是否成立,
具体为:
a)理论正畸弓丝曲线弯制点复杂度验证;
如果iCr>Crb成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的复杂度大于复杂度分界值Crb,则跳转至步骤五a);
如果iCr>Crb不成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的复杂度不大于复杂度分界值Crb,则对理论正畸弓丝曲线归一化弯制点角距比进行验证,跳转至步骤四b);
b)理论正畸弓丝曲线归一化弯制点角距比验证;
如果成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比大于归一化弯制点角距比分界值(TE*)b,跳转至步骤五a);
如果不成立,说明处于最终位姿下的理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比不大于归一化弯制点角距比分界值(TE*)b,则跳转至步骤五b);
步骤五、实际正畸弓丝曲线弯制点误差波动度评价:
根据iσd、iσa和iσc的计算公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度、角误差波动度和曲率误差波动度,i的初始值为i=1;
a)实际正畸弓丝曲线弯制点曲率误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度,判断iσc≤(σc)max是否成立,
具体为:
如果iσc≤(σc)max成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度在允许范围内,则跳转至步骤五b);
如果iσc≤(σc)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,正畸弓丝评价结束;
b)实际正畸弓丝曲线弯制点线误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度,判断iσd≤(σd)max是否成立,
具体为:
如果iσd≤(σd)max成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度在允许范围内,则跳转至步骤五c);
如果iσd≤(σd)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,正畸弓丝评价结束;
c)实际正畸弓丝曲线弯制点角误差波动度评价;
根据公式计算m条实际正畸弓丝曲线第i个弯制点的角误差波动度,判断iσa≤(σa)max是否成立,
具体为:
如果iσa≤(σa)max成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度在允许范围内,则跳转至步骤六;
如果iσa≤(σa)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的角误差波动度超出许可范围,正畸弓丝评价结束;
步骤六、判断m条实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对m条实际正畸弓丝曲线上所有弯制点进行评价,则令i=i+1,即表示对下一组实际正畸弓丝曲线弯制点进行评价,跳转至步骤四;
如果i=n成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行评价,且所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,则输出所有实际正畸弓丝曲线弯制点其相应评价的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,正畸弓丝评价结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210402124.7A CN115024838B (zh) | 2022-04-18 | 2022-04-18 | 一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210402124.7A CN115024838B (zh) | 2022-04-18 | 2022-04-18 | 一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115024838A CN115024838A (zh) | 2022-09-09 |
CN115024838B true CN115024838B (zh) | 2023-09-05 |
Family
ID=83119257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210402124.7A Active CN115024838B (zh) | 2022-04-18 | 2022-04-18 | 一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115024838B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160082842A (ko) * | 2014-12-29 | 2016-07-11 | 한라대학교산학협력단 | 컴퓨터를 이용한 치열교정용 아치와이어의 굽힘점 및 굽힘각도 산출 방법 및 그 프로그램 |
CN111588499A (zh) * | 2020-06-01 | 2020-08-28 | 哈尔滨理工大学 | 一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法 |
CN111588502A (zh) * | 2020-06-01 | 2020-08-28 | 哈尔滨理工大学 | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 |
-
2022
- 2022-04-18 CN CN202210402124.7A patent/CN115024838B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20160082842A (ko) * | 2014-12-29 | 2016-07-11 | 한라대학교산학협력단 | 컴퓨터를 이용한 치열교정용 아치와이어의 굽힘점 및 굽힘각도 산출 방법 및 그 프로그램 |
CN111588499A (zh) * | 2020-06-01 | 2020-08-28 | 哈尔滨理工大学 | 一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法 |
CN111588502A (zh) * | 2020-06-01 | 2020-08-28 | 哈尔滨理工大学 | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 |
Non-Patent Citations (1)
Title |
---|
Modeling and Experimentation of the Unidirectional Orthodontic Force of Second Sequential Loop Orthodontic Archwire;Jin-Gang Jiang等;Applied bionics and biomechanics;第1-11页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115024838A (zh) | 2022-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114972184B (zh) | 一种基于权值比例法的正畸弓丝误差评价方法 | |
CN114943685B (zh) | 一种基于轮廓降维法的正畸弓丝误差评价方法 | |
CN116644558B (zh) | 一种基于误差评定域的正畸弓丝误差评价方法 | |
CN114983593B (zh) | 一种基于向量共线的正畸弓丝弯制点误差评价方法 | |
EP2056734B1 (en) | System and method for modeling of interproximal reduction of teeth | |
CN115024840B (zh) | 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 | |
CN115040275B (zh) | 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法 | |
CN111588502B (zh) | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 | |
CN114972183B (zh) | 一种基于残差平方和区间划分的正畸弓丝误差评价方法 | |
CN114943058B (zh) | 一种基于位置误差判断的正畸弓丝误差波动度评价方法 | |
CN111588499A (zh) | 一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法 | |
CN115024838B (zh) | 一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 | |
CN115024839B (zh) | 一种基于归一化弯制点角距比的正畸弓丝误差率评价方法 | |
CN117216967B (zh) | 一种基于矩阵范数的正畸弓丝误差评价方法 | |
CN115024841B (zh) | 一种基于复杂度区分法的正畸弓丝误差评价方法 | |
CN115035196A (zh) | 一种基于弯制点复杂度判断的正畸弓丝误差率评价方法 | |
CN111588500B (zh) | 一种正畸弓丝弯制顺序规划等角度划分角度确定方法 | |
CN114983594B (zh) | 一种基于弯制点价值度的正畸弓丝评价方法 | |
CN117408943B (zh) | 一种基于曲线拟合法的正畸弓丝误差评价方法 | |
CN115035197A (zh) | 一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 | |
CN114998190B (zh) | 一种基于点对区间的正畸弓丝评价方法 | |
CN111588497B (zh) | 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法 | |
CN117218073B (zh) | 一种基于空间等距分割面的正畸弓丝误差评价方法 | |
CN117219256A (zh) | 一种基于弯制点曲率和降维角距偏差域的正畸弓丝误差评价方法 | |
CN115018245A (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 |