CN115024840B - 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 - Google Patents
一种基于归一化弯制点密度的正畸弓丝误差率评价方法 Download PDFInfo
- Publication number
- CN115024840B CN115024840B CN202210402396.7A CN202210402396A CN115024840B CN 115024840 B CN115024840 B CN 115024840B CN 202210402396 A CN202210402396 A CN 202210402396A CN 115024840 B CN115024840 B CN 115024840B
- Authority
- CN
- China
- Prior art keywords
- orthodontic archwire
- curve
- bending point
- 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 481
- 238000011156 evaluation Methods 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 claims description 54
- 230000009191 jumping Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 8
- 230000001131 transforming effect Effects 0.000 claims description 6
- 210000004513 dentition Anatomy 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 230000036346 tooth eruption Effects 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 8
- 238000011158 quantitative evaluation Methods 0.000 abstract description 4
- 238000012937 correction Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 2
- 208000025157 Oral disease Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 208000030194 mouth disease Diseases 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 230000007704 transition Effects 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Computational Mathematics (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Epidemiology (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Dentistry (AREA)
- General Engineering & Computer Science (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (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之间连线的长度;
以根据理论正畸弓丝曲线弯制出的具有n个弯制点的实际正畸弓丝曲线,计算并输入实际正畸弓丝曲线弯制点信息集RP'={Rp′1,Rp'2,Rp'3,...,Rp′i,...,Rp'n},Rp′i=(Rα′i,Rβ'i,Rγ′i,Rd′i)为实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Rα′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Rβ′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Rγ′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Rd′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;实际正畸弓丝曲线左端点为p's,实际正畸弓丝曲线右端点为p'f,p's和p'f之间连线的中点为Ro',对实际正畸弓丝曲线进行空间变换:令点Ro'与三维正畸弓丝误差标定坐标系w的原点o重合,实际弓丝曲线左端点p's位于y轴负半轴,实际正畸弓丝曲线右端点p'f位于y轴正半轴,且实际正畸弓丝曲线与x轴正半轴无交点,再令实际正畸弓丝曲线沿y轴正方向顺时针旋转,直至实际正畸弓丝曲线与x轴正半轴出现交点,将实际正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系w中的最终位姿,计算并输入处于最终位姿下的实际正畸弓丝曲线弯制点信息集RP={Rp1,Rp2,Rp3,...,Rpi,...,Rpn},Rpi=(Rαi,Rβi,Rγi,Rdi)为处于最终位姿下的实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Rαi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Rβi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Rγi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Rdi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
步骤二、计算并判断实际正畸弓丝曲线端距误差:
定义实际正畸弓丝曲线端距误差,用符号Δ表示,由于正畸弓丝弯制过程中,各弯制点的成形误差会累计到正畸弓丝两端点之间的距离上,因此需要对端距误差Δ进行计算和判断,规定Δ=|Ta-Ra|,Ta表示理论正畸弓丝曲线左端点ps和理论正畸弓丝曲线右端点pf之间的直线距离,Ra表示实际正畸弓丝曲线左端点p's和实际正畸弓丝右端点p'f之间的距离,规定端距误差Δ的上限值为Δmax,判断Δ≤Δmax是否成立,
具体为:
如果Δ≤Δmax成立,说明该实际正畸弓丝曲线端距误差在允许范围内,则跳转至步骤三;
如果Δ≤Δmax不成立,说明该实际正畸弓丝曲线端距误差超出允许范围,则输出实际正畸弓丝曲线端距误差过大,正畸弓丝评价结束;
步骤三、理论正畸弓丝曲线弯制点复杂度计算:
定义理论正畸弓丝曲线弯制点的复杂度,用符号Cr表示,Cr是对理论正畸弓丝弯制点的弯制难易程度的综合量化描述,理论正畸弓丝曲线弯制点的Cr值越高,即该弯制点在弯制时越困难,说明该弯制点误差率的大小对于实际正畸弓丝评价的影响越大,规定理论正畸弓丝曲线第i个弯制点的复杂度表示为表示理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比,规定/> TEi表示理论正畸弓丝曲线第i个弯制点的弯制点角距比,弯制点角距比是对正畸弓丝弓丝曲线单个弯制点的弯制复杂程度的量化描述,规定/> Tθi为作用在理论正畸弓丝曲线弯制点Tpi处的弯制角度,表示作用在理论正畸弓丝曲线第i个弯制点处的弯制距离,即理论正畸弓丝曲线弯制点Tpi-1与Tpi之间的曲线段的长度,对于理论正畸弓丝曲线第1个弯制点Tp1,/>表示弯制点Tp1到理论正畸弓丝曲线左端点ps之间的曲线段长度,TEmin为理论正畸弓丝曲线弯制点角距比的最小值,TEmax为理论正畸弓丝曲线弯制点角距比的最大值,规定/>的上限值为/>表示理论正畸弓丝曲线第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为理论正畸弓丝曲线弯制点密度的最大值;A表示在计算复杂度时考虑的参数个数;
步骤四、理论正畸弓丝曲线归一化弯制点角距比验证:
按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比/>即/>表示理论正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,通过比较取出/>中的最大值/>对条件/>进行验证,
具体为:
如果成立,说明处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpn}内的最大归一化弯制点角距比不大于所设定的归一化弯制点角距比上限值/>则理论正畸弓丝曲线上每个弯制点的归一化角距比都小于等于归一化弯制点角距比上限值/>可知在所取的包含n个弯制点的理论正畸弓丝曲线上,所有弯制点都能够满足评价系统对理论正畸弓丝曲线归一化弯制点角距比的要求,即因理论正畸弓丝曲线归一化弯制点角距比而导致的弯制难度较小,所以在对依据此类具有特殊归一化弯制点角距比的理论正畸弓丝曲线而弯制出的实际正畸弓丝曲线的弯制点进行复杂度计算时,无需考虑理论正畸弓丝曲线归一化弯制点角距比的影响,则本方法仅以理论正畸弓丝曲线归一化弯制点密度为参数计算理论正畸弓丝曲线弯制点的复杂度,此时A=1,/>跳转至步骤五;
如果不成立,说明该评价方法不适用于此正畸弓丝曲线,则输出该评价方法不适用于此正畸弓丝曲线,正畸弓丝误差率评价结束;
步骤五、实际正畸弓丝曲线弯制点加权线误差率和实际正畸弓丝曲线弯制点加权平均偏置误差率的设定:
定义实际正畸弓丝曲线弯制点的加权线误差率,用符号表示,加权线误差率/>为理论正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系原点o之间直线距离和与理论正畸弓丝曲线弯制点相对应的实际正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系原点o之间直线距离的误差的量化描述,规定实际正畸弓丝曲线第i个弯制点的加权线误差率表示为/> ied表示实际正畸弓丝曲线第i个弯制点的线误差率,规定规定实际正畸弓丝弯制点加权线误差率/>的上限值为/>定义实际正畸弓丝曲线弯制点的加权平均偏置误差率,用符号/>表示,加权平均偏置误差率/>为理论正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系各坐标轴之间的夹角和与理论正畸弓丝曲线弯制点相对应的实际正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系各坐标轴之间夹角的平均误差的量化描述,规定实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率表示为/> iea表示实际正畸弓丝曲线第i个弯制点的平均偏置误差率,规定/>其中ieα为理论正畸弓丝曲线第i个弯制点的角Tαi与实际正畸弓丝曲线第i个弯制点的角Rαi之间的误差率,规定/> ieβ为理论正畸弓丝曲线第i个弯制点的角Tβi与实际正畸弓丝曲线第i个弯制点的角Rβi之间的误差率,规定 ieγ为理论正畸弓丝曲线第i个弯制点的角Tγi与实际正畸弓丝曲线第i个弯制点的角Rγi之间的误差率,规定/>规定实际正畸弓丝曲线弯制点加权平均偏置误差率/>的上限为/>
步骤六、实际正畸弓丝曲线误差率评估:
根据公式计算实际正畸弓丝曲线第i个弯制点的加权线误差率,根据公式/>计算实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率,i的初始值为i=1;
a)实际正畸弓丝曲线弯制点加权线误差率评价
根据公式计算出的实际正畸弓丝曲线第i个弯制点的线误差率,判断是否成立,
具体为:
如果成立,说明实际正畸弓丝曲线第i个弯制点的加权线误差率在允许范围内,则跳转至步骤六b);
如果不成立,说明实际正畸弓丝曲线第i个弯制点的加权线误差率超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的加权线误差率超出许可范围,正畸弓丝评价结束;
b)实际正畸弓丝弯制点加权平均偏置误差率评价
根据公式计算出的实际正畸弓丝曲线第i个弯制点的平均偏置误差率,判断/>是否成立,
具体为:
如果成立,说明实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率在允许范围内,则跳转至步骤七;
如果不成立,说明实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率超出许可范围,正畸弓丝评价结束;
步骤七、判断实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对所有实际正畸弓丝曲线弯制点进行评价,则令i=i+1,即表示对下一个实际正畸弓丝曲线弯制点进行评价,跳转至步骤六a);
如果i=n成立,说明已对所有实际正畸弓丝曲线弯制点进行评价,且所有实际正畸弓丝曲线弯制点的加权线误差率和加权平均偏置误差率均在允许范围内,则输出所有实际正畸弓丝曲线弯制点的加权线误差率和加权平均偏置误差率均在允许范围内,正畸弓丝评价结束。
本发明的有益效果为:
1、本发明针对正畸弓丝评价方法,提出将归一化弯制点角距比作为正畸弓丝评价的预先判断参数,设定正畸弓丝曲线归一化弯制点角距比的上限值为在正畸弓丝评价前预先验证正畸弓丝曲线上弯制点的最大归一化弯制点角距比不大于归一化弯制点角距比的上限值/>可得出正畸弓丝曲线上每个弯制点的归一化弯制点角距比/>都符合要求,从而为本方法提供了前提条件的约束,确定了正畸弓丝曲线弯制点的复杂度计算方法,避免了在对正畸弓丝曲线评价前需要对每个正畸弓丝曲线的弯制点进行判断从而确定正畸弓丝曲线弯制点复杂度计算方法的繁琐步骤,提高了评价效率。
2、本发明针对的是一类具有归一化弯制点角距比较小这一特殊属性的正畸弓丝曲线,而正畸弓丝曲线归一化弯制点角距比能够反映正畸弓丝的形状特征;正畸弓丝曲线的归一化弯制点角距比较小,说明作用在正畸弓丝曲线单位曲线段上的弯制角度较小,即正畸弓丝曲线的形状较为简单,因此对具有此特殊属性的正畸弓丝曲线不需要再进行额外的正畸弓丝曲线弯制形状准确性的评价;通过计算实际正畸弓丝曲线弯制点在三维正畸弓丝误差标定坐标系w中的位置与理论正畸弓丝曲线弯制点在三维正畸弓丝误差标定坐标系w中的位置的偏差,结合此类具有特殊弯制点角距比的理论正畸弓丝曲线所确定的理论正畸弓丝曲线弯制点复杂度,完成对实际正畸弓丝弯制点弯制位置准确性的量化评价,即可实现对正畸弓丝整体弯制准确性的评价。
3、同本发明人同日申报的发明专利《一种基于归一化弯制点角距比的正畸弓丝误差率评价方法》相比,虽然两种方法都适用于具有特殊属性的一类个性化正畸弓丝曲线,但对于《一种基于归一化弯制点角距比的正畸弓丝误差率评价方法》中所提及的方法所侧重的前提是每个理论正畸弓丝曲线弯制点的归一化弯制点密度都满足设定的要求,进而仅以理论正畸弓丝曲线弯制点的归一化弯制点角距比作为计算理论正畸弓丝曲线弯制点复杂度时考虑的参数,本方法所侧重的前提是每个理论正畸弓丝曲线弯制点的归一化弯制点角距比均满足设定要求,进而仅以理论正畸弓丝曲线弯制点的归一化弯制点密度/>作为计算理论正畸弓丝曲线弯制点复杂度时的参数,且本方法不需要再对正畸弓丝曲线进行额外的弯制形状评估;两种方法在进行实际正畸弓丝评价时应用的情况不同,因此本方法的提出与另种方法相互补偿,进而完善了实际正畸弓丝评价的系列方法。
4、同本发明人同日申报的发明专利《一种基于弯制点复杂度判断的正畸弓丝误差率评价方法》相比,的提出增加了在理论正畸弓丝曲线弯制点复杂度计算时的约束条件,该方法以一类个性化正畸弓丝曲线上的弯制点具有归一化弯制点角距比较小的特殊属性为前提,不再对理论正畸弓丝曲线上每个弯制点进行复杂度判断,且由于本发明针对的是具有归一化弯制点角距比较小这一特征的正畸弓丝曲线,因此本方法不需要再对实际正畸弓丝曲线弯制点进行加权曲率误差率评价,简化了评价步骤;由于在进行正畸弓丝评价前对弯制点的归一化弯制点角距比/>进行判定,因此在理论正畸弓丝曲线弯制点复杂度计算的过程中不再考虑理论正畸弓丝曲线弯制点归一化弯制点角距比,从而在理论正畸弓丝曲线弯制点复杂度的计算中仅考虑弯制点的归一化密度,不仅满足了实际正畸弓丝曲线的误差评价,而且使得正畸弓丝评价算法的复杂程度降低,提高了评价效率。
附图说明
为了易于说明,本发明由下述的具体实施及附图作以详细描述。
图1为一种基于归一化弯制点密度的正畸弓丝误差率评价方法流程图;
图2为理论正畸弓丝曲线在空间变换前的位置示意图;
图3为实际正畸弓丝曲线在空间变换前的位置示意图;
图4为处于最终位姿下的理论正畸弓丝曲线以及存在弯制点误差超出上限值的实际正畸弓丝曲线的示意图;
图5为处于最终位姿下的理论正畸弓丝曲线以及误差均在允许范围内的实际正畸弓丝曲线的示意图;
具体实施方式
为使本发明专利的目的、技术方案和优点更加清楚明了,下面通过附图中示出的具体实施例来描述本发明专利,但是应该理解,这些描述只是示例性的,而并非要限制本发明专利的范围,此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明专利的概念。
实施实例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之间连线的长度;
以根据理论正畸弓丝曲线弯制出的具有n个弯制点的实际正畸弓丝曲线,计算并输入实际正畸弓丝曲线弯制点信息集RP'={Rp′1,Rp'2,Rp'3,...,Rp′i,...,Rp'n},Rp′i=(Rα′i,Rβ′i,Rγ′i,Rd′i)为实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Rα′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Rβ′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Rγ′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Rd′i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;实际正畸弓丝曲线左端点为p's,实际正畸弓丝曲线右端点为p'f,p's和p'f之间连线的中点为Ro',对实际正畸弓丝曲线进行空间变换:令点Ro'与三维正畸弓丝误差标定坐标系w的原点o重合,实际弓丝曲线左端点p's位于y轴负半轴,实际正畸弓丝曲线右端点p'f位于y轴正半轴,且实际正畸弓丝曲线与x轴正半轴无交点,再令实际正畸弓丝曲线沿y轴正方向顺时针旋转,直至实际正畸弓丝曲线与x轴正半轴出现交点,将实际正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系w中的最终位姿,计算并输入处于最终位姿下的实际正畸弓丝曲线弯制点信息集RP={Rp1,Rp2,Rp3,...,Rpi,...,Rpn},Rpi=(Rαi,Rβi,Rγi,Rdi)为处于最终位姿下的实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Rαi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Rβi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Rγi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Rdi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
步骤二、计算并判断实际正畸弓丝曲线端距误差:
定义实际正畸弓丝曲线端距误差,用符号Δ表示,规定Δ=|Ta-Ra|,Ta表示理论正畸弓丝曲线左端点ps和理论正畸弓丝曲线右端点pf之间的直线距离,Ra表示实际正畸弓丝曲线左端点p's和实际正畸弓丝右端点p'f之间的距离,规定端距误差Δ的上限值为Δmax,判断Δ≤Δmax是否成立,
具体为:
如果Δ≤Δmax成立,说明该实际正畸弓丝曲线端距误差在允许范围内,则跳转至步骤三;
如果Δ≤Δmax不成立,说明该实际正畸弓丝曲线端距误差超出允许范围,则输出实际正畸弓丝曲线端距误差过大,正畸弓丝评价结束;
步骤三、理论正畸弓丝曲线弯制点复杂度计算:
定义理论正畸弓丝曲线弯制点的复杂度,用符号Cr表示,Cr是对理论正畸弓丝弯制点的弯制难易程度的综合量化描述,规定理论正畸弓丝曲线第i个弯制点的复杂度表示为表示理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比,规定 TEi表示理论正畸弓丝曲线第i个弯制点的弯制点角距比,弯制点角距比是对正畸弓丝曲线上单个弯制点的弯制复杂程度的量化描述,规定/> Tθi为作用在理论正畸弓丝曲线弯制点Tpi处的弯制角度,/>表示作用在理论正畸弓丝曲线第i个弯制点处的弯制距离,即理论正畸弓丝曲线弯制点Tpi-1与Tpi之间的曲线段的长度,对于理论正畸弓丝曲线第1个弯制点Tp1,/>表示弯制点Tp1到理论正畸弓丝曲线左端点ps之间的曲线段长度,TEmin为理论正畸弓丝曲线弯制点角距比的最小值,TEmax为理论正畸弓丝曲线弯制点角距比的最大值,规定/>的上限值为/>表示理论正畸弓丝曲线第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为理论正畸弓丝曲线弯制点密度的最大值;A表示在计算复杂度时考虑的参数个数;
步骤四、理论正畸弓丝曲线归一化弯制点角距比验证:
按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比/>即/>表示理论正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,通过比较取出/>中的最大值/>对条件/>进行验证,
具体为:
如果成立,说明处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpn}内的最大归一化弯制点角距比不大于所设定的归一化弯制点角距比上限值/>则理论正畸弓丝曲线上每个弯制点的归一化角距比都小于等于归一化弯制点角距比上限值/>此时A=1,/>跳转至步骤五;
如果不成立,说明该评价方法不适用于此正畸弓丝曲线,则输出该评价方法不适用于此正畸弓丝曲线,正畸弓丝误差率评价结束;
步骤五、实际正畸弓丝曲线弯制点加权线误差率和实际正畸弓丝曲线弯制点加权平均偏置误差率的设定:
定义实际正畸弓丝曲线弯制点的加权线误差率,用符号表示,加权线误差率/>为理论正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系原点o之间直线距离和与理论正畸弓丝曲线弯制点相对应的实际正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系原点o之间直线距离的误差的量化描述,规定实际正畸弓丝曲线第i个弯制点的加权线误差率表示为/> ied表示实际正畸弓丝曲线第i个弯制点的线误差率,规定规定实际正畸弓丝曲线弯制点加权线误差率/>的上限值为/>定义实际正畸弓丝曲线弯制点的加权平均偏置误差率,用符号/>表示,加权平均偏置误差率/>为理论正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系各坐标轴之间的夹角和与理论正畸弓丝曲线弯制点相对应的实际正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系各坐标轴之间夹角的平均误差的量化描述,规定实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率表示为/> iea表示实际正畸弓丝曲线第i个弯制点的平均偏置误差率,规定/>其中ieα为理论正畸弓丝曲线第i个弯制点的角Tαi与实际正畸弓丝曲线第i个弯制点的角Rαi之间的误差率,规定/> ieβ为理论正畸弓丝曲线第i个弯制点的角Tβi与实际正畸弓丝曲线第i个弯制点的角Rβi之间的误差率,规定 ieγ为理论正畸弓丝曲线第i个弯制点的角Tγi与实际正畸弓丝曲线第i个弯制点的角Rγi之间的误差率,规定/>规定实际正畸弓丝曲线弯制点加权平均偏置误差率/>的上限为/>
步骤六、实际正畸弓丝曲线误差率评估:
根据公式计算实际正畸弓丝曲线第i个弯制点的加权线误差率,根据公式/>计算实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率,i的初始值为i=1;
a)实际正畸弓丝曲线弯制点加权线误差率评价
根据公式计算出的实际正畸弓丝曲线第i个弯制点的线误差率,判断是否成立,
具体为:
如果成立,说明实际正畸弓丝曲线第i个弯制点的加权线误差率在允许范围内,则跳转至步骤六b);
如果不成立,说明实际正畸弓丝曲线第i个弯制点的加权线误差率超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的加权线误差率超出许可范围,正畸弓丝评价结束;
b)实际正畸弓丝弯制点加权平均偏置误差率评价
根据公式计算出的实际正畸弓丝曲线第i个弯制点的平均偏置误差率,判断/>是否成立,
具体为:
如果成立,说明实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率在允许范围内,则跳转至步骤七;
如果不成立,说明实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率超出许可范围,正畸弓丝评价结束;
步骤七、判断实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对所有实际正畸弓丝曲线弯制点进行评价,则令i=i+1,即表示对下一个实际正畸弓丝曲线弯制点进行评价,跳转至步骤六a);
如果i=n成立,说明已对所有实际正畸弓丝曲线弯制点进行评价,且所有实际正畸弓丝曲线弯制点的加权线误差率和加权平均偏置误差率均在允许范围内,则输出所有实际正畸弓丝曲线弯制点的加权线误差率和加权平均偏置误差率均在允许范围内,正畸弓丝评价结束。
实施实例2:如图4所示,在一条包含n=16弯制点的实际正畸弓丝曲线上进行基于归一化弯制点密度的正畸弓丝误差率评价过程中,i的初始值为i=1,设定实际正畸弓丝曲线端距误差的上限值Δmax=4,计算该实际正畸弓丝曲线端距误差Δ=2.4,通过验证该实际正畸弓丝曲线的端距误差Δ≤Δmax成立,则跳转至步骤三;设定理论正畸弓丝曲线弯制点的归一化弯制点角距比的上限值通过计算实际正畸弓丝曲线16个弯制点的归一化弯制点角距比,比较得出其中的最大值/>通过判断可得该实际正畸弓丝曲线弯制点角距比存在/>则该评价方法适用于此实际正畸弓丝曲线;设定实际正畸弓丝曲线弯制点的加权线误差率的上限/>加权平均偏置误差率的上限根据步骤三中计算理论正畸弓丝曲线弯制点复杂度值的公式计算理论正畸弓丝曲线第1个弯制点Tp1的复杂度1Cr=0.65,根据步骤五计算实际正畸弓丝曲线第1个弯制点Rp1的加权线误差率/>和加权平均偏置误差率/>并跳转至步骤六a)判断得出/>成立,跳转至步骤六b)判断得出/>成立,即实际正畸弓丝曲线第一个弯制点的加权线误差率/>和加权平均偏置误差率/>均在允许范围内,则跳转至步骤七,通过步骤七验证此时1=16不成立,说明此时还未对实际正畸弓丝曲线的所有弯制点进行评价,因此令i=i+1并跳转至步骤六a)进行际正畸弓丝曲线第2个弯制点的评价,通过重复该步骤,可完成对后续实际正畸弓丝曲线弯制点的评价,假设实际正畸弓丝曲线的弯制点Rp2,Rp3,Rp4,Rp5,Rp6,Rp7,Rp8的加权线误差率和加权平均偏置误差率均在允许范围内,但是实际正畸弓丝曲线第9个弯制点存在/>和/>即实际正畸弓丝曲线第9个弯制点的加权线误差率超出上限,则输出实际正畸弓丝曲线第9个弯制点的加权线误差率过大,不再对后续实际正畸弓丝弯制点进行评价,正畸弓丝评价结束。
实施实例3:如图5所示,在一条包含n=16弯制点的实际正畸弓丝曲线上进行基于归一化弯制点密度的正畸弓丝误差率评价过程中,i的初始值为i=1,设定实际正畸弓丝曲线端距误差的上限值Δmax=4,计算该实际正畸弓丝曲线端距误差Δ=3.2,通过验证该实际正畸弓丝曲线的端距误差Δ≤Δmax成立,设定理论正畸弓丝曲线弯制点的归一化弯制点角距比的上限值通过计算实际正畸弓丝曲线16个弯制点的归一化弯制点角距比,比较得出其中的最大值/>通过判断该实际正畸弓丝曲线弯制点角距比存在则该评价方法适用于此实际正畸弓丝曲线;设定实际正畸弓丝曲线弯制点的加权线误差率的上限/>加权平均偏置误差率的上限/>根据步骤三中计算理论正畸弓丝曲线弯制点复杂度值的公式计算理论正畸弓丝曲线第1个弯制点Tp1的复杂度1Cr=0.35,根据步骤五计算实际正畸弓丝曲线第1个弯制点Rp1的加权线误差率和加权平均偏置误差率/>并跳转至步骤六a)判断得出/>成立,跳转至步骤六b)判断得出/>成立,即实际正畸弓丝曲线第一个弯制点的加权线误差率/>和加权平均偏置误差率/>均在允许范围内,则跳转至步骤七,通过步骤七验证此时1=16不成立,说明此时还未对实际正畸弓丝曲线的所有弯制点进行评价,因此令i=i+1并跳转至步骤六a)进行际正畸弓丝曲线第2个弯制点的评价,通过重复该步骤,可完成对后续实际正畸弓丝曲线弯制点的评价,通过计算得到后续实际正畸弓丝曲线的弯制点Rp2,Rp3,Rp4,Rp5,Rp6,Rp7,Rp8,Rp9,Rp10,Rp11,Rp12,Rp13,Rp14,Rp15,Rp16的加权线误差率和加权平均偏置误差率均未超出上限值,则输出实际正畸弓丝曲线弯制点的加权线误差率和加权平均偏置误差率均在允许范围内,正畸弓丝评价结束。/>
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之间连线的长度;
以根据理论正畸弓丝曲线弯制出的具有n个弯制点的实际正畸弓丝曲线,计算并输入实际正畸弓丝曲线弯制点信息集RP'={Rp'1,Rp'2,Rp'3,...,Rp'i,...,Rp'n},Rp'i=(Rα'i,Rβ'i,Rγ'i,Rd'i)为实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Rα'i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Rβ'i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Rγ'i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Rd'i为实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;实际正畸弓丝曲线左端点为p's,实际正畸弓丝曲线右端点为p'f,p's和p'f之间连线的中点为Ro',对实际正畸弓丝曲线进行空间变换:令点Ro'与三维正畸弓丝误差标定坐标系w的原点o重合,实际弓丝曲线左端点p's位于y轴负半轴,实际正畸弓丝曲线右端点p'f位于y轴正半轴,且实际正畸弓丝曲线与x轴正半轴无交点,再令实际正畸弓丝曲线沿y轴正方向顺时针旋转,直至实际正畸弓丝曲线与x轴正半轴出现交点,将实际正畸弓丝曲线经空间变换后的位姿设定为在三维正畸弓丝误差标定坐标系w中的最终位姿,计算并输入处于最终位姿下的实际正畸弓丝曲线弯制点信息集RP={Rp1,Rp2,Rp3,...,Rpi,...,Rpn},Rpi=(Rαi,Rβi,Rγi,Rdi)为处于最终位姿下的实际正畸弓丝曲线第i个弯制点相对于三维正畸弓丝误差标定坐标系w的位姿信息,Rαi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与x轴所成的夹角,Rβi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与y轴所成的夹角,Rγi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间的连线与z轴所成的夹角,Rdi为处于最终位姿下的实际正畸弓丝曲线第i个弯制点和三维正畸弓丝误差标定坐标系w的原点o之间连线的长度;
步骤二、计算并判断实际正畸弓丝曲线端距误差:
定义实际正畸弓丝曲线端距误差,用符号Δ表示,规定Δ=|Ta-Ra|,Ta表示理论正畸弓丝曲线左端点ps和理论正畸弓丝曲线右端点pf之间的直线距离,Ra表示实际正畸弓丝曲线左端点p's和实际正畸弓丝右端点p'f之间的距离,规定端距误差Δ的上限值为Δmax,判断Δ≤Δmax是否成立,
具体为:
如果Δ≤Δmax成立,说明该实际正畸弓丝曲线端距误差在允许范围内,则跳转至步骤三;
如果Δ≤Δmax不成立,说明该实际正畸弓丝曲线端距误差超出允许范围,则输出实际正畸弓丝曲线端距误差过大,正畸弓丝误差率评价结束;
步骤三、理论正畸弓丝曲线弯制点复杂度计算:
定义理论正畸弓丝曲线弯制点的复杂度,用符号Cr表示,Cr是对理论正畸弓丝弯制点的弯制难易程度的综合量化描述,规定理论正畸弓丝曲线第i个弯制点的复杂度表示为 表示理论正畸弓丝曲线第i个弯制点的归一化弯制点角距比,规定 TEi表示理论正畸弓丝曲线第i个弯制点的弯制点角距比,弯制点角距比是对正畸弓丝曲线上单个弯制点的弯制复杂程度的量化描述,规定/> Tθi为作用在理论正畸弓丝曲线弯制点Tpi处的弯制角度,Tpi-1 Tpi表示作用在理论正畸弓丝曲线第i个弯制点处的弯制距离,即理论正畸弓丝曲线弯制点Tpi-1与Tpi之间的曲线段的长度,对于理论正畸弓丝曲线第1个弯制点Tp1,Tp0 Tp1表示弯制点Tp1到理论正畸弓丝曲线左端点ps之间的曲线段长度,TEmin为理论正畸弓丝曲线弯制点角距比的最小值,TEmax为理论正畸弓丝曲线弯制点角距比的最大值,规定/>的上限值为/> 表示理论正畸弓丝曲线第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为理论正畸弓丝曲线弯制点密度的最大值;A表示在计算复杂度时考虑的参数个数;
步骤四、理论正畸弓丝曲线归一化弯制点角距比验证:
按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比/>即/>表示理论正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,通过比较取出/>中的最大值/>对条件/>进行验证,
具体为:
如果成立,说明处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpn}内的最大归一化弯制点角距比不大于所设定的归一化弯制点角距比上限值/>则理论正畸弓丝曲线上每个弯制点的归一化角距比都小于等于归一化弯制点角距比上限值/>此时A=1,/>跳转至步骤五;
如果不成立,说明该评价方法不适用于此正畸弓丝曲线,则输出该评价方法不适用于此正畸弓丝曲线,正畸弓丝误差率评价结束;
步骤五、实际正畸弓丝曲线弯制点加权线误差率和实际正畸弓丝曲线弯制点加权平均偏置误差率的设定:
定义实际正畸弓丝曲线弯制点的加权线误差率,用符号表示,加权线误差率/>为理论正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系原点o之间直线距离和与理论正畸弓丝曲线弯制点相对应的实际正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系原点o之间直线距离的误差的量化描述,规定实际正畸弓丝曲线第i个弯制点的加权线误差率表示为/> ied表示实际正畸弓丝曲线第i个弯制点的线误差率,规定规定实际正畸弓丝曲线弯制点加权线误差率/>的上限值为/>定义实际正畸弓丝曲线弯制点的加权平均偏置误差率,用符号/>表示,加权平均偏置误差率/>为理论正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系各坐标轴之间的夹角和与理论正畸弓丝曲线弯制点相对应的实际正畸弓丝曲线弯制点与三维正畸弓丝误差标定坐标系各坐标轴之间夹角的平均误差的量化描述,规定实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率表示为/> iea表示实际正畸弓丝曲线第i个弯制点的平均偏置误差率,规定/>其中ieα为理论正畸弓丝曲线第i个弯制点的角Tαi与实际正畸弓丝曲线第i个弯制点的角Rαi之间的误差率,规定/> ieβ为理论正畸弓丝曲线第i个弯制点的角Tβi与实际正畸弓丝曲线第i个弯制点的角Rβi之间的误差率,规定ieγ为理论正畸弓丝曲线第i个弯制点的角Tγi与实际正畸弓丝曲线第i个弯制点的角Rγi之间的误差率,规定/>规定实际正畸弓丝曲线弯制点加权平均偏置误差率/>的上限为/>
步骤六、实际正畸弓丝曲线误差率评价:
根据公式计算实际正畸弓丝曲线第i个弯制点的加权线误差率,根据公式计算实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率,i的初始值为i=1;
a)实际正畸弓丝曲线弯制点加权线误差率评价
根据公式计算出的实际正畸弓丝曲线第i个弯制点的线误差率,判断是否成立,
具体为:
如果成立,说明实际正畸弓丝曲线第i个弯制点的加权线误差率在允许范围内,则跳转至步骤六b);
如果不成立,说明实际正畸弓丝曲线第i个弯制点的加权线误差率超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的加权线误差率超出许可范围,正畸弓丝误差率评价结束;
b)实际正畸弓丝曲线弯制点加权平均偏置误差率评价
根据公式计算出的实际正畸弓丝曲线第i个弯制点的平均偏置误差率,判断/>是否成立,
具体为:
如果成立,说明实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率在允许范围内,则跳转至步骤七;
如果不成立,说明实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的加权平均偏置误差率超出许可范围,正畸弓丝误差率评价结束;
步骤七、判断实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对所有实际正畸弓丝曲线弯制点进行评价,则令i=i+1,即表示对下一个实际正畸弓丝曲线弯制点进行评价,跳转至步骤六a);
如果i=n成立,说明已对所有实际正畸弓丝曲线弯制点进行评价,且所有实际正畸弓丝曲线弯制点的加权线误差率和加权平均偏置误差率均在允许范围内,则输出所有实际正畸弓丝曲线弯制点的加权线误差率和加权平均偏置误差率均在允许范围内,正畸弓丝误差率评价结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210402396.7A CN115024840B (zh) | 2022-04-18 | 2022-04-18 | 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210402396.7A CN115024840B (zh) | 2022-04-18 | 2022-04-18 | 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115024840A CN115024840A (zh) | 2022-09-09 |
CN115024840B true CN115024840B (zh) | 2023-10-24 |
Family
ID=83119645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210402396.7A Active CN115024840B (zh) | 2022-04-18 | 2022-04-18 | 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115024840B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117297806B (zh) * | 2023-04-27 | 2024-06-07 | 哈尔滨理工大学 | 一种基于阈值区间划分的正畸弓丝误差评价包络面半径确定方法 |
CN117218073B (zh) * | 2023-09-01 | 2024-05-14 | 哈尔滨理工大学 | 一种基于空间等距分割面的正畸弓丝误差评价方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003116882A (ja) * | 2001-10-12 | 2003-04-22 | Saburo Yamamoto | 歯列矯正用アーチワイヤ、ワイヤ部材およびワイヤ切断工具 |
CN101988182A (zh) * | 2010-12-08 | 2011-03-23 | 沈阳工业大学 | 耐磨蚀牙科正畸弓丝表面改性方法及所得的耐磨蚀弓丝 |
KR20160082842A (ko) * | 2014-12-29 | 2016-07-11 | 한라대학교산학협력단 | 컴퓨터를 이용한 치열교정용 아치와이어의 굽힘점 및 굽힘각도 산출 방법 및 그 프로그램 |
WO2018235146A1 (ja) * | 2017-06-19 | 2018-12-27 | 好哉 嘉数 | 歯列矯正用アーチワイヤ、歯列矯正用アーチワイヤの製造方法、及び歯列矯正用アーチワイヤの取付方法 |
CN111588502A (zh) * | 2020-06-01 | 2020-08-28 | 哈尔滨理工大学 | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 |
-
2022
- 2022-04-18 CN CN202210402396.7A patent/CN115024840B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003116882A (ja) * | 2001-10-12 | 2003-04-22 | Saburo Yamamoto | 歯列矯正用アーチワイヤ、ワイヤ部材およびワイヤ切断工具 |
CN101988182A (zh) * | 2010-12-08 | 2011-03-23 | 沈阳工业大学 | 耐磨蚀牙科正畸弓丝表面改性方法及所得的耐磨蚀弓丝 |
KR20160082842A (ko) * | 2014-12-29 | 2016-07-11 | 한라대학교산학협력단 | 컴퓨터를 이용한 치열교정용 아치와이어의 굽힘점 및 굽힘각도 산출 방법 및 그 프로그램 |
WO2018235146A1 (ja) * | 2017-06-19 | 2018-12-27 | 好哉 嘉数 | 歯列矯正用アーチワイヤ、歯列矯正用アーチワイヤの製造方法、及び歯列矯正用アーチワイヤの取付方法 |
EP3456285A1 (en) * | 2017-06-19 | 2019-03-20 | Yoshiya Kakazu | Orthodontic archwire, method for manufacturing orthodontic archwire, and method for attaching orthodontic archwire |
CN111588502A (zh) * | 2020-06-01 | 2020-08-28 | 哈尔滨理工大学 | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115024840A (zh) | 2022-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115024840B (zh) | 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 | |
CN114943685B (zh) | 一种基于轮廓降维法的正畸弓丝误差评价方法 | |
CN114972184B (zh) | 一种基于权值比例法的正畸弓丝误差评价方法 | |
CN115040275B (zh) | 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法 | |
CN114972183B (zh) | 一种基于残差平方和区间划分的正畸弓丝误差评价方法 | |
CN116644558A (zh) | 一种基于误差评定域的正畸弓丝误差评价方法 | |
CN111588502B (zh) | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 | |
CN114943058B (zh) | 一种基于位置误差判断的正畸弓丝误差波动度评价方法 | |
CN111588499A (zh) | 一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法 | |
CN115024841B (zh) | 一种基于复杂度区分法的正畸弓丝误差评价方法 | |
CN115035196A (zh) | 一种基于弯制点复杂度判断的正畸弓丝误差率评价方法 | |
CN115024839B (zh) | 一种基于归一化弯制点角距比的正畸弓丝误差率评价方法 | |
CN111588494B (zh) | 一种基于弯制点密度的正畸弓丝变角度划分方法 | |
CN115024838B (zh) | 一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 | |
CN117297806B (zh) | 一种基于阈值区间划分的正畸弓丝误差评价包络面半径确定方法 | |
CN115035197A (zh) | 一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 | |
CN117408943B (zh) | 一种基于曲线拟合法的正畸弓丝误差评价方法 | |
CN117219256A (zh) | 一种基于弯制点曲率和降维角距偏差域的正畸弓丝误差评价方法 | |
CN117218073B (zh) | 一种基于空间等距分割面的正畸弓丝误差评价方法 | |
Jiang et al. | Digital expression and interactive adjustment method of personalized orthodontic archwire for robotic bending | |
CN117217992B (zh) | 一种基于曲率波动率及步长误差判断的牙体预备轨迹插补精度评价方法 | |
CN115018245A (zh) | 一种基于弯制域节点的正畸弓丝弯制域评价方法 | |
CN111588503B (zh) | 一种基于弯制点密度的正畸弓丝变半径圆域划分方法 | |
CN115690009A (zh) | 一种基于共面等角向量的正畸弓丝评价方法 | |
CN111588495B (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 |