Automated generation of curved planar reformations from MR images of the spine

Volumetric Parameterization of the Placenta to a Flattened Template

We present a volumetric mesh-based algorithm for parameterizing the placenta to a flattened template to enable effective visualization of local anatomy and function. MRI shows potential as a research tool as it provides signals directly related to placental function. However, due to the curved and highly variable in vivo shape of the placenta, interpreting and visualizing these images is difficult. We address interpretation challenges by mapping the placenta so that it resembles the familiar ex vivo shape. We formulate the parameterization as an optimization problem for mapping the placental shape represented by a volumetric mesh to a flattened template. We employ the symmetric Dirichlet energy to control local distortion throughout the volume. Local injectivity in the mapping is enforced by a constrained line search during the gradient descent optimization. We validate our method using a research study of 111 placental shapes extracted from BOLD MRI images. Our mapping achieves sub-voxel accuracy in matching the template while maintaining low distortion throughout the volume. We demonstrate how the resulting flattening of the placenta improves visualization of anatomy and function. Our code is freely available at https://github.com/mabulnaga/placenta-flattening.

The third dimension of scoliosis: The forgotten axial plane

Idiopathic scoliosis is a three-dimensional (3D) deformity of the spine. In clinical practice, however, the diagnosis and treatment of scoliosis consider only two dimensions (2D) as they rely solely on postero-anterior (PA) and lateral radiographs. Thus, the projections of the deformity are evaluated in only the coronal and sagittal planes, whereas those in the axial plane are disregarded, precluding an accurate assessment of the 3D deformity. A universal dogma in engineering is that designing a 3D object requires drawing projections of the object in all three planes. Similarly, when dealing with a 3D deformity, knowledge of the abnormalities in all three planes is crucial, as each plane is as important as the other two planes. This article reviews the chronological development of axial plane imaging and spinal deformity measurement.

VolHOG: a volumetric object recognition approach based on bivariate histograms of oriented gradients for vertebra detection in cervical spine MRI

PURPOSE

The automatic recognition of vertebrae in volumetric images is an important step toward automatic spinal diagnosis and therapy support systems. There are many applications such as the detection of pathologies and segmentation which would benefit from automatic initialization by the detection of vertebrae. One possible application is the initialization of local vertebral segmentation methods, eliminating the need for manual initialization by a human operator. Automating the initialization process would optimize the clinical workflow. However, automatic vertebra recognition in magnetic resonance (MR) images is a challenging task due to noise in images, pathological deformations of the spine, and image contrast variations.

METHODS

This work presents a fully automatic algorithm for 3D cervical vertebra detection in MR images. We propose a machine learning method for cervical vertebra detection based on new features combined with a linear support vector machine for classification. An algorithm for bivariate gradient orientation histogram generation from three-dimensional raster image data is introduced which allows us to describe three-dimensional objects using the authors' proposed bivariate histograms.

RESULTS

A detailed performance evaluation on 21 T2-weighted MR images of the cervical vertebral region is given. A single model for cervical vertebrae C3-C7 is generated and evaluated. The results show that the generic model performs equally well for each of the cervical vertebrae C3-C7. The algorithm's performance is also evaluated on images containing various levels of artificial noise. The results indicate that the proposed algorithm achieves good results despite the presence of severe image noise.

CONCLUSIONS

The proposed detection method delivers accurate locations of cervical vertebrae in MR images which can be used in diagnosis and therapy. In order to achieve absolute comparability with the results of future work, the authors are following an open data approach by making the image dataset used in their performance evaluation available to the public.

Automated detection, 3D segmentation and analysis of high resolution spine MR images using statistical shape models

Recent advances in high resolution magnetic resonance (MR) imaging of the spine provide a basis for the automated assessment of intervertebral disc (IVD) and vertebral body (VB) anatomy. High resolution three-dimensional (3D) morphological information contained in these images may be useful for early detection and monitoring of common spine disorders, such as disc degeneration. This work proposes an automated approach to extract the 3D segmentations of lumbar and thoracic IVDs and VBs from MR images using statistical shape analysis and registration of grey level intensity profiles. The algorithm was validated on a dataset of volumetric scans of the thoracolumbar spine of asymptomatic volunteers obtained on a 3T scanner using the relatively new 3D T2-weighted SPACE pulse sequence. Manual segmentations and expert radiological findings of early signs of disc degeneration were used in the validation. There was good agreement between manual and automated segmentation of the IVD and VB volumes with the mean Dice scores of 0.89 ± 0.04 and 0.91 ± 0.02 and mean absolute surface distances of 0.55 ± 0.18 mm and 0.67 ± 0.17 mm respectively. The method compares favourably to existing 3D MR segmentation techniques for VBs. This is the first time IVDs have been automatically segmented from 3D volumetric scans and shape parameters obtained were used in preliminary analyses to accurately classify (100% sensitivity, 98.3% specificity) disc abnormalities associated with early degenerative changes.

