Fully automatic segmentation of the brain in MRI

Functional magnetic resonance imaging: comparison between activation maps and computation pipelines in a clinical context

In this study new evaluation strategies for comparing different Statistical Parametric Maps computed from fMRI time-series analysis software tools are proposed. The aim of our work is to assess and quantitatively evaluate the statistical agreement of activation maps. Some pre-processing steps are necessary to compare SPMs (Statistical Parametric Maps), including segmentation and co-registration. The study of the statistical agreement is carried out following two ways. The first way considers SPMs as the result of two classification processes and extracts confusion matrix and Cohen's kappa index to assess agreement. Some considerations will be made on the statistical dependence of classes and a new formulation of kappa index will be used for overcoming this problem. The second way considers SPMs as two 3D images, and computes the similarity of SPMs images with a fuzzy formulation of the Jaccard Index. Several experiments were conducted both to assess the performance of the proposed evaluation tools and to compare activation maps computation pipelines from two widely used software tools in a clinical context.

A 3-D visualization method for image-guided brain surgery

This paper deals with a 3D methodology for brain tumor image-guided surgery. The methodology is based on development of a visualization process that mimics the human surgeon behavior and decision-making. In particular, it originally constructs a 3D representation of a tumor by using the segmented version of the 2D MRI images. Then it develops an optimal path for the tumor extraction based on minimizing the surgical effort and penetration area. A cost function, incorporated in this process, minimizes the damage surrounding healthy tissues taking into consideration the constraints of a new snake-like surgical tool proposed here. The tumor extraction method presented in this paper is compared with the ordinary method used on brain surgery, which is based on a straight-line based surgical tool. Illustrative examples based on real simulations present the advantages of the 3D methodology proposed here.

Patient specific prostate segmentation in 3-d magnetic resonance images

Accurate localization of the prostate and its surrounding tissue is essential in the treatment of prostate cancer. This paper presents a novel approach to fully automatically segment the prostate, including its seminal vesicles, within a few minutes of a magnetic resonance (MR) scan acquired without an endorectal coil. Such MR images are important in external beam radiation therapy, where using an endorectal coil is highly undesirable. The segmentation is obtained using a deformable model that is trained on-the-fly so that it is specific to the patient's scan. This case specific deformable model consists of a patient specific initialized triangulated surface and image feature model that are trained during its initialization. The image feature model is used to deform the initialized surface by template matching image features (via normalized cross-correlation) to the features of the scan. The resulting deformations are regularized over the surface via well established simple surface smoothing algorithms, which is then made anatomically valid via an optimized shape model. Mean and median Dice's similarity coefficients (DSCs) of 0.85 and 0.87 were achieved when segmenting 3T MR clinical scans of 50 patients. The median DSC result was equal to the inter-rater DSC and had a mean absolute surface error of 1.85 mm. The approach is showed to perform well near the apex and seminal vesicles of the prostate.

Cerebellar grey-matter deficits, cannabis use and first-episode schizophrenia in adolescents and young adults

Epidemiological data link adolescent cannabis use to psychosis and schizophrenia, but its contribution to schizophrenia neuropathology remains controversial. First-episode schizophrenia (FES) patients show regional cerebral grey- and white-matter changes as well as a distinct pattern of regional grey-matter loss in the vermis of the cerebellum. The cerebellum possesses a high density of cannabinoid type 1 receptors involved in the neuronal diversification of the developing brain. Cannabis abuse may interfere with this process during adolescent brain maturation leading to 'schizophrenia-like' cerebellar pathology. Magnetic resonance imaging and cortical pattern matching techniques were used to investigate cerebellar grey and white matter in FES patients with and without a history of cannabis use and non-psychiatric cannabis users. In the latter group we found lifetime dose-dependent regional reduction of grey matter in the right cerebellar lobules and a tendency for more profound grey-matter reduction in lobule III with younger age at onset of cannabis use. The overall regional grey-matter differences in cannabis users were within the normal variability of grey-matter distribution. By contrast, FES subjects had lower total cerebellar grey-matter:total cerebellar volume ratio and marked grey-matter loss in the vermis, pedunculi, flocculi and lobules compared to pair-wise matched healthy control subjects. This pattern and degree of grey-matter loss did not differ from age-matched FES subjects with comorbid cannabis use. Our findings indicate small dose-dependent effects of juvenile cannabis use on cerebellar neuropathology but no evidence of an additional effect of cannabis use on FES cerebellar grey-matter pathology.

