Adaptive segmentation of MRI data

Middle and inferior temporal gyrus gray matter volume abnormalities in chronic schizophrenia: an MRI study

OBJECTIVE

The middle temporal gyrus and inferior temporal gyrus subserve language and semantic memory processing, visual perception, and multimodal sensory integration. Functional deficits in these cognitive processes have been well documented in patients with schizophrenia. However, there have been few in vivo structural magnetic resonance imaging (MRI) studies of the middle temporal gyrus and inferior temporal gyrus in schizophrenia.

METHOD

Middle temporal gyrus and inferior temporal gyrus gray matter volumes were measured in 23 male patients diagnosed with chronic schizophrenia and 28 healthy male subjects by using high-spatial-resolution MRI. For comparison, superior temporal gyrus and fusiform gyrus gray matter volumes were also measured. Correlations between these four regions and clinical symptoms were also investigated.

RESULTS

Relative to healthy subjects, the patients with chronic schizophrenia showed gray matter volume reductions in the left middle temporal gyrus (13% difference) and bilateral inferior temporal gyrus (10% difference in both hemispheres). In addition, the patients showed gray matter volume reductions in the left superior temporal gyrus (13% difference) and bilateral fusiform gyrus (10% difference in both hemispheres). More severe hallucinations were significantly correlated with smaller left hemisphere volumes in the superior temporal gyrus and middle temporal gyrus.

CONCLUSIONS

These results suggest that patients with schizophrenia evince reduced gray matter volume in the left middle temporal gyrus and bilateral reductions in the inferior temporal gyrus. In conjunction with findings of left superior temporal gyrus reduction and bilateral fusiform gyrus reductions, these data suggest that schizophrenia may be characterized by left hemisphere-selective dorsal pathophysiology and bilateral ventral pathophysiology in temporal lobe gray matter.

Robust Image Segmentation Using FCM With Spatial Constraints Based on New Kernel-Induced Distance Measure

Fuzzy c-means clustering (FCM) with spatial constraints (FCM_S) is an effective algorithm suitable for image segmentation. Its effectiveness contributes not only to the introduction of fuzziness for belongingness of each pixel but also to exploitation of spatial contextual information. Although the contextual information can raise its insensitivity to noise to some extent, FCM_S still lacks enough robustness to noise and outliers and is not suitable for revealing non-Euclidean structure of the input data due to the use of Euclidean distance (L2 norm). In this paper, to overcome the above problems, we first propose two variants, FCM_S1 and FCM_S2, of FCM_S to aim at simplifying its computation and then extend them, including FCM_S, to corresponding robust kernelized versions KFCM_S, KFCM_S1 and KFCM_S2 by the kernel methods. Our main motives of using the kernel methods consist in: inducing a class of robust non-Euclidean distance measures for the original data space to derive new objective functions and thus clustering the non-Euclidean structures in data; enhancing robustness of the original clustering algorithms to noise and outliers, and still retaining computational simplicity. The experiments on the artificial and real-world datasets show that our proposed algorithms, especially with spatial constraints, are more effective.

Method for bias field correction of brain T1-weighted magnetic resonance images minimizing segmentation error

