3D lumbar spine intervertebral disc segmentation and compression simulation from MRI using shape-aware models.

Purpose: More accurate and robust image segmentations are needed for identification of spine pathologies and to assist with spine surgery planning and simulation. A framework for 3D segmentation of healthy and herniated intervertebral discs from T2-weighted magnetic resonance imaging was developed that exploits weak shape priors encoded in simplex mesh active surface models.

Methods: Weak shape priors inherent in simplex mesh deformable models have been exploited to automatically segment intervertebral discs. An ellipsoidal simplex template mesh was initialized within the disc image boundary through affine landmark-based registration and was allowed to deform according to image gradient forces. Coarse-to-fine multi-resolution approach was adopted in conjunction with decreasing shape memory forces to accurately capture the disc boundary. User intervention is allowed to turn off the shape feature and guide model deformation when the internal simplex shape memory influence hinders detection of pathology. A resulting surface mesh was utilized for disc compression simulation under gravitational and weight loads using Simulation Open Framework Architecture. For testing, 16 healthy discs were automatically segmented, and five pathological discs were segmented with minimal supervision.

Results: Segmentation results were validated against expert guided segmentation and demonstrate mean absolute shape distance error of <1 mm. Healthy intervertebral disc compression simulation resulted in a bulging disc under vertical pressure of 100 N/cm(2).

Conclusion: This study presents the application of a simplex active surface model featuring weak shape priors for 3D segmentation of healthy as well as herniated discs. A framework was developed that enables the application of shape priors in the healthy part of disc anatomy, with user intervention when the priors were inapplicable. The surface-mesh-based segmentation method is part of a processing pipeline for anatomical modelling to support interactive surgery simulation.