Jerry L. Prince

Snakes, shapes, and gradient vector flow

Snakes, or active contours, are used extensively in computer vision and image processing applications, particularly to locate object boundaries. Problems associated with initialization and poor convergence to boundary concavities, however, have limited their utility. This paper presents a new external force for active contours, largely solving both problems. This external force, which we call gradient vector flow (GVF), is computed as a diffusion of the gradient vectors of a gray-level or binary edge map derived from the image. It differs fundamentally from traditional snake external forces in that it cannot be written as the negative gradient of a potential function, and the corresponding snake is formulated directly from a force balance condition rather than a variational formulation. Using several two-dimensional (2-D) examples and one three-dimensional (3-D) example, we show that GVF has a large capture range and is able to move snakes into boundary concavities.

Adaptive fuzzy segmentation of magnetic resonance images

An algorithm is presented for the fuzzy segmentation of two-dimensional (2-D) and three-dimensional (3-D) multispectral magnetic resonance (MR) images that have been corrupted by intensity inhomogeneities, also known as shading artifacts. The algorithm is an extension of the 2-D adaptive fuzzy C-means algorithm (2-D AFCM) presented in previous work by the authors. This algorithm models the intensity inhomogeneities as a gain field that causes image intensities to smoothly and slowly vary through the image space. It iteratively adapts to the intensity inhomogeneities and is completely automated. In this paper, the authors fully generalize 2-D AFCM to three-dimensional (3-D) multispectral images. Because of the potential size of 3-D image data, they also describe a new faster multigrid-based algorithm for its implementation. They show, using simulated MR data, that 3-D AFCM yields lower error rates than both the standard fuzzy C-means (FCM) algorithm and two other competing methods, when segmenting corrupted images. Its efficacy is further demonstrated using real 3-D scalar and multispectral MR brain images.

Gradient vector flow: a new external force for snakes

Snakes, or active contours, are used extensively in computer vision and image processing applications, particularly to locate object boundaries. Problems associated with initialization and poor convergence to concave boundaries, however, have limited their utility. This paper develops a new external force for active contours, largely solving both problems. This external force, which we call gradient vector flow (GVF) is computed as a diffusion of the gradient vectors of a gray-level or binary edge map derived from the image. The resultant field has a large capture range and forces active contours into concave regions. Examples on simulated images and one real image are presented.

Cardiac Motion Tracking Using CINE Harmonic Phase (HARP) Magnetic Resonance Imaging

This article introduces a new image processing technique for rapid analysis of tagged cardiac magnetic resonance image sequences. The method uses isolated spectral peaks in SPAMM‐tagged magnetic resonance images, which contain information about cardiac motion. The inverse Fourier transform of a spectral peak is a complex image whose calculated angle is called a harmonic phase (HARP) image. It is shown how two HARP image sequences can be used to automatically and accurately track material points through time. A rapid, semiautomated procedure to calculate circumferential and radial Lagrangian strain from tracked points is described. This new computational approach permits rapid analysis and visualization of myocardial strain within 5–10 min after the scan is complete. Its performance is demonstrated on MR image sequences reflecting both normal and abnormal cardiac motion. Results from the new method are shown to compare very well with a previously validated tracking algorithm. Magn Reson Med 42:1048–1060, 1999.

Generalized gradient vector flow external forces for active contours

