KR102238483B1 - Calculation Method of Body Similarity Index, and Method and Apparatus for Surgery Planning Simulation Using Thereof - Google Patents

Abstract

The present invention relates to a method for calculating a body similarity index and a method and apparatus for simulating a surgery plan using the same. The method for calculating a body similarity index according to the present invention comprises segmenting a body region of interest from 3D medical image data acquired by an image acquisition unit. The step of mirroring the first body region located on one side of the reference plane among the segmented body regions based on the reference plane to overlap the second body region located on the other side of the reference plane, and corresponding to the first body region The similarity between the first body region and the second body region is calculated by using the first volume to be determined, the second volume corresponding to the second body region, and the third volume in which the second body region and the mirrored first body region are overlapped. It includes the step of.

Description

Calculation Method of Body Similarity Index, and Method and Apparatus for Surgery Planning Simulation Using Thereof}

The present invention relates to a method of calculating a body similarity index and a method and apparatus for simulating a surgical plan using the same.

For patients with facial asymmetry, it is important to evaluate facial asymmetry in order to plan the operation or evaluate the surgical outcome. However, most conventional methods using a 3D laser or optical scan to evaluate facial asymmetry depend on several markers of soft tissue, and are not sufficient to intuitively and quantitatively evaluate the degree of facial asymmetry by mainly analyzing only the surface morphology.

U.S. Patent No. 9,808,326 (Registration Date: November 7, 2017)

The technical problem to be solved by the present invention is to provide a method for calculating a body similarity index and a method and apparatus for simulating a surgical plan using the same.

The surgical plan simulation method according to the present invention for solving the above technical problem includes the steps of: (a) segmenting the maxillary region and the mandible region from the 3D medical image data acquired by the image acquisition unit, (b) the segmented Dividing the mandibular region into an anterior fracture region and a left and right posterior fracture region, (c) passing through the midpoint (M) of the line segment connecting the left and right mental foramen extracted from the segmented mandible region. Creating a modified AMP (absolute mandibular midsagittal plane) in which a center point (T) is set in a plane perpendicular to the line segment, (d) FHP (Frankfort horizontal plane), rotate the modified AMP as much as it is rotated parallel to the horizontal plane, and rotate the modified AMP parallel to the midsagittal plane (MSP) based on the midpoint (M) on an axial plane, The step of moving and positioning the anterior fracture segment region as much as the modified AMP is moved in parallel so that the center point T is located at the MSP, and (e) the anterior fracture segment region has a maximum similarity on the left and right sides based on the MSP. It includes the step of positioning to be.

The method may further include (f) adjusting the anterior bone fracture region to a final position by using the normal average value of the head measurements obtained in advance.

The method includes: (g) positioning the condyle center point of the left posterior fracture region on the left vertical line passing through the uppermost point of the left articular fossa of the maxillary region, and placing the condyle center point of the right posterior fracture region on the maxillary region. Positioning the right joint and the uppermost point on the right vertical line, and (h) while moving the left and right posterior fracture region in a vertical direction in a state in which the condyle center points of the left and right posterior fracture regions are respectively located on the left and right vertical lines. It may further include the step of finally positioning the left and right posterior fracture region to a position where the similarity is maximized.

The central point (T) is set as the midpoint of the anterior fracture segment located in the plane, or is the midpoint of a Menton, point B, or genial tubercle projected perpendicularly to the plane. Can be set.

The method may further include: (i) finally positioning the anterior fracture region and the left and right posterior fracture regions, and then removing a region where the anterior fracture region and the left and right posterior fracture regions overlap. .

Based on the MSP, the left and right similarity of the anterior fracture region, the similarity of the left and right posterior fracture region, and the left and right similarity of the entire mandible may be obtained, respectively.

The method may further include performing a virtual surgery on the maxillary region before step (b).

The method may further include generating data for a final wafer based on positional data of the maxillary bone region, the anterior fracture region, and the left and right posterior fracture region after the step (i).

_{The similarity (SI) obtained based on the MSP is a first volume (V 1} ) corresponding to the left region and a first volume (V 1) corresponding to the right region by mirroring and overlapping the left region and the right region based on the MSP. It can be calculated by the equation SI = 2×V ₃ /(V ₁ +V ₂ ) by using the 2 volume (V _{2 ) and the third volume (V 3} ) in which the left area and the right area overlap.

