CN110811653A - Method for measuring tooth three-dimensional movement in orthodontic process - Google Patents

Method for measuring tooth three-dimensional movement in orthodontic process Download PDF

Info

Publication number
CN110811653A
CN110811653A CN201911000420.9A CN201911000420A CN110811653A CN 110811653 A CN110811653 A CN 110811653A CN 201911000420 A CN201911000420 A CN 201911000420A CN 110811653 A CN110811653 A CN 110811653A
Authority
CN
China
Prior art keywords
tooth
model
dimensional
treatment
axis
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
Application number
CN201911000420.9A
Other languages
Chinese (zh)
Inventor
刘鑫
汤文成
刘懋
严斌
吴斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201911000420.9A priority Critical patent/CN110811653A/en
Publication of CN110811653A publication Critical patent/CN110811653A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • A61B6/51
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50Clinical applications
    • A61B6/501Clinical applications involving diagnosis of head, e.g. neuroimaging, craniography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61CDENTISTRY; APPARATUS OR METHODS FOR ORAL OR DENTAL HYGIENE
    • A61C7/00Orthodontics, i.e. obtaining or maintaining the desired position of teeth, e.g. by straightening, evening, regulating, separating, or by correcting malocclusions

Abstract

The invention discloses a method for measuring tooth three-dimensional movement in an orthodontic process, which is characterized in that a CBCT tooth model and an intraoral scanning model at the initial treatment are superposed to construct a complete tooth model; matching and overlapping the intraoral scanning model after a period of treatment with the model at the initial moment, and comparing; isolating a single tooth for analysis, and constructing a local three-dimensional coordinate system of the target tooth; three reference points are marked at the positions of the dental crown supporting grooves on the front and back intraoral models respectively, a space transformation matrix is calculated according to the coordinates of the reference points, and then the three-dimensional translation and rotation motion of the teeth are calculated, so that the three-dimensional tooth movement is obtained. The same method analyzes all teeth and obtains the displacement of the entire dentition. The model is continuously collected for a plurality of times of intraoral scanning of the patient, and the continuous complete dentition three-dimensional displacement in the orthodontic process can be measured.

Description

