Quantitating Age-Related BMD Textural Variation from DXA Region-Free-Analysis: A Study of Hip Fracture Prediction in Three Cohorts.

The risk of osteoporotic fracture is inversely related to bone mineral density (BMD), but how spatial BMD pattern influences fracture risk remains incompletely understood. This study used a pixel-level spatiotemporal atlas of proximal femoral BMD in 13,338 white European women (age 20-97 years) to quantitate age-related texture variation in BMD maps and generate a "reference" map of bone aging. We introduce a new index, called Densitometric Bone Age (DBA), as the age at which an individual site-specific BMD map (the proximal femur is studied here) best matches the median aging trajectory at that site in terms of the root mean squared error (RMSE). The ability of DBA to predict incident hip fracture and hip fracture pattern over 5 years following baseline BMD was compared against conventional region-based BMD analysis in a subset of 11,899 women (age 45-97 years), for which follow-up fracture records exist. There were 208 subsequent incident hip fractures in the study populations (138 femoral necks [FNs], 52 trochanteric [TR], 18 sites unspecified). DBA had modestly better performance compared to the conventional FN-BMD, TR-BMD, and total hip (TOT)-BMD in identifying hip fractures measured as the area under the curve (AUC) using receiver operating characteristics (ROC) curve analysis by 2% (95% confidence interval [CI], -0.5% to 3.5%), 3% (95% CI, 1.0% to 4.0%), and 1% (95% CI, 0.4% to 1.6%), respectively. Compared to FN-BMD T-score, DBA improved the ROC-AUC for predicting TR fractures by ~5% (95% CI, 1.1% to 9.8%) with similar performance in identifying FN fractures. Compared to TR-BMD T-score, DBA improved the ROC-AUC for the prediction of FN fractures by ~3% (95% CI, 1.1% to 4.9%), with similar performance in identifying TR fractures. Our findings suggest that DBA may provide a spatially sensitive measure of proximal femoral fragility that is not captured by FN-BMD or TR-BMD alone. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Origami: Single-cell 3D shape dynamics oriented along the apico-basal axis of folding epithelia from fluorescence microscopy data.

A common feature of morphogenesis is the formation of three-dimensional structures from the folding of two-dimensional epithelial sheets, aided by cell shape changes at the cellular-level. Changes in cell shape must be studied in the context of cell-polarised biomechanical processes within the epithelial sheet. In epithelia with highly curved surfaces, finding single-cell alignment along a biological axis can be difficult to automate in silico. We present 'Origami', a MATLAB-based image analysis pipeline to compute direction-variant cell shape features along the epithelial apico-basal axis. Our automated method accurately computed direction vectors denoting the apico-basal axis in regions with opposing curvature in synthetic epithelia and fluorescence images of zebrafish embryos. As proof of concept, we identified different cell shape signatures in the developing zebrafish inner ear, where the epithelium deforms in opposite orientations to form different structures. Origami is designed to be user-friendly and is generally applicable to fluorescence images of curved epithelia.

A Spatio-Temporal Ageing Atlas of the Proximal Femur.

Osteoporosis is an age-associated disease characterised by low bone mineral density (BMD) and micro-architectural deterioration leading to enhanced fracture risk. Conventional dual-energy X-ray absorptiometry (DXA) analysis has facilitated our understanding of BMD reduction in specific regions of interest (ROIs) within the femur, but cannot resolve spatial BMD patterns nor reflect age-related changes in bone microarchitecture due to its inherent averaging of pixel BMD values into large ROIs. To address these limitations and develop a comprehensive model of involutional bone loss, this paper presents a fully automatic pipeline to build a spatio-temporal atlas of ageing bone in the proximal femur. A new technique, termed DXA region free analysis (DXA RFA), is proposed to eliminate morphological variation between DXA scans by warping each image into a reference template. To construct the atlas, we use unprocessed DXA data from Caucasian women aged 20-97 years participating in three cohort studies in Western Europe ( ,000). A novel calibration procedure, termed quantile matching regression, is proposed to integrate data from different DXA manufacturers. Pixel-wise BMD evolution with ageing was modelled using smooth quantile curves. This technique enables characterisation of spatially-complex BMD change patterns with ageing, visualised using heat-maps. Furthermore, quantile curves plotted at different pixel coordinates showed consistently different rates of bone loss at different regions within the femoral neck. Given the close relationship between spatio-temporal bone loss and osteoporotic fracture, improved understanding of the bone ageing process could lead to enhanced prognostic, preventive and therapeutic strategies for the disease.

