Martin D. Fox

Level set evolution without re-initialization: a new variational formulation

In this paper, we present a new variational formulation for geometric active contours that forces the level set function to be close to a signed distance function, and therefore completely eliminates the need of the costly re-initialization procedure. Our variational formulation consists of an internal energy term that penalizes the deviation of the level set function from a signed distance function, and an external energy term that drives the motion of the zero level set toward the desired image features, such as object boundaries. The resulting evolution of the level set function is the gradient flow that minimizes the overall energy functional. The proposed variational level set formulation has three main advantages over the traditional level set formulations. First, a significantly larger time step can be used for numerically solving the evolution partial differential equation, and therefore speeds up the curve evolution. Second, the level set function can be initialized with general functions that are more efficient to construct and easier to use in practice than the widely used signed distance function. Third, the level set evolution in our formulation can be easily implemented by simple finite difference scheme and is computationally more efficient. The proposed algorithm has been applied to both simulated and real images with promising results.

Distance Regularized Level Set Evolution and Its Application to Image Segmentation

Level set methods have been widely used in image processing and computer vision. In conventional level set formulations, the level set function typically develops irregularities during its evolution, which may cause numerical errors and eventually destroy the stability of the evolution. Therefore, a numerical remedy, called reinitialization, is typically applied to periodically replace the degraded level set function with a signed distance function. However, the practice of reinitialization not only raises serious problems as when and how it should be performed, but also affects numerical accuracy in an undesirable way. This paper proposes a new variational level set formulation in which the regularity of the level set function is intrinsically maintained during the level set evolution. The level set evolution is derived as the gradient flow that minimizes an energy functional with a distance regularization term and an external energy that drives the motion of the zero level set toward desired locations. The distance regularization term is defined with a potential function such that the derived level set evolution has a unique forward-and-backward (FAB) diffusion effect, which is able to maintain a desired shape of the level set function, particularly a signed distance profile near the zero level set. This yields a new type of level set evolution called distance regularized level set evolution (DRLSE). The distance regularization effect eliminates the need for reinitialization and thereby avoids its induced numerical errors. In contrast to complicated implementations of conventional level set formulations, a simpler and more efficient finite difference scheme can be used to implement the DRLSE formulation. DRLSE also allows the use of more general and efficient initialization of the level set function. In its numerical implementation, relatively large time steps can be used in the finite difference scheme to reduce the number of iterations, while ensuring sufficient numerical accuracy. To demonstrate the effectiveness of the DRLSE formulation, we apply it to an edge-based active contour model for image segmentation, and provide a simple narrowband implementation to greatly reduce computational cost.

Noninvasive optical polarimetric glucose sensing using a true phase measurement technique

The development and testing of a noninvasive true phase optical polarimetry sensing system to monitor in vivo glucose concentrations is described. To demonstrate the applicability of this optical sensor for glucose movement, the authors calibrate the system and then test it in vitro using both a glass test cell filled with glucose solution in the physiologic range, with a path length of 0.9 cm to approximate the 1-cm path length present in the anterior chamber of the eye, and then on an excised human eye. The technique used helium neon laser light which was coupled through a rotating linear polarizer along with two stationary linear polarizers and two detectors to produce reference and signal outputs whose amplitudes varied sinusoidally with a frequency of twice the angular velocity of the rotating polarizer, and whose phase was proportional to the rotation of the linear polarization vector passing through the glucose solution.<<ETX>>

Classifying mammographic lesions using computerized image analysis

The classification of 3 common breast lesions, fibroadenomas, cysts, and cancers, was achieved using computerized image analysis of tumor shape in conjunction with patient age. The process involved the digitization of 69 mammographic images using a video camera and a commercial frame grabber on a PC-based computer system. An interactive segmentation procedure identified the tumor boundary using a thresholding technique which successfully segmented 57% of the lesions. Several features were chosen based on the gross and fine shape describing properties of the tumor boundaries as seen on the radiographs. Patient age was included as a significant feature in determining whether the tumor was a cyst, fibroadenoma, or cancer and was the only patient history information available for this study. The concept of a radial length measure provided a basis from which 6 of the 7 shape describing features were chosen, the seventh being tumor circularity. The feature selection process was accomplished using linear discriminant analysis and a Euclidean distance metric determined group membership. The effectiveness of the classification scheme was tested using both the apparent and the leaving-one-out test methods. The best results using the apparent test method resulted in correctly classifying 82% of the tumors segmented using the entire feature space and the highest classification rate using the leaving-one-out test method was 69% using a subset of the feature space. The results using only the shape descriptors, and excluding patient age resulted in correctly classifying 72% using the entire feature space (except age), and 51% using a subset of the feature space.

