Pulmonary Nodule Classification Using a Multiview Residual Selective Kernel Network.

Lung cancer is one of the leading causes of death worldwide and early detection is crucial to reduce the mortality. A reliable computer-aided diagnosis (CAD) system can help facilitate early detection of malignant nodules. Although existing methods provide adequate classification accuracy, there is still room for further improvement. This study is dedicated to investigating a new CAD scheme for predicting the malignant likelihood of lung nodules in computed tomography (CT) images in light of a deep learning strategy. Conceived from the residual learning and selective kernel, we investigated an efficient residual selective kernel (RSK) block to handle the diversity of lung nodules with various shapes and obscure structures. Founded on this RSK block, we established a multiview RSK network (MRSKNet), to which three anatomical planes in the axial, coronal, and sagittal directions were fed. To reinforce the classification efficiency, seven handcrafted texture features with a filter-like computation strategy were explored, among which the homogeneity (HOM) feature maps are combined with the corresponding intensity CT images for concatenation input, leading to an improved network architecture. Evaluated on the public benchmark Lung Image Database Consortium and Image Database Resource Initiative (LIDC-IDRI) challenge database with ten-fold cross validation of binary classification, our experimental results indicated high area under receiver operating characteristic (AUC) and accuracy scores. A better compromise between recall and specificity was struck using the suggested concatenation strategy comparing to many state-of-the-art approaches. The proposed pulmonary nodule classification framework exhibited great efficacy and achieved a higher AUC of 0.9711. The association of handcrafted texture features with deep learning models is promising in advancing the classification performance. The developed pulmonary nodule CAD network architecture is of potential in facilitating the diagnosis of lung cancer for further image processing applications.

Automated Stroke Lesion Segmentation in Rat Brain MR Images Using an Encoder-Decoder Framework.

Stroke is a leading cause of serious long-term disability and the major cause of mortality worldwide. Experimental ischemic stroke models play an important role in realizing the mechanism of cerebral ischemia and evaluating the development of pathological extent. An accurate and reliable image segmentation tool to automatically identify the stroke lesion is important in the subsequent processes. However, the intensity distribution of the infarct region in the diffusion weighted imaging (DWI) images is usually nonuniform with blurred boundaries. A deep learning-based infarct region segmentation framework is developed in this paper to address the segmentation difficulties. The proposed solution is an encoder-decoder network that includes a hybrid block model for efficient multiscale feature extraction. An in-house DWI image dataset was created to evaluate this automated stroke lesion segmentation scheme. Through massive experiments, accurate segmentation results were obtained, which outperformed many competitive methods both qualitatively and quantitatively. Our stroke lesion segmentation system is potential in providing a decent tool to facilitate preclinical stroke investigation using DWI images.Clinical Relevance- This facilitates neuroscientists the investigation of a new scoring system with less examination time and better inter-rater reliability, which helps to understand the function of specific brain areas underlying neuroimaging signatures clinically.

Multimodal Image Registration Using a Viscous Fluid Model with the Bhattacharyya Distance.

Image registration is an elementary task in medical image processing and analysis, which can be divided into monomodal and multimodal. Direct 3D multimodal registration in volumetric medical images can provide more insight into the interpretation of subsequent image processing applications than 2D methods. This paper is dedicated to the development of a 3D multimodal image registration algorithm based on a viscous fluid model associated with the Bhattacharyya distance. In our approach, a modified Navier-Stoke's equation is exploited as the foundation of the multimodal image registration framework. The hopscotch method is numerically implemented to solve the velocity field, whose values at the explicit locations are first computed and the values at the implicit positions are solved by transposition. The differential of the Bhattacharyya distance is incorporated into the body force function, which is the main driving force for deformation, to enable multimodal registration. A variety of simulated and real brain MR images were utilized to assess the proposed 3D multimodal image registration system. Preliminary experimental results indicated that our algorithm produced high registration accuracy in various registration scenarios and outperformed other competing methods in many multimodal image registration tasks.Clinical Relevance- This facilitates the disease diagnosis and treatment planning that requires accurate 3D multimodal image registration without massive image data and extensive training regardless of the imaging modality.

RU-Net: skull stripping in rat brain MR images after ischemic stroke with rat U-Net.

