Automated analysis of digital fundus autofluorescence images of geographic atrophy in advanced age-related macular degeneration using confocal scanning laser ophthalmoscopy (cSLO)
© Deckert et al; licensee BioMed Central Ltd. 2005
Received: 28 October 2004
Accepted: 06 April 2005
Published: 06 April 2005
Fundus autofluorescence (AF) imaging using confocal scanning laser ophthalmoscopy (cSLO) provides an accurate delineation of areas of geographic atrophy (GA).
Automated computer-assisted methods for detecting and removing interfering vessels are needed to support the GA quantification process in longitudinal studies and in reading centres.
A test tool was implemented that uses region-growing techniques to segment GA areas. An algorithm for illuminating shadows can be used to process low-quality images. Agreement between observers and between three different methods was evaluated by two independent readers in a pilot study. Agreement and objectivity were assessed using the Bland-Altman approach.
The new method (C) identifies vascular structures that interfere with the delineation of GA. Results are comparable to those of two commonly used procedures (A, B), with a mean difference between C and A of -0.67 mm2 (95% CI [-0.99, -0.36]), between B and A of -0.81 mm2, (95% CI [-1.08, -0.53]), and between C and B of 0.15 mm2 (95% CI [-0.12, 0.41]). Objectivity of a method is quantified by the mean difference between observers: A 0.30 mm2 (95% CI [0.02, 0.57]), B -0.11 mm2 (95% CI [-0.28, 0.10]), and C 0.12 mm2 (95% CI [0.02, 0.22]).
The novel procedure is comparable with regard to objectivity and inter-reader agreement to established methods of quantifying GA. It considerably speeds up the lengthy measurement process in AF with well defined GA zones.
Age-related macular degeneration (AMD) is the most common cause of legal blindness among industrialized nations in the population aged 50 years and above [1–4]. Besides choroidal neovascularization and detachments of the retinal pigment epithelium, geographic atrophy (GA) of retinal pigment epithelium (RPE) is a common cause of severe visual loss in patients with AMD [5–7]. Changes in time can be documented by fundus autofluorescence images (AF) mediated by RPE lipofuscin accumulations and its spatial distribution over retinal areas, obtained in vivo using a confocal scanning laser ophthalmosocope (Heidelberg Retina Angiograph (HRA), Dossenheim, Germany) [8–10]. Areas of GA are usually associated with a well-defined zone of decreased autofluorescence due to the absence of fluorophores residing in RPE lipofuscin granules [11, 12]. The deduction of clinically relevant information from these pictures is a complex process which should be optimised and automated, especially in the context of multicenter studies.
Generally, AF are recorded using a confocal scanning laser ophthalmoscope (Heidelberg Retina Angiograph, HRA, Heidelberg Engineering, Germany; which includes the Heidelberg Eye Explorer (HEE) software package). The images are immediately digitised and processed using a flexible frame processor and subsequently displayed on a computer screen. Corresponding with funduscopically visible atrophic areas, fundus intensity of AF is markedly decreased .
With method A, atrophic areas are outlined on the screen using the mouse-driven cursor of the HEE software program. The areas are then measured and the data exported manually to an Excel spreadsheet by cut and paste. This completely manual, mouse-driven method A is time-consuming and can exaggerate subjective impressions. Mistakes can occur as a result of the error-prone interface between the user's hand and the computer mouse . This implies that the accuracy of mouse-driven contour painting depends not only on subjective impressions, but also on the user's dexterity.
With method B, the images are exported as bitmap files from the HEE program. Interfering vascular structures, which appear as dark atrophic areas, are manually repainted white using the mouse-driven paintbrush of Microsoft Paint. The modified images are then transferred to Global Lab Image 2 and the remaining dark areas are measured using a threshold procedure tool. The resulting data are exported to Microsoft Excel manually by cut and paste . This semi-automated method B only allows an interpretation using high-quality AF images. It requires the manual, mouse-driven whitewashing of interfering vessels that touch upon, or extend into, the GA area. However, as the vessels exhibit grey levels that are similar to atrophic patches, automated segmentation is not possible . Both methods rely on the circumstantial handling of different software tools that have not been specifically adapted to the problems of GA measurements.
We have developed a novel customized image analysis test tool that includes an adapted algorithm for automated identification of interfering vascular structures, and compared this new method with the previous ones.
The test tool combines all steps of GA area measurement, including the automatic export of data into an Excel spreadsheet for further analysis. The method is based on region-growing, instead of segmentation with a threshold value similar to that of method B, allowing the reader to sort out non-atrophic areas with similarly diminished grey values as in actual GA areas.
After selecting the GA area by moving the mouse cursor and clicking on the region of interest, a first segmentation is started using a default parameter value calculated by the image's mean grey value. Automatic correction of the generated contour is possible by adjusting the associated parameter value with a relocatable graphical element. Holes within the detected GA area can be identified and further GA areas in the same image can be integrated.
For segmentation of GA areas with interfering vascular structures, the tool includes an optional algorithm for detecting and whitewashing such structures . The number of detected vessels depends on the individual settings of the parameters for vessels diameter, length and cross-linking. Incorrect contour segments of the GA caused by whitewashed vessel stubs can reach into the GA area and are corrected, either automatically with a default parameter or by user interaction.
The algorithm eliminates vessel stubs within the GA area and corrects the contour segments in these regions. One disadvantage of this process is that small contour jags which were correctly detected a step before now disappear, and instead the contour is minimally smoothed and widened. As a result the process produces correct contour segments in passages from GA area to vessels and incorrectly widened contour segments otherwise. A new segmentation step produces the correct contour without interfering vessels.
Using well-aligned relocatable convex and radial hulls facilitates fine-tuning of the actual delineation and measurement of GA areas. These tools enable alternative contour finding in difficult contour segments with tangentially interfering vessels. Poor quality images often contain large shady outer areas causing erupting segmentation of GA areas before identifying the true borders of the GA. This can be countered by the method's option for illuminating such areas of shade. The original image is shown in a separate window. An electronic magnifier is integrated to support contour finding in difficult sections as well as for small GA areas.
Validation is based on the material collected in a previous study , with published data from the right eye only. Data for both eyes are available  and are used for validation. The same readers evaluated the same material, but using the new method. The images were evaluated in random order and the readers had no access to the previous results based on methods A and B. With the unpublished material from 2001, our validation sample consists of 10 left and 24 right eyes. The Bland-Altman Design for method comparison studies  was applied. The limits of agreement (REF) were calculated for each comparison. Nonparametric tests for paired and unpaired observations as well as ANOVA (Analysis of Variance) were used for statistical assessment. The significance level was set to α = 0.05.
The study was a sub-study of the FAM study which followed the tenets of the Declaration of Helsinki and was approved by the Ethics Committee of the University of Heidelberg. Informed consent was obtained from the patients prior to recruitment into the study.
The complete GA measurement data set
reader 1 (mm2)
reader 2 (mm2)
ID right eyes
ID left eyes
ID in 2001 excluded eyes
Quality measures to assess method comparison
Mean difference (mm2)
95% confidence interval
Limits of agreement
Comparability, repeatability, and objectivity are crucial factors to consider when developing a new method. Comparability with at least two methods has to be determined; objectivity calls for a multitude of readers; and repeatability requires multiple measurements by the same reader using each method. An appropriate study would randomly allocate image and reader in a cross-over design.
Bland-Altman plots visualize comparability, repeatability, and objectivity in terms of the limits of agreement in which 95% of data should appear. Furthermore, the data within these borders must lie within the region of clinically irrelevant differences. Therewith, verification of a clinically comparable quality can be obtained. The measurements for method comparison have to be performed at the same time to exclude subjective effects caused by changing user criteria rather than by the methods themselves.
Additionally, improved effectiveness must be demonstrated to justify replacement of methods . Effectiveness is defined here as the number of successfully assessed AF images and the time required to complete the assessment procedure. In 2001  no exact measurements of reading duration were taken. However, both readers had the impression that method C speeds up the reading process. Furthermore, low-quality images were excluded from the readings in 2001 . Some of these pictures could be handled with method C. Again both readers agreed that method C was better suited to evaluate low-quality images. A further study has to put these subjective impressions into an appropriate objective framework to demonstrate effectiveness.
We conclude that method C is not inferior to the two commonly methods used in measuring GA areas. This has been shown by the comparability of the measured values in the Bland-Altman design. Thereby, the inter-methods comparisons align with the degree of decision freedom for each method. Method A  enables the setting of each contour pixel individually with no relation between contour pixels. Method B  has only one degree of decision freedom using one threshold value for the whole contour, and method C has fewer degrees of decision freedom than A, but more than B: it allows for the individual exclusion of non-atrophic holes or a fine-tuning of critical contour-sections. In accordance with the degrees of decision freedom, the mean values of the new method C lie between those of method A and B (data not shown). Altogether, method C is similar more to method B than to method A (see table 2, A-C, B-C). Increasing the number of degrees of decision freedom increases the influence of subjectivity. It also gives a larger impact of user's competence.
There are subjective inter-reader effects and effects caused by changing assessment criteria within the same reader over a period of years due to increasing experience or changing criteria, but not due to the method used. This influences the bias in method comparison and therewith the statements about measurement error. Considering this confounder, method comparison can only show that method C is no worse than methods A and B. The effect of the influence by the subjectivity of readers is most distinct in method A (see table 2).
Heterogeneity in agreement between the readers over the methods was demonstrated by a significant interaction term (reader by method) in an ANOVA model. There is a loss of histological information due to the process of producing AF images from the retina by confocal laser scanning ophthalmoscopy. Several underlying factors for alterations in grey values in the AF images, mainly in the border region of GA areas, force the reader to interpret by his subjective impressions, more or less supported by the interaction between him and the method.
The study is a pilot study, yet it has provided important information for the further design and testing of method C. Method C has been developed and launched in consideration of several viewpoints. One aim has been to reduce the influence of manual skills and the number of procedural steps, while another important goal was not to restrict the competence of the medical user.
Method C can be implemented to quickly assess both unproblematic AF images and, with the additional accessory tools, difficult AF images. For example, the presence of shadows within the marginal areas of the images appears to be the main reason for poor quality. Presenting a processed image with illuminated shadows together with the original image facilitates segmentation by different algorithms. Furthermore, robust algorithms have been integrated to rectify segmentation and to eliminate interfering vessels in an effective way. Hence, a large proportion of GAs can be measured in a short time using vessel detection and default parameters.
The individual nature of each AF image and the wide variation of possible combinations of features suggest that it would be impossible to develop a clearly defined method that could measure GA areas with only a few procedural steps. Thus, it seems to be recommendable to represent in the future a combination of several methods, with adaptations depending on the quality of AF images, in one tool. However, in attempting to decide upon a particular combination of method parts, some basic questions arose during the process of developing and testing method C, which have still to be clarified.
The selection of further algorithms and methods will depend on the definition of the clinical relevance of differences between GA areas. Consequently, care should be taken that no clinical relevance is attached to artificial differences due solely to repeated use of a method or to different users. Furthermore, a method comparison using the Bland-Altman design will only be meaningful if there are well-defined limits of agreement and clinically relevant bounds . For a method to be useful, the probability of differences between repeated measures and objectivity must not transgress this bound significantly. If two methods produce absolute differences below this limit with non-significant probabilities, then the method with low dispersion but nearer the bound should be preferred. It should be favoured over a method with less bias but larger dispersion. The reason is that a uniform bias could be eliminated by training or considered in the statistical evaluation. With regard to objectivity, it is also important to know whether high objectivity by one threshold value like in method B limits precise contours and facilitates medically correct detection of GA borders. If the precision in imaging the GA borders is important for quantification of disease progression, a segmentation algorithm like region-growing with more degrees of decision freedom should be used. For the decision, clinical relevance in method comparison should be the main criteria, too.
Supported by the Deutsche Forschungsgemeinschaft (DFG) Research Priority Program SPP 1088: MA 1723/1-1 and HO 1926/1-2
- Bressler NM, Bressler SB, Fine SL: Age-related macular degeneration. Surv Ophthalmol. 1988, 32: 375-413. 10.1016/0039-6257(88)90052-5.View ArticlePubMedGoogle Scholar
- Klein R, Wang Q, Klein BE, Moss SE, Meuer SM: The relationship of age-related maculopathy, cataract and glaucoma to visual acuity. Invest Ophthalmol Vis Sci. 1995, 36: 182-191.PubMedGoogle Scholar
- Leibowitz HM, Krueger DE, Maunder LR, Milton RC, Kini MM, Kahn HA, Nickerson RJ, Pool J, Colton TL, Ganley JP, Loewenstein JI, Dawber TR: The Framingham Eye Study monograph: an ophthalmological and epidemiological study of cataract, glaucoma, diabetic retinopathy, macular degeneration, and visual acuity in a general population of 2631 adults, 1973–1975. Surv Ophthalmol. 1980, 24 (Suppl): 335-610.PubMedGoogle Scholar
- Wormald R: Assessing the prevalence of eye disease in the community. Eye. 1995, 9: 674-6.View ArticlePubMedGoogle Scholar
- Holz FG, Wolfensberger TJ, Piguet B, Gross-Jendroska M, Wells JA, Minassian DC, Chisholm IH, Bird AC: Bilateral drusen in age-related macular degeneration – prognosis and risk factors. Ophthalmology. 1994, 101: 1522-1528.View ArticlePubMedGoogle Scholar
- Maguire P, Vine AK: Geographic atrophy of the retinal pigment epithelium. Am J Ophthalmol. 1986, 102: 621-625.View ArticlePubMedGoogle Scholar
- Sunness JS, Bressler NM, Tian Y, Alexander J, Applegate CA: Measuring geographic atrophy in advanced age-related macular degeneration. Invest Ophththalmol Vis Sci. 1999, 40: 1761-1769.Google Scholar
- Rückmann AV, Fitzke FW, Bird AC: Distribution of fundus autofluorescence with a scanning laser ophthalmoscope. BR J Ophthalmol. 1995, 79: 407-412.View ArticleGoogle Scholar
- Rückmann AV, Fitzke FW, Bird AC: Fundus autofluorescence in age-related macular disease imaged with a laser scanning ophthalmoscope. Invest Ophthalmol Vis Sci. 1997, 38: 478-486.Google Scholar
- Solbach U, Keilhauer C, Knabben H, Wolf S: Imaging of retinal autofluorescence in patients with age-related macular degeneration. Retina. 1997, 17: 385-389.View ArticlePubMedGoogle Scholar
- Holz FG, Bellmann C, Margaritidis M, Schutt F, Otto TP, Volcker HE: Patterns of increased in vivo fundus autofluorescence in the junctional zone of geographic atrophy of the retinal pigment epithelium associated with age-related macular degeneration. Graefes Arch Clin Exp Ophthalmol. 1999, 237: 145-152. 10.1007/s004170050209.View ArticlePubMedGoogle Scholar
- Rückmann AV, Fitzke FW, Bird AC: In vivo fundus autofluorescence in macular dystrophies. Arch Ophthalmol. 1997, 115: 609-615.View ArticleGoogle Scholar
- Schmitz-Valckenberg S, Jorzik J, Unnebrink C, Holz FG: Analysis of digital scanning laser ophthalmoscopy fundus autofluorescence images of geographic atrophy in advanced age-related macular degeneration. Graefes Arch Clin Exp Ophthalmol. 2002, 240: 73-78.View ArticlePubMedGoogle Scholar
- Steger C: An unbiased detector of curvilinear structures. Technical Report FGBV-96-03, Technical University Munich, research team image processing. 1996Google Scholar
- Bland JM, Altman DG: Measuring agreement in method comparison studies. Stat Methods Med Res. 1986, 8: 135-160. 10.1191/096228099673819272.View ArticleGoogle Scholar
- Sachs L: Angewandte Statistik. 2004, Berlin: Springer-Verlag, 11View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2415/5/8/prepub