Hierarchical statistical shape models of multiobject anatomical structures: application to brain MRI

The accurate segmentation of subcortical brain structures in magnetic resonance (MR) images is of crucial importance in the interdisciplinary field of medical imaging. Although statistical approaches such as active shape models (ASMs) have proven to be particularly useful in the modeling of multiobject shapes, they are inefficient when facing challenging problems. Based on the wavelet transform, the fully generic multiresolution framework presented in this paper allows us to decompose the interobject relationships into different levels of detail. The aim of this hierarchical decomposition is twofold: to efficiently characterize the relationships between objects and their particular localities. Experiments performed on an eight-object structure defined in axial cross sectional MR brain images show that the new hierarchical segmentation significantly improves the accuracy of the segmentation, and while it exhibits a remarkable robustness with respect to the size of the training set.

An accurate skull stripping method based on simplex meshes and histogram analysis for magnetic resonance images

Skull stripping methods are designed to eliminate the non-brain tissue in magnetic resonance (MR) brain images. Removal of non-brain tissues is a fundamental step in enabling the processing of brain MR images. The aim of this study is to develop an automatic accurate skull stripping method based on deformable models and histogram analysis. A rough-segmentation step is used to find the optimal starting point for the deformation and is based on thresholds and morphological operators. Thresholds are computed using comparisons with an atlas, and modeling by Gaussians. The deformable model is based on a simplex mesh and its deformation is controlled by the image local gray levels and the information obtained on the gray level modeling of the rough-segmentation. Our Simplex Mesh and Histogram Analysis Skull Stripping (SMHASS) method was tested on the following international databases commonly used in scientific articles: BrainWeb, Internet Brain Segmentation Repository (IBSR), and Segmentation Validation Engine (SVE). A comparison was performed against three of the best skull stripping methods previously published: Brain Extraction Tool (BET), Brain Surface Extractor (BSE), and Hybrid Watershed Algorithm (HWA). Performance was measured using the Jaccard index (J) and Dice coefficient (κ). Our method showed the best performance and differences were statistically significant (p<0.05): J=0.904 and κ=0.950 on BrainWeb; J=0.905 and κ=0.950 on IBSR; J=0.946 and κ=0.972 on SVE.

A review of atlas-based segmentation for magnetic resonance brain images

Normal and abnormal brains can be segmented by registering the target image with an atlas. Here, an atlas is defined as the combination of an intensity image (template) and its segmented image (the atlas labels). After registering the atlas template and the target image, the atlas labels are propagated to the target image. We define this process as atlas-based segmentation. In recent years, researchers have investigated registration algorithms to match atlases to query subjects and also strategies for atlas construction. In this paper we present a review of the automated approaches for atlas-based segmentation of magnetic resonance brain images. We aim to point out the strengths and weaknesses of atlas-based methods and suggest new research directions. We use two different criteria to present the methods. First, we refer to the algorithms according to their atlas-based strategy: label propagation, multi-atlas methods, and probabilistic techniques. Subsequently, we classify the methods according to their medical target: the brain and its internal structures, tissue segmentation in healthy subjects, tissue segmentation in fetus, neonates and elderly subjects, and segmentation of damaged brains. A quantitative comparison of the results reported in the literature is also presented.

Automatic brain extraction methods for T1 magnetic resonance images using region labeling and morphological operations

