CN105139442A - Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) - Google Patents
Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) Download PDFInfo
- Publication number
- CN105139442A CN105139442A CN201510437619.3A CN201510437619A CN105139442A CN 105139442 A CN105139442 A CN 105139442A CN 201510437619 A CN201510437619 A CN 201510437619A CN 105139442 A CN105139442 A CN 105139442A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- knee joint
- mri
- model
- select
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The invention relates to a method for establishing a human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging), and belongs to the technical field of medical images. The method comprises six steps of CT image acquisition, MRI image acquisition, data storage, CT two-dimensional image importing and three-dimensional model establishment, MRI two-dimensional image importing and three-dimensional model establishment, and knee joint anatomy simulation model registration and modification. The method combines the advantages of CT and MRI scanning, so that the scanning is convenient and fast, the images are clear, the dissection relation is definite, the density resolution is very high, and the scanning range is very wide. According to the method, two sets of images are registered through a three-dimensional image registration technology, three-dimensional coordinate systems of two types of images are unified by fitting point cloud contours and anatomical feature points of femurs, tibias and fibulas, and registration is performed through a registration module to obtain a double-lower-limb full-length knee joint three-dimensional digital model with structures of bones, cartilages, meniscuses, ligaments, muscle tendons and the like.
Description
Technical field
The invention belongs to Medical Imaging Technology field, be specifically related to a kind of method setting up complete human body's knee joint geometry dissection realistic model in conjunction with CT and MRI two dimensional image.
Background technology
CT and MRI is Imaging Technology conventional clinically at present, and the two has respective relative merits.CT full name computed tomography (computedtomography, CT), launches by scanner the transverse scan that x-ray carries out successively to a certain scope of human body, obtains information, obtain the cross-section internal anatomy of human body as calculated after machine process.CT scan is convenient, rapidly, clear image, anatomy relationship is clear and definite, has very high density resolution, can distinguish the different tissues that density only differs 5-6Hu.Especially bone tissue is observed and there is obvious advantage, the density of different bone tissue and bone tissue and soft tissue interface can be separated in circle of good definition, but not good enough to knee joint soft tissue resolving effect.MRI full name Magnetic resonance imaging (NuclearMagneticResonanceImaging, MRI), utilize nuclear magnetic resonance (nuclearmagneticresonnance, NMR) principle, according to the decay that the energy discharged is different in the inner different structure environment of material, detected the electromagnetic wave launched by additional gradient magnetic, can learn and form the nuclear position of this object and kind, the structural images of interior of articles can be depicted as accordingly.MRI adopts multiparameter imaging, high to soft tissue resolution, object to different parts different structure soft-tissue imaging can be reached by regulating imaging parameters, also be free to by regulating magnetic field select desired profile, can obtain other imaging technique institute can not the close or image that is difficult to close to position.But MRI is limited by the size of scanning yoke, and single scan range is limited.
Present inventor studies rear discovery, sets up in the process of three-dimensional model at two dimensional image, and the model obtained based on different scan modes can be variant to the precision of different tissues.The change of cartilage stress cannot be analyzed with CT to the complete finite element model that knee joint surface direct digital measurement is set up, especially can not carry out modeling satisfactorily to soft tissues such as meniscuss.Although the three-dimensional finite element model built based on MRI two dimensional image is building the deficiency that compensate for CT in meniscus, ligament, cartilage etc., but the gray scale difference value of organizational boundary is less, a large amount of artificial corrections is needed in the process set up, considerably increase workload, there is certain error, to the structure less effective of cortex of bone, articular cavity structure cannot be set up exactly, and the three-dimensional model constructing lower limb total length that cannot be complete, there is limitation for the structure lower limb line of force and dissection axis.
Therefore, we are on the basis of large quantity research, to same volunteer simultaneously row CT and MRI scan.Accurately construct the bone tissue three-dimensional digitalization model of lower limb total length based on CT data, accurately construct the three-dimensional digitalization model of knee-joint anatomy structure based on MRI data.In computer software, by the technology of 3-D view registration, two groups of images are carried out registration, by the position of matching femur, shin bone, fibula, just obtains and comprise bone, cartilage, meniscus, ligament and tendon isostructural pair of lower limb total length knee joint three-dimensional digitalization model.This model, except accurately reflecting kneed anatomical structure, accurately can also draw the lower limb line of force, ossis center line, articular cavity axis, rotation etc.Find a kind of method setting up complete human body's knee joint geometry dissection realistic model, for the biomechanics Research of knee joint and joint replacement provides good basic model, this in the past researcher build knee joint model by independent a kind of scan method and cannot accomplish.
Summary of the invention
The object of the invention is to for the deficiencies in the prior art, overcome single utilization CT or MRI two dimensional image set up knee joint three-dimensional model cannot the shortcoming of complete model knee-joint anatomy structure, in conjunction with both respective advantage, by the position of the bone tissues such as matching femur, shin bone, fibula, in conjunction with soft tissue anatomical data, obtain and comprise bone, cartilage, meniscus, ligament and tendon isostructural pair of lower limb total length knee joint three-dimensional digital and dissect realistic model.
The technical solution used in the present invention is as follows:
Set up a method for human body knee joint three-dimensional simulation model in conjunction with CT and MRI two dimensional image, comprise the steps:
Step (1), CT image acquisition:
Adopt the two lower limb of CT scan, scanning position is that knee joint stretches naturally and outward turning 10 ° ~ 15 ° is fixed; Sweep limit be top to midpelvic plane, below comprises foot completely;
Step (2), MRI image acquisition:
Adopt the two knee of magnetic resonance machine scanning, magnetic resonance machine is using head coil as receiving coil; Scanning position be knee joint naturally stretch and outward turning 10 ° ~ 15 ° degree angle fix; Sweep limit is centered by knee space, respectively scans 10cm up and down;
Step (3), data are preserved:
The scan image data of step (1) and step (2) gained is stored with DICOM3.0 (DigitalImagingandCommunicationinMedicine) form and is burnt on CD-ROM on computer workstation;
Step (4), the importing of CT two dimensional image and three-dimension modeling:
On computer workstation, knee joint CT scanning view data is imported Mimicsl6.0 in dicom format, behind definition up, down, left, right, before and after direction, in Mimicsl6.0, demonstrate the two dimensional image of sagittal plain, Coronal, volume shape position;
A kind of green (Green) is newly set up in " Masks ", select " Thresholding " in " SegmentationMenu ", covered the bone tissue of structure at all levels completely to green " Masks " by lower threshold value in " Thresholding " adjustment, the CT value of bone tissue is between 226-1542Hu;
Select " Erase " in " EditMasks " menu to wipe green unnecessary around femur, define bone tissue edge; Click GreenMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different bone;
Select " Calculate3DfromMasks " in " SegmentationMenu ", adopt " HighQuality " computing method, after operation, set up out knee joint three-dimensional digitalization model;
Only there is bone tissue model with the knee joint three-dimensional digitalization model that CT two dimensional image is set up, comprise femur, shin bone, fibula, kneecap.
Step (5), the importing of MRI two dimensional image and three-dimension modeling:
The first step, on computer workstation, imports Mimicsl6.0 in dicom format by knee joint MRI scan image data, behind definition up, down, left, right, before and after direction, demonstrates the two dimensional image of sagittal plain, Coronal, volume shape position in Mimicsl6.0; A kind of green (Green) is newly set up in " Masks ", select " Thresholding " in " SegmentationMenu ", first covered the bone tissue of structure at all levels by lower threshold value in " Thresholding " adjustment to green " Masks " completely; Select " Erase " in " EditMasks " menu to wipe green unnecessary around femur, define bone tissue edge; Click GreenMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different bone;
Second step sets up a kind of blueness (Blue) newly in " Masks ", and the upper lower threshold value of adjustment " Thresholding " covers the ligament of structure at all levels completely to blue " Masks "; Select " Erase " in " EditMasks " menu to wipe blueness unnecessary around ligament, define all ligament tissue edges; Click BlueMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different ligament;
3rd step sets up a kind of redness (Red) newly in " Masks ", and the upper lower threshold value of adjustment " Thresholding " covers cartilage and the meniscus of structure at all levels completely to red " Masks "; Select " Erase " in " EditMasks " menu to wipe the unnecessary redness of periarticular, define all cartilages and meniscal tissue edge; In RedMasks, select " RegionGrowing ", click different cartilages or meniscal tissue, " NewMasks ", complete different cartilage and meniscal division; Click RedMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete different cartilage and meniscal division;
Select " Calculate3DfromMasks " option in " SegmentationMenu " menu bar, adopt in dialog box and click " HighQuality " computing method, after operation, set up out knee joint three-dimensional digitalization model;
Have bone tissue model and soft tissue model two parts with the knee joint three-dimensional digitalization model that MRI two dimensional image is set up, wherein bone tissue model comprises femur, shin bone, fibula, kneecap; Soft tissue model comprises medial meniscus, lateral meniscus, anterior cruciate ligament, posterior cruciate ligament of knee, MCL, lateral collateral ligament, femoral cartilage, tibial cartilage, fibula cartilage, Patella Cartilage, quadriceps muscle of thigh, kneecap tendon.
MRI image is higher to soft tissue resolution, meniscus, cartilage and ligament model can be set up out, but its gray scale difference is less than CT data, when adopting computing machine automatically to calculate, noise is on the high side, need the method by manually deleting benefit, fill the cavity in tissue regions, remove the burr of segmenting edge and the noise of other positions introducing.Especially need according to anatomical structure, carry out careful defining to organization edge, defining of edge causes very large impact by the really degree of institute's Modling model.
Step (6), knee-joint anatomy realistic model registration and modification:
By the three-dimensional digitalization model obtained by CT and MRI scan image data, namely the knee joint three-dimensional digitalization model Imageware13.0 software that step (4) the knee joint three-dimensional digitalization model that adopts CT scan view data to obtain and step (5) adopt MRI scan image data to obtain is converted to point cloud model, by in registration module in Imageware13.0 software of the some cloud profile of the bone tissue obtained by MRI and the bone tissue point cloud profile obtained by CT by optimum matching computing, the two is made to put cloud profile registration, and finely tune based on clear and definite anatomical features point, the space coordinates of the three-dimensional digitalization model that the three-dimensional digitalization model that CT scan view data is obtained and MRI scan image data obtain are unified,
Then again the some cloud profile of the soft tissue obtained by MRI and bone tissue is imported in these space coordinates, again by registration module by soft tissue and bone tissue registration, and finely tune based on clear and definite anatomical features point, namely obtain human body knee joint geometry and dissect realistic model.
In technique scheme, the CT described in step (1) is that SIEMENS128 arranges double source spiral CT.
Further, preferably, the sweep parameter of CT scan is set as: tube voltage 120kV and 70kV of double source CT, tube current 100mA and 60mA, thickness 0.6mm, interlayer every 0.6mm, pitch 0.5.
In technique scheme, the magnetic resonance machine described in step (2) is GE1.5T superconduct magnetic resonance machine.
Further, preferably, the sweep parameter of magnetic resonance machine scanning is set as: sagittal plain 3D proton density weighted imaging sequence, TR11000ms, TE25ms; Thickness 1.0mm; Interlamellar spacing 0.2mm; Echo train 14; Encourage 2 times; Matrix 192/320; FOV18.
Compared with prior art, its beneficial effect is in the present invention:
1. the inventive method combines the advantage that both CT and MRI scan, and scanning is convenient, rapidly, clear image, anatomy relationship is clear and definite, has very high density resolution and sweep limit.
2. by 3-D view registration technology, two groups of images are carried out registration, by matching femur, shin bone, the point cloud profile of fibula and dissection unique point, the three-dimensional system of coordinate of unified two kinds of images, on this basis by soft tissue model by registration module and bone tissue Model registration, [soft tissue model refer to set up in the middle of step (5) comprise meniscus, anterior cruciate ligament, posterior cruciate ligament of knee, MCL, lateral collateral ligament, articular cartilage, quadriceps muscle of thigh, kneecap tendon is at interior three-dimensional model, what bone tissue model referred to that CT and MRI set up comprises femur, shin bone, fibula, kneecap is at interior three-dimensional model], just obtain and comprise bone, cartilage, meniscus, ligament and tendon isostructural pair of lower limb total length knee joint three-dimensional digitalization model.
3. the model set up of this method is except accurately reflecting kneed anatomical structure, accurately can also draw the lower limb line of force, ossis center line, articular cavity axis, rotation etc., for the biomechanics Research of knee joint and joint replacement provides good basic model.
Accompanying drawing explanation
Fig. 1 is knee joint CT scanning two dimensional image, and wherein a represents sagittal plain; B represents Coronal; C represents cross-section position, can find out it is two lower limb total length positions.
Fig. 2 is that knee joint MRI scans two dimensional image, and wherein a represents sagittal plain; B represents Coronal; C represents cross-section position.
CT two dimensional image imports in Mimicsl6.0 software by Fig. 3.
Fig. 4 is the two lower limb three-dimensional digitalization models set up based on CT two dimensional image, and wherein a represents Coronal, and b represents cross-section position, and c represents sagittal plain, and d represents three-dimensional digitalization model.
Fig. 5 imports based on MRI two dimensional image in Mimicsl6.0 software.
Fig. 6 is the knee joint three-dimensional digitalization model set up based on MRI two dimensional image, and wherein a represents Coronal, and b represents cross-section position, and c represents sagittal plain, and d represents three-dimensional digitalization model.
Fig. 7 dissects realistic model to carrying out the human body knee joint geometry after registration based on CT and MRI two kinds of data in Imgeware13.0 software.
Fig. 8 is the full knee system of PFCSigma stationary platform, and wherein a represents femoral prosthesis; B represents tibial prosthesis; C represents the rear prosthese of assembling.
Fig. 9 carries out 3-D scanning to knee-joint prosthesis, and wherein a represents that being coated with bletilla to prosthetic surface locates; B represents the scanning process of three-dimensional laser scanner.
Figure 10 is the three-dimensional digitalization model being set up knee-joint prosthesis by GeomagicStudio12 software.
Figure 11 simulates femoral head geometric center by femoral head point cloud model.
Figure 12 draws the lower limb line of force by Hypermesh11.0 software.
Figure 13 is normal person's lower limb line of force and osteotomy schematic diagram, and wherein a represents normal person's lower limb line of force, and b represents by lower limb line of force osteotomy schematic diagram.
Figure 14 simulates thighbone end prosthesis cross section by Hypermesh11.0 software.
Figure 15 is femoral-posterior condyles osteotomy.
Figure 16 is femur osteotomy surface.
Figure 17 is Tibial osteotomy face.
Figure 18 is three-dimensional digitalization model after knee joint virtual osteotomy, and wherein a represents that side is seen, and b represents that front is seen.
Figure 19 is three-dimensional digitalization model after TKA, and wherein a represents that side is seen, and b represents that front is seen.
Figure 20 is three-dimensional finite element model after computer virtual total knee replacement.
Figure 21 be Three-Dimensional Dynamic finite element model difference go down on one's knees angle time femur rotating relative to shin bone.
Figure 22 be Three-Dimensional Dynamic finite element model difference go down on one's knees angle time femur relative shin bone displacement.
Embodiment
Below in conjunction with embodiment, the present invention is described in further detail.
It will be understood to those of skill in the art that the following example only for illustration of the present invention, and should not be considered as limiting scope of the present invention.Unreceipted concrete technology or condition person in embodiment, according to the technology described by the document in this area or condition or carry out according to product description.Agents useful for same or the unreceipted production firm person of instrument, being can by buying the conventional products obtained.
Experimental example 1
Human body knee joint 3 D anatomical realistic model is set up based on CT and MRI two dimensional image
Step (1), CT image acquisition
SIEMENS128 is adopted to arrange the two lower limb of double source Spiral CT scan.Scanning position: knee joint stretches naturally and outward turning 10 ° ~ 15 ° angles are fixed.Sweep limit: top is to midpelvic plane, and below comprises foot completely.CT image be mainly used in observe bone tissue, sweep parameter is set as: tube voltage 120kV and 70kV of double source CT, tube current 100mA and 60mA, thickness 0.6mm, interlayer every 0.6mm, pitch 0.5.
Step (2), MRI image acquisition
Adopt the two knee of GE1.5T superconduct magnetic resonance machine scanning, magnetic resonance machine is using head coil as receiving coil.Scanning position be knee joint naturally stretch and outward turning 10 ° ~ 15 ° degree angle fix.Sweep limit: centered by knee space, respectively scans 10cm up and down.MRI image is mainly used in observing cartilage and ligament tissue, and retouching setting parameter is: sagittal plain 3D proton density weighted imaging sequence, TR11000ms, TE25ms; Thickness 1.0mm; Interlamellar spacing 0.2mm; Echo train 14; Encourage 2 times; Matrix 192/320; FOV18.
Step (3), data are preserved
Step (1) and step (2) gained scan image data are stored with DICOM3.0 (DigitalImagingandCommunicationinMedicine) form on a workstation and be burnt on CD-ROM, as shown in Figure 1, MRI scan image as shown in Figure 2 for CT scan image.
Dicom standard belongs to medical information system field, and it is mainly used in exchanging numerical information between medical image equipment.The DICOM compatible equipment that different vendor produces can carry out interconnected easily, can be used for the Image Communication between two Medical Devices and the interface between image acquisition equipment and image processing workstations.
Step (4), the importing of CT two dimensional image and three-dimension modeling
On computer workstation, knee joint CT scanning image is imported Mimicsl6.0 in dicom format, behind definition up, down, left, right, before and after direction, demonstrate the two dimensional image of sagittal plain, Coronal, volume shape position in Mimicsl6.0, as shown in Figure 3;
A kind of green (Green) is newly set up in " Masks ", select " Thresholding " in " SegmentationMenu ", covered the bone tissue (CT value is between 226-1542Hu) of structure at all levels by lower threshold value in " Thresholding " adjustment to green " Masks " completely.Select " Erase " in " EditMasks " menu to wipe green unnecessary around femur, define bone tissue edge.Click GreenMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different bone;
Select " Calculate3DfromMasks " in " SegmentationMenu ", adopt " HighQuality " computing method, set up out knee joint three-dimensional digitalization model after operation, as shown in Figure 4.
Step (5), the importing of MRI two dimensional image and three-dimension modeling
The first step, on computer workstation, knee joint MRI scan image data is imported Mimicsl6.0 in dicom format, behind definition up, down, left, right, before and after direction, the two dimensional image of sagittal plain, Coronal, volume shape position is demonstrated, as shown in Figure 5 in Mimicsl6.0; In " Masks ", set up a kind of green newly, select " Thresholding " in " SegmentationMenu ", first covered the bone tissue of structure at all levels by lower threshold value in " Thresholding " adjustment to green " Masks " completely; Select " Erase " in " EditMasks " menu to wipe green unnecessary around femur, define bone tissue edge; Click GreenMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different bone;
Second step sets up a kind of blueness newly in " Masks ", and the upper lower threshold value of adjustment " Thresholding " covers the ligament of structure at all levels completely to blue " Masks "; Select " Erase " in " EditMasks " menu to wipe blueness unnecessary around ligament, define all ligament tissue edges; Click BlueMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different ligament;
3rd step sets up a kind of redness newly in " Masks ", and the upper lower threshold value of adjustment " Thresholding " covers cartilage and the meniscus of structure at all levels completely to red " Masks "; Select " Erase " in " EditMasks " menu to wipe the unnecessary redness of periarticular, define all cartilages and meniscal tissue edge; Click RedMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete different cartilage and meniscal division;
Select " Calculate3DfromMasks " option in " SegmentationMenu " menu bar, adopt in dialog box and click " HighQuality " computing method, set up out knee joint three-dimensional digitalization model after operation, as shown in Figure 6;
Set up the process of three-dimensional model based on MRI two dimensional image and CT image method basically identical, but it should be noted, MRI image is higher to soft tissue resolution, meniscus, cartilage and ligament model can be set up out, but its gray scale difference is less than CT data, when adopting computing machine automatically to calculate, noise is on the high side, needs the method by manually deleting benefit, fill the cavity in tissue regions, remove the burr of segmenting edge and the noise of other positions introducing.Especially need according to anatomical structure, carry out careful defining to organization edge, defining of edge causes very large impact by the really degree of institute's Modling model.
Step (6), knee joint model registration and modification
By the three-dimensional digitalization model obtained by CT and MRI scan image data, namely the knee joint three-dimensional digitalization model Imageware13.0 software that step (4) the knee joint three-dimensional digitalization model that adopts CT scan view data to obtain and step (5) adopt MRI scan image data to obtain is converted to point cloud model, so just can obtain entirety point cloud profile and the clear and definite anatomical features point of bone tissue, therefore we adopt the principle of 3-D view registration, by in registration module in Imageware13.0 software of the some cloud profile of the bone tissue obtained by MRI and the bone tissue point cloud profile obtained by CT by optimum matching computing, the two is made to put cloud profile registration, and finely tune based on clear and definite anatomical features point, so just the space coordinates of the three-dimensional digitalization model obtained by different pieces of information are unified.
Then again the some cloud profile of the soft tissue obtained by MRI and bone tissue is imported in these space coordinates, again by registration module by soft tissue and bone tissue registration, and finely tune based on clear and definite anatomical features point, namely obtain comparatively accurate and complete human body knee joint geometry and dissect realistic model, as shown in Figure 7.
Experimental example 2
The foundation of three-dimensional finite element model after TKA
1. prosthese scans and sets up three-dimensional digitalization model
Carry out white after being connected with register mark point by complete for PFCSigma stationary platform knee system prosthese to spray, scan with three-dimensional laser scanner, as shown in Figure 8 and Figure 9.Scan-data is imported GeomagicStudio12 reverse engineering software, by anchor point registration three-dimensional point cloud model, and carry out curved surface optimization.The cavity of model surface is filled by the principle first based on continual curvature, removes the non-characteristic impression of model surface, and smooth release surface, prevents the generation in non-distinctive higher curvature and self intersection face.Then the change of descriptive model surface curvature is carried out with outline line, adopt the U-V mesh lines of the every block occluding contour of parametrization regulative mode condition, thus regulate the fineness of NURBS (non-uniform rational B-spline) curved surface, realize coming model of fit surface with smooth curved surface, finally in software each 3D solid module set up is carried out fitting assembling, modules is assembled into many bodies, forms the three-dimensional digitalization model of a total knee prostheses, as shown in Figure 10.
2. computer virtual osteotomy
Knee joint three-dimensional model is imported Hypermesh11.0 software with STL form, the geometric center of femoral head point cloud model is simulated by software, as shown in figure 11, a line is done to ankle-joint center from late-segmental collapse, on coronal-plane, this line, inside fossa intercondylaris femoris center deviation, is femur axial alignment, as shown in figure 12.This line becomes 3 ° to turn up with centre of body weight line, is 6 ° with femur anatomical axis wire clamp angle.The line at both centers of tibial plateau and astragalus forms shank line of force axis, also has the valgus angle of 3 ° with centre of body weight line.On coronal-plane, femur mechanics axle overlaps with shin bone axial alignment, as shown in figure 13.Prosthese three-dimensional model imports Hypermesh11.0 software with STL form, simulates knee-joint prosthesis femur end section, as femur osteotomy instrument, as shown in figure 14 by software.
Conveniently replacement knee in arthroplasty requirement, femur, Tibial osteotomy face are perpendicular to the lower limb line of force.Postoperative at knee prosthesis, when two lower limb are stood naturally, knee joint transverse axis is parallel to ankle-joint and ground, has so just recovered the normal forces credit cloth of articular surface.Femoral prosthesis Turning matched is with femur surgery epicondyle line for reference axis, and it is closest to stretch of knee joint axle.DF cross section orthogonal is in the femur line of force, femoral-posterior condyles cross section is parallel with surgery epicondyle line (depression of internal epicondyle of femur vicinity and the line of external epicondyle of femur), and perpendicular to front and back axis (i.e. Whiteside line, the line on femoral bone pulley minimum point and fossa intercondyloidea summit), place relative to the 3 ° of outward turnings of femoral-posterior condyles axle, as shown in figure 15.DF osteotomy and postartis osteotomy site carry out osteotomy with knee-joint prosthesis cross section to femur after determining, obtain femur osteotomy surface, as shown in figure 16.Proximal tibia osteotomy direction is that coronal-plane and sagittal plane are all hung down as Tibia Anatomy axle, without hypsokinesis, as shown in figure 17.Respectively in " solid segmentation " module of computer software, simulation osteotomy is carried out to femur and shin bone according to above-mentioned osteotomy surface, obtain three-dimensional digitalization model after knee joint virtual osteotomy, as shown in figure 18.Measure osteotomy surface size, select No. 3 femoral prosthesis and No. 3 tibial prosthesis, polyethylene liner mat thickness 8mm.
3. joint prosthesis knee-joint prosthesis Model Mounting
Respectively shifting function is carried out to each prosthesis assembly model in software, make femoral prosthesis model be placed in knee joint femoral far-end.The rotation carrying out tibial prosthesis with Tibial osteotomy face transverse axis and antero posterior axis intersection point for rotating neutral position is placed.High-molecular polythene shim models and tibia support model fit like a glove, and measure its concave surface and femoral prosthesis perigee distance 1mm, three-dimensional digitalization model after acquisition TKA, as shown in figure 19.
4. the postoperative three-dimensional model gridding of knee prosthesis divides
In Hypermesh11.0 software to knee replacements after model carry out stress and strain model, use tetrahedral grid partitioning model, like this can reserving model geological information preferably.Three-dimensional finite element model complete after comprising femur, shin bone, fibula, interior lateral collateral ligament, femoral prosthesis, tibial prosthesis and tygon pad total knee replacement is finally built into, as shown in figure 20 after stress and strain model.
Experimental example 3
Dynamic Finite Element analysis after TKA
1. the coordinate of knee joint Three-Dimensional Dynamic finite element model calculates
Stretch relative to the flexing of shin bone according to knee sprung motion analysis femur, interior outward turning, adduction abduction angle description, adopt the define method of card red angular coordinate conversion to set up knee joint Three-Dimensional Dynamic model.The coordinate of definition knee joint tibial when to stretch is reference frame 1, knee joint femoral coordinate system is defined as fixed coordinate system, other phase place knee joint femoral coordinate systems are defined as moving coordinate system 2, then knee joint femoral stretches from 0 ° the Z-Y-X Eulerian angle that a coordinate system 1 transforms to other phase place coordinate systems 2 and is respectively φ, θ, ψ, wherein angle φ represents the interior outward turning of femur relative to shin bone, θ represents the adduction abduction of the relative shin bone of femur, and ψ represents flexing or the stretching, extension of the relative shin bone of femur.Computing formula is as follows:
ψ＝Atan2(r
32，r
33)
For the front and back of knee joint femoral relative to shin bone, the description of inside and outside and upper and lower translation, similarly, it is reference frame 3 that a shin bone coordinate is stretched in definition, 0 ° of femur coordinate is fixed coordinate system 1, other phase place femur coordinates are moving coordinate system 2, then knee joint femoral stretches a coordinate system 1 translation transformation from 0 ° and is respectively along shin bone coordinate system 3 to the front and back of other phase place coordinate systems 2, inside and outside and upper and lower translation.Represent translation inside and outside the relative shin bone of femur with X, Y represents the anterior-posterior translation of the relative shin bone of femur, and Z represents the upper and lower translation of the relative shin bone of femur.Computing formula is as follows:
3T(X,Y,Z)＝
3O
1(X
1,Y
1,Z
1)-
3O
2(X
2,Y
2,Z
2)
2. Definition Model material properties
Each structured material attribute of model is defined as follows: (table 1).
The each structured material attribute of three-dimensional model after table 1 TKA
3. the load of knee joint Three-Dimensional Dynamic model and boundary condition
Three-dimensional finite element model after total knee arthroplasty is applied to the quadriceps muscle of thigh power of 400N, the action direction of quadriceps muscle of thigh, along being parallel to femoral shaft direction, points to quadriceps muscle of thigh starting point.And apply the power of 300N with simulated body weight along the knee joint line of force.Femur is under applied muscular force controls, and relative shin bone carries out the motion of full degree of freedom.Motion is not initiatively applied to shin bone in a model.The motion of shin bone is determined by other the loading environment be distributed on bone, articular cartilage, muscle and ligament.For femur, high-molecular polythene pad, tibial plateau, patellar prosthesis and other soft tissue, define 8 surface area contact.The friction factor of high-molecular polythene and cobalt chrome molybdenum material is defined as 0.04.Contact definition have employed band weighting factor penalty function method, therefore, contact force be defined as interarea invade auxiliary identity distance from function.Select weighting factor so that interarea (femoral prosthesis surface) is larger for the impact of contact intrusion calculating.Penalty function contact allows definition contact between rigid body.
Experimental example 4
Knee prosthesis postoperative Three-Dimensional Dynamic finite element model motion conditions is analyzed
In Abaqus6.10 software, load load and the boundary condition of knee joint Three-Dimensional Dynamic model, obtained the data of knee joint relative motion by computing machine.The relative shin bone of femur obtains according to the Z-Y-X Eulerian angle method around moving coordinate system, the wherein interior outward turning of femur, outreach adduction and range of flexion represent femur coordinate system under the different range of flexion of knee joint respectively successively around the angle of stretching reference frame Z, Y and X-axis and rotating, defining around reference frame Z, Y and X-axis is negative just clockwise, is just counterclockwise.In table, femoral rotation data definition inward turning is negative, and outward turning is just, abduction is negative, and adduction is just, as shown in figure 21.Translation data represent postartis central point X, Y along shin bone reference frame and the translation of Z axis inside and outside knee joint femoral under different flexion angle respectively, define along each axle of coordinate system forward on the occasion of, as shown in figure 22.
Artificial joint replacement patella in the process of 0 ° of-90 ° of flexing along with the intensification of knee sprung, the outward turning of the relative shin bone of femur constantly increases, maximum outward turning about 13 degree, the relative outward turning of flexing femur 90 °-120 ° time starts to reduce between 30 ° of-90 ° of flexings, there is the abduction of about 2.5 ° in femur, after more than 90 °, femur starts adduction relative to shin bone, to adduction about 8 degree during flexing about 130 degree.Axially along with flexing is deepened, during displacement patella 0 °-90 °, femur translation backward increases, and when 90 °-120 °, translation backward reduces.Distally there is the minute movement of about 3mm in condylus lateralis femoris, departs from tibial plateau surfaces be lifted away from when high flexing to flexing about 60 ° from stretching position.In whole flexing process, interior lateral direction occurs in the less translation within 5mm.
The embodiment of the present invention 1 is the basic model of embodiment 2-4, and embodiment 2-4 is the further research done on basis of the present invention.The present invention establishes the three-dimensional model of Normal Knee, on the basis of this model, is established the model of total knee replacement by embodiment 2.In the research in future, can also the present invention be passed through, set up the three-dimensional model etc. of single condyle knee replacements, can follow-up experimental analysis be carried out.
More than show and describe ultimate principle of the present invention, principal character and advantage of the present invention.The technician of the industry should understand; the present invention is not restricted to the described embodiments; what describe in above-described embodiment and instructions just illustrates principle of the present invention; without departing from the spirit and scope of the present invention; the present invention also has various changes and modifications, and these changes and improvements all fall in the claimed scope of the invention.Application claims protection domain is defined by appending claims and equivalent thereof.
Claims (5)
1. set up a method for human body knee joint three-dimensional simulation model in conjunction with CT and MRI two dimensional image, it is characterized in that, comprise the steps:
Step (1), CT image acquisition:
Adopt the two lower limb of CT scan, scanning position is that knee joint stretches naturally and outward turning 10 ° ~ 15 ° is fixed; Sweep limit be top to midpelvic plane, below comprises foot completely;
Step (2), MRI image acquisition:
Adopt the two knee of magnetic resonance machine scanning, magnetic resonance machine is using head coil as receiving coil; Scanning position be knee joint naturally stretch and outward turning 10 ° ~ 15 ° degree angle fix; Sweep limit is centered by knee space, respectively scans 10cm up and down;
Step (3), data are preserved:
The scan image data of step (1) and step (2) gained is stored with DICOM3.0 form and is burnt on CD-ROM on computer workstation;
Step (4), the importing of CT two dimensional image and three-dimension modeling:
On computer workstation, knee joint CT scanning view data is imported Mimicsl6.0 in dicom format, behind definition up, down, left, right, before and after direction, in Mimicsl6.0, demonstrate the two dimensional image of sagittal plain, Coronal, volume shape position;
A kind of green is newly set up in " Masks ", select " Thresholding " in " SegmentationMenu ", covered the bone tissue of structure at all levels completely to green " Masks " by lower threshold value in " Thresholding " adjustment, the CT value of bone tissue is between 226-1542Hu;
Select " Erase " in " EditMasks " menu to wipe green unnecessary around femur, define bone tissue edge; Click GreenMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different bone;
Select " Calculate3DfromMasks " in " SegmentationMenu ", adopt " HighQuality " computing method, after operation, set up out knee joint three-dimensional digitalization model;
Step (5), the importing of MRI two dimensional image and three-dimension modeling:
The first step, on computer workstation, imports Mimicsl6.0 in dicom format by knee joint MRI scan image data, behind definition up, down, left, right, before and after direction, demonstrates the two dimensional image of sagittal plain, Coronal, volume shape position in Mimicsl6.0; In " Masks ", set up a kind of green newly, select " Thresholding " in " SegmentationMenu ", first covered the bone tissue of structure at all levels by lower threshold value in " Thresholding " adjustment to green " Masks " completely; Select " Erase " in " EditMasks " menu to wipe green unnecessary around femur, define bone tissue edge; Click GreenMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different bone;
Second step sets up a kind of blueness newly in " Masks ", and the upper lower threshold value of adjustment " Thresholding " covers the ligament of structure at all levels completely to blue " Masks "; Select " Erase " in " EditMasks " menu to wipe blueness unnecessary around ligament, define all ligament tissue edges; Click BlueMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete the division of different ligament;
3rd step sets up a kind of redness newly in " Masks ", and the upper lower threshold value of adjustment " Thresholding " covers cartilage and the meniscus of structure at all levels completely to red " Masks "; Select " Erase " in " EditMasks " menu to wipe the unnecessary redness of periarticular, define all cartilages and meniscal tissue edge; Click RedMasks by right key, select " RegionGrowing " in a menu, click the different bone tissues in two dimensional image, in dialog box, select " NewMasks ", complete different cartilage and meniscal division;
Select " Calculate3DfromMasks " option in " SegmentationMenu " menu bar, adopt in dialog box and click " HighQuality " computing method, after operation, set up out knee joint three-dimensional digitalization model;
Step (6), knee-joint anatomy realistic model registration and modification:
By the three-dimensional digitalization model obtained by CT and MRI scan image data, namely the knee joint three-dimensional digitalization model Imageware13.0 software that step (4) the knee joint three-dimensional digitalization model that adopts CT scan view data to obtain and step (5) adopt MRI scan image data to obtain is converted to point cloud model, by in registration module in Imageware13.0 software of the some cloud profile of the bone tissue obtained by MRI and the bone tissue point cloud profile obtained by CT by optimum matching computing, the two is made to put cloud profile registration, and finely tune based on clear and definite anatomical features point, the space coordinates of the three-dimensional digitalization model that the three-dimensional digitalization model that CT scan view data is obtained and MRI scan image data obtain are unified,
Then again the some cloud profile of the soft tissue obtained by MRI and bone tissue is imported in these space coordinates, again by registration module by soft tissue and bone tissue registration, and finely tune based on clear and definite anatomical features point, namely obtain human body knee joint geometry and dissect realistic model.
2. the method setting up human body knee joint three-dimensional simulation model in conjunction with CT and MRI two dimensional image according to claim 1, is characterized in that, the CT described in step (1) is that SIEMENS128 arranges double source spiral CT.
3. the method setting up human body knee joint three-dimensional simulation model in conjunction with CT and MRI two dimensional image according to claim 2, it is characterized in that, the sweep parameter of CT scan is set as: tube voltage 120kV and 70kV of double source CT, tube current 100mA and 60mA, thickness 0.6mm, interlayer every 0.6mm, pitch 0.5.
4. the method setting up human body knee joint three-dimensional simulation model in conjunction with CT and MRI two dimensional image according to claim 1, is characterized in that, the magnetic resonance machine described in step (2) is GE1.5T superconduct magnetic resonance machine.
5. the method setting up human body knee joint three-dimensional simulation model in conjunction with CT and MRI two dimensional image according to claim 4, it is characterized in that, the sweep parameter of magnetic resonance machine scanning is set as: sagittal plain 3D proton density weighted imaging sequence, TR11000ms, TE25ms; Thickness 1.0mm; Interlamellar spacing 0.2mm; Echo train 14; Encourage 2 times; Matrix 192/320; FOV18.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510437619.3A CN105139442A (en) | 2015-07-23 | 2015-07-23 | Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510437619.3A CN105139442A (en) | 2015-07-23 | 2015-07-23 | Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105139442A true CN105139442A (en) | 2015-12-09 |
Family
ID=54724774
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510437619.3A Pending CN105139442A (en) | 2015-07-23 | 2015-07-23 | Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105139442A (en) |
Cited By (48)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930617A (en) * | 2016-05-17 | 2016-09-07 | 南方医科大学 | Method for designing and forming stiffness-controllable bone tumor defect repair implant |
| CN106055848A (en) * | 2016-07-14 | 2016-10-26 | 哈尔滨理工大学 | Microstructure parameter-based individualized articular cartilage simulation method |
| CN106175802A (en) * | 2016-08-29 | 2016-12-07 | 吉林大学 | A kind of in body osteoarthrosis stress distribution detection method |
| CN106264731A (en) * | 2016-10-11 | 2017-01-04 | 昆明医科大学第附属医院 | A kind of method based on point-to-point registration technique virtual knee joint single condyle replacement model construction |
| CN106963487A (en) * | 2017-04-17 | 2017-07-21 | 四川大学华西医院 | Knee joint disc-shaped meniscus simulation operation method |
| CN107016666A (en) * | 2017-02-28 | 2017-08-04 | 四川大学华西医院 | Image post-processing method for pelvis CT three-dimensional reconstruction based on coordinate system |
| CN107913100A (en) * | 2017-12-21 | 2018-04-17 | 成都真实维度科技有限公司 | A kind of definite method of late-segmental collapse point |
| CN108042217A (en) * | 2017-12-21 | 2018-05-18 | 成都真实维度科技有限公司 | A kind of definite method of three dimensions lower-limbs biology force-line |
| CN108143487A (en) * | 2017-12-21 | 2018-06-12 | 成都真实维度科技有限公司 | A kind of determining method of ankle-joint central point |
| CN108175378A (en) * | 2017-12-21 | 2018-06-19 | 成都真实维度科技有限公司 | A kind of determining method of knee joint central point |
| CN108711187A (en) * | 2018-08-03 | 2018-10-26 | 华侨大学 | The method that registration fusion CT and MRI signal establish human body lumbar vertebrae three-dimensional simulation model |
| CN109389676A (en) * | 2018-10-22 | 2019-02-26 | 南方医科大学 | A kind of creation method based on Medical Imaging Technology rehabilitation brace Morphological Model |
| CN109528195A (en) * | 2018-11-08 | 2019-03-29 | 殷晓亮 | In a kind of MRI/CT image bone boundary and method for distinguishing is known with reference to bone mark object |
| CN110097546A (en) * | 2019-04-30 | 2019-08-06 | 北京大学第三医院 | A kind of method and device for assessing knee cartilage tissue defect |
| CN110169820A (en) * | 2019-04-24 | 2019-08-27 | 艾瑞迈迪科技石家庄有限公司 | A kind of joint replacement surgery pose scaling method and device |
| CN110232732A (en) * | 2019-05-08 | 2019-09-13 | 广州市红十字会医院（暨南大学医学院附属广州红十字会医院） | Method for making computer model of defect filling block of knee joint cartilage |
| CN110443839A (en) * | 2019-07-22 | 2019-11-12 | 艾瑞迈迪科技石家庄有限公司 | A kind of skeleton model spatial registration method and device |
| CN110613469A (en) * | 2019-09-18 | 2019-12-27 | 北京理工大学 | Automatic leg bone and lower limb force line detection method and device |
| WO2020087257A1 (en) * | 2018-10-30 | 2020-05-07 | 西安大医集团有限公司 | Image guidance method and device, and medical equipment and computer readable storage medium |
| CN111292363A (en) * | 2020-02-13 | 2020-06-16 | 张逸凌 | Joint image processing method and device and computing equipment |
| CN111292306A (en) * | 2020-02-04 | 2020-06-16 | 北京航空航天大学 | Knee joint CT and MR image fusion method and device |
| CN111311738A (en) * | 2020-03-04 | 2020-06-19 | 杭州市第三人民医院 | Ureter 3D digital-analog establishing method adopting imaging and data acquisition device thereof |
| CN111557732A (en) * | 2020-05-06 | 2020-08-21 | 重庆市九龙坡区中医院 | Foot arch image model reconstruction method |
| CN111640093A (en) * | 2020-05-20 | 2020-09-08 | 上海联影智能医疗科技有限公司 | Quality control method for medical image and computer-readable storage medium |
| CN111724468A (en) * | 2020-06-11 | 2020-09-29 | 青岛聚铄智能精准医疗有限公司 | Preparation method of conch holographic digital modeling |
| CN111951241A (en) * | 2020-08-06 | 2020-11-17 | 暨南大学 | Method for measuring and displaying muscle deformation of aquatic animals in movement process |
| CN111973270A (en) * | 2020-08-03 | 2020-11-24 | 上海交通大学 | Method for establishing individual knee joint bionic ligament biomechanical model |
| CN112155732A (en) * | 2020-09-29 | 2021-01-01 | 苏州微创畅行机器人有限公司 | Readable storage medium, bone modeling and registering system and bone surgery system |
| CN112200780A (en) * | 2020-09-29 | 2021-01-08 | 上海联影医疗科技股份有限公司 | Bone tissue positioning method, bone tissue positioning device, computer equipment and storage medium |
| CN112258640A (en) * | 2020-10-30 | 2021-01-22 | 李艳 | Skull model establishing method and device, storage medium and electronic equipment |
| CN112258638A (en) * | 2020-10-30 | 2021-01-22 | 李艳 | Human body model modeling method and device, storage medium and electronic equipment |
| CN112270705A (en) * | 2020-10-30 | 2021-01-26 | 李艳 | Skull slicing method, storage medium and electronic device |
| CN112288797A (en) * | 2020-10-30 | 2021-01-29 | 李艳 | Skull correction scheme generation system, construction method, acquisition method and device |
| CN112294352A (en) * | 2020-10-30 | 2021-02-02 | 李艳 | Skull slice acquisition method and device, storage medium and electronic equipment |
| CN112515764A (en) * | 2020-12-01 | 2021-03-19 | 江门市中心医院 | Intracranial tumor three-dimensional printing method, three-dimensional printing device and readable storage medium |
| CN113112900A (en) * | 2021-05-28 | 2021-07-13 | 北京大学第三医院（北京大学第三临床医学院） | Method for manufacturing meniscus 3D printing model for preoperative drilling |
| CN113192191A (en) * | 2021-05-28 | 2021-07-30 | 北京大学第三医院（北京大学第三临床医学院） | Meniscus three-dimensional model reconstruction method based on CT scanning |
| CN113205542A (en) * | 2021-04-30 | 2021-08-03 | 深圳高性能医疗器械国家研究院有限公司 | Three-dimensional image registration method using machine vision and blood brain barrier opening device |
| WO2021159519A1 (en) * | 2020-02-14 | 2021-08-19 | 西安大医集团股份有限公司 | Image guidance method and apparatus, radiotherapy device, and computer storage medium |
| CN113515875A (en) * | 2021-07-20 | 2021-10-19 | 北京工业大学 | Bone biomechanical modeling method, system and device based on multi-modal images |
| CN113616273A (en) * | 2021-08-04 | 2021-11-09 | 长安大学 | Positioning block manufacturing method and system for precise replacement of artificial knee joint |
| CN113689406A (en) * | 2021-08-24 | 2021-11-23 | 北京长木谷医疗科技有限公司 | Knee joint femoral posterior condylar point identification method and system based on motion simulation algorithm |
| CN113850810A (en) * | 2021-12-01 | 2021-12-28 | 杭州柳叶刀机器人有限公司 | Method and surgical system for correcting femur, storage medium, and electronic device |
| CN113952083A (en) * | 2021-09-29 | 2022-01-21 | 北京科技大学 | Biomechanically high-adaptation bionic artificial meniscus, molding system and molding method |
| CN114796901A (en) * | 2022-05-30 | 2022-07-29 | 北京大学第一医院 | Automatic drawing method, equipment and storage medium for lumbosacral nerve root |
| CN115272615A (en) * | 2022-08-10 | 2022-11-01 | 新乡学院 | Medical semantic segmentation method for femoral point cloud model and application of medical semantic segmentation method |
| CN115753795A (en) * | 2022-11-17 | 2023-03-07 | 上海智能制造功能平台有限公司 | Image denoising method and system based on optical measuring instrument and industrial CT machine and application |
| CN116306326A (en) * | 2023-05-25 | 2023-06-23 | 南方医科大学珠江医院 | Joint contact mechanics simulation model building method and device and electronic equipment |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102940530A (en) * | 2012-11-16 | 2013-02-27 | 昆明医科大学第一附属医院 | Method for virtually building anterior cruciate ligament on femur and tibia tunnels |
| CN102968565A (en) * | 2012-11-26 | 2013-03-13 | 昆明医科大学第一附属医院 | Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm |
-
2015
- 2015-07-23 CN CN201510437619.3A patent/CN105139442A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102940530A (en) * | 2012-11-16 | 2013-02-27 | 昆明医科大学第一附属医院 | Method for virtually building anterior cruciate ligament on femur and tibia tunnels |
| CN102968565A (en) * | 2012-11-26 | 2013-03-13 | 昆明医科大学第一附属医院 | Method for obtaining central axial line of anterior cruciate ligament based on linear fitting algorithm |
Non-Patent Citations (5)
| Title |
|---|
| 仝路 等: "三维数字化模型在膝关节修复重建中的应用进展", 《中国修复重建外科杂志》 * |
| 何川: "前交叉韧带虚拟重建系统的应用及术后膝关节内结构的三维有限元研究", 《中国博士学位论文全文数据库》 * |
| 王建平 等: "人体膝关节几何解剖模型仿真", 《系统仿真学报》 * |
| 陈欣 等: "基于二维轮廓序列的膝关节三维重建", 《计算机应用》 * |
| 黄赞 等: "基于MRI和CT二维图像重建膝关节三维数字化模型的股骨髁扭转角差异性研究", 《中国修复重建外科杂志》 * |
Cited By (76)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930617A (en) * | 2016-05-17 | 2016-09-07 | 南方医科大学 | Method for designing and forming stiffness-controllable bone tumor defect repair implant |
| CN106055848A (en) * | 2016-07-14 | 2016-10-26 | 哈尔滨理工大学 | Microstructure parameter-based individualized articular cartilage simulation method |
| CN106055848B (en) * | 2016-07-14 | 2017-10-13 | 哈尔滨理工大学 | One kind is based on microstructural parameter individuation articular cartilage emulation mode |
| CN106175802A (en) * | 2016-08-29 | 2016-12-07 | 吉林大学 | A kind of in body osteoarthrosis stress distribution detection method |
| CN106264731A (en) * | 2016-10-11 | 2017-01-04 | 昆明医科大学第附属医院 | A kind of method based on point-to-point registration technique virtual knee joint single condyle replacement model construction |
| CN107016666A (en) * | 2017-02-28 | 2017-08-04 | 四川大学华西医院 | Image post-processing method for pelvis CT three-dimensional reconstruction based on coordinate system |
| CN107016666B (en) * | 2017-02-28 | 2020-09-18 | 四川大学华西医院 | Image post-processing method for pelvis CT three-dimensional reconstruction based on coordinate system |
| CN106963487B (en) * | 2017-04-17 | 2021-08-13 | 四川大学华西医院 | Knee joint disc-shaped meniscus simulation operation method |
| CN106963487A (en) * | 2017-04-17 | 2017-07-21 | 四川大学华西医院 | Knee joint disc-shaped meniscus simulation operation method |
| CN108175378A (en) * | 2017-12-21 | 2018-06-19 | 成都真实维度科技有限公司 | A kind of determining method of knee joint central point |
| CN108143487A (en) * | 2017-12-21 | 2018-06-12 | 成都真实维度科技有限公司 | A kind of determining method of ankle-joint central point |
| CN107913100A (en) * | 2017-12-21 | 2018-04-17 | 成都真实维度科技有限公司 | A kind of definite method of late-segmental collapse point |
| CN108042217A (en) * | 2017-12-21 | 2018-05-18 | 成都真实维度科技有限公司 | A kind of definite method of three dimensions lower-limbs biology force-line |
| CN108711187A (en) * | 2018-08-03 | 2018-10-26 | 华侨大学 | The method that registration fusion CT and MRI signal establish human body lumbar vertebrae three-dimensional simulation model |
| CN108711187B (en) * | 2018-08-03 | 2021-12-31 | 华侨大学 | Method for establishing human lumbar vertebra three-dimensional simulation model by registering and fusing CT (computed tomography) and MRI (magnetic resonance imaging) signals |
| CN109389676A (en) * | 2018-10-22 | 2019-02-26 | 南方医科大学 | A kind of creation method based on Medical Imaging Technology rehabilitation brace Morphological Model |
| CN109389676B (en) * | 2018-10-22 | 2023-04-07 | 南方医科大学 | Medical image technology-based rehabilitation brace morphological model creation method |
| WO2020087257A1 (en) * | 2018-10-30 | 2020-05-07 | 西安大医集团有限公司 | Image guidance method and device, and medical equipment and computer readable storage medium |
| US11961245B2 (en) | 2018-10-30 | 2024-04-16 | Our United Corporation | Method and apparatus for performing image guidance, medical equipment, and computer-readable storage medium |
| CN111386555A (en) * | 2018-10-30 | 2020-07-07 | 西安大医集团股份有限公司 | Image guidance method and device, medical equipment and computer readable storage medium |
| CN111386555B (en) * | 2018-10-30 | 2024-03-29 | 西安大医集团股份有限公司 | Image guidance method and device, medical equipment and computer readable storage medium |
| CN109528195A (en) * | 2018-11-08 | 2019-03-29 | 殷晓亮 | In a kind of MRI/CT image bone boundary and method for distinguishing is known with reference to bone mark object |
| CN110169820A (en) * | 2019-04-24 | 2019-08-27 | 艾瑞迈迪科技石家庄有限公司 | A kind of joint replacement surgery pose scaling method and device |
| CN110097546A (en) * | 2019-04-30 | 2019-08-06 | 北京大学第三医院 | A kind of method and device for assessing knee cartilage tissue defect |
| CN110232732A (en) * | 2019-05-08 | 2019-09-13 | 广州市红十字会医院（暨南大学医学院附属广州红十字会医院） | Method for making computer model of defect filling block of knee joint cartilage |
| CN110443839A (en) * | 2019-07-22 | 2019-11-12 | 艾瑞迈迪科技石家庄有限公司 | A kind of skeleton model spatial registration method and device |
| CN110613469A (en) * | 2019-09-18 | 2019-12-27 | 北京理工大学 | Automatic leg bone and lower limb force line detection method and device |
| CN111292306A (en) * | 2020-02-04 | 2020-06-16 | 北京航空航天大学 | Knee joint CT and MR image fusion method and device |
| CN111292363B (en) * | 2020-02-13 | 2022-02-22 | 张逸凌 | Joint image processing method and device and computing equipment |
| CN111292363A (en) * | 2020-02-13 | 2020-06-16 | 张逸凌 | Joint image processing method and device and computing equipment |
| WO2021159519A1 (en) * | 2020-02-14 | 2021-08-19 | 西安大医集团股份有限公司 | Image guidance method and apparatus, radiotherapy device, and computer storage medium |
| US11883684B2 (en) | 2020-02-14 | 2024-01-30 | Our United Corporation | Image-guided method, radio therapy device, and computer storage medium |
| CN111311738B (en) * | 2020-03-04 | 2023-08-11 | 杭州市第三人民医院 | Ureter 3D digital-analog establishing method adopting imaging and data acquisition device thereof |
| CN111311738A (en) * | 2020-03-04 | 2020-06-19 | 杭州市第三人民医院 | Ureter 3D digital-analog establishing method adopting imaging and data acquisition device thereof |
| CN111557732B (en) * | 2020-05-06 | 2023-05-26 | 重庆市九龙坡区中医院 | Arch image model reconstruction method |
| CN111557732A (en) * | 2020-05-06 | 2020-08-21 | 重庆市九龙坡区中医院 | Foot arch image model reconstruction method |
| CN111640093A (en) * | 2020-05-20 | 2020-09-08 | 上海联影智能医疗科技有限公司 | Quality control method for medical image and computer-readable storage medium |
| CN111724468B (en) * | 2020-06-11 | 2024-05-28 | 青岛聚铄智能精准医疗有限公司 | Preparation method of conch holographic digital modeling |
| CN111724468A (en) * | 2020-06-11 | 2020-09-29 | 青岛聚铄智能精准医疗有限公司 | Preparation method of conch holographic digital modeling |
| CN111973270B (en) * | 2020-08-03 | 2022-02-22 | 上海涛影医疗科技有限公司 | Method for establishing individual knee joint bionic ligament biomechanical model |
| CN111973270A (en) * | 2020-08-03 | 2020-11-24 | 上海交通大学 | Method for establishing individual knee joint bionic ligament biomechanical model |
| CN111951241B (en) * | 2020-08-06 | 2023-07-18 | 暨南大学 | Method for measuring and displaying muscle deformation in aquatic animal exercise process |
| CN111951241A (en) * | 2020-08-06 | 2020-11-17 | 暨南大学 | Method for measuring and displaying muscle deformation of aquatic animals in movement process |
| CN112200780B (en) * | 2020-09-29 | 2023-09-29 | 上海联影医疗科技股份有限公司 | Bone tissue positioning method, device, computer equipment and storage medium |
| CN112155732A (en) * | 2020-09-29 | 2021-01-01 | 苏州微创畅行机器人有限公司 | Readable storage medium, bone modeling and registering system and bone surgery system |
| CN112155732B (en) * | 2020-09-29 | 2022-05-17 | 苏州微创畅行机器人有限公司 | Readable storage medium, bone modeling and registering system and bone surgery system |
| CN112200780A (en) * | 2020-09-29 | 2021-01-08 | 上海联影医疗科技股份有限公司 | Bone tissue positioning method, bone tissue positioning device, computer equipment and storage medium |
| WO2022089058A1 (en) * | 2020-10-30 | 2022-05-05 | 李艳 | Modeling method and apparatus for human body model, storage medium and electronic device |
| CN112294352B (en) * | 2020-10-30 | 2021-12-28 | 李艳 | Skull slice acquisition method and device, storage medium and electronic equipment |
| CN112288797A (en) * | 2020-10-30 | 2021-01-29 | 李艳 | Skull correction scheme generation system, construction method, acquisition method and device |
| CN112294352A (en) * | 2020-10-30 | 2021-02-02 | 李艳 | Skull slice acquisition method and device, storage medium and electronic equipment |
| CN112270705A (en) * | 2020-10-30 | 2021-01-26 | 李艳 | Skull slicing method, storage medium and electronic device |
| CN112258640B (en) * | 2020-10-30 | 2024-03-22 | 李艳 | Skull model building method and device, storage medium and electronic equipment |
| WO2022089055A1 (en) * | 2020-10-30 | 2022-05-05 | 李艳 | Method, apparatus, storage medium, and electronic device for obtaining cranial section |
| WO2022089054A1 (en) * | 2020-10-30 | 2022-05-05 | 李艳 | Method and apparatus for establishing skull model, storage medium, and electronic device |
| CN112258638A (en) * | 2020-10-30 | 2021-01-22 | 李艳 | Human body model modeling method and device, storage medium and electronic equipment |
| CN112270705B (en) * | 2020-10-30 | 2021-08-31 | 李艳 | Skull slicing method, storage medium and electronic device |
| CN112258640A (en) * | 2020-10-30 | 2021-01-22 | 李艳 | Skull model establishing method and device, storage medium and electronic equipment |
| CN112515764A (en) * | 2020-12-01 | 2021-03-19 | 江门市中心医院 | Intracranial tumor three-dimensional printing method, three-dimensional printing device and readable storage medium |
| CN113205542A (en) * | 2021-04-30 | 2021-08-03 | 深圳高性能医疗器械国家研究院有限公司 | Three-dimensional image registration method using machine vision and blood brain barrier opening device |
| CN113112900A (en) * | 2021-05-28 | 2021-07-13 | 北京大学第三医院（北京大学第三临床医学院） | Method for manufacturing meniscus 3D printing model for preoperative drilling |
| CN113192191A (en) * | 2021-05-28 | 2021-07-30 | 北京大学第三医院（北京大学第三临床医学院） | Meniscus three-dimensional model reconstruction method based on CT scanning |
| CN113192191B (en) * | 2021-05-28 | 2023-07-28 | 北京大学第三医院（北京大学第三临床医学院） | CT scanning-based meniscus three-dimensional model reconstruction method |
| CN113515875A (en) * | 2021-07-20 | 2021-10-19 | 北京工业大学 | Bone biomechanical modeling method, system and device based on multi-modal images |
| CN113515875B (en) * | 2021-07-20 | 2023-12-29 | 北京工业大学 | Bone biomechanics modeling method, system and device based on multi-modal image |
| CN113616273A (en) * | 2021-08-04 | 2021-11-09 | 长安大学 | Positioning block manufacturing method and system for precise replacement of artificial knee joint |
| CN113689406A (en) * | 2021-08-24 | 2021-11-23 | 北京长木谷医疗科技有限公司 | Knee joint femoral posterior condylar point identification method and system based on motion simulation algorithm |
| CN113952083B (en) * | 2021-09-29 | 2023-01-24 | 北京科技大学 | Biomechanically high-adaptation bionic artificial meniscus, molding system and molding method |
| CN113952083A (en) * | 2021-09-29 | 2022-01-21 | 北京科技大学 | Biomechanically high-adaptation bionic artificial meniscus, molding system and molding method |
| CN113850810A (en) * | 2021-12-01 | 2021-12-28 | 杭州柳叶刀机器人有限公司 | Method and surgical system for correcting femur, storage medium, and electronic device |
| CN114796901A (en) * | 2022-05-30 | 2022-07-29 | 北京大学第一医院 | Automatic drawing method, equipment and storage medium for lumbosacral nerve root |
| CN115272615A (en) * | 2022-08-10 | 2022-11-01 | 新乡学院 | Medical semantic segmentation method for femoral point cloud model and application of medical semantic segmentation method |
| CN115272615B (en) * | 2022-08-10 | 2023-09-08 | 新乡学院 | Femur point cloud model medical semantic segmentation method and application thereof |
| CN115753795A (en) * | 2022-11-17 | 2023-03-07 | 上海智能制造功能平台有限公司 | Image denoising method and system based on optical measuring instrument and industrial CT machine and application |
| CN116306326B (en) * | 2023-05-25 | 2023-09-15 | 南方医科大学珠江医院 | Joint contact mechanics simulation model building method and device and electronic equipment |
| CN116306326A (en) * | 2023-05-25 | 2023-06-23 | 南方医科大学珠江医院 | Joint contact mechanics simulation model building method and device and electronic equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105139442A (en) | Method for establishing human knee joint three-dimensional simulation model in combination with CT (Computed Tomography) and MRI (Magnetic Resonance Imaging) | |
| CN106264731B (en) | A method of based on the virtual knee joint single condyle displacement technique model construction of point-to-point registration technique | |
| US9681956B2 (en) | Acquiring and utilizing kinematic information for patient-adapted implants, tools and surgical procedures | |
| Fellows et al. | Magnetic resonance imaging for in vivo assessment of three-dimensional patellar tracking | |
| Bae et al. | Biomechanical analysis of the effects of medial meniscectomy on degenerative osteoarthritis | |
| Patel et al. | A three‐dimensional MRI analysis of knee kinematics | |
| Mizuno et al. | Q‐angle influences tibiofemoral and patellofemoral kinematics | |
| Penrose et al. | Development of an accurate three-dimensional finite element knee model | |
| CN108478250A (en) | Femur localization method, individuation osteotomy guide plate and the prosthese of total knee arthroplasty | |
| Jorge et al. | Finite element simulations of a hip joint with femoroacetabular impingement | |
| Monk et al. | The shape of the distal femur: a geometrical study using MRI | |
| CN112733404B (en) | Method for establishing human knee joint finite element model in motion process | |
| CA2753485A1 (en) | Customized orthopaedic implants and related methods | |
| CN107822745A (en) | The method of precision customized knee-joint prosthesis | |
| CN106963487A (en) | Knee joint disc-shaped meniscus simulation operation method | |
| Incavo et al. | Total knee arthroplasty using anatomic alignment can produce mid-flexion laxity | |
| CN108836457A (en) | A kind of 3D printing personalization Healing in Anterior Cruciate Ligament Reconstruction femoral bone tunnel locator and preparation method thereof | |
| Yamada et al. | In vivo movement analysis of the patella using a three-dimensional computer model | |
| Fernandez et al. | Integrating modelling, motion capture and x-ray fluoroscopy to investigate patellofemoral function during dynamic activity | |
| CN108056800A (en) | knee joint osteotomy tool and manufacturing system and manufacturing method thereof | |
| Maderbacher et al. | Reproducibility of navigation based kinematic analysis of the knee–A cadaveric investigation | |
| Xu et al. | Biomechanical strengthening design for limb articulation based on reconstructed skeleton kinesthetics | |
| Fitzpatrick et al. | A tibial-based coordinate system for three-dimensional data | |
| Zahra et al. | 3D designing and imaging process of the human knee joint: a review | |
| Lane | Human knee FEA model for transtibial amputee tibial cartilage pressure in gait and cycling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151209 |