Ales Neubert

Automated 3D Segmentation of Vertebral Bodies and Intervertebral Discs from MRI

Recent developments in high resolution MRI scanning of the human spine are providing increasing opportunities for the development of accurate automated approaches for pathoanatomical assessment of intervertebral discs and vertebrae. We are developing a fully automated 3D segmentation approach for MRI scans of the human spine based on statistical shape analysis and template matching of grey level intensity profiles. The algorithm reported in the present study was validated on a dataset of high resolution volumetric scans of lower thoracic and lumbar spine obtained on a 3T scanner using the relatively new 3D SPACE (T2-weighted) pulse sequence, and on a dataset of axial T1-weighted scans of lumbar spine obtained on a 1.5T system. A 3D spine curve is initially extracted and used to position the statistical shape models for final segmentation. Initial validating experiments show promising results on both MRI datasets.

Three-dimensional morphological and signal intensity features for detection of intervertebral disc degeneration from magnetic resonance images

BACKGROUND AND OBJECTIVES Advances in MRI hardware and sequences are continually increasing the amount and complexity of data such as those generated in high-resolution three-dimensional (3D) scanning of the spine. Efficient informatics tools offer considerable opportunities for research and clinically based analyses of magnetic resonance studies. In this work, we present and validate a suite of informatics tools for automated detection of degenerative changes in lumbar intervertebral discs (IVD) from both 3D isotropic and routine two-dimensional (2D) clinical T2-weighted MRI. MATERIALS AND METHODS An automated segmentation approach was used to extract morphological (traditional 2D radiological measures and novel 3D shape descriptors) and signal appearance (extracted from signal intensity histograms) features. The features were validated against manual reference, compared between 2D and 3D MRI scans and used for quantification and classification of IVD degeneration across magnetic resonance datasets containing IVD with early and advanced stages of degeneration. RESULTS AND CONCLUSIONS Combination of the novel 3D-based shape and signal intensity features on 3D (area under receiver operating curve (AUC) 0.984) and 2D (AUC 0.988) magnetic resonance data deliver a significant improvement in automated classification of IVD degeneration, compared to the combination of previously used 2D radiological measurement and signal intensity features (AUC 0.976 and 0.983, respectively). Further work is required regarding the usefulness of 2D and 3D shape data in relation to clinical scores of lower back pain. The results reveal the potential of the proposed informatics system for computer-aided IVD diagnosis from MRI in large-scale research studies and as a possible adjunct for clinical diagnosis.

Validity and reliability of computerized measurement of lumbar intervertebral disc height and volume from magnetic resonance images.

BACKGROUND CONTEXT Magnetic resonance (MR) examinations of morphologic characteristics of intervertebral discs (IVDs) have been used extensively for biomechanical studies and clinical investigations of the lumbar spine. Traditionally, the morphologic measurements have been performed using time- and expertise-intensive manual segmentation techniques not well suited for analyses of large-scale studies.. PURPOSE The purpose of this study is to introduce and validate a semiautomated method for measuring IVD height and mean sagittal area (and volume) from MR images to determine if it can replace the manual assessment and enable analyses of large MR cohorts. STUDY DESIGN/SETTING This study compares semiautomated and manual measurements and assesses their reliability and agreement using data from repeated MR examinations. METHODS Seven healthy asymptomatic males underwent 1.5-T MR examinations of the lumbar spine involving sagittal T2-weighted fast spin-echo images obtained at baseline, pre-exercise, and postexercise conditions. Measures of the mean height and the mean sagittal area of lumbar IVDs (L1-L2 to L4-L5) were compared for two segmentation approaches: a conventional manual method (10-15 minutes to process one IVD) and a specifically developed semiautomated method (requiring only a few mouse clicks to process each subject). RESULTS Both methods showed strong test-retest reproducibility evaluated on baseline and pre-exercise examinations with strong intraclass correlations for the semiautomated and manual methods for mean IVD height (intraclass correlation coefficient [ICC]=0.99, 0.98) and mean IVD area (ICC=0.98, 0.99), respectively. A bias (average deviation) of 0.38 mm (4.1%, 95% confidence interval 0.18-0.59 mm) was observed between the manual and semiautomated methods for the IVD height, whereas there was no statistically significant difference for the mean IVD area (0.1%±3.5%). The semiautomated and manual methods both detected significant exercise-induced changes in IVD height (0.20 and 0.28 mm) and mean IVD area (5.7 and 8.3 mm(2)), respectively. CONCLUSIONS The presented semiautomated method provides an alternative to time- and expertise-intensive manual procedures for analysis of larger, cross-sectional, interventional, and longitudinal MR studies for morphometric analyses of lumbar IVDs.

