Automatic section thickness determination using an absolute gradient focus function.

Quantitative analysis of microstructures using computerized stereology systems is an essential tool in many disciplines of bioscience research. Section thickness determination in current nonautomated approaches requires manual location of upper and lower surfaces of tissue sections. In contrast to conventional autofocus functions that locate the optimally focused optical plane using the global maximum on a focus curve, this study identified by two sharp 'knees' on the focus curve as the transition from unfocused to focused optical planes. Analysis of 14 grey-scale focus functions showed, the thresholded absolute gradient function, was best for finding detectable bends that closely correspond to the bounding optical planes at the upper and lower tissue surfaces. Modifications to this function generated four novel functions that outperformed the original. The 'modified absolute gradient count' function outperformed all others with an average error of 0.56 µm on a test set of images similar to the training set; and, an average error of 0.39 µm on a test set comprised of images captured from a different case, that is, different staining methods on a different brain region from a different subject rat. We describe a novel algorithm that allows for automatic section thickness determination based on just out-of-focus planes, a prerequisite for fully automatic computerized stereology.

Comprehensive processing, display and analysis for in vivo MR spectroscopic imaging.

Image reconstruction for magnetic resonance spectroscopic imaging (MRSI) requires specialized spatial and spectral data processing methods and benefits from the use of several sources of prior information that are not commonly available, including MRI-derived tissue segmentation, morphological analysis and spectral characteristics of the observed metabolites. In addition, incorporating information obtained from MRI data can enhance the display of low-resolution metabolite images and multiparametric and regional statistical analysis methods can improve detection of altered metabolite distributions. As a result, full MRSI processing and analysis can involve multiple processing steps and several different data types. In this paper, a processing environment is described that integrates and automates these data processing and analysis functions for imaging of proton metabolite distributions in the normal human brain. The capabilities include normalization of metabolite signal intensities and transformation into a common spatial reference frame, thereby allowing the formation of a database of MR-measured human metabolite values as a function of acquisition, spatial and subject parameters. This development is carried out under the MIDAS project (Metabolite Imaging and Data Analysis System), which provides an integrated set of MRI and MRSI processing functions. It is anticipated that further development and distribution of these capabilities will facilitate more widespread use of MRSI for diagnostic imaging, encourage the development of standardized MRSI acquisition, processing and analysis methods and enable improved mapping of metabolite distributions in the human brain.

A generic knowledge-guided image segmentation and labeling system using fuzzy clustering algorithms.

Segmentation of an image into regions and the labeling of the regions is a challenging problem. In this paper, an approach that is applicable to any set of multifeature images of the same location is derived. Our approach applies to, for example, medical images of a region of the body; repeated camera images of the same area; and satellite images of a region. The segmentation and labeling approach described here uses a set of training images and domain knowledge to produce an image segmentation system that can be used without change on images of the same region collected over time. How to obtain training images, integrate domain knowledge, and utilize learning to segment and label images of the same region taken under any condition for which a training image exists is detailed. It is shown that clustering in conjunction with image processing techniques utilizing an iterative approach can effectively identify objects of interest in images. The segmentation and labeling approach described here is applied to color camera images and two other image domains are used to illustrate the applicability of the approach.

Automatic segmentation of non-enhancing brain tumors in magnetic resonance images.

Tumor segmentation from magnetic resonance (MR) images may aid in tumor treatment by tracking the progress of tumor growth and/or shrinkage. In this paper we present the first automatic segmentation method which separates non-enhancing brain tumors from healthy tissues in MR images to aid in the task of tracking tumor size over time. The MR feature images used for the segmentation consist of three weighted images (T1, T2 and proton density (PD)) for each axial slice through the head. An initial segmentation is computed using an unsupervised fuzzy clustering algorithm. Then, integrated domain knowledge and image processing techniques contribute to the final tumor segmentation. They are applied under the control of a knowledge-based system. The system knowledge was acquired by training on two patient volumes (14 images). Testing has shown successful tumor segmentations on four patient volumes (31 images). Our results show that we detected all six non-enhancing brain tumors, located tumor tissue in 35 of the 36 ground truth (radiologist labeled) slices containing tumor and successfully separated tumor regions from physically connected CSF regions in all the nine slices. Quantitative measurements are promising as correspondence ratios between ground truth and segmented tumor regions ranged between 0.368 and 0.871 per volume, with percent match ranging between 0.530 and 0.909 per volume.

Report of research activities in fuzzy AI and medicine at USF CSE.