The surgical planning simulation apparatus according to the present invention for solving the above technical problem comprises segmenting a maxillary region and a mandible region from 3D medical image data acquired by an image acquisition unit, and using the segmented mandible region as an anterior fracture segment. , The center point (T) is divided into left and right posterior fracture regions, and a center point (T) is placed on a plane perpendicular to the line segment while passing the midpoint (M) of the line segment connecting the left and right mental foramen extracted from the segmented mandibular region. A set modified AMP (absolute mandibular midsagittal plane) is created, and the modified AMP is rotated by rotating the line segment parallel to FHP (Frankfort horizontal plane) on a coronal plane based on the mid point (M). And, after rotating the modified AMP on the axial plane to be parallel to the midsagittal plane (MSP) based on the midpoint (M), the modified AMP so that the center point (T) is located at the MSP. And a data processing unit for moving and positioning the anterior bone fragment region as much as is parallelly moved, and positioning the anterior bone fragment region such that the left and right similarity is maximized with respect to the MSP.

The method for calculating a body similarity index according to the present invention for solving the above technical problem includes the steps of segmenting a body region of interest from 3D medical image data obtained by an image acquisition unit, and a reference plane set in advance among the segmented body regions. Mirroring a first body region located on one side of the reference plane based on the reference plane to overlap a second body region located on the other side of the reference plane, and a first volume corresponding to the first body region, the second Calculating a similarity between the first body region and the second body region using a second volume corresponding to the body region and a third volume in which the second body region and the mirrored first body region overlap. Includes.

The degree of similarity between the first body region and the second body region may be calculated by the following equation.

SI = 2×V ₃ /(V ₁ +V ₂ )

Here, SI may be a degree of similarity between the first body region and the second body region, V ₁ may be the first volume, V ₂ may be the second volume, and V ₃ may be the third volume.

A computer-readable recording medium in which a program for executing the method described above is recorded in a computer according to the present invention for solving the above technical problem may be further included.

Using the body similarity index according to the present invention, a user, such as a clinician, can intuitively and quantitatively evaluate the degree of body asymmetry including facial asymmetry. By providing a simulation of a surgery plan using a body similarity index according to the present invention, it may help to establish a surgery plan.

1 is a block diagram of a surgical plan simulation apparatus according to an embodiment of the present invention.
2 is a flow chart for explaining the operation of the surgery plan simulation apparatus according to an embodiment of the present invention.
3 to 12 are views for explaining the procedure of performing a virtual mandible surgery according to the present invention.
13 is a diagram illustrating a method of calculating a body region similarity index for body asymmetry evaluation according to an embodiment of the present invention.
14 is a diagram for describing a method of calculating a body region similarity index for facial asymmetry evaluation according to an embodiment of the present invention.
15 is a view for explaining a method of calculating a body region similarity index for facial asymmetry evaluation according to another embodiment of the present invention.

Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings so that those of ordinary skill in the art may easily implement the present invention.

1 is a block diagram of a surgical plan simulation apparatus according to an embodiment of the present invention.

Referring to FIG. 1, the apparatus for simulating an operation plan according to the present invention may include an image acquisition unit 100, a data processing unit 200, and a user interface unit 300.

The image acquisition unit 100 is a device that acquires a three-dimensional (3D) medical image of a patient and may be a computed tomography (Computed Tomography) device. Here, the 3D medical image is 3D acquired in types such as CT (Computed Tomography), CBCT (Cone Beam Computed Tomography), MRI (Magnetic Resonance Imaging), PET (Positron Emmission Tomography), and SPECT (Single-Photon Emmission Computed Tomography). It may be image data.

The data processing unit 200 performs data processing such as calculating a similarity index between the body region segment of interest and the body regions divided based on the reference plane from the medical image data acquired by the image acquisition unit 100, and performing a virtual operation on the body region. For this, one or more processors and memory modules may be included.

The user interface unit 300 may output various types of information related to the operation of the operation plan simulation apparatus. For example, the user interface unit 300 may display the contents of a virtual surgery performed on the screen based on the 3D medical image data acquired by the image acquisition unit 100 as an image. In addition, the user interface unit 300 may receive various commands related to the operation of the operation plan simulation apparatus from the user. The user interface unit 300 may include an input/output module such as a monitor, a speaker, a touch panel, a mouse, and a keyboard.

