Adaptive fuzzy segmentation of magnetic resonance images

Taste dysfunction in multiple sclerosis

Empirical studies of taste function in multiple sclerosis (MS) are rare. Moreover, a detailed assessment of whether quantitative measures of taste function correlate with the punctate and patchy myelin-related lesions found throughout the CNS of MS patients has not been made. We administered a 96-trial test of sweet (sucrose), sour (citric acid), bitter (caffeine) and salty (NaCl) taste perception to the left and right anterior (CN VII) and posterior (CN IX) tongue regions of 73 MS patients and 73 matched controls. The number and volume of lesions were assessed using quantitative MRI in 52 brain regions of 63 of the MS patients. Taste identification scores were significantly lower in the MS patients for sucrose (p = 0.0002), citric acid (p = 0.0001), caffeine (p = 0.0372) and NaCl (p = 0.0004) and were present in both anterior and posterior tongue regions. The percent of MS patients with identification scores falling below the 5th percentile of controls was 15.07 % for caffeine, 21.9 % for citric acid, 24.66 % for sucrose, and 31.50 % for NaCl. Such scores were inversely correlated with lesion volumes in the temporal, medial frontal, and superior frontal lobes, and with the number of lesions in the left and right superior frontal lobes, right anterior cingulate gyrus, and left parietal operculum. Regardless of the subject group, women outperformed men on the taste measures. These findings indicate that a sizable number of MS patients exhibit taste deficits that are associated with MS-related lesions throughout the brain.

Effects of Hormone Therapy on Brain Volumes Changes of Postmenopausal Women Revealed by Optimally-Discriminative Voxel-Based Morphometry

Backgrounds

The Women's Health Initiative Memory Study Magnetic Resonance Imaging (WHIMS-MRI) provides an opportunity to evaluate how menopausal hormone therapy (HT) affects the structure of older women’s brains. Our earlier work based on region of interest (ROI) analysis demonstrated potential structural changes underlying adverse effects of HT on cognition. However, the ROI-based analysis is limited in statistical power and precision, and cannot provide fine-grained mapping of whole-brain changes.

Methods

We aimed to identify local structural differences between HT and placebo groups from WHIMS-MRI in a whole-brain refined level, by using a novel method, named Optimally-Discriminative Voxel-Based Analysis (ODVBA). ODVBA is a recently proposed imaging pattern analysis approach for group comparisons utilizing a spatially adaptive analysis scheme to accurately locate areas of group differences, thereby providing superior sensitivity and specificity to detect the structural brain changes over conventional methods.

Results

Women assigned to HT treatments had significant Gray Matter (GM) losses compared to the placebo groups in the anterior cingulate and the adjacent medial frontal gyrus, and the orbitofrontal cortex, which persisted after multiple comparison corrections. There were no regions where HT was significantly associated with larger volumes compared to placebo, although a trend of marginal significance was found in the posterior cingulate cortical area. The CEE-Alone and CEE+MPA groups, although compared with different placebo controls, demonstrated similar effects according to the spatial patterns of structural changes.

Conclusions

HT had adverse effects on GM volumes and risk for cognitive impairment and dementia in older women. These findings advanced our understanding of the neurobiological underpinnings of HT effects.

3D multimodal spatial fuzzy segmentation of intramuscular connective and adipose tissue from ultrashort TE MR images of calf muscle

PURPOSE

To develop and evaluate an automated algorithm to segment intramuscular adipose (IMAT) and connective (IMCT) tissue from musculoskeletal MRI images acquired with a dual echo Ultrashort TE (UTE) sequence.

THEORY AND METHODS

The dual echo images and calculated structure tensor images are the inputs to the multichannel fuzzy cluster mean (MCFCM) algorithm. Modifications to the basic multichannel fuzzy cluster mean include an adaptive spatial term and bias shading correction. The algorithm was tested on digital phantoms simulating IMAT/IMCT tissue under varying conditions of image noise and bias and on ten subjects with varying amounts of IMAT/IMCT.

RESULTS

The MCFCM including the adaptive spatial term and bias shading correction performed better than the original MCFCM and adaptive spatial MCFCM algorithms. IMAT/IMCT was segmented from the unsmoothed simulated phantom data with a mean Dice coefficient of 0.933 ±0.001 when contrast-to-noise (CNR) was 140 and bias was varied between 30% and 65%. The algorithm yielded accurate in vivo segmentations of IMAT/IMCT with a mean Dice coefficient of 0.977 ±0.066.

CONCLUSION

The proposed algorithm is completely automated and yielded accurate segmentation of intramuscular adipose and connective tissue in the digital phantom and in human calf data. Magn Reson Med 77:870-883, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Temporal filtering of longitudinal brain magnetic resonance images for consistent segmentation

Longitudinal analysis of magnetic resonance images of the human brain provides knowledge of brain changes during both normal aging as well as the progression of many diseases. Previous longitudinal segmentation methods have either ignored temporal information or have incorporated temporal consistency constraints within the algorithm. In this work, we assume that some anatomical brain changes can be explained by temporal transitions in image intensities. Once the images are aligned in the same space, the intensities of each scan at the same voxel constitute a temporal (or 4D) intensity trend at that voxel. Temporal intensity variations due to noise or other artifacts are corrected by a 4D intensity-based filter that smooths the intensity values where appropriate, while preserving real anatomical changes such as atrophy. Here smoothing refers to removal of sudden changes or discontinuities in intensities. Images processed with the 4D filter can be used as a pre-processing step to any segmentation method. We show that such a longitudinal pre-processing step produces robust and consistent longitudinal segmentation results, even when applying 3D segmentation algorithms. We compare with state-of-the-art 4D segmentation algorithms. Specifically, we experimented on three longitudinal datasets containing 4-12 time-points, and showed that the 4D temporal filter is more robust and has more power in distinguishing between healthy subjects and those with dementia, mild cognitive impairment, as well as different phenotypes of multiple sclerosis.

Rough Sets and Stomped Normal Distribution for Simultaneous Segmentation and Bias Field Correction in Brain MR Images

