Mark P. Wilson

Digital stereo image analyzer for generating automated 3-D measures of optic disc deformation in glaucoma

The major limitations of precise evaluation of retinal structures in present clinical situations are the lack of standardization, the inherent subjectivity involved in the interpretation of retinal images, and intra- as well as interobserver variability. While evaluating optic disc deformation in glaucoma, these limitations could be overcome by using advanced digital image analysis techniques to generate precise metrics; from stereo optic disc image pairs. A digital stereovision system for visualizing the topography of the optic nerve head from stereo optic disc images is presented. We have developed an algorithm, combining power cepstrum and zero-mean-normalized cross correlation techniques, which extracts depth information using coarse-to-fine disparity between corresponding windows in a stereo pair. The gray level encoded sparse disparity matrix is subjected to a cubic B-spline operation to generate smooth representations of the optic cup/disc surfaces and new three-dimensional (3-D) metrics from isodisparity contours. Despite the challenges involved in 3-D surface recovery, the robustness of our algorithm in finding disparities within the constraints used has been validated using stereo pairs with known disparities. In a preliminary longitudinal study of glaucoma patients, a strong correlation is found between the computer-generated quantitative cup/disc volume metrics and manual metrics commonly used in a clinic. The computer generated new metrics, however, eliminate the subjective variability and greatly reduce the time and cost involved in manual metric generation in follow-up studies of glaucoma.

Computer-aided methods for quantitative assessment of longitudinal changes in retinal images presenting with maculopathy

This paper presents the results from applying a computer- based methodology for making precise measurements of longitudinal changes in a patients digital retinal images presenting with age-related macular degeneration. The digital retinal image analysis system applies recognized principles in automatic image segmentation and integrates the automation with a graphical user interface. Drusen, retinal lesions associated with age-related macular degeneration (ARMD), were segmented using a region-growing algorithm. The algorithm calculates the 76 percentile intensity in a region to provide seed points for the neighborhood-growing algorithm. Twenty-one cases were analyzed. Agreement statistics (kappa) were determined by comparing the automated results with those provided from manually derived measurements. Agreement statistics ranged from 0.49 to 0.71 for different regions of the retina. The manual analysis ground truth was performed by trained graders from the University of Wisconsin Reading Center using guidelines found in the Wisconsin Age-Related Maculopathy Degeneration Grading Scheme (WARMGS). Because of the time required, the ophthalmic graders can only grade (size, area, type) the most prominent drusen in specific regions, resulting in a small sampling of drusen lesions in the retina. The computer-based approach allows one to efficiently and comprehensively grade all of the lesions for larger numbers of images. The additional advantage, however, is in the precision and total area that can be graded with the computer-aided technology. Computer-registered longitudinal images produced a precise determination of the temporal changes in the individual lesions. This study has demonstrated a robust segmentation and registration methodology for automatic and semiautomatic detection and measurement of abnormal regions in longitudinal retinal images.

Automated detection of microcalcifications in mammograms through application of image pixel remapping and statistical filter

Mammography and finding of suspicious masses during self-examinations and clinical breast examinations form the primary screening tools for early detection of breast cancer. Mammography is essential for early detection of cancer, prior to any other means. Many abnormalities other than cancer are shown by screening mammography, leading to biopsy of the suspicious defect. Any enhancement that reduces the number of unnecessary biopsies is welcomed. Digital enhancement may aid in early detection of some patterns, such as microcalcification clusters, indicating the onset of DCIS (ductal carcinoma in situ) that accounts for 20% of all mammographically detected breast cancers, and early detection may enable a complete cure. The individual calcifications are hard to detect, due to their size and shape variability and the inhomogeneous background structure. Our study addresses only early detection of microcalcifications. We present an algorithm which locates microcalcifications based on local gray-scale variability, tissue structures and image statistics. The mammographs are digitally enhanced to accent textures and a previously developed threshold filter creates a binary mask of the calcifications spatial location.

Full automation of morphological segmentation of retinal images: a comparison with human-based analysis

Age-Related Macular Degeneration (ARMD) is the leading cause of irreversible visual loss among the elderly in the US and Europe. A computer-based system has been developed to provide the ability to track the position and margin of the ARMD associated lesion; drusen. Variations in the subjects retinal pigmentation, size and profusion of the lesions, and differences in image illumination and quality present significant challenges to most segmentation algorithms. An algorithm is presented that first classifies the image to optimize the variables of a mathematical morphology algorithm. A binary image is found by applying Otsus method to the reconstructed image. Lesion size and area distribution statistics are then calculated. For training and validation, the University of Wisconsin provided longitudinal images of 22 subjects from their 10 year Beaver Dam Study. Using the Wisconsin Age-Related Maculopathy Grading System, three graders classified the retinal images according to drusen size and area of involvement. The percentages within the acceptable error between the three graders and the computer are as follows: Grader-A: Area: 84% Size: 81%; Grader-B: Area: 63% Size: 76%; Grader-C: Area: 81% Size: 88%. To validate the segmented position and boundary one grader was asked to digitally outline the drusen boundary. The average accuracy based on sensitivity and specificity was 0.87 for thirty four marked regions.

A JAVA-based system for segmentation and analysis of retinal images

