BackgroundClinical diagnosis and therapy for the lumbar disc herniation requires accurate vertebra segmentation. The complex anatomical structure and the degenerative deformations of the vertebrae makes its segmentation challenging.MethodsAn improved level set method, namely edge- and region-based level set method (ERBLS), is proposed for vertebra CT images segmentation. By considering the gradient information and local region characteristics of images, the proposed model can efficiently segment images with intensity inhomogeneity and blurry or discontinuous boundaries. To reduce the dependency on manual initialization in many active contour models and for an automatic segmentation, a simple initialization method for the level set function is built, which utilizes the Otsu threshold. In addition, the need of the costly re-initialization procedure is completely eliminated.ResultsExperimental results on both synthetic and real images demonstrated that the proposed ERBLS model is very robust and efficient. Compared with the well-known local binary fitting (LBF) model, our method is much more computationally efficient and much less sensitive to the initial contour. The proposed method has also applied to 56 patient data sets and produced very promising results.ConclusionsAn improved level set method suitable for vertebra CT images segmentation is proposed. It has the flexibility of segmenting the vertebra CT images with blurry or discontinuous edges, internal inhomogeneity and no need of re-initialization.

Nucleus replacement was deemed to have therapeutic potential for patients with intervertebral disc herniation. However, whether a patient would benefit from nucleus replacement is technically unclear. This study aimed to investigate the influence of nucleus pulposus (NP) removal on the biomechanical behavior of a lumbar motion segment and to further explore a computational method of biomechanical characteristics of NP removal, which can evaluate the mechanical stability of pulposus replacement. We, respectively, reconstructed three types of models for a mildly herniated disc and three types of models for a severely herniated disc based on a L4-L5 segment finite element model with computed tomography image data from a healthy adult. First, the NP was removed from the herniated disc models, and the biomechanical behavior of NP removal was simulated. Second, the NP cavities were filled with an experimental material (Poisson's ratio = 0.3; elastic modulus = 3 MPa), and the biomechanical behavior of pulposus replacement was simulated. The simulations were carried out under the five loadings of axial compression, flexion, lateral bending, extension, and axial rotation. The changes of the four biomechanical characteristics, i.e. the rotation degree, the maximum stress in the annulus fibrosus (AF), joint facet contact forces, and the maximum disc deformation, were computed for all models. Experimental results showed that the rotation range, the maximum AF stress, and joint facet contact forces increased, and the maximum disc deformation decreased after NP removal, while they changed in the opposite way after the nucleus cavities were filled with the experimental material.

Abstract. In this paper, we propose an approach to calculation and visualization of mechanical features for human lumbar L4-L5 segment. A finite-element analysis model of an intact lumbar L4-L5 segment is formed from a series of CT-scanning images. Axial compressive loads, which simulate the pressure from above, are applied to the FEA, and 10Nm moment are applied in flexion, extension, lateral bending and axial rotation to simulate bending. The simulation calculation shows the stress-strain distribution and deformation of the spine. The finite element method model can help surgeons analyze biomechanical characteristics of lumbar effectively and optimize individual therapy for orthopedic clinical.