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

Automated medical image segmentation techniques

Accurate segmentation of medical images is a key step in contouring during radiotherapy planning. Computed topography (CT) and Magnetic resonance (MR) imaging are the most widely used radiographic techniques in diagnosis, clinical studies and treatment planning. This review provides details of automated segmentation methods, specifically discussed in the context of CT and MR images. The motive is to discuss the problems encountered in segmentation of CT and MR images, and the relative merits and limitations of methods currently available for segmentation of medical images.

Image segmentation using fuzzy region competition and spatial/frequency information

This paper presents a multiphase fuzzy region competition model that takes into account spatial and frequency information for image segmentation. In the proposed energy functional, each region is represented by a fuzzy membership function and a data fidelity term that measures the conformity of spatial and frequency data within each region to (generalized) gaussian densities whose parameters are determined jointly with the segmentation process. Compared with the classical region competition model, our approach gives soft segmentation results via the fuzzy membership functions, and moreover, the use of frequency data provides additional region information that can improve the overall segmentation result. To efficiently solve the minimization of the energy functional, we adopt an alternate minimization procedure and make use of Chambolle's fast duality projection algorithm. We apply the proposed method to synthetic and natural textures as well as real-world natural images. Experimental results show that our proposed method has very promising segmentation performance compared with the current state-of-the-art approaches.

A MR Brain Classification Method Based on Multiscale and Multiblock Fuzzy C-means

A fully automatic, multiscale and multiblock fuzzy C-means (MsbFCM) classification method with intensity correction for MR images is presented in this paper. We use a bilateral filter to process MR images and to build a multiscale image series by increasing the standard deviation of spatial function and reducing the standard deviation of range function. We separate every scale image into multiple blocks and for every block a multiscale fuzzy C-means classification method is applied along the scales from the coarse to fine levels to overcome the effect of intensity inhomogeneity. The method is robust for noise MR images with intensity inhomogeneity because of its multiscale and multiblock bilateral filtering scheme. Our method was compared with the conventional FCM, a modified FCM (MFCM) and multiscale FCM (MsFCM) method on synthesized images, simulated brain MR images, and real MR images. The MsbFCM method achieved an overlap ratio of greater than 91% as validated by the ground truth even if original images have 9% noise and 40% intensity inhomogeneity. Experimental results using real MR images demonstrate the effectiveness of the proposed method. Our MsbFCM classification method is accurate and robust for various MR images.

A modified possibilistic fuzzy c-means clustering algorithm for bias field estimation and segmentation of brain MR image

A modified possibilistic fuzzy c-means clustering algorithm is presented for fuzzy segmentation of magnetic resonance (MR) images that have been corrupted by intensity inhomogeneities and noise. By introducing a novel adaptive method to compute the weights of local spatial in the objective function, the new adaptive fuzzy clustering algorithm is capable of utilizing local contextual information to impose local spatial continuity, thus allowing the suppression of noise and helping to resolve classification ambiguity. To estimate the intensity inhomogeneity, the global intensity is introduced into the coherent local intensity clustering algorithm and takes the local and global intensity information into account. The segmentation target therefore is driven by two forces to smooth the derived optimal bias field and improve the accuracy of the segmentation task. The proposed method has been successfully applied to 3 T, 7 T, synthetic and real MR images with desirable results. Comparisons with other approaches demonstrate the superior performance of the proposed algorithm. Moreover, the proposed algorithm is robust to initialization, thereby allowing fully automatic applications.

Intensity inhomogeneity correction for the breast sonogram: constrained fuzzy cell-based bipartitioning and polynomial surface modeling

PURPOSE

To develop an intensity inhomogeneity algorithm for breast sonograms in order to assist visual identification and automatic delineation of lesion boundaries.

METHODS