2 is a flow chart for explaining the operation of the surgery plan simulation apparatus according to an embodiment of the present invention.

Referring to FIGS. 1 and 2, the image acquisition unit 100 may acquire three-dimensional (3D) medical image data photographing a patient's face (S205).

Next, the data processing unit 200 may segment the maxillary region and the mandible region from the 3D medical image data acquired by the image acquisition unit 100 (S210 ). In step S210, the data processing unit 200 may automatically segment the maxillary region and the mandible region from the 3D medical image data using a neural network model such as deep learning.

The neural network model may be in the form of a machine learning algorithm such as a convolution neural network (CNN). Convolutional neural networks can be implemented with deep CNN models such as AlexNet, VGGNet, GoogLeNet, ResNet, and the like. A technique for segmenting an ROI in an input image using a trained neural network model is well known to those skilled in the art, and thus a detailed description thereof will be omitted.

The data processing unit 200 may extract a 3D cephalometric reference point and an additional reference point from the 3D medical image data (S215). 3D head measurement reference points include Porion, Orbitale, Nasion, Sella, Menton, point B, and genial tubercle. I can. In addition, the additional reference point may include an uppermost point of the left and right joints (articular fossa), a predetermined reference point in the maxilla and mandible.

Reference points are automatically extracted from the data processing unit 200 using the neural network model trained in advance in step S215, or 3D medical image data is displayed through the user interface unit 300, and reference points from a user such as a clinician It can be implemented to select and extract them. In step S215, the data processing unit 200 may extract 2D head measurement reference points by using a neural network model trained in advance on the medical image data acquired by the image acquisition unit 100.

Thereafter, the data processing unit 200 may set a reference plane using the reference points extracted in step S215 (S220). The reference plane set in step S220 may be a Frankfort horizontal plane (FHP), a midsagittal plane (MSP), an absolute mandibular midsagittal plane (AMP), and the like. The FHP is the plane passing through the left and right porions and the ophthalmic point. The MSP is the plane perpendicular to the FHP, passing through the non-peripheral and cellar. AMP is the plane through the menton, point B, and nodule.

According to an embodiment, steps S215 and S220 may be performed before step S210 or may be performed simultaneously.

The data processing unit 200 may perform a virtual operation on the maxilla (S225). For example, the data processing unit 200 may virtually perform Le Fort I maxillary bone amputation on the segmented maxillary region. The segmented maxillary region can be moved in 3D based on FHP and MSP. The optimal position of the maxillary A point, ANS (Anterior nasal spine), PNS (Posterior nasal spine), etc. can be determined using 2D or 3D head measurement analysis. In addition, the segmented maxillary region can be translated back and forth on the MSP, and a pitch and an upward movement are performed based on point A to move to an optimal position. When moving the maxilla, the reference plane is the plane connecting the mesiobuccal cusp tips of the maxillary bilateral central incisors and the maxillary bilateral first molars. Reference points for determining these planes may be determined manually or through deep learning. The yaw movement can be performed so that the center line of the maxillary reference plane coincides with the MSP. The roll movement can be performed so that the maxillary reference plane is horizontal with the FHP around point A. In addition, after the maxillary bone is moved in three dimensions, the upper part of the maxillary bone, that is, the part overlapping with the skull may be removed (trimming).

In step S225, the virtual surgery on the maxilla may be performed by applying a maxillary surgery method different from that described above. Meanwhile, in the case of a patient who does not need maxillary surgery, step S225 may be omitted.

Next, the data processing unit 200 may generate data for an intermediate wafer based on the final position data of the maxillary region on which the virtual surgery was performed (S227). Intermediate wafer refers to a wafer that holds the patient's maxilla.

In addition, the data processing unit 200 may perform a virtual operation on the mandible (S230).

3 to 13 are views for explaining the procedure of performing a virtual mandible surgery according to the present invention.