Method for measuring tooth three-dimensional movement in orthodontic process
Technical Field
The invention relates to a technology for monitoring and measuring three-dimensional space movement of orthodontic teeth, and belongs to the technical field of orthodontic treatment detection.
Background
Orthodontic treatment is the adjustment of the coordination of the dentition itself, the dentition and the jaw bone by means of an orthodontic device. The ultimate goal is to achieve balance, stability and aesthetics of the oromandibular system.
In the orthodontic process, the position and the movement condition of the teeth of a patient can be detected and controlled in time, and the orthodontic treatment method has great significance for determining and adjusting the treatment scheme. While the morphology of the entire tooth location is often not well reflected by viewing only the visible crown portion. Therefore, prior to orthodontic treatment, physicians often perform CBCT scans of the patient to visualize the position of the teeth (including the root of the tooth). Since CBCT scanning has radiation to the patient and the patient needs to receive as little radiation as possible during the treatment, the partial motion of the teeth during the orthodontic procedure is usually measured clinically by using a ruler, cephalogram measurement and other methods. However, the movement of the tooth in three-dimensional space involves a complex process of rotation and translation. These simple measurement methods do not completely reflect the tooth movement, and in some complicated cases, such as embedded teeth and abnormal tooth movement, the doctor may still need to understand the tooth displacement and position by means of CBCT scanning of the patient during the operation to make and adjust the treatment measures in the next stage.
To date, no convenient method for measuring the three-dimensional space movement of teeth is proposed and applied to orthodontic clinics, so as to detect the tooth movement and reduce the need of CBCT scanning in operation.
Disclosure of Invention
The invention aims to solve the problems, provides a method for measuring and evaluating the complex three-dimensional movement of teeth in the clinical orthodontic treatment process, has simple and convenient digital operation, does not need a patient to bear additional radiation, and solves the problems that the displacement of the teeth (including tooth roots) cannot be mastered in time in the traditional orthodontic treatment process, the tooth movement is evaluated by excessively depending on the position of an observed tooth crown and the experience of a doctor, and the like.
In order to solve the above problems, the present invention provides a method for measuring three-dimensional movement of teeth in orthodontic process, comprising the steps of:
(1) superimposing the intraoral scan model of treatment initiation (T0) and the CBCT dental model;
(2) matching and superposing the intraoral scanning model after a period of treatment (T1) and the model at the initial moment;
(3) analyzing the single tooth, and establishing a local three-dimensional coordinate system for the target tooth on the obtained complete tooth model;
(4) marking three corresponding reference points of the target dental crown on the front model and the rear model respectively, and calculating a space transformation matrix of the target tooth position before and after the treatment at the stage according to the coordinates of the reference points;
(5) calculating the three-dimensional translation amount and the rotation angle of the tooth according to the transformation matrix;
(6) repeatedly analyzing each tooth to obtain displacement information of the complete dentition at the T1 moment;
(7) and continuously acquiring intraoral scanning models of the patient at different moments in the orthodontic process for analysis, and finally obtaining the displacement change information of all teeth of the patient in the whole treatment period.
As an improvement of the invention, the specific method of the step (1) is as follows: respectively preprocessing the intraoral scanning model and the CT tooth model, and matching and overlapping the preprocessed intraoral scanning model and the preprocessed CBCT tooth model.
As a modification of the present invention, the specific method in step (2) is to acquire the characteristics of the jaw and the third palatine fold of the anterior and posterior intraoral scan models, respectively, and to register the jaw and the third palatine fold with reference to each other.
As an improvement of the invention, the specific method of the step (3) is as follows: a local three-dimensional coordinate system of a single tooth is established by taking the central position of the tooth impedance as a coordinate origin, the mesial-distal direction as an X axis, the lingual-labial direction as a Y axis and the central axis of the tooth as a Z axis.
As an improvement of the present invention, the method for marking three corresponding reference points of the target crown on the anterior-posterior model in the step (4) marks three reference points on the initial model is as follows: and matching and marking the reference points on the initial model to the corresponding positions of the subsequent model to obtain three reference points on the second model.
As an improvement of the invention, in the step (4), a space transformation matrix TraM of the tooth body transformed from the position before the treatment to the position after the treatment for a period of time is solved by using a method of three-dimensional space coordinate system transformation in space geometry, as shown in the following formula
Figure BDA0002241142260000021
As a modification of the present invention, in step (5), the tooth three-dimensional translation amount and rotation angle are calculated, including extracting a translation vector from the spatial transformation matrix and calculating a rotation angle from the spatial transformation matrix; the translation vector is X, Y and Z, and the three terms respectively represent the translation distances of the impedance center of the tooth along the X-axis direction, the Y-axis direction and the Z-axis direction;
the rotation angle is calculated from the space transformation matrix, and the rotation angle comprises the steps of extracting a rotation matrix from the space transformation matrix and converting the rotation matrix into Euler angles according to a ZYX rotation sequence;
extracting a rotation matrix from the space transformation matrix, wherein the rotation matrix is formed by extracting three-dimensional matrix elements;
converting the rotation matrix into Euler angles in ZYX rotation order, comprising: according to the formula
Figure BDA0002241142260000022
Calculating the rotation angle around the X-axisAccording to the formula
Figure BDA0002241142260000024
Calculating a rotation angle ω about the Y-axis; according to the formula theta-atan 2 (r)21,r11) The rotation angle θ about the Z-axis is calculated.
Has the advantages that:
through the digital processing of the three-dimensional image, the three-dimensional movement of teeth (including tooth roots) in the orthodontic process is detected and measured in time, the mastering of orthodontic effects by clinicians and the adjustment of a later treatment scheme can be assisted, and the requirements of CBCT scanning in an operation can be effectively reduced.
Drawings
FIG. 1 is a flow chart of a method of monitoring and measuring tooth movement.
FIG. 2 is a front view of a CBCT dental model at T0.
Fig. 3 is a front view of the intraoral scan model at T0.
FIG. 4 is a schematic overlay of the intraoral scan model and CBCT dental model at T0.
Fig. 5 is a schematic diagram showing the comparison of the T0 model and the intraoral scan model at T1.
Fig. 6 is a schematic diagram of a local three-dimensional coordinate system of the patient's upper right cuspid teeth.
Fig. 7 is a schematic illustration of three reference points of the crown at the time of the marker T0.
Detailed Description
Illustrating according to what is contained in the claims
Example 1:
a method of measuring three-dimensional movement of teeth during orthodontic procedures, as shown in fig. 1, comprising the steps of:
(1) performing CBCT oral cavity scanning and intraoral laser scanning on the patient at the initial time of orthodontic treatment (T0), and reconstructing a CBCT tooth model (figure 2) and an intraoral scanning model (figure 3) of the patient; and (3) introducing the tooth model and the oral model into reverse software Geomagic, and performing segmentation, hole filling and smooth finishing to remove unnecessary structures such as gingiva, lips or cheeks.
The trimming aims to: the reconstructed model may have defects of concave-convex, noise point and the like, and the surface of the digital model is more real and smooth after the reconstruction.
(2) And the tooth model and the intraoral scan model are superimposed with the crown portion of each tooth as a reference, as shown in fig. 4.
(3) After the patient is treated in one stage, or when a doctor wants to pay attention to the current tooth displacement and tooth (including tooth root) shape of the patient, setting the time as T1, carrying out intraoral scanning on the patient, and obtaining an intraoral scanning model at the time; the intraoral model after a period of treatment is superimposed in registration with the initial intraoral model, with reference to the models at the jaw and third palatal fold, as shown in fig. 5.
(4) Respectively carrying out displacement analysis on the single tooth, taking the upper right cuspid as an example, as shown in fig. 6, taking the empirical tooth impedance center position as the origin of coordinates, the mesial-distal direction as the X axis, the lingual-labial direction as the Y axis, and the central axis of the tooth as the Z axis to establish a local three-dimensional coordinate system of the single tooth for describing tooth movement.
The tooth movement in the coordinate system can be expressed as
Figure BDA0002241142260000032
The corresponding physical and clinical meanings of each item are shown in table 1.
Figure BDA0002241142260000031
(5) And selecting three characteristic points as reference points on a bracket of a target tooth of the initial intraoral model, as shown in fig. 7, after the target dental crown is separately segmented, registering the segmented target dental crown with the target dental crown of the intraoral model after a period of treatment, and acquiring the positions of the three reference points of the target dental crown.
(6) And calculating a space transformation matrix of the positions of the target teeth before and after the treatment at the stage by using a three-dimensional space coordinate system transformation method in space geometry according to the coordinates of six reference points before and after the treatment, wherein the space transformation matrix is shown in the following formula.
Figure BDA0002241142260000041
The front three-dimensional matrix of the transformation matrix is a rotation matrix, and the fourth column is a translation vector.
(7) And calculating the tooth displacement according to the space transformation matrix.
The following formula is used to extract the movement terms (X, Y, Z) from the translation vector, and the following formula can be used to calculate the rotation terms from the rotation matrix
Figure BDA0002241142260000042
(X Y Z)=(txtytz)
Figure BDA0002241142260000043
Figure BDA0002241142260000044
θ=atan2(r21,r11)
(8) And repeating the steps 4-7, measuring all teeth of the patient and obtaining the displacement information of the whole dentition at the time T1.
(9) And repeating the steps 3-8, continuously collecting the intraoral scanning model of the patient in the treatment process for analysis, and finally obtaining all tooth displacement change information in the whole treatment period of the patient.