Excitation and geometrically matched local encoding of curved slices

In this work, the concept of excitation and geometrically matched local in-plane encoding of curved slices (ExLoc) is introduced. ExLoc is based on a set of locally near-orthogonal spatial encoding magnetic fields, thus maintaining a local rectangular shape of the individual voxels and avoiding potential problems arising due to highly irregular voxel shapes. Unlike existing methods for exciting curved slices based on multidimensional radiofrequency-pulses, excitation and geometrically matched local encoding of curved slices does not require long duration or computationally expensive radiofrequency-pulses. As each encoding field consists of a superposition of potentially arbitrary (spatially linear or nonlinear) magnetic field components, the resulting field shape can be adapted with high flexibility to the specific region of interest. For extended nonplanar structures, this results in improved relevant volume coverage for fewer excited slices and thus increased efficiency. In addition to the mathematical description for the generation of dedicated encoding fields and data reconstruction, a verification of the ExLoc concept in phantom experiments and examples for in vivo curved single and multislice imaging are presented.

Classification of three-dimensional thoracic deformities in adolescent idiopathic scoliosis from a multivariate analysis

PURPOSE

Understanding how to classify and quantify three-dimensional (3D) spinal deformities remains an open question in adolescent idiopathic scoliosis. The objective of this study was to perform a 3D manifold characterization of scoliotic spines demonstrating thoracic deformations using a novel geometric and intuitive statistical tool to determine patterns in pathological cases.

METHODS

Personalized 3D reconstructions of thoracic (T)/lumbar (L) spines from a cohort of 170 Lenke Type-1 patients were analyzed with a non-linear manifold embedding algorithm in order to reduce the high-dimensionality of the data, using statistical properties of neighbouring spine models. We extracted sub-groups of the data from the underlying manifold structure using an unsupervised clustering algorithm to understand the inherent distribution and determine classes of pathologies which appear from the low-dimensional space.

RESULTS

For Lenke Type-1 patients, four clusters were detected from the low-dimensional manifold of 3D models: (1) normal kyphosis (T) with hyper-lordosis (L) and high Cobb angles (37 cases), (2) low kyphosis (T) and normal lordosis (L), with high rotation of plane of maximum curvature (55 cases), (3) hypo-kyphotic (T) and hyper-lordosis (L) (21 cases) and (4) hyper-kyphotic (T) with strong vertebral rotation (57 cases). Results show the manifold representation can potentially be useful for classification of 3D spinal pathologies such as idiopathic scoliosis and serve as a tool for understanding the progression of deformities in longitudinal studies.

CONCLUSIONS

Quantitative evaluation illustrates that the complex space of spine variability can be modeled by a low-dimensional manifold and shows the existence of an additional hyper-kyphotic subgroup from the cohort of 3D spine reconstructions of Lenke Type-1 patients when compared with previous findings on the 3D classification of spinal deformities.

Automated detection of spinal centrelines, vertebral bodies and intervertebral discs in CT and MR images of lumbar spine

We propose a completely automated algorithm for the detection of the spinal centreline and the centres of vertebral bodies and intervertebral discs in images acquired by computed tomography (CT) and magnetic resonance (MR) imaging. The developed methods are based on the analysis of the geometry of spinal structures and the characteristics of CT and MR images and were evaluated on 29 CT and 13 MR images of lumbar spine. The overall mean distance between the obtained and the ground truth spinal centrelines and centres of vertebral bodies and intervertebral discs were 1.8 +/- 1.1 mm and 2.8 +/- 1.9 mm, respectively, and no considerable differences were detected among the results for CT, T(1)-weighted MR and T(2)-weighted MR images. The knowledge of the location of the spinal centreline and the centres of vertebral bodies and intervertebral discs is valuable for the analysis of the spine. The proposed method may therefore be used to initialize the techniques for labelling and segmentation of vertebrae.

Personalized X-ray 3-D reconstruction of the scoliotic spine from hybrid statistical and image-based models

This paper presents a novel 3-D reconstruction method of the scoliotic spine using prior vertebra models with image-based information taken from biplanar X-ray images. We first propose a global modeling approach by exploiting the 3-D scoliotic curve reconstructed from a coronal and sagittal X-ray image in order to generate an approximate statistical model from a 3-D database of scoliotic patients based on a transformation algorithm which incorporates intuitive geometrical properties. The personalized 3-D reconstruction of the spine is then achieved with a novel segmentation method which takes into account the variable appearance of scoliotic vertebrae (rotation, wedging) from standard quality images in order to segment and isolate individual vertebrae on the radiographic planes. More specifically, it uses prior 3-D models regulated from 2-D image level set functionals to identify and match corresponding bone structures on the biplanar X-rays. An iterative optimization procedure integrating similarity measures such as deformable vertebral contours regulated from high-level anatomical primitives, morphological knowledge and epipolar constraints is then applied to globally refine the 3-D anatomical landmarks on each vertebra level of the spine. This method was validated on twenty scoliotic patients by comparing results to a standard manual approach. The qualitative evaluation of the retro-projection of the vertebral contours confirms that the proposed method can achieve better consistency to the X-ray image's natural content. A comparison to synthetic models and real patient data also yields good accuracy on the localization of low-level primitives such as anatomical landmarks identified by an expert on each vertebra. The experiments reported in this paper demonstrate that the proposed method offers a better matching accuracy on a set of landmarks from biplanar views when compared to a manual technique for each evaluated cases, and its precision is comparable to 3-D models generated from magnetic resonance images, thus suitable for routine 3-D clinical assessment of spinal deformities.

