|
|
||||||||
1From the Glaucoma Research Unit, Moorfields Eye Hospital, London, United Kingdom; 2Department of Information Systems and Computing, Brunel University, London, United Kingdom; 3Discoveries In Sight, Devers Eye Institute, Portland, Oregon; and the 4Department of Optometry and Visual Science, City University, London, United Kingdom.
| Abstract |
|---|
|
|
|---|
METHODS. Previously, interpoint correlations were generated for all possible pairs of VF test points in a dataset of 98,821 Humphrey VF test results taken from the Moorfields Eye Hospital archive. The relationship between these correlations and the physical distance between the VF test point pairs was evaluated by Pearsons correlation coefficient and multiple regression analysis. The distance between the pairs of VF test points was calculated in two ways. First, the anatomic map was used to estimate the angular distance at the optic nerve head (ONH), between the RNFL bundles corresponding to the VF test points in each pair (ONHd). Second, the retinal distance between pairs of test points was calculated from the Humphrey VF template (RETd). A best-fit model for predicting functional correlation (FC) from ONHd and RETd was constructed and used to formulate a filter incorporating the anatomic-functional correlation data.
RESULTS. All scatterplots showed a negative association between interpoint retinal sensitivity correlation values and distance between points: ONHd (R2 = 0.60) and RETd (R2 = 0.33). The raw sensitivity correlation values could be predicted from a multiple regression model using ONHd, RETd, and a combined interaction of ONHd and RETd (R2 = 0.75, P < 0.00001). The construction of a new filter was based on the equation FC = 0.9325 (0.0029 · ONHd) (0.0077 · RETd) + (0.0001 · ONHd · RETd).
CONCLUSIONS. A good level of association was observed between the strength of correlation between points in the VF and the relative location of those test points in the peripheral retina and in corresponding RNFL bundles at the ONH. These results help to validate the relationship between structure and function and may be of use in the further refinement of physiologically derived VF filters to reduce measurement noise.
One approach to reducing measurement variability is the post hoc application of a spatial filter. Spatial filtering is a technique adapted from digital imaging processing whereby the measured sensitivity of a particular test point is adjusted based on the sensitivity of other points in the vicinity. A novel spatial filter has been designed with the intention that it closely mimics the physiological relationship between VF test points.8 This filter predicts the sensitivity of each test point based on the sensitivities at other locations in the field, weighted according to the predictive power of each of those locations. Predictive power was determined by examining the correlations and covariances between sensitivities among all pairs of test locations, within a large database of 98,821 predominantly glaucomatous VFs. The derivation of this filter has generated an array of interpoint absolute correlations for the entire field. These values effectively constitute a mathematically derived "functional map" based on physiological data. The relationship between points demonstrated by the filter should therefore closely mirror the structural pattern of the nerve fiber layer, although it may be influenced by idiosyncrasies of the Humphrey Visual Field Analyzer (HFA; Carl Zeiss Meditec, Inc., Dublin, CA) testing algorithm used.
In this study, we examined the relationship between this functional map and an anatomic map relating the RNFL distribution to the ONH.9 The purpose of the study was to identify whether the magnitude of functional correlation between points is related to the relative proximity of the points at the ONH and in the retinal periphery. This information may give insight to the anatomic organization of the ONH, the glaucomatous disease process, and enable the refinement of filters applied to the VF series to reduce measurement noise.
| Methods |
|---|
|
|
|---|
|
Description and Application of the Functional Map
The derivation of the novel spatial filter has been described in detail elsewhere.8 Briefly, the filter was derived from VFs from the Moorfields Eye Hospital database. This contains 98,821 VFs, taken from 14,675 individual patients with suspected glaucoma. The tests performed were all standard white-on-white, full-threshold tests; only complete 24-2 tests were used in the generation of the filter, although all levels of test reliability were included. The relationship between individual test points was elucidated by examining the covariances and correlations between the sensitivities of all test-point pairings throughout the database. These interpoint correlations (52 x 51 = 2652 in total) were used in this study as a gauge of "functional relationship" between points, and the basis of the comparison with the physical distances obtained from the anatomic map.
Comparison of the Anatomic Map and the Functional Map
Interpoint correlations were compared to RETd and ONHd using Pearson correlation coefficient and multiple regression analyses. Comparisons were performed for the entire template, for the upper and lower hemispheres and between hemispheres. Linear models for predicting VF correlations (FCs) from ONHd and RETd were assessed. The application of the prediction model to derive a "physiological" filter was illustrated for one VF test point. First, FCs predicted for that test point were generated; R2 was generated by squaring the predicted FCs. An R2 cutoff of 0.7 was selected to identify which test point pairings should be included in the filter for that test point. All statistical analyses were performed using S (AT & T Bell Laboratories, Murray Hill, NJ). This study adhered to the tenets of the Declaration of Helsinki.
| Results |
|---|
|
|
|---|
|
|
|
|
|
|
|
![]() |
A new filter was derived for test point 49 by using this predictive model; its derivation is illustrated in Figure 7 . The ONHd (Fig. 7a) and RETd (Fig. 7b) for point 49 were calculated. Using the regression equation, we predicted FCs for each test point in relation to point 49. The FC results were squared to generate R2 (Fig. 7c) . The test points included in the filter were those with R2
0.7 (highlighted in bold in Fig. 7c ). Finally, the weightings by which each test points sensitivity influences the sensitivity of point 49 were calculated by dividing each R2 value by the sum of all R2 values included in the filter (Fig. 7d) .
|
| Discussion |
|---|
|
|
|---|
A more useful approach to constructing a spatial filter would be to exploit the functional or anatomic relationship between test points. Initially, the point by point spatial dependence was determined by multiple regression analysis of sensitivity values for each test point in a dataset of 440 Humphrey VFs.13 A similar investigation reports the relationship between sensitivities of test points using the 32 program of Octopus 1-2-3 (Interzeag, Schlieren-Zurich, Switzerland).14 In this study, linear regression analysis among each of the locations and the rest of the points in the field was performed. The methodology used in the construction of our filter was similar, although the mathematical relationship between sensitivities was assessed using covariances and correlations and the number of fields assessed was much larger. In particular, all available VF data were used in the construction of the filter, so as to be truly representative of a glaucoma clinic population. It may therefore not be suitable for use in normal subjects or subjects with nonglaucomatousfor example, neurologicalfield defects. With simulated progressing VF data, the novel filter was found to improve both specificity and sensitivity.8 When used on a 50-patient sample of longitudinal field data, the filter has been shown to reduce variability, and it does not reduce detection of loss by total deviation maps (Artes PH et al., IOVS 2005;46:ARVO E-Abstract 3732). This method represents a clear improvement in the performance of the Gaussian filter, although the effect of the filter has yet to be fully assessed on prospective clinical data. An additional observation from this study is that the filter improved the "pattern" of progression compared with unfiltered VFs, so that it more closely resembled the defect appearance expected in glaucoma. This result would be expected if the physiological relationships exploited in the construction of the filter are valid.
An encouraging level of agreement between the magnitude of functional correlation between points and the relative location of the points at the ONH and the retinal periphery was observed in this study. There was a negative association between functional correlation and both ONHd and RETd. Using a multiple regression model with the product of ONHd and RETd, we were able to predict interpoint correlations. It should be noted that the model continues to predict FC well with the interaction term removed, which may indicate that the dependent association between ONHd/RETd and interpoint sensitivity correlation may not be large. However, although the coefficient for the interaction term is small (0.0001), it cannot be dismissed completely, as an increase in R2 from 0.6 to 0.75 was observed when the interaction was included. The interaction term is intended to account for the nonlinearity observed in the models shown in Figures 2 and 3 . The minimal impact of the interaction term on the predictive model may suggest that the nonlinearity observed has a negligible influence. The regression equation used to construct the example filter was therefore derived from a predictive model that included the interaction between ONHd and RETd. In the construction of the filter, only predicted FCs with ONHd correlations >0.84 (R2
0.7) were included, and the ONHd/FC relation is clearly linear over this range (Fig. 5) . However, as should be expected, the relationship between ONHd and the VF correlations appears to be wholly valid (and linear) within the same hemisphere but not between hemispheres (Fig. 6) .
As glaucomatous damage is believed to manifest at the ONH, it seems logical that VF locations that correspond to similar regions of the ONH should be well correlated; damage to that area of the ONH affects all such points. In this study, retinal proximity was also found to be a predictor of the strength of correlation between two points. This observation has implications for both disease process and anatomic organization, although with the caveat that it is unknown whether the finding is real or spurious. If the observation is "real," it may support the hypothesis that RNFL bundles from similar peripheral eccentricities are closely located at the optic nerve head. Experimental studies in different species of the macaque monkey have generally suggested that a degree of retinotopic organization exists with respect to the eccentricity of axonal origin, although they tend to differ in terms of exact detail.15 16 17 18 To date, there has been little by way of clinical observation to support this hypothesis,19 although the results of this study may support such a finding in the context of a disease model where glaucomatous damage occurs at the ONH. An alternative explanation applies to a model in which damage occurs primarily in the retina. In this situation, damage may propagate from dysfunctional or dying retinal ganglion cells locally within the retina. ONHd has a much higher coefficient of determination than RETd, suggesting that the glaucomatous process more likely occurs at the ONH, although the ONH and retinal models are not mutually exclusive. The FC/RETd relationship, however, may in part be spurious, resulting from measurement error. The error may be systematic, perhaps related to inaccuracies in the anatomic map. Random error may relate to interindividual variation of ONH position in relation to the fovea,9 or to fixation losses occurring during visual field testing.
The comparison between interpoint functional correlation and the anatomic map is dependent on the assumption that the relationships described by the anatomic map are valid. Alternative maps have been described that were developed with similar techniques.20 21 The map used in the present study has already been used in structure-function studies in glaucoma.22 23 24 The adoption of an alternative map, such as that developed by Junemann et al.22 has been based on the simplicity of use, as opposed to any perceived greater integrity compared to the map used in our study.25 Recently, a map has been described that was developed using both static automated perimetry and Heidelberg Retina Tomograph (HRT) data.26 This newer map therefore differs from the map used in our study, in that it incorporated both structural and functional information in its development, although the result is similar to the map used in the present study.
Structural data have been incorporated into the construction of a spatial filter, based on the multiple regression predictive model that incorporates the angular distance between test points at the optic nerve head, the angular distance between test points in the retinal periphery, and the interaction between the two distances. In the example used in this study, which is for a single test point, the filter has a similar distribution of test-point associations compared with those generated using the "physiological" filter of Gardiner et al.8 Both filters follow an "arcuate" pattern, in keeping with what might be expected given the distribution of the retinal nerve fiber layer. The newly developed "structural" filter does include fewer points, however, that have more similar weightings relative to each other, compared with the physiological filter. This may be explained by the fact that all the points, bar one, included in the structural filter are directly adjacent to the point of interest and as such may be expected to have a similar relation, according to a linear model. The method of constructing the physiological filter downplayed points that strongly covaried with, but had lower predictive value than, other predicting points. This method was not used in the new structural filter. By not downplaying strongly associated points, measurement noise reduction may be improved through increased signal averaging; however, whether this confers an advantage in the detection of signal should be tested and will be the subject of further work. The use of the predictive equation developed in this study enables the construction of filters that may be customized on a point-wise basis. A "bespoke" spatial filter is particularly useful if one wants to exclude test points from a longitudinal series if they are consistently depressed by a mechanism other than glaucomafor example, by lid artifact or chorioretinal scarring. Likewise, the physiological filter is limited as it is designed for use with the 24-2 program of the HFA. A point-wise customizable filter may be adopted in the context of different Humphrey programs (such as 10-2) and may also be used in alternative proprietary perimeters.
The associations identified in our study help to validate the structure-function relationship in glaucoma and give insight into the anatomic organization of the ONH and glaucomatous disease process. The incorporation of structural data may be of benefit in the development of more refined spatial filters to reduce measurement noise in VF testing.
| Footnotes |
|---|
Submitted for publication December 29, 2005; revised May 17 and August 23, 2006; accepted October 11, 2006.
Disclosure: N.G. Strouthidis, Heidelberg Engineering (F); V. Vinciotti, None; A.J. Tucker, None; S.K. Gardiner, None; D.P. Crabb, None; D.F. Garway-Heath, Carl Zeiss Meditec (C)
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: David F. Garway-Heath, Glaucoma Research Unit, Moorfields Eye Hospital, 162 City Road, London, EC1V 2PD, UK; david.garway-heath{at}moorfields.nhs.uk.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Wall, C. A. Johnson, R. H. Kardon, and D. P. Crabb Use of a Continuous Probability Scale to Display Visual Field Damage Arch Ophthalmol, June 1, 2009; 127(6): 749 - 756. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Ferreras, L. E. Pablo, D. F. Garway-Heath, P. Fogagnolo, and J. Garcia-Feijoo Mapping Standard Automated Perimetry to the Peripapillary Retinal Nerve Fiber Layer in Glaucoma Invest. Ophthalmol. Vis. Sci., July 1, 2008; 49(7): 3018 - 3025. [Abstract] [Full Text] [PDF] |
||||
![]() |
B C Chauhan, D F Garway-Heath, F J Goni, L Rossetti, B Bengtsson, A C Viswanathan, and A Heijl Practical recommendations for measuring rates of visual field change in glaucoma Br J Ophthalmol, April 1, 2008; 92(4): 569 - 573. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Cheng, M. Laron, J. S. Schiffman, R. A. Tang, and L. J. Frishman The Relationship between Visual Field and Retinal Nerve Fiber Layer Measurements in Patients with Multiple Sclerosis Invest. Ophthalmol. Vis. Sci., December 1, 2007; 48(12): 5798 - 5805. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |