Normalization of series of fundus images to monitor the geographic atrophy growth in dry age-related macular degeneration.

BACKGROUND AND OBJECTIVE

Age-related macular degeneration (ARMD) is a degenerative disease that affects the retina, and the leading cause of visual loss. In its dry form, the pathology is characterized by the progressive, centrifugal expansion of retinal lesions, called geographic atrophy (GA). In infrared eye fundus images, the GA appears as localized bright areas and its growth can be observed in series of images acquired at regular time intervals. However, illumination distortions between the images make impossible the direct comparison of intensities in order to study the GA progress. Here, we propose a new method to compensate for illumination distortion between images.

METHODS

We process all images of the series so that any two images have comparable gray levels. Our approach relies on an illumination/reflectance model. We first estimate the pixel-wise illumination ratio between any two images of the series, in a recursive way; then we correct each image against all the others, based on those estimates. The algorithm is applied on a sliding temporal window to cope with large changes in reflectance. We also propose morphological processing to suppress illumination artefacts.

RESULTS

The corrected illumination function is homogeneous in the series, enabling the direct comparison of grey-levels intensities in each pixel, and so the detection of the GA growth between any two images. To demonstrate that, we present numerous experiments performed on a dataset of 18 series (328 images), manually segmented by an ophthalmologist. First, we show that the normalization preprocessing dramatically increases the contrast of the GA growth areas. Secondly, we apply segmentation algorithms derived from Otsu's thresholding to detect automatically the GA total growth and the GA progress between consecutive images. We demonstrate qualitatively and quantitatively that these algorithms, although fully automatic, unsupervised and basic, already lead to interesting segmentation results when applied to the normalized images. Colored maps representing the GA evolution can be derived from the segmentations.

CONCLUSION

To our knowledge, the proposed method is the first one which corrects automatically and jointly the illumination inhomogeneity in a series of fundus images, regardless of the number of images, the size, shape and progression of lesion areas. This algorithm greatly facilitates the visual interpretation by the medical expert. It opens up the possibility of treating automatically each series as a whole (not just in pairs of images) to model the GA growth.