This work presents a new algorithm (nonuniform intensity correction; NIC) for correction of intensity inhomogeneities in T1-weighted magnetic resonance (MR) images. The bias field and a bias-free image are obtained through an iterative process that uses brain tissue segmentation. The algorithm was validated by means of realistic phantom images and a set of 24 real images. The first evaluation phase was based on a public domain phantom dataset, used previously to assess bias field correction algorithms. NIC performed similar to previously described methods in removing the bias field from phantom images, without introduction of degradation in the absence of intensity inhomogeneity. The real image dataset was used to compare the performance of this new algorithm to that of other widely used methods (N3, SPM'99, and SPM2). This dataset included both low and high bias field images from two different MR scanners of low (0.5 T) and medium (1.5 T) static fields. Using standard quality criteria for determining the goodness of the different methods, NIC achieved the best results, correcting the images of the real MR dataset, enabling its systematic use in images from both low and medium static field MR scanners. A limitation of our method is that it might fail if the bias field is so high that the initial histogram does not show bimodal distribution for white and gray matter.

A method for the automatic segmentation of autoradiographic image stacks and spatial normalization of functional cortical activity patterns

This paper introduces two new methods for the automatic anatomical and functional analysis of neurobiological autoradiographic image stacks, such as 2-fluoro-deoxyglucose (2FDG) images. The difficulty in the evaluation of these "2(1/2)D" datasets is that they do not inherently represent a continuous 3D data volume (as generated by MRI or CT), but consist of a stack of images from single tissue slices, suffering from unavoidable preparation artifacts. In the first part of the paper, a semi-automatic segmentation method is presented which generates a 3D surface model of certain brain structures and which is robust against different cutting directions with respect to the brain coordinate system. The method saves man-hours compared to manual segmentation and the results are highly reproducible. In the second part, a fully automatic method for the extraction, analysis and 3D visualization of functional information is described, which allows not only a more accurate localization of activation sites, but also greatly enhances the comparability of different individuals. Results are shown for 2FDG autoradiographs from rat brains under acoustical stimulation.

Probabilistic segmentation of white matter lesions in MR imaging

A new method has been developed for fully automated segmentation of white matter lesions (WMLs) in cranial MR imaging. The algorithm uses information from T1-weighted (T1-w), inversion recovery (IR), proton density-weighted (PD), T2-weighted (T2-w) and fluid attenuation inversion recovery (FLAIR) scans. It is based on the K-Nearest Neighbor (KNN) classification technique that builds a feature space from voxel intensities and spatial information. The technique generates images representing the probability per voxel being part of a WML. By application of thresholds on these probability maps, binary segmentations can be obtained. ROC curves show that the segmentations achieve both high sensitivity and specificity. A similarity index (SI), overlap fraction (OF) and extra fraction (EF) are calculated for additional quantitative analysis of the result. The SI is also used for determination of the optimal probability threshold for generation of the binary segmentation. Using probabilistic equivalents of the SI, OF and EF, the probability maps can be evaluated directly, providing a powerful tool for comparison of different classification results. This method for automated WML segmentation reaches an accuracy that is comparable to methods for multiple sclerosis (MS) lesion segmentation and is suitable for detection of WMLs in large and longitudinal population studies.

Automatic brain tumor segmentation by subject specific modification of atlas priors

RATIONALE AND OBJECTIVES

Manual segmentation of brain tumors from magnetic resonance images is a challenging and time-consuming task. An automated system has been developed for brain tumor segmentation that will provide objective, reproducible segmentations that are close to the manual results. Additionally, the method segments white matter, grey matter, cerebrospinal fluid, and edema. The segmentation of pathology and healthy structures is crucial for surgical planning and intervention.

MATERIALS AND METHODS

The method performs the segmentation of a registered set of magnetic resonance images using an expectation-maximization scheme. The segmentation is guided by a spatial probabilistic atlas that contains expert prior knowledge about brain structures. This atlas is modified with the subject-specific brain tumor prior that is computed based on contrast enhancement.

RESULTS

Five cases with different types of tumors are selected for evaluation. The results obtained from the automatic segmentation program are compared with results from manual and semi-automated methods. The automated method yields results that have surface distances at roughly 1-4 mm compared with the manual results.

CONCLUSION

The automated method can be applied to different types of tumors. Although its performance is below that of the semi-automated method, it has the advantage of requiring no user supervision.

Effects of heavy drinking, binge drinking, and family history of alcoholism on regional brain metabolites

BACKGROUND

The main goals are to investigate the effects of chronic active heavy drinking on N-acetylaspartate (NAA) and other metabolites throughout the brain and to determine whether they are affected by family history (FH) of alcoholism and long-term drinking pattern.

METHODS

Forty-six chronic heavy drinkers (HD) and 52 light drinkers (LD) were recruited from the community and compared on measures of regional brain structure using magnetic resonance imaging and measures of common brain metabolites in gray matter (GM) and white matter (WM) of the major lobes, subcortical nuclei, brainstem, and cerebellum using short-echo time magnetic resonance spectroscopic imaging. Regional atrophy-corrected levels of NAA, myoinositol (mI), creatine, and choline-containing metabolites were compared as a function of group, FH of alcoholism, and bingeing.

RESULTS

Frontal WM NAA was lower in FH-negative HD than FH-positive HD and tended to be lower in women than men. Creatine-containing metabolites in parietal GM were higher in HD than LD. FH-negative compared with FH-positive HD also had more mI in the brainstem and tended to have lower NAA and more mI in frontal GM. Although parietal GM NAA was not significantly lower in HD than LD, it was lower in non-binge drinkers than bingers. Frontal WM NAA was lower in HD than LD, with the difference driven by a small number of women, FH-negative HD, and older age. Lower frontal WM NAA in HD was associated with lower executive and working memory functions and with lower P3b amplitudes at frontal electrodes.

CONCLUSIONS

Community-dwelling HD who are not in alcoholism treatment have brain metabolite changes that are associated with lower brain function and are likely of behavioral significance. Age, FH, and binge drinking modulate brain metabolite abnormalities. Metabolite changes in active HD are less pronounced and present with a different spatial and metabolite pattern than reported in abstinent alcoholics.

Improved watershed transform for medical image segmentation using prior information

The watershed transform has interesting properties that make it useful for many different image segmentation applications: it is simple and intuitive, can be parallelized, and always produces a complete division of the image. However, when applied to medical image analysis, it has important drawbacks (oversegmentation, sensitivity to noise, poor detection of thin or low signal to noise ratio structures). We present an improvement to the watershed transform that enables the introduction of prior information in its calculation. We propose to introduce this information via the use of a previous probability calculation. Furthermore, we introduce a method to combine the watershed transform and atlas registration, through the use of markers. We have applied our new algorithm to two challenging applications: knee cartilage and gray matter/white matter segmentation in MR images. Numerical validation of the results is provided, demonstrating the strength of the algorithm for medical image segmentation.

Automated segmentation of human brain MR images using a multi-agent approach

Image interpretation consists in finding a correspondence between radiometric information and symbolic labeling with respect to specific spatial constraints. It is intrinsically a distributed process in terms of goals to be reached, zones in the image to be processed and methods to be applied. To cope the the difficulty inherent in this process, several information processing steps are required to gradually extract information. In this paper we advocate the use of situated cooperative agents as a framework for managing such steps. Dedicated agent behaviors are dynamically adapted depending on their position in the image, of their topographic relationships and of the radiometric information available. The information collected by the agents is gathered, shared via qualitative maps, or used as soon as available by acquaintances. Incremental refinement of interpretation is obtained through a coarse to fine strategy. Our work is essentially focused on radiometry-based tissue interpretation where knowledge is introduced or extracted at several levels to estimate models for tissue-intensity distribution and to cope with noise, intensity non-uniformity and partial volume effect. Several experiments on phantom and real images were performed. A complete volume can be segmented in less than 5 min with about 0.84% accuracy of the segmented reference. Comparison with other techniques demonstrates the potential interest of our approach for magnetic resonance imaging (MRI) brain scan interpretation.

Prefrontal cortical thickness in first-episode psychosis: a magnetic resonance imaging study

BACKGROUND

Findings from postmortem studies suggest reduced prefrontal cortical thickness in schizophrenia; however, cortical thickness in first-episode schizophrenia has not been evaluated using magnetic resonance imaging (MRI).

METHODS

Prefrontal cortical thickness was measured using MRI in first-episode schizophrenia patients (n = 17), first-episode affective psychosis patients (n = 17), and normal control subjects (n = 17); subjects were age-matched within 2 years and within a narrow age range (18-29 years). A previous study using the same subjects reported reduced prefrontal gray matter volume in first-episode schizophrenia. Manual editing was performed on those prefrontal segmentations before cortical thickness was measured.

RESULTS

Prefrontal cortical thickness was not significantly different among groups. Prefrontal gray matter volume and thickness were, however, positively correlated in both schizophrenia and control subjects. The product of boundary complexity and thickness, an alternative measure of volume, was positively correlated with volume for all three groups. Finally, age and age at first medication were negatively correlated with prefrontal cortical thickness only in first-episode schizophrenia.

CONCLUSIONS

This study demonstrates the potential usefulness of MRI for the study of cortical thickness abnormalities in schizophrenia. Correlations between cortical thickness and age and between cortical thickness and age at first medication suggest that the longer the schizophrenic process has been operative, the thinner the prefrontal cortex, although this needs confirmation in a longitudinal study.

Coupling Statistical Segmentation and PCA Shape Modeling

This paper presents a novel segmentation approach featuring shape constraints of multiple structures. A framework is developed combining statistical shape modeling with a maximum a posteriori segmentation problem. The shape is characterized by signed distance maps and its modes of variations are generated through principle component analysis. To solve the maximum a posteriori segmentation problem a robust Expectation Maximization implementation is used. The Expectation Maximization segmenter generates a label map, calculates image intensity inhomogeneities, and considers shape constraints for each structure of interest. Our approach enables high quality segmentations of structures with weak image boundaries which is demonstrated by automatically segmenting 32 brain MRIs into right and left thalami.

Accurate template-based correction of brain MRI intensity distortion with application to dementia and aging

This paper examines an alternative approach to separating magnetic resonance imaging (MRI) intensity inhomogeneity from underlying tissue-intensity structure using a direct template-based paradigm. This permits the explicit spatial modeling of subtle intensity variations present in normal anatomy which may confound common retrospective correction techniques using criteria derived from a global intensity model. A fine-scale entropy driven spatial normalisation procedure is employed to map intensity distorted MR images to a tissue reference template. This allows a direct estimation of the relative bias field between template and subject MR images, from the ratio of their low-pass filtered intensity values. A tissue template for an aging individual is constructed and used to correct distortion in a set of data acquired as part of a study on dementia. A careful validation based on manual segmentation and correction of nine datasets with a range of anatomies and distortion levels is carried out. This reveals a consistent improvement in the removal of global intensity variation in terms of the agreement with a global manual bias estimate, and in the reduction in the coefficient of intensity variation in manually delineated regions of white matter.

Statistics-based approach for aneurysm volume measurements

PURPOSE

To evaluate the ability of high-resolution MRA to monitor changes in intracranial aneurysm volume, and devise a highly reliable technique for obtaining these measurements.

MATERIALS AND METHODS

To obtain a baseline estimate of the repeatability of MRA scans and validate the statistics-based technique for aneurysm volume measurement, multiple scans were obtained on individual subjects over a period of up to 1 year. These 3D MRA data sets were coregistered and then analyzed using the volumetric analysis of segmented data and the proposed statistical method.

RESULTS

It was shown that high-resolution MRA provides highly repeatable data sets. Both methods used for the aneurysm volume measurements showed consistent results. However, the proposed statistical method had lower error and was much less sensitive to the choice of segmentation parameter than the volumetric analysis of segmented data. A change of 1 mm in the average radius of the aneurysm was detectable with the statistics-based technique.

CONCLUSIONS

This study demonstrates that the statistical method of aneurysm volume measurement in high-resolution MRA allows reliable and accurate assessments of aneurysm volume changes.

Sequence-independent segmentation of magnetic resonance images

We present a set of techniques for embedding the physics of the imaging process that generates a class of magnetic resonance images (MRIs) into a segmentation or registration algorithm. This results in substantial invariance to acquisition parameters, as the effect of these parameters on the contrast properties of various brain structures is explicitly modeled in the segmentation. In addition, the integration of image acquisition with tissue classification allows the derivation of sequences that are optimal for segmentation purposes. Another benefit of these procedures is the generation of probabilistic models of the intrinsic tissue parameters that cause MR contrast (e.g., T1, proton density, T2*), allowing access to these physiologically relevant parameters that may change with disease or demographic, resulting in nonmorphometric alterations in MR images that are otherwise difficult to detect. Finally, we also present a high band width multiecho FLASH pulse sequence that results in high signal-to-noise ratio with minimal image distortion due to B0 effects. This sequence has the added benefit of allowing the explicit estimation of T2* and of reducing test-retest intensity variability.

A fully automatic and robust brain MRI tissue classification method

A novel, fully automatic, adaptive, robust procedure for brain tissue classification from 3D magnetic resonance head images (MRI) is described in this paper. The procedure is adaptive in that it customizes a training set, by using a 'pruning' strategy, such that the classification is robust against anatomical variability and pathology. Starting from a set of samples generated from prior tissue probability maps (a 'model') in a standard, brain-based coordinate system ('stereotaxic space'), the method first reduces the fraction of incorrectly labeled samples in this set by using a minimum spanning tree graph-theoretic approach. Then, the corrected set of samples is used by a supervised kNN classifier for classifying the entire 3D image. The classification procedure is robust against variability in the image quality through a non-parametric implementation: no assumptions are made about the tissue intensity distributions. The performance of this brain tissue classification procedure is demonstrated through quantitative and qualitative validation experiments on both simulated MRI data (10 subjects) and real MRI data (43 subjects). A significant improvement in output quality was observed on subjects who exhibit morphological deviations from the model due to aging and pathology.

Dynamic programming generation of boundaries of local coordinatized submanifolds in the neocortex: application to the planum temporale

Dynamic programming is used to define boundaries of cortical submanifolds with focus on the planum temporale (PT) of the superior temporal gyrus (STG), which has been implicated in a variety of neuropsychiatric disorders. To this end, automated methods are used to generate the PT manifold from 10 high-resolution MRI subvolumes ROI masks encompassing the STG. A procedure to define the subvolume ROI masks from original MRI brain scans is developed. Bayesian segmentation is then used to segment the subvolumes into cerebrospinal fluid, gray matter (GM), and white matter (WM). 3D isocontouring using the intensity value at which there is equal probability of GM and WM is used to reconstruct the triangulated graph representing the STG cortical surface, enabling principal curvature at each point on the graph to be computed. Dynamic programming is used to delineate the PT manifold by tracking principal curves from the retro-insular end of the Heschl's gyrus (HG) to the STG, along the posterior STG up to the start of the ramus and back to the retro-insular end of the HG. A coordinate system is then defined on the PT manifold. The origin is defined by the retro-insular end of the HG and the y-axis passes through the point on the posterior STG where the ramus begins. Automated labeling of GM in the STG is robust with L(1) distances between Bayesian and manual segmentation in the range 0.001-0.12 (n = 20). PT reconstruction is also robust with 90% of the vertices of the reconstructed PT within about 1 voxel (n = 20) from semiautomated contours. Finally, the reliability index (based on interrater intraclass correlation) for the surface area derived from repeated reconstructions is 0.96 for the left PT and 0.94 for the right PT, thus demonstrating the robustness of dynamic programming in defining a coordinate system on the PT. It provides a method with potential significance in the study of neuropsychiatric disorders.

Abnormal CNV in chronic heavy drinkers

OBJECTIVE

We used the contingent negative variation (CNV), a slow negative shift in the human electroencephalogram, to investigate the effects of heavy chronic alcohol use on frontal lobe function.

METHODS

Event-related potentials (ERPs) were recorded from 30 heavy drinkers (HD) and 30 age-, sex-, and education-matched light or non-drinkers (LD), using a classical two-stimulus reaction time (RT) paradigm. Structural magnetic resonance images and neuropsychological tests were also administered.

RESULTS

The amplitude of the late CNV was significantly reduced in HD relative to light drinkers. Moreover, diminished CNV amplitudes in HD appear to be closely related to the amount of recent alcohol consumption. There were no significant differences in neuropsychological measures of frontal lobe function and frontal lobe volume between light and HD. However, in HD, reduced late CNV amplitudes were associated with decreased frontal lobe gray matter volume and poor performance on the Trail Making Test B. In LD but not in HD, late CNV amplitude correlated positively with RT, suggesting that the late CNV reflects some aspect of motor and cognitive preparation.

CONCLUSIONS

The inverse relationships between frontal lobe gray matter volume, performance on the Trail Making Test B, and late CNV amplitude in HD suggest that the ERP abnormalities observed in the current study may be indices of alcohol-related damage to the frontal lobe. The lack of a significant relationship between CNV amplitude and RT in HD suggests that chronic heavy alcohol use may disrupt response preparation.

Time-series analysis of MRI intensity patterns in multiple sclerosis

In progressive neurological disorders, such as multiple sclerosis (MS), magnetic resonance imaging (MRI) follow-up is used to monitor disease activity and progression and to understand the underlying pathogenic mechanisms. This article presents image postprocessing methods and validation for integrating multiple serial MRI scans into a spatiotemporal volume for direct quantitative evaluation of the temporal intensity profiles. This temporal intensity signal and its dynamics have thus far not been exploited in the study of MS pathogenesis and the search for MRI surrogates of disease activity and progression. The integration into a four-dimensional data set comprises stages of tissue classification, followed by spatial and intensity normalization and partial volume filtering. Spatial normalization corrects for variations in head positioning and distortion artifacts via fully automated intensity-based registration algorithms, both rigid and nonrigid. Intensity normalization includes separate stages of correcting intra- and interscan variations based on the prior tissue class segmentation. Different approaches to image registration, partial volume correction, and intensity normalization were validated and compared. Validation included a scan-rescan experiment as well as a natural-history study on MS patients, imaged in weekly to monthly intervals over a 1-year follow-up. Significant error reduction was observed by applying tissue-specific intensity normalization and partial volume filtering. Example temporal profiles within evolving multiple sclerosis lesions are presented. An overall residual signal variance of 1.4% +/- 0.5% was observed across multiple subjects and time points, indicating an overall sensitivity of 3% (for axial dual echo images with 3-mm slice thickness) for longitudinal study of signal dynamics from serial brain MRI.

Three-dimensional analysis of the geometry of individual multiple sclerosis lesions: detection of shape changes over time using spherical harmonics

PURPOSE

To suggest a quantitative method for assessing the temporal changes in the geometry of individual multiple sclerosis (MS) lesions in follow-up studies of MS patients.

MATERIALS AND METHODS

Computer simulated and in vivo magnetic resonance (MR) imaged MS lesions were studied. Ten in vivo MS lesions were identified from sets of axial MR images acquired from a patient scanned consecutively for 24 times during a one-year period. Each of the lesions was segmented and its three-dimensional surface approximated using spherical harmonics (SH). From the obtained SH polynomial coefficients, indices of shape were defined, and analysis of the temporal changes in each lesion's geometry throughout the year was performed by determining the mean discrete total variation of the shape indices.

RESULTS

The results demonstrate that most of the studied lesions undergo notable geometrical changes with time. These changes are not necessarily associated with similar changes in size/volume. Furthermore, it was found that indices corresponding to changes in lesion shape could be 1.4 to 8.0 times higher than those corresponding to changes in the lesion size/volume.

CONCLUSION

Quantitative three-dimensional shape analysis can serve as a new tool for monitoring MS lesion activity and study patterns of MS lesion evolution over time.

An adaptive spatial fuzzy clustering algorithm for 3-D MR image segmentation

An adaptive spatial fuzzy c-means clustering algorithm is presented in this paper for the segmentation of three-dimensional (3-D) magnetic resonance (MR) images. The input images may be corrupted by noise and intensity nonuniformity (INU) artifact. The proposed algorithm takes into account the spatial continuity constraints by using a dissimilarity index that allows spatial interactions between image voxels. The local spatial continuity constraint reduces the noise effect and the classification ambiguity. The INU artifact is formulated as a multiplicative bias field affecting the true MR imaging signal. By modeling the log bias field as a stack of smoothing B-spline surfaces, with continuity enforced across slices, the computation of the 3-D bias field reduces to that of finding the B-spline coefficients, which can be obtained using a computationally efficient two-stage algorithm. The efficacy of the proposed algorithm is demonstrated by extensive segmentation experiments using both simulated and real MR images and by comparison with other published algorithms.

Progressive decrease of left Heschl gyrus and planum temporale gray matter volume in first-episode schizophrenia: a longitudinal magnetic resonance imaging study

BACKGROUND

The Heschl gyrus and planum temporale have crucial roles in auditory perception and language processing. Our previous investigation using magnetic resonance imaging (MRI) indicated smaller gray matter volumes bilaterally in the Heschl gyrus and in left planum temporale in patients with first-episode schizophrenia but not in patients with first-episode affective psychosis. We sought to determine whether there are progressive decreases in anatomically defined MRI gray matter volumes of the Heschl gyrus and planum temporale in patients with first-episode schizophrenia and also in patients with first-episode affective psychosis.

METHODS

At a private psychiatric hospital, we conducted a prospective high spatial resolution MRI study that included initial scans of 28 patients at their first hospitalization (13 with schizophrenia and 15 with affective psychosis, 13 of whom had a manic psychosis) and 22 healthy control subjects. Follow-up scans occurred, on average, 1.5 years after the initial scan.

RESULTS

Patients with first-episode schizophrenia showed significant decreases in gray matter volume over time in the left Heschl gyrus (6.9%) and left planum temporale (7.2%) compared with patients with first-episode affective psychosis or control subjects.

CONCLUSIONS

These findings demonstrate a left-biased progressive volume reduction in the Heschl gyrus and planum temporale gray matter in patients with first-episode schizophrenia in contrast to patients with first-episode affective psychosis and control subjects. Schizophrenia but not affective psychosis seems to be characterized by a postonset progression of neocortical gray matter volume loss in the left superior temporal gyrus and thus may not be developmentally fixed.

A unified variational approach to denoising and bias correction in MR

We propose a novel bias correction method for magnetic resonance (MR) imaging that uses complementary body coil and surface coil images. The former are spatially homogeneous but have low signal intensity; the latter provide excellent signal response but have large bias fields. We present a variational framework where we optimize an energy functional to estimate the bias field and the underlying image using both observed images. The energy functional contains smoothness-enforcing regularization for both the image and the bias field. We present extensions of our basic framework to a variety of imaging protocols. We solve the optimization problem using a computationally efficient numerical algorithm based on coordinate descent, preconditioned conjugate gradient, half-quadratic regularization, and multigrid techniques. We show qualitative and quantitative results demonstrating the effectiveness of the proposed method in producing debiased and denoised MR images.

Automatic landmark extraction from image data using modified growing neural gas network

A new method for automatic landmark extraction from MR brain images is presented. In this method, landmark extraction is accomplished by modifying growing neural gas (GNG), which is a neural-network-based cluster-seeking algorithm. Using modified GNG (MGNG) corresponding dominant points of contours extracted from two corresponding images are found. These contours are borders of segmented anatomical regions from brain images. The presented method is compared to: 1) the node splitting-merging Kohonen model and 2) the Teh-Chin algorithm (a well-known approach for dominant points extraction of ordered curves). It is shown that the proposed algorithm has lower distortion error, ability of extracting landmarks from two corresponding curves simultaneously, and also generates the best match according to five medical experts.