Several projects involving the use of fuzzy and neuro-fuzzy methods in medical applications, developed by members of the Department of Computer Science and Engineering, University of South Florida, Tampa, Florida, are briefly reviewed. The successful applications are emphasized.

Feature extraction for MRI segmentation.

Magnetic resonance images (MRIs) of the brain are segmented to measure the efficacy of treatment strategies for brain tumors. To date, no reproducible technique for measuring tumor size is available to the clinician, which hampers progress of the search for good treatment protocols. Many segmentation techniques have been proposed, but the representation (features) of the MRI data has received little attention. A genetic algorithm (GA) search was used to discover a feature set from multi-spectral MRI data. Segmentations were performed using the fuzzy c-means (FCM) clustering technique. Seventeen MRI data sets from five patients were evaluated. The GA feature set produces a more accurate segmentation. The GA fitness function that achieves the best results is the Wilks's lambda statistic when applied to FCM clusters. Compared to linear discriminant analysis, which requires class labels, the same or better accuracy is obtained by the features constructed from a GA search without class labels, allowing fully operator independent segmentation. The GA approach therefore provides a better starting point for the measurement of the response of a brain tumor to treatment.

Automatic tumor segmentation using knowledge-based techniques.

A system that automatically segments and labels glioblastoma-multiforme tumors in magnetic resonance images (MRI's) of the human brain is presented. The MRI's consist of T1-weighted, proton density, and T2-weighted feature images and are processed by a system which integrates knowledge-based (KB) techniques with multispectral analysis. Initial segmentation is performed by an unsupervised clustering algorithm. The segmented image, along with cluster centers for each class are provided to a rule-based expert system which extracts the intracranial region. Multispectral histogram analysis separates suspected tumor from the rest of the intracranial region, with region analysis used in performing the final tumor labeling. This system has been trained on three volume data sets and tested on thirteen unseen volume data sets acquired from a single MRI system. The KB tumor segmentation was compared with supervised, radiologist-labeled "ground truth" tumor volumes and supervised k-nearest neighbors tumor segmentations. The results of this system generally correspond well to ground truth, both on a per slice basis and more importantly in tracking total tumor volume during treatment over time.

Medical image analysis with fuzzy models.

This paper updates several recent surveys on the use of fuzzy models for segmentation and edge detection in medical image data. Our survey is divided into methods based on supervised and unsupervised learning (that is, on whether there are or are not labelled data available for supervising the computations), and is organized first and foremost by groups (that we know of!) that are active in this area. Our review is aimed more towards 'who is doing it' rather than 'how good it is'. This is partially dictated by the fact that direct comparisons of supervised and unsupervised methods is somewhat akin to comparing apples and oranges. There is a further subdivision into methods for two- and three-dimensional data and/or problems. We do not cover methods based on neural-like networks or fuzzy reasoning systems. These topics are covered in a recently published companion survey by keller et al.

Normal brain volume measurements using multispectral MRI segmentation.

The performance of a supervised k-nearest neighbor (kNN) classifier and a semisupervised fuzzy c-means (SFCM) clustering segmentation method are evaluated for reproducible measurement of the volumes of normal brain tissues and cerebrospinal fluid. The stability of the two segmentation methods is evaluated for (a) operator selection of training data, (b) reproducibility during repeat imaging sessions to determine any variations in the sensor performance over time, (c) variations in the measured volumes between different subjects, and (d) variability with different imaging parameters. The variations were found to be dependent on the type of measured tissue and the operator performing the segmentations. The variability during repeat imaging sessions for the SFCM method was < 3%. The absolute volumes of the brain matter and cerebrospinal fluid between subjects varied quite large, ranging from 9% to 13%. The intraobserver and interobserver reproducibility for SFCM were < 4% for the soft tissues and 6% for cerebrospinal fluid. The corresponding results for the kNN segmentation method were higher compared to the SFCM method.

Monitoring brain tumor response to therapy using MRI segmentation.

The performance evaluation of a semi-supervised fuzzy c-means (SFCM) clustering method for monitoring brain tumor volume changes during the course of routine clinical radiation-therapeutic and chemo-therapeutic regimens is presented. The tumor volume determined using the SFCM method was compared with the volume estimates obtained using three other methods: (a) a k nearest neighbor (kNN) classifier, b) a grey level thresholding and seed growing (ISG-SG) method and c) a manual pixel labeling (GT) method for ground truth estimation. The SFCM and kNN methods are applied to the multispectral, contrast enhanced T1, proton density, and T2 weighted, magnetic resonance images (MRI) whereas the ISG-SG and GT methods are applied only to the contrast enhanced T1 weighted image. Estimations of tumor volume were made on eight patient cases with follow-up MRI scans performed over a 32 week interval during treatment. The tumor cases studied include one meningioma, two brain metastases and five gliomas. Comparisons with manually labeled ground truth estimations showed that there is a limited agreement between the segmentation methods for absolute tumor volume measurements when using images of patients after treatment. The average intraobserver reproducibility for the SFCM, kNN and ISG-SG methods was found to be 5.8%, 6.6% and 8.9%, respectively. The average of the interobserver reproducibility of these methods was found to be 5.5%, 6.5% and 11.4%, respectively. For the measurement of relative change of tumor volume as required for the response assessment, the multi-spectral methods kNN and SFCM are therefore preferred over the seedgrowing method.