In this work we propose two brain extraction methods (BEM) that solely depend on the brain anatomy and its intensity characteristics. Our methods are simple, unsupervised and knowledge based. Using an adaptive intensity thresholding method on the magnetic resonance images of head scans, a binary image is obtained. The binary image is labeled using the anatomical facts that the scalp is the boundary between head and background, and the skull is the boundary separating brain and scalp. A run length scheme is applied on the labeled image to get a rough brain mask. Morphological operations are then performed to obtain the fine brain on the assumption that brain is the largest connected component (LCC). But the LCC concept failed to work on some slices where brain is composed of more than one connected component. To solve this problem a 3-D approach is introduced in the BEM. Experimental results on 61 sets of T1 scans taken from MRI scan center and neuroimage web services showed that our methods give better results than the popular methods, FSL's Brain Extraction Tool (BET), BrainSuite's Brain Surface Extractor (BSE) gives results comparable to that of Model-based Level Sets (MLS) and works well even where MLS failed. The average Dice similarity index computed using the "Gold standard" and the specificity values are 0.938 and 0.992, respectively, which are higher than that for BET, BSE and MLS. The average processing time by one of our methods is ≈1s/slice, which is smaller than for MLS, which is ≈4s/slice. One of our methods produces the lowest false positive rate of 0.075, which is smaller than that for BSE, BET and MLS. It is independent of imaging orientation and works well for slices with abnormal features like tumor and lesion in which the existing methods fail in certain cases.

Skull-stripping method for brain MRI using a 3D level set with a speedup operator

PURPOSE

To extract the brain region from brain magnetic resonance (MR) images using a fast 3D level set method and a refinement process.

MATERIALS AND METHODS

The proposed method introduces a speedup operator to the conventional 3D level set method in order to accelerate the level set evolution. While the processing time for brain extraction is reduced by the speedup operator, the accuracy of brain extraction is also improved by adopting a refinement process.

RESULTS

The speedup operator yielded a 75% reduction in the total iteration numbers for the synthesized volume. The proposed method was applied to several datasets and compared with previous methods, ie, BrainVisa, BET, and FreeSurfer. The proposed method provided a Jaccard index of 0.971 ± 0.0114 for the BrainWeb dataset, 0.864 ± 0.035 for the IBSR dataset, and 0.9414 ± 0.0517 for a self-produced dataset acquired with a 3T MRI system.

CONCLUSION

Utilizing a speedup operator, the proposed method reduced the evolution time. Robust and accurate results for various datasets were obtained in experiments.

Cerebral Contour Extraction with Particle Method in Neonatal MR Images

Cerebral surface extraction from neonatal MR images is the basic work of quantifying the deformation of the cerebrum. Although there are many conventional methods of segmenting the cerebral region, only the rough area is given by counting the number of surface voxels in the segmented region. This article proposes a new method of extraction that is based on the particle method. The method introduces three kinds of particles that correspond to cerebrospinal fluid, gray matter, and white matter; it converts the brain MR images into the set of particles. The proposed method was applied to neonatal magnetic resonance images, and the experimental results showed that the cerebral contour was extracted with a root-mean-square-error of 0.51 mm compared with the ground truth contour given by a physician.

Fully automatic brain extraction algorithm for axial T2-weighted magnetic resonance images

In this paper we propose two brain extraction algorithms (BEA) for T2-weighted magnetic resonance imaging (MRI) scans. The T2-weighted image is first filtered with a low pass filter (LPF) to remove or subdue the background noise. Then the image is diffused to enhance the brain boundaries. Using Ridler's method a threshold value for intensity is obtained. Using the threshold value a rough binary brain image is obtained. By performing morphological operations and using the largest connected component (LCC) analysis, a brain mask is obtained from which the brain is extracted. This method uses only 2D information of slices and is named as 2D-BEA. The concept of LCC failed in few slices. To overcome this problem, 3D information available in adjacent slices is used which resulted in 3D-BEA. Experimental results on 20 MRI data sets show that the proposed 3D-BEA gave excellent results. The performance of this 3D-BEA is better than 2D-BEA and other popular methods, brain extraction tool (BET) and brain surface extractor (BSE).

Semi-automated segmentation to assess the lateral meniscus in normal and osteoarthritic knees

OBJECTIVE

The goal of this study was to develop an algorithm to semi-automatically segment the meniscus in a series of magnetic resonance (MR) images to use for normal knees and those with moderate osteoarthritis (OA).

METHOD

The segmentation method was developed then evaluated on 10 baseline MR images obtained from subjects with no evidence, symptoms, or risk factors of knee (OA), and 14 from subjects with established knee OA enrolled in the Osteoarthritis Initiative (OAI). After manually choosing a seed point within the meniscus, a threshold level was calculated through a Gaussian fit model. Under anatomical, intensity, and range constraints, a threshold operation was completed followed by conditional dilation and post-processing. The post-processing operation reevaluates the pixels included and excluded in the area surrounding the meniscus to improve accuracy. The developed method was evaluated for both normal and degenerative menisci by comparing the segmentation algorithm results with manual segmentations from five human readers.