BACKGROUND: Experimental ischemic stroke models play a fundamental role in interpreting the mechanism of cerebral ischemia and appraising the development of pathological extent. An accurate and automatic skull stripping tool for rat brain image volumes with magnetic resonance imaging (MRI) are crucial in experimental stroke analysis. Due to the deficiency of reliable rat brain segmentation methods and motivated by the demand for preclinical studies, this paper develops a new skull stripping algorithm to extract the rat brain region in MR images after stroke, which is named Rat U-Net (RU-Net). METHODS: Based on a U-shape like deep learning architecture, the proposed framework integrates batch normalization with the residual network to achieve efficient end-to-end segmentation. A pooling index transmission mechanism between the encoder and decoder is exploited to reinforce the spatial correlation. Two different modalities of diffusion-weighted imaging (DWI) and T2-weighted MRI (T2WI) corresponding to two in-house datasets with each consisting of 55 subjects were employed to evaluate the performance of the proposed RU-Net. RESULTS: Extensive experiments indicated great segmentation accuracy across diversified rat brain MR images. It was suggested that our rat skull stripping network outperformed several state-of-the-art methods and achieved the highest average Dice scores of 98.04% (p < 0.001) and 97.67% (p < 0.001) in the DWI and T2WI image datasets, respectively. CONCLUSION: The proposed RU-Net is believed to be potential for advancing preclinical stroke investigation and providing an efficient tool for pathological rat brain image extraction, where accurate segmentation of the rat brain region is fundamental.

Segmentation of Rat Brains and Cerebral Hemispheres in Triphenyltetrazolium Chloride-Stained Images after Stroke.

Ischemic stroke is one of the leading causes of death among the aged population in the world. Experimental stroke models with rodents play a fundamental role in the investigation of the mechanism and impairment of cerebral ischemia. For its celerity and veracity, the 2,3,5-triphenyltetrazolium chloride (TTC) staining of rat brains has been extensively adopted to visualize the infarction, which is subsequently photographed for further processing. Two important tasks are to segment the brain regions and to compute the midline that separates the brain. This paper investigates automatic brain extraction and hemisphere segmentation algorithms in camera-based TTC-stained rat images. For rat brain extraction, a saliency region detection scheme on a superpixel image is exploited to extract the brain regions from the raw complicated image. Subsequently, the initial brain slices are refined using a parametric deformable model associated with color image transformation. For rat hemisphere segmentation, open curve evolution guided by the gradient vector flow in a medial subimage is developed to compute the midline. A wide variety of TTC-stained rat brain images captured by a smartphone were produced and utilized to evaluate the proposed segmentation frameworks. Experimental results on the segmentation of rat brains and cerebral hemispheres indicated that the developed schemes achieved high accuracy with average Dice scores of 92.33% and 97.15%, respectively. The established segmentation algorithms are believed to be potential and beneficial to facilitate experimental stroke study with TTC-stained rat brain images.

Automatic brain extraction and hemisphere segmentation in rat brain MR images after stroke using deformable models.

PURPOSE: Experimental ischemic stroke models play an essential role in understanding the mechanisms of cerebral ischemia and evaluating the development of pathological extent. An important precursor to the investigation of ischemic strokes associated with rodents is the brain extraction and hemisphere segmentation in rat brain diffusion-weighted imaging (DWI) and T2-weighted MRI (T2WI) images. Accurate and reliable image segmentation tools for extracting the rat brain and hemispheres in the MR images are critical in subsequent processes, such as lesion identification and injury analysis. This study is an attempt to investigate rat brain extraction and hemisphere segmentation algorithms that are practicable in both DWI and T2WI images. METHODS: To automatically perform brain extraction, the proposed framework is based on an efficient geometric deformable model. By introducing an additional image force in response to the rat brain characteristics into the skull stripping model, we establish a unique rat brain extraction scheme in DWI and T2WI images. For the subsequent hemisphere segmentation, we develop an efficient brain feature detection algorithm to approximately separate the rat brain. A refinement process is enforced by constructing a gradient vector flow in the proximity of the midsurface, where a parametric active contour is attracted to achieve hemisphere segmentation. RESULTS: Extensive experiments with 55 DWI and T2WI subjects were executed in comparison with the state-of-the-art methods. Experimental results indicated that our rat brain extraction and hemisphere segmentation schemes outperformed the competitive methods and exhibited high performance both qualitatively and quantitatively. For rat brain extraction, the average Dice scores were 97.13% and 97.42% in DWI and T2WI image volumes, respectively. Rat hemisphere segmentation results based on the Hausdorff distance metric revealed average values of 0.17 and 0.15 mm for DWI and T2WI subjects, respectively. CONCLUSIONS: We believe that the established frameworks are advantageous to facilitate preclinical stroke investigation and relevant neuroscience research that requires accurate brain extraction and hemisphere segmentation using rat DWI and T2WI images.