The proposed algorithm was composed of two essential ideas. One was decomposing the region of interest (ROI) into foreground and background regions by a cell-based segmentation algorithm, called constrained fuzzy cell-based bipartition-EM (CFCB-EM) algorithm. The CFCB-EM algorithm deformed the contour in a fuzzy cell-based deformation fashion with the cell structures derived by the fuzzy cell competition (FCC) algorithm as the deformation unit and the boundary estimated by the normalized cut (NC) algorithm as the reference contour. The other was modeling the intensity inhomogeneity in an ROI as a spatially variant normal distribution with a constant variance and spatially variant means, which formed a polynomial surface of order n. The proposed algorithm was formulated as a nested EM algorithm comprising the outer-layer EM algorithm, i.e., the intensity inhomogeneity correction-EM (IIC-EM) algorithm, and the inner-layer EM algorithm, i.e., the CFCB-EM algorithm. The E step of the IIC-EM algorithm was to provide a reasonably good bipartition separating the ROI into foreground and background regions, which included three major component algorithms, namely, the FCC, the NC, and the CFCB-EM. The M step of the IIC-EM algorithm was to estimate and correct the intensity inhomogeneity field by least-squared fitting the intensity inhomogeneity to an nth order polynomial surface. Forty-nine breast sonograms with intensity inhomogeneity, each from a different subject, were randomly selected for performance analysis. Three assessments were carried out to evaluate the effectiveness of the proposed algorithm.

RESULTS

Based on the visual evaluation of two experienced radiologists, in the first assessment, 46 out of 49 breast lesions were considered to have better contrasts on the inhomogeneity-corrected images by both radiologists. The interrater reliability for the radiologists was found to be kappa = 0.479 (p = 0.001). In the second assessment, the mean gradients of the low-gradient boundary points before and after correction of the intensity inhomogeneity were compared by the paired t-test, yielding a p-value of 0.000, which suggested the proposed intensity inhomogeneity algorithm may enhance the mean gradient of the low-gradient boundary points. By using the paired t-test, the third assessment further showed that the Chan and Vese level set method could derive a much better lesion boundary on the inhomogeneity-corrected image than on the original image (p = 0.000).

CONCLUSIONS

The proposed intensity inhomogeneity correction algorithm could not only augment the visibility of lesion boundary but also improve the segmentation result on a breast sonogram.

Fast algorithm for calculation of inhomogeneity gradient in magnetic resonance imaging data

PURPOSE

To develop and implement a new approach for correcting the intensity inhomogeneity in magnetic resonance imaging (MRI) data.

MATERIALS AND METHODS

The algorithm is based on the assumption that intensity inhomogeneity in MR data is multiplicative and smoothly varying. Using a statistically stable method, the algorithm first calculates the partial derivative of the inhomogeneity gradient across the data. The algorithm then solves for the gradient field and fits it to a parametric surface. It was tested on both simulated and real human and animal MRI data.

RESULTS

The algorithm is shown to restore the homogeneity in all images that were tested. On real human brain images the algorithm demonstrated superior or comparable performance relative to some of the commonly used intensity inhomogeneity correction methods such as SPM, BrainSuite, and N3.

CONCLUSION

The proposed algorithm provides an alternative method for correcting the intensity inhomogeneity in MR images. It is shown to be fast and its performance is superior or comparable to algorithms described in the published literature. Due to its generality, this algorithm is applicable to MR images of both humans and animals.

Development of a virtual reality guided diagnostic tool based on magnetic resonance imaging

Computed tomography (CT) and virtual reality (VR) made it possible to create internal views of the human body without actual penetration. During the last two decades, several endoscopic diagnosis procedures have received virtual counter candidates. This paper presents an own concept of a virtual reality guided diagnostic tool, based on magnetic resonance images representing parallel cross-sections of the investigated organ. A series of image processing methods are proposed for image quality enhancement, accurate segmentation in two dimensions, and three-dimensional reconstruction of detected surfaces. These techniques provide improved accuracy in image segmentation, and thus they represent excellent support for three dimensional imaging. The implemented software system allows interactive navigation within the investigated volume, and provides several facilities to quantify important physical properties including distances, areas, and volumes.