Diffusion MRI for Assessment of Bone Quality; A Review of Findings in Healthy Aging and Osteoporosis.

Diffusion MRI (dMRI) is a growing imaging technique with the potential to provide biomarkers of tissue variation, such as cellular density, tissue anisotropy, and microvascular perfusion. However, the role of dMRI in characterizing different aspects of bone quality, especially in aging and osteoporosis, has not yet been fully established, particularly in clinical applications. The reason lies in the complications accompanied with implementation of dMRI in assessment of human bone structure, in terms of acquisition and quantification. Bone is a composite tissue comprising different elements, each contributing to the overall quality and functional competence of bone. As diffusion is a critical biophysical process in biological tissues, early changes of tissue microstructure and function can affect diffusive properties of the tissue. While there are multiple MRI methods to detect variations of individual properties of bone quality due to aging and osteoporosis, dMRI has potential to serve as a superior method for characterizing different aspects of bone quality within the same framework but with higher sensitivity to early alterations. This is mainly because several properties of the tissue including directionality and anisotropy of trabecular bone and cell density can be collected using only dMRI. In this review article, we first describe components of human bone that can be potentially detected by their diffusivity properties and contribute to variations in bone quality during aging and osteoporosis. Then we discuss considerations and challenges of dMRI in bone imaging, current status, and suggestions for development of dMRI in research studies and clinics to segregate different contributing components of bone quality in an integrated acquisition. Level of Evidence: 5 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2020;51:975-992.

Histological data of axons, astrocytes, and myelin in deep subcortical white matter populations.

This immunohistochemistry dataset contains the main structures in deep subcortical white matter (axons, astrocytes, and myelinated axons) in a representative cohort of an ageing population. A set of samples from 90 subjects of the Cognitive Function and Ageing Study (CFAS) were analysed, stratified into three groups of 30 subjects each, in relation to the presence of age-associated deep subcortical lesions. High-resolution microscopy data enables the extraction of valuable information, such as volume fractions, for the construction and validation of diffusion MRI (dMRI) models. The dataset provided here was used in Coelho et al. [1].

Quantitative histomorphometry of capillary microstructure in deep white matter.

White matter lesions represent a major risk factor for dementia in elderly people. Magnetic Resonance Imaging (MRI) studies have demonstrated cerebral blood flow reduction in age-related white matter lesions, indicating that vascular alterations are involved in developing white matter lesions. Hypoperfusion and changes in capillary morphology are generally linked to dementia. However, a quantitative study describing these microvascular alterations in white matter lesions is missing in the literature; most previous microvascular studies being on the cortex. The aim of this work is to identify and quantify capillary microstructural changes involved in the appearance of deep subcortical lesions (DSCL). We characterize the distribution of capillary diameter, thickness, and density in the deep white matter in a population of 75 elderly subjects, stratified into three equal groups according to DSCL: Control (subject without DSCL), Lesion (sample presenting DSCL), and Normal Appearing White Matter (NAWM, the subject presented DSCL but not at the sampled tissue location). Tissue samples were selected from the Cognitive Function and Aging Study (CFAS), a cohort representative of an aging population, from which immunohistochemically-labeled histological images were acquired. To accurately estimate capillary diameters and thicknesses from the 2D histological images, we also introduce a novel semi-automatic method robust to non-perpendicular incidence angle of capillaries into the imaging plane, and to non-circular deformations of capillary cross sections. Subjects with DSCL presented a significant increase in capillary wall thickness, a decrease in the diameter intra-subject variability (but not in the mean), and a decrease in capillary density. No significant difference was observed between controls and NAWM. Both capillary wall thickening and reduction in capillary density contribute to the reduction of cerebral blood flow previously reported for white matter lesions. The obtained distributions provide reliable statistics of capillary microstructure useful to inform the modeling of human cerebral blood flow, for instance to define microcirculation models for their estimation from MRI or to perform realistic cerebral blood flow simulations.