Claims (7)

1. A method of measuring three-dimensional movement of teeth during orthodontic procedures, comprising the steps of:
(1) superposing the intraoral scanning model of the treatment initial T0 and the CBCT tooth model;
(2) matching and superposing the intraoral scanning model of T1 after a period of treatment and the model at the initial moment;
(3) analyzing the single tooth, and establishing a local three-dimensional coordinate system for the target tooth on the obtained complete tooth model;
(4) marking three corresponding reference points of the target dental crown on the front model and the rear model respectively, and calculating a space transformation matrix of the target tooth position before and after the treatment at the stage according to the coordinates of the reference points;
(5) calculating the three-dimensional translation amount and the rotation angle of the tooth according to the transformation matrix;
(6) repeatedly analyzing each tooth to obtain displacement information of the complete dentition at the T1 moment;
(7) and continuously acquiring intraoral scanning models of the patient at different moments in the orthodontic process for analysis, and finally obtaining the displacement change information of all teeth of the patient in the whole treatment period.
2. The method for measuring the three-dimensional movement of teeth in orthodontic treatment according to claim 1, wherein the specific method in the step (1) is as follows: and respectively preprocessing the intraoral scanning model and the CBCT tooth model, and matching and overlapping the preprocessed intraoral scanning model and the preprocessed CBCT tooth model.
3. The method of claim 1, wherein the step (2) comprises obtaining the characteristics of the jaw and the third palatal fold of the anterior and posterior intraoral scan models, respectively, and registering the jaw and the third palatal fold with reference thereto.
4. The method for measuring the three-dimensional movement of teeth in orthodontic treatment according to claim 1, wherein the specific method in the step (3) is as follows: a local three-dimensional coordinate system of a single tooth is established by taking the central position of the tooth impedance as a coordinate origin, the mesial-distal direction as an X axis, the lingual-labial direction as a Y axis and the central axis of the tooth as a Z axis.
5. The method of claim 1, wherein the step (4) of marking three corresponding reference points of the target crown on the anterior-posterior model marks three reference points on the initial model by: and matching and marking the reference points on the initial model to the corresponding positions of the subsequent model to obtain three reference points on the second model.
6. The method of claim 1, wherein the transformation of the three-dimensional space coordinate system in the space geometry is used in step (4) to solve a spatial transformation matrix TraM for transforming the tooth from the pre-treatment position to the post-treatment position over time, as shown in
Figure FDA0002241142250000011
7. The method for measuring three-dimensional tooth movement in orthodontic treatment according to claim 1, wherein in the step (5), the three-dimensional tooth translation amount and the rotation angle are calculated, including extracting a translation vector from a spatial transformation matrix and calculating the rotation angle from the spatial transformation matrix; the translation vector is X, Y and Z, and the three terms respectively represent the translation distances of the impedance center of the tooth along the X-axis direction, the Y-axis direction and the Z-axis direction;
the rotation angle is calculated from the space transformation matrix, and the rotation angle comprises the steps of extracting a rotation matrix from the space transformation matrix and converting the rotation matrix into Euler angles according to a ZYX rotation sequence;
extracting a rotation matrix from the space transformation matrix, wherein the rotation matrix is formed by extracting three-dimensional matrix elements;
converting the rotation matrix into Euler angles in ZYX rotation order, comprising: according to the formula
Figure FDA0002241142250000021
Calculating the rotation angle around the X-axis
Figure FDA0002241142250000022
According to the formula
Figure FDA0002241142250000023
Calculating a rotation angle ω about the Y-axis; according to the formula theta-atan 2 (r)21,r11) The rotation angle θ about the Z-axis is calculated.
CN201911000420.9A 2019-10-21 2019-10-21 Method for measuring tooth three-dimensional movement in orthodontic process Pending CN110811653A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911000420.9A CN110811653A (en) 2019-10-21 2019-10-21 Method for measuring tooth three-dimensional movement in orthodontic process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911000420.9A CN110811653A (en) 2019-10-21 2019-10-21 Method for measuring tooth three-dimensional movement in orthodontic process