Multicontext fuzzy clustering for separation of brain tissues in magnetic resonance images

A local image model is proposed to eliminate the adverse impact of both artificial and inherent intensity inhomogeneities in magnetic resonance imaging on intensity-based image segmentation methods. The estimation and correction procedures for intensity inhomogeneities are no longer indispensable because the highly convoluted spatial distribution of different tissues in the brain is taken into consideration. On the basis of the local image model, multicontext fuzzy clustering (MCFC) is proposed for classifying 2D and 3D MR data into tissues of white matter, gray matter, and cerebral spinal fluid automatically. In MCFC, multiple clustering contexts are generated for each pixel, and fuzzy clustering is independently performed in each context to calculate the degree of membership of a pixel to each tissue class. To maintain the statistical reliability and spatial continuity of membership distributions, a fusion strategy is adopted to integrate the clustering outcomes from different contexts. The fusion result is taken as the final membership value of the pixel. Experimental results on both real MR images and simulated volumetric MR data show that MCFC outperforms the classic fuzzy c-means (FCM) as well as other segmentation methods that deal with intensity inhomogeneities.

A unifying framework for partial volume segmentation of brain MR images

Accurate brain tissue segmentation by intensity-based voxel classification of magnetic resonance (MR) images is complicated by partial volume (PV) voxels that contain a mixture of two or more tissue types. In this paper, we present a statistical framework for PV segmentation that encompasses and extends existing techniques. We start from a commonly used parametric statistical image model in which each voxel belongs to one single tissue type, and introduce an additional downsampling step that causes partial voluming along the borders between tissues. An expectation-maximization approach is used to simultaneously estimate the parameters of the resulting model and perform a PV classification. We present results on well-chosen simulated images and on real MR images of the brain, and demonstrate that the use of appropriate spatial prior knowledge not only improves the classifications, but is often indispensable for robust parameter estimation as well. We conclude that general robust PV segmentation of MR brain images requires statistical models that describe the spatial distribution of brain tissues more accurately than currently available models.

