Adaptive segmentation of MRI data

Finite element modeling of human head from MRI

The FEM geometry modeling of realistic head is a key issue for the research on FEM-based EEG/MEG. In this paper, a methodology is developed to construct this kind of model. By using this method, a five-layer realistic head FEM model is obtained, and with its application in FEM-based EEG, a satisfying result shows the reliability of the model.

AdaPT: An adaptive preterm segmentation algorithm for neonatal brain MRI

Advances in neonatal care have improved the survival of infants born prematurely although these infants remain at increased risk of adverse neurodevelopmental outcome. The measurement of white matter structure and features of the cortical surface can help define biomarkers that predict this risk. The measurement of these structures relies upon accurate automated segmentation routines, but these are often confounded by neonatal-specific imaging difficulties including poor contrast, low resolution, partial volume effects and the presence of significant natural and pathological anatomical variability. In this work we develop and evaluate an adaptive preterm multi-modal maximum a posteriori expectation-maximisation segmentation algorithm (AdaPT) incorporating an iterative relaxation strategy that adapts the tissue proportion priors toward the subject data. Also incorporated are intensity non-uniformity correction, a spatial homogeneity term in the form of a Markov random field and furthermore, the proposed method explicitly models the partial volume effect specifically mitigating the neonatal specific grey and white matter contrast inversion. Spatial priors are iteratively relaxed, enabling the segmentation of images with high anatomical disparity from a normal population. Experiments performed on a clinical cohort of 92 infants are validated against manual segmentation of normal and pathological cortical grey matter, cerebellum and ventricular volumes. Dice overlap scores increase significantly when compared to a widely-used maximum likelihood expectation maximisation algorithm for pathological cortical grey matter, cerebellum and ventricular volumes. Adaptive maximum a posteriori expectation maximisation is shown to be a useful tool for accurate and robust neonatal brain segmentation.

3-T MR-guided brachytherapy for gynecologic malignancies

Gynecologic malignancies are a leading cause of death in women worldwide. Standard treatment for many primary and recurrent gynecologic cancer cases includes external-beam radiation followed by brachytherapy. Magnetic resonance (MR) imaging is beneficial in diagnostic evaluation, in mapping the tumor location to tailor radiation dose and in monitoring the tumor response to treatment. Initial studies of MR guidance in gynecologic brachytherapy demonstrate the ability to optimize tumor coverage and reduce radiation dose to normal tissues, resulting in improved outcomes for patients. In this article, we describe a methodology to aid applicator placement and treatment planning for 3 Tesla (3-T) MR-guided brachytherapy that was developed specifically for gynecologic cancers. This methodology has been used in 18 cases from September 2011 to May 2012 in the Advanced Multimodality Image Guided Operating (AMIGO) suite at Brigham and Women's Hospital. AMIGO comprises state-of-the-art tools for MR imaging, image analysis and treatment planning. An MR sequence using three-dimensional (3D)-balanced steady-state free precession in a 3-T MR scanner was identified as the best sequence for catheter identification with ballooning artifact at the tip. 3D treatment planning was performed using MR images. Items in development include software designed to support virtual needle trajectory planning that uses probabilistic bias correction, graph-based segmentation and image registration algorithms. The results demonstrate that 3-T MR image guidance has a role in gynecologic brachytherapy. These novel developments have the potential to improve targeted treatment to the tumor while sparing the normal tissues.

An open source multivariate framework for n-tissue segmentation with evaluation on public data