The segmentation of brain MR images into different tissue classes is an important task for automatic image analysis technique, particularly due to the presence of intensity inhomogeneity artifact in MR images. In this regard, this paper presents a novel approach for simultaneous segmentation and bias field correction in brain MR images. It integrates judiciously the concept of rough sets and the merit of a novel probability distribution, called stomped normal (SN) distribution. The intensity distribution of a tissue class is represented by SN distribution, where each tissue class consists of a crisp lower approximation and a probabilistic boundary region. The intensity distribution of brain MR image is modeled as a mixture of finite number of SN distributions and one uniform distribution. The proposed method incorporates both the expectation-maximization and hidden Markov random field frameworks to provide an accurate and robust segmentation. The performance of the proposed approach, along with a comparison with related methods, is demonstrated on a set of synthetic and real brain MR images for different bias fields and noise levels.

A Robust Energy Minimization Algorithm for MS-Lesion Segmentation

The detection of multiple sclerosis lesion is important for many neuroimaging studies. In this paper, a new automatic robust algorithm for lesion segmentation based on MR images is proposed. This method takes full advantage of the decomposition of MR images into the true image that characterizes a physical property of the tissues and the bias field that accounts for the intensity inhomogeneity. An energy function is defined in term of the property of true image and bias field. The energy minimization is proposed for seeking the optimal segmentation result of lesions and white matter. Then postprocessing operations is used to select the most plausible lesions in the obtained hyperintense signals. The experimental results show that our approach is effective and robust for the lesion segmentation.

Automatic Contrast Enhancement of Brain MR Images Using Hierarchical Correlation Histogram Analysis

Parkinson’s disease is a progressive neurodegenerative disorder that has a higher probability of occurrence in middle-aged and older adults than in the young. With the use of a computer-aided diagnosis (CAD) system, abnormal cell regions can be identified, and this identification can help medical personnel to evaluate the chance of disease. This study proposes a hierarchical correlation histogram analysis based on the grayscale distribution degree of pixel intensity by constructing a correlation histogram, that can improves the adaptive contrast enhancement for specific objects. The proposed method produces significant results during contrast enhancement preprocessing and facilitates subsequent CAD processes, thereby reducing recognition time and improving accuracy. The experimental results show that the proposed method is superior to existing methods by using two estimation image quantitative methods of PSNR and average gradient values. Furthermore, the edge information pertaining to specific cells can effectively increase the accuracy of the results.

Improved R2* liver iron concentration assessment using a novel fuzzy c-mean clustering scheme

BACKGROUND

In thalassemia patients, R2* liver iron concentration (LIC) measurement is a common clinical tool for assessing iron overload and for determining necessary chelator dose and evaluating its efficacy. Despite the importance of accurate LIC measurement, existing methods suffer from LIC variability, especially at the severe iron overload range due to inclusion of vessel parts in LIC calculation. In this study, we build upon previous Fuzzy C-Mean (FCM) clustering work to formulate a scheme with superior performance in segmenting vessel pixels from the parenchyma. Our method (MIX-FCM) combines our novel 2D-FCM with the existing 1D-FCM algorithm. This study further assessed possible optimal clustering parameters (OP scheme) and proposed a semi-automatic (SA) scheme for routine clinical application.

METHODS

Segmentation of liver parenchyma and vessels was performed on T2* images and their LIC maps in 196 studies from 147 thalassemia major patients. We used manual segmentation as the reference. 1D-FCM clustering was performed on the acquired image alone and 2D-FCM used both the acquired image and its LIC data. To execute the MIX-FCM method, the best outcome (OP-MIX-FCM) was selected from the aforementioned methods and was compared to the SA-MIX-FCM scheme. We used the percent value of the normalized interquartile range (nIQR) to its median to evaluate the variability of all methods.

RESULTS

2D-FCM clustering is more effective than 1D-FCM clustering at the severe overload range only, but inferior for other ranges (where 1D-FCM provides suitable results). This complementary performance between the two methods allows MIX-FCM to improve results for all ranges. OP-MIX-FCM clustering error was 2.1 ± 2.3%, compared with 10.3 ± 9.9% and 7.0 ± 11.9% from 1D- and 2D-FCM clustering, respectively. SA-MIX-FCM result was comparable to OP-MIX-FCM result, with both schemes showing ability to decrease overall nIQR by approximately 30%.

CONCLUSION

Our proposed 2D-FCM algorithm is not as superior to 1D-FCM as hypothesized. In contrast, our MIX-FCM method benefits from the best of both methods to obtain the highest segmentation accuracy at all ranges. Moreover, segmentation accuracy of the practical scheme (SA-MIX-FCM) is comparable to segmentation accuracy of the reference scheme (OP-MIX-FCM). Finally, we confirmed that segmentation is crucial to improving LIC assessments, especially at the severe iron overload range.

Improved Fuzzy C-Means based Particle Swarm Optimization (PSO) initialization and outlier rejection with level set methods for MR brain image segmentation

In this paper, a new image segmentation method based on Particle Swarm Optimization (PSO) and outlier rejection combined with level set is proposed. A traditional approach to the segmentation of Magnetic Resonance (MR) images is the Fuzzy C-Means (FCM) clustering algorithm. The membership function of this conventional algorithm is sensitive to the outlier and does not integrate the spatial information in the image. The algorithm is very sensitive to noise and in-homogeneities in the image, moreover, it depends on cluster centers initialization. To improve the outlier rejection and to reduce the noise sensitivity of conventional FCM clustering algorithm, a novel extended FCM algorithm for image segmentation is presented. In general, in the FCM algorithm the initial cluster centers are chosen randomly, with the help of PSO algorithm the clusters centers are chosen optimally. Our algorithm takes also into consideration the spatial neighborhood information. These a priori are used in the cost function to be optimized. For MR images, the resulting fuzzy clustering is used to set the initial level set contour. The results confirm the effectiveness of the proposed algorithm.

Resource atlases for multi-atlas brain segmentations with multiple ontology levels based on T1-weighted MRI

