Craig Engstrom

MR Morphometry of Posterior Tibialis Muscle in Adult Acquired Flat Foot

We conducted magnetic resonance imaging of the posterior tibial (PT) and flexor digitorum longus (FDL) muscle bellies in 12 patients undergoing surgical treatment for unilateral posterior tibial tendon (PTT) dysfunction. All patients had atrophy of the PT muscle compared to the normal leg (mean 10.7%, p=0.008). In those patients with a complete rupture of PTT there was replacement of the PT muscle by fatty infiltration. Conversely, the FDL muscle showed a compensatory hypertrophy (mean 17.2%, p<0.002). We support the use of FDL as an appropriate tendon for augmentation of PTT in stage II disease. This study also demonstrates that in the presence of a complete rupture, excision of the PTT is a reasonable surgical procedure and pure tenodesis will be more likely to fail because the PT muscle belly undergoes fatty infiltration. In patients with a diseased but intact PTT there was no fatty infiltration and the muscle volume was at least 83% of the normal side in all cases. We therefore suggest that in the presence of an intact PTT the PT muscle belly may provide some useful function if used to augment the FDL transfer when the diseased tendon is excised.

Focused shape models for hip joint segmentation in 3D magnetic resonance images

Deformable models incorporating shape priors have proved to be a successful approach in segmenting anatomical regions and specific structures in medical images. This paper introduces weighted shape priors for deformable models in the context of 3D magnetic resonance (MR) image segmentation of the bony elements of the human hip joint. The fully automated approach allows the focusing of the shape model energy to a priori selected anatomical structures or regions of clinical interest by preferentially ordering the shape representation (or eigen-modes) within this type of model to the highly weighted areas. This focused shape model improves accuracy of the shape constraints in those regions compared to standard approaches. The proposed method achieved femoral head and acetabular bone segmentation mean absolute surface distance errors of 0.55±0.18mm and 0.75±0.20mm respectively in 35 3D unilateral MR datasets from 25 subjects acquired at 3T with different limited field of views for individual bony components of the hip joint.

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.

Automatic hip cartilage segmentation from 3D MR images using arc-weighted graph searching

Accurate segmentation of hip joint cartilage from magnetic resonance (MR) images offers opportunities for quantitative investigations of pathoanatomical conditions such as osteoarthritis. In this paper, we present a fully automatic scheme for the segmentation of the individual femoral and acetabular cartilage plates in the human hip joint from high-resolution 3D MR images. The developed scheme uses an improved optimal multi-object multi-surface graph search framework with an arc-weighted graph representation that incorporates prior morphological knowledge as a basis for segmentation of the individual femoral and acetabular cartilage plates despite weak or incomplete boundary interfaces. This automated scheme was validated against manual segmentations from 3D true fast imaging with steady-state precession (TrueFISP) MR examinations of the right hip joints in 52 asymptomatic volunteers. Compared with expert manual segmentations of the combined, femoral and acetabular cartilage volumes, the automatic scheme obtained mean (± standard deviation) Dices similarity coefficients of 0.81 (± 0.03), 0.79 (± 0.03) and 0.72 (± 0.05). The corresponding mean absolute volume difference errors were 8.44% (± 6.36), 9.44% (± 7.19) and 9.05% (± 8.02). The mean absolute differences between manual and automated measures of cartilage thickness for femoral and acetabular cartilage plates were 0.13 mm (± 0.12) and 0.11 mm (± 0.11), respectively.

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.

Segmentation of the quadratus lumborum muscle using statistical shape modeling

To compare automated segmentation of the quadratus lumborum (QL) based on statistical shape modeling (SSM) with conventional manual processing of magnetic resonance (MR) images for segmentation of this paraspinal muscle.

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.

Automated MR Hip Bone Segmentation

The accurate segmentation of the bone and articular cartilages from magnetic resonance (MR) images of the hip is important for clinical studies and drug trials into conditions like osteoarthritis. In current studies, segmentations are obtained using time-consuming manual or semi-automatic algorithms which have high inter- and intra-observer variabilities. This paper presents an important step towards obtaining automatic and accurate segmentations of the hip cartilages, namely an approach to automatically segment the bones. The segmentation is performed using three-dimensional active shape models, which are initialized using an affine registration to an atlas. The accuracy and robustness of the approach was experimentally validated using an MR database of we VIBE, we DESS and MEDIC MR images. The (left, right) femoral and ace tabular bone segmentation had a median Dice similarity coefficient of (0.921, 0.926) and (0.830, 0.813).