CN115035197A - 一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 - Google Patents
一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 Download PDFInfo
- Publication number
- CN115035197A CN115035197A CN202210415167.9A CN202210415167A CN115035197A CN 115035197 A CN115035197 A CN 115035197A CN 202210415167 A CN202210415167 A CN 202210415167A CN 115035197 A CN115035197 A CN 115035197A
- Authority
- CN
- China
- Prior art keywords
- arch wire
- orthodontic arch
- curve
- bending point
- theoretical
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computational Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Dentistry (AREA)
- Epidemiology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Dental Tools And Instruments Or Auxiliary Dental Instruments (AREA)
Abstract
本发明公开了一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法,它涉及正畸弓丝弯制评价领域,本发明针对具有弯制点复杂度较大或归一化弯制点角距比较大特征的正畸弓丝曲线,基于空间变换后的理论正畸弓丝曲线弯制点信息集和实际正畸弓丝曲线信息集,通过对理论正畸弓丝曲线弯制点复杂度和归一化弯制点角距比判断和评价标定值的设定、结合实际正畸弓丝曲线弯制点的曲率误差波动度、线误差波动度和角误差波动度,建立一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法。本发明通过前置条件约束,避免了对每个正畸弓丝曲线弯制点的复杂度或归一化弯制点角距比进行计算分类的繁琐步骤,提高了评价效率。通过计算m条实际正畸弓丝曲线弯制点的曲率误差波动度、线误差波动度和角误差波动度,实现了对具有此特征正畸弓丝弯制稳定性的量化评估。
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之间连线的长度;
步骤二、实际正畸弓丝曲线曲率误差波动度和评价标定值的设定:
定义实际正畸弓丝曲线的曲率误差波动度,用符号σc表示,σc是对正畸弓丝曲线弯制形状稳定性的量化描述,规定实际正畸弓丝曲线第i个弯制点的曲率误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的曲率误差率,规定 TKi表示理论正畸弓丝曲线第i个弯制点的曲率,表示第j条实际正畸弓丝曲线第i个弯制点的曲率,表示m条实际正畸弓丝曲线的第i个弯制点的曲率误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度的上限值为(σc)max;定义评价标定值,用符号ξ表示,评价标定值ξ是对进行的误差波动度评价是否为曲率误差波动度的量化描述,当进行的误差波动度评价是曲率误差波动度时ξ=1,当进行的误差波动度评价不是曲率误差波动度时ξ=0,ξ的初始值为ξ=0;
步骤三、理论正畸弓丝曲线弯制点复杂度的设定:
定义理论正畸弓丝曲线弯制点复杂度,用符号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;
步骤四、理论正畸弓丝曲线弯制点复杂度及归一化弯制点角距比验证:
按照公式计算理论正畸弓丝曲线上第i个弯制点的复杂度,即1Cr表示理论正畸弓丝曲线上第1个弯制点Tp1的复杂度,通过比较取出iCr中的最小值iminCr,对条件iminCr>Crb进行验证;按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比,即表示理论正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,通过比较取出中的最小值TminE*,对条件TminE*>(TE*)b进行验证,具体为:
a)理论正畸弓丝曲线弯制点复杂度验证;
如果iminCr>Crb成立,说明处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpi,...,Tpn}内的最小弯制点复杂度仍大于所设定的复杂度分界值Crb,可知在所取的包含n个弯制点的理论正畸弓丝曲线上每个弯制点的复杂度都大于弯制点复杂度的分界值Crb,则此类正畸弓丝曲线复杂度大,即弯制难度大;且由于在实际正畸弓丝曲线的弯制过程中,对于正畸弓丝曲线形状的弯制要比正畸弓丝弯制点位置的弯制更复杂,在进行正畸弓丝形状弯制的过程中更容易产生误差,而弯制点曲率误差即是反映正畸弓丝形状弯制准确性的指标,因此在评价依据此类具有复杂度较大这一特征的理论正畸弓丝曲线而弯制的实际正畸弓丝曲线的误差波动度时,预先对实际正畸弓丝曲线弯制带你的曲率误差波动度进行评价,跳转至步骤五;
如果iminCr>Crmax不成立,说明该此类理论正畸弓丝曲线弯制点的复杂度不大,但仍存在理论正畸弓丝曲线弯制点归一化弯制点角距比较大,而归一化弯制点密度较小这一情况;归一化弯制点角距比包含实际正畸弓丝曲线的形状特征,因此需要对理论正畸弓丝曲线归一化弯制点角距比进行验证,跳转至步骤四b);
b)理论正畸弓丝曲线归一化弯制点角距比验证;
如果TminE*>(TE*)b成立,说明处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpi,...,Tpn}内弯制点的最小归一化弯制点角距比仍大于所设定的归一化弯制点分界值(TE*)b,可知在所取的包含n个弯制点的理论正畸弓丝曲线上每个弯制点的归一化弯制点角距比都大于归一化弯制点的上限值(TE*)max,则此类正畸弓丝曲线的形状较为复杂,在评价依据此类具有归一化弯制点角距比较大这一特征的理论正畸弓丝曲线而弯制的实际正畸弓丝曲线的误差波动度时,同样需要预先对实际正畸弓丝曲线弯制点的曲率误差波动度进行评价,跳转至步骤五;
如果TminE*>(TE*)b不成立,说明此类理论正畸弓丝曲线弯制点的复杂度和归一化弯制点角距比均不大,该评价方法不适用于此正畸弓丝曲线,则输出该评价方法不适用于此正畸弓丝曲线,正畸弓丝误差波动度评价结束;
步骤五、实际正畸弓丝曲线曲率误差波动度评价
根据iσc的计算公式计算m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度,判断iσc≤(σc)max是否成立,
具体为:
如果iσc≤(σc)max成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度在允许范围内,则令ξ=1,跳转至步骤八;
如果iσc≤(σc)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出许可范围,正畸弓丝评价结束;
步骤六、实际正畸弓丝曲线线误差波动度、角误差波动度的设定:
定义实际正畸弓丝曲线的线误差波动度,用符号σ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σd和iσa的计算公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度和角误差波动度,i的初始值为i=1;
a)实际正畸弓丝曲线线误差波动度评价
具体为:
如果iσd≤(σd)max成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度在允许范围内,则跳转至步骤四b);
如果iσd≤(σd)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的线误差波动度超出许可范围,正畸弓丝评价结束;
b)实际正畸弓丝曲线角误差波动度评价
具体为:
如果iσa≤(σa)max成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度在允许范围内,则令ξ=0,跳转至步骤八;
如果iσa≤(σa)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的角误差波动度超出许可范围,正畸弓丝评价结束;
步骤八、判断m条实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对m条实际正畸弓丝曲线上所有弯制点进行评价,则令i=i+1,即表示对下一组实际正畸弓丝曲线弯制点进行评价,判断ξ=1是否成立,
具体为:
如果ξ=1成立,说明进行的是曲率误差波动度的评价,则跳转至步骤五;
如果ξ=1不成立,说明进行的不是曲率误差波动度的评价,则跳转至步骤七a);
如果i=n成立,则判断ξ=1是否成立,
具体为:
如果ξ=1成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行了曲率误差波动度评价,则可进行后续实际正畸弓丝曲线弯制点线误差波动度和角误差波动度评价,跳转至步骤六;
如果ξ=1不成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行了曲率误差波动度、线误差波动度和角误差波动度评价,且所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,则输出所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,正畸弓丝评价结束。
本发明的有益效果为:
1、本发明针对正畸弓丝评价方法,提出了以弯制点复杂度或以弯制点的归一化弯制点角距比为正畸弓丝曲线评价的预先判断参数,设定正畸弓丝曲线弯制点复杂度的分界值为Crb,正畸弓丝曲线弯制点的归一化弯制点角距比的分界值为(TE*)b,在正畸弓丝曲线评价前预先验证正畸弓丝曲线上弯制点的最小复杂度仍大于弯制点复杂度的分界值Crb或正畸弓丝曲线上弯制点的归一化弯制点角距比的最小值仍大于归一化弯制点角距比的分界值(TE*)b,可得出正畸弓丝曲线上弯制点的复杂度iCr或正畸弓丝曲线弯制点的归一化弯制点角距比TEi都符合要求,从而为本方法提供了前提条件的约束,确定了正畸弓丝曲线上弯制点的复杂度或归一化弯制点角距比特征,避免了在对正畸弓丝曲线评价前需要对每个正畸弓丝曲线弯制点的复杂或归一化弯制点角距比进行计算分类并应用不同评价方案的繁琐步骤,提高了评价效率。
2、本发明针对具有弯制点复杂度较大或归一化弯制点角距比较大这种特殊属性的正畸弓丝曲线,提出了正畸弓丝曲线弯制点的曲率误差波动度σc的概念,并设定σc的上限值为(σc)max;考虑到对于弯制点复杂度较大或归一化弯制点角距比较大的正畸弓丝曲线,形状复杂,在进行正畸弓丝曲线形状弯制的过程中更容易产生误差,误差的波动度也更大,而弯制点曲率误差即是反映正畸弓丝形状弯制准确性的指标,因此本发明首先对正畸弓丝曲线弯制点的曲率误差波动度进行评价,完成对正畸弓丝曲线弯制点形状弯制稳定性的量化评估;当正畸弓丝曲线弯制点的曲率误差波动度超出许可范围时,则无需再对后续的正畸弓丝曲线弯制点线误差波动度和正畸弓丝曲线弯制点角误差波动度进行评价,简化了评价步骤。
3、本发明针对进行的误差波动度评价是否为曲率误差波动度的判断方法,提出了评价标定值ξ的概念,实现了对进行的误差波动度评价是否为曲率误差波动度的量化描述;由于在正畸弓丝曲线弯制点误差评价过程中,预先对弯制点的曲率误差波动度进行评价,而后在对弯制点的线误差波动度和弯制点的角误差波动度进行评价,在对正畸弓丝曲线弯制点的循环判断中,无法判断进行误差波动度评价是否为曲率误差波动度,进而无法确定需要跳转的步骤;而通过定义评价标定值ξ,则能在误差波动度的评价中,标定进行误差波动度评价的类型,进而为后续的循环跳转步骤的确定提供了判断依据,简化了评价步骤。
4、同本发明人同日申报的发明专利《一种基于位置误差判断的正畸弓丝误差波动度评价方法》相比,虽然两种方法都适用于具有特殊属性的一类个性化实际正畸弓丝曲线,但对于《一种基于位置误差判断的正畸弓丝误差波动度评价方法》中所提及的方法所侧重的前提是每个实际正畸弓丝曲线弯制点的复杂度和归一化弯制点角距比均小于相应设定的分界值,进而仅验证实际正畸弓丝曲线弯制点的线误差波动度和角误差波动度,本方法所侧重的前提是每个实际正畸弓丝曲线弯制点的复杂度或归一化弯制点角距比满足设定要求,进而预先判断实际正畸弓丝曲线弯制点的曲率误差波动度,在实际正畸弓丝曲线所有弯制点的曲率误差波动度均满足要求后,再进行实际正畸弓丝曲线弯制点的线误差波动度和角误差波动度的评价;两种方法在进行实际正畸弓丝曲线评价时应用的情况不同,因此本方法的提出与另种方法相互补偿,进而完善了实际正畸弓丝曲线评价的系列方法。
5、同本发明人同日申报的发明专利《一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法》相比,该方法以一类个性化正畸弓丝曲线上的弯制点具有复杂度较大或归一化弯制点角距比较大的特殊属性为前提,不再对正畸弓丝曲线上每个弯制点进行复杂度和归一化弯制点角距比判断进而确定是否需要对弯制点进行加权曲率误差率评价,不仅满足了实际正畸弓丝曲线的误差波动度评价,而且使得正畸弓丝评价算法的复杂程度降低,提高了评价效率。
附图说明
为了易于说明,本发明由下述的具体实施及附图作以详细描述。
图1为一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法流程图;
图2为理论正畸弓丝曲线在空间变换前的位姿示意图;
图3为m条实际正畸弓丝曲线在空间变换前的位姿示意图;
图4为第9个实际正畸弓丝曲线弯制点线误差波动度大于上限值的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之间连线的长度;
步骤二、实际正畸弓丝曲线弯制点曲率误差波动度和评价标定值的设定:
定义实际正畸弓丝曲线弯制点的曲率误差波动度,用符号σc表示,σc是对正畸弓丝曲线弯制点弯制形状稳定性的量化描述,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的曲率误差率,规定 TKi表示理论正畸弓丝曲线第i个弯制点的曲率,表示第j条实际正畸弓丝曲线第i个弯制点的曲率,表示m条实际正畸弓丝曲线的第i个弯制点的曲率误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度的上限值为(σc)max;定义评价标定值,用符号ξ表示,评价标定值ξ是对进行的误差波动度评价是否为曲率误差波动度的量化描述,当进行的误差波动度评价是曲率误差波动度时ξ=1,当进行的误差波动度评价不是曲率误差波动度时ξ=0,ξ的初始值为ξ=0;
步骤三、理论正畸弓丝曲线弯制点复杂度的设定:
定义理论正畸弓丝曲线弯制点复杂度,用符号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;
步骤四、理论正畸弓丝曲线弯制点复杂度及归一化弯制点角距比验证:
按照公式计算理论正畸弓丝曲线上第i个弯制点的复杂度,即1Cr表示理论正畸弓丝曲线上第1个弯制点Tp1的复杂度,i的初始值为i=1,通过比较取出iCr中的最小值(iCr)min,对条件(iCr)min>Crb进行验证;按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比,即表示理论正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,通过比较取出中的最小值对条件进行验证,具体为:
a)理论正畸弓丝曲线弯制点复杂度验证;
如果(iCr)min>Crb成立,说明处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpi,...,Tpn}内的最小弯制点复杂度仍大于所设定的复杂度分界值Crb,可知在所取的包含n个弯制点的理论正畸弓丝曲线上每个弯制点的复杂度都大于弯制点复杂度的分界值Crb,跳转至步骤五;
如果(iCr)min>Crb不成立,说明该此类理论正畸弓丝曲线弯制点的复杂度不大,但仍存在理论正畸弓丝曲线弯制点归一化弯制点角距比较大,而归一化弯制点密度较小这一情况,因此对理论正畸弓丝曲线归一化弯制点角距比进行验证,跳转至步骤四b);
b)理论正畸弓丝曲线归一化弯制点角距比验证;
如果成立,说明处于最终位姿下的理论正畸弓丝曲线弯制点信息集PT={Tp1,Tp2,Tp3,...,Tpi,...,Tpn}内弯制点的最小归一化弯制点角距比仍大于所设定的归一化弯制点分界值(TE*)b,可知在所取的包含n个弯制点的理论正畸弓丝曲线上每个弯制点的归一化弯制点角距比都大于归一化弯制点的分界值(TE*)b,跳转至步骤五;
如果TminE*>(TE*)b不成立,说明此类理论正畸弓丝曲线弯制点的复杂度和归一化弯制点角距比均不大于其相应的分界值,该评价方法不适用于此正畸弓丝曲线,则输出该评价方法不适用于此正畸弓丝曲线,正畸弓丝误差波动度评价结束;
步骤五、实际正畸弓丝曲线弯制点曲率误差波动度评价
根据iσc的计算公式计算m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度,判断iσc≤(σc)max是否成立,
具体为:
如果iσc≤(σc)max成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度在允许范围内,令ξ=1,跳转至步骤八;
如果iσc≤(σc)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,正畸弓丝评价结束;
步骤六、实际正畸弓丝曲线弯制点线误差波动度、角误差波动度的设定:
定义实际正畸弓丝曲线弯制点的线误差波动度,用符号σd表示,σd是对实际正畸弓丝曲线弯制点弯制距离稳定性的量化描述,规定m条实际正畸弓丝曲线第i个弯制点的线误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的线误差率,规定 表示m条实际正畸弓丝曲线第i个弯制点的线误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的线误差波动度的上限值为(σd)max;定义实际正畸弓丝曲线弯制点的角误差波动度,用符号σa表示,σa是对实际正畸弓丝曲线弯制点弯制角度稳定性的量化描述,规定m条实际正畸弓丝曲线第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σd和iσa的计算公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度和角误差波动度,i的初始值为i=1;
a)实际正畸弓丝曲线弯制点线误差波动度评价
具体为:
如果iσd≤(σd)max成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度在允许范围内,则跳转至步骤七b);
如果iσd≤(σd)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,正畸弓丝评价结束;
b)实际正畸弓丝曲线弯制点角误差波动度评价
具体为:
如果iσa≤(σa)max成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度在允许范围内,则令ξ=0,跳转至步骤八;
如果iσa≤(σa)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,正畸弓丝评价结束;
步骤八、判断m条实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对m条实际正畸弓丝曲线上所有弯制点进行评价,则令i=i+1,即表示对下一组实际正畸弓丝曲线弯制点进行评价,判断ξ=1是否成立,
具体为:
如果ξ=1成立,说明进行的是曲率误差波动度的评价,则跳转至步骤五;
如果ξ=1不成立,说明进行的不是曲率误差波动度的评价,则跳转至步骤七a);
如果i=n成立,则判断ξ=1是否成立,
具体为:
如果ξ=1成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行了曲率误差波动度评价,则可进行后续实际正畸弓丝曲线弯制点线误差波动度和角误差波动度评价,跳转至步骤六;
如果ξ=1不成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行了曲率误差波动度、线误差波动度和角误差波动度评价,且所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,则输出所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,正畸弓丝评价结束。
实施实例2:如图4所示,在m=20条包含n=16弯制点的实际正畸弓丝曲线上进行基于曲率误差预先判断的正畸弓丝误差波动度评价过程中,i的初始值为i=1,设定弯制点复杂度的分界值为Crb=0.5,归一化弯制点角距比的分界值(TE*)b=0.4,曲率误差波动度的上限值为(σc)max=0.07,根据步骤四中计算理论正畸弓丝弯制点复杂度的公式计算16个理论正畸弓丝曲线弯制点的复杂度,从中取出最小值为(12Cr)min=12Cr=0.57,通过比较得出(12Cr)min>Crb,则跳转至步骤五进行实际正畸弓丝曲线曲率误差波动度评价,根据iσc的计算公式,得出20条实际正畸弓丝曲线的第1个弯制点的曲率误差波动度1σc=0.05<(σc)max,令ξ=1后跳转至步骤八,根据判断此时1=16不成立且ξ=1,则令i=i+1,跳转至步骤五进行后续弯制点的误差波动度评价,根据此步骤计算得出后续弯制点的曲率误差波动度均小于上限值(σc)max,此时n=16,ξ=1,跳转至步骤八,根据判断此时16=16成立且ξ=1,则跳转至步骤六进行弯制点线误差波动度和角误差波动度的验证;设定弯制点线误差波动度的上限值(σd)max=0.06,弯制点角误差波动度的上限值(σa)max=0.09,在步骤七a)中,i的初始值为i=1,根据iσd的计算公式得出1σd=0.02<(σd)max,跳转至步骤七b),根据iσa的计算公式得出1σa=0.02<(σa)max,则令ξ=0并跳转至步骤八,根据判断此时1=16不成立且ξ=0,则令i=i+1并跳转至七a)进行后续弯制点的线误差波动度和角误差波动度评价,根据此步骤,计算得出20条实际正畸弓丝曲线的第2、第3、第4、第5、第6、第7、第8个弯制点的线误差波动度和角误差波动度均小于相应的上限值,在进行20条实际正畸弓丝曲线的第9个弯制点线误差波动度计算时,得出9σd=0.1>(σd)max,则输出实际正畸弓丝曲线第9个弯制点的线误差波动度超出允许范围,不再对后续实际正畸弓丝曲线弯制点进行评价,正畸弓丝评价结束。
实施实例3:如图5所示,在m=20条包含n=16弯制点的实际正畸弓丝曲线上进行基于曲率误差预先判断的正畸弓丝误差波动度评价过程中,i的初始值为i=1,设定弯制点复杂度的分界值为Crb=0.5,归一化弯制点角距比的分界值(TE*)b=0.4,曲率误差波动度的上限值为(σc)max=0.07,根据步骤四中计算理论正畸弓丝弯制点复杂度的公式计算16个理论正畸弓丝曲线弯制点的复杂度,从中取出最小值为(7Cr)min=7Cr=0.42,通过比较得出(7Cr)min<Crb,计算16个理论正畸弓丝曲线弯制点的归一化弯制点角距比,从中取出最小值为通过比较得出则跳转至步骤五进行实际正畸弓丝曲率误差波动度评价,根据iσc的计算公式,得出20条实际正畸弓丝曲线的第1个弯制点的曲率误差波动度1σc=0.06<(σc)max,令ξ=1后跳转至步骤八,根据判断此时1=16不成立且ξ=1,则令i=i+1,跳转至步骤五进行后续弯制点的误差波动度评价,根据此步骤计算得出后续弯制点的曲率误差波动度均小于上限值(σc)max,此时n=16,ξ=1,跳转至步骤八,根据判断此时16=16成立且ξ=1,则跳转至步骤六进行弯制点线误差波动度和角误差波动度的验证;设定弯制点线误差波动度的上限值(σd)max=0.06,弯制点角误差波动度的上限值(σa)max=0.09,在步骤七a)中,i的初始值为i=1,根据iσd的计算公式得出1σd=0.04<(σd)max,跳转至步骤七b),根据iσa的计算公式得出1σa=0.07<(σa)max,则令ξ=0并跳转至步骤八,根据判断此时1=16不成立且ξ=0,则令i=i+1并跳转至七a)进行后续弯制点的线误差波动度和角误差波动度评价,根据此步骤,计算得出20条实际正畸弓丝曲线的第2、第3、第4、第5、第6、第7、第8、第9、第10、第11、第12、第13、第14、第15、第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之间连线的长度;
步骤二、实际正畸弓丝曲线弯制点曲率误差波动度和评价标定值的设定:
定义实际正畸弓丝曲线弯制点的曲率误差波动度,用符号σc表示,σc是对正畸弓丝曲线弯制点弯制形状稳定性的量化描述,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的曲率误差率,规定 TKi表示理论正畸弓丝曲线第i个弯制点的曲率,表示第j条实际正畸弓丝曲线第i个弯制点的曲率,表示m条实际正畸弓丝曲线的第i个弯制点的曲率误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度的上限值为(σc)max;定义评价标定值,用符号ξ表示,评价标定值ξ是对进行的误差波动度评价是否为曲率误差波动度的量化描述,当进行的误差波动度评价是曲率误差波动度时ξ=1,当进行的误差波动度评价不是曲率误差波动度时ξ=0,ξ的初始值为ξ=0;
步骤三、理论正畸弓丝曲线弯制点复杂度的设定:
定义理论正畸弓丝曲线弯制点复杂度,用符号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;
步骤四、理论正畸弓丝曲线弯制点复杂度及归一化弯制点角距比验证:
按照公式计算理论正畸弓丝曲线上第i个弯制点的复杂度,即1Cr表示理论正畸弓丝曲线上第1个弯制点Tp1的复杂度,i的初始值为i=1,通过比较取出iCr中的最小值(iCr)min,对条件(iCr)min>Crb进行验证;按照公式计算理论正畸弓丝曲线上第i个弯制点的归一化弯制点角距比,即表示理论正畸弓丝曲线上第1个弯制点Tp1的归一化弯制点角距比,通过比较取出中的最小值对条件进行验证,具体为:
a)理论正畸弓丝曲线弯制点复杂度验证;
如果(iCr)min>Crb成立,跳转至步骤五;
如果(iCr)min>Crb不成立,则对理论正畸弓丝曲线归一化弯制点角距比进行验证,跳转至步骤四b);
b)理论正畸弓丝曲线归一化弯制点角距比验证;
步骤五、实际正畸弓丝曲线弯制点曲率误差波动度评价
根据iσc的计算公式计算m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度,判断iσc≤(σc)max是否成立,
具体为:
如果iσc≤(σc)max成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度在允许范围内,令ξ=1,跳转至步骤八;
如果iσc≤(σc)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的曲率误差波动度超出允许范围,正畸弓丝评价结束;
步骤六、实际正畸弓丝曲线弯制点线误差波动度、角误差波动度的设定:
定义实际正畸弓丝曲线弯制点的线误差波动度,用符号σd表示,σd是对实际正畸弓丝曲线弯制点弯制距离稳定性的量化描述,规定m条实际正畸弓丝曲线第i个弯制点的线误差波动度表示为 表示第j条实际正畸弓丝曲线第i个弯制点的线误差率,规定 表示m条实际正畸弓丝曲线第i个弯制点的线误差率的平均值,规定m条实际正畸弓丝曲线第i个弯制点的线误差波动度的上限值为(σd)max;定义实际正畸弓丝曲线弯制点的角误差波动度,用符号σa表示,σa是对实际正畸弓丝曲线弯制点弯制角度稳定性的量化描述,规定m条实际正畸弓丝曲线第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σd和iσa的计算公式计算m条实际正畸弓丝曲线第i个弯制点的线误差波动度和角误差波动度,i的初始值为i=1;
a)实际正畸弓丝曲线弯制点线误差波动度评价
具体为:
如果iσd≤(σd)max成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度在允许范围内,则跳转至步骤七b);
如果iσd≤(σd)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的线误差波动度超出允许范围,正畸弓丝评价结束;
b)实际正畸弓丝曲线弯制点角误差波动度评价
具体为:
如果iσa≤(σa)max成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度在允许范围内,则令ξ=0,跳转至步骤八;
如果iσa≤(σa)max不成立,说明m条实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,则输出实际正畸弓丝曲线第i个弯制点的角误差波动度超出允许范围,正畸弓丝评价结束;
步骤八、判断m条实际正畸弓丝曲线弯制点是否评价完毕:
判断i与实际正畸弓丝曲线弯制点的个数n是否相等,
具体为:
如果i=n不成立,说明未对m条实际正畸弓丝曲线上所有弯制点进行评价,则令i=i+1,即表示对下一组实际正畸弓丝曲线弯制点进行评价,判断ξ=1是否成立,
具体为:
如果ξ=1成立,说明进行的是曲率误差波动度的评价,则跳转至步骤五;
如果ξ=1不成立,说明进行的不是曲率误差波动度的评价,则跳转至步骤七a);
如果i=n成立,则判断ξ=1是否成立,
具体为:
如果ξ=1成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行了曲率误差波动度评价,则进行实际正畸弓丝曲线弯制点线误差波动度和角误差波动度评价,跳转至步骤六;
如果ξ=1不成立,说明已对m条实际正畸弓丝曲线上所有弯制点进行了曲率误差波动度、线误差波动度和角误差波动度评价,且所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,则输出所有实际正畸弓丝曲线弯制点的线误差波动度、角误差波动度和曲率误差波动度均在允许范围内,正畸弓丝评价结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210415167.9A CN115035197A (zh) | 2022-04-18 | 2022-04-18 | 一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210415167.9A CN115035197A (zh) | 2022-04-18 | 2022-04-18 | 一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115035197A true CN115035197A (zh) | 2022-09-09 |
Family
ID=83119386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210415167.9A Pending CN115035197A (zh) | 2022-04-18 | 2022-04-18 | 一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115035197A (zh) |
-
2022
- 2022-04-18 CN CN202210415167.9A patent/CN115035197A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114943685B (zh) | 一种基于轮廓降维法的正畸弓丝误差评价方法 | |
CN111588502A (zh) | 一种基于弯制点角距比和的正畸弓丝变半径圆域划分方法 | |
CN116644558A (zh) | 一种基于误差评定域的正畸弓丝误差评价方法 | |
CN114943058A (zh) | 一种基于位置误差判断的正畸弓丝误差波动度评价方法 | |
CN111588499A (zh) | 一种基于正畸弓丝弯制点密度的平面等半径圆域划分半径确定方法 | |
CN115024840B (zh) | 一种基于归一化弯制点密度的正畸弓丝误差率评价方法 | |
CN115024841B (zh) | 一种基于复杂度区分法的正畸弓丝误差评价方法 | |
CN114972184B (zh) | 一种基于权值比例法的正畸弓丝误差评价方法 | |
CN115035196A (zh) | 一种基于弯制点复杂度判断的正畸弓丝误差率评价方法 | |
CN114972183B (zh) | 一种基于残差平方和区间划分的正畸弓丝误差评价方法 | |
CN115035197A (zh) | 一种基于曲率误差预先判断的正畸弓丝误差波动度评价方法 | |
CN115040275B (zh) | 一种基于空间平移子坐标系卦限判断的正畸弓丝评价方法 | |
CN115024839B (zh) | 一种基于归一化弯制点角距比的正畸弓丝误差率评价方法 | |
CN115024838A (zh) | 一种基于弯制点复杂度判断的正畸弓丝误差波动度评价方法 | |
CN111588497B (zh) | 一种基于正畸弓丝弯制点角距比和的平面等半径圆域划分半径确定方法 | |
CN111588494B (zh) | 一种基于弯制点密度的正畸弓丝变角度划分方法 | |
CN115018245A (zh) | 一种基于弯制域节点的正畸弓丝弯制域评价方法 | |
CN117408943B (zh) | 一种基于曲线拟合法的正畸弓丝误差评价方法 | |
CN117218073B (zh) | 一种基于空间等距分割面的正畸弓丝误差评价方法 | |
CN114983593B (zh) | 一种基于向量共线的正畸弓丝弯制点误差评价方法 | |
CN111588503B (zh) | 一种基于弯制点密度的正畸弓丝变半径圆域划分方法 | |
CN117297806B (zh) | 一种基于阈值区间划分的正畸弓丝误差评价包络面半径确定方法 | |
CN117219256A (zh) | 一种基于弯制点曲率和降维角距偏差域的正畸弓丝误差评价方法 | |
CN114998190A (zh) | 一种基于点对区间的正畸弓丝评价方法 | |
CN114983594B (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 |