CN111968236A - K-D tree method-based contraposition to line fracture reduction system - Google Patents
K-D tree method-based contraposition to line fracture reduction system Download PDFInfo
- Publication number
- CN111968236A CN111968236A CN202010688900.5A CN202010688900A CN111968236A CN 111968236 A CN111968236 A CN 111968236A CN 202010688900 A CN202010688900 A CN 202010688900A CN 111968236 A CN111968236 A CN 111968236A
- Authority
- CN
- China
- Prior art keywords
- fracture
- femoral
- shaft
- points
- point
- 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
- 230000009467 reduction Effects 0.000 title claims abstract description 49
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000012634 fragment Substances 0.000 claims abstract description 112
- 230000000694 effects Effects 0.000 claims abstract description 8
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 210000003484 anatomy Anatomy 0.000 claims abstract description 5
- 238000005516 engineering process Methods 0.000 claims abstract description 3
- 208000010392 Bone Fractures Diseases 0.000 claims description 179
- 206010017076 Fracture Diseases 0.000 claims description 168
- 208000008924 Femoral Fractures Diseases 0.000 claims description 24
- 239000013598 vector Substances 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 17
- 239000013256 coordination polymer Substances 0.000 claims description 12
- 238000013519 translation Methods 0.000 claims description 10
- 210000000988 bone and bone Anatomy 0.000 claims description 8
- 210000001930 leg bone Anatomy 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 201000008181 benign familial infantile epilepsy Diseases 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000004804 winding Methods 0.000 claims description 6
- 208000006670 Multiple fractures Diseases 0.000 claims description 4
- 238000012952 Resampling Methods 0.000 claims description 4
- 210000002436 femur neck Anatomy 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 208000024779 Comminuted Fractures Diseases 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 206010016454 Femur fracture Diseases 0.000 claims description 2
- 238000004971 IR microspectroscopy Methods 0.000 claims description 2
- 210000000845 cartilage Anatomy 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000009499 grossing Methods 0.000 claims description 2
- 238000012060 immune response imaging Methods 0.000 claims description 2
- 230000003993 interaction Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 claims description 2
- 210000000689 upper leg Anatomy 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 claims description 2
- 238000012800 visualization Methods 0.000 claims description 2
- 238000009616 inductively coupled plasma Methods 0.000 claims 2
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000011946 reduction process Methods 0.000 description 6
- 239000008267 milk Substances 0.000 description 3
- 210000004080 milk Anatomy 0.000 description 3
- 235000013336 milk Nutrition 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 210000003414 extremity Anatomy 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 208000032912 Local swelling Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000740 bleeding effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000839 emulsion Substances 0.000 description 1
- 230000035876 healing Effects 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Computer Graphics (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Surgical Instruments (AREA)
Abstract
The invention discloses a K-D tree method-based alignment and alignment fracture reduction system, which combines medical science and computer technology to perform alignment and alignment iterative splicing on fracture fragments, restore the native anatomical structure of limbs and complete the primary steps of treating fracture. The system comprises a fracture fragment three-dimensional geometric modeling module, a fracture fragment information extraction module, a fracture fragment position determining and adjusting module, a registration splicing module and a display module. Obtaining the correct position before fracture reduction through a fracture fragment information extraction module and a fracture fragment position determining and adjusting module; judging whether the fragment matching degree meets the requirement or not according to the result of the matching and splicing module; and obtaining the effect after splicing the fragments through the display module. The invention guides doctors to complete the fracture fragment reduction operation, reduces the operation risk and time consumption, avoids the interference of human factors, is not limited by time, can repeatedly carry out simulation exercise, and provides technical reference for formulating the fracture reduction operation scheme.
Description
Technical Field
The invention belongs to the field of biomedical engineering, and particularly relates to a K-D tree method-based contraposition alignment fracture reduction system.
Background
Bones are the mechanism connecting the whole limbs of the human body, and after fracture, local swelling, bone exposure, bleeding, infection, deformity and obvious pain symptoms occur. After the fracture of most patients occurs, the daily life is influenced, and the psychology of the patients is also harmed. Despite the consistent interest in research relating to fracture reduction, there are various reasons why fracture reduction may not be well-established.
The fracture reduction process mainly comprises the following steps: the position of the fracture fragment is adjusted, the bone length is recovered, the anatomical structure of the bone is restored to the original geometric form, the fracture restoration result is not good, the later fracture healing process is influenced, and therefore limb deformation and abnormal growth can be caused.
Research on fracture reduction has been ongoing, and at present, fracture reduction for treating fracture fragments mainly depends on the alignment splicing performed by the clinical experience of doctors, reduces the displacement between the fracture fragments, increases the distance between fracture surfaces of the fracture fragments, but still has many disadvantages and shortcomings:
1. only by the clinical experience of doctors, the fracture reduction operation time is long, and no more accurate fracture reduction method exists;
2. the fracture reduction is only carried out with contraposition splicing, and the study on line splicing is not carried out, thus not meeting the medical standard on the fracture reduction;
3. the fracture reduction process is complex, and the condition that one-to-many and false splicing exist in the splicing of fracture fragments is not considered.
Disclosure of Invention
The invention aims to provide a K-D tree method-based alignment-alignment fracture reduction system which is used for predicting a complex fracture reduction process so as to obtain an optimal fracture reduction operation planning scheme.
The purpose of the invention is realized by the following technical scheme: a system for contraposition and alignment fracture reduction based on a K-D tree method is characterized by comprising a fracture fragment three-dimensional geometric modeling module, a fracture fragment information extraction module, a fracture fragment position determining and adjusting module, a matching and splicing module and a display module;
the fracture fragment three-dimensional geometric modeling module is used for importing two-dimensional tomographic image data obtained by CT medical scanning into medical image processing software in a DICOM format, and establishing a three-dimensional grid geometric model representing fracture fragments after a series of image processing so as to provide a useful model for fracture reduction;
the fracture fragment information extraction module is used for extracting the geometrical information such as points, lines, surfaces and the like which can represent femoral fractures and are required by fracture reduction from the obtained fracture fragment three-dimensional grid geometrical model, and establishing an information structure body of the femoral fracture fragments to provide useful fracture fragment information for fracture reduction;
the fracture fragment position determining and adjusting module is used for positioning and evaluating the displaced femoral fracture fragments to provide accurate positions for fracture reduction;
the matching and splicing module is used for aligning femoral fracture fragments, and restoring the fracture fragments to the original anatomical structure position by finding and determining the corresponding relation between different fracture fragments, so that the fracture reduction stability is ensured;
the display module is used for observing and displaying the comminuted fracture alignment and line reduction effect and realizing the interactivity between a user and the fracture reduction system;
compared with the prior art, the invention has the beneficial effects that:
1. the K-D tree method-based alignment and alignment fracture reduction system is developed on the basis of a dialog box form under a Windows platform, the simulation of the fracture reduction process is realized through autonomous programming, the system is simple and easy to operate, and the system flow can be rapidly mastered in a short time;
2. the pair wires are spliced and added into a fracture reduction simulation system, so that the fracture reduction process can be reflected more truly, the fracture reduction result meets the medical clinical actual condition, and the splicing result is more accurate;
3. the fracture reduction process is divided into primary splicing and accurate splicing, compared with a primary splicing process, false splicing can be deleted in the primary splicing process, and the condition that fracture fragments are one-to-many before accurate splicing is avoided;
4. by constructing the fracture reduction simulation method, guidance can be provided for a doctor to designate an optimal operation scheme, so that the operation success rate is improved, the fracture reduction quality is improved, and the condition of fracture reduction deformity is reduced;
5. by constructing the fracture reduction simulation method, repeated experiments can be carried out on the established simulation model for many times, so that real biological experiments are reduced, time and resources are saved, efficiency is improved, and humanitarian disputes are avoided.
In summary, the simulation method of the present invention overcomes the disadvantages and shortcomings of the prior art.
Drawings
FIG. 1 is a schematic diagram showing the relationship between modules of the fracture reduction simulation system of the present invention;
FIG. 2 is a flow chart of a three-dimensional geometric modeling module for establishing fracture fragments;
FIG. 3 is a flow chart of a module for establishing and extracting fracture fragment information;
FIG. 4 is a flow chart of a module for establishing the determination and adjustment of the fracture fragment position;
FIG. 5 is a flow chart of building a matched tiling module;
FIG. 6 is a flow chart for establishing a preliminary splice;
FIG. 7 is a flow chart for establishing an accurate splice;
FIG. 8 is a flowchart illustrating the process of creating a K _ D tree;
fig. 9 is a flow chart of the alignment-to-alignment stitching algorithm.
Detailed Description
The first embodiment is as follows: as shown in fig. 1, the system for reduction of alignment-alignment fracture based on the K-D tree according to the present embodiment includes: the fracture fragment three-dimensional geometric modeling module, the fracture fragment information extraction module, the fracture fragment position determining and adjusting module, the matching and splicing module and the display module are sequentially connected;
the three-dimensional geometric modeling module of the fracture fragments is used for importing two-dimensional tomographic image data obtained by CT medical scanning into medical image processing software in a DICOM format, and establishing a three-dimensional grid geometric model representing the fracture fragments after a series of image processing, so that a useful model is provided for fracture reduction.
The fracture fragment information extracting module is used for extracting geometrical information such as points, lines, surfaces and the like which can represent femoral fractures and are required by fracture reduction from the obtained fracture fragment three-dimensional grid geometrical model, establishing an information structure body of the femoral fracture fragments and providing useful fracture fragment information for fracture reduction.
The module for determining and adjusting the position of the fracture fragment is used for positioning and evaluating the displaced femoral fracture fragment to provide an accurate position for fracture reduction.
The matching and splicing module is used for aligning femoral fracture fragments, and restoring the fracture fragments to the original anatomical structure position by finding and determining the corresponding relation between different fracture fragments, so that the stability of fracture restoration is ensured.
The display module is used for observing and displaying the comminuted fracture alignment and line resetting effect, and realizes the interactivity between the user and the fracture resetting system.
The second embodiment is as follows: as shown in fig. 2, in the present embodiment, the specific process of the three-dimensional geometric modeling module for fracture fragments to achieve its function is as follows:
cartilage tissues and bones are separated through a proper threshold value, the bone tissues are marked or identified by a segmentation method based on region growth, fracture fragments are extracted, then a three-dimensional grid geometric model of the fracture fragments required for reduction is obtained through the grid division and three-dimensional reconstruction processes of an MC algorithm, and optimization technologies such as smoothing, simplification or regridding and the like are carried out on the obtained three-dimensional grid geometric model, so that the characteristics of the model are improved, and the visualization effect and the manageability of the model are further improved.
Other components and connection relationships of the present embodiment are the same as those of the first embodiment.
The third concrete implementation mode: as shown in fig. 3, in the present embodiment, the specific process of the module for extracting fracture fragment information to implement its function is as follows:
a. extracting the geometric information of the three-dimensional broken bone model, and establishing an information structure body of the femoral fracture fragments, wherein the information structure body is expressed as follows:
in the formula, fracture fragment information B ═ { B ═ B1,B2,…Bn},BiThe ith fracture fragment information in the fracture fragment information structure body; gMGeometric elements of a three-dimensional model of a fracture fragment, including the geometric shape GFSide GEPoint GP;IMFor fracture patient information, including patient name INSex ISAge IALeft and right leg bones IL;
b. Traversing the established information structure BFIS of the femoral fracture fragment to obtain a femoral fracture fragment model;
obtaining the total bounding box X of all broken bones and obtaining the point set at the lower part of the bounding box XFinding corresponding point set from BFISWherein the content of the first and second substances,the point set of the ith broken bone in the fracture fragment information structure body is shown; b isiObtaining the ith fracture fragment information in the fracture fragment information structure body to obtain a femoral shaft fragment model Mshaft;
Acquiring the total bounding box X of all broken bones according to the left and right leg bones I in the fracture fragment information BLInformation, if ILIs the left leg bone, the point set around the minimum value of the bounding box X is obtainedIf IRIs the right leg bone, the point set around the maximum value of the bounding box X is obtainedFinding corresponding point set from BFISWherein the content of the first and second substances,the point set of the ith fracture fragment in the fracture fragment information structure body; b isiObtaining the information of the ith fracture fragment in the information structure of the fracture fragments to obtain a femoral head fracture model Mhead;
c. Extracting the axis of the femur from the obtained femur fracture fragment model;
according to the femoral head bone-breaking model MheadResampling the femoral head model to obtain a three-dimensional coordinate point set S of the femoral headhead;
According to the obtained point set SheadThe sphere is fitted by the least square method, and a function H (x) of the sum of squares of the errors of the least square method is constructed as shown below0,y0,z0R) to find the center coordinates of the fitting sphere and the sphere radius, can be expressed as follows:
in the formula, xi、yi、ziSet of points SheadCoordinate value of the ith point; n is a set of points SheadThe number of middle coordinate points; x is the number of0、y0、z0To the centre P of the fitted sphere0The coordinates of (a); r is the radius of the fitting sphere;
set of points SheadCenter of sphere P of median fitting sphere0(x0,y0,z0) The far points form a femur neck three-dimensional coordinate point set SneckAccording to the set of points SneckThe function F (x ') was constructed as follows'0,y′0,z′0) Represents:
thereby obtaining the coordinate of the central point of the femoral neck as PN(x′0,y′0,z′0)xj、yj、zjSet of points SneckCoordinate value of j point; n is a set of points SneckThe number of middle coordinate points;
connection point PN(x′0,y′0,z′0) And P0(x0,y0,z0) Obtaining the femoral head axis HL;
broken bone model M according to femoral shaftshaft,Resampling the femoral shaft model to obtain a femoral shaft three-dimensional coordinate point set Sshaft;
According to the obtained point set SshaftSet points SshaftThe N points are sorted according to the sequence of Z coordinates from small to large, the points within the range of 30-60% of the maximum value of the Z coordinates are extracted, and a point set S 'is constructed'shaft(ii) a Will click set S'shaftThe points in the formula are sorted from small to large according to Z coordinates, and are sequentially divided into 15 parts, wherein each part has n coordinate points, and the formula is as follows:
calculating the center point P of each coordinate point seti(xi,yi,zi) And will point Pi(xi,yi,zi) Put in set SshaftFitting a space straight line according to a least square method to obtain a femoral shaft axis line emulsion;
d. extracting a femoral fracture surface point set;
wherein x isj、xkPoints of the information structure BFIS of femoral fracture fragments; x is the number ofkIs xjPoints in the left and right σ ranges; n is a radical ofσIs f (x)j,xk) The number of medium non-zero points; f (x)j,xk) Is a point xj、xkThe angle between the normal vectors of (a); n isj、nkIs xj、xkThe normal vector of (a); σ is xkAnd xjThe distance between the two, usually 2-4 mm;
function W (x)j) Record the point xjThe included angle of other points must be larger than the threshold value t2(ii) a If point xjSatisfy the requirement ofThen point xjTo a point on the fracture surface, otherwise not on the fracture surface.
Other components and connection relationships of this embodiment are the same as those of the first to second embodiments.
The fourth concrete implementation mode: as shown in fig. 4, in this embodiment, the specific process of determining and adjusting the position of the fracture fragment by the module to achieve its function is as follows:
e. determining an axial plane AP, a coronal plane CP and a sagittal plane SP of a femoral shaft fracture surface according to the obtained femoral shaft axial lead, and adjusting the positions of the coronal plane CP and the axial plane AP;
the axial line milk of the femoral shaft is vertical to the axial plane AP; the coronal plane CP and the sagittal plane SP are perpendicular to the axial plane AP and both pass through the axis milk of the femoral shaft; after determining the positions of the coronal plane CP and the sagittal plane SP, manually adjusting the coronal plane CP and the sagittal plane SP to the appropriate positions according to medical common knowledge;
f. calculating and adjusting the neck trunk angle and the anteversion angle of the femoral shaft;
the neck angle θ, can be expressed as follows:
in the formula, N1Is the vector of the femoral shaft axis milk; n is a radical of2Is the vector of the femoral head axis HL; theta is an included angle between the two axial lines;
in the formula, N3Is a projection vector of a femoral head axis HL on an axial plane AP; n is a radical of4Is the intersection vector of the coronal plane CP and the axial plane AP;is the included angle between the two vectors;
if the neck angle theta is in the range of 120-135 DEG and the anteversion angleThe range of 12-15 degrees meets the requirements of medical clinical specifications; otherwise, readjust the angle theta of the shaft and the anteversion angleUntil the requirements of medical clinical specifications are met.
Other components and connection relationships of this embodiment are the same as those of the first to third embodiments.
The fifth concrete implementation mode: as shown in fig. 5, in this embodiment, the specific process of the matching and splicing module to implement its function is as follows:
g. completing the initial splicing by using a central point coincidence method, as shown in FIG. 6;
in order to perform center coincidence on the femoral head and the fracture surface point set of the femoral shaft, reduce the translation dislocation between the two point sets and provide initial values of translation and rotation parameters for accurate splicing;
calculating the centroid of the set of femoral head and femoral shaft fracture surface points can be expressed as follows:
in the formula (I), the compound is shown in the specification,fracture surface point sets of the femoral head and the femoral shaft respectively; n is the number of corresponding fracture surface point sets;respectively as the centroids of the femoral head and femoral shaft fracture surface point sets;
the offset vector of the two centroids is recorded and can be expressed as follows:
in the formula (I), the compound is shown in the specification,is a translation vector;respectively as the centroids of the femoral head and femoral shaft fracture surface point sets;
the coordinates of each point in the femoral head fracture surface point set are translated by taking the femoral shaft as a reference object, and can be expressed as follows:
in the formula (I), the compound is shown in the specification,is a fracture surface point set of femoral head;is a translation vector;the coordinates of the femoral head fracture surface point set after the center of gravity is formed;
h. an ICP algorithm is used for completing accurate splicing of fracture fragments of the femoral head and the femoral shaft, and the splicing is shown in figure 7;
in order to improve the accuracy of splicing, the false splicing needs to be deleted through the processes of rotation, movement and the like after the initial splicing, and the shape of the fracture fragment is finally adjusted and aligned;
using the femoral shaft as a reference object, searching by using a K _ D tree to find the corresponding closest point of the points in the femoral shaft fracture surface point set in the femoral head fracture surface point set, as shown in FIG. 8;
the Euclidean distance can be expressed as follows:
in the formula, d is the Euclidean distance between the corresponding point vectors of the femoral shaft and the femoral head; (x)1,,y1,z1),(x2,,y2,z2) Coordinates of corresponding points of the femoral head and the femoral shaft are respectively;
solving a rotation matrix and a translation matrix;
the coordinate system of the femoral head model is OcXcYcZcThe transformed coordinate system is expressed asWherein d is the distance between the femoral head model and the position after restoration; omega is femoral head model winding XwThe angle of rotation; theta is the femoral head model winding YwThe angle of rotation;for the femoral head model winding ZwThe angle of rotation; the transformation matrix R of the femoral head is three matricesThe product of (c) can be expressed as follows:
as shown in fig. 9, verifying whether the distance between the two sets of fracture surface points of the corresponding femoral head and femoral shaft is smaller than 4mm, if so, adjusting the three-dimensional spatial position of the femoral head according to the transformation matrix R of the femoral head, obtaining the information of the three-dimensional fractured model according to the adjusted spatial position coordinates of the femoral head, replacing the information structure BFIS of the femoral fracture fragment, and finishing the precise splicing; if not, returning to the step of searching the closest point by using the K _ D tree.
Other components and connection relationships of the present embodiment are the same as those of the first to fourth embodiments.
The sixth specific implementation mode: in this embodiment, the specific process of the display module to realize its function is as follows:
and connecting the output display equipment of the computer to the computer for visual interaction, so that the effect of fracture reduction can be previewed.
Other components and connection relationships of the present embodiment are the same as those of one of the first to fifth embodiments.
Claims (1)
1. An alignment-to-alignment fracture reduction system based on a K-D tree method is characterized by comprising:
the fracture fragment three-dimensional geometric modeling module (1), the fracture fragment information extraction module (2), the fracture fragment position determining and adjusting module (3), the matching and splicing module (4) and the display module (5);
the three-dimensional geometric modeling module (1) of fracture fragments is used for importing two-dimensional tomographic image data obtained by CT medical scanning into medical image processing software in a DICOM format, and establishing a three-dimensional grid geometric model representing the fracture fragments after a series of image processing so as to provide a useful model for fracture reduction;
the three-dimensional geometric modeling module (1) for the fracture fragments has the specific process of realizing the functions as follows:
cartilage tissues and bones are separated through a proper threshold value, the bone tissues are marked or identified by a segmentation method based on region growth, so that fracture fragments are extracted, then a three-dimensional grid geometric model of the fracture fragments required for reduction is obtained by adopting an MC algorithm to carry out grid division and three-dimensional reconstruction processes, and optimization technologies such as smoothing, simplification or regridding and the like are carried out on the obtained three-dimensional grid geometric model, so that the characteristics of the model are improved, and the visualization effect and the manageability of the model are further improved;
the fracture fragment information extraction module (2) is used for extracting the geometrical information such as points, lines, surfaces and the like which can represent femoral fractures and are required by fracture reduction from the obtained fracture fragment three-dimensional grid geometrical model, and establishing an information structure body of the femoral fracture fragments to provide useful fracture fragment information for fracture reduction;
the module (2) for extracting the fracture fragment information realizes the specific process of the function as follows:
a. extracting the geometric information of the three-dimensional broken bone model, and establishing an information structure body of the femoral fracture fragments, wherein the information structure body is expressed as follows:
in the formula, fracture fragment information B ═ { B ═ B1,B2,…Bn},BiThe ith fracture fragment information in the fracture fragment information structure body; gMGeometric elements of a three-dimensional model of a fracture fragment, including the geometric shape GFSide GEPoint GP;IMFor fracture patient information, including patient name INSex ISAge IALeft and right leg bones IL;
b. Traversing the established information structure BFIS of the femoral fracture fragment to obtain a femoral fracture fragment model;
obtaining the total bounding box X of all broken bones and obtaining the point set at the lower part of the bounding box XFinding corresponding point set from BFISWherein the content of the first and second substances,the point set of the ith broken bone in the fracture fragment information structure body is shown; b isiObtaining the ith fracture fragment information in the fracture fragment information structure body to obtain a femoral shaft fragment model Mshaft;
Acquiring the total bounding box X of all broken bones according to the left and right leg bones I in the fracture fragment information BLInformation, if ILIs the left leg bone, the point set around the minimum value of the bounding box X is obtainedIf IRIs the right leg bone, the point set around the maximum value of the bounding box X is obtainedFinding corresponding point set from BFISWherein the content of the first and second substances,the point set of the ith fracture fragment in the fracture fragment information structure body; b isiObtaining the information of the ith fracture fragment in the information structure of the fracture fragments to obtain a femoral head fracture model Mhead;
c. Extracting the axis of the femur from the obtained femur fracture fragment model;
according to the femoral head bone-breaking model MheadResampling the femoral head model to obtain a three-dimensional coordinate point set S of the femoral headhead;
According to the obtained point set SheadThe sphere is fitted by the least square method, and a function H (x) of the sum of squares of the errors of the least square method is constructed as shown below0,y0,z0R) to find the center coordinates of the fitting sphere and the sphere radius, can be expressed as follows:
in the formula, xi、yi、ziSet of points SheadCoordinate value of the ith point; n is a set of points SheadThe number of middle coordinate points; x is the number of0、y0、z0To the centre P of the fitted sphere0The coordinates of (a); r is the radius of the fitting sphere;
set of points SheadCenter of sphere P of median fitting sphere0(x0,y0,z0) The far points form a femur neck three-dimensional coordinate point set SneckAccording to the set of points SneckThe function F (x ') was constructed as follows'0,y′0,z′0) Represents:
thereby obtaining the coordinate of the central point of the femoral neck as PN(x′0,y′0,z′0)xj、yj、zjSet of points SneckCoordinate value of j point; n is a set of points SneckThe number of middle coordinate points;
connection point PN(x′0,y′0,z′0) And P0(x0,y0,z0) Obtaining the femoral head axis HL;
broken bone model M according to femoral shaftshaftResampling the femoral shaft model to obtain a femoral shaft three-dimensional coordinate point set Sshaft;
According to the obtained point set SshaftSet points SshaftThe N points are sorted according to the sequence of Z coordinates from small to large, the points within the range of 30-60% of the maximum value of the Z coordinates are extracted, and a point set S 'is constructed'shaft(ii) a Will click set S'shaftThe points in the formula are sorted from small to large according to Z coordinates, and are sequentially divided into 15 parts, wherein each part has n coordinate points, and the formula is as follows:
calculating the center point P of each coordinate point seti(xi,yi,zi) And will point Pi(xi,yi,zi) Put in set SshaftFitting a space straight line according to a least square method to obtain a femoral shaft axis SL;
d. extracting a femoral fracture surface point set;
wherein x isj、xkPoints of the information structure BFIS of femoral fracture fragments; x is the number ofkIs xjPoints in the left and right σ ranges; n is a radical ofσIs f (x)j,xk) The number of medium non-zero points; f (x)j,xk) Is a point xj、xkThe angle between the normal vectors of (a); n isj、nkIs xj、xkThe normal vector of (a); σ is xkAnd xjThe distance between the two, usually 2-4 mm;
function W (x)j) Record the point xjThe included angle of other points must be larger than the threshold value t2(ii) a If point xjSatisfy the requirement ofThen point xjPoints belonging to the fracture surface, otherwise not points on the fracture surface;
the fracture fragment position determining and adjusting module (3) is used for positioning and evaluating the displaced femoral fracture fragments to provide accurate positions for fracture reduction;
the module (3) for determining and adjusting the position of the fracture fragment has the specific process of realizing the functions as follows:
e. determining an axial plane AP, a coronal plane CP and a sagittal plane SP of a femoral shaft fracture surface according to the obtained femoral shaft axial lead, and adjusting the positions of the coronal plane CP and the axial plane AP;
the axial line SL of the femoral shaft is vertical to the axial plane AP; the coronal plane CP and the sagittal plane SP are perpendicular to the axial plane AP and both pass through an axis SL of the femoral shaft; after determining the positions of the coronal plane CP and the sagittal plane SP, manually adjusting the coronal plane CP and the sagittal plane SP to the appropriate positions according to medical common knowledge;
f. calculating and adjusting the neck trunk angle and the anteversion angle of the femoral shaft;
the neck angle θ, can be expressed as follows:
in the formula, N1Is the vector of the femoral shaft axis SL; n is a radical of2Is the vector of the femoral head axis HL; theta isThe included angle between the two axial lines;
in the formula, N3Is a projection vector of a femoral head axis HL on an axial plane AP; n is a radical of4Is the intersection vector of the coronal plane CP and the axial plane AP;is the included angle between the two vectors;
if the neck angle theta is in the range of 120-135 DEG and the anteversion angleThe range of 12-15 degrees meets the requirements of medical clinical specifications; otherwise, readjust the angle theta of the shaft and the anteversion angleUntil the requirements of medical clinical specifications are met;
the matching and splicing module (4) is used for aligning femoral fracture fragments, and restoring the fracture fragments to the original anatomical structure position by finding and determining the corresponding relation among different fracture fragments, so that the stability of fracture restoration is ensured;
the matching splicing module (4) realizes the functions by the following specific processes:
g. completing primary splicing by using a central point coincidence method;
in order to perform center coincidence on the femoral head and the fracture surface point set of the femoral shaft, reduce the translation dislocation between the two point sets and provide initial values of translation and rotation parameters for accurate splicing;
calculating the centroid of the set of femoral head and femoral shaft fracture surface points can be expressed as follows:
in the formula (I), the compound is shown in the specification,fracture surface point sets of the femoral head and the femoral shaft respectively; n is the number of corresponding fracture surface point sets;respectively as the centroids of the femoral head and femoral shaft fracture surface point sets;
the offset vector of the two centroids is recorded and can be expressed as follows:
in the formula (I), the compound is shown in the specification,is a translation vector;respectively as the centroids of the femoral head and femoral shaft fracture surface point sets;
the coordinates of each point in the femoral head fracture surface point set are translated by taking the femoral shaft as a reference object, and can be expressed as follows:
in the formula (I), the compound is shown in the specification,is a fracture surface point set of femoral head;is a translation vector;the coordinates of the femoral head fracture surface point set after the center of gravity is formed;
h. accurate splicing of fracture fragments of the femoral head and the femoral shaft is completed by applying an ICP (inductively coupled plasma) algorithm;
in order to improve the accuracy of splicing, the false splicing needs to be deleted through the processes of rotation, movement and the like after the initial splicing, and the shape of the fracture fragment is finally adjusted and aligned;
using the femoral shaft as a reference object, and searching and finding out the corresponding nearest point of the point in the femoral shaft fracture surface point set in the femoral head fracture surface point set by using a K _ D tree;
the Euclidean distance can be expressed as follows:
in the formula, d is the Euclidean distance between the corresponding point vectors of the femoral shaft and the femoral head; (x)1,,y1,z1),(x2,,y2,z2) Coordinates of corresponding points of the femoral head and the femoral shaft are respectively;
solving a rotation matrix and a translation matrix;
the coordinate system of the femoral head model is OcXcYcZcThe transformed coordinate system is expressed asWherein d is the distance between the femoral head model and the position after restoration; omega is femoral head model winding XwThe angle of rotation; theta is the femoral head model winding YwThe angle of rotation;for the femoral head model winding ZwThe angle of rotation; the transformation matrix R of the femoral head is three matrices Rx(ω),Ry(θ),The product of (c) can be expressed as follows:
verifying whether the distance between the corresponding femoral head and femoral shaft fracture surface point sets is smaller than 4mm, if so, adjusting the three-dimensional space position of the femoral head according to the femoral head transformation matrix R, obtaining the information of a three-dimensional bone fracture model according to the adjusted space position coordinates of the femoral head, replacing an information structure body BFIS of femoral fracture fragments, and finishing accurate splicing; if not, returning to the step of searching the closest point by using the K _ D tree;
the display module (5) is used for observing and displaying the comminuted fracture alignment and line reduction effect and realizing the interactivity between a user and the fracture reduction system;
the display module (5) realizes the functions of the display module by the following specific processes:
and connecting the output display equipment of the computer to the computer for visual interaction, so that the effect of fracture reduction can be previewed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010688900.5A CN111968236A (en) | 2020-07-17 | 2020-07-17 | K-D tree method-based contraposition to line fracture reduction system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010688900.5A CN111968236A (en) | 2020-07-17 | 2020-07-17 | K-D tree method-based contraposition to line fracture reduction system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111968236A true CN111968236A (en) | 2020-11-20 |
Family
ID=73361909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010688900.5A Pending CN111968236A (en) | 2020-07-17 | 2020-07-17 | K-D tree method-based contraposition to line fracture reduction system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111968236A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113450435A (en) * | 2021-06-25 | 2021-09-28 | 网易(杭州)网络有限公司 | Vertex animation processing method and device |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102961176A (en) * | 2012-12-13 | 2013-03-13 | 哈尔滨理工大学 | Method for selecting optimal inner fixing modes for femur and neck bone fractures |
US20130114866A1 (en) * | 2011-11-08 | 2013-05-09 | Mako Surgical Corporation | Computer-Aided Planning with Dual Alpha Angles in Femoral Acetabular Impingement Surgery |
CN103632371A (en) * | 2013-12-06 | 2014-03-12 | 河海大学常州校区 | Compatibility mesh segmentation based skeleton parameter computation method |
WO2015142291A1 (en) * | 2014-03-20 | 2015-09-24 | National University Of Singapore | Computer-aided planning of craniomaxillofacial and orthopedic surgery |
US20180168731A1 (en) * | 2015-05-29 | 2018-06-21 | The Penn State Research Foundation | Individualized preoperative planning system and method |
US20180325599A1 (en) * | 2017-05-12 | 2018-11-15 | Coreline Soft Co., Ltd. | Computer assistance system and method for image-guided reduction of fracture |
CN109528305A (en) * | 2018-10-12 | 2019-03-29 | 河海大学常州校区 | Computer assisted intertrochanter fracture of femur contraposition is to the semi-automatic repositioning method of line |
US20190298452A1 (en) * | 2018-03-27 | 2019-10-03 | Virginia Commonwealth University | Automated hip analysis methods and devices |
-
2020
- 2020-07-17 CN CN202010688900.5A patent/CN111968236A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130114866A1 (en) * | 2011-11-08 | 2013-05-09 | Mako Surgical Corporation | Computer-Aided Planning with Dual Alpha Angles in Femoral Acetabular Impingement Surgery |
CN102961176A (en) * | 2012-12-13 | 2013-03-13 | 哈尔滨理工大学 | Method for selecting optimal inner fixing modes for femur and neck bone fractures |
CN103632371A (en) * | 2013-12-06 | 2014-03-12 | 河海大学常州校区 | Compatibility mesh segmentation based skeleton parameter computation method |
WO2015142291A1 (en) * | 2014-03-20 | 2015-09-24 | National University Of Singapore | Computer-aided planning of craniomaxillofacial and orthopedic surgery |
US20180168731A1 (en) * | 2015-05-29 | 2018-06-21 | The Penn State Research Foundation | Individualized preoperative planning system and method |
US20180325599A1 (en) * | 2017-05-12 | 2018-11-15 | Coreline Soft Co., Ltd. | Computer assistance system and method for image-guided reduction of fracture |
US20190298452A1 (en) * | 2018-03-27 | 2019-10-03 | Virginia Commonwealth University | Automated hip analysis methods and devices |
CN109528305A (en) * | 2018-10-12 | 2019-03-29 | 河海大学常州校区 | Computer assisted intertrochanter fracture of femur contraposition is to the semi-automatic repositioning method of line |
Non-Patent Citations (2)
Title |
---|
王虹;丁焕文;刘宝;涂强;刘辉亮;沈健坚;: "计算机辅助逆行交锁髓内钉治疗股骨髁间髁上粉碎性骨折", 中国骨与关节损伤杂志, no. 08 * |
邓子越;蒋俊锋;陈正鸣;黄瑞;: "对位对线的粗隆间骨折计算机辅助复位方法", 计算机辅助设计与图形学学报, no. 05 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113450435A (en) * | 2021-06-25 | 2021-09-28 | 网易(杭州)网络有限公司 | Vertex animation processing method and device |
CN113450435B (en) * | 2021-06-25 | 2023-05-26 | 网易(杭州)网络有限公司 | Vertex animation processing method and device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10937250B2 (en) | Methods of reconstructing skulls | |
Hacihaliloglu et al. | Local phase tensor features for 3-D ultrasound to statistical shape+ pose spine model registration | |
Deng et al. | A novel skull registration based on global and local deformations for craniofacial reconstruction | |
JP2003530177A (en) | Computer aided bone lengthening | |
CN110891488B (en) | Sagittal rotation determination | |
US20220071708A1 (en) | Patient positioning using a skeleton model | |
CN112509022A (en) | Non-calibration object registration method for preoperative three-dimensional image and intraoperative perspective image | |
CN109528305A (en) | Computer assisted intertrochanter fracture of femur contraposition is to the semi-automatic repositioning method of line | |
CN112164443B (en) | Interactive simulation fracture reduction method | |
Irwansyah et al. | Algorithm for segmentation and reduction of fractured bones in computer-aided preoperative surgery | |
Liu et al. | Virtual plate pre-bending for the long bone fracture based on axis pre-alignment | |
CN109350059B (en) | Combined steering engine and landmark engine for elbow auto-alignment | |
CN111968236A (en) | K-D tree method-based contraposition to line fracture reduction system | |
Morooka et al. | A survey on statistical modeling and machine learning approaches to computer assisted medical intervention: Intraoperative anatomy modeling and optimization of interventional procedures | |
WO2022133442A1 (en) | Systems and methods for generating a three-dimensional model of a joint from two-dimensional images | |
US9639977B2 (en) | Processing model data representing a biological body section having a defect portion | |
Li et al. | A fully automatic surgical registration method for percutaneous abdominal puncture surgical navigation | |
Berti et al. | Medical simulation services via the grid | |
CN114283179A (en) | Real-time fracture far-near end space pose acquisition and registration system based on ultrasonic images | |
Xie et al. | A small-scaled intraoperative 3d visualization navigation system for femoral head repair surgery | |
US20200197096A1 (en) | Pre-operative planning for reorientation surgery: surface-model-free approach using simulated x-rays | |
Irwansyah et al. | Repositioning bone fragments using registration of paired-points and assisted-constraints in virtual bone reduction surgery | |
CN112085698A (en) | Method and device for automatically analyzing left and right breast ultrasonic images | |
Wang et al. | Virtual reassembly of fractured bones for orthopedic surgery | |
Zhang et al. | Development and model form assessment of an automatic subject-specific vertebra reconstruction method |
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 |