An automatic restoration framework based on GPU-accelerated collateral filtering in brain MR images.

BACKGROUND: Image restoration is one of the fundamental and essential tasks within image processing. In medical imaging, developing an effective algorithm that can automatically remove random noise in brain magnetic resonance (MR) images is challenging. The collateral filter has been shown a more powerful algorithm than many existing methods. However, the computation of the collateral filter is more time-consuming and the selection of the filter parameters is also laborious. This paper proposes an automatic noise removal system based on the accelerated collateral filter for brain MR images. METHODS: To solve these problems, we first accelerated the collateral filter with parallel computing using the graphics processing unit (GPU) architecture. We adopted the compute unified device architecture (CUDA), an application programming interface for the GPU by NVIDIA, to hasten the computation. Subsequently, the optimal filter parameters were selected and the automation was achieved by artificial neural networks. Specifically, an artificial neural network system associated with image feature analysis was adopted to establish the automatic image restoration framework. The best feature combination was selected by the paired t-test and the sequential forward floating selection (SFFS) methods. RESULTS: Experimental results indicated that not only did the proposed automatic image restoration algorithm perform dramatically faster than the traditional collateral filter, but it also effectively removed the noise in a wide variety of brain MR images. A speed up gain of 34 was attained to process an MR image, which completed within 0.1 s. Representative illustrations of brain tumor images demonstrated the capability of identifying lesion boundaries, which outperformed many existing methods. CONCLUSIONS: We believe that our accelerated and automated restoration framework is promising for achieving robust filtering in many brain MR image restoration applications.

An Automatic Parameter Decision System of Bilateral Filtering with GPU-Based Acceleration for Brain MR Images.

Bilateral filters have been extensively utilized in a number of image denoising applications such as segmentation, registration, and tissue classification. However, it requires burdensome adjustments of the filter parameters to achieve the best performance for each individual image. To address this problem, this paper proposes a computer-aided parameter decision system based on image texture features associated with neural networks. In our approach, parallel computing with the GPU architecture is first developed to accelerate the computation of the conventional bilateral filter. Subsequently, a back propagation network (BPN) scheme using significant image texture features as the input is established to estimate the GPU-based bilateral filter parameters and its denoising process. The k-fold cross validation method is exploited to evaluate the performance of the proposed automatic restoration framework. A wide variety of T1-weighted brain MR images were employed to train and evaluate this parameter-free decision system with GPU-based bilateral filtering, which resulted in a speed-up factor of 208 comparing to the CPU-based computation. The proposed filter parameter prediction system achieved a mean absolute percentage error (MAPE) of 6% and was classified as "high accuracy". Our automatic denoising framework dramatically removed noise in numerous brain MR images and outperformed several state-of-the-art methods based on the peak signal-to-noise ratio (PSNR). The usage of image texture features associated with the BPN to estimate the GPU-based bilateral filter parameters and to automate the denoising process is feasible and validated. It is suggested that this automatic restoration system is advantageous to various brain MR image-processing applications.

Brain MR Image Restoration Using an Automatic Trilateral Filter With GPU-Based Acceleration.

OBJECTIVE: Noise reduction in brain magnetic resonance (MR) images has been a challenging and demanding task. This study develops a new trilateral filter that aims to achieve robust and efficient image restoration. METHODS: Extended from the bilateral filter, the proposed algorithm contains one additional intensity similarity funct-ion, which compensates for the unique characteristics of noise in brain MR images. An entropy function adaptive to intensity variations is introduced to regulate the contributions of the weighting components. To hasten the computation, parallel computing based on the graphics processing unit (GPU) strategy is explored with emphasis on memory allocations and thread distributions. To automate the filtration, image texture feature analysis associated with machine learning is investigated. Among the 98 candidate features, the sequential forward floating selection scheme is employed to acquire the optimal texture features for regularization. Subsequently, a two-stage classifier that consists of support vector machines and artificial neural networks is established to predict the filter parameters for automation. RESULTS: A speedup gain of 757 was reached to process an entire MR image volume of 256 × 256 × 256 pixels, which completed within 0.5 s. Automatic restoration results revealed high accuracy with an ensemble average relative error of 0.53 ± 0.85% in terms of the peak signal-to-noise ratio. CONCLUSION: This self-regulating trilateral filter outperformed many state-of-the-art noise reduction methods both qualitatively and quantitatively. SIGNIFICANCE: We believe that this new image restoration algorithm is of potential in many brain MR image processing applications that require expedition and automation.