Automated segmentation of multispectral brain MR images

This work presents a robust and comprehensive approach for the in vivo automated segmentation and quantitative tissue volume measurement of normal brain composition from multispectral magnetic resonance imaging (MRI) data. Statistical pattern recognition methods based on a finite mixture model are used to partition the intracranial volume into gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF) spaces. A masking algorithm initially extracts the brain volume from surrounding extrameningeal tissue. Radio frequency (RF) field inhomogeneity effects in the images are then removed using a recursive method that adapts to the intrinsic local tissue contrast. Our technique supports heterogeneous data with multispectral MR images of different contrast and intensity weighting acquired at varying spatial resolution and orientation. The proposed image segmentation methods have been tested using multispectral T1-, proton density-, and T2-weighted MRI data from young and aged non-human primates as well as from human subjects.

Automatic brain tissue extraction method using erosion-dilation treatment (BREED) from three-dimensional magnetic resonance imaging T1-weighted data

To improve the efficiency of brain image analysis, we propose a full-automatic method for extracting brain tissue from three-dimensional magnetic resonance imaging of T1-weighted data on the human head (brain tissue extraction method using erosion-dilation treatment [BREED]). The extraction processing is realized by combining signal intensity thresholding by means of the discriminant analysis method and an erosion-dilation treatment of the image. The accuracy of BREED is evaluated using both simulated and subject data. BREED can extract brain tissues with high accuracy (approximately 97%) for either simulated or subject data.