Segmentation of external force field for automatic initialization and splitting of snakes

Active contours or snakes have been extensively utilized in handling image segmentation and classification problems. In traditional active contour models, snake initialization is performed manually by users, and topological changes, such as splitting of the snake, cannot be automatically handled. In this paper, we introduce a new method to solve the snake initialization and splitting problem, based on an area segmentation approach: the external force field is segmented first, and then the snake initialization and splitting can be automatically performed by using the segmented external force field. Such initialization and splitting produces multiple snakes, each of which is within the capture range associated to an object and will be evolved to the object boundary. The external force used in this paper is a gradient vector flow with an edge-preserving property (EPGVF), which can prevent the snakes from passing over weak boundaries. To segment the external force field, we represent it with a graph, and a graph-theory approach can be taken to determine the membership of each pixel. Experimental results establish the effectiveness of the proposed approach.

Segmentation of edge preserving gradient vector flow: an approach toward automatically initializing and splitting of snakes

Active contours or snakes have been extensively utilized in handling image segmentation and classification problems. In traditional active contour models, snake initialization is performed manually by users, and topological changes, such as splitting of the snake, can not be automatically handled. In this paper, we introduce a new method to solve the snake initialization and splitting problem, based on an area segmentation approach: the external force field is segmented first, and then the snake initialization and splitting can be automatically performed by using the segmented external force field. Such initialization and splitting produces multiple snakes, each of which is within the capture range associated to an object and evolved to the object boundary. The external force used in this paper is a gradient vector flow with an edge-preserving property (EPGVF), which can prevent the snakes from passing over weak boundaries. To segment the external force field, we represent it with a graph, and a graph-theory approach can be taken to determine the membership of each pixel. Experimental results establish the effectiveness of the proposed approach.

Fast Hartley transforms for image processing

The fast Hartley transform (FHT) is used to transform two-dimensional image data. Because the Hartley transform is real-valued, it does not require complex operations. Both spectra and autocorrelations of two-dimensional ultrasound images of normal and abnormal livers were computed.

Ultrasound image enhancement: A review

Abstract Medical ultrasound imaging uses pulsed acoustic waves that are transmitted and received by a hand-held transducer. This is a mature technology that it is widely used around the world. Among its advantages are that it is cost-effective, flexible, and does not require ionizing radiation. However, the image quality is affected by degradation of ultrasound signals when propagating through biological tissues. Many efforts have been done in the last three decades to improve the quality of the images. This paper reviews some of the most important methods for ultrasound enhancement. We classified these techniques into two groups: preprocessing and post-processing, analyzed their benefits and limitations, and presented our beliefs about where ultrasound research could be directed to, in order to improve its effectiveness and broaden its applications.

Three-dimensional Doppler velocimetry of flow jets

A closed-form solution is derived of the Doppler equation for the magnitude and angle to the three-dimensional velocity vector. The resultant solution is more general than previous formulations and, since it is based on multiple conventional transmit-receive Doppler probes, it can be readily applied by adapting existing Doppler units. Turntable experiments demonstrated that the three-dimensional angle-independent theory correctly predicted velocity within average 32% over a 26. degrees -range of Doppler angles. Experiments with a flow phantom showed general agreement with the angle-independent theory in the more demanding setting of an actual flow stream. Experiments utilizing a jet stream which modeled the clinically important conditions that exist in stenotic vessels, valvular abnormalities, and septal defects within 1.4-3.8% of the actual velocity, over an 18 degrees range of Doppler angles.<<ETX>>

Enhancement of chest radiographs with gradient operators

Reference is made to the Sobel and Roberts gradient operators used to enhance image edges. Overall, the Sobel operator was found to be superior to the Roberts operator in edge enhancement. A theoretical explanation for the superior performance of the Sobel operator was developed based on the concept of analyzing the x and y Sobel masks as linear filters. By applying pill-box, Gaussian, or median filtering prior to applying a gradient operator, noise was reduced. The pill-box and Gaussian filters were more computationally efficient than the median filter with approximately equal effectiveness in noise reduction.