Population-based Bayesian regularization for microstructural diffusion MRI with NODDIDA.

PURPOSE: Information on the brain microstructure can be probed by Diffusion Magnetic Resonance Imaging (dMRI). Neurite Orientation Dispersion and Density Imaging with Diffusivities Assessment (NODDIDA) is one of the simplest microstructural model proposed. However, the estimation of the NODDIDA parameters from clinically plausible dMRI acquisition is ill-posed, and different parameter sets can describe the same measurements equally well. A few approaches to resolve this problem focused on developing better optimization strategies for this non-convex optimization. However, this fundamentally does not resolve ill-posedness. This article introduces a Bayesian estimation framework, which is regularized through knowledge from an extensive dMRI measurement set on a population of healthy adults (henceforth population-based prior). METHODS: We reformulate the problem as a Bayesian maximum a posteriori estimation, which includes as a special case previous approach using non-informative uniform priors. A population-based prior is estimated from 35 subjects of the MGH Adult Diffusion data (Human Connectome Project), acquired with an extensive acquisition protocol including high b-values. The accuracy and robustness of different approaches with and without the population-based prior is tested on subsets of the MGH dataset, and an independent dataset from a clinically comparable scanner, with only clinically plausible dMRI measurements. RESULTS: The population-based prior produced substantially more accurate and robust parameter estimates, compared to the conventional uniform priors, for clinically feasible protocols, without introducing any evident bias. CONCLUSIONS: The use of the proposed Bayesian population-based prior can lead to clinically feasible and robust estimation of NODDIDA parameters without changing the acquisition protocol.

Resolving degeneracy in diffusion MRI biophysical model parameter estimation using double diffusion encoding.

PURPOSE: Biophysical tissue models are increasingly used in the interpretation of diffusion MRI (dMRI) data, with the potential to provide specific biomarkers of brain microstructural changes. However, it has been shown recently that, in the general Standard Model, parameter estimation from dMRI data is ill-conditioned even when very high b-values are applied. We analyze this issue for the Neurite Orientation Dispersion and Density Imaging with Diffusivity Assessment (NODDIDA) model and demonstrate that its extension from single diffusion encoding (SDE) to double diffusion encoding (DDE) resolves the ill-posedness for intermediate diffusion weightings, producing an increase in accuracy and precision of the parameter estimation. METHODS: We analyze theoretically the cumulant expansion up to fourth order in b of SDE and DDE signals. Additionally, we perform in silico experiments to compare SDE and DDE capabilities under similar noise conditions. RESULTS: We prove analytically that DDE provides invariant information non-accessible from SDE, which makes the NODDIDA parameter estimation injective. The in silico experiments show that DDE reduces the bias and mean square error of the estimation along the whole feasible region of 5D model parameter space. CONCLUSIONS: DDE adds additional information for estimating the model parameters, unexplored by SDE. We show, as an example, that this is sufficient to solve the previously reported degeneracies in the NODDIDA model parameter estimation.

Iba-1-/CD68+ microglia are a prominent feature of age-associated deep subcortical white matter lesions.

Deep subcortical lesions (DSCL) of the brain, are present in ~60% of the ageing population, and are linked to cognitive decline and depression. DSCL are associated with demyelination, blood brain barrier (BBB) dysfunction, and microgliosis. Microglia are the main immune cell of the brain. Under physiological conditions microglia have a ramified morphology, and react to pathology with a change to a more rounded morphology as well as showing protein expression alterations. This study builds on previous characterisations of DSCL and radiologically 'normal-appearing' white matter (NAWM) by performing a detailed characterisation of a range of microglial markers in addition to markers of vascular integrity. The Cognitive Function and Ageing Study (CFAS) provided control white matter (WM), NAWM and DSCL human post mortem tissue for immunohistochemistry using microglial markers (Iba-1, CD68 and MHCII), a vascular basement membrane marker (collagen IV) and markers of BBB integrity (fibrinogen and aquaporin 4). The immunoreactive profile of CD68 increased in a stepwise manner from control WM to NAWM to DSCL. This correlated with a shift from small, ramified cells, to larger, more rounded microglia. While there was greater Iba-1 immunoreactivity in NAWM compared to controls, in DSCL, Iba-1 levels were reduced to control levels. A prominent feature of these DSCL was a population of Iba-1-/CD68+ microglia. There were increases in collagen IV, but no change in BBB integrity. Overall the study shows significant differences in the immunoreactive profile of microglial markers. Whether this is a cause or effect of lesion development remains to be elucidated. Identifying microglia subpopulations based on their morphology and molecular markers may ultimately help decipher their function and role in neurodegeneration. Furthermore, this study demonstrates that Iba-1 is not a pan-microglial marker, and that a combination of several microglial markers is required to fully characterise the microglial phenotype.

