GB2623864A - Method for bionic design of interbody fusion cage based on lumbar statistical shape model - Google Patents
Method for bionic design of interbody fusion cage based on lumbar statistical shape model Download PDFInfo
- Publication number
- GB2623864A GB2623864A GB2312651.9A GB202312651A GB2623864A GB 2623864 A GB2623864 A GB 2623864A GB 202312651 A GB202312651 A GB 202312651A GB 2623864 A GB2623864 A GB 2623864A
- Authority
- GB
- United Kingdom
- Prior art keywords
- lumbar
- shape model
- shape
- model
- statistical
- 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
- 230000004927 fusion Effects 0.000 title claims abstract description 48
- 238000013461 design Methods 0.000 title claims abstract description 38
- 239000011664 nicotinic acid Substances 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000010224 classification analysis Methods 0.000 claims abstract description 14
- 238000002591 computed tomography Methods 0.000 claims abstract description 9
- 238000007781 pre-processing Methods 0.000 claims abstract description 5
- 238000012549 training Methods 0.000 claims description 25
- 208000007623 Lordosis Diseases 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 6
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims description 2
- 230000000717 retained effect Effects 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 2
- 230000011218 segmentation Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 7
- 208000008035 Back Pain Diseases 0.000 description 3
- 208000008930 Low Back Pain Diseases 0.000 description 3
- 238000000605 extraction Methods 0.000 description 2
- 208000002193 Pain Diseases 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008468 bone growth Effects 0.000 description 1
- 230000007850 degeneration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 238000007639 printing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F2/00—Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
- A61F2/02—Prostheses implantable into the body
- A61F2/30—Joints
- A61F2/44—Joints for the spine, e.g. vertebrae, spinal discs
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61F—FILTERS IMPLANTABLE INTO BLOOD VESSELS; PROSTHESES; DEVICES PROVIDING PATENCY TO, OR PREVENTING COLLAPSING OF, TUBULAR STRUCTURES OF THE BODY, e.g. STENTS; ORTHOPAEDIC, NURSING OR CONTRACEPTIVE DEVICES; FOMENTATION; TREATMENT OR PROTECTION OF EYES OR EARS; BANDAGES, DRESSINGS OR ABSORBENT PADS; FIRST-AID KITS
- A61F2/00—Filters implantable into blood vessels; Prostheses, i.e. artificial substitutes or replacements for parts of the body; Appliances for connecting them with the body; Devices providing patency to, or preventing collapsing of, tubular structures of the body, e.g. stents
- A61F2/02—Prostheses implantable into the body
- A61F2/30—Joints
- A61F2/44—Joints for the spine, e.g. vertebrae, spinal discs
- A61F2/4455—Joints for the spine, e.g. vertebrae, spinal discs for the fusion of spinal bodies, e.g. intervertebral fusion of adjacent spinal bodies, e.g. fusion cages
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/02—Affine transformations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/772—Determining representative reference patterns, e.g. averaging or distorting patterns; Generating dictionaries
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/101—Computer-aided simulation of surgical operations
- A61B2034/105—Modelling of the patient, e.g. for ligaments or bones
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/108—Computer aided selection or customisation of medical implants or cutting guides
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/16—Customisation or personalisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/10—Additive manufacturing, e.g. 3D printing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
- G06T2207/30012—Spine; Backbone
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/03—Recognition of patterns in medical or anatomical images
- G06V2201/033—Recognition of patterns in medical or anatomical images of skeletal patterns
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Multimedia (AREA)
- Artificial Intelligence (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Heart & Thoracic Surgery (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Veterinary Medicine (AREA)
- Neurology (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- Computer Hardware Design (AREA)
- Transplantation (AREA)
- Surgery (AREA)
- Cardiology (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Molecular Biology (AREA)
- Mathematical Analysis (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
Abstract
A method for bionic design of an interbody fusion cage based on a lumbar statistical shape model includes medical computer tomography (CT) data segmentation and preprocessing, automated construction of a lumbar statistical shape model database, classification analysis of the statistical shape model database to form shape model sub-databases in order to obtain an average shape model for each sub-databse, and bionic design of an interbody fusion cage. According to the method, the lumbar statistical shape model database can be automatically generated, and a large number of 3D lumbar shape models can be automatically and quickly generated. An interbody fusion cage that can meet the requirements of a larger group range and better conforms to a shape of an individual lumbar vertebra is obtained by extracting lumbar shape features and performing the classification analysis and the bionic design.
Description
METHOD FOR BIONIC DESIGN OF INTERBODY FUSION CAGE BASED ON
LUMBAR STATISTICAL SHAPE MODEL
TECHNICAL FIELD
[0001] The present disclosure relates to the field of bionic implant device design, and in particular, to a method for bionic design of an interbody fusion cage based on a lumbar statistical shape model.
BACKGROUND
[0002] Low back pain is a common lumbar disease, which is generally caused by lumbar degeneration and changes in lumbar anatomical structure and mechanical characteristics caused by specific occupational habits.
[0003] Interbody fusion is currently the most common and effective way to treat low back pain. An interbody fusion cage is used in the interbody fusion to perfomi vertebral fusion to maintain lumbar stability and relieve pain of patients with low back pain, and provides an important function in supporting bone growth and maintaining a normal physiological height of a lumbar vertebra [0004] Literature research shows that current development of the interbody fusion cage focuses on reducing the subsidence of the interbody fusion cage after implantation by adjusting a porosity and a material of the interbody fusion cage, so as to achieve a better fusion effect [0005] Clinically, an interbody fusion cage with a suitable dimension is usually selected from a limited number of commercially available interbody fusion cages by an experienced orthopedic surgeon according to a patient's medical images.
[0006] Due to great individual differences, interbody fusion cages clinically widely used at present cannot well match upper and lower vertebral endplates and vertebrae of patients in shape, which leads to uneven stress on the upper and lower endplates of the patients' lumbar vertebrae, thereby affecting the fusion effect and easily leading to problems such as subsidence after the fusion [0007] In recent years, some researchers have customized interbody fusion cages based on specific patients by combining a three-dimensional (3D) printing technology, so as to meet requirements that an interbody fusion cage fits with upper and lower vertebral endplates of a specific patient in shape. However, the steps of customizing an interbody fusion cage in this method are complicated and labor and time costs are high, so that it is difficult to apply the method in a large scale [0008] The statistical shape model is an effective method for describing shape changes of an object. The statistical shape model is established based on the principle that shape changes of an object satisfy Gaussian distribution when there are large samples of the object, and a changing law of the object shape may be described by combining an average shape model with a changing mode.
[0009] The statistical shape model enables a large number of object shape models to be effectively and quickly formed, and facilitates extraction of characteristics of object shape changes, thereby providing a new design idea for designing an interbody fusion cage that meets requirements of large-scale groups and conforms to the shape of an individual lumbar vertebra
SUMMARY
100101 Against the technical problem that it is difficult to use a current interbody cage design to meet requirements of a larger group and conform to a shape of an individual lumbar vertebra, the present disclosure provides a method for bionic design of an interbody fusion cage based on a lumbar statistical shape model.
100111 According to the method, a lumbar statistical shape model database can be effectively and automatically constructed, and bionic design of an interbody fusion cage can be performed based on the lumbar statistical shape model database, so as to design a bionic interbody fusion cage that meets requirements of a larger group and better conforms to a shape of an individual lumbar vertebra.
100121 The specific technical solutions of the present disclosure are as follow: 10013] As shown in FIG. 1, a method for bionic design of an interbody fusion cage based on a lumbar statistical shape model includes the following steps: 100141 1) data preprocessing: segmenting N human lumbar vertebrae from medical computer tomography (CT) data, reconstructing a 3D lumbar shape model, outputting a standard template library (STL) file, and simplifying and optimizing a grid of the 3D lumbar shape model as a lumbar training set; [0015] 2) automatically building a lumbar statistical shape model database by applying a relevant algorithm and writing program, specifically including the following substeps: [0016] 2.01) randomly selecting one 3D lumbar shape model from the training set as a template model, and aligning the remaining 3D lumbar shape models with the template model by means of an iterative closest point algorithm; 100171 2.02) causing two point sets to coincide to a largest extent by rotation, translation and other transformation in the algorithm, so as to reduce shape errors caused by the rotation and the translation; f [0018] 2.03) assuming that a to-be-aligned point set s and a template point set is calculating whether a result of an objective function reaches a threshold as an index for determining whether to stop iteration, where M and X are point sets, and are points in the point sets, T is a translation matrix, R is a rotation matrix, and 1\ta1 is a number of points in the to-be-aligned point set; 100191 2.04) in order to unify 3D lumbar shape models with different numbers of points into 3D lumbar shape models with the same number of points and establish a corresponding relationship between points, registering the template model to the aligned 3D lumbar shape model by means of a non-rigid iterative closest point algorithm; [0020] 2.05) adding affine transformation to rigid transformation in the algorithm, so that the template model is deformed onto the to-be-registered 3D lumbar shape model to complete the registration process, establishing a corresponding relationship between the template model and points of the to-be-registered 3D lumbar shape model, and establishing a new lumbar training set; [0021] 2.06) analyzing a principal component of the registered lumbar training set; 100221 2.07) calculating an average shape model of the lumbar training set: , where n represents the number of 3D shape models in the lumbar training set, and X represents the 3D lumbar shape models; [0023], 2.08) calculating a covariance matrix of the lumbar training set: [0024] 2.09) obtaining an eigenvector and eigenvalue Sy, =1,(p, of the covariance matrix, where Ai represents the eigenvalue, (pi represents the eigenvector, the eigenvector represents a main pattern of a shape change, the eigenvalue represents a variance of the principal component, and a larger eigenvalue represents a greater amount of data retained in a direction corresponding to the eigenvector;
X
[0025] 2.10) expressing the lumbar statistical shape model database as where t represents sorting the eigenvalues in descending order, and for first t eigenyalues, a shape parameter b, is independent, obeys Gaussian distribution of (0, X.,), is used to control a change range of the shape, and is generally; and [0026] 2.11) executing a compiled automated program to automatically generate the lumbar statistical shape model database, 100271 3) performing classification analysis on the lumbar statistical shape model database to form lumbar statistical shape model sub-databases, and obtaining an average shape model for each of the sub-databases: 100281 3.01) obtaining a median sagittal plane of each shape model in the lumbar statistical shape model database by means of a mirror image method; [0029] 3.02) mirroring the 3D lumbar shape model in any sagittal plane in the mirror image method to obtain a new model, aligning the new model with the original model by means of the iterative closest point algorithm, and obtaining a symmetry plane with the aligned two models as a whole, that is, a median sagittal plane; [0030] 3.03) extracting required lumbar anatomical features, specifically, measuring an intervertebral disc height and a segmental lordosis angle on the sagittal plane of each 3D lumbar shape model in the lumbar statistical shape model database; [0031] 3.04) performing classification analysis on the lumbar statistical shape model database based on the intervertebral disc height and the segmental lordosis angle, and classifying 3D lumbar shape models with different intervertebral disc heights and different segmental lordosis angles; and [0032] 3.05) after the classification, forming a plurality of sub-databases, and obtaining an r----average shape model for each of the sub-databases, that is, , as a 3D lumbar shape model for bionic design, where j is a number of 3D lumbar shape models in each sub-database, and Xi represents the 3D lumbar shape models in the sub-databases; and [0033] 4) performing bionic design of an interbody fusion cage according to the average shape model for each of the sub-databases: 100341 4.01) extracting a design surface required for design, and extracting an upper endplate surface and a lower endplate surface of two adjacent vertebrae of the average shape model for each of the sub-databases; and [0035] 4.02) performing bionic design according to the extracted design surface to form a series of interbody fusion cages with bionic shapes.
[0036] The method according to the present disclosure mainly includes four processes: medical CT data segmentation and preprocessing, automated construction of a lumbar statistical shape model database, classification analysis of the statistical shape model database, and bionic design of an interbody fusion cage.
[0037] According to the present disclosure, collected medical CT images are first preprocessed, and the lumbar statistical shape model database is automatically established by means of alignment, registration, principal component analysis and other processes [0038] Classification analysis is performed on the established lumbar statistical shape model database to form a plurality of sub-databases, and an average shape model is obtained for each of the sub-databases.
100391 The bionic design of an interbody fusion cage is performed according to the obtained lumbar average shape model for each of the sub-databases.
[0040] The present disclosure has following beneficial effects: [0041] 1. The lumbar statistical shape model database is automatically established, and a large number of 3D lumbar shape models can be automatically and quickly generated, which facilitates extraction of lumbar shape features, and effectively reduces time and labor costs. 100421 2. Based on the lumbar statistical shape model, by classification analysis and bionic design, an interbody fusion cage that can meet requirements of a larger group range and better conforms to a shape of an individual lumbar vertebra is designed, which significantly improves fit between the fusion cage and a patient's vertebra, and further helps to reduce the risk of loosening and subsidence.
BRIEF DESCRIPTION OF THE DRAWINGS
[0043] FIG. t is a flowchart of the present disclosure; [0044] FIG. 2 is an effect diagram after alignment of two L4-L5 lumbar vertebrae by means of an iterative closest point algorithm and registration thereof by means of a non-rigid iterative closest point algorithm; [0045] FIG. 3 is an effect diagram of changes under a first principal component of L4-L5 lumbar vertebrae; [0046] FIG. 4 is a schematic diagram of measuring an intervertebral disc height and a segmental lordosis angle on a median sagittal plane of an L4-L5 3D lumbar shape model, and [0047] FIG. 5 is a schematic diagram of designing an anterior interbody fusion cage with a bionic shape based on a sub-database.
DETAILED DESCRIPTION OF THE EMBODIMENTS
[0048] The following further describes specific implementations of the present disclosure in detail with reference to accompanying drawings [0049] In a specific implementation, taking L4-L5 vertebrae as an example, a statistical shape model thereof is established, and an anterior interbody fusion cage is designed based on the statistical shape model [0050] Medical CT images of 30 volunteers' lumbar vertebrae were segmented, an L4-L5 3D lumbar shape model was reconstructed, and an STL format file was outputted.
[0051] The L4-L5 3D lumbar shape model was further processed, and grids were simplified and optimized to keep a number of vertices between 4000 and 5000, so as to improve efficiency of a subsequent algorithm.
100521 Through preprocessing, a training set containing 30 3D shape model samples of L4-L5 lumbar vertebrae was obtained.
[0053] The 3D shape model samples of L4-L5 lumbar vertebrae in the training set might be expressed as follows: [0054] X = (xii, vii, zii, xin, yin, zit), where xin, yin and zin represent coordinates of points of the samples.
100551 The L4-L5 lumbar training set might be expressed as follows: X = (Xi, X2, ..., X:40). 100561 One L4-L5 3D lumbar shape model Xi was randomly selected from the L4-L5 3D lumbar shape model training set as a template model, and the remaining L4-L5 3D lumbar shape models were aligned with the template model by means of an iterative closest point algorithm. [0057] This process was to minimize position errors caused by different people and different photographing environments.
[0058] The training set of the aligned L4-L5 3D lumbar shape models was updated to: X = [0059] After the alignment, the template model Xi was registered to the training set of the aligned L4-L5 3D lumbar shape models by means of a non-rigid iterative closest point algorithm. [0060] In this process, by registration by means of the non-rigid iterative closest point algorithm, the training set of the L4-L5 3D lumbar shape models with the same number of points and corresponding points was obtained.
[0061] After the registration, the training set of the L4-L5 3D lumbar shape models was updated to: X = (Xi, X2, ..., X30).
[0062] FIG. 2 is an effect diagram showing alignment and registration of one L4-L5 3D lumbar shape model and a template model.
[0063] A principal component of the training set of the L4-L5 3D lumbar shape models is analyzed: [0064] An average shape model of the training set of the L4-L5 3D lumbar shape models is calculated as follows.
100651 A covariance matrix of the training set of the L4-L5 3D lumbar shape models is calculated as follows: [0066] [0067] An eigenvalue and an eigenvector of the covariance matrix are calculated as follows: Scpi = ki(pi, where X represents the eigenvalue of the covariance matrix, eigenvalues are sequentially arranged from large to small, and (pi represents the eigenvector.
[00681 [0069] where m is a number of change patterns, k represents first k change patterns, and P represents a probability of the first k change patterns in the total change patterns, with a value of P 95%.
100701 In this example, when k = 8, the probability was 96.7%, which was greater than 95%, so it was set that k = 8.
100711 An L4-L5 lumbar statistical shape model database might be expressed as follows: = [0072] A program for automatically establishing a lumbar statistical shape database is run to automatically generate an L4-L5 lumbar statistical shape model database.
[0073] FIG. 3 is an effect diagram of changes under a first principal component of L4-L5 lumbar vertebrae.
[0074] Classification analysis is performed on the obtained L4-L5 lumbar statistical shape model database.
[0075] This process involves measurement of an intervertebral disc height and a segmental lordosis angle on a median sagittal plane of each L4-L5 3D lumbar shape model.
[0076] A median sagittal plane of the L4-L5 3D lumbar shape model is obtained by means of a mirror image method.
100771 The L4-L5 3D lumbar shape model needs to be mirrored in any sagittal plane in the mirror image method to obtain a new model, the new model is aligned with the original model by means of the iterative closest point algorithm, and a symmetry plane with the aligned two models as a whole, that is, a median sagittal plane is obtained.
100781 Required L4-L5 lumbar anatomical features are extracted, specifically, an intervertebral disc height and a segmental lordosis angle are measured on the median sagittal plane of each L4-L5 3D lumbar shape model, and classification analysis is performed on the measured L4-L5 lumbar statistical shape model database based on the intervertebral disc height and the segmental lordosis angle, to obtain different sub-databases.
[0079] The different sub-databases differ in different intervertebral disc heights and segmental lordosis angles, and an average shape model is obtained for each sub-database.
[0080] FIG. 4 is a schematic diagram of measuring an intervertebral disc height and a segmental lordosis angle on a median sagittal plane of an L4-L5 3D lumbar shape model [0081] The bionic design of an interbody fusion cage is performed according to the obtained average shape model for the sub-database, and a lower endplate surface of an L4 vertebra and an upper endplate surface of an LS vertebra are extracted as design surfaces of the interbody fusion cage 100821 Based on the design surfaces, the bionic design of the interbody fusion cage is performed, and an anterior interbody fusion cage that can meet requirements of a larger group range and better conforms to a shape of an individual lumbar vertebra is designed.
[0083] FIG. 5 shows an anterior interbody fusion cage with a bionic shape designed based on an average shape model of a sub-database.
[0084] In conclusion, according to the present disclosure, the lumbar statistical shape model can be automatically established, the bionic design of the interbody fusion cage is performed based on the lumbar statistical shape model through classification analysis, and specific examples illustrate the principle and effects of the present disclosure.
Claims (4)
- WHAT IS CLAIMED IS: I. A method for bionic design of an interbody fusion cage based on a lumbar statistical shape model, comprising the following steps: 1) data preprocessing: segmenting N human lumbar vertebrae from medical computer tomography (CT) data, reconstructing a three-dimensional (3D) lumbar shape model, outputting a standard template library (STL) file, and simplifying and optimizing a grid of the 3D lumbar shape model as a lumbar training set; 2) automatically building a lumbar statistical shape model database by applying a relevant algorithm and writing program; 3) performing classification analysis on the lumbar statistical shape model database to form lumbar statistical shape model sub-databases, and obtaining an average shape model for each of the sub-databases; and 4) performing bionic design of an interbody fusion cage according to the average shape model for each of the sub-databases.2. The method for bionic design of an interbody fusion cage based on a lumbar statistical shape model according to claim 1, wherein the automatically building a lumbar statistical shape model database in step 2) comprises the following substeps: 2.01) randomly selecting one 3D lumbar shape model from the training set as a template model, and aligning the remaining 3D lumbar shape models with the template model by means of an iterative closest point algorithm; 2.02) causing two point sets to coincide to a largest extent by rotation, translation and other transformation in the algorithm, so as to reduce shape errors caused by the rotation and the translation; 2.04 assuming that a to-be-aligned point set is and a template point set is calculating whether a result of an objective function reaches a threshold as an index for determining whether to stop iteration, wherein M and X are point sets, and xj are points in the point sets, T is a translation matrix, R is a rotation matrix, and N. is a number of points in the template point set; 2.04) in order to unify 3D lumbar shape models with different numbers of points into 3D lumbar shape models with the same number of points and establish a corresponding relationship between points, registering the template model to the aligned 3D lumbar shape model by means of a non-rigid iterative closest point algorithm, 2 05) adding affine transformation to rigid transformation in the algorithm, so that the template model is deformed onto the to-be-registered 3D lumbar shape model to complete the registration process, establishing a corresponding relationship between the template model and points of the to-be-registered 3D lumbar shape model, and establishing a new lumbar training set; 2.06) analyzing a principal component of the registered lumbar training set; 2.07) calculating an average shape model of the lumbar training set: , wherein n represents the number of 3D shape models in the lumbar training set, and X represents the 3D lumbar shape models; 2.08) calculating a covariance matrix of the lumbar training set.
- 2.09) obtaining an eigenvector and eigenvalue ST; = kipi of the covariance matrix, wherein Xi represents the eigenvalue, (pi represents the eigenvector, the eigenvector represents a main pattern of a shape change, the eigenvalue represents a variance of the principal component, and a larger eigenvalue represents a greater amount of data retained in a direction corresponding to the eigenvector; 2.10) expressing the lumbar statistical shape model database as wherein t represents sorting the eigenvalues in descending order, and for first t eigenvalues, a shape parameter bi is independent, obeys Gaussian distribution of (0, Ai), is used to control a change range of the shape, and is generally and 2.11) executing a compiled automated program to automatically generate the lumbar statistical shape model database.
- 3. The method for bionic design of an interbody fusion cage based on a lumbar statistical shape model according to claim 1, wherein the performing classification analysis on the lumbar statistical shape model database to form lumbar statistical shape model sub-databases, and obtaining an average shape model for each of the sub-databases in step 3) comprises the following sub step s: 3.01) obtaining a median sagittal plane of each shape model in the lumbar statistical shape model database by means of a mirror image method; 3.02) mirroring the 3D lumbar shape model in any sagittal plane in the mirror image method to obtain a new model, aligning the new model with the original model by means of the iterative closest point algorithm, and obtaining a symmetry plane with the aligned two models as a whole, that is, a median sagittal plane; 3.03) extracting required lumbar anatomical features, specifically, measuring an 10 intervertebral disc height and a segmental lordosis angle on the sagittal plane of each 3D lumbar shape model in the lumbar statistical shape model database; 3.04) performing classification analysis on the lumbar statistical shape model database based on the intervertebral disc height and the segmental lordosis angle, and classifying 3D lumbar shape models with different intervertebral disc heights and different segmental lordosis angles; and 3.05) after the classification analysis, forming a plurality of sub-databases, and obtaining an average shape model for each of the sub-databases, that is, , as a 3D lumbar shape model for bionic design, wherein j is a number of 3D lumbar shape models in each sub-database, and X, represents the 3D lumbar shape models in the sub-databases
- 4. The method for bionic design of an interbody fusion cage based on a lumbar statistical shape model according to claim 1, wherein the performing bionic design of an interbody fusion cage according to the average shape model for each of the sub-databases in step 4) comprises the following steps: 4.01) extracting a design surface required for design, and extracting an upper endplate surface and a lower endplate surface of two adjacent vertebrae of the average shape model for each of the sub-databases; and 4.02) performing bionic design according to the extracted design surface to form a series of interbody fusion cages with bionic shapes.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211306517.4A CN115630455A (en) | 2022-10-25 | 2022-10-25 | Intervertebral fusion cage bionic design method based on lumbar vertebra statistical shape model |
Publications (2)
Publication Number | Publication Date |
---|---|
GB202312651D0 GB202312651D0 (en) | 2023-10-04 |
GB2623864A true GB2623864A (en) | 2024-05-01 |
Family
ID=84906412
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
GB2312651.9A Pending GB2623864A (en) | 2022-10-25 | 2023-08-18 | Method for bionic design of interbody fusion cage based on lumbar statistical shape model |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN115630455A (en) |
GB (1) | GB2623864A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102208117A (en) * | 2011-05-04 | 2011-10-05 | 西安电子科技大学 | Method for constructing vertebral three-dimensional geometry and finite element mixture model |
CN113284243A (en) * | 2021-03-02 | 2021-08-20 | 上海交通大学 | Method and system for establishing statistical shape model of three-dimensional surface model, storage medium and terminal |
AU2021107289A4 (en) * | 2021-08-25 | 2021-12-09 | 3Dmorphic Pty Ltd | A Method and Medical Implant |
US20220054197A1 (en) * | 2019-05-09 | 2022-02-24 | Materialise N.V. | Surgery planning system with automated defect quantification |
CN114648612A (en) * | 2022-04-15 | 2022-06-21 | 南京大学 | Spine three-dimensional reconstruction method and system |
-
2022
- 2022-10-25 CN CN202211306517.4A patent/CN115630455A/en active Pending
-
2023
- 2023-08-18 GB GB2312651.9A patent/GB2623864A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102208117A (en) * | 2011-05-04 | 2011-10-05 | 西安电子科技大学 | Method for constructing vertebral three-dimensional geometry and finite element mixture model |
US20220054197A1 (en) * | 2019-05-09 | 2022-02-24 | Materialise N.V. | Surgery planning system with automated defect quantification |
CN113284243A (en) * | 2021-03-02 | 2021-08-20 | 上海交通大学 | Method and system for establishing statistical shape model of three-dimensional surface model, storage medium and terminal |
AU2021107289A4 (en) * | 2021-08-25 | 2021-12-09 | 3Dmorphic Pty Ltd | A Method and Medical Implant |
CN114648612A (en) * | 2022-04-15 | 2022-06-21 | 南京大学 | Spine three-dimensional reconstruction method and system |
Non-Patent Citations (1)
Title |
---|
CLOGENSON MARINE; DUFF JOHN M.; LUETHI MARCEL; LEVIVIER MARC; MEULI RETO; BAUR CHARLES; HENEIN SIMON: "A statistical shape model of the human second cervical vertebra", INTERNATIONAL JOURNAL OF COMPUTER ASSISTED RADIOLOGY AND SURGERY, SPRINGER, DE, vol. 10, no. 7, 30 October 2014 (2014-10-30), DE , pages 1097 - 1107, XP035506227, ISSN: 1861-6410, DOI: 10.1007/s11548-014-1121-x * |
Also Published As
Publication number | Publication date |
---|---|
CN115630455A (en) | 2023-01-20 |
GB202312651D0 (en) | 2023-10-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10068671B2 (en) | Methods and systems for producing an implant | |
US20180189421A1 (en) | Producing a three-dimensional model of an implant | |
WO2023142956A1 (en) | Total hip replacement preoperative planning system based on deep learning | |
US20130018497A1 (en) | Methods and systems for producing an implant | |
CN112037200A (en) | Method for automatically identifying anatomical features and reconstructing model in medical image | |
CN114119872B (en) | Method for analyzing 3D printing intraspinal plants based on artificial intelligence big data | |
CN108711187A (en) | The method that registration fusion CT and MRI signal establish human body lumbar vertebrae three-dimensional simulation model | |
CN114998301B (en) | Vertebral body sub-region segmentation method and device and storage medium | |
Bruner et al. | Landmark‐based shape analysis of the archaic Homo calvarium from Ceprano (Italy) | |
Urciuoli et al. | A comparative analysis of the vestibular apparatus in Epipliopithecus vindobonensis: Phylogenetic implications | |
GB2623864A (en) | Method for bionic design of interbody fusion cage based on lumbar statistical shape model | |
CN116228720A (en) | Medical image processing method, system, equipment and medium | |
CN109035404A (en) | A kind of human foot model positioning transformations method and system | |
CN113506301A (en) | Tooth image segmentation method and device | |
Qian et al. | Attention-based Shape-Deformation Networks for Artifact-Free Geometry Reconstruction of Lumbar Spine from MR Images | |
CN113506302B (en) | Interactive object updating method, device and processing system | |
Güldeniz et al. | A Semi-automatic Parameter Extraction Process for Modeling of Lumbar Vertebral Body Segments | |
Li et al. | Toward patient-specific computational study of aortic diseases: a population based shape modeling approach | |
CN117338488A (en) | Construction method of cervical vertebra model and printing method of cervical vertebra prosthesis | |
CN118493858A (en) | 3D printing modeling design method and preparation method for bone defect filling body | |
Gorrila | Data Driven Dense 3D Facial Reconstruction From 3D Skull Shape | |
CN113936100A (en) | Extraction and reconstruction method for human knee joint cruciate ligament insertion points | |
CN118691632A (en) | Automatic airway segmentation result post-processing method based on oral CBCT image | |
CN118135313A (en) | Spine medical image key point detection system and method based on deep learning | |
CN118397216A (en) | Bone material non-uniform assignment method based on slope difference distribution |