Computational simulation of breast compression based on segmented breast and fibroglandular tissues on magnetic resonance images

This study presents a finite element-based computational model to simulate the three-dimensional deformation of a breast and fibroglandular tissues under compression. The simulation was based on 3D MR images of the breast, and craniocaudal and mediolateral oblique compression, as used in mammography, was applied. The geometry of the whole breast and the segmented fibroglandular tissues within the breast were reconstructed using triangular meshes by using the Avizo 6.0 software package. Due to the large deformation in breast compression, a finite element model was used to simulate the nonlinear elastic tissue deformation under compression, using the MSC.Marc software package. The model was tested in four cases. The results showed a higher displacement along the compression direction compared to the other two directions. The compressed breast thickness in these four cases at a compression ratio of 60% was in the range of 5-7 cm, which is a typical range of thickness in mammography. The projection of the fibroglandular tissue mesh at a compression ratio of 60% was compared to the corresponding mammograms of two women, and they demonstrated spatially matched distributions. However, since the compression was based on magnetic resonance imaging (MRI), which has much coarser spatial resolution than the in-plane resolution of mammography, this method is unlikely to generate a synthetic mammogram close to the clinical quality. Whether this model may be used to understand the technical factors that may impact the variations in breast density needs further investigation. Since this method can be applied to simulate compression of the breast at different views and different compression levels, another possible application is to provide a tool for comparing breast images acquired using different imaging modalities--such as MRI, mammography, whole breast ultrasound and molecular imaging--that are performed using different body positions and under different compression conditions.

A novel fuzzy C-means algorithm for unsupervised heterogeneous tumor quantification in PET

PURPOSE

Accurate and robust image segmentation was identified as one of the most challenging issues facing PET quantification in oncological imaging. This difficulty is compounded by the low spatial resolution and high noise characteristics of PET images. The fuzzy C-means (FCM) clustering algorithm was largely used in various medical image segmentation approaches. However, the algorithm is sensitive to both noise and intensity heterogeneity since it does not take into account spatial contextual information.

METHODS

To overcome this limitation, a new fuzzy segmentation technique adapted to typical noisy and low resolution oncological PET data is proposed. PET images smoothed using a nonlinear anisotropic diffusion filter are added as a second input to the proposed FCM algorithm to incorporate spatial information (FCM-S). In addition, a methodology was developed to integrate the a trous wavelet transform in the standard FCM algorithm (FCM-SW) to allow handling of heterogeneous lesions' uptake. The algorithm was applied to the simulated data of the NCAT phantom, incorporating heterogeneous lesions in the lung and clinical PET/CT images of 21 patients presenting with histologically proven nonsmall-cell lung cancer (NSCLC) and 7 patients presenting with laryngeal squamous cell carcinoma (LSCC) to assess its performance for segmenting tumors with arbitrary size, shape, and tracer uptake. For NSCLC patients, the maximal tumor diameters measured from the macroscopic examination of the surgical specimen served as the ground truth for comparison with the maximum diameter estimated by the segmentation technique, whereas for LSCC patients, the 3D macroscopic tumor volume was considered as the ground truth for comparison with the corresponding PET-based volume. The proposed algorithm was also compared to the classical FCM segmentation technique.

RESULTS

There is a good correlation (R2 = 0.942) between the actual maximal diameter of primary NSCLC tumors estimated using the proposed PET segmentation procedure and those measured from the macroscopic examination, and the regression line agreed well with the line of identity (slope = 1.08) for the group analysis of the clinical data. The standard FCM algorithm seems to underestimate actual maximal diameters of the clinical data, resulting in a mean error of -4.6 mm (relative error of -10.8 +/- 23.1%) for all data sets. The mean error of maximal diameter estimation was reduced to 0.1 mm (0.9 +/- 14.4%) using the proposed FCM-SW algorithm. Likewise, the mean relative error on the estimated volume for LSCC patients was reduced from 21.7 +/- 22.0% for FCM to 8.6 +/- 28.3% using the proposed FCM-SW technique.

