CN102968565A - Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm - Google Patents
Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm Download PDFInfo
- Publication number
- CN102968565A CN102968565A CN2012104870245A CN201210487024A CN102968565A CN 102968565 A CN102968565 A CN 102968565A CN 2012104870245 A CN2012104870245 A CN 2012104870245A CN 201210487024 A CN201210487024 A CN 201210487024A CN 102968565 A CN102968565 A CN 102968565A
- Authority
- CN
- China
- Prior art keywords
- cruciate ligament
- anterior cruciate
- software
- femur
- section
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
The invention discloses a method for obtaining a central axial line of an anterior cruciate ligament based on a linear fitting algorithm. The method comprises the steps of: first establishing a complete three-dimensional solid digital model of two knee joints of an STL (standard template library) format, leading in finite element analysis software for calculation to obtain a cross section of the ACL (anterior cruciate ligament) and the coordinates of the center of the cross section, fitting to obtain a center curve of the anterior cruciate ligament, solving the equation, calculating the coordinates of the center of a femur and tibia attachment point, and calculating the specific value of a planar angle by using a trigonometric function and the coordinates of the geometric center of a femur and tibia attachment surface. Accurate positioning of femur and tibia tunnels of the ACL can be realized. Preoperative planning, surgical simulation exercise, surgery teaching, guidance of actual surgical operations, post-operative rehabilitation and the like can be performed by adopting information of a virtual environment. The surgical accuracy can be improved, the approximately normal functions of the ACL can be reestablished, and the method is an effective new path for repair and reestablishment of the efficacy of the ACL.
Description
Technical field
The invention belongs to the computer-aided medical science technical field, the preoperative planning virtual technology while relating to Healing in Anterior Cruciate Ligament Reconstruction, be specifically related to a kind ofly calculate the method that obtains the anterior cruciate ligament axis stretching under the knee state by linear fit.
Background technology
Anterior cruciate ligament (Anterior cruciate ligament, ACL) damage is clinical comparatively common and have a strong impact on the damage of knee Stability, its incidence of disease is higher, all occupies suitable high specific gravity in athletic contest motion, traffic accident injury and military training.In China, along with extensively carrying out and the increase of traffic hazard of nationwide fitness programs, cross ligament damage increases day by day, has become one of common frdquently encountered disease clinically.U.S.'s cross ligament damage incidence of disease is annual 38/10 ten thousand people at population at present, at football and skiing, is respectively annual 60/10 ten thousand people and annual 70/10 ten thousand people.China there is no the report of whole man group's incidence of disease, but the investigation that has the author to carry out China active service training sportsman is found, the overall incidence of disease 0.47% of anterior cruciate ligament, in China professional athlete cross ligament damage person, women's the rate of injury is 0.71%, the male sex is 0.29%, the women is the male sex 2.37 times.
Cross ligament damage is difficult to self-healing, and malpractice can cause instability of knee, causes knee joint secondary meniscus, the too early regression of cartilage damage and joint, has a strong impact on function of joint and patient's live and work quality, must repair in time.Therefore the Healing in Anterior Cruciate Ligament Reconstruction operation is one of orthopaedics and the modal operation of sports medical science.Although under arthroscope, Healing in Anterior Cruciate Ligament Reconstruction operation has very high success ratio, yet still have 10%～15% mortality in bibliographical information and up to 10%～20% revision rate.The researchs such as Navy Medicine center Khafayan show, there is 10%～40% Healing in Anterior Cruciate Ligament Reconstruction tunnel location incorrect, and the accuracy of anterior cruciate ligament osseous tunnel position is the key of operation, it will directly affect postoperative effect, be to cause the knee joint unstability, Range of motion of knee joint changes, gonarthromeningitis, Traumatic arthritis, the key factor that ligamentaum cruciatum clashes into.The Kodali of Northwestern Univ USA etc. thinks and rebuilds by the osseous tunnel positional fault again and cause up to 40% anterior cruciate ligament.Therefore, how one of study hotspot that accurate location is Healing in Anterior Cruciate Ligament Reconstruction and difficult problem are carried out in anterior cruciate ligament osseous tunnel position.
Summary of the invention
The object of the invention is to for the deficiencies in the prior art, by three-dimensional digital model the computer ancillary technique of setting up anterior cruciate ligament of knee joint, computing is carried out in the axis of stretching the anterior cruciate ligament under the knee state, linear fit calculates the axis of anterior cruciate ligament, obtain the accurate location of Healing in Anterior Cruciate Ligament Reconstruction osseous tunnel by this axis, further instruct reconstruction of anterior cruciate ligament by arthroscopy, for the accurate location of anterior cruciate ligament osseous tunnel, the success ratio that improves Healing in Anterior Cruciate Ligament Reconstruction lays the foundation.
Purpose of the present invention is achieved through the following technical solutions.
A kind of method that obtains the anterior cruciate ligament axis based on the linear fit algorithm comprises the following steps:
(1) software preparation: the Mimics of medical finite elements simulation software, reverse engineering software and finite element analysis software ANSYS;
(2) foundation of three-dimensional digital model
To import the Mimics of medical finite elements simulation software 10.01 by the software input port on computers according to the knee joint MRI fault image data of dicom standard storage, after locating picture, organizing picture and interior difference to process, automatically generate suitable threshold value; Software generates the profile of the structure of every aspect automatically, carries out at every tomographic image by hand that marginality is cut apart, selective editing, filling-up hole process, and removes artifact and miscellaneous data; Generate shin bone, femur, kneecap, fibula, anterior cruciate ligament, posterior cruciate ligament of knee and meniscus how much facial mask files separately by region growing; Finally set up the complete double knee joint 3D solid digital model that comprises DF, proximal tibia, anterior cruciate ligament, meniscus, kneecap and fibula, with the STL form, derive and preserve;
(3) importing of anterior cruciate ligament three-dimensional digital model:
The anterior cruciate ligament model part of STL form is imported in ANSYS13.0 finite element analysis software working interface, and show separately anterior cruciate ligament;
(4) the anterior cruciate ligament cross section obtains and centre coordinate numerical value
Obtain anterior cruciate ligament cross section and this kernel of section coordinate in the ANSYS13.0 finite element analysis software, and be parallel to the anterior cruciate ligament face and up or down the anterior cruciate ligament model carried out to cross section respectively, and obtain its center and coordinate values thereof with same method;
(5) matching of anterior cruciate ligament kernel of section curve and curve equation thereof solves
By the coordinate values of obtained anterior cruciate ligament kernel of section, simulate the center curve of anterior cruciate ligament and solve its equation after importing ANASYS13.0 software, equation is:
Z=0.2412+0.7879x-6.5738y+?2.9976x
2+?1.8441y
2+2.5504xy
(6) femur and Tibial attachment centre coordinate value tries to achieve
In the ANSYS13.0 finite element analysis software, try to achieve the spatial value of anterior cruciate ligament axis and femur intersection point by the method for high linear algebra; In like manner, draw the spatial value of anterior cruciate ligament axis and shin bone intersection point;
(7) the anterior cruciate ligament central axis is in the acquisition of each Plane Angle of coordinate system
By the method for trigonometric function, apply the concrete value of femur and shin bone attachment surface Geometric center coordinates numerical evaluation out-of-plane angle again.
With respect to prior art, the present invention has following useful result:
The present invention has realized the accurate location in anterior cruciate ligament tunnel on femur and shin bone, and the method can allow elementary surgeon carry out the work such as preoperative planning, surgical simulation rehearsal, operation teaching, guided operation practical operation and postoperative rehabilitation by the information in virtual environment.Contributing to improve the operation precision, rebuild and approach normal acl feature, is a fruitful new way that improves ACL reconstruction curative effect.Its assistance application in ACL rebuilds can be improved preoperative assessment and planning, greatly improves the accuracy of tunnel location, reduces incidence, operating time and surgeon's training time of postoperative complication.
The accompanying drawing explanation
The knee joint that Fig. 1 obtains and anterior cruciate ligament MRI two dimensional image;
Schematic diagram when Fig. 2 sample MRI imports Mimics10.01;
Fig. 3 sample MRI imports the schematic diagram after Mimics10.01;
Fig. 4 moves the schematic diagram before cutting apart three-dimensional reconstruction Mimics10.01 expert;
Fig. 5 sample three-dimensional reconstruction digital model schematic diagram;
Fig. 6 three-dimensional reconstruction sample double knee joint digital model normotopia;
Position after Fig. 7 three-dimensional reconstruction sample double knee joint digital model;
Fig. 8 anterior cruciate ligament angle is lost the shape position;
Fig. 9 anterior cruciate ligament angle Coronal;
Figure 10 anterior cruciate ligament angle horizontal position.
Embodiment
Below in conjunction with accompanying drawing, with implementing example, the present invention is described in further detail, but they are not limitation of the invention.
1. experimental subjects: get the normal double knee joint sample of 20 routine fresh adult corpse (corpse age at death 25～50 years old) as research object, experiment is front-seat except exogenous injury, tumour, regression, deformity etc.
2. software and computer run environment:
Mimics 10.01(Materialise company of medical finite elements simulation software, Belgium), reverse engineering software (Korea S INUS company), finite element analysis software ANSYS13.0(U.S. ANSYS company).
3. allocation of computer:
Intel (R) Core (TM) i5CPU M 430@2.27GHz 2.27GHz, figure video card ATI Mobility Radeon
tMhD 5650 GRAPHICS 1G, hard disk 320G, internal memory 8G.7,64 of operating system Windows.
4. the foundation of three-dimensional digital model
Will be upper by software input port Mimics 10.01(Materialise company of importing medical finite elements simulation software at computer (operating system Windows 7) according to the knee joint nuclear magnetic resonance scanning MRI fault image data of dicom standard storage, Belgium), through the location picture, organize picture, interior difference to process, set suitable threshold value (226～2782Hounsfield unit).Software generates the profile of the structure of every aspect automatically.Carry out at every tomographic image by hand that marginality is cut apart, selective editing, filling-up hole process, and removes artifact and miscellaneous data.Generate shin bone, femur, kneecap, fibula, anterior cruciate ligament, posterior cruciate ligament of knee and meniscus how much facial mask (Mask) files separately by region growing (region growing).Finally set up the complete double knee joint 3D solid digital model (seeing Fig. 1～Fig. 7) that comprises DF, proximal tibia, anterior cruciate ligament, meniscus, kneecap and fibula etc., and derive and preserve with the STL form.
5. the Healing in Anterior Cruciate Ligament Reconstruction model is imported in ANSYS13.0 finite element analysis software working interface with the STL form and shows separately anterior cruciate ligament.
6. the anterior cruciate ligament cross section obtains and centre coordinate numerical value
Obtaining cross section and this kernel of section coordinate in the ANSYS13.0 finite element analysis software is XC=0.23815E-06 YC=-0.21219E-06 ZC=-0.32095E-06.And be parallel to this plane and up or down the anterior cruciate ligament model carried out to cross section respectively, and obtain its center and coordinate values thereof with same method, and finally altogether obtaining 40 central points and coordinate values thereof, concrete numerical value is arranged as follows in order.
Central point and the coordinate values thereof in 40 cross sections that table 1. obtains
7. the matching of anterior cruciate ligament kernel of section curve and curve equation thereof solves
By the coordinate values of the above anterior cruciate ligament kernel of section obtained, after importing ANASYS13.0 software, simulate the center curve of anterior cruciate ligament and solve its equation.
8. the anterior cruciate ligament center curve general equation simulated is:
Z=0.2412+0.7879x-6.5738y+?2.9976x
2+?1.8441y
2+2.5504xy
9. femur and Tibial attachment centre coordinate value tries to achieve
In the ANSYS13.0 finite element analysis software, the spatial value of trying to achieve anterior cruciate ligament axis and femur intersection point by the method for high linear algebra is (x=-4.9497; Y=7.5249; Z=-23.9527).
In like manner, the spatial value that draws anterior cruciate ligament axis and shin bone intersection point is (x=-4.6631; Y=-1.8740; Z=18.7653).
10. the anterior cruciate ligament center curve is in the acquisition of each Plane Angle of coordinate system
By the method for trigonometric function, the concrete value of application femur and shin bone attachment surface Geometric center coordinates numerical evaluation out-of-plane angle, result of calculation is as follows:
Anterior cruciate ligament is 68 in sagittal angle
.; (see Fig. 8, the anterior cruciate ligament angle is lost the shape position)
Anterior cruciate ligament is 22 in the angle of coronal-plane
.; (seeing Fig. 9, anterior cruciate ligament angle Coronal)
Anterior cruciate ligament is 19 in the angle of surface level.(seeing Figure 10, anterior cruciate ligament angle horizontal position)
11. result
11.1. the space curve general equation of anterior cruciate ligament axis is:
Z=0.2412+0.7879x-6.5738y+?2.9976x
2+?1.8441y
2+2.5504xy
11.2. the volume coordinate at the center of anterior cruciate ligament femur and shin bone attachment surface is (x=-4.9497; Y=7.5249; Z=-23.9527); (x=-4.6631; Y=-1.8740; Z=18.7653).
11.3. the angle of anterior cruciate ligament sagittal plane, surface level and coronal-plane is respectively 68
.、22
。with 19
.。
Claims (1)
1. a method that obtains the anterior cruciate ligament axis based on the linear fit algorithm comprises the following steps:
(1) software preparation: the Mimics of medical finite elements simulation software, reverse engineering software and finite element analysis software ANSYS;
(2) foundation of three-dimensional digital model
To import the Mimics of medical finite elements simulation software 10.01 by the software input port on computers according to the knee joint MRI fault image data of dicom standard storage, after locating picture, organizing picture and interior difference to process, automatically generate suitable threshold value; Software generates the profile of the structure of every aspect automatically, carries out at every tomographic image by hand that marginality is cut apart, selective editing, filling-up hole process, and removes artifact and miscellaneous data; Generate shin bone, femur, kneecap, fibula, anterior cruciate ligament, posterior cruciate ligament of knee and meniscus how much facial mask files separately by region growing; Finally set up the complete double knee joint 3D solid digital model that comprises DF, proximal tibia, anterior cruciate ligament, meniscus, kneecap and fibula, with the STL form, derive and preserve;
(3) importing of anterior cruciate ligament three-dimensional digital model:
The anterior cruciate ligament model part of STL form is imported in ANSYS13.0 finite element analysis software working interface, and show separately anterior cruciate ligament;
(4) the anterior cruciate ligament cross section obtains and centre coordinate numerical value
Obtain anterior cruciate ligament cross section and this kernel of section coordinate in the ANSYS13.0 finite element analysis software, and be parallel to the anterior cruciate ligament face and up or down the anterior cruciate ligament model carried out to cross section respectively, and obtain its center and coordinate values thereof with same method;
(5) matching of anterior cruciate ligament kernel of section curve and curve equation thereof solves
By the coordinate values of obtained anterior cruciate ligament kernel of section, simulate the center curve of anterior cruciate ligament and solve its equation after importing ANASYS13.0 software, equation is:
Z=0.2412+0.7879x-6.5738y+?2.9976x
2+?1.8441y
2+2.5504xy
(6) femur and Tibial attachment centre coordinate value tries to achieve
In the ANSYS13.0 finite element analysis software, try to achieve the spatial value of anterior cruciate ligament axis and femur intersection point by the method for high linear algebra; In like manner, draw the spatial value of anterior cruciate ligament axis and shin bone intersection point;
(7) the anterior cruciate ligament central axis again by the method for trigonometric function, is applied the concrete value of femur and shin bone attachment surface Geometric center coordinates numerical evaluation out-of-plane angle in the acquisition of each Plane Angle of coordinate system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2012104870245A CN102968565A (en) | 2012-11-26 | 2012-11-26 | Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2012104870245A CN102968565A (en) | 2012-11-26 | 2012-11-26 | Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102968565A true CN102968565A (en) | 2013-03-13 |
Family
ID=47798702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2012104870245A Pending CN102968565A (en) | 2012-11-26 | 2012-11-26 | Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102968565A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104166762A (en) * | 2014-08-15 | 2014-11-26 | 段少银 | Construction method of finite element model of atlanto-axial joint, ligaments and vertebral artery |
| CN105139442A (en) * | 2015-07-23 | 2015-12-09 | 昆明医科大学第一附属医院 | Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) |
| CN103854288B (en) * | 2014-03-11 | 2017-04-19 | 深圳市旭东数字医学影像技术有限公司 | Cruciate ligament segmentation method and system |
| CN108053470A (en) * | 2017-12-28 | 2018-05-18 | 昆明医科大学第附属医院 | One kind is based on MRI difficult airways three-dimensional rebuilding method and application |
| CN116012437A (en) * | 2023-03-27 | 2023-04-25 | 山东捷瑞数字科技股份有限公司 | Method, device and equipment for extracting three-dimensional model center line based on three-dimensional engine |
-
2012
- 2012-11-26 CN CN2012104870245A patent/CN102968565A/en active Pending
Non-Patent Citations (3)
| Title |
|---|
| 李彦林等: "基于MRI二维影像下股骨髁间窝的三维可视化研究", 《中国修复重建外科杂志》 * |
| 王疆: "逆向工程软件结合Micro_CT技术在牙体组织精细三维有限元模型建立中的应用", 《第四军医大学学报》 * |
| 陈文栋等: "基于MRI建立膝关节前交叉韧带三维数字化模型", 《中国组织工程研究与临床康复》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103854288B (en) * | 2014-03-11 | 2017-04-19 | 深圳市旭东数字医学影像技术有限公司 | Cruciate ligament segmentation method and system |
| CN104166762A (en) * | 2014-08-15 | 2014-11-26 | 段少银 | Construction method of finite element model of atlanto-axial joint, ligaments and vertebral artery |
| CN105139442A (en) * | 2015-07-23 | 2015-12-09 | 昆明医科大学第一附属医院 | Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) |
| CN108053470A (en) * | 2017-12-28 | 2018-05-18 | 昆明医科大学第附属医院 | One kind is based on MRI difficult airways three-dimensional rebuilding method and application |
| CN116012437A (en) * | 2023-03-27 | 2023-04-25 | 山东捷瑞数字科技股份有限公司 | Method, device and equipment for extracting three-dimensional model center line based on three-dimensional engine |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | Application of 3D printing in the surgical planning of trimalleolar fracture and doctor-patient communication | |
| Carbone et al. | TLEM 2.0–A comprehensive musculoskeletal geometry dataset for subject-specific modeling of lower extremity | |
| Chabanas et al. | Patient specific finite element model of the face soft tissues for computer-assisted maxillofacial surgery | |
| Tucker et al. | Comparison of actual surgical outcomes and 3-dimensional surgical simulations | |
| CN108478250A (en) | Femur localization method, individuation osteotomy guide plate and the prosthese of total knee arthroplasty | |
| CN104546060B (en) | Patella osteotomy navigation template for knee arthroplasty and manufacture method | |
| CN105608741A (en) | Computer simulation method for predicting soft tissue appearance change after maxillofacial bone plastic surgery | |
| CN105361883A (en) | Method for determining lower limb biological force line in three-dimensional space for total knee arthroplasty | |
| CN102968565A (en) | Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm | |
| Chen et al. | Research status and application prospects of digital technology in orthopaedics | |
| CN102940530A (en) | Method for virtually building anterior cruciate ligament on femur and tibia tunnels | |
| CN107106239A (en) | Surgery is planned and method | |
| CN106963487A (en) | A kind of discoid meniscus simulates operation method | |
| Sinclair et al. | Three-dimensional kinematic differences between the preferred and non-preferred limbs during maximal instep soccer kicking | |
| CN104112384A (en) | Neurosurgery virtual surgery training system | |
| Liu et al. | Cloud computing-enabled IIOT system for neurosurgical simulation using augmented reality data access | |
| Li et al. | Research and development of the regulation training system of minimally invasive knee surgery based on VR technology | |
| Tang et al. | Development of software for human muscle force estimation | |
| Bakke et al. | Shape-model scaled gait models can neglect segment markers without consequential change to inverse kinematics results | |
| ÖZKADİF et al. | Contribution of virtual anatomic models to medical education | |
| Xie et al. | 3D visualization and simulation in surgical planning system of orbital hypertelorism | |
| Rosen | Advanced surgical technologies for plastic and reconstructive surgery | |
| Chao et al. | Virtual Interactive Musculoskeletal System (VIMS) in orthopaedic translational research | |
| Niu et al. | Reconstruction and optimization of the 3D geometric anatomy structure model for subject-specific human knee joint based on CT and MRI images | |
| Gargiulo et al. | Use of 3 dimensional models and navigation system in preparation for brain surgery |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130313 |