Abtin Rasoulian

Lumbar Spine Segmentation Using a Statistical Multi-Vertebrae Anatomical Shape+Pose Model

Segmentation of the spinal column from computed tomography (CT) images is a preprocessing step for a range of image-guided interventions. One intervention that would benefit from accurate segmentation is spinal needle injection. Previous spinal segmentation techniques have primarily focused on identification and separate segmentation of each vertebra. Recently, statistical multi-object shape models have been introduced to extract common statistical characteristics between several anatomies. These models can be used for segmentation purposes because they are robust, accurate, and computationally tractable. In this paper, we develop a statistical multi-vertebrae shape+pose model and propose a novel registration-based technique to segment the CT images of spine. The multi-vertebrae statistical model captures the variations in shape and pose simultaneously, which reduces the number of registration parameters. We validate our technique in terms of accuracy and robustness of multi-vertebrae segmentation of CT images acquired from lumbar vertebrae of 32 subjects. The mean error of the proposed technique is below 2 mm, which is sufficient for many spinal needle injection procedures, such as facet joint injections.

Group-Wise Registration of Point Sets for Statistical Shape Models

This paper presents a novel, fast, group-wise registration technique based on establishing soft correspondences between groups of point sets. The registration approach is used to create a statistical shape model, capable of learning the shape variations within a training set. The shape model consists of a mean shape and its transformations to all training shapes. We decouple the procedure into two steps: updating the mean shape and registering it to the training shapes. The algorithm alternates between these two steps until convergence. Following the generation of the statistical shape model, we use the soft correspondence approach to register the model to a new observation. We perform extensive experiments on two data sets: lumbar spine and hippocampi. We compare our algorithm to available state-of-the-art group-wise registration algorithms including feature-based and volume-based approaches. We demonstrate improved generalization, specificity and compactness compared to these algorithms.

A multi-center milestone study of clinical vertebral CT segmentation☆

A multiple center milestone study of clinical vertebra segmentation is presented in this paper. Vertebra segmentation is a fundamental step for spinal image analysis and intervention. The first half of the study was conducted in the spine segmentation challenge in 2014 International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI) Workshop on Computational Spine Imaging (CSI 2014). The objective was to evaluate the performance of several state-of-the-art vertebra segmentation algorithms on computed tomography (CT) scans using ten training and five testing dataset, all healthy cases; the second half of the study was conducted after the challenge, where additional 5 abnormal cases are used for testing to evaluate the performance under abnormal cases. Dice coefficients and absolute surface distances were used as evaluation metrics. Segmentation of each vertebra as a single geometric unit, as well as separate segmentation of vertebra substructures, was evaluated. Five teams participated in the comparative study. The top performers in the study achieved Dice coefficient of 0.93 in the upper thoracic, 0.95 in the lower thoracic and 0.96 in the lumbar spine for healthy cases, and 0.88 in the upper thoracic, 0.89 in the lower thoracic and 0.92 in the lumbar spine for osteoporotic and fractured cases. The strengths and weaknesses of each method as well as future suggestion for improvement are discussed. This is the first multi-center comparative study for vertebra segmentation methods, which will provide an up-to-date performance milestone for the fast growing spinal image analysis and intervention.

Feature-based multibody rigid registration of CT and ultrasound images of lumbar spine.

PURPOSE Fusion of intraprocedure ultrasound and preprocedure CT data is proposed for guidance in percutaneous spinal injections, a common procedure for pain management. CT scan of the lumbar spine is usually collected in a supine position, whereas spinal injections are performed in prone or sitting positions. This leads to a difference in the spine curvature between CT and ultrasound images; as such, a single-body rigid registration approach cannot be used for the whole lumbar vertebrae. METHODS To compensate for the difference in the spinal curvature between the two imaging modalities, a multibody rigid registration algorithm is presented. The approach utilizes a point-based registration technique based on the unscented Kalman filter (UKF), taking as input segmented vertebrae surfaces in both CT and ultrasound data. Ultrasound images are automatically segmented using a dynamic programming method, while the CT images are semiautomatically segmented using thresholding. The registration approach is designed to simultaneously align individual groups of points segmented from each vertebra in the two imaging modalities. A biomechanical model is developed to constrain the vertebrae transformation parameters during the registration and to ensure convergence. RESULTS The proposed methodology is evaluated on human spine phantoms and a sheep cadaver. Registrations on phantom data have a mean target registration error (TRE) of 1.99 mm in the region of interest and converged in 87% of cases. Registrations on sheep cadaver have a mean target registration error of 2.2 mm and converged in 82% of cases. CONCLUSIONS The proposed technique can robustly and simultaneously register several vertebrae extracted from CT images to the ultrasound volumes. The registration error below 2.2 mm is sufficient for most spinal injections.