RESULTS

The semi-automated segmentation method produces results similar to those of trained observers, with an average similarity index over 0.80 for normal participants and 0.75, 0.67, and 0.64 for participants with established knee OA with Osteoarthritis Research Society International (OARSI) joint space narrowing (JSN) scores of 0, one, and two respectively.

CONCLUSION

The semi-automatic segmentation method produced accurate and consistent segmentations of the meniscus when compared to manual segmentations in the assessment of normal menisci in mild to moderate OA. Future studies will examine the change in volume, thickness, and intensity characteristics at different stages of OA.

A multidimensional segmentation evaluation for medical image data

Evaluation of segmentation methods is a crucial aspect in image processing, especially in the medical imaging field, where small differences between segmented regions in the anatomy can be of paramount importance. Usually, segmentation evaluation is based on a measure that depends on the number of segmented voxels inside and outside of some reference regions that are called gold standards. Although some other measures have been also used, in this work we propose a set of new similarity measures, based on different features, such as the location and intensity values of the misclassified voxels, and the connectivity and the boundaries of the segmented data. Using the multidimensional information provided by these measures, we propose a new evaluation method whose results are visualized applying a Principal Component Analysis of the data, obtaining a simplified graphical method to compare different segmentation results. We have carried out an intensive study using several classic segmentation methods applied to a set of MRI simulated data of the brain with several noise and RF inhomogeneity levels, and also to real data, showing that the new measures proposed here and the results that we have obtained from the multidimensional evaluation, improve the robustness of the evaluation and provides better understanding about the difference between segmentation methods.

Joint brain parametric T1-map segmentation and RF inhomogeneity calibration

We propose a constrained version of Mumford and Shah's (1989) segmentation model with an information-theoretic point of view in order to devise a systematic procedure to segment brain magnetic resonance imaging (MRI) data for parametric T(1)-Map and T(1)-weighted images, in both 2-D and 3D settings. Incorporation of a tuning weight in particular adds a probabilistic flavor to our segmentation method, and makes the 3-tissue segmentation possible. Moreover, we proposed a novel method to jointly segment the T(1)-Map and calibrate RF Inhomogeneity (JSRIC). This method assumes the average T(1) value of white matter is the same across transverse slices in the central brain region, and JSRIC is able to rectify the flip angles to generate calibrated T(1)-Maps. In order to generate an accurate T(1)-Map, the determination of optimal flip-angles and the registration of flip-angle images are examined. Our JSRIC method is validated on two human subjects in the 2D T(1)-Map modality and our segmentation method is validated by two public databases, BrainWeb and IBSR, of T(1)-weighted modality in the 3D setting.

Accurate and robust extraction of brain regions using a deformable model based on radial basis functions

Brain extraction from head magnetic resonance (MR) images is a classification problem of segmenting image volumes into brain and non-brain regions. It is a difficult task due to the convoluted brain surface and the inapparent brain/non-brain boundaries in images. This paper presents an automated, robust, and accurate brain extraction method which utilizes a new implicit deformable model to well represent brain contours and to segment brain regions from MR images. This model is described by a set of Wendland's radial basis functions (RBFs) and has the advantages of compact support property and low computational complexity. Driven by the internal force for imposing the smoothness constraint and the external force for considering the intensity contrast across boundaries, the deformable model of a brain contour can efficiently evolve from its initial state toward its target by iteratively updating the RBF locations. In the proposed method, brain contours are separately determined on 2D coronal and sagittal slices. The results from these two views are generally complementary and are thus integrated to obtain a complete 3D brain volume. The proposed method was compared to four existing methods, Brain Surface Extractor, Brain Extraction Tool, Hybrid Watershed Algorithm, and Model-based Level Set, by using two sets of MR images as well as manual segmentation results obtained from the Internet Brain Segmentation Repository. Our experimental results demonstrated that the proposed approach outperformed these four methods when jointly considering extraction accuracy and robustness.

Skull stripping based on region growing for magnetic resonance brain images