We introduce Atropos, an ITK-based multivariate n-class open source segmentation algorithm distributed with ANTs ( http://www.picsl.upenn.edu/ANTs). The Bayesian formulation of the segmentation problem is solved using the Expectation Maximization (EM) algorithm with the modeling of the class intensities based on either parametric or non-parametric finite mixtures. Atropos is capable of incorporating spatial prior probability maps (sparse), prior label maps and/or Markov Random Field (MRF) modeling. Atropos has also been efficiently implemented to handle large quantities of possible labelings (in the experimental section, we use up to 69 classes) with a minimal memory footprint. This work describes the technical and implementation aspects of Atropos and evaluates its performance on two different ground-truth datasets. First, we use the BrainWeb dataset from Montreal Neurological Institute to evaluate three-tissue segmentation performance via (1) K-means segmentation without use of template data; (2) MRF segmentation with initialization by prior probability maps derived from a group template; (3) Prior-based segmentation with use of spatial prior probability maps derived from a group template. We also evaluate Atropos performance by using spatial priors to drive a 69-class EM segmentation problem derived from the Hammers atlas from University College London. These evaluation studies, combined with illustrative examples that exercise Atropos options, demonstrate both performance and wide applicability of this new platform-independent open source segmentation tool.

Reconstruction of fetal brain MRI with intensity matching and complete outlier removal

We propose a method for the reconstruction of volumetric fetal MRI from 2D slices, comprising super-resolution reconstruction of the volume interleaved with slice-to-volume registration to correct for the motion. The method incorporates novel intensity matching of acquired 2D slices and robust statistics which completely excludes identified misregistered or corrupted voxels and slices. The reconstruction method is applied to motion-corrupted data simulated from MRI of a preterm neonate, as well as 10 clinically acquired thick-slice fetal MRI scans and three scan-sequence optimized thin-slice fetal datasets. The proposed method produced high quality reconstruction results from all the datasets to which it was applied. Quantitative analysis performed on simulated and clinical data shows that both intensity matching and robust statistics result in statistically significant improvement of super-resolution reconstruction. The proposed novel EM-based robust statistics also improves the reconstruction when compared to previously proposed Huber robust statistics. The best results are obtained when thin-slice data and the correct approximation of the point spread function is used. This paper addresses the need for a comprehensive reconstruction algorithm of 3D fetal MRI, so far lacking in the scientific literature.

3D Slicer as an image computing platform for the Quantitative Imaging Network

Quantitative analysis has tremendous but mostly unrealized potential in healthcare to support objective and accurate interpretation of the clinical imaging. In 2008, the National Cancer Institute began building the Quantitative Imaging Network (QIN) initiative with the goal of advancing quantitative imaging in the context of personalized therapy and evaluation of treatment response. Computerized analysis is an important component contributing to reproducibility and efficiency of the quantitative imaging techniques. The success of quantitative imaging is contingent on robust analysis methods and software tools to bring these methods from bench to bedside. 3D Slicer is a free open-source software application for medical image computing. As a clinical research tool, 3D Slicer is similar to a radiology workstation that supports versatile visualizations but also provides advanced functionality such as automated segmentation and registration for a variety of application domains. Unlike a typical radiology workstation, 3D Slicer is free and is not tied to specific hardware. As a programming platform, 3D Slicer facilitates translation and evaluation of the new quantitative methods by allowing the biomedical researcher to focus on the implementation of the algorithm and providing abstractions for the common tasks of data communication, visualization and user interface development. Compared to other tools that provide aspects of this functionality, 3D Slicer is fully open source and can be readily extended and redistributed. In addition, 3D Slicer is designed to facilitate the development of new functionality in the form of 3D Slicer extensions. In this paper, we present an overview of 3D Slicer as a platform for prototyping, development and evaluation of image analysis tools for clinical research applications. To illustrate the utility of the platform in the scope of QIN, we discuss several use cases of 3D Slicer by the existing QIN teams, and we elaborate on the future directions that can further facilitate development and validation of imaging biomarkers using 3D Slicer.

Unsupervised Neural Techniques Applied to MR Brain Image Segmentation

The primary goal of brain image segmentation is to partition a given brain image into different regions representing anatomical structures. Magnetic resonance image (MRI) segmentation is especially interesting, since accurate segmentation in white matter, grey matter and cerebrospinal fluid provides a way to identify many brain disorders such as dementia, schizophrenia or Alzheimer’s disease (AD). Then, image segmentation results in a very interesting tool for neuroanatomical analyses. In this paper we show three alternatives to MR brain image segmentation algorithms, with the Self-Organizing Map (SOM) as the core of the algorithms. The procedures devised do not use any a priori knowledge about voxel class assignment, and results in fully-unsupervised methods for MRI segmentation, making it possible to automatically discover different tissue classes. Our algorithm has been tested using the images from the Internet Brain Image Repository (IBSR) outperforming existing methods, providing values for the average overlap metric of 0.7 for the white and grey matter and 0.45 for the cerebrospinal fluid. Furthermore, it also provides good results for high-resolution MR images provided by the Nuclear Medicine Service of the “Virgen de las Nieves” Hospital (Granada, Spain).

Simultaneous segmentation and inhomogeneity correction in magnetic resonance images

In Magnetic Resonance Imaging (MRI), intensity inhomogeneity has been an issue affecting the quality of post processing. In this paper, we present a simultaneous segmentation and inhomogeneity correction (IC) method based on active contour algorithm. It uses a generative model which is a modified Mumford-Shah functional proposed by Chan and Vese. The piecewise constant image model in the functional is multiplied by an underlying intensity inhomogeneity field. The inhomogeneity field and piecewise constant function are jointly estimated in an iterative way including solving the associated contour evolution equation and updating corresponding parameters. The algorithm is implemented using the level set framework. Test on MRI leg data shows our method achieves more accurate segmentation and IC results than other related methods in MR images with strong intensity inhomogeneity.

Ray-based approach to skeletal muscle perfusion measurement on interventional x-ray systems

PURPOSE

Periprocedural assessment of tissue perfusion by imaging methods could improve outcome control during treatment of peripheral vascular disease. Currently, endovascular revascularization treatments are assessed by planar angiography which only allows for qualitative inspection of blood flow in vessels. In this paper, we present a method for periprocedural perfusion estimation based on temporal attenuation curves in skeletal muscles using angiographic C-arm systems.

METHODS

The proposed method tackles the loss of spatial depth information which occurs in conventional angiography by combining the acquired angiograms with two additional C-arm rotational soft tissue scans. The area subject to contrast propagation is segmented from the two images that are tomographically reconstructed from the rotational scans and is then used to estimate the spatially averaged temporal contrast attenuation along the x-ray directions from the angiograms. A segmentation method which is tailored to the estimation procedure is applied to limit inaccuracies in the estimation. The accuracy of the method in estimating tissue blood flow in muscular tissue is evaluated in a simulation study using experimental data from CT perfusion acquisitions.

RESULTS

Results show that perfusion estimation accuracy is limited owing to spatial inhomogeneity of contrast in muscular tissue and to the presence of vessels along the x-ray directions. Nonetheless, the spatially averaged perfusion quantification allows for improved visual differentiation of normal and hypoperfused tissue in comparison with conventional digital subtraction angiography.

CONCLUSIONS

Periprocedural assessment of muscle perfusion through digital subtraction angiography is challenging due to lack of longitudinal information in the planar projections. By including additional 3D information on the anatomy retrieved from rotational soft tissue scans, the visualization and differentiation of normal and hypoperfused areas can be improved.

Simultaneous Segmentation and Statistical Label Fusion

Labeling or segmentation of structures of interest in medical imaging plays an essential role in both clinical and scientific understanding. Two of the common techniques to obtain these labels are through either fully automated segmentation or through multi-atlas based segmentation and label fusion. Fully automated techniques often result in highly accurate segmentations but lack the robustness to be viable in many cases. On the other hand, label fusion techniques are often extremely robust, but lack the accuracy of automated algorithms for specific classes of problems. Herein, we propose to perform simultaneous automated segmentation and statistical label fusion through the reformulation of a generative model to include a linkage structure that explicitly estimates the complex global relationships between labels and intensities. These relationships are inferred from the atlas labels and intensities and applied to the target using a non-parametric approach. The novelty of this approach lies in the combination of previously exclusive techniques and attempts to combine the accuracy benefits of automated segmentation with the robustness of a multi-atlas based approach. The accuracy benefits of this simultaneous approach are assessed using a multi-label multi- atlas whole-brain segmentation experiment and the segmentation of the highly variable thyroid on computed tomography images. The results demonstrate that this technique has major benefits for certain types of problems and has the potential to provide a paradigm shift in which the lines between statistical label fusion and automated segmentation are dramatically blurred.

Fuzzy local Gaussian mixture model for brain MR image segmentation

Accurate brain tissue segmentation from magnetic resonance (MR) images is an essential step in quantitative brain image analysis. However, due to the existence of noise and intensity inhomogeneity in brain MR images, many segmentation algorithms suffer from limited accuracy. In this paper, we assume that the local image data within each voxel's neighborhood satisfy the Gaussian mixture model (GMM), and thus propose the fuzzy local GMM (FLGMM) algorithm for automated brain MR image segmentation. This algorithm estimates the segmentation result that maximizes the posterior probability by minimizing an objective energy function, in which a truncated Gaussian kernel function is used to impose the spatial constraint and fuzzy memberships are employed to balance the contribution of each GMM. We compared our algorithm to state-of-the-art segmentation approaches in both synthetic and clinical data. Our results show that the proposed algorithm can largely overcome the difficulties raised by noise, low contrast, and bias field, and substantially improve the accuracy of brain MR image segmentation.

Machine Learning for Brain Image Segmentation

In this chapter, the authors review a variety of algorithms developed by different groups for automatically segmenting structures in medical images, such as brain MRI scans. Some of the simpler methods, based on active contours, deformable image registration, and anisotropic Markov random fields, have known weaknesses, which can be largely overcome by learning methods that better encode knowledge on anatomical variability. The authors show how the anatomical segmentation problem may be re-cast in a Bayesian framework. They then present several different learning techniques increasing in complexity until they derive two algorithms recently proposed by the authors. The authors show how these automated algorithms are validated empirically, by comparison with segmentations by experts, which serve as independent ground truth, and in terms of their power to detect disease effects in Alzheimer’s disease. They show how these methods can be used to investigate factors that influence disease progression in databases of thousands of images. Finally the authors indicate some promising directions for future work.

Magnetic Resonance Image Analysis for Brain CAD Systems with Machine Learning

This chapter describes the image analysis for brain Computer-Aided Diagnosis (CAD) systems with machine learning techniques, which could assist radiologists in the detection of such brain diseases as asymptomatic unruptured aneurysms, Alzheimer’s Disease (AD), vascular dementia, and Multiple Sclerosis (MS) by magnetic resonance imaging. Image analysis in CAD systems consists of image enhancement, initial detection, and image feature extraction, including segmentation. In addition, the authors review the classification of true and false positives using machine learning techniques, as well as the evaluation methods and development cycle for CAD systems.

A fast and fully automatic method for cerebrovascular segmentation on time-of-flight (TOF) MRA image

The precise three-dimensional (3-D) segmentation of cerebral vessels from magnetic resonance angiography (MRA) images is essential for the detection of cerebrovascular diseases (e.g., occlusion, aneurysm). The complex 3-D structure of cerebral vessels and the low contrast of thin vessels in MRA images make precise segmentation difficult. We present a fast, fully automatic segmentation algorithm based on statistical model analysis and improved curve evolution for extracting the 3-D cerebral vessels from a time-of-flight (TOF) MRA dataset. Cerebral vessels and other tissue (brain tissue, CSF, and bone) in TOF MRA dataset are modeled by Gaussian distribution and combination of Rayleigh with several Gaussian distributions separately. The region distribution combined with gradient information is used in edge-strength of curve evolution as one novel mode. This edge-strength function is able to determine the boundary of thin vessels with low contrast around brain tissue accurately and robustly. Moreover, a fast level set method is developed to implement the curve evolution to assure high efficiency of the cerebrovascular segmentation. Quantitative comparisons with 10 sets of manual segmentation results showed that the average volume sensitivity, the average branch sensitivity, and average mean absolute distance error are 93.6%, 95.98%, and 0.333 mm, respectively. By applying the algorithm to 200 clinical datasets from three hospitals, it is demonstrated that the proposed algorithm can provide good quality segmentation capable of extracting a vessel with a one-voxel diameter in less than 2 min. Its accuracy and speed make this novel algorithm more suitable for a clinical computer-aided diagnosis system.

Longitudinally guided level sets for consistent tissue segmentation of neonates

Quantification of brain development as well as disease-induced pathologies in neonates often requires precise delineation of white matter, grey matter and cerebrospinal fluid. Unlike adults, tissue segmentation in neonates is significantly more challenging due to the inherently lower tissue contrast. Most existing methods take a voxel-based approach and are limited to working with images from a single time-point, even though longitudinal scans are available. We take a different approach by taking advantage of the fact that the pattern of the major sulci and gyri are already present in the neonates and generally preserved but fine-tuned during brain development. That is, the segmentation of late-time-point image can be used to guide the segmentation of neonatal image. Accordingly, we propose a novel longitudinally guided level-sets method for consistent neonatal image segmentation by combining local intensity information, atlas spatial prior, cortical thickness constraint, and longitudinal information into a variational framework. The minimization of the proposed energy functional is strictly derived from a variational principle. Validation performed on both simulated and in vivo neonatal brain images shows promising results.

A NOVEL FUZZY C-MEANS ALGORITHM AND ITS APPLICATION

The accurate and effective algorithm for segmenting image is very useful in many fields, especially in medical images. In this paper we introduced a novel method that focuses on segmenting the brain MR Image that is important for neural diseases. Because of many noises embedded in the acquiring procedure, such as eddy currents, susceptibility artifacts, rigid body motion and intensity inhomogeneity, segmenting the brain MR image is a difficult work. In this algorithm, we overcame the inhomogeneity shortage, by modifying the objective function by compensating its immediate neighborhood effect using Gaussian smooth method for decreasing the influence of the inhomogeneity and increasing the segmenting accuracy. Using simulate image and clinical MRI data, experiments show that our proposed algorithm is effective.

NON-RIGID REGISTRATION AND CORRESPONDENCE FINDING IN MEDICAL IMAGE ANALYSIS USING MULTIPLE-LAYER FLEXIBLE MESH TEMPLATE MATCHING

In this paper we present a novel technique for non-rigid medical image registration and correspondence finding based on a multiple-layer flexible mesh template matching technique. A statistical anatomical model is built in the form of a tetrahedral mesh, which incorporates both shape and density properties of the anatomical structure. After the affine transformation and global deformation of the model are computed by optimizing an energy function, a multiple-layer flexible mesh template matching is applied to find the vertex correspondence and achieve local deformation. The multiple-layer structure of the template can be used to describe different scale of anatomical features; furthermore, the template matching is flexible which makes the correspondence finding robust. A leave-one-out validation has been conducted to demonstrate the effectiveness and accuracy of our method.

HYBRID APPROACH FOR BRAIN TUMOR SEGMENTATION IN MAGNETIC RESONANCE IMAGES USING CELLULAR NEURAL NETWORKS AND OPTIMIZATION TECHNIQUES

Tumor segmentation from brain magnetic resonance image data is an important but time consuming task performed manually by medical experts. Automating this process is challenging due to the high diversity in appearance of tumor tissue among different patients and in many cases, similarity between tumor and normal tissue. This paper deals with an efficient segmentation algorithm for extracting brain tumors in magnetic resonance images using Cellular Neural Networks (CNN). Learning CNN templates values are formulated as an optimization problem. The template coefficients (weights) of an CNN which will give a desired performance, can be derived by learning genetic algorithm and simulated annealing optimization techniques. The objective of this work is to compare the performance of genetic algorithm (GA) and simulated annealing (SA) for finding the optimum template values in the CNN which is used for segmenting the tumor region in the abnormal MR images. The method is applied on real data of MRI images of thirty patients with four different types of tumors. The results are compared with radiologist labeled ground truth. Quantitative analysis between ground truth and segmented tumor is presented in terms of segmentation efficiency. From the analysis and performance measures like segmentation accuracy, it is inferred that the brain tumor segmentation is best done using CNN with genetic algorithm template optimization than CNN with simulated annealing template optimization. An average accuracy rate of above 95% was obtained using this segmentation algorithm.

Adaptive neonate brain segmentation

Babies born prematurely are at increased risk of adverse neurodevelopmental outcomes. Recent advances suggest that measurement of brain volumes can help in defining biomarkers for neurodevelopmental outcome. These techniques rely on an accurate segmentation of the MRI data. However, due to lack of contrast, partial volume (PV) effect, the existence of both hypo- and hyper-intensities and significant natural and pathological anatomical variability, the segmentation of neonatal brain MRI is challenging. We propose a pipeline for image segmentation that uses a novel multi-model Maximum a posteriori Expectation Maximisation (MAP-EM) segmentation algorithm with a prior over both intensities and the tissue proportions, a B0 inhomogeneity correction, and a spatial homogeneity term through the use of a Markov Random Field. This robust and adaptive technique enables the segmentation of images with high anatomical disparity from a normal population. Furthermore, the proposed method implicitly models Partial Volume, mitigating the problem of neonatal white/grey matter intensity inversion. Experiments performed on a clinical cohort show expected statistically significant correlations with gestational age at birth and birthweight. Furthermore, the proposed method obtains statistically significant improvements in Dice scores when compared to the a Maximum Likelihood EM algorithm.

Learning likelihoods for labeling (L3): a general multi-classifier segmentation algorithm

PURPOSE

To develop an MRI segmentation method for brain tissues, regions, and substructures that yields improved classification accuracy. Current brain segmentation strategies include two complementary strategies. Multi-spectral classification techniques generate excellent segmentations for tissues with clear intensity contrast, but fail to identify structures defined largely by location, such as lobar parcellations and certain subcortical structures. Conversely, multi-template label fusion methods are excellent for structures defined largely by location, but perform poorly when segmenting structures that cannot be accurately identified through a consensus of registered templates.

METHODS

We propose here a novel multi-classifier fusion algorithm with the advantages of both types of segmentation strategy. We illustrate and validate this algorithm using a group of 14 expertly hand-labeled images.

RESULTS

Our method generated segmentations of cortical and subcortical structures that were more similar to hand-drawn segmentations than majority vote label fusion or a recently published intensity/label fusion method.

CONCLUSIONS

We have presented a novel, general segmentation algorithm with the advantages of both statistical classifiers and label fusion techniques.

A joint shape evolution approach to medical image segmentation using expectation-maximization algorithm

This study proposes an expectation-maximization (EM)-based curve evolution algorithm for segmentation of magnetic resonance brain images. In the proposed algorithm, the evolution curve is constrained not only by a shape-based statistical model but also by a hidden variable model from image observation. The hidden variable model herein is defined by the local voxel labeling, which is unknown and estimated by the expected likelihood function derived from the image data and prior anatomical knowledge. In the M-step, the shapes of the structures are estimated jointly by encoding the hidden variable model and the statistical prior model obtained from the training stage. In the E-step, the expected observation likelihood and the prior distribution of the hidden variables are estimated. In experiments, the proposed automatic segmentation algorithm is applied to multiple gray nuclei structures such as caudate, putamens and thalamus of three-dimensional magnetic resonance imaging in volunteers and patients. As for the robustness and accuracy of the segmentation algorithm, the results of the proposed EM-joint shape-based algorithm outperformed those obtained using the statistical shape model-based techniques in the same framework and a current state-of-the-art region competition level set method.

Automatic segmentation of meningioma from non-contrasted brain MRI integrating fuzzy clustering and region growing

Background

In recent years, magnetic resonance imaging (MRI) has become important in brain tumor diagnosis. Using this modality, physicians can locate specific pathologies by analyzing differences in tissue character presented in different types of MR images.

This paper uses an algorithm integrating fuzzy-c-mean (FCM) and region growing techniques for automated tumor image segmentation from patients with menigioma. Only non-contrasted T1 and T2 -weighted MR images are included in the analysis. The study's aims are to correctly locate tumors in the images, and to detect those situated in the midline position of the brain.

Methods

The study used non-contrasted T1- and T2-weighted MR images from 29 patients with menigioma. After FCM clustering, 32 groups of images from each patient group were put through the region-growing procedure for pixels aggregation. Later, using knowledge-based information, the system selected tumor-containing images from these groups and merged them into one tumor image. An alternative semi-supervised method was added at this stage for comparison with the automatic method. Finally, the tumor image was optimized by a morphology operator. Results from automatic segmentation were compared to the "ground truth" (GT) on a pixel level. Overall data were then evaluated using a quantified system.

Results

The quantified parameters, including the "percent match" (PM) and "correlation ratio" (CR), suggested a high match between GT and the present study's system, as well as a fair level of correspondence. The results were compatible with those from other related studies. The system successfully detected all of the tumors situated at the midline of brain.

Six cases failed in the automatic group. One also failed in the semi-supervised alternative. The remaining five cases presented noticeable edema inside the brain. In the 23 successful cases, the PM and CR values in the two groups were highly related.

Conclusions

Results indicated that, even when using only two sets of non-contrasted MR images, the system is a reliable and efficient method of brain-tumor detection. With further development the system demonstrates high potential for practical clinical use.

Mapping fetal brain development in utero using magnetic resonance imaging: the Big Bang of brain mapping

The development of tools to construct and investigate probabilistic maps of the adult human brain from magnetic resonance imaging (MRI) has led to advances in both basic neuroscience and clinical diagnosis. These tools are increasingly being applied to brain development in adolescence and childhood, and even to neonatal and premature neonatal imaging. Even earlier in development, parallel advances in clinical fetal MRI have led to its growing use as a tool in challenging medical conditions. This has motivated new engineering developments encompassing optimal fast MRI scans and techniques derived from computer vision, the combination of which allows full 3D imaging of the moving fetal brain in utero without sedation. These promise to provide a new and unprecedented window into early human brain growth. This article reviews the developments that have led us to this point, examines the current state of the art in the fields of fast fetal imaging and motion correction, and describes the tools to analyze dynamically changing fetal brain structure. New methods to deal with developmental tissue segmentation and the construction of spatiotemporal atlases are examined, together with techniques to map fetal brain growth patterns.

Trimmed-likelihood estimation for focal lesions and tissue segmentation in multisequence MRI for multiple sclerosis

We present a new automatic method for segmentation of multiple sclerosis (MS) lesions in magnetic resonance images. The method performs tissue classification using a model of intensities of the normal appearing brain tissues. In order to estimate the model, a trimmed likelihood estimator is initialized with a hierarchical random approach in order to be robust to MS lesions and other outliers present in real images. The algorithm is first evaluated with simulated images to assess the importance of the robust estimator in presence of outliers. The method is then validated using clinical data in which MS lesions were delineated manually by several experts. Our method obtains an average Dice similarity coefficient (DSC) of 0.65, which is close to the average DSC obtained by raters (0.66).

Segmentation of multiple sclerosis lesions in MR images: a review

INTRODUCTION

Multiple sclerosis (MS) is an inflammatory demyelinating disease that the parts of the nervous system through the lesions generated in the white matter of the brain. It brings about disabilities in different organs of the body such as eyes and muscles. Early detection of MS and estimation of its progression are critical for optimal treatment of the disease.

METHODS

For diagnosis and treatment evaluation of MS lesions, they may be detected and segmented in Magnetic Resonance Imaging (MRI) scans of the brain. However, due to the large amount of MRI data to be analyzed, manual segmentation of the lesions by clinical experts translates into a very cumbersome and time consuming task. In addition, manual segmentation is subjective and prone to human errors. Several groups have developed computerized methods to detect and segment MS lesions. These methods are not categorized and compared in the past.

RESULTS

This paper reviews and compares various MS lesion segmentation methods proposed in recent years. It covers conventional methods like multilevel thresholding and region growing, as well as more recent Bayesian methods that require parameter estimation algorithms. It also covers parameter estimation methods like expectation maximization and adaptive mixture model which are among unsupervised techniques as well as kNN and Parzen window methods that are among supervised techniques.

CONCLUSIONS

Integration of knowledge-based methods such as atlas-based approaches with Bayesian methods increases segmentation accuracy. In addition, employing intelligent classifiers like Fuzzy C-Means, Fuzzy Inference Systems, and Artificial Neural Networks reduces misclassified voxels.

Segmentation of multiple sclerosis lesions in MR images: a review

INTRODUCTION

Multiple sclerosis (MS) is an inflammatory demyelinating disease that the parts of the nervous system through the lesions generated in the white matter of the brain. It brings about disabilities in different organs of the body such as eyes and muscles. Early detection of MS and estimation of its progression are critical for optimal treatment of the disease.

METHODS

For diagnosis and treatment evaluation of MS lesions, they may be detected and segmented in Magnetic Resonance Imaging (MRI) scans of the brain. However, due to the large amount of MRI data to be analyzed, manual segmentation of the lesions by clinical experts translates into a very cumbersome and time consuming task. In addition, manual segmentation is subjective and prone to human errors. Several groups have developed computerized methods to detect and segment MS lesions. These methods are not categorized and compared in the past.

RESULTS

This paper reviews and compares various MS lesion segmentation methods proposed in recent years. It covers conventional methods like multilevel thresholding and region growing, as well as more recent Bayesian methods that require parameter estimation algorithms. It also covers parameter estimation methods like expectation maximization and adaptive mixture model which are among unsupervised techniques as well as kNN and Parzen window methods that are among supervised techniques.

CONCLUSIONS

Integration of knowledge-based methods such as atlas-based approaches with Bayesian methods increases segmentation accuracy. In addition, employing intelligent classifiers like Fuzzy C-Means, Fuzzy Inference Systems, and Artificial Neural Networks reduces misclassified voxels.

On the convergence of EM-like algorithms for image segmentation using Markov random fields

Inference of Markov random field images segmentation models is usually performed using iterative methods which adapt the well-known expectation-maximization (EM) algorithm for independent mixture models. However, some of these adaptations are ad hoc and may turn out numerically unstable. In this paper, we review three EM-like variants for Markov random field segmentation and compare their convergence properties both at the theoretical and practical levels. We specifically advocate a numerical scheme involving asynchronous voxel updating, for which general convergence results can be established. Our experiments on brain tissue classification in magnetic resonance images provide evidence that this algorithm may achieve significantly faster convergence than its competitors while yielding at least as good segmentation results.

A discriminative model-constrained EM approach to 3D MRI brain tissue classification and intensity non-uniformity correction

We describe a fully automated method for tissue classification, which is the segmentation into cerebral gray matter (GM), cerebral white matter (WM), and cerebral spinal fluid (CSF), and intensity non-uniformity (INU) correction in brain magnetic resonance imaging (MRI) volumes. It combines supervised MRI modality-specific discriminative modeling and unsupervised statistical expectation maximization (EM) segmentation into an integrated Bayesian framework. While both the parametric observation models and the non-parametrically modeled INUs are estimated via EM during segmentation itself, a Markov random field (MRF) prior model regularizes segmentation and parameter estimation. Firstly, the regularization takes into account knowledge about spatial and appearance-related homogeneity of segments in terms of pairwise clique potentials of adjacent voxels. Secondly and more importantly, patient-specific knowledge about the global spatial distribution of brain tissue is incorporated into the segmentation process via unary clique potentials. They are based on a strong discriminative model provided by a probabilistic boosting tree (PBT) for classifying image voxels. It relies on the surrounding context and alignment-based features derived from a probabilistic anatomical atlas. The context considered is encoded by 3D Haar-like features of reduced INU sensitivity. Alignment is carried out fully automatically by means of an affine registration algorithm minimizing cross-correlation. Both types of features do not immediately use the observed intensities provided by the MRI modality but instead rely on specifically transformed features, which are less sensitive to MRI artifacts. Detailed quantitative evaluations on standard phantom scans and standard real-world data show the accuracy and robustness of the proposed method. They also demonstrate relative superiority in comparison to other state-of-the-art approaches to this kind of computational task: our method achieves average Dice coefficients of 0.93 ± 0.03 (WM) and 0.90 ± 0.05 (GM) on simulated mono-spectral and 0.94 ± 0.02 (WM) and 0.92 ± 0.04 (GM) on simulated multi-spectral data from the BrainWeb repository. The scores are 0.81 ± 0.09 (WM) and 0.82 ± 0.06 (GM) and 0.87 ± 0.05 (WM) and 0.83 ± 0.12 (GM) for the two collections of real-world data sets-consisting of 20 and 18 volumes, respectively-provided by the Internet Brain Segmentation Repository.

Subregional 6-[18F]fluoro-ʟ-m-tyrosine uptake in the striatum in Parkinson's disease

BACKGROUND

In idiopathic Parkinson's disease (PD) the clinical features are heterogeneous and include different predominant symptoms. The aim of the present study was to determine the relationship between subregional aromatic l-amino acid decarboxylase (AADC) activity in the striatum and the cardinal motor symptoms of PD using high-resolution positron emission tomography (PET) with an AADC tracer, 6-[18F]fluoro-ʟ-m-tyrosine (FMT).

METHODS

We assessed 101 patients with PD and 19 healthy volunteers. PD was diagnosed based on the UK Brain Bank criteria by two experts on movement disorders. Motor symptoms were measured with the Unified Parkinson's Disease Rating Scale (UPDRS). FMT uptake in the subregions of the striatum was analyzed using semi-automated software for region-of-interest demarcation on co-registered magnetic resonance images.

RESULTS

In all PD patients, FMT uptake was decreased in the posterior putamen regardless of predominant motor symptoms and disease duration. Smaller uptake values were found in the putamen contralateral to the side with more affected limbs. The severity of bradykinesia, rigidity, and axial symptoms was correlated with the decrease of FMT uptake in the putamen, particularly in the anterior part. No significant correlation was observed between tremors and FMT uptake.

CONCLUSIONS

Decrease of FMT uptake in the posterior putamen appears to be most sensitive in mild PD and uptake in the anterior putamen may reflect the severity of main motor symptoms, except for tremor.