CONCLUSIONS

A novel unsupervised PET image segmentation technique that allows the quantification of lesions in the presence of heterogeneity of tracer uptake was developed and evaluated. The technique is being further refined and assessed in clinical setting to delineate treatment volumes for the purpose of PET-guided radiation therapy treatment planning but could find other applications in clinical oncology such as the assessment of response to treatment.

Effective fuzzy c-means based kernel function in segmenting medical images

The objective of this paper is to develop an effective robust fuzzy c-means for a segmentation of breast and brain magnetic resonance images. The widely used conventional fuzzy c-means for medical image segmentations has limitations because of its squared-norm distance measure to measure the similarity between centers and data objects of medical images which are corrupted by heavy noise, outliers, and other imaging artifacts. To overcome the limitations this paper develops a novel objective function based standard objective function of fuzzy c-means that incorporates the robust kernel-induced distance for clustering the corrupted dataset of breast and brain medical images. By minimizing the novel objective function this paper obtains effective equation for optimal cluster centers and equation to achieve optimal membership grades for partitioning the given dataset. In order to solve the problems of clustering performance affected by initial centers of clusters, this paper introduces a specialized center initialization method for executing the proposed algorithm in segmenting medical images. Experiments are performed with synthetic, real breast and brain images to assess the performance of the proposed method. Further the validity of clustering results is obtained using silhouette method and this paper compares the results with the results of other recent reported fuzzy c-means methods. The experimental results show the superiority of the proposed clustering results.

A robust fuzzy local information C-Means clustering algorithm

This paper presents a variation of fuzzy c-means (FCM) algorithm that provides image clustering. The proposed algorithm incorporates the local spatial information and gray level information in a novel fuzzy way. The new algorithm is called fuzzy local information C-Means (FLICM). FLICM can overcome the disadvantages of the known fuzzy c-means algorithms and at the same time enhances the clustering performance. The major characteristic of FLICM is the use of a fuzzy local (both spatial and gray level) similarity measure, aiming to guarantee noise insensitiveness and image detail preservation. Furthermore, the proposed algorithm is fully free of the empirically adjusted parameters (a, ¿(g), ¿(s), etc.) incorporated into all other fuzzy c-means algorithms proposed in the literature. Experiments performed on synthetic and real-world images show that FLICM algorithm is effective and efficient, providing robustness to noisy images.

Local bone enhancement fuzzy clustering for segmentation of MR trabecular bone images

PURPOSE

Segmentation of trabecular bone from magnetic resonance (MR) images is a challenging task due to spatial resolution limitations, signal-to-noise ratio constraints, and signal intensity inhomogeneities. This article examines an alternative approach to trabecular bone segmentation using partial membership segmentation termed fuzzy C-means clustering incorporating local second order features for bone enhancement (BE-FCM) at multiple scales. This approach is meant to allow for a soft segmentation that accounts for partial volume effects while suppressing the influence of noise.

METHODS

A soft segmentation method was developed and evaluated on three different sets of data; interscan reproducibility was evaluated on six test-retest in vivo MR scans of the proximal femur, correlation between MR and HR-pQCT measurements was evaluated on 49 in vivo scans from the distal tibia, and the potential for fracture discrimination was evaluated using MR scans of calcaneus specimens from 15 participants with and 15 participants without vertebral fracture. The algorithm was compared to fuzzy clustering using the intensity as the only feature (I-FCM) and a dual thresholding algorithm. The metric evaluated was bone volume over total volume (BV/TV) within user-defined regions of interest.

RESULTS