CUDA-based acceleration and BPN-assisted automation of bilateral filtering for brain MR image restoration.

PURPOSE: Bilateral filters have been substantially exploited in numerous magnetic resonance (MR) image restoration applications for decades. Due to the deficiency of theoretical basis on the filter parameter setting, empirical manipulation with fixed values and noise variance-related adjustments has generally been employed. The outcome of these strategies is usually sensitive to the variation of the brain structures and not all the three parameter values are optimal. This article is in an attempt to investigate the optimal setting of the bilateral filter, from which an accelerated and automated restoration framework is developed. METHODS: To reduce the computational burden of the bilateral filter, parallel computing with the graphics processing unit (GPU) architecture is first introduced. The NVIDIA Tesla K40c GPU with the compute unified device architecture (CUDA) functionality is specifically utilized to emphasize thread usages and memory resources. To correlate the filter parameters with image characteristics for automation, optimal image texture features are subsequently acquired based on the sequential forward floating selection (SFFS) scheme. Subsequently, the selected features are introduced into the back propagation network (BPN) model for filter parameter estimation. Finally, the k-fold cross validation method is adopted to evaluate the accuracy of the proposed filter parameter prediction framework. RESULTS: A wide variety of T1-weighted brain MR images with various scenarios of noise levels and anatomic structures were utilized to train and validate this new parameter decision system with CUDA-based bilateral filtering. For a common brain MR image volume of 256 × 256 × 256 pixels, the speed-up gain reached 284. Six optimal texture features were acquired and associated with the BPN to establish a "high accuracy" parameter prediction system, which achieved a mean absolute percentage error (MAPE) of 5.6%. Automatic restoration results on 2460 brain MR images received an average relative error in terms of peak signal-to-noise ratio (PSNR) less than 0.1%. In comparison with many state-of-the-art filters, the proposed automation framework with CUDA-based bilateral filtering provided more favorable results both quantitatively and qualitatively. CONCLUSIONS: Possessing unique characteristics and demonstrating exceptional performances, the proposed CUDA-based bilateral filter adequately removed random noise in multifarious brain MR images for further study in neurosciences and radiological sciences. It requires no prior knowledge of the noise variance and automatically restores MR images while preserving fine details. The strategy of exploiting the CUDA to accelerate the computation and incorporating texture features into the BPN to completely automate the bilateral filtering process is achievable and validated, from which the best performance is reached.

Automatic brain MR image denoising based on texture feature-based artificial neural networks.

Noise is one of the main sources of quality deterioration not only for visual inspection but also in computerized processing in brain magnetic resonance (MR) image analysis such as tissue classification, segmentation and registration. Accordingly, noise removal in brain MR images is important for a wide variety of subsequent processing applications. However, most existing denoising algorithms require laborious tuning of parameters that are often sensitive to specific image features and textures. Automation of these parameters through artificial intelligence techniques will be highly beneficial. In the present study, an artificial neural network associated with image texture feature analysis is proposed to establish a predictable parameter model and automate the denoising procedure. In the proposed approach, a total of 83 image attributes were extracted based on four categories: 1) Basic image statistics. 2) Gray-level co-occurrence matrix (GLCM). 3) Gray-level run-length matrix (GLRLM) and 4) Tamura texture features. To obtain the ranking of discrimination in these texture features, a paired-samples t-test was applied to each individual image feature computed in every image. Subsequently, the sequential forward selection (SFS) method was used to select the best texture features according to the ranking of discrimination. The selected optimal features were further incorporated into a back propagation neural network to establish a predictable parameter model. A wide variety of MR images with various scenarios were adopted to evaluate the performance of the proposed framework. Experimental results indicated that this new automation system accurately predicted the bilateral filtering parameters and effectively removed the noise in a number of MR images. Comparing to the manually tuned filtering process, our approach not only produced better denoised results but also saved significant processing time.