Active contours, or snakes, are used extensively in computer vision and image processing applications, particularly to locate object boundaries. A new type of external force for active contours, called gradient vector flow (GVF) was introduced recently to address problems associated with initialization and poor convergence to boundary concavities. GVF is computed as a di⁄usion of the gradient vectors of a gray-level or binary edge map derived from the image. In this paper, we generalize the GVF formulation to include two spatially varying weighting functions. This improves active contour convergence to long, thin boundary indentations, while maintaining other desirable properties of GVF, such as an extended capture range. The original GVF is a special case of this new generalized GVF (GGVF) model. An error analysis for active contour results on simulated test images is also presented. ( 1998 Elsevier Science B.V. All rights reserved.

Measurement of radiotracer concentration in brain gray matter using positron emission tomography: MRI-based correction for partial volume effects.

Accuracy in in vivo quantitation of brain function with positron emission tomography (PET) has often been limited by partial volume effects. This limitation becomes prominent in studies of aging and degenerative brain diseases where partial volume effects vary with different degrees of atrophy. The present study describes how the actual gray matter (GM) tracer concentration can be estimated using an algorithm that relates the regional fraction of GM to partial volume effects. The regional fraction of GM was determined by magnetic resonance imaging (MRI). The procedure is designated as GM PET. In computer simulations and phantom studies, the GM PET algorithm permitted a 100% recovery of the actual tracer concentration in neocortical GM and hippocampus, irrespective of the GM volume. GM PET was applied in a test case of temporal lobe epilepsy revealing an increase in radiotracer activity in GM that was undetected in the PET image before correction for partial volume effects. In computer simulations, errors in the segmentation of GM and errors in registration of PET and MRI images resulted in less than 15% inaccuracy in the GM PET image. In conclusion, GM PET permits accurate determination of the actual radiotracer concentration in human brain GM in vivo. The method differentiates whether a change in the apparent radiotracer concentration reflects solely an alteration in GM volume or rather a change in radiotracer concentration per unit volume of GM.

A topology preserving level set method for geometric deformable models

Active contour and surface models, also known as deformable models, are powerful image segmentation techniques. Geometric deformable models implemented using level set methods have advantages over parametric models due to their intrinsic behavior, parameterization independence, and ease of implementation. However, a long claimed advantage of geometric deformable models-the ability to automatically handle topology changes-turns out to be a liability in applications where the object to be segmented has a known topology that must be preserved. We present a new class of geometric deformable models designed using a novel topology-preserving level set method, which achieves topology preservation by applying the simple point concept from digital topology. These new models maintain the other advantages of standard geometric deformable models including subpixel accuracy and production of nonintersecting curves or surfaces. Moreover, since the topology-preserving constraint is enforced efficiently through local computations, the resulting algorithm incurs only nominal computational overhead over standard geometric deformable models. Several experiments on simulated and real data are provided to demonstrate the performance of this new deformable model algorithm.

Imaging heart motion using harmonic phase MRI

Describes a new image processing technique for rapid analysis and visualization of tagged cardiac magnetic resonance (MR) images. The method is based on the use of isolated spectral peaks in spatial modulation of magnetization (SPAMM)-tagged magnetic resonance images. The authors call the calculated angle of the complex image corresponding to one of these peaks a harmonic phase (HARP) image and show that HARP images can be used to synthesize conventional tag lines, reconstruct displacement fields for small motions, and calculate two-dimensional (2-D) strain. The performance of this new approach is demonstrated using both real and simulated tagged MR images. Potential for use of HARP images in fast imaging techniques and three-dimensional (3-D) analyses are discussed.

A computerized approach for morphological analysis of the corpus callosum

OBJECTIVE A new technique for analyzing the morphology of the corpus callosum is presented, and it is applied to a group of elderly subjects. MATERIALS AND METHODS The proposed approach normalizes subject data into the Talairach space using an elastic deformation transformation. The properties of this transformation are used as a quantitative description of the callosal shape with respect to the Talairach atlas, which is treated as a standard. In particular, a deformation function measures the enlargement/shrinkage associated with this elastic deformation. Intersubject comparisons are made by comparing deformation functions. RESULTS This technique was applied to eight male and eight female subjects. Based on the average deformation functions of each group, the posterior region of the female corpus callosum was found to be larger than its corresponding region in the males. The average callosal shape of each group was also found, demonstrating visually the callosal shape differences between the two groups in this sample. CONCLUSION The proposed methodology utilizes the full resolution of the data, rather than relying on global descriptions such as area measurements. The application of this methodology to an elderly group indicated sex-related differences in the callosal shape and size.

Tag and contour detection in tagged MR images of the left ventricle

Tracking magnetic resonance tags in myocardial tissue promises to be an effective tool for the assessment of myocardial motion. The authors describe a hierarchy of image processing steps which rapidly detects both the contours of the myocardial boundaries of the left ventricle and the tags within the myocardium. The method works on both short axis and long axis images containing radial and parallel tag patterns, respectively. Left ventricular boundaries are detected by first removing the tags using morphological closing and then selecting candidate edge points. The best inner and outer boundaries are found using a dynamic program that minimizes a nonlinear combination of several local cost functions. Tags are tracked by matching a template of their expected profile using a least squares estimate. Since blood pooling, contiguous and adjacent tissue, and motion artifacts sometimes cause detection errors, a graphical user interface was developed to allow user correction of anomalous points. The authors present results on several tagged images of a human. A fully automated run generally finds the endocardial boundary and the tag lines extremely well, requiring very little manual correction. The epicardial boundary sometimes requires more intervention to obtain an acceptable result. These methods are currently being used in the analysis of cardiac strain and as a basis for the analysis of alternate tag geometries.