Technologies for multi-atlas brain segmentation of T1-weighted MRI images have rapidly progressed in recent years, with highly promising results. This approach, however, relies on a large number of atlases with accurate and consistent structural identifications. Here, we introduce our atlas inventories (n=90), which cover ages 4-82years with unique hierarchical structural definitions (286 structures at the finest level). This multi-atlas library resource provides the flexibility to choose appropriate atlases for various studies with different age ranges and structure-definition criteria. In this paper, we describe the details of the atlas resources and demonstrate the improved accuracy achievable with a dynamic age-matching approach, in which atlases that most closely match the subject's age are dynamically selected. The advanced atlas creation strategy, together with atlas pre-selection principles, is expected to support the further development of multi-atlas image segmentation.

A Variational Approach to Simultaneous Image Segmentation and Bias Correction

This paper presents a novel variational approach for simultaneous estimation of bias field and segmentation of images with intensity inhomogeneity. We model intensity of inhomogeneous objects to be Gaussian distributed with different means and variances, and then introduce a sliding window to map the original image intensity onto another domain, where the intensity distribution of each object is still Gaussian but can be better separated. The means of the Gaussian distributions in the transformed domain can be adaptively estimated by multiplying the bias field with a piecewise constant signal within the sliding window. A maximum likelihood energy functional is then defined on each local region, which combines the bias field, the membership function of the object region, and the constant approximating the true signal from its corresponding object. The energy functional is then extended to the whole image domain by the Bayesian learning approach. An efficient iterative algorithm is proposed for energy minimization, via which the image segmentation and bias field correction are simultaneously achieved. Furthermore, the smoothness of the obtained optimal bias field is ensured by the normalized convolutions without extra cost. Experiments on real images demonstrated the superiority of the proposed algorithm to other state-of-the-art representative methods.

Interpreting support vector machine models for multivariate group wise analysis in neuroimaging

Machine learning based classification algorithms like support vector machines (SVMs) have shown great promise for turning a high dimensional neuroimaging data into clinically useful decision criteria. However, tracing imaging based patterns that contribute significantly to classifier decisions remains an open problem. This is an issue of critical importance in imaging studies seeking to determine which anatomical or physiological imaging features contribute to the classifier's decision, thereby allowing users to critically evaluate the findings of such machine learning methods and to understand disease mechanisms. The majority of published work addresses the question of statistical inference for support vector classification using permutation tests based on SVM weight vectors. Such permutation testing ignores the SVM margin, which is critical in SVM theory. In this work we emphasize the use of a statistic that explicitly accounts for the SVM margin and show that the null distributions associated with this statistic are asymptotically normal. Further, our experiments show that this statistic is a lot less conservative as compared to weight based permutation tests and yet specific enough to tease out multivariate patterns in the data. Thus, we can better understand the multivariate patterns that the SVM uses for neuroimaging based classification.

Patient-specific biomechanical model as whole-body CT image registration tool

Whole-body computed tomography (CT) image registration is important for cancer diagnosis, therapy planning and treatment. Such registration requires accounting for large differences between source and target images caused by deformations of soft organs/tissues and articulated motion of skeletal structures. The registration algorithms relying solely on image processing methods exhibit deficiencies in accounting for such deformations and motion. We propose to predict the deformations and movements of body organs/tissues and skeletal structures for whole-body CT image registration using patient-specific non-linear biomechanical modelling. Unlike the conventional biomechanical modelling, our approach for building the biomechanical models does not require time-consuming segmentation of CT scans to divide the whole body into non-overlapping constituents with different material properties. Instead, a Fuzzy C-Means (FCM) algorithm is used for tissue classification to assign the constitutive properties automatically at integration points of the computation grid. We use only very simple segmentation of the spine when determining vertebrae displacements to define loading for biomechanical models. We demonstrate the feasibility and accuracy of our approach on CT images of seven patients suffering from cancer and aortic disease. The results confirm that accurate whole-body CT image registration can be achieved using a patient-specific non-linear biomechanical model constructed without time-consuming segmentation of the whole-body images.

Modeling clustered activity increase in amyloid-beta positron emission tomographic images with statistical descriptors

Background

Amyloid-beta (Aβ) imaging with positron emission tomography (PET) holds promise for detecting the presence of Aβ plaques in the cortical gray matter. Many image analyses focus on regional average measurements of tracer activity distribution; however, considerable additional information is available in the images. Metrics that describe the statistical properties of images, such as the two-point correlation function (S₂), have found wide applications in astronomy and materials science. S₂ provides a detailed characterization of spatial patterns in images typically referred to as clustering or flocculence. The objective of this study was to translate the two-point correlation method into Aβ-PET of the human brain using 11C-Pittsburgh compound B (11C-PiB) to characterize longitudinal changes in the tracer distribution that may reflect changes in Aβ plaque accumulation.

Methods

We modified the conventional S₂ metric, which is primarily used for binary images and formulated a weighted two-point correlation function (wS₂) to describe nonbinary, real-valued PET images with a single statistical function. Using serial 11C-PiB scans, we calculated wS₂ functions from two-dimensional PET images of different cortical regions as well as three-dimensional data from the whole brain. The area under the wS₂ functions was calculated and compared with the mean/median of the standardized uptake value ratio (SUVR). For three-dimensional data, we compared the area under the wS₂ curves with the subjects’ cerebrospinal fluid measures.

Results

Overall, the longitudinal changes in wS₂ correlated with the increase in mean SUVR but showed lower variance. The whole brain results showed a higher inverse correlation between the cerebrospinal Aβ and wS₂ than between the cerebrospinal Aβ and SUVR mean/median. We did not observe any confounding of wS₂ by region size or injected dose.

Conclusion

The wS₂ detects subtle changes and provides additional information about the binding characteristics of radiotracers and Aβ accumulation that are difficult to verify with mean SUVR alone.

Lung segmentation on standard and mobile chest radiographs using oriented Gaussian derivatives filter

Background