Local volume fraction distributions of axons, astrocytes, and myelin in deep subcortical white matter.

This study aims to statistically describe histologically stained white matter brain sections to subsequently inform and validate diffusion MRI techniques. For the first time, we characterise volume fraction distributions of three of the main structures in deep subcortical white matter (axons, astrocytes, and myelinated axons) in a representative cohort of an ageing population for which well-characterized neuropathology data is available. We analysed a set of samples from 90 subjects of the Cognitive Function and Ageing Study (CFAS), stratified into three groups of 30 subjects each, in relation to the presence of age-associated deep subcortical lesions. This provides volume fraction distributions in different scenarios relevant to brain diffusion MRI in dementia. We also assess statistically significant differences found between these groups. In agreement with previous literature, our results indicate that white matter lesions are related with a decrease in the myelinated axons fraction and an increase in astrocytic fraction, while no statistically significant changes occur in axonal mean fraction. In addition, we introduced a framework to quantify volume fraction distributions from 2D immunohistochemistry images, which is validated against in silico simulations. Since a trade-off between precision and resolution emerged, we also performed an assessment of the optimal scale for computing such distributions.

Simultaneous magnetic resonance diffusion and pseudo-diffusion tensor imaging.

PURPOSE: An emerging topic in diffusion magnetic resonance is imaging blood microcirculation alongside water diffusion using the intravoxel incoherent motion (IVIM) model. Recently, a combined IVIM diffusion tensor imaging (IVIM-DTI) model was proposed, which accounts for both anisotropic pseudo-diffusion due to blood microcirculation and anisotropic diffusion due to tissue microstructures. In this article, we propose a robust IVIM-DTI approach for simultaneous diffusion and pseudo-diffusion tensor imaging. METHODS: Conventional IVIM estimation methods can be broadly divided into two-step (diffusion and pseudo-diffusion estimated separately) and one-step (diffusion and pseudo-diffusion estimated simultaneously) methods. Here, both methods were applied on the IVIM-DTI model. An improved one-step method based on damped Gauss-Newton algorithm and a Gaussian prior for the model parameters was also introduced. The sensitivities of these methods to different parameter initializations were tested with realistic in silico simulations and experimental in vivo data. RESULTS: The one-step damped Gauss-Newton method with a Gaussian prior was less sensitive to noise and the choice of initial parameters and delivered more accurate estimates of IVIM-DTI parameters compared to the other methods. CONCLUSION: One-step estimation using damped Gauss-Newton and a Gaussian prior is a robust method for simultaneous diffusion and pseudo-diffusion tensor imaging using IVIM-DTI model. Magn Reson Med 79:2367-2378, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Robustness of common hemodynamic indicators with respect to numerical resolution in 38 middle cerebral artery aneurysms.

BACKGROUND: Using computational fluid dynamics (CFD) to compute the hemodynamics in cerebral aneurysms has received much attention in the last decade. The usability of these methods depends on the quality of the computations, highlighted in recent discussions. The purpose of this study is to investigate the convergence of common hemodynamic indicators with respect to numerical resolution. METHODS: 38 middle cerebral artery bifurcation aneurysms were studied at two different resolutions (one comparable to most studies, and one finer). Relevant hemodynamic indicators were collected from two of the most cited studies, and were compared at the two refinements. In addition, correlation to rupture was investigated. RESULTS: Most of the hemodynamic indicators were very well resolved at the coarser resolutions, correlating with the finest resolution with a correlation coefficient >0.95. The oscillatory shear index (OSI) had the lowest correlation coefficient of 0.83. A logarithmic Bland-Altman plot revealed noticeable variations in the proportion of the aneurysm under low shear, as well as in spatial and temporal gradients not captured by the correlation alone. CONCLUSION: Statistically, hemodynamic indicators agree well across the different resolutions studied here. However, there are clear outliers visible in several of the hemodynamic indicators, which suggests that special care should be taken when considering individual assessment.