In this paper, we propose a new skull stripping method for T1-weighted magnetic resonance (MR) brain images. Skull stripping has played an important role in neuroimage research because it is a basic preliminary step in many clinical applications. The process of skull stripping can be challenging due to the complexity of the human brain, variable parameters of MR scanners, individual characteristics, etc. In this paper, we aim to develop a computationally efficient and robust method. In the proposed algorithm, after eliminating the background voxels with histogram analysis, two seed regions of the brain and non-brain regions were automatically identified using a mask produced by morphological operations. Then we expanded these seed regions with a 2D region growing algorithm based on general brain anatomy information. The proposed algorithm was validated using 56 volumes of human brain data and simulated phantom data with manually segmented masks. It was compared with two popular automated skull stripping methods: the brain surface extractor (BSE) and the brain extraction tool (BET). The experimental results showed that the proposed algorithm produced accurate and stable results against data sets acquired from various MR scanners and effectively addressed difficult problems such as low contrast and large anatomical connections between the brain and surrounding tissues. The proposed method was also robust against noise, RF, and intensity inhomogeneities.

Pervasive access to MRI bias artifact suppression service on a grid

Bias artifact corrupts MRIs in such a way that the image is afflicted by illumination variations. Some of the authors proposed the exponential entropy-driven homomorphic unsharp masking ( E(2)D-HUM) algorithm that corrects this artifact without any a priori hypothesis about the tissues or the MRI modality. Moreover, E(2)D-HUM does not care about the body part under examination and does not require any particular training task. People who want to use this algorithm, which is Matlab-based, have to set their own computers in order to execute it. Furthermore, they have to be Matlab-skilled to exploit all the features of the algorithm. In this paper, we propose to make such algorithm available as a service on a grid infrastructure, so that people can use it almost from everywhere, in a pervasive fashion, by means of a suitable user interface running on smartphones. The proposed solution allows physicians to use the E(2)D-HUM algorithm (or any other kind of algorithm, given that it is available as a service on the grid), being it remotely executed somewhere in the grid, and the results are sent back to the user's device. This way, physicians do not need to be aware of how to use Matlab to process their images. The pervasive service provision for medical image enhancement is presented, along with some experimental results obtained using smartphones connected to an existing Globus-based grid infrastructure.

A Scalable Framework For Segmenting Magnetic Resonance Images

A fast, accurate and fully automatic method of segmenting magnetic resonance images of the human brain is introduced. The approach scales well allowing fast segmentations of fine resolution images. The approach is based on modifications of the soft clustering algorithm, fuzzy c-means, that enable it to scale to large data sets. Two types of modifications to create incremental versions of fuzzy c-means are discussed. They are much faster when compared to fuzzy c-means for medium to extremely large data sets because they work on successive subsets of the data. They are comparable in quality to application of fuzzy c-means to all of the data. The clustering algorithms coupled with inhomogeneity correction and smoothing are used to create a framework for automatically segmenting magnetic resonance images of the human brain. The framework is applied to a set of normal human brain volumes acquired from different magnetic resonance scanners using different head coils, acquisition parameters and field strengths. Results are compared to those from two widely used magnetic resonance image segmentation programs, Statistical Parametric Mapping and the FMRIB Software Library (FSL). The results are comparable to FSL while providing significant speed-up and better scalability to larger volumes of data.

Image filtering via generalized scale

In medical imaging, low signal-to-noise ratio (SNR) and/or contrast-to-noise ratio (CNR) often cause many image processing algorithms to perform poorly. Postacquisition image filtering is an important off-line image processing approach widely employed to enhance the SNR and CNR. A major drawback of many filtering techniques is image degradation by diffusing/blurring edges and/or fine structures. In this paper, we introduce a scale-based filtering method that employs scale-dependent diffusion conductance to perform filtering. This approach utilizes novel object scale information via a concept called generalized scale, which imposes no shape, size, or anisotropic constraints unlike previously published ball scale-based filtering strategies. The object scale allows us to better control the filtering process by constraining smoothing in regions with fine details and in the vicinity of boundaries while permitting effective smoothing in the interior of homogeneous regions. A new quantitative evaluation strategy that captures the SNR to CNR trade-off behavior of filtering methods is presented. The evaluations based on the Brainweb data sets show superior performance of generalized scale-based diffusive filtering over two existing methods, namely, ball scale-based and nonlinear complex diffusion processes. Qualitative experiments based on both phantom and patient magnetic resonance images demonstrate that the generalized scale-based approach leads to better preservation of fine details and edges.