Unsupervised lung segmentation method is one of the mandatory processes in order to develop a Content Based Medical Image Retrieval System (CBMIRS) of CXR. The purpose of the study is to present a robust solution for lung segmentation of standard and mobile chest radiographs using fully automated unsupervised method.

Methods

The novel method is based on oriented Gaussian derivatives filter with seven orientations, combined with Fuzzy C-Means (FCM) clustering and thresholding to refine the lung region. In addition, a new algorithm to automatically generate a threshold value for each Gaussian response is also proposed. The algorithms are applied to both PA and AP chest radiographs from both public JSRT dataset and our private datasets from collaborative hospital. Two pre-processing blocks are introduced to standardize the images from different machines. Comparisons with the previous works found in the literature on JSRT dataset shows that our method gives a reasonably good result. We also compare our algorithm with other unsupervised methods to provide fairly comparative measures on the performances for all datasets.

Results

Performance measures (accuracy, F-score, precision, sensitivity and specificity) for the segmentation of lung in public JSRT dataset are above 0.90 except for the overlap measure is 0.87. The standard deviations for all measures are very low, from 0.01 to 0.06. The overlap measure for the private image database is 0.81 (images from standard machine) and 0.69 (images from two mobile machines). The algorithm is fully automated and fast, with the average execution time of 12.5 s for 512 by 512 pixels resolution.

Conclusions

Our proposed method is fully automated, unsupervised, with no training or learning stage is necessary to segment the lungs taken using both a standard machine and two different mobile machines. The proposed pre-processing blocks are significantly useful to standardize the radiographs from mobile machines. The algorithm gives good performance measures, robust, and fast for the application of the CBMIRS.

Patterns of myeloarchitecture in lower limb amputees: an MRI study

Functional studies of cortical plasticity in humans suggest that the motor cortex reorganizes when the descending motor output pathway is disrupted as a result of limb amputation. The question thus arises if the underlying anatomical organization of the motor cortex is also altered in limb amputation. Owing to challenges involved in imaging the thin cerebral cortex in vivo, there is limited data available on the anatomical or morphological plasticity of the motor cortex in amputation. In this paper, we study the morphology of the primary motor cortex in four lower limb amputees with 37 or more years of amputation and four age and gender-matched controls using 0.7 mm isotropic, T₁-weighted MRI optimized to produce enhanced intracortical contrast based on myelin content. We segment the cortex into myelinated and unmyelinated gray matter. We determine the myelinated thickness which is the thickness of the well-myelinated tissue in the deeper layers of the cortex. We compare the bilateral differences in the myelinated thickness between amputees and controls. We also compare bilateral differences in cortical thickness between the two groups. Our measurements show no statistically significant difference between the amputees and controls in the myelinated thickness and in cortical thickness, in the region of the primary motor cortex representing the lower leg.

Spatial patterns of structural brain changes in type 2 diabetic patients and their longitudinal progression with intensive control of blood glucose

OBJECTIVE

Understanding the effect of diabetes as well as of alternative treatment strategies on cerebral structure is critical for the development of targeted interventions against accelerated neurodegeneration in type 2 diabetes. We investigated whether diabetes characteristics were associated with spatially specific patterns of brain changes and whether those patterns were affected by intensive versus standard glycemic treatment.

RESEARCH DESIGN AND METHODS

Using baseline MRIs of 488 participants with type 2 diabetes from the Action to Control Cardiovascular Risk in Diabetes-Memory in Diabetes (ACCORD-MIND) study, we applied a new voxel-based analysis methodology to identify spatially specific patterns of gray matter and white matter volume loss related to diabetes duration and HbA1c. The longitudinal analysis used 40-month follow-up data to evaluate differences in progression of volume loss between intensive and standard glycemic treatment arms.

RESULTS

Participants with longer diabetes duration had significantly lower gray matter volumes, primarily in certain regions in the frontal and temporal lobes. The longitudinal analysis of treatment effects revealed a heterogeneous pattern of decelerated loss of gray matter volume associated with intensive glycemic treatment. Intensive treatment decelerated volume loss, particularly in regions adjacent to those cross-sectionally associated with diabetes duration. No significant relationship between low versus high baseline HbA1c levels and brain changes was found. Finally, regions in which cognitive change was associated with longitudinal volume loss had only small overlap with regions related to diabetes duration and to treatment effects.

CONCLUSIONS

Applying advanced quantitative image pattern analysis methods on longitudinal MRI data of a large sample of patients with type 2 diabetes, we demonstrate that there are spatially specific patterns of brain changes that vary by diabetes characteristics and that the progression of gray matter volume loss is slowed by intensive glycemic treatment, particularly in regions adjacent to areas affected by diabetes.

Magnetic resonance image tissue classification using an automatic method

Background

Brain segmentation in magnetic resonance images (MRI) is an important stage in clinical studies for different issues such as diagnosis, analysis, 3-D visualizations for treatment and surgical planning. MR Image segmentation remains a challenging problem in spite of different existing artifacts such as noise, bias field, partial volume effects and complexity of the images. Some of the automatic brain segmentation techniques are complex and some of them are not sufficiently accurate for certain applications. The goal of this paper is proposing an algorithm that is more accurate and less complex).

Methods

In this paper we present a simple and more accurate automated technique for brain segmentation into White Matter, Gray Matter and Cerebrospinal fluid (CSF) in three-dimensional MR images. The algorithm’s three steps are histogram based segmentation, feature extraction and final classification using SVM. The integrated algorithm has more accurate results than what can be obtained with its individual components. To produce much more efficient segmentation method our framework captures different types of features in each step that are of special importance for MRI, i.e., distributions of tissue intensities, textural features, and relationship with neighboring voxels or spatial features.

Results

Our method has been validated on real images and simulated data, with desirable performance in the presence of noise and intensity inhomogeneities.

Conclusions

The experimental results demonstrate that our proposed method is a simple and accurate technique to define brain tissues with high reproducibility in comparison with other techniques.

Virtual Slides

The virtual slide(s) for this article can be found here: http://www.diagnosticpathology.diagnomx.eu/vs/13000_2014_207