Automatic "pipeline" analysis of 3-D MRI data for clinical trials: application to multiple sclerosis

The quantitative analysis of magnetic resonance imaging (MRI) data has become increasingly important in both research and clinical studies aiming at human brain development, function, and pathology. Inevitably, the role of quantitative image analysis in the evaluation of drug therapy will increase, driven in part by requirements imposed by regulatory agencies. However, the prohibitive length of time involved and the significant intraand inter-rater variability of the measurements obtained from manual analysis of large MRI databases represent major obstacles to the wider application of quantitative MRI analysis. We have developed a fully automatic "pipeline" image analysis framework and have successfully applied it to a number of large-scale, multicenter studies (more than 1,000 MRI scans). This pipeline system is based on robust image processing algorithms, executed in a parallel, distributed fashion. This paper describes the application of this system to the automatic quantification of multiple sclerosis lesion load in MRI, in the context of a phase III clinical trial. The pipeline results were evaluated through an extensive validation study, revealing that the obtained lesion measurements are statistically indistinguishable from those obtained by trained human observers. Given that intra- and inter-rater measurement variability is eliminated by automatic analysis, this system enhances the ability to detect small treatment effects not readily detectable through conventional analysis techniques. While useful for clinical trial analysis in multiple sclerosis, this system holds widespread potential for applications in other neurological disorders, as well as for the study of neurobiology in general.