Quantifying Pelvic Periprosthetic Bone Remodeling Using Dual-Energy X-Ray Absorptiometry Region-Free Analysis.

The gold standard tool for measuring periprosthetic bone mineral density (BMD) is dual-energy X-ray absorptiometry (DXA). However, resolution of the method is limited due to the aggregation of pixel data into large regions of interest for clinical and statistical analysis. We have previously validated a region-free analysis method (DXA-RFA) for quantitating BMD change at the pixel level around femoral prostheses. Here, we applied the DXA-RFA method to the pelvis, and quantitated its precision in this setting using repeated DXA scans taken on the same day after repositioning in 29 patients after total hip arthroplasty. Scans were semiautomatically segmented using edge detection, intensity thresholding, and morphologic operations, and elastically registered to a common template generated through generalized Procrustes analysis. Pixel-wise BMD precision between repeated scans was expressed as a coefficient of variation %. Longitudinal BMD change was assessed in an independent group of 24 patients followed up for 260 wk. DXA-RFA spatial resolution of 0.31 mm2 provided approximately 12,500 data points per scan. The median data-point precision was 17.8% (interquartile range 14.3%-22.7%). The anatomic distribution of the precision errors showed poorer precision at the bone borders and superior precision to the obturator foramen. Evaluation of longitudinal BMD showed focal BMD change at 260 wk of -26.8% adjacent to the prosthesis-bone interface (1% of bone map area). In contrast, BMD change of +39.0% was observed at the outer aspect of the ischium (3% of bone map area). Pelvic DXA-RFA is less precise than conventional DXA analysis. However, it is sensitive for detecting local BMD change events in groups of patients, and provides a novel tool for quantitating local bone mass after joint replacement. Using this method, we were able to resolve BMD change over small areas adjacent to the implant-bone interface and in the ischial region over 260 wk after total hip arthroplasty.

Quantitating the effect of prosthesis design on femoral remodeling using high-resolution region-free densitometric analysis (DXA-RFA).

Dual energy X-ray absorptiometry (DXA) is the reference standard method used to study bone mineral density (BMD) after total hip arthroplasty (THA). However, the subtle, spatially complex changes in bone mass due to strain-adaptive bone remodeling relevant to different prosthesis designs are not readily resolved using conventional DXA analysis. DXA region free analysis (DXA RFA) is a novel computational image analysis technique that provides a high-resolution quantitation of periprosthetic BMD. Here, we applied the technique to quantitate the magnitude and areal size of periprosthetic BMD changes using scans acquired during two previous randomized clinical trials (2004 to 2009); one comparing three cemented prosthesis design geometries, and the other comparing a hip resurfacing versus a conventional cementless prosthesis. DXA RFA resolved subtle differences in magnitude and area of bone remodeling between prosthesis designs not previously identified in conventional DXA analyses. A mean bone loss of 10.3%, 12.1%, and 11.1% occurred for the three cemented prostheses within a bone area fraction of 14.8%, 14.4%, and 6.2%, mostly within the lesser trochanter (p < 0.001). For the cementless prosthesis, a diffuse pattern of bone loss (-14.3%) was observed at the shaft of femur in a small area fraction of 0.6% versus no significant bone loss for the hip resurfacing prosthesis (p < 0.001). BMD increases were observed consistently at the greater trochanter for all prostheses except the hip-resurfacing prosthesis, where BMD increase was widespread across the metaphysis (p < 0.001). DXA RFA provides high-resolution insights into the effect of prosthesis design on the local strain environment in bone. © 2017 The Authors Journal of Orthopaedic Research published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society. J Orthop Res 35:2203-2210, 2017.