Finding imaging patterns of structural covariance via Non-Negative Matrix Factorization

In this paper, we investigate the use of Non-Negative Matrix Factorization (NNMF) for the analysis of structural neuroimaging data. The goal is to identify the brain regions that co-vary across individuals in a consistent way, hence potentially being part of underlying brain networks or otherwise influenced by underlying common mechanisms such as genetics and pathologies. NNMF offers a directly data-driven way of extracting relatively localized co-varying structural regions, thereby transcending limitations of Principal Component Analysis (PCA), Independent Component Analysis (ICA) and other related methods that tend to produce dispersed components of positive and negative loadings. In particular, leveraging upon the well known ability of NNMF to produce parts-based representations of image data, we derive decompositions that partition the brain into regions that vary in consistent ways across individuals. Importantly, these decompositions achieve dimensionality reduction via highly interpretable ways and generalize well to new data as shown via split-sample experiments. We empirically validate NNMF in two data sets: i) a Diffusion Tensor (DT) mouse brain development study, and ii) a structural Magnetic Resonance (sMR) study of human brain aging. We demonstrate the ability of NNMF to produce sparse parts-based representations of the data at various resolutions. These representations seem to follow what we know about the underlying functional organization of the brain and also capture some pathological processes. Moreover, we show that these low dimensional representations favorably compare to descriptions obtained with more commonly used matrix factorization methods like PCA and ICA.

Technical considerations for functional magnetic resonance imaging analysis

Clinical application of functional magnetic resonance imaging (fMRI) based on blood oxygenation level-dependent (BOLD) effect has increased over the past decade because of its ability to map regional blood flow in response to brain stimulation. This mapping is primarily achieved by exploiting the BOLD effect precipitated by changes in the magnetic properties of hemoglobin. BOLD fMRI has utility in neurosurgical planning and mapping neuronal functional connectivity. Conventional echo planar imaging techniques are used to acquire stimulus-driven fMR imaging BOLD data. This article highlights technical aspects of fMRI data analysis to make it more accessible in clinical settings.

Partial volume effect modeling for segmentation and tissue classification of brain magnetic resonance images: A review

Quantitative analysis of magnetic resonance (MR) brain images are facilitated by the development of automated segmentation algorithms. A single image voxel may contain of several types of tissues due to the finite spatial resolution of the imaging device. This phenomenon, termed partial volume effect (PVE), complicates the segmentation process, and, due to the complexity of human brain anatomy, the PVE is an important factor for accurate brain structure quantification. Partial volume estimation refers to a generalized segmentation task where the amount of each tissue type within each voxel is solved. This review aims to provide a systematic, tutorial-like overview and categorization of methods for partial volume estimation in brain MRI. The review concentrates on the statistically based approaches for partial volume estimation and also explains differences to other, similar image segmentation approaches.

Partial volume estimation in brain MRI revisited

We propose a fast algorithm to estimate brain tissue concentrations from conventional T1-weighted images based on a Bayesian maximum a posteriori formulation that extends the "mixel" model developed in the 90's. A key observation is the necessity to incorporate additional prior constraints to the "mixel" model for the estimation of plausible concentration maps. Experiments on the ADNI standardized dataset show that global and local brain atrophy measures from the proposed algorithm yield enhanced diagnosis testing value than with several widely used soft tissue labeling methods.

Self-assessed performance improves statistical fusion of image labels

PURPOSE

Expert manual labeling is the gold standard for image segmentation, but this process is difficult, time-consuming, and prone to inter-individual differences. While fully automated methods have successfully targeted many anatomies, automated methods have not yet been developed for numerous essential structures (e.g., the internal structure of the spinal cord as seen on magnetic resonance imaging). Collaborative labeling is a new paradigm that offers a robust alternative that may realize both the throughput of automation and the guidance of experts. Yet, distributing manual labeling expertise across individuals and sites introduces potential human factors concerns (e.g., training, software usability) and statistical considerations (e.g., fusion of information, assessment of confidence, bias) that must be further explored. During the labeling process, it is simple to ask raters to self-assess the confidence of their labels, but this is rarely done and has not been previously quantitatively studied. Herein, the authors explore the utility of self-assessment in relation to automated assessment of rater performance in the context of statistical fusion.

METHODS

The authors conducted a study of 66 volumes manually labeled by 75 minimally trained human raters recruited from the university undergraduate population. Raters were given 15 min of training during which they were shown examples of correct segmentation, and the online segmentation tool was demonstrated. The volumes were labeled 2D slice-wise, and the slices were unordered. A self-assessed quality metric was produced by raters for each slice by marking a confidence bar superimposed on the slice. Volumes produced by both voting and statistical fusion algorithms were compared against a set of expert segmentations of the same volumes.

RESULTS

Labels for 8825 distinct slices were obtained. Simple majority voting resulted in statistically poorer performance than voting weighted by self-assessed performance. Statistical fusion resulted in statistically indistinguishable performance from self-assessed weighted voting. The authors developed a new theoretical basis for using self-assessed performance in the framework of statistical fusion and demonstrated that the combined sources of information (both statistical assessment and self-assessment) yielded statistically significant improvement over the methods considered separately.

CONCLUSIONS

The authors present the first systematic characterization of self-assessed performance in manual labeling. The authors demonstrate that self-assessment and statistical fusion yield similar, but complementary, benefits for label fusion. Finally, the authors present a new theoretical basis for combining self-assessments with statistical label fusion.

A modified method for MRF segmentation and bias correction of MR image with intensity inhomogeneity

Markov random field (MRF) model is an effective method for brain tissue classification, which has been applied in MR image segmentation for decades. However, it falls short of the expected classification in MR images with intensity inhomogeneity for the bias field is not considered in the formulation. In this paper, we propose an interleaved method joining a modified MRF classification and bias field estimation in an energy minimization framework, whose initial estimation is based on k-means algorithm in view of prior information on MRI. The proposed method has a salient advantage of overcoming the misclassifications from the non-interleaved MRF classification for the MR image with intensity inhomogeneity. In contrast to other baseline methods, experimental results also have demonstrated the effectiveness and advantages of our algorithm via its applications in the real and the synthetic MR images.