Incorporating Non-rigid Registration into Expectation Maximization Algorithm to Segment MR Images

The paper introduces an algorithm which allows the automatic segmentation of multi channel magnetic resonance images. We extended the Expectation Maximization-Mean Field Approximation Segmenter, to include Local Prior Probability Maps. Thereby our algorithm estimates the bias field in the image while simultaneously assigning voxels to different tissue classes under prior probability maps. The probability maps were aligned to the subject using nonrigid registration. This allowed the parcellation of cortical sub-structures including the superior temporal gyrus. To our knowledge this is the first description of an algorithm capable of automatic cortical parcellation incorporating strong noise reduction and image intensity correction.

An accurate and efficient bayesian method for automatic segmentation of brain MRI

Automatic three-dimensional (3-D) segmentation of the brain from magnetic resonance (MR) scans is a challenging problem that has received an enormous amount of attention lately. Of the techniques reported in the literature, very few are fully automatic. In this paper, we present an efficient and accurate, fully automatic 3-D segmentation procedure for brain MR scans. It has several salient features; namely, the following. 1) Instead of a single multiplicative bias field that affects all tissue intensities, separate parametric smooth models are used for the intensity of each class. 2) A brain atlas is used in conjunction with a robust registration procedure to find a nonrigid transformation that maps the standard brain to the specimen to be segmented. This transformation is then used to: segment the brain from nonbrain tissue; compute prior probabilities for each class at each voxel location and find an appropriate automatic initialization. 3) Finally, a novel algorithm is presented which is a variant of the expectation-maximization procedure, that incorporates a fast and accurate way to find optimal segmentations, given the intensity models along with the spatial coherence assumption. Experimental results with both synthetic and real data are included, as well as comparisons of the performance of our algorithm with that of other published methods.