The data processing unit 200 may divide the segmented mandible region into an anterior fracture region, a left fracture region, and a right posterior fracture region (S231). For example, as illustrated in FIG. 3, by virtually performing intraoral vertical ramus osteotomy (IVRO), the mandible region is an anterior fracture region 10, a left fracture region 21, and a right posterior fracture region. It can be divided into (22). Mandibular vertical fracture is a surgery that vertically divides both mandibles based on the line connecting the left and right antegonial notch and sigmoid notch of the mandible. Of course, as illustrated in FIG. 4, a mandibular sagittal split ramus osteotomy (SSRO) is virtually performed to determine the segmented mandible region in an anterior fracture region 10, a left fracture region 21, and a right posterior region. It can also be divided into a fracture section 22. In addition, it is possible to divide the segmented mandibular region into an anterior fracture region, a left fracture region, and a right posterior fracture region by performing a virtual operation according to another surgical method in step S231.

Next, referring to FIG. 2 again, the data processing unit 200 may generate a modified AMP (absolute mandibular midsagittal plane) (S232).

5, while passing through the midpoint (M) of the line segment (L) connecting the left and right teeth (mental foramen) (P _mf, P _{mf') extracted from the segmented mandible region, perpendicular to the line segment (L).} A center point T (not shown) can be set on the plane A. Hereinafter, the plane (A) in which the center point (T) is set is defined as a modified AMP.

The central point T may be set as the midpoint of the anterior fracture segment 10 located in the plane A. Alternatively, the central point T may be set as the midpoint of the menton, point B, and genial tubercle projected perpendicular to the plane A.

Referring back to FIG. 2, the data processing unit 200 may move and position the anterior bone fragment region 10 by rotating and parallel moving the modified AMP to overlap the MSP (S233). A detailed description of step S235 is as follows.

First, as illustrated in FIG. 6, the modified AMP(A) is rotated by rotating the line segment L in parallel with the Frankfort horizontal plane (FHP) on the coronal plane based on the midpoint M. The modified AMP is not in a state parallel to the midsagittal plane (MSP) as illustrated in FIG. 7. Next, the modified AMP on the axial plane is rotated to be parallel to the midsagittal plane (MSP) based on the midpoint M (see FIG. 8). Here, the axial plane is a horizontal (transverse) plane, which is a plane parallel to the FHP. Then, the modified AMP is moved in parallel so that the center point T set in the modified AMP (A) is located in the MSP (see FIG. 9).

6 to 9 by rotating and translating the modified AMP (A) so as to move to overlap with the MSP rotates and translates the anterior bone fragment region 10.

Referring to FIG. 2 again, after the step S233, the data processing unit 200 may rotate the anterior bone fracture region 10 within a certain angle and position the left and right sides to maximize the similarity of the left and right sides based on the MSP (S234). . For example, the anterior fracture segment 10 may be positioned at a position where the left and right similarity is maximized based on the MSP while yawing and rolling within the range of ± 15°. The yaw axis may be an axis perpendicular to a horizontal (transverse) plane while passing through the center point T, and the rolling axis may be an axis perpendicular to a coronal plane while passing through the center point T.

Next, the data processing unit 200 may adjust the anterior bone fracture region 10 to the final position by using the normal average of the head measurements obtained in advance (S235). For example lower point (pogonion) the horizontal distance between (also in 10 P _O) and Nperp plane (handy that the passing each vertical tubular (coronal) orientation plane in MSP and FHP) (FHP in step (S235) After moving the anterior fracture section 10 anteriorly or posteriorly so that the parallel distance, D) is in the range of approximately 0-8 mm (the exact amount of movement may vary from person to person), the upper and lower anterior teeth horizontal overjet and vertical Each of the overbite is 2-3 mm (refer to FIG. 11), and the position of the anterior fracture segment 10 may be finally adjusted and positioned so that the maxillary posterior and the lower posterior teeth are in contact. The exact amount of movement can be determined using a normal value according to a two-dimensional or three-dimensional profile or anteroposterior standard measurement analysis method, and may vary according to the characteristics of an individual's skeletal malocclusion. Of course, the range may be set or adjusted according to the medical judgment of the clinician who is the user.

Next, the data processing unit 200 may move the left and right posterior bone fracture regions 21 and 22 to a designated position (S236). In the case of performing mandibular vertical fracture (IVRO) on the mandible, it may not be necessary to specify an exact position in step S236, since it is possible to perform an osteotomy, but do not perform a separate internal fixation, and immediately perform intermaxillary fixation. However, in the case of mandibular sagittal fracture (SSRO), the left and right posterior fracture regions (21, 22) and the anterior fracture region (10) need to be fixed and fixed together, so the left and right posterior fracture regions (21, 22) It is necessary to specify the exact location of.

Position designation for the left and right posterior fracture regions 21 and 22 in step S236 may be performed using the following method. As illustrated in FIG. 12, first, the condyle center point P21 of the left posterior fracture region 21 is positioned on the left vertical line L3 passing through the uppermost point P31 of the left articular fossa of the maxillary region, and the right posterior The condylar central point P22 of the bone fracture region 22 may be positioned on a right vertical line L4 passing through the right joint and the uppermost point P32 of the maxillary region. And as illustrated in FIG. 13, in a state in which the condyle center points P21 and P22 of the left and right posterior fracture regions 21 and 22 are respectively positioned on the left and right vertical lines L3 and L4, the left and right posterior fracture regions 21 and 22 ) May be moved in the vertical direction within a certain range, and the similarity of the left and right posterior fracture regions 21 and 22 may be finally positioned at a position where the degree of similarity is maximized based on the MSP.

The range in which the left posterior fracture region 21 and the right posterior fracture region 22 can move in the vertical direction may be preset by the user. For example, it can be set by using the average range of the condyle and the gap between the joint and the condyle examined for a normal person. Of course, the range may be set or adjusted according to the medical judgment of the clinician who is the user. In addition, while moving the left and right posterior fracture regions 21 and 22 in the vertical direction, yawing, pitching, and rolling are performed within a certain angle (for example, ±15°), and a position at which the degree of similarity is maximized based on the MSP may be designated.

Referring to FIG. 2 again, the data processing unit 200 finally positions the anterior fracture region 10 and the left and right posterior fracture regions 21 and 22, and then the anterior fracture region 10 and the left and right posterior fracture regions. The area where (21, 22) overlaps may be removed (S237).

Finally, the data processing unit 200 may generate data for the final wafer based on the positional data of the maxillary bone region, the anterior fracture region, and the left and right posterior fracture region (S240). The final wafer refers to a wafer that holds the mandible position.

14 is a diagram for describing a method of calculating a body region similarity index for body asymmetry evaluation according to an embodiment of the present invention.

Referring to FIG. 14, when the body region to be evaluated for body asymmetry is the mandible, the similarity index of the left and right regions may be calculated by dividing the mandible based on an absolute mandibular midsagittal plane (AMP). Here, AMP may be defined as a plane consisting of a menton of the mandible, point B, and a genial tubercle.

14, the mandible is divided into a first area (A) and a second area (B) based on the AMP, and the area (B) located on one side is mirrored, and the area (A) located on the other side. When superimposed, a non-overlapping volume is created in which the first region A and the second region B are overlapped with each other (A∩B).

As the first region (A) and the second region (B) become symmetrical based on the AMP, the overlapped region increases and the non-overlapping region decreases.

Therefore, the similarity index of the two regions (A, B) of the mandible divided on the basis of the AMP can be obtained using Equation 1 below.

[Equation 1]

SI = 2×V ₃ /(V ₁ +V ₂ )

Here, SI is the similarity index of the two regions (A, B) of the mandible divided based on AMP, V ₁ is the volume of the first region (A), V ₂ is the volume of the second region (B), V ₃ may be the volume of the overlapped area (A∩B).

When the two areas (A, B) are perfectly symmetrical based on AMP, the SI will be '1', and if it is completely asymmetric, the SI will be '0' because there is no overlapping area.

15 is a view for explaining a method of calculating a body region similarity index for facial asymmetry evaluation according to another embodiment of the present invention.

Referring to FIG. 15, the similarity index may be calculated by applying Equation 1 to the two regions C and D of the mandible divided based on the MSP.

Of course, for a body region other than the mandible, the degree of similarity between the two body regions divided based on a predetermined reference plane may be calculated using Equation 1.

Meanwhile, while simulating the operation plan as in the embodiment of FIG. 2, not only can the body similarity index obtained according to the present invention be used, but also the body similarity index using Equation 1 can be obtained and compared before and after the operation to evaluate the operation result. have.

For example, in the surgical simulation process as in the previous steps (S234, S236), in addition to obtaining the left and right similarity of the anterior fracture region based on the MSP and the similarity between the left posterior fracture region and the right posterior fracture region based on the MSP, surgery In order to evaluate the surgical results before or after surgery, it is also possible to obtain and compare the left and right similarities of the anterior fracture region and the similarity between the left posterior fracture region and the right posterior fracture region based on the MSP.

The embodiments described above may be implemented as a hardware component, a software component, and/or a combination of a hardware component and a software component. For example, the devices, methods, and components described in the embodiments are, for example, a processor, a controller, an arithmetic logic unit (ALU), a digital signal processor, a microcomputer, a field programmable gate (FPGA). array), programmable logic unit (PLU), microprocessor, or any other device capable of executing and responding to instructions. The processing device may execute an operating system (OS) and one or more software applications executed on the operating system. Further, the processing device may access, store, manipulate, process, and generate data in response to the execution of software. For the convenience of understanding, although it is sometimes described that one processing device is used, one of ordinary skill in the art, the processing device is a plurality of processing elements and/or a plurality of types of processing elements. It can be seen that it may include. For example, the processing device may include a plurality of processors or one processor and one controller. In addition, other processing configurations are possible, such as a parallel processor.