Automatic segmentation of breast MR images through a Markov random field statistical model

An algorithm dedicated to automatic segmentation of breast magnetic resonance images is presented in this paper. Our approach is based on a pipeline that includes a denoising step and statistical segmentation. The noise removal preprocessing relies on an anisotropic diffusion scheme, whereas the statistical segmentation is conducted through a Markov random field model. The continuous updating of all parameters governing the diffusion process enables automatic denoising, and the partial volume effect is also addressed during the labeling step. To assess the relevance, the Jaccard similarity coefficient was computed. Experiments were conducted on synthetic data and breast magnetic resonance images extracted from a high-risk population. The relevance of the approach for the dataset is highlighted, and we demonstrate accuracy superior to that of traditional clustering algorithms. The results emphasize the benefits of both denoising guided by input data and the inclusion of spatial dependency through a Markov random field. For example, the Jaccard coefficient for the clinical data was increased by 114%, 109%, and 140% with respect to a K-means algorithm and, respectively, for the adipose, glandular and muscle and skin components. Moreover, the agreement between the manual segmentations provided by an experienced radiologist and the automatic segmentations performed with this algorithm was good, with Jaccard coefficients equal to 0.769, 0.756, and 0.694 for the above-mentioned classes.

Analysis of spatio-temporal brain imaging patterns by Hidden Markov Models and serial MRI images

Brain changes due to development and maturation, normal aging, or degenerative disease are continuous, gradual, and variable across individuals. To quantify the individual progression of brain changes, we propose a spatio-temporal methodology based on Hidden Markov Models (HMM), and apply it on four-dimensional structural brain magnetic resonance imaging series of older individuals. First, regional brain features are extracted in order to reduce image dimensionality. This process is guided by the objective of the study or the specific imaging patterns whose progression is of interest, for example, the evaluation of Alzheimer-like patterns of brain change in normal individuals. These regional features are used in conjunction with HMMs, which aim to measure the dynamic association between brain structure changes and progressive stages of disease over time. A bagging framework is used to obtain models with good generalization capability, since in practice the number of serial scans is limited. An application of the proposed methodology was to detect individuals with the risk of developing MCI, and therefore it was tested on modeling the progression of brain atrophy patterns in older adults. With HMM models, the state-transition paths corresponding to longitudinal brain changes were constructed from two completely independent datasets, the Alzheimer Disease Neuroimaging Initiative and the Baltimore Longitudinal Study of Aging. The statistical analysis of HMM-state paths among the normal, progressive MCI, and MCI groups indicates that, HMM-state index 1 is likely to be a predictor of the conversion from cognitively normal to MCI, potentially many years before clinical symptoms become measurable.

Development of new population-averaged standard templates for spatial normalization and segmentation of MR images for postnatal piglet brains

PROPOSE

To design a set of brain templates for postnatal piglet brains based on high-resolution T1-weighted imaging for voxel-based morphometric analysis.

MATERIALS AND METHODS

Using a 3.0 T magnetic resonance (MR) scanner, a population-based whole brain template was developed by averaging forty T1 images in the brains of postnatal piglets at 38 days of age. The templates for gray and white matter, and cerebrospinal fluid were designed based on the corresponding probability maps by adapting individual data sets using statistical parametric mapping. Anatomical labeling maps were generated from labeling propagation derived from the established Pig Brain Atlas. Differences in the coordinates from four significant structural landmarks in the template, plus an additional 12 normalized images and anatomical labeling maps were measured to validate the accuracy of the registration of the template.

RESULTS

A whole brain template, a set of tissue-specific probability and anatomical labeling maps were developed. The location deviation of the four significant structural landmarks, including the anterior and posterior regions in the corpus callosum, and the left and right caudate nucleus, was found to be <0.25 cm, validating the sensitivity and resolution of the template.

CONCLUSION

A whole brain template map and a set of tissue-specific probability and anatomical labeling maps were developed to analyze the morphometric imaging of the postnatal piglet brain, an animal model of the human infant.

Sub-millimeter imaging of brain-free water for rapid volume assessment in atrophic brains

INTRODUCTION

Cerebral atrophy occurs in healthy aging, and in disease processes such as multiple sclerosis (MS), it correlates with disability accumulation. Imaging measurements of brain atrophy are commonly based on tissue segmentation, which is susceptible to classification errors and inconsistencies. High-resolution imaging techniques with strong contrast between brain parenchyma and cerebrospinal fluid (CSF) might allow fully automated, rapid, threshold-based determination of the free water in the brain. We hypothesized that total brain-free-water (BFW) volume and BFW volume expressed as a normalized fraction of the intracranial volume ("BFW fraction"), determined from heavily T2-weighted images, would be useful surrogates for cerebral atrophy and therefore would correlate with clinical measures of disability in MS.

METHODS

Whole brains of 83 MS cases and 7 healthy volunteers were imaged with a 4.7-min, heavily T2-weighted sequence on a 3T MRI scanner, acquiring 650-μm isotropic voxels. MS cases were clinically assessed on the Expanded Disability Status Scale (EDSS), Scripps Neurological Rating Scale (SNRS), Paced Auditory Serial Addition Test (PASAT), 9-Hole Peg Test (9HPT), Symbol Digit Modalities Test (SDMT), and 25-Foot Timed Walk. Twelve of the MS cases were rescanned within an average of 1.8 months to assess reproducibility. Automated calculations of BFW volume and BFW fraction were correlated with clinical measures of disability upon adjusting for age and sex. Results were compared to data from T1-based approaches (SIENAX and Lesion-TOADS).

RESULTS AND DISCUSSION

BFW volume was automatically derived from heavily T2-weighted images with no need for separate skull stripping. BFW volume and fraction had mean scan-rescan coefficients of variation of 1.5% and 1.9%, respectively, similar to the T1-based approaches tested here. BFW fraction more strongly correlated with clinical measures than T1-derived results. Among those clinical measures, modality-specific disability scores, such as SDMT and 9HPT, were more strongly associated with BFW fraction than composite measures, such as EDSS and SNRS.