MRI study of caudate nucleus volume and its cognitive correlates in neuroleptic-naive patients with schizotypal personality disorder

OBJECTIVE

"Cognitive" circuits anatomically link the frontal lobe to subcortical structures; therefore, pathology in any of the core components of these circuits, such as in the caudate nucleus, may result in neurobehavioral syndromes similar to those of the frontal lobe. Neuroleptic medication, however, affects the size of the caudate nucleus. For this reason, individuals diagnosed with schizotypal personality disorder offer an ideal group for the measurement of the caudate nucleus because they may be genetically related to individuals with schizophrenia but do not require neuroleptic treatment because of their less severe symptoms.

METHOD

Magnetic resonance imagining (MRI) scans obtained on a 1.5-T magnet with 1.5-mm contiguous slices were used to measure the caudate nucleus and lateral ventricles in 15 right-handed male subjects with schizotypal personality disorder who had no previous neuroleptic exposure and in 14 normal comparison subjects. Subjects were group matched for parental socioeconomic status, handedness, and gender.

RESULTS

First, the authors found significantly lower left and right absolute (13.1%, 13.2%) and relative (9.1%, 9.2%) caudate nucleus volumes in never-medicated subjects with schizotypal personality disorder than in normal subjects. Second, they found significant, inverse correlations between caudate nucleus volume and the severity of perseveration in two distinct working memory tasks in these neuroleptic-naive subjects with schizotypal personality disorder.

CONCLUSIONS

These data are consistent with the findings of reduced caudate nucleus volume reported in studies of neuroleptic-naive patients experiencing their first episode of schizophrenia and support the association of intrinsic pathology in the caudate nucleus with abnormalities in working memory in the schizophrenia spectrum.

BrainSuite: an automated cortical surface identification tool

We describe a new magnetic resonance (MR) image analysis tool that produces cortical surface representations with spherical topology from MR images of the human brain. The tool provides a sequence of low-level operations in a single package that can produce accurate brain segmentations in clinical time. The tools include skull and scalp removal, image nonuniformity compensation, voxel-based tissue classification, topological correction, rendering, and editing functions. The collection of tools is designed to require minimal user interaction to produce cortical representations. In this paper we describe the theory of each stage of the cortical surface identification process. We then present classification validation results using real and phantom data. We also present a study of interoperator variability.

Shape recovery algorithms using level sets in 2-D/3-D medical imagery: a state-of-the-art review