Publications (1)

Publication Number Publication Date
CN110811653A true CN110811653A (en) 2020-02-21

Family

ID=69549879

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911000420.9A Pending CN110811653A (en) 2019-10-21 2019-10-21 Method for measuring tooth three-dimensional movement in orthodontic process

Country Status (1)

Country Link
CN (1) CN110811653A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113034557A (en) * 2021-03-22 2021-06-25 西安交通大学口腔医院 Method and equipment for registering orthodontic anterior and posterior palatine fold
CN113052902A (en) * 2020-12-29 2021-06-29 上海银马科技有限公司 Dental treatment monitoring method
WO2021174479A1 (en) * 2020-03-05 2021-09-10 深圳先进技术研究院 Method and apparatus for calculating dental movement amount before and after oral orthodontic treatment
CN113712587A (en) * 2021-09-06 2021-11-30 吉林大学 Invisible orthodontic progress monitoring method, system and device based on oral scanning model
CN115153909A (en) * 2022-08-11 2022-10-11 南昌大学附属口腔医院(江西省口腔医院) Method for calculating tooth movement angle variation and visualizing tooth position angle value

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106327535A (en) * 2016-08-16 2017-01-11 苏州迪凯尔医疗科技有限公司 CBCT tooth root and intraoral scanning dental crown fusion method
CN108320325A (en) * 2018-01-04 2018-07-24 华夏天宇(北京)科技发展有限公司 The generation method and device of dental arch model
CN108765474A (en) * 2018-04-17 2018-11-06 天津工业大学 A kind of efficient method for registering for CT and optical scanner tooth model
CN109259875A (en) * 2018-09-06 2019-01-25 西安增材制造国家研究院有限公司 A kind of deviation detecting method of correction stage design model

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106327535A (en) * 2016-08-16 2017-01-11 苏州迪凯尔医疗科技有限公司 CBCT tooth root and intraoral scanning dental crown fusion method
CN108320325A (en) * 2018-01-04 2018-07-24 华夏天宇(北京)科技发展有限公司 The generation method and device of dental arch model
CN108765474A (en) * 2018-04-17 2018-11-06 天津工业大学 A kind of efficient method for registering for CT and optical scanner tooth model
CN109259875A (en) * 2018-09-06 2019-01-25 西安增材制造国家研究院有限公司 A kind of deviation detecting method of correction stage design model

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
张东霞等: "基于口腔计算机断层扫描图像与激光扫描图像融合的牙齿三维模型重构" *
张东霞等: "基于口腔计算机断层扫描图像与激光扫描图像融合的牙齿三维模型重构", 《生物医学工程学杂志》 *
李占利等: "虚拟正畸治疗中的错位牙齿自动排列方法", 《图学学报》 *
李林峰等: "三次均匀B样条在工业机器人轨迹规划中的应用研究", 《科学技术与工程》 *
焦培峰等: "基于激光扫描与CT建立带牙根的三维数字化牙颌模型", 《中国临床解剖学杂志》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021174479A1 (en) * 2020-03-05 2021-09-10 深圳先进技术研究院 Method and apparatus for calculating dental movement amount before and after oral orthodontic treatment
CN113052902A (en) * 2020-12-29 2021-06-29 上海银马科技有限公司 Dental treatment monitoring method
CN113034557A (en) * 2021-03-22 2021-06-25 西安交通大学口腔医院 Method and equipment for registering orthodontic anterior and posterior palatine fold
CN113034557B (en) * 2021-03-22 2023-12-05 西安交通大学口腔医院 Registration method and equipment for orthodontic anterior and posterior palate folds
CN113712587A (en) * 2021-09-06 2021-11-30 吉林大学 Invisible orthodontic progress monitoring method, system and device based on oral scanning model
CN113712587B (en) * 2021-09-06 2023-07-18 吉林大学 Invisible orthodontic progress monitoring method, system and device based on oral cavity scanning model
CN115153909A (en) * 2022-08-11 2022-10-11 南昌大学附属口腔医院(江西省口腔医院) Method for calculating tooth movement angle variation and visualizing tooth position angle value