The software may include a computer program, code, instructions, or a combination of one or more of these, configuring the processing unit to behave as desired or processed independently or collectively. You can command the device. Software and/or data may be interpreted by a processing device or, to provide instructions or data to a processing device, of any type of machine, component, physical device, virtual equipment, computer storage medium or device. It can be permanently or temporarily embody. The software may be distributed over networked computer systems and stored or executed in a distributed manner. Software and data may be stored on one or more computer-readable recording media.

The method according to the embodiment may be implemented in the form of program instructions that can be executed through various computer means and recorded in a computer-readable medium. The computer-readable medium may include program instructions, data files, data structures, etc. alone or in combination. The program instructions recorded on the medium may be specially designed and configured for the embodiment, or may be known and usable to those skilled in computer software. Examples of computer-readable recording media include magnetic media such as hard disks, floppy disks, and magnetic tapes, optical media such as CD-ROMs and DVDs, and magnetic media such as floptical disks. -A hardware device specially configured to store and execute program instructions such as magneto-optical media, and ROM, RAM, flash memory, and the like. Examples of program instructions include not only machine language codes such as those produced by a compiler, but also high-level language codes that can be executed by a computer using an interpreter or the like. The hardware device described above may be configured to operate as one or more software modules to perform the operation of the embodiment, and vice versa.

