A cross-center smoothness prior for variational Bayesian brain tissue segmentation

Suppose one is faced with the challenge of tissue segmentation in MR images, without annotators at their center to provide labeled training data. One option is to go to another medical center for a trained classifier. Sadly, tissue classifiers do not generalize well across centers due to voxel intensity shifts caused by center-specific acquisition protocols. However, certain aspects of segmentations, such as spatial smoothness, remain relatively consistent and can be learned separately. Here we present a smoothness prior that is fit to segmentations produced at another medical center. This informative prior is presented to an unsupervised Bayesian model. The model clusters the voxel intensities, such that it produces segmentations that are similarly smooth to those of the other medical center. In addition, the unsupervised Bayesian model is extended to a semi-supervised variant, which needs no visual interpretation of clusters into tissues.Comment: 12 pages, 2 figures, 1 table. Accepted to the International Conference on Information Processing in Medical Imaging (2019

Asymmetric Image-Template Registration

Authors Manuscript received: 2010 May 4. 12th International Conference, London, UK, September 20-24, 2009, Proceedings, Part IA natural requirement in pairwise image registration is that the resulting deformation is independent of the order of the images. This constraint is typically achieved via a symmetric cost function and has been shown to reduce the effects of local optima. Consequently, symmetric registration has been successfully applied to pairwise image registration as well as the spatial alignment of individual images with a template. However, recent work has shown that the relationship between an image and a template is fundamentally asymmetric. In this paper, we develop a method that reconciles the practical advantages of symmetric registration with the asymmetric nature of image-template registration by adding a simple correction factor to the symmetric cost function. We instantiate our model within a log-domain diffeomorphic registration framework. Our experiments show exploiting the asymmetry in image-template registration improves alignment in the image coordinates.NAMIC (NIH NIBIB NAMIC U54-EB005149)NAC (NIH NCRR NAC P41- RR13218)mBIRN (NIH NCRR mBIRN U24-RR021382)NIH NINDS (R01-NS051826 Grant)National Science Foundation (U.S.) (CAREER Grant 0642971)NIBIB (R01 EB001550)NIBIB (R01EB006758)NCRR (R01 RR16594-01A1)NCRR (P41-RR14075)NINDS (R01 NS052585-01)Singapore. Agency for Science, Technology and Researc

Symmetric Log-Domain Diffeomorphic Registration: A Demons-based Approach

pmid 18979814International audienceModern morphometric studies use non-linear image registration to compare anatomies and perform group analysis. Recently, log-Euclidean approaches have contributed to promote the use of such computational anatomy tools by permitting simple computations of statistics on a rather large class of invertible spatial transformations. In this work, we propose a non-linear registration algorithm perfectly fit for log-Euclidean statistics on diffeomorphisms. Our algorithm works completely in the log-domain, i.e. it uses a stationary velocity field. This implies that we guarantee the invertibility of the deformation and have access to the true inverse transformation. This also means that our output can be directly used for log-Euclidean statistics without relying on the heavy computation of the log of the spatial transformation. As it is often desirable, our algorithm is symmetric with respect to the order of the input images. Furthermore, we use an alternate optimization approach related to Thirion's demons algorithm to provide a fast non-linear registration algorithm. First results show that our algorithm outperforms both the demons algorithm and the recently proposed diffeomorphic demons algorithm in terms of accuracy of the transformation while remaining computationally efficient

Post traumatic brain perfusion SPECT analysis using reconstructed ROI maps of radioactive microsphere derived cerebral blood flow and statistical parametric mapping

<p>Abstract</p> <p>Background</p> <p>Assessment of cerebral blood flow (CBF) by SPECT could be important in the management of patients with severe traumatic brain injury (TBI) because changes in regional CBF can affect outcome by promoting edema formation and intracranial pressure elevation (with cerebral hyperemia), or by causing secondary ischemic injury including post-traumatic stroke. The purpose of this study was to establish an improved method for evaluating regional CBF changes after TBI in piglets.</p> <p>Methods</p> <p>The focal effects of moderate traumatic brain injury (TBI) on cerebral blood flow (CBF) by SPECT cerebral blood perfusion (CBP) imaging in an animal model were investigated by parallelized statistical techniques. Regional CBF was measured by radioactive microspheres and by SPECT 2 hours after injury in sham-operated piglets versus those receiving severe TBI by fluid-percussion injury to the left parietal lobe. Qualitative SPECT CBP accuracy was assessed against reference radioactive microsphere regional CBF measurements by map reconstruction, registration and smoothing. Cerebral hypoperfusion in the test group was identified at the voxel level using statistical parametric mapping (SPM).</p> <p>Results</p> <p>A significant area of hypoperfusion (P < 0.01) was found as a response to the TBI. Statistical mapping of the reference microsphere CBF data confirms a focal decrease found with SPECT and SPM.</p> <p>Conclusion</p> <p>The suitability of SPM for application to the experimental model and ability to provide insight into CBF changes in response to traumatic injury was validated by the SPECT SPM result of a decrease in CBP at the left parietal region injury area of the test group. Further study and correlation of this characteristic lesion with long-term outcomes and auxiliary diagnostic modalities is critical to developing more effective critical care treatment guidelines and automated medical imaging processing techniques.</p

Preconditioned ADMM with nonlinear operator constraint

We are presenting a modification of the well-known Alternating Direction Method of Multipliers (ADMM) algorithm with additional preconditioning that aims at solving convex optimisation problems with nonlinear operator constraints. Connections to the recently developed Nonlinear Primal-Dual Hybrid Gradient Method (NL-PDHGM) are presented, and the algorithm is demonstrated to handle the nonlinear inverse problem of parallel Magnetic Resonance Imaging (MRI)

Maternal prepregnancy body mass index and offspring white matter microstructure: results from three birth cohorts

Prepregnancy maternal obesity is a global health problem and has been associated with offspring metabolic and mental ill-health. However, there is a knowledge gap in understanding potential neurobiological factors related to these associations. This study explored the relation between maternal prepregnancy body mass index (BMI) and offspring brain white matter microstructure at the age of 6, 10, and 26 years in three independent cohorts. Maternal BMI was associated with higher FA and lower MD in multiple brain tracts in offspring aged 10 and 26 years, but not at 6 years of age. Future studies should examine whether our observations can be replicated and explore the potential causal nature of the findings.This work was supported by the European Union’s Horizon 2020 research and innovation program [grant agreement no. 633595 DynaHEALTH] and no. 733206 LifeCycle], the Netherlands Organization for Health Research and Development [ZONMW Vici project 016.VICI.170.200]. The PREOBE cohort was funded by Spanish Ministry of Innovation and Science. Junta de Andalucía: Excellence Projects (P06-CTS-02341) and Spanish Ministry of Economy and Competitiveness (BFU2012-40254-C03-01). The first phase of the Generation R Study is made possible by financial support from the Erasmus Medical Centre, the Erasmus University, and the Netherlands Organization for Health Research and Development (ZonMW, grant ZonMW Geestkracht 10.000.1003). The Northern Finland Birth Cohort 1986 is funded by University of Oulu, University Hospital of Oulu, Academy of Finland (EGEA), Sigrid Juselius Foundation, European Commission (EURO-BLCS, Framework 5 award QLG1-CT-2000-01643), NIH/NIMH (5R01MH63706:02

Use of brain MRI atlases to determine boundaries of age-related pathology: the importance of statistical method

Neurodegenerative disease diagnoses may be supported by the comparison of an individual patient's brain magnetic resonance image (MRI) with a voxel-based atlas of normal brain MRI. Most current brain MRI atlases are of young to middle-aged adults and parametric, e.g., mean ± standard deviation (SD); these atlases require data to be Gaussian. Brain MRI data, e.g., grey matter (GM) proportion images, from normal older subjects are apparently not Gaussian. We created a nonparametric and a parametric atlas of the normal limits of GM proportions in older subjects and compared their classifications of GM proportions in Alzheimer's disease (AD) patients.Using publicly available brain MRI from 138 normal subjects and 138 subjects diagnosed with AD (all 55-90 years), we created: a mean ± SD atlas to estimate parametrically the percentile ranks and limits of normal ageing GM; and, separately, a nonparametric, rank order-based GM atlas from the same normal ageing subjects. GM images from AD patients were then classified with respect to each atlas to determine the effect statistical distributions had on classifications of proportions of GM in AD patients.The parametric atlas often defined the lower normal limit of the proportion of GM to be negative (which does not make sense physiologically as the lowest possible proportion is zero). Because of this, for approximately half of the AD subjects, 25-45% of voxels were classified as normal when compared to the parametric atlas; but were classified as abnormal when compared to the nonparametric atlas. These voxels were mainly concentrated in the frontal and occipital lobes.To our knowledge, we have presented the first nonparametric brain MRI atlas. In conditions where there is increasing variability in brain structure, such as in old age, nonparametric brain MRI atlases may represent the limits of normal brain structure more accurately than parametric approaches. Therefore, we conclude that the statistical method used for construction of brain MRI atlases should be selected taking into account the population and aim under study. Parametric methods are generally robust for defining central tendencies, e.g., means, of brain structure. Nonparametric methods are advisable when studying the limits of brain structure in ageing and neurodegenerative disease

Unsupervised clustering using diffusion maps for local shape modelling

Understanding the biological variability of anatomical objects, is essential for statistical shape analysis and to distinguish between healthy and pathological structures. Statistical Shape Modelling (SSM) can be used to analyse the shapes of sub-structures aiming to describe their variation across individual objects and between groups of them [1]. However, when the shapes exhibit; self-similarity or are intrinsically fractal, such as often encountered in biomedical problems, global shape models result in highly non-linear shape spaces and it can be difficult; to determine a compact set; of modes of variation. In this work, we present, a method for local shape, modelling and analysis that uses Diffusion Maps [2] for non-linear, spectral clustering to build a set of linear shape spaces for such analysis. The method uses a curvature scale-space (CSS) description of shape to partition them into sets of self-similar parts and these are then linearly mixed to more compactly model the global shape

Non-parametric Diffeomorphic Image Registration with the Demons Algorithm: Non-parametric Diffeomorphic Image Registration with the Demons Algorithm

International audienceWe propose a non-parametric diffeomorphic image registration algorithm based on Thirion's demons algorithm. The demons algorithm can be seen as an optimization procedure on the entire space of displacement fields. The main idea of our algorithm is to adapt this procedure to a space of diffeomorphic transformations. In contrast to many diffeomorphic registration algorithms, our solution is computationally efficient since in practice it only replaces an addition of free form deformations by a few compositions. Our experiments show that in addition to being diffeomorphic, our algorithm provides results that are similar to the ones from the demons algorithm but with transformations that are much smoother and closer to the true ones in terms of Jacobians