Local Phase Tensor Features for 3-D Ultrasound to Statistical Shape+Pose Spine Model Registration

Most conventional spine interventions are performed under X-ray fluoroscopy guidance. In recent years, there has been a growing interest to develop nonionizing imaging alternatives to guide these procedures. Ultrasound guidance has emerged as a leading alternative. However, a challenging problem is automatic identification of the spinal anatomy in ultrasound data. In this paper, we propose a local phase-based bone feature enhancement technique that can robustly identify the spine surface in ultrasound images. The local phase information is obtained using a gradient energy tensor filter. This information is used to construct local phase tensors in ultrasound images, which highlight the spine surface. We show that our proposed approach results in a more distinct enhancement of the bone surfaces compared to recently proposed techniques based on monogenic scale-space filters and logarithmic Gabor filters. We also demonstrate that registration accuracy of a statistical shape+pose model of the spine to 3-D ultrasound images can be significantly improved, using the proposed method, compared to those obtained using monogenic scale-space filters and logarithmic Gabor filters.

Augmentation of paramedian 3D ultrasound images of the spine

The blind placement of an epidural needle is among the most difficult regional anesthetic techniques. The challenge is to insert the needle in the mid-sagittal plane and to avoid overshooting the needle into the spinal cord. Prepuncture 2D ultrasound scanning has been introduced as a reliable tool to localize the target and facilitate epidural needle placement. Ideally, real-time ultrasound should be used during needle insertion. However, several issues inhibit the use of standard 2D ultrasound, including the obstruction of the puncture site by the ultrasound probe, low visibility of the target in ultrasound images, and increased pain due to longer needle trajectory. An alternative is to use 3D ultrasound imaging, where the needle and target could be visible within the same reslice of a 3D volume; however, novice ultrasound users (i.e., many anesthesiologists) still have difficulty interpreting ultrasound images of the spine and identifying the target epidural space. In this paper, we propose to augment 3D ultrasound images by registering a multi-vertebrae statistical shape+pose model. We use such augmentation for enhanced interpretation of the ultrasound and identification of the mid-sagittal plane for the needle insertion. Validation is performed on synthetic data derived from the CT images, and 64 in vivo ultrasound volumes.

Statistical Shape Model to 3D Ultrasound Registration for Spine Interventions Using Enhanced Local Phase Features

Accurate registration of ultrasound images to statistical shape models is a challenging problem in percutaneous spine injection procedures due to the typical imaging artifacts inherent to ultrasound. In this paper we propose a robust and accurate registration method that matches local phase bone features extracted from ultrasound images to a statistical shape model. The local phase information for enhancing the bone surfaces is obtained using a gradient energy tensor filter, which combines advantages of the monogenic scale-space and Gaussian scale-space filters, resulting in an improved simultaneous estimation of phase and orientation information. A novel statistical shape model was built by separating the pose statistics from the shape statistics. This model is then registered to the local phase bone surfaces using an iterative expectation maximization registration technique. Validation on 96 in vivo clinical scans obtained from eight patients resulted in a root mean square registration error of 2 mm (SD: 0.4 mm), which is below the clinically acceptable threshold of 3.5 mm. The improvement achieved in registration accuracy using the new features was also significant (p < 0.05) compared to state of the art local phase image processing methods.

Biomechanically Constrained Surface Registration: Application to MR-TRUS Fusion for Prostate Interventions