The class of geometric deformable models, also known as level sets, has brought tremendous impact to medical imagery due to its capability of topology preservation and fast shape recovery. In an effort to facilitate a clear and full understanding of these powerful state-of-the-art applied mathematical tools, this paper is an attempt to explore these geometric methods, their implementations and integration of regularizers to improve the robustness of these topologically independent propagating curves/surfaces. This paper first presents the origination of level sets, followed by the taxonomy of level sets. We then derive the fundamental equation of curve/surface evolution and zero-level curves/surfaces. The paper then focuses on the first core class of level sets, known as "level sets without regularizers." This class presents five prototypes: gradient, edge, area-minimization, curvature-dependent and application driven. The next section is devoted to second core class of level sets, known as "level sets with regularizers." In this class, we present four kinds: clustering-based, Bayesian bidirectional classifier-based, shape-based and coupled constrained-based. An entire section is dedicated to optimization and quantification techniques for shape recovery when used in the level set framework. Finally, the paper concludes with 22 general merits and four demerits on level sets and the future of level sets in medical image segmentation. We present applications of level sets to complex shapes like the human cortex acquired via MRI for neurological image analysis.

Automatic MRI adipose tissue mapping using overlapping mosaics

This paper presents an automatic method of correcting non-uniform RF coil response for the classification of body composition using MR imaging. By linear mosaic modelling, the smoothly but non-linearly varying bias field, which modulates tissue intensities within the image, was corrected. The overlapping between adjacent mosaics ensured consistent segmentation of body fat content and the effectiveness of the technique was validated by both phantom and in vivo experiments. Ten whole body composition data sets, each with 39 trans-axial slices, were acquired. Automatic segmentation results using the proposed technique were compared with those from manual delineations. The automatic segmentation method was found to be highly accurate and the mean percentage error between the two methods was less than 1.5%.

A modified fuzzy C-means algorithm for bias field estimation and segmentation of MRI data

In this paper, we present a novel algorithm for fuzzy segmentation of magnetic resonance imaging (MRI) data and estimation of intensity inhomogeneities using fuzzy logic. MRI intensity inhomogeneities can be attributed to imperfections in the radio-frequency coils or to problems associated with the acquisition sequences. The result is a slowly varying shading artifact over the image that can produce errors with conventional intensity-based classification. Our algorithm is formulated by modifying the objective function of the standard fuzzy c-means (FCM) algorithm to compensate for such inhomogeneities and to allow the labeling of a pixel (voxel) to be influenced by the labels in its immediate neighborhood. The neighborhood effect acts as a regularizer and biases the solution toward piecewise-homogeneous labelings. Such a regularization is useful in segmenting scans corrupted by salt and pepper noise. Experimental results on both synthetic images and MR data are given to demonstrate the effectiveness and efficiency of the proposed algorithm.

Quantitative analysis of MRI signal abnormalities of brain white matter with high reproducibility and accuracy

PURPOSE

To assess the reproducibility and accuracy compared to radiologists of three automated segmentation pipelines for quantitative magnetic resonance imaging (MRI) measurement of brain white matter signal abnormalities (WMSA).

MATERIALS AND METHODS

WMSA segmentation was performed on pairs of whole brain scans from 20 patients with multiple sclerosis (MS) and 10 older subjects who were positioned and imaged twice within 30 minutes. Radiologist outlines of WMSA on 20 sections from 16 patients were compared with the corresponding results of each segmentation method.

RESULTS

The segmentation method combining expectation-maximization (EM) tissue segmentation, template-driven segmentation (TDS), and partial volume effect correction (PVEC) demonstrated the highest accuracy (the absolute value of the Z-score was 0.99 for both groups of subjects), as well as high interscan reproducibility (repeatability coefficient was 0.68 mL in MS patients and 1.49 mL in aging subjects).

CONCLUSION

The addition of TDS to the EM segmentation and PVEC algorithms significantly improved the accuracy of WMSA volume measurements, while also improving measurement reproducibility.

Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain

We present a technique for automatically assigning a neuroanatomical label to each voxel in an MRI volume based on probabilistic information automatically estimated from a manually labeled training set. In contrast to existing segmentation procedures that only label a small number of tissue classes, the current method assigns one of 37 labels to each voxel, including left and right caudate, putamen, pallidum, thalamus, lateral ventricles, hippocampus, and amygdala. The classification technique employs a registration procedure that is robust to anatomical variability, including the ventricular enlargement typically associated with neurological diseases and aging. The technique is shown to be comparable in accuracy to manual labeling, and of sufficient sensitivity to robustly detect changes in the volume of noncortical structures that presage the onset of probable Alzheimer's disease.

Robotics for surgery

Robotic technology is enhancing surgery through improved precision, stability, and dexterity. In image-guided procedures, robots use magnetic resonance and computed tomography image data to guide instruments to the treatment site. This requires new algorithms and user interfaces for planning procedures; it also requires sensors for registering the patient's anatomy with the preoperative image data. Minimally invasive procedures use remotely controlled robots that allow the surgeon to work inside the patient's body without making large incisions. Specialized mechanical designs and sensing technologies are needed to maximize dexterity under these access constraints. Robots have applications in many surgical specialties. In neurosurgery, image-guided robots can biopsy brain lesions with minimal damage to adjacent tissue. In orthopedic surgery, robots are routinely used to shape the femur to precisely fit prosthetic hip joint replacements. Robotic systems are also under development for closed-chest heart bypass, for microsurgical procedures in ophthalmology, and for surgical training and simulation. Although results from initial clinical experience is positive, issues of clinician acceptance, high capital costs, performance validation, and safety remain to be addressed.