CONCLUSION

The BFW method robustly estimates cerebral atrophy in an automated, fast, and reliable manner, and as such may prove a useful addition to imaging protocols for clinical practice and trials.

Multi-contrast multi-atlas parcellation of diffusion tensor imaging of the human brain

In this paper, we propose a novel method for parcellating the human brain into 193 anatomical structures based on diffusion tensor images (DTIs). This was accomplished in the setting of multi-contrast diffeomorphic likelihood fusion using multiple DTI atlases. DTI images are modeled as high dimensional fields, with each voxel exhibiting a vector valued feature comprising of mean diffusivity (MD), fractional anisotropy (FA), and fiber angle. For each structure, the probability distribution of each element in the feature vector is modeled as a mixture of Gaussians, the parameters of which are estimated from the labeled atlases. The structure-specific feature vector is then used to parcellate the test image. For each atlas, a likelihood is iteratively computed based on the structure-specific vector feature. The likelihoods from multiple atlases are then fused. The updating and fusing of the likelihoods is achieved based on the expectation-maximization (EM) algorithm for maximum a posteriori (MAP) estimation problems. We first demonstrate the performance of the algorithm by examining the parcellation accuracy of 18 structures from 25 subjects with a varying degree of structural abnormality. Dice values ranging 0.8-0.9 were obtained. In addition, strong correlation was found between the volume size of the automated and the manual parcellation. Then, we present scan-rescan reproducibility based on another dataset of 16 DTI images - an average of 3.73%, 1.91%, and 1.79% for volume, mean FA, and mean MD respectively. Finally, the range of anatomical variability in the normal population was quantified for each structure.

A fuzzy-rule-based approach for single frame super resolution

In this paper, a novel fuzzy rule-based prediction framework is developed for high-quality image zooming. In classical interpolation-based image zooming, resolution is increased by inserting pixels using certain interpolation techniques. Here, we propose a patch-based image zooming technique, where each low-resolution (LR) image patch is replaced by an estimated high-resolution (HR) patch. Since an LR patch can be generated from any of the many possible HR patches, it would be natural to develop rules to find different possible HR patches and then to combine them according to rule strength to get the estimated HR patch. Here, we generate a large number of LR–HR patch pairs from a collection of natural images, group them into different clusters, and then generate a fuzzy rule for each of these clusters. The rule parameters are also learned from these LR-HR patch pairs. As a result, an efficient mapping from LR patch space to HR patch space can be formulated. The performance of the proposed method is tested on different images, and is also compared with other representative as well as state-of-the-art image zooming techniques. Experimental results show that the proposed method is better than the competing methods and is capable of reconstructing thin lines, edges, fine details, and textures in the image efficiently.

Multiplicative intrinsic component optimization (MICO) for MRI bias field estimation and tissue segmentation

This paper proposes a new energy minimization method called multiplicative intrinsic component optimization (MICO) for joint bias field estimation and segmentation of magnetic resonance (MR) images. The proposed method takes full advantage of the decomposition of MR images into two multiplicative components, namely, the true image that characterizes a physical property of the tissues and the bias field that accounts for the intensity inhomogeneity, and their respective spatial properties. Bias field estimation and tissue segmentation are simultaneously achieved by an energy minimization process aimed to optimize the estimates of the two multiplicative components of an MR image. The bias field is iteratively optimized by using efficient matrix computations, which are verified to be numerically stable by matrix analysis. More importantly, the energy in our formulation is convex in each of its variables, which leads to the robustness of the proposed energy minimization algorithm. The MICO formulation can be naturally extended to 3D/4D tissue segmentation with spatial/sptatiotemporal regularization. Quantitative evaluations and comparisons with some popular softwares have demonstrated superior performance of MICO in terms of robustness and accuracy.

Non-locally regularized segmentation of multiple sclerosis lesion from multi-channel MRI data

Segmentation of multiple sclerosis (MS) lesion is important for many neuroimaging studies. In this paper, we propose a novel algorithm for automatic segmentation of MS lesions from multi-channel MR images (T1W, T2W and FLAIR images). The proposed method is an extension of Li et al.'s algorithm in [1], which only segments the normal tissues from T1W images. The proposed method is aimed to segment MS lesions, while normal tissues are also segmented and bias field is estimated to handle intensity inhomogeneities in the images. Another contribution of this paper is the introduction of a nonlocal means technique to achieve spatially regularized segmentation, which overcomes the influence of noise. Experimental results have demonstrated the effectiveness and advantages of the proposed algorithm.

A new multistage medical segmentation method based on superpixel and fuzzy clustering

The medical image segmentation is the key approach of image processing for brain MRI images. However, due to the visual complex appearance of image structures and the imaging characteristic, it is still challenging to automatically segment brain MRI image. A new multi-stage segmentation method based on superpixel and fuzzy clustering (MSFCM) is proposed to achieve the good brain MRI segmentation results. The MSFCM utilizes the superpixels as the clustering objects instead of pixels, and it can increase the clustering granularity and overcome the influence of noise and bias effectively. In the first stage, the MRI image is parsed into several atomic areas, namely, superpixels, and a further parsing step is adopted for the areas with bigger gray variance over setting threshold. Subsequently, designed fuzzy clustering is carried out to the fuzzy membership of each superpixel, and an iterative broadcast method based on the Butterworth function is used to redefine their classifications. Finally, the segmented image is achieved by merging the superpixels which have the same classification label. The simulated brain database from BrainWeb site is used in the experiments, and the experimental results demonstrate that MSFCM method outperforms the traditional FCM algorithm in terms of segmentation accuracy and stability for MRI image.

Multimodal MRI-based Imputation of the Aβ+ in Early Mild Cognitive Impairment

Objective