A divide and conquer deformable contour method with a model based searching algorithm

A divide and conquer deformable contour method is presented with an initial inside closed contour being divided into arbitrary segments, and these segments are allowed to deform separately preserving the segments' connectivity. A maximum area threshold, A/sub max/, is used to stop these outward contour segments' marching. Clear and blur contour points are then identified to partition the whole contour into clear and blur segments. A bi-directional searching method is then recursively applied to each blur segment including a search for contour-within-contour segment to reach a final close contour. Further improvements are provided by a model based searching algorithm. It is a two-step process with step 1 being a linked contour model matching operation where landmarks are extracted, and step 2 being a posteriori probability model matching and correction operation where large error segments are fine tuned to obtain the final results. The experiments include ultrasound images of pig heart, MRI brain images, MRI knee images having complex shapes with or without gaps, and inhomogeneous interior and contour region brightness distributions. These experiments have shown that the method has the capability of moving a contour into the neighboring region of the desired boundary by overcoming inhomogeneous interior, and by adapting each contour segment searching operation to different local difficulties, through a contour partition and repartition scheme in searching for a final solution.

A dual-stage method for lesion segmentation on digital mammograms

Mass lesion segmentation on mammograms is a challenging task since mass lesions are usually embedded and hidden in varying densities of parenchymal tissue structures. In this article, we present a method for automatic delineation of lesion boundaries on digital mammograms. This method utilizes a geometric active contour model that minimizes an energy function based on the homogeneities inside and outside of the evolving contour. Prior to the application of the active contour model, a radial gradient index (RGI)-based segmentation method is applied to yield an initial contour closer to the lesion boundary location in a computationally efficient manner. Based on the initial segmentation, an automatic background estimation method is applied to identify the effective circumstance of the lesion, and a dynamic stopping criterion is implemented to terminate the contour evolution when it reaches the lesion boundary. By using a full-field digital mammography database with 739 images, we quantitatively compare the proposed algorithm with a conventional region-growing method and an RGI-based algorithm by use of the area overlap ratio between computer segmentation and manual segmentation by an expert radiologist. At an overlap threshold of 0.4, 85% of the images are correctly segmented with the proposed method, while only 69% and 73% of the images are correctly delineated by our previous developed region-growing and RGI methods, respectively. This resulting improvement in segmentation is statistically significant.

Segmentation of brain parenchyma using bilateral filtering and region growing

When the non-diffusion weighted images (non-DWIs) and the diffusion weighted images (DWIs) are acquired by a fast imaging sequence, they suffer from several artifacts such as N/2 ghost, subject motion, eddy current, etc. These artifacts act as a noise in the background area of the human brain. To extract the brain region from the noisy background, brain parenchyma segmentation has been used. Several segmentation methods presented so far cannot address this problem well. In this study, we propose a novel segmentation method of brain contour in non-DWIs using bilateral filtering, which can reduce the background noise while edge-preserving, and region growing. We compare the segmentation results from various methods, and the proposed method shows better segmentation results than those from other schemes.

Virtual definition of neuronal tissue by cluster analysis of multi-parametric imaging (virtual-dot-com imaging)

Individual mapping of cerebral, morphological, functionally related structures using MRI was carried out using a new multi-contrast acquisition and analysis framework, called virtual-dot-com imaging. So far, conventional anatomical MRI has been able to provide gross segmentation of gray/white matter boundaries and a few sub-cortical structures. By combining a handful of imaging contrasts mechanisms (T1, T2, magnetization transfer, T2* and proton density), we were able to further segment sub-cortical tissue to its sub-nuclei arrangement, a segmentation that is difficult based on conventional, single-contrast MRI. Using an automatic four-step image and signal processing algorithm, we segmented the thalamus to at least 7 sub-nuclei with high similarity across subjects and high statistical significance within subjects (p<0.0001). The identified sub-nuclei resembled the known anatomical arrangement of the thalamus given in various atlases. Each cluster was characterized by a unique MRI contrast fingerprint. With this procedure, the weighted proportions of the different cellular compartments could be estimated, a property available to date only by histological analysis. Each sub-nucleus could be characterized in terms of normalized MRI contrast and compared to other sub-nuclei. The different weights of the contrasts (T1/T2/T2*/PD/MT, etc.) for each sub-nuclei cluster might indicate the intra-cluster morphological arrangement of the tissue that it represents. The implications of this methodology are far-ranging, from non-invasive, in vivo, individual mapping of histologically distinct brain areas to automatic identification of pathological processes.