Unsupervised measurement of brain tumor volume on MR images.

We examined unsupervised methods of segmentation of MR images of the brain for measuring tumor volume in response to treatment. Two clustering methods were used: fuzzy c-means and a nonfuzzy clustering algorithm. Results were compared with volume segmentations by two supervised methods, k-nearest neighbors and region growing, and all results were compared with manual labelings. Results of individual segmentations are presented as well as comparisons on the application of the different methods with 10 data sets of patients with brain tumors. Unsupervised segmentation is preferred for measuring tumor volumes in response to treatment, as it eliminates operator dependency and may be adequate for delineation of the target volume in radiation therapy. Some obstacles need to be overcome, in particular regarding the detection of anatomically relevant tissue classes. This study shows that these improvements are possible.

Comparison of supervised MRI segmentation methods for tumor volume determination during therapy.

Two different multispectral pattern recognition methods are used to segment magnetic resonance images (MRI) of the brain for quantitative estimation of tumor volume and volume changes with therapy. A supervised k-nearest neighbor (kNN) rule and a semi-supervised fuzzy c-means (SFCM) method are used to segment MRI slice data. Tumor volumes as determined by the kNN and SFCM segmentation methods are compared with two reference methods, based on image grey scale, as a basis for an estimation of ground truth, namely: (a) a commonly used seed growing method that is applied to the contrast enhanced T1-weighted image, and (b) a manual segmentation method using a custom-designed graphical user interface applied to the same raw image (T1-weighted) dataset. Emphasis is placed on measurement of intra and inter observer reproducibility using the proposed methods. Intra- and interobserver variation for the kNN method was 9% and 5%, respectively. The results for the SFCM method was a little better at 6% and 4%, respectively. For the seed growing method, the intra-observer variation was 6% and the interobserver variation was 17%, significantly larger when compared with the multispectral methods. The absolute tumor volume determined by the multispectral segmentation methods was consistently smaller than that observed for the reference methods. The results of this study are found to be very patient case-dependent. The results for SFCM suggest that it should be useful for relative measurements of tumor volume during therapy, but further studies are required. This work demonstrates the need for minimally supervised or unsupervised methods for tumor volume measurements.

MRI segmentation: methods and applications.

The current literature on MRI segmentation methods is reviewed. Particular emphasis is placed on the relative merits of single image versus multispectral segmentation, and supervised versus unsupervised segmentation methods. Image pre-processing and registration are discussed, as well as methods of validation. The application of MRI segmentation for tumor volume measurements during the course of therapy is presented here as an example, illustrating problems associated with inter- and intra-observer variations inherent to supervised methods.

Application of fuzzy c-means segmentation technique for tissue differentiation in MR images of a hemorrhagic glioblastoma multiforme.

The application of a raw data-based, operator-independent MR segmentation technique to differentiate boundaries of tumor from edema or hemorrhage is demonstrated. A case of a glioblastoma multiforme with gross and histopathologic correlation is presented. The MR image data set was segmented into tissue classes based on three different MR weighted image parameters (T1-, proton density-, and T2-weighted) using unsupervised fuzzy c-means (FCM) clustering algorithm technique for pattern recognition. A radiological examination of the MR images and correlation with fuzzy clustering segmentations was performed. Results were confirmed by gross and histopathology which, to the best of our knowledge, reports the first application of this demanding approach. Based on the results of neuropathologic correlation, the application of FCM MR image segmentation to several MR images of a glioblastoma multiforme represents a viable technique for displaying diagnostically relevant tissue contrast information used in 3D volume reconstruction. With this technique, it is possible to generate segmentation images that display clinically important neuroanatomic and neuropathologic tissue contrast information from raw MR image data.