Evaluation and comparison of 3D intervertebral disc localization and segmentation methods for 3D T2 MR data: A grand challenge.

&NA; The evaluation of changes in Intervertebral Discs (IVDs) with 3D Magnetic Resonance (MR) Imaging (MRI) can be of interest for many clinical applications. This paper presents the evaluation of both IVD localization and IVD segmentation methods submitted to the Automatic 3D MRI IVD Localization and Segmentation challenge, held at the 2015 International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI2015) with an on‐site competition. With the construction of a manually annotated reference data set composed of 25 3D T2‐weighted MR images acquired from two different studies and the establishment of a standard validation framework, quantitative evaluation was performed to compare the results of methods submitted to the challenge. Experimental results show that overall the best localization method achieves a mean localization distance of 0.8 mm and the best segmentation method achieves a mean Dice of 91.8%, a mean average absolute distance of 1.1 mm and a mean Hausdorff distance of 4.3 mm, respectively. The strengths and drawbacks of each method are discussed, which provides insights into the performance of different IVD localization and segmentation methods. HighlightsEstablish a standard framework with 25 manually annotated 3D T2 MRI data for an objective comparison of intervertebral disc (IVD) localization and segmentation methods.Investigate strengths and limitations of a representative selection of the state‐of‐the‐art IVD localization and segmentation methods with a challenge setup.Results achieved by the best algorithms in this study set new frontiers for IVD localization and segmentation from MR data. Graphical abstract Figure. No caption available.

Automated segmentation and analysis of normal and osteoarthritic knee menisci from magnetic resonance images: data from the Osteoarthritis Initiative

OBJECTIVE To validate an automatic scheme for the segmentation and quantitative analysis of the medial meniscus (MM) and lateral meniscus (LM) in magnetic resonance (MR) images of the knee. METHOD We analysed sagittal water-excited double-echo steady-state MR images of the knee from a subset of the Osteoarthritis Initiative (OAI) cohort. The MM and LM were automatically segmented in the MR images based on a deformable model approach. Quantitative parameters including volume, subluxation and tibial-coverage were automatically calculated for comparison (Wilcoxon tests) between knees with variable radiographic osteoarthritis (rOA), medial and lateral joint space narrowing (mJSN, lJSN) and pain. Automatic segmentations and estimated parameters were evaluated for accuracy using manual delineations of the menisci in 88 pathological knee MR examinations at baseline and 12 months time-points. RESULTS The median (95% confidence-interval (CI)) Dice similarity index (DSI) (2 ∗|Auto ∩ Manual|/(|Auto|+|Manual|)∗ 100) between manual and automated segmentations for the MM and LM volumes were 78.3% (75.0-78.7), 83.9% (82.1-83.9) at baseline and 75.3% (72.8-76.9), 83.0% (81.6-83.5) at 12 months. Pearson coefficients between automatic and manual segmentation parameters ranged from r = 0.70 to r = 0.92. MM in rOA/mJSN knees had significantly greater subluxation and smaller tibial-coverage than no-rOA/no-mJSN knees. LM in rOA knees had significantly greater volumes and tibial-coverage than no-rOA knees. CONCLUSION Our automated method successfully segmented the menisci in normal and osteoarthritic knee MR images and detected meaningful morphological differences with respect to rOA and joint space narrowing (JSN). Our approach will facilitate analyses of the menisci in prospective MR cohorts such as the OAI for investigations into pathophysiological changes occurring in early osteoarthritis (OA) development.

Automatic bone segmentation and bone-cartilage interface extraction for the shoulder joint from magnetic resonance images

We present a statistical shape model approach for automated segmentation of the proximal humerus and scapula with subsequent bone-cartilage interface (BCI) extraction from 3D magnetic resonance (MR) images of the shoulder region. Manual and automated bone segmentations from shoulder MR examinations from 25 healthy subjects acquired using steady-state free precession sequences were compared with the Dice similarity coefficient (DSC). The mean DSC scores between the manual and automated segmentations of the humerus and scapula bone volumes surrounding the BCI region were 0.926  ±  0.050 and 0.837  ±  0.059, respectively. The mean DSC values obtained for BCI extraction were 0.806  ±  0.133 for the humerus and 0.795  ±  0.117 for the scapula. The current model-based approach successfully provided automated bone segmentation and BCI extraction from MR images of the shoulder. In future work, this framework appears to provide a promising avenue for automated segmentation and quantitative analysis of cartilage in the glenohumeral joint.

Automatic segmentation of the glenohumeral cartilages from magnetic resonance images