Similar Documents

Publication Publication Date Title
CN110811653A (en) Method for measuring tooth three-dimensional movement in orthodontic process
US20180153659A1 (en) System and method for effective planning, visualization, and optimization of dental restorations
EP2134290B1 (en) Computer-assisted creation of a custom tooth set-up using facial analysis
Lin et al. Artifact-resistant superimposition of digital dental models and cone-beam computed tomography images
US8731280B2 (en) Virtual cephalometric imaging
US7080979B2 (en) Method and workstation for generating virtual tooth models from three-dimensional tooth data
US20060127854A1 (en) Image based dentition record digitization
CN113509281A (en) Historical scan reference for intraoral scanning
CN104715099B (en) Virtual articulator design method for tooth crown bridge working bite relation face
CN106570859B (en) Tooth image analysis system
JPWO2006033483A1 (en) Human body information extraction device, human body photographing information reference plane conversion method, and cross-section information detection device
WO2006065955A2 (en) Image based orthodontic treatment methods
CN103908352B (en) For generating the method and system of digital virtual jaw frame
CN112451151B (en) Orthodontic model establishing method utilizing mixed reality technology
CN104867131A (en) Dental crown data extraction method based on digital model
Chen et al. Quantification of tooth displacement from cone-beam computed tomography images
TWI397402B (en) An integration method of dental implant positioning and implant guide planning
Barone et al. Geometrical modeling of complete dental shapes by using panoramic X-ray, digital mouth data and anatomical templates
CN116958169A (en) Tooth segmentation method for three-dimensional dental model
CN115690045A (en) Quantitative assessment method for bone increment before and after periodontitis treatment based on curved surface fault slice
Lin et al. A new method for the integration of digital dental models and cone-beam computed tomography images
CN113397585B (en) Tooth body model generation method and system based on oral CBCT and oral scan data
Anes et al. Mounting digital casts on a virtual articulator by using two-dimensional facial photographs with a facebow: A technique
Zou et al. Research of dental arch curve extraction and application
Enciso et al. 3D visualization of the craniofacial patient: Volume segmentation, data integration and animation

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200221