Review of MR image segmentation techniques using pattern recognition.

This paper has reviewed, with somewhat variable coverage, the nine MR image segmentation techniques itemized in Table II. A wide array of approaches have been discussed; each has its merits and drawbacks. We have also given pointers to other approaches not discussed in depth in this review. The methods reviewed fall roughly into four model groups: c-means, maximum likelihood, neural networks, and k-nearest neighbor rules. Both supervised and unsupervised schemes require human intervention to obtain clinically useful results in MR segmentation. Unsupervised techniques require somewhat less interaction on a per patient/image basis. Maximum likelihood techniques have had some success, but are very susceptible to the choice of training region, which may need to be chosen slice by slice for even one patient. Generally, techniques that must assume an underlying statistical distribution of the data (such as LML and UML) do not appear promising, since tissue regions of interest do not usually obey the distributional tendencies of probability density functions. The most promising supervised techniques reviewed seem to be FF/NN methods that allow hidden layers to be configured as examples are presented to the system. An example of a self-configuring network, FF/CC, was also discussed. The relatively simple k-nearest neighbor rule algorithms (hard and fuzzy) have also shown promise in the supervised category. Unsupervised techniques based upon fuzzy c-means clustering algorithms have also shown great promise in MR image segmentation. Several unsupervised connectionist techniques have recently been experimented with on MR images of the brain and have provided promising initial results. A pixel-intensity-based edge detection algorithm has recently been used to provide promising segmentations of the brain. This is also an unsupervised technique, older versions of which have been susceptible to oversegmenting the image because of the lack of clear boundaries between tissue types or finding uninteresting boundaries between slightly different types of the same tissue. To conclude, we offer some remarks about improving MR segmentation techniques. The better unsupervised techniques are too slow. Improving speed via parallelization and optimization will improve their competitiveness with, e.g., the k-nn rule, which is the fastest technique covered in this review. Another area for development is dynamic cluster validity. Unsupervised methods need better ways to specify and adjust c, the number of tissue classes found by the algorithm. Initialization is a third important area of research. Many of the schemes listed in Table II are sensitive to good initialization, both in terms of the parameters of the design, as well as operator selection of training data.(ABSTRACT TRUNCATED AT 400 WORDS)

Knowledge-based classification and tissue labeling of MR images of human brain.

Presents a knowledge-based approach to automatic classification and tissue labeling of 2D magnetic resonance (MR) images of the human brain. The system consists of 2 components: an unsupervised clustering algorithm and an expert system. MR brain data is initially segmented by the unsupervised algorithm, then the expert system locates a landmark tissue or cluster and analyzes it by matching it with a model or searching in it for an expected feature. The landmark tissue location and its analysis are repeated until a tumor is found or all tissues are labeled. The knowledge base contains information on cluster distribution in feature space and tissue models. Since tissue shapes are irregular, their models and matching are specially designed: 1) qualitative tissue models are defined for brain tissues such as white matter; 2) default reasoning is used to match a model with an MR image; that is, if there is no mismatch between a model and an image, they are taken as matched. The system has been tested with 53 slices of MR images acquired at different times by 2 different scanners. It accurately identifies abnormal slices and provides a partial labeling of the tissues. It provides an accurate complete labeling of all normal tissues in the absence of large amounts of data nonuniformity, as verified by radiologists. Thus the system can be used to provide automatic screening of slices for abnormality. It also provides a first step toward the complete description of abnormal images for use in automatic tumor volume determination.

A comparison of neural network and fuzzy clustering techniques in segmenting magnetic resonance images of the brain.

Magnetic resonance (MR) brain section images are segmented and then synthetically colored to give visual representations of the original data with three approaches: the literal and approximate fuzzy c-means unsupervised clustering algorithms, and a supervised computational neural network. Initial clinical results are presented on normal volunteers and selected patients with brain tumors surrounded by edema. Supervised and unsupervised segmentation techniques provide broadly similar results. Unsupervised fuzzy algorithms were visually observed to show better segmentation when compared with raw image data for volunteer studies. For a more complex segmentation problem with tumor/edema or cerebrospinal fluid boundary, where the tissues have similar MR relaxation behavior, inconsistency in rating among experts was observed, with fuzz-c-means approaches being slightly preferred over feedforward cascade correlation results. Various facets of both approaches, such as supervised versus unsupervised learning, time complexity, and utility for the diagnostic process, are compared.