In surface-based registration for image-guided interventions, the presence of missing data can be a significant issue. This often arises with real-time imaging modalities such as ultrasound, where poor contrast can make tissue boundaries difficult to distinguish from surrounding tissue. Missing data poses two challenges: ambiguity in establishing correspondences; and extrapolation of the deformation field to those missing regions. To address these, we present a novel non-rigid registration method. For establishing correspondences, we use a probabilistic framework based on a Gaussian mixture model (GMM) that treats one surface as a potentially partial observation. To extrapolate and constrain the deformation field, we incorporate biomechanical prior knowledge in the form of a finite element model (FEM). We validate the algorithm, referred to as GMM-FEM, in the context of prostate interventions. Our method leads to a significant reduction in target registration error (TRE) compared to similar state-of-the-art registration algorithms in the case of missing data up to 30%, with a mean TRE of 2.6 mm. The method also performs well when full segmentations are available, leading to TREs that are comparable to or better than other surface-based techniques. We also analyze robustness of our approach, showing that GMM-FEM is a practical and reliable solution for surface-based registration.

Temporal-lobe morphology differs between healthy adolescents and those with early-onset of depression

Major depressive disorder (MDD) has previously been linked to structural changes in several brain regions, particularly in the medial temporal lobes (Bellani, Baiano, Brambilla, 2010; Bellani, Baiano, Brambilla, 2011). This has been determined using voxel-based morphometry, segmentation algorithms, and analysis of shape deformations (Bell-McGinty et al., 2002; Bergouignan et al., 2009; Posener et al., 2003; Vasic et al., 2008; Zhao et al., 2008): these are methods in which information related to the shape and the pose (the size, and anatomical position and orientation) of structures is lost. Here, we incorporate information about shape and pose to measure structural deformation in adolescents and young adults with and without depression (as measured using the Beck Depression Inventory and Diagnostic and Statistical Manual of Mental Disorders criteria). As a hypothesis-generating study, a significance level of p < 0.05, uncorrected for multiple comparisons, was used, so that subtle morphological differences in brain structures between adolescent depressed individuals and control participants could be identified. We focus on changes in cortical and subcortical temporal structures, and use a multi-object statistical pose and shape model to analyze imaging data from 16 females (aged 16–21) and 3 males (aged 18) with early-onset MDD, and 25 female and 1 male normal control participants, drawn from the same age range. The hippocampus, parahippocampal gyrus, putamen, and superior, inferior and middle temporal gyri in both hemispheres of the brain were automatically segmented using the LONI Probabilistic Brain Atlas (Shattuck et al., 2008) in MNI space. Points on the surface of each structure in the atlas were extracted and warped to each participants structural MRI. These surface points were analyzed to extract the pose and shape features. Pose differences were detected between the two groups, particularly in the left and right putamina, right hippocampus, and left and right inferior temporal gyri. Shape differences were detected between the two groups, particularly in the left hippocampus and in the left and right parahippocampal gyri. Furthermore, pose measures were significantly correlated with BDI score across the whole (clinical and control) sample. Since the clinical participants were experiencing their very first episodes of MDD, morphological alteration in the medial temporal lobe appears to be an early sign of MDD, and is unlikely to result from treatment with antidepressants. Pose and shape measures of morphology, which are not usually analyzed in neuromorphometric studies, appear to be sensitive to depressive symptomatology.

Deep learning for automatic localization, identification, and segmentation of vertebral bodies in volumetric MR images

This paper proposes an automatic method for vertebra localization, labeling, and segmentation in multi-slice Magnetic Resonance (MR) images. Prior work in this area on MR images mostly requires user interaction while our method is fully automatic. Cubic intensity-based features are extracted from image voxels. A deep learning approach is used for simultaneous localization and identification of vertebrae. The localized points are refined by local thresholding in the region of the detected vertebral column. Thereafter, a statistical multi-vertebrae model is initialized on the localized vertebrae. An iterative Expectation Maximization technique is used to register the vertebral body of the model to the image edges and obtain a segmentation of the lumbar vertebral bodies. The method is evaluated by applying to nine volumetric MR images of the spine. The results demonstrate 100% vertebra identification and a mean surface error of below 2.8 mm for 3D segmentation. Computation time is less than three minutes per high-resolution volumetric image.