PURPOSE Magnetic resonance (MR) imaging plays a key role in investigating early degenerative disorders and traumatic injuries of the glenohumeral cartilages. Subtle morphometric and biochemical changes of potential relevance to clinical diagnosis, treatment planning, and evaluation can be assessed from measurements derived from in vivo MR segmentation of the cartilages. However, segmentation of the glenohumeral cartilages, using approaches spanning manual to automated methods, is technically challenging, due to their thin, curved structure and overlapping intensities of surrounding tissues. Automatic segmentation of the glenohumeral cartilages from MR imaging is not at the same level compared to the weight-bearing knee and hip joint cartilages despite the potential applications with respect to clinical investigation of shoulder disorders. In this work, the authors present a fully automated segmentation method for the glenohumeral cartilages using MR images of healthy shoulders. METHODS The method involves automated segmentation of the humerus and scapula bones using 3D active shape models, the extraction of the expected bone-cartilage interface, and cartilage segmentation using a graph-based method. The cartilage segmentation uses localization, patient specific tissue estimation, and a model of the cartilage thickness variation. The accuracy of this method was experimentally validated using a leave-one-out scheme on a database of MR images acquired from 44 asymptomatic subjects with a true fast imaging with steady state precession sequence on a 3 T scanner (Siemens Trio) using a dedicated shoulder coil. The automated results were compared to manual segmentations from two experts (an experienced radiographer and an experienced musculoskeletal anatomist) using the Dice similarity coefficient (DSC) and mean absolute surface distance (MASD) metrics. RESULTS Accurate and precise bone segmentations were achieved with mean DSC of 0.98 and 0.93 for the humeral head and glenoid fossa, respectively. Mean DSC scores of 0.74 and 0.72 were obtained for the humeral and glenoid cartilage volumes, respectively. The manual interobserver reliability evaluated by DSC was 0.80 ± 0.03 and 0.76 ± 0.04 for the two cartilages, implying that the automated results were within an acceptable 10% difference. The MASD between the automatic and the corresponding manual cartilage segmentations was less than 0.4 mm (previous studies reported mean cartilage thickness of 1.3 mm). CONCLUSIONS This work shows the feasibility of volumetric segmentation and separation of the glenohumeral cartilages from MR images. To their knowledge, this is the first fully automated algorithm for volumetric segmentation of the individual glenohumeral cartilages from MR images. The approach was validated against manual segmentations from experienced analysts. In future work, the approach will be validated on imaging datasets acquired with various MR contrasts in patients.

Constrained reverse diffusion for thick slice interpolation of 3D volumetric MRI images.

Due to physical limitations inherent in magnetic resonance imaging scanners, three dimensional volumetric scans are often acquired with anisotropic voxel resolution. We investigate several interpolation approaches to reduce the anisotropy and present a novel approach - constrained reverse diffusion for thick slice interpolation. This technique was compared to common methods: linear and cubic B-Spline interpolation and a technique based on non-rigid registration of neighboring slices. The methods were evaluated on artificial MR phantoms and real MR scans of human brain. The constrained reverse diffusion approach delivered promising results and provides an alternative for thick slice interpolation, especially for higher anisotropy factors.

Local contrast-enhanced MR images via high dynamic range processing

To develop a local contrast‐enhancing and feature‐preserving high dynamic range (HDR) image processing algorithm for multichannel and multisequence MR images of multiple body regions and tissues, and to evaluate its performance for structure visualization, bias field (correction) mitigation, and automated tissue segmentation.

Comparison of 3D bone models of the knee joint derived from CT and 3T MR imaging

PURPOSE To examine whether magnetic resonance (MR) imaging can offer a viable alternative to computed tomography (CT) based 3D bone modeling. METHODS CT and MR (SPACE, TrueFISP, VIBE) images were acquired from the left knee joint of a fresh-frozen cadaver. The distal femur, proximal tibia, proximal fibula and patella were manually segmented from the MR and CT examinations. The MR bone models obtained from manual segmentations of all three sequences were compared to CT models using a similarity measure based on absolute mesh differences. RESULTS The average absolute distance between the CT and the various MR-based bone models were all below 1mm across all bones. The VIBE sequence provided the best agreement with the CT model, followed by the SPACE, then the TrueFISP data. The most notable difference was for the proximal tibia (VIBE 0.45mm, SPACE 0.82mm, TrueFISP 0.83mm). CONCLUSIONS The study indicates that 3D MR bone models may offer a feasible alternative to traditional CT-based modeling. A single radiological examination using the MR imaging would allow simultaneous assessment of both bones and soft-tissues, providing anatomically comprehensive joint models for clinical evaluation, without the ionizing radiation of CT imaging.