BE-FCM had a higher interscan reproducibility (rms CV: 2.0%) compared to I-FCM (5.6%) and thresholding (4.2%), and expressed higher correlation to HR-pQCT data (r = 0.79, p < 10(-11)) compared to I-FCM (r = 0.74, p < 10(-8)) and thresholding (r = 0.70, p < 10(-6)). BE-FCM was also the method that was best able to differentiate between a control and a vertebral fracture group at a 95% significance level.

CONCLUSIONS

The results suggest that trabecular bone segmentation by BE-FCM can provide a precise BV/TV measurement that is sensitive to pathology. The segmentation method may become useful in MR imaging-based quantification of bone microarchitecture.

Automated classification of multi-spectral MR images using Linear Discriminant Analysis

Magnetic resonance imaging (MRI) is a valuable instrument in medical science owing to its capabilities in soft tissue characterization and 3D visualization. A potential application of MRI in clinical practice is brain parenchyma classification. This work proposes a novel approach called "Unsupervised Linear Discriminant Analysis (ULDA)" to classify and segment the three major tissues, i.e. gray matter (GM), white matter (WM) and cerebral spinal fluid (CSF), from a multi-spectral MR image of the human brain. The ULDA comprises two processes, namely Target Generation Process (TGP) and Linear Discriminant Analysis (LDA) classification. TGP is a fuzzy-set process that generates a set of potential targets from unknown information, and applies these targets to train the optimal division boundary by LDA, such that three tissues GM, WM and CSF are separated. Finally, two sets of images, namely computer-generated phantom images and real MR images are used in the experiments to evaluate the effectiveness of ULDA. Experiment results reveal that UDLA segments a multi-spectral MR image much more effectively than either FMRIB's Automated Segmentation Tool (FAST) or Fuzzy C-means (FC).

Pattern Recognition of Longitudinal Trial Data with Nonignorable Missingness: An Empirical Case Study

Methods for identifying meaningful growth patterns of longitudinal trial data with both nonignorable intermittent and drop-out missingness are rare. In this study, a combined approach with statistical and data mining techniques is utilized to address the nonignorable missing data issue in growth pattern recognition. First, a parallel mixture model is proposed to model the nonignorable missing information from a real-world patient-oriented study and concurrently to estimate the growth trajectories of participants. Then, based on individual growth parameter estimates and their auxiliary feature attributes, a fuzzy clustering method is incorporated to identify the growth patterns. This case study demonstrates that the combined multi-step approach can achieve both statistical gener ality and computational efficiency for growth pattern recognition in longitudinal studies with nonignorable missing data.

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.

A robust hidden Markov Gauss mixture vector quantizer for a noisy source

Noise is ubiquitous in real life and changes image acquisition, communication, and processing characteristics in an uncontrolled manner. Gaussian noise and Salt and Pepper noise, in particular, are prevalent in noisy communication channels, camera and scanner sensors, and medical MRI images. It is not unusual for highly sophisticated image processing algorithms developed for clean images to malfunction when used on noisy images. For example, hidden Markov Gauss mixture models (HMGMM) have been shown to perform well in image segmentation applications, but they are quite sensitive to image noise. We propose a modified HMGMM procedure specifically designed to improve performance in the presence of noise. The key feature of the proposed procedure is the adjustment of covariance matrices in Gauss mixture vector quantizer codebooks to minimize an overall minimum discrimination information distortion (MDI). In adjusting covariance matrices, we expand or shrink their elements based on the noisy image. While most results reported in the literature assume a particular noise type, we propose a framework without assuming particular noise characteristics. Without denoising the corrupted source, we apply our method directly to the segmentation of noisy sources. We apply the proposed procedure to the segmentation of aerial images with Salt and Pepper noise and with independent Gaussian noise, and we compare our results with those of the median filter restoration method and the blind deconvolution-based method, respectively. We show that our procedure has better performance than image restoration-based techniques and closely matches to the performance of HMGMM for clean images in terms of both visual segmentation results and error rate.

Improvement of brain segmentation accuracy by optimizing non-uniformity correction using N3