Adaptive registration of magnetic resonance images based on a viscous fluid model.

This paper develops a new viscous fluid registration algorithm that makes use of a closed incompressible viscous fluid model associated with mutual information. In our approach, we treat the image pixels as the fluid elements of a viscous fluid governed by the nonlinear Navier-Stokes partial differential equation (PDE) that varies in both temporal and spatial domains. We replace the pressure term with an image-based body force to guide the transformation that is weighted by the mutual information between the template and reference images. A computationally efficient algorithm with staggered grids is introduced to obtain stable solutions of this modified PDE for transformation. The registration process of updating the body force, the velocity and deformation fields is repeated until the mutual information reaches a prescribed threshold. We have evaluated this new algorithm in a number of synthetic and medical images. As consistent with the theory of the viscous fluid model, we found that our method faithfully transformed the template images into the reference images based on the intensity flow. Experimental results indicated that the proposed scheme achieved stable registrations and accurate transformations, which is of potential in large-scale medical image deformation applications.

Automatic segmentation of brain MR images using an adaptive balloon snake model with fuzzy classification.

Skull-stripping in magnetic resonance (MR) images is one of the most important preprocessing steps in medical image analysis. We propose a hybrid skull-stripping algorithm based on an adaptive balloon snake (ABS) model. The proposed framework consists of two phases: first, the fuzzy possibilistic c-means (FPCM) is used for pixel clustering, which provides a labeled image associated with a clean and clear brain boundary. At the second stage, a contour is initialized outside the brain surface based on the FPCM result and evolves under the guidance of an adaptive balloon snake model. The model is designed to drive the contour in the inward normal direction to capture the brain boundary. The entire volume is segmented from the center slice toward both ends slice by slice. Our ABS algorithm was applied to numerous brain MR image data sets and compared with several state-of-the-art methods. Four similarity metrics were used to evaluate the performance of the proposed technique. Experimental results indicated that our method produced accurate segmentation results with higher conformity scores. The effectiveness of the ABS algorithm makes it a promising and potential tool in a wide variety of skull-stripping applications and studies.

Active shape modeling with electric flows.

Physics-based particle systems are an effective tool for shape modeling. Also, there has been much interest in the study of shape modeling using deformable contour approaches. In this paper, we describe a new deformable model with electric flows based upon computer simulations of a number of charged particles embedded in an electrostatic system. Making use of optimized numerical techniques, the electric potential associated with the electric field in the simulated system is rapidly calculated using the finite-size particle (FSP) method. The simulation of deformation evolves based upon the vector sum of two interacting forces: one from the electric fields and the other from the image gradients. Inspired by the concept of the signed distance function associated with the entropy condition in the level set framework, we efficiently handle topological changes at the interface. In addition to automatic splitting and merging, the evolving contours enable simultaneous detection of various objects with varying intensity gradients at both interior and exterior boundaries. This electric flows approach for shape modeling allows one to connect electric properties in electrostatic equilibrium and classical active contours based upon the theory of curve evolution. Our active contours can be applied to model arbitrarily complicated objects including shapes with sharp corners and cusps, and to situations where no a priori knowledge about the object's topology and geometry is made. We demonstrate the capabilities of this new algorithm in recovering a wide variety of structures on simulated and real images in both 2D and 3D.

Performance measure characterization for evaluating neuroimage segmentation algorithms.

Characterizing the performance of segmentation algorithms in brain images has been a persistent challenge due to the complexity of neuroanatomical structures, the quality of imagery and the requirement of accurate segmentation. There has been much interest in using the Jaccard and Dice similarity coefficients associated with Sensitivity and Specificity for evaluating the performance of segmentation algorithms. This paper addresses the essential characteristics of the fundamental performance measure coefficients adopted in evaluation frameworks. While exploring the properties of the Jaccard, Dice and Specificity coefficients, we propose new measure coefficients Conformity and Sensibility for evaluating image segmentation techniques. It is indicated that Conformity is more sensitive and rigorous than Jaccard and Dice in that it has better discrimination capabilities in detecting small variations in segmented images. Comparing to Specificity, Sensibility provides consistent and reliable evaluation scores without the incorporation of image background properties. The merits of the proposed coefficients are illustrated by extracting neuroanatomical structures in a wide variety of brain images using various segmentation techniques.