This paper presents a platform independent system, RetinaView, that is designed for manual analysis as well as automatic segmentation and characterization of multi-modality retinal images and videos. A user-friendly interface that allows display of individual retinal images in four or more windows has been designed for simultaneous viewing of multiple modalities in research-based retinal studies. A variety of tools have been developed for RetinaView, including automatic segmentation of retinal vasculature, segmentation of small vascular abnormalities found in early stages of diabetic retinopathy (microaneurysms), segmentation of drusen, lesions pathoneumonic of age related macular degeneration (ARMD), analysis of fluorescein and indocyanine green videos, and manual annotation capabilities and fovea location. This system will enable ophthalmological researchers, clinicians, ophthalmologists to study, diagnose, and monitor the progression of pathologies in retinal images with much greater precision than is possible with totally manual techniques. This paper presents an overview of the RetinaView system and gives examples of its application in segmentation of normal anatomical features as well as pathological lesions.

Digital stereo-optic disc image analyzer for monitoring progression of glaucoma

This paper describes an automated 3-D surface recovery algorithm for consistent and quantitative evaluation of the deformation in the ONH (optic nerve head). Additional measures, such as the changes in the volume of the cup and the disc as an improvement to the traditional cup to disc ratios, can thus be developed for longitudinal follow-up study of a patient. We propose an automated computerized technique for stereo pair registration and surface visualization of the ONH. Power cepstrum and zero mean cross correlation are embedded in the registration and a 3-D surface recovery technique is proposed. Preprocessing, as well as an overall registration, is performed upon stereo pairs. Then a coarse to fine feature matching strategy is used to reduce the ambiguity in finding the conjugate pair of the same point within the constraints of the epipolar plane. A cubic B-spline interpolation smooths the representation of the ONH obtained, while superimposition of features such as blood vessels is added. Studies show high correlation between traditional cup/disc measures derived from manual segmentation by ophthalmologists and computer generated cup/disc volume ratio. Such longitudinal studies over a large population of glaucoma patients are currently in progress for validation of the surface recovery algorithm.

Optimal scanning, display, and segmentation of the International Labor Organization (ILO) X-ray images set for pneumoconiosis

A method for scanning and displaying chest radiographs (X-rays) is presented. The new method treats the scanning as independent of the image display, allowing for maximum information content to be captured during the scanning process, and then for this information to be optimally displayed using a new function for maximizing the contrast variation throughout the image. The quality of the digitized X-ray images was compared against the original X-ray films and was found to be of comparable visualization quality. The rib parenchyma is then segmented by an active shape model.

Optimizing retinal image digitization for improved digital processing and visualization

Determines if a relationship exists between quantitative measures of digital retinal image quality, such as contrast and spatial resolution, and perceptual preferences of human graders. Ad-hoc digitization of retinal 35-mm color slides is less than optimal due to the unusual light scattering properties of the retina. Quantitative and qualitative measures of contrast and entropy were calculated for two test images. The resolution was determined by a scanning procedure. A pairwise comparison technique was adopted for measuring the medical technicians visual preference, and the results were analyzed using a formal quantitative methodology. The results showed that measures such as entropy and contrast generally correlated with the perceived image quality. A single, universally applicable set of scanning parameters was not found that would consistently produce the most preferred image for all image quality criteria.

Performance of multiresolution pattern classifiers in medical image encoding from wavelet coefficient distributions

The fidelity of the reconstructed image in an image coding/decoding scheme and the lowest transmission bit rate from rate-distortion theory can be predicted provided the image statistics are known. Currently popular subband image coding assumes Gaussian source with memory for optimal performance. However, most images do not follow the ideal distribution. The advantage of subband coding lies in the fact that the wavelet coefficients in decomposed subimages have probability distribution functions (pdfs) that can be modeled as a generalized Gaussian when proper parameters are chosen experimentally. However, the filter length chosen for digital implementation of a specific wavelet is crucial in shaping the pdf characteristics and hence in the ability to predict the achievable bit rate at minimum distortion in a quantization scheme. We have analyzed the pdfs of a number of wavelets and chosen filter lengths providing the best fit to a generalized Gaussian distribution for encoding an image by vector quantization of multiresolution wavelet subimages using an adaptive clustering. Our results demonstrate that the performance of the adaptive vector quantizer improves significantly when wavelet filter lengths are chosen to fit the generalized Gaussian distribution.

Feature extraction and segmentation in medical images by statistical optimization and point operation approaches

Feature extraction is a critical preprocessing step, which influences the outcome of the entire process of developing significant metrics for medical image evaluation. The purpose of this paper is firstly to compare the effect of an optimized statistical feature extraction methodology to a well designed combination of point operations for feature extraction at the preprocessing stage of retinal images for developing useful diagnostic metrics for retinal diseases such as glaucoma and diabetic retinopathy. Segmentation of the extracted features allow us to investigate the effect of occlusion induced by these features on generating stereo disparity mapping and 3-D visualization of the optic cup/disc. Segmentation of blood vessels in the retina also has significant application in generating precise vessel diameter metrics in vascular diseases such as hypertension and diabetic retinopathy for monitoring progression of retinal diseases.