Smoothly varying and multiplicative intensity variations within MR images that are artifactual, can reduce the accuracy of automated brain segmentation. Fortunately, these can be corrected. Among existing correction approaches, the nonparametric non-uniformity intensity normalization method N3 (Sled, J.G., Zijdenbos, A.P., Evans, A.C., 1998. Nonparametric method for automatic correction of intensity nonuniformity in MRI data. IEEE Trans. Med. Imag. 17, 87-97.) is one of the most frequently used. However, at least one recent study (Boyes, R.G., Gunter, J.L., Frost, C., Janke, A.L., Yeatman, T., Hill, D.L.G., Bernstein, M.A., Thompson, P.M., Weiner, M.W., Schuff, N., Alexander, G.E., Killiany, R.J., DeCarli, C., Jack, C.R., Fox, N.C., 2008. Intensity non-uniformity correction using N3 on 3-T scanners with multichannel phased array coils. NeuroImage 39, 1752-1762.) suggests that its performance on 3 T scanners with multichannel phased-array receiver coils can be improved by optimizing a parameter that controls the smoothness of the estimated bias field. The present study not only confirms this finding, but additionally demonstrates the benefit of reducing the relevant parameter values to 30-50 mm (default value is 200 mm), on white matter surface estimation as well as the measurement of cortical and subcortical structures using FreeSurfer (Martinos Imaging Centre, Boston, MA). This finding can help enhance precision in studies where estimation of cerebral cortex thickness is critical for making inferences.

MR image segmentation using a power transformation approach

This study proposes a segmentation method for brain MR images using a distribution transformation approach. The method extends traditional Gaussian mixtures expectation-maximization segmentation to a power transformed version of mixed intensity distributions, which includes Gaussian mixtures as a special case. As MR intensities tend to exhibit non-Gaussianity due to partial volume effects, the proposed method is designed to fit non-Gaussian tissue intensity distributions. One advantage of the method is that it is intuitively appealing and computationally simple. To avoid performance degradation caused by intensity inhomogeneity, different methods for correcting bias fields were applied prior to image segmentation, and their correction effects on the segmentation results were examined in the empirical study. The partitions of brain tissues (i.e., gray and white matter) resulting from the method were validated and evaluated against manual segmentation results based on 38 real T1-weighted image volumes from the internet brain segmentation repository, and 18 simulated image volumes from BrainWeb. The Jaccard and Dice similarity indexes were computed to evaluate the performance of the proposed approach relative to the expert segmentations. Empirical results suggested that the proposed segmentation method yielded higher similarity measures for both gray matter and white matter as compared with those based on the traditional segmentation using the Gaussian mixtures approach.

Nonrigid registration of joint histograms for intensity standardization in magnetic resonance imaging

A major disadvantage of magnetic resonance imaging (MRI) compared to other imaging modalities like computed tomography is the fact that its intensities are not standardized. Our contribution is a novel method for MRI signal intensity standardization of arbitrary MRI scans, so as to create a pulse sequence dependent standard intensity scale. The proposed method is the first approach that uses the properties of all acquired images jointly (e.g., T1- and T2-weighted images). The image properties are stored in multidimensional joint histograms. In order to normalize the probability density function (pdf) of a newly acquired data set, a nonrigid image registration is performed between a reference and the joint histogram of the acquired images. From this matching a nonparametric transformation is obtained, which describes a mapping between the corresponding intensity spaces and subsequently adapts the image properties of the newly acquired series to a given standard. As the proposed intensity standardization is based on the probability density functions of the data sets only, it is independent of spatial coherence or prior segmentations of the reference and current images. Furthermore, it is not designed for a particular application, body region or acquisition protocol. The evaluation was done using two different settings. First, MRI head images were used, hence the approach can be compared to state-of-the-art methods. Second, whole body MRI scans were used. For this modality no other normalization algorithm is known in literature. The Jeffrey divergence of the pdfs of the whole body scans was reduced by 45%. All used data sets were acquired during clinical routine and thus included pathologies.

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.