Comparison between shifted Spearman rank correlation test and SPM in fMRI

A new improving method to compute Spearman rank correlation test has been developed. This method is based on the use of a response to stimulation delay referred by many authors. Visual results tend to prove the method efficiency which is confirmed by obtained overlap rates; moreover this method can easily be adapted in a clinical context.

Statistical morphological skull stripping of adult and infant MRI data

This paper describes a novel automatic statistical morphology skull stripper (SMSS) that uniquely exploits a statistical self-similarity measure and a 2-D brain mask to delineate the brain. The result of applying SMSS to 20 MRI data set volumes, including scans of both adult and infant subjects is also described. Quantitative performance assessment was undertaken with the use of brain masks provided by a brain segmentation expert. The performance is compared with an alternative technique known as brain extraction tool. The results suggest that SMSS is capable of skull-stripping neurological data with small amounts of over- and under-segmentation.

MRI tissue classification with neighborhood statistics: a nonparametric, entropy-minimizing approach

We introduce a novel approach for magnetic resonance image (MRI) brain tissue classification by learning image neighborhood statistics from noisy input data using nonparametric density estimation. The method models images as random fields and relies on minimizing an entropy-based metric defined on high dimensional probability density functions. Combined with an atlas-based initialization, it is completely automatic. Experiments on real and simulated data demonstrate the advantages of the method in comparison to other approaches.

New time-shifted Z-score and Student's test in fMRI

A new approach to compute z-score and Student's test in functional MRI has been developed. This approach tends to involve standard z-score and Student's test computation. This approach is based on the delay of the response compared to the stimulation introduced by many authors. The results obtained prove the methods efficiency; moreover these methods can be easily adapted in a clinical context. This paper presents the new computation and the validation.

Atlas-driven lung lobe segmentation in volumetric X-ray CT images

High-resolution X-ray computed tomography (CT) imaging is routinely used for clinical pulmonary applications. Since lung function varies regionally and because pulmonary disease is usually not uniformly distributed in the lungs, it is useful to study the lungs on a lobe-by-lobe basis. Thus, it is important to segment not only the lungs, but the lobar fissures as well. In this paper, we demonstrate the use of an anatomic pulmonary atlas, encoded with a priori information on the pulmonary anatomy, to automatically segment the oblique lobar fissures. Sixteen volumetric CT scans from 16 subjects are used to construct the pulmonary atlas. A ridgeness measure is applied to the original CT images to enhance the fissure contrast. Fissure detection is accomplished in two stages: an initial fissure search and a final fissure search. A fuzzy reasoning system is used in the fissure search to analyze information from three sources: the image intensity, an anatomic smoothness constraint, and the atlas-based search initialization. Our method has been tested on 22 volumetric thin-slice CT scans from 12 subjects, and the results are compared to manual tracings. Averaged across all 22 data sets, the RMS error between the automatically segmented and manually segmented fissures is 1.96 +/- 0.71 mm and the mean of the similarity indices between the manually defined and computer-defined lobe regions is 0.988. The results indicate a strong agreement between the automatic and manual lobe segmentations.

Intrathecal immunoglobulin G synthesis and brain injury by quantitative MRI in multiple sclerosis

OBJECTIVES

It was the aim of this study to evaluate if the quantitative intrathecal immunoglobulin G (IgG) synthesis correlates with the brain atrophy and the total lesion volume (TLV) in brain magnetic resonance imaging (MRI) of multiple sclerosis (MS) patients.

METHODS

A total of 50 patients with relapsing-remitting MS were included in this study. MRIs were performed and cerebrospinal fluid samples were collected during the diagnostic determination when patients were in remission without treatment.

RESULTS