The primary goal of this study was to identify brain atrophy from structural MRI (magnetic resonance imaging) and cerebral blood flow (CBF) patterns from arterial spin labeling perfusion MRI that are best predictors of the Aβ-burden, measured as composite ¹⁸F-AV45-PET (positron emission tomography) uptake, in individuals with early mild cognitive impairment (MCI). Furthermore, another objective was to assess the relative importance of imaging modalities in classification of Aβ+/Aβ− early MCI.

Methods

Sixty-seven Alzheimer's Disease Neuroimaging Initiative (ADNI)-GO/2 participants with early MCI were included. Voxel-wise anatomical shape variation measures were computed by estimating the initial diffeomorphic mapping momenta from an unbiased control template. CBF measures normalized to average motor cortex CBF were mapped onto the template space. Using partial least squares regression, we identified the structural and CBF signatures of Aβ after accounting for normal cofounding effects of age, gender, and education.

Results

¹⁸F-AV45-positive early MCIs could be identified with 83% classification accuracy, 87% positive predictive value, and 84% negative predictive value by multidisciplinary classifiers combining demographics data, ApoE ε4-genotype, and a multimodal MRI-based Aβ score.

Interpretation

Multimodal MRI can be used to predict the amyloid status of early-MCI individuals. MRI is a very attractive candidate for the identification of inexpensive and noninvasive surrogate biomarkers of Aβ deposition. Our approach is expected to have value for the identification of individuals likely to be Aβ+ in circumstances where cost or logistical problems prevent Aβ detection using cerebrospinal fluid analysis or Aβ-PET. This can also be used in clinical settings and clinical trials, aiding subject recruitment and evaluation of treatment efficacy. Imputation of the Aβ-positivity status could also complement Aβ-PET by identifying individuals who would benefit the most from this assessment.

Quantitative tumor segmentation for evaluation of extent of glioblastoma resection to facilitate multisite clinical trials

Standard-of-care therapy for glioblastomas, the most common and aggressive primary adult brain neoplasm, is maximal safe resection, followed by radiation and chemotherapy. Because maximizing resection may be beneficial for these patients, improving tumor extent of resection (EOR) with methods such as intraoperative 5-aminolevulinic acid fluorescence-guided surgery (FGS) is currently under evaluation. However, it is difficult to reproducibly judge EOR in these studies due to the lack of reliable tumor segmentation methods, especially for postoperative magnetic resonance imaging (MRI) scans. Therefore, a reliable, easily distributable segmentation method is needed to permit valid comparison, especially across multiple sites. We report a segmentation method that combines versatile region-of-interest blob generation with automated clustering methods. We applied this to glioblastoma cases undergoing FGS and matched controls to illustrate the method's reliability and accuracy. Agreement and interrater variability between segmentations were assessed using the concordance correlation coefficient, and spatial accuracy was determined using the Dice similarity index and mean Euclidean distance. Fuzzy C-means clustering with three classes was the best performing method, generating volumes with high agreement with manual contouring and high interrater agreement preoperatively and postoperatively. The proposed segmentation method allows tumor volume measurements of contrast-enhanced T 1-weighted images in the unbiased, reproducible fashion necessary for quantifying EOR in multicenter trials.

Imaging-based biomarkers of cognitive performance in older adults constructed via high-dimensional pattern regression applied to MRI and PET

In this study, we used high-dimensional pattern regression methods based on structural (gray and white matter; GM and WM) and functional (positron emission tomography of regional cerebral blood flow; PET) brain data to identify cross-sectional imaging biomarkers of cognitive performance in cognitively normal older adults from the Baltimore Longitudinal Study of Aging (BLSA). We focused on specific components of executive and memory domains known to decline with aging, including manipulation, semantic retrieval, long-term memory (LTM), and short-term memory (STM). For each imaging modality, brain regions associated with each cognitive domain were generated by adaptive regional clustering. A relevance vector machine was adopted to model the nonlinear continuous relationship between brain regions and cognitive performance, with cross-validation to select the most informative brain regions (using recursive feature elimination) as imaging biomarkers and optimize model parameters. Predicted cognitive scores using our regression algorithm based on the resulting brain regions correlated well with actual performance. Also, regression models obtained using combined GM, WM, and PET imaging modalities outperformed models based on single modalities. Imaging biomarkers related to memory performance included the orbito-frontal and medial temporal cortical regions with LTM showing stronger correlation with the temporal lobe than STM. Brain regions predicting executive performance included orbito-frontal, and occipito-temporal areas. The PET modality had higher contribution to most cognitive domains except manipulation, which had higher WM contribution from the superior longitudinal fasciculus and the genu of the corpus callosum. These findings based on machine-learning methods demonstrate the importance of combining structural and functional imaging data in understanding complex cognitive mechanisms and also their potential usage as biomarkers that predict cognitive status.

Mapping of the internal structure of human habenula with ex vivo MRI at 7T

The habenula is a small but important nucleus located next to the third ventricle in front of the pineal body. It helps to control the human reward system and is considered to play a key role in emotion, showing increased activation in major depressive disorders. Its dysfunction may underlie several neurological and psychiatric disorders. It is now possible to visualize the habenula and its anatomical subdivisions—medial habenula (MHB) and lateral habenula (LHB)—using MR techniques. The aim of this study was to further differentiate substructures within human lateral habenula (LHB) using ex vivo ultra-high field MR structural imaging, distinguishing between a medial part (m-LHB) and a lateral part (l-LHB). High resolution T1w images with 0.3-mm isotropic resolution and T2^*w images with 60-micrometer isotropic resolution were acquired on a 7T MR scanner and quantitative maps of T1 and T2^* were calculated. Cluster analysis of image intensity was performed using the Fuzzy and Noise Tolerant Adaptive Segmentation Method (FANTASM) tool. Ultra-high resolution structural MRI of ex vivo brain tissue at 7T provided sufficient SNR and contrast to discriminate the medial and lateral habenular nuclei. Heterogeneity was observed in the lateral habenula (LHB) nuclei, with clear distinctions between lateral and medial parts (m-LHB, l-LHB) and with the neighboring medial habenula (MHB). Clustering analysis based on the T1 and T2^* maps strongly showed 4–6 clusters as subcomponents of lateral and medial habenula.