As described above, although the embodiments have been described by the limited drawings, a person of ordinary skill in the art can apply various technical modifications and variations based on the above. For example, the described techniques are performed in a different order from the described method, and/or components such as systems, structures, devices, circuits, etc. described are combined or combined in a form different from the described method, or other components Alternatively, even if substituted or substituted by an equivalent, an appropriate result can be achieved.

Claims (19)

(a) segmenting the maxillary region and the mandible region from the 3D medical image data obtained by the data processing unit by the image acquisition unit,
(b) dividing the segmented mandible region into an anterior fracture region and a left and right posterior fracture region, by the data processing unit,
(c) A modified AMP in which the data processing unit sets the center point T on a plane perpendicular to the line segment while passing the midpoint M of the line segment connecting the left and right side foramen extracted from the segmented mandible region. creating an absolute mandibular midsagittal plane),
(d) The data processing unit rotates the modified AMP by rotating the line segment in parallel with a Frankfort horizontal plane (FHP) on a coronal plane with respect to the midpoint M, and an axial plane ), the modified AMP is rotated to be parallel to the midsagittal plane (MSP) based on the midpoint (M), and then the modified AMP is moved in parallel so that the center point (T) is located at the MSP. Moving and positioning the fracture section, and
(e) positioning the data processing unit so that the left and right similarity is maximized based on the MSP, by the data processing unit. In claim 1,
(f) The data processing unit further comprises the step of adjusting the anterior fracture region to a final position by using the normal average value of the head measurements obtained in advance. In claim 2,
(g) The data processing unit locates the condyle center point of the left posterior fracture region on the left vertical line passing through the uppermost point of the left articular fossa of the maxillary region, and the condyle center point of the right posterior fracture region is the maxillary region. Positioning on the right vertical line passing through the right joint and the top point of the
(h) In a state in which the data processing unit positions the condyle center points of the left and right posterior fracture regions on the left and right vertical lines, the similarity of the left and right posterior fracture regions is maximized while moving the left and right posterior fracture regions in a vertical direction. Surgical planning simulation method further comprising the step of finally positioning the position to be. In claim 3,
The center point (T) is,
It is set as the midpoint of the anterior fracture segment located in the plane, or
A method of simulating a surgical plan that is set as the midpoint of a menton, point B, and a genial tubercle projected perpendicular to the plane. In claim 3,
(i) after the data processing unit finally positions the anterior fracture region and the left and right posterior fracture regions, the operation plan further comprises removing an area where the anterior fracture region and the left and right posterior fracture regions overlap Simulation method. In claim 5,
Surgery planning simulation method in which the data processing unit obtains the left and right similarity of the anterior fracture region, the similarity of the left and right posterior fracture region, and the left and right similarity of the entire mandible based on the MSP. In claim 5,
The operation planning simulation method further comprising the step of performing a virtual operation on the maxillary region before the data processing unit (b) step. In claim 7,
The operation plan simulation method further comprising the step of generating data for the final wafer based on the position data of the maxillary bone region, the anterior fracture region, and the left and right posterior fracture region after the data processing step (i). In claim 6,
The similarity (SI) calculated based on the MSP is,
The left area and the right area are mirrored and overlapped based on the MSP, and a corresponding first volume (V ₁ ) corresponding to the left area, a second volume corresponding to the right area (V ₂ ), and the left area Surgery plan simulation method calculated by the equation SI = 2 × V ₃ / (V ₁ + V ₂ ) using the third volume (V _{3) overlapping the right area.} The maxillary region and the mandible region are segmented from the 3D medical image data acquired by the image acquisition unit, and the segmented mandibular region is divided into an anterior fracture segment region and a left and right posterior fracture segment region, and in the segmented mandible region A plane perpendicular to the line segment is created while passing through the midpoint (M) of the line segment connecting the extracted left and right teeth (mental foramen), and a center point (T) obtained based on the anterior fracture segment is set in the plane. A modified AMP (absolute mandibular midsagittal plane) is created, and the modified AMP is rotated by rotating the line segment in parallel with the Frankfort horizontal plane (FHP) on a coronal plane based on the midpoint (M). , After rotating the modified AMP on an axial plane to be parallel to a midsagittal plane (MSP) based on the midpoint (M), the modified AMP is positioned so that the center point (T) is located at the MSP. A surgical planning simulation apparatus comprising a data processing unit that moves and locates the anterior fracture segment as much as it is moved in parallel, and positions the anterior fracture segment region such that the left and right similarity is maximized based on the MSP. In claim 10,
The data processing unit is a surgical planning simulation device that adjusts the anterior bone fracture region to a final position by using a normal average value of the head measurements obtained in advance. In claim 11,
The data processing unit locates the condyle center point of the left posterior fracture region on a left vertical line passing through the uppermost point of the left articular fossa of the maxillary region, and locates the condyle center point of the right posterior fracture region to the right of the maxillary region. In a state in which the condyle center points of the left and right posterior fracture regions are positioned on the left and right vertical lines, respectively, the left and right posterior fracture regions are moved in a vertical direction while moving the left and right posterior fracture regions in a vertical direction. Surgery planning simulation device that finally positions the position where the degree of similarity is maximized. In claim 12,
The central point (T), a menton (Menton), point B (point B) projected perpendicular to the first plane, a surgical planning simulation device that is set as the midpoint of a genial tubercle. In claim 12,
The data processing unit, after finally positioning the anterior fracture region and the left and right posterior fracture region, the surgical planning simulation device for removing the area where the anterior fracture region and the left and right posterior fracture region overlap. In claim 14,
The data processing unit, based on the MSP, a surgery planning simulation device for obtaining the left and right similarity of the anterior fracture region, the similarity of the left and right posterior fracture region, and the left and right similarity of the entire mandible. In claim 15,
The similarity (SI) obtained based on the MSP is mirrored and overlapped with the left area and the right area based on the MSP, and a corresponding first volume (V ₁ ) corresponding to the left area, corresponding to the right area Surgery calculated by the equation SI = 2×V ₃ /(V ₁ +V ₂ ) using the second volume (V _{2 ) and the third volume (V 3} ) in which the left area and the right area are overlapped Planning simulation device. delete delete A computer-readable recording medium in which a program for executing any one of the above methods 1 to 9 is recorded on a computer.

Images