A review of methods for quantitative evaluation of spinal curvature

The aim of this paper is to provide a complete overview of the existing methods for quantitative evaluation of spinal curvature from medical images, and to summarize the relevant publications, which may not only assist in the introduction of other researchers to the field, but also be a valuable resource for studying the existing methods or developing new methods and evaluation strategies. Key evaluation issues and future considerations, supported by the results of the overview, are also discussed.

A review of methods for quantitative evaluation of axial vertebral rotation

Quantitative evaluation of axial vertebral rotation is essential for the determination of reference values in normal and pathological conditions and for understanding the mechanisms of the progression of spinal deformities. However, routine quantitative evaluation of axial vertebral rotation is difficult and error-prone due to the limitations of the observer, characteristics of the observed vertebral anatomy and specific imaging properties. The scope of this paper is to review the existing methods for quantitative evaluation of axial vertebral rotation from medical images along with all relevant publications, which may provide a valuable resource for studying the existing methods or developing new methods and evaluation strategies. The reviewed methods are divided into the methods for evaluation of axial vertebral rotation in 2D images and the methods for evaluation of axial vertebral rotation in 3D images. Key evaluation issues and future considerations, supported by the results of the overview, are also discussed.

Automated model-based vertebra detection, identification, and segmentation in CT images

For many orthopaedic, neurological, and oncological applications, an exact segmentation of the vertebral column including an identification of each vertebra is essential. However, although bony structures show high contrast in CT images, the segmentation and labelling of individual vertebrae is challenging. In this paper, we present a comprehensive solution for automatically detecting, identifying, and segmenting vertebrae in CT images. A framework has been designed that takes an arbitrary CT image, e.g., head-neck, thorax, lumbar, or whole spine, as input and provides a segmentation in form of labelled triangulated vertebra surface models. In order to obtain a robust processing chain, profound prior knowledge is applied through the use of various kinds of models covering shape, gradient, and appearance information. The framework has been tested on 64 CT images even including pathologies. In 56 cases, it was successfully applied resulting in a final mean point-to-surface segmentation error of 1.12+/-1.04mm. One key issue is a reliable identification of vertebrae. For a single vertebra, we achieve an identification success of more than 70%. Increasing the number of available vertebrae leads to an increase in the identification rate reaching 100% if 16 or more vertebrae are shown in the image.

Quantitative analysis of spinal curvature in 3D: application to CT images of normal spine

The purpose of this study is to present a framework for quantitative analysis of spinal curvature in 3D. In order to study the properties of such complex 3D structures, we propose two descriptors that capture the characteristics of spinal curvature in 3D. The descriptors are the geometric curvature (GC) and curvature angle (CA), which are independent of the orientation and size of spine anatomy. We demonstrate the two descriptors that characterize the spinal curvature in 3D on 30 computed tomography (CT) images of normal spine and on a scoliotic spine. The descriptors are determined from 3D vertebral body lines, which are obtained by two different methods. The first method is based on the least-squares technique that approximates the manually identified vertebra centroids, while the second method searches for vertebra centroids in an automated optimization scheme, based on computer-assisted image analysis. Polynomial functions of the fourth and fifth degree were used for the description of normal and scoliotic spinal curvature in 3D, respectively. The mean distance to vertebra centroids was 1.1 mm (+/-0.6 mm) for the first and 2.1 mm (+/-1.4 mm) for the second method. The distributions of GC and CA values were obtained along the 30 images of normal spine at each vertebral level and show that maximal thoracic kyphosis (TK), thoracolumbar junction (TJ) and maximal lumbar lordosis (LL) on average occur at T3/T4, T12/L1 and L4/L5, respectively. The main advantage of GC and CA is that the measurements are independent of the orientation and size of the spine, thus allowing objective intra- and inter-subject comparisons. The positions of maximal TK, TJ and maximal LL can be easily identified by observing the GC and CA distributions at different vertebral levels. The obtained courses of the GC and CA for the scoliotic spine were compared to the distributions of GC and CA for the normal spines. The significant difference in values indicates that the descriptors of GC and CA may be used to detect and quantify scoliotic spinal curvatures. The proposed framework may therefore improve the understanding of spine anatomy and aid in the clinical quantitative evaluation of spinal deformities.