MRI tissue classification and bias field estimation based on coherent local intensity clustering: a unified energy minimization framework

This paper presents a new energy minimization method for simultaneous tissue classification and bias field estimation of magnetic resonance (MR) images. We first derive an important characteristic of local image intensities--the intensities of different tissues within a neighborhood form separable clusters, and the center of each cluster can be well approximated by the product of the bias within the neighborhood and a tissue-dependent constant. We then introduce a coherent local intensity clustering (CLIC) criterion function as a metric to evaluate tissue classification and bias field estimation. An integration of this metric defines an energy on a bias field, membership functions of the tissues, and the parameters that approximate the true signal from the corresponding tissues. Thus, tissue classification and bias field estimation are simultaneously achieved by minimizing this energy. The smoothness of the derived optimal bias field is ensured by the spatially coherent nature of the CLIC criterion function. As a result, no extra effort is needed to smooth the bias field in our method. Moreover, the proposed algorithm is robust to the choice of initial conditions, thereby allowing fully automatic applications. Our algorithm has been applied to high field and ultra high field MR images with promising results.

An artificial ant colonies approach to medical image segmentation

The success of image analysis depends heavily upon accurate image segmentation algorithms. This paper presents a novel segmentation algorithm based on artificial ant colonies (AC). Recent studies show that the self-organization of ants is similar to neurons in the human brain in many respects. Therefore, it has been used successfully for understanding biological systems. It is also widely used in many applications in robotics, computer graphics, etc. Considering the features of artificial ant colonies, we present an extended model for image segmentation. In our model, each ant can memorize a reference object, which will be refreshed when it finds a new target. A fuzzy connectedness measure is adopted to evaluate the similarity between target and the reference object. The behavior of an ant is affected by the neighbors and the cooperation between ants is performed by exchanging information through pheromone updating. Experimental results show that the new algorithm can preserve the detail of the object and is also insensitive to noise.

Segmentation of brain magnetic resonance images for measurement of gray matter atrophy in multiple sclerosis patients

Multiple sclerosis (MS) affects both white matter and gray matter (GM). Measurement of GM volumes is a particularly useful method to estimate the total extent of GM tissue damage because it can be done with conventional magnetic resonance images (MRI). Many algorithms exist for segmentation of GM, but none were specifically designed to handle issues associated with MS, such as atrophy and the effects that MS lesions may have on the classification of GM. A new GM segmentation algorithm has been developed specifically for calculation of GM volumes in MS patients. The new algorithm uses a combination of intensity, anatomical, and morphological probability maps. Several validation tests were performed to evaluate the algorithm in terms of accuracy, reproducibility, and sensitivity to MS lesions. The accuracy tests resulted in error rates of 1.2% and 3.1% for comparisons to BrainWeb and manual tracings, respectively. Similarity indices indicated excellent agreement with the BrainWeb segmentation (0.858-0.975, for various levels of noise and rf inhomogeneity). The scan-rescan reproducibility test resulted in a mean coefficient of variation of 1.1% for GM fraction. Tests of the effects of varying the size of MS lesions revealed a moderate and consistent dependence of GM volumes on T2 lesion volume, which suggests that GM volumes should be corrected for T2 lesion volumes using a simple scale factor in order to eliminate this technical artifact. The new segmentation algorithm can be used for improved measurement of GM volumes in MS patients, and is particularly applicable to retrospective datasets.

Bias artifact suppression on MR volumes

RF-inhomogeneity correction is a relevant research topic in the field of magnetic resonance imaging (MRI). A volume corrupted by this artifact exhibits nonuniform illumination both inside a single slice and between adjacent ones. In this work a bias correction technique is presented, which suppresses this artifact on MR volumes scanned from different body parts without any a priori hypothesis on the artifact model. Theoretical foundations of the method are reported together with experimental results and a comparison is presented with both the 2D version of the algorithm and other techniques that are widely used in MRI literature.