Computer-assisted extraction of intracranial aneurysms on 3D rotational angiograms for computational fluid dynamics modeling.

PURPOSE: Three-dimensional rotational angiography (3DRA) is an evolving imaging procedure from traditional digital subtraction angiography and is gaining much interest for detecting intracranial aneurysms. Computational fluid dynamics (CFD) modeling plays an important role in understanding the biomechanical properties and in facilitating the prediction of aneurysm rupture. A successful computational study relies on an accurate description of the vascular geometry that is obtained from volumetric images. METHODS: The authors propose a new aneurysm segmentation algorithm to facilitate the study of CFD. This software combines a region-growing segmentation method with the 3D extension of a deformable contour based on a charged fluid model. A charged fluid model essentially consists of a set of charged elements that are governed by the nature of electrostatics. The approach requires no prior knowledge of anatomic structures and automatically segments the vasculature after the end-user selects a vessel section in a plane image. RESULTS: Experimental results on 15 cases indicate that aneurysm structures were effectively segmented and in good agreement with manual delineation outcomes. In comparison with the existing methods, the algorithm provided a much higher overlap index with respect to the ground truth. Furthermore, the outcomes of the proposed approach achieved a clean representation of vascular structures that is advantageous for hemodynamics analyses. CONCLUSIONS: A new aneurysm segmentation framework in an attempt to automatically segment vascular structures in 3DRA image volumes has been developed. The proposed algorithm demonstrated promising performance and unique characteristics to adequately segment aneurysms in 3DRA image volumes for further study in computational fluid dynamics.

Computer simulations of isolated conductors in electrostatic equilibrium.

A computer simulation model is introduced to study the characteristics of isolated conductors in electrostatic equilibrium. Drawing an analogy between electrons and how they move to the surface of isolated conductors, we randomly initialize a large number of particles inside a small region at the center of simulated conductors and advance them according to their forces of repulsion. By use of optimized numerical techniques of the finite-size particle method associated with Poisson's equation, the particles are quickly advanced using a fast Fourier transform and their charge is efficiently shared using the clouds-in-cells method. The particle populations in the simulations range from 50x10;{3} to 1x10;{6} that move in various computation domains equal to 128x128 , 256x256 , and 512x512 grids. When the particles come to an electrostatic equilibrium, they lie on the boundaries of the simulated conductors, from which the equilibrium properties are obtained. Consistent with the theory of electrostatics and charged conductors, we found that the particles move in response to the conductor geometry in such a way that the electrostatic energy is minimized. Good approximation results for the equilibrium properties were obtained using the proposed computer simulation model.

An electrostatic deformable model for medical image segmentation.

A new deformable model, the charged fluid model (CFM), that uses the simulation of charged elements was used to segment medical images. Poisson's equation was used to guide the evolution of the CFM in two steps. In the first step, the elements of the charged fluid were distributed along the propagating interface until electrostatic equilibrium was achieved. In the second step, the propagating front of the charged fluid was deformed in response to the image gradient. The CFM provided sub-pixel precision, required only one parameter setting, and required no prior knowledge of the anatomy of the segmented object. The characteristics of the CFM were compared to existing deformable models using CT and MR images. The results indicate that the CFM is a promising approach for the segmentation of anatomic structures in a wide variety of medical images across different modalities.

Segmentation of brain MR images using a charged fluid model.

In this paper, we developed a new deformable model, the charged fluid model (CFM), that uses the simulation of a charged fluid to segment anatomic structures in magnetic resonance (MR) images of the brain. Conceptually, the charged fluid behaves like a liquid such that it flows through and around different obstacles. The simulation evolves in two steps governed by Poisson's equation. The first step distributes the elements of the charged fluid within the propagating interface until an electrostatic equilibrium is achieved. The second step advances the propagating front of the charged fluid such that it deforms into a new shape in response to the image gradient. This approach required no prior knowledge of anatomic structures, required the use of only one parameter, and provided subpixel precision in the region of interest. We demonstrated the performance of this new algorithm in the segmentation of anatomic structures on simulated and real brain MR images of different subjects. The CFM was compared to the level-set-based methods [Caselles et al. (1993) and Malladi et al (1995)] in segmenting difficult objects in a variety of brain MR images. The experimental results in different types of MR images indicate that the CFM algorithm achieves good segmentation results and is of potential value in brain image processing applications.