Multiresolution eXtended Free-Form Deformations (XFFD) for non-rigid registration with discontinuous transforms.

Image registration is an essential technique to obtain point correspondences between anatomical structures from different images. Conventional non-rigid registration methods assume a continuous and smooth deformation field throughout the image. However, the deformation field at the interface of different organs is not necessarily continuous, since the organs may slide over or separate from each other. Therefore, imposing continuity and smoothness ubiquitously would lead to artifacts and increased errors near the discontinuity interface. In computational mechanics, the eXtended Finite Element Method (XFEM) was introduced to handle discontinuities without using computational meshes that conform to the discontinuity geometry. Instead, the interpolation bases themselves were enriched with discontinuous functional terms. Borrowing this concept, we propose a multiresolution eXtented Free-Form Deformation (XFFD) framework that seamlessly integrates within and extends the standard Free-Form Deformation (FFD) approach. Discontinuities are incorporated by enriching the B-spline basis functions coupled with extra degrees of freedom, which are only introduced near the discontinuity interface. In contrast with most previous methods, restricted to sliding motion, no ad hoc penalties or constraints are introduced to reduce gaps and overlaps. This allows XFFD to describe more general discontinuous motions. In addition, we integrate XFFD into a rigorously formulated multiresolution framework by introducing an exact parameter upsampling method. The proposed method has been evaluated in two publicly available datasets: 4D pulmonary CT images from the DIR-Lab dataset and 4D CT liver datasets. The XFFD achieved a Target Registration Error (TRE) of 1.17 ± 0.85 mm in the DIR-lab dataset and 1.94 ± 1.01 mm in the liver dataset, which significantly improves on the performance of the state-of-the-art methods handling discontinuities.

Intervertebral disc classification by its degree of degeneration from T2-weighted magnetic resonance images.

PURPOSE: The primary goal of this article is to achieve an automatic and objective method to compute the Pfirrmann's degeneration grade of intervertebral discs (IVD) from MRI. This grading system is used in the diagnosis and management of patients with low back pain (LBP). In addition, biomechanical models, which are employed to assess the treatment on patients with LBP, require this grading value to compute proper material properties. MATERIALS AND METHODS: T2-weighted MR images of 48 patients were employed in this work. The 240 lumbar IVDs were divided into a training set (140) and a testing set (100). Three experts manually classified the whole set of IVDs using the Pfirrmann's grading system and the ground truth was selected as the most voted value among them. The developed method employs active contour models to delineate the boundaries of the IVD. Subsequently, the classification is achieved using a trained Neural Network (NN) with eight designed features that contain shape and intensity information of the IVDs. RESULTS: The classification method was evaluated using the testing set, resulting in a mean specificity (95.5 %) and sensitivity (87.3 %) comparable to those of every expert with respect to the ground truth. CONCLUSIONS: Our results show that the automatic method and humans perform equally well in terms of the classification accuracy. However, human annotations have inherent inter- and intra-observer variabilities, which lead to inconsistent assessments. In contrast, the proposed automatic method is objective, being only dependent on the input MRI.

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.

Statistical Interspace Models (SIMs): Application to Robust 3D Spine Segmentation.

Statistical shape models (SSM) are used to introduce shape priors in the segmentation of medical images. However, such models require large training datasets in the case of multi-object structures, since it is required to obtain not only the individual shape variations but also the relative position and orientation among objects. A solution to overcome this limitation is to model each individual shape independently. However, this approach does not take into account the relative position, orientations and shapes among the parts of an articulated object, which may result in unrealistic geometries, such as with object overlaps. In this article, we propose a new Statistical Model, the Statistical Interspace Model (SIM), which provides information about the interaction of all the individual structures by modeling the interspace between them. The SIM is described using relative position vectors between pair of points that belong to different objects that are facing each other. These vectors are divided into their magnitude and direction, each of these groups modeled as independent manifolds. The SIM was included in a segmentation framework that contains an SSM per individual object. This framework was tested using three distinct types of datasets of CT images of the spine. Results show that the SIM completely eliminated the inter-process overlap while improving the segmentation accuracy.