At study baseline, IgG index values were elevated in 36 patients (72%), and oligoclonal IgG bands were positive in 42 of 50 patients (84%). Brain MRI was abnormal in 94% of patients, and, compared with healthy controls, brain atrophy was observed in MS patients. A positive correlation among IgG index, cerebrospinal fluid leukocyte count and TLV was observed; the Expanded Disability Status Scale correlated positively with TLV and the number of lesions, although a significant relationship between disability and brain atrophy was not demonstrated.

CONCLUSIONS

Although new parameters will be necessary in longitudinal studies to characterize the axonal injury in various stages of the disease, the data suggest that the high intrathecal IgG synthesis may predict a greater brain lesion burden.

STACS: new active contour scheme for cardiac MR image segmentation

The paper presents a novel stochastic active contour scheme (STACS) for automatic image segmentation designed to overcome some of the unique challenges in cardiac MR images such as problems with low contrast, papillary muscles, and turbulent blood flow. STACS minimizes an energy functional that combines stochastic region-based and edge-based information with shape priors of the heart and local properties of the contour. The minimization algorithm solves, by the level set method, the Euler-Lagrange equation that describes the contour evolution. STACS includes an annealing schedule that balances dynamically the weight of the different terms in the energy functional. Three particularly attractive features of STACS are: 1) ability to segment images with low texture contrast by modeling stochastically the image textures; 2) robustness to initial contour and noise because of the utilization of both edge and region-based information; 3) ability to segment the heart from the chest wall and the undesired papillary muscles due to inclusion of heart shape priors. Application of STACS to a set of 48 real cardiac MR images shows that it can successfully segment the heart from its surroundings such as the chest wall and the heart structures (the left and right ventricles and the epicardium.) We compare STACS' automatically generated contours with manually-traced contours, or the "gold standard," using both area and edge similarity measures. This assessment demonstrates very good and consistent segmentation performance of STACS.

An advanced boundary element method (BEM) implementation for the forward problem of electromagnetic source imaging

The forward problem of electromagnetic source imaging has two components: a numerical model to solve the related integral equations and a model of the head geometry. This study is on the boundary element method (BEM) implementation for numerical solutions and realistic head modelling. The use of second-order (quadratic) isoparametric elements and the recursive integration technique increase the accuracy in the solutions. Two new formulations are developed for the calculation of the transfer matrices to obtain the potential and magnetic field patterns using realistic head models. The formulations incorporate the use of the isolated problem approach for increased accuracy in solutions. If a personal computer is used for computations, each transfer matrix is calculated in 2.2 h. After this pre-computation period, solutions for arbitrary source configurations can be obtained in milliseconds for a realistic head model. A hybrid algorithm that uses snakes, morphological operations, region growing and thresholding is used for segmentation. The scalp, skull, grey matter, white matter and eyes are segmented from the multimodal magnetic resonance images and meshes for the corresponding surfaces are created. A mesh generation algorithm is developed for modelling the intersecting tissue compartments, such as eyes. To obtain more accurate results quadratic elements are used in the realistic meshes. The resultant BEM implementation provides more accurate forward problem solutions and more efficient calculations. Thus it can be the firm basis of the future inverse problem solutions.

Robust unsupervised segmentation of infarct lesion from diffusion tensor MR images using multiscale statistical classification and partial volume voxel reclassification

Manual region tracing method for segmentation of infarction lesions in images from diffusion tensor magnetic resonance imaging (DT-MRI) is usually used in clinical works, but it is time consuming. A new unsupervised method has been developed, which is a multistage procedure, involving image preprocessing, calculation of tensor field and measurement of diffusion anisotropy, segmentation of infarction volume based on adaptive multiscale statistical classification (MSSC), and partial volume voxel reclassification (PVVR). The method accounts for random noise, intensity overlapping, partial volume effect (PVE), and intensity shading artifacts, which always appear in DT-MR images. The proposed method was applied to 20 patients with clinically diagnosed brain infarction by DT-MRI scans. The accuracy and reproducibility in terms of identifying the infarction lesion have been confirmed by clinical experts. This automatic segmentation method is promising not only in detecting the location and the size of infarction lesion in stroke patient but also in quantitatively analyzing diffusion anisotropy of lesion to guide clinical diagnoses and therapy.