|
|
||||||||
From the Glaucoma Research Unit, Moorfields Eye Hospital, London, United Kingdom.
| Abstract |
|---|
|
|
|---|
METHODS. The experimental reference plane was positioned so that (1) it always lay entirely below the margin of the optic nerve head (ONH), (2) it remained at a set z-axis distance below the ONH in images of each eye, and (3) it was at a level where variability in rim area is least. Twenty normal control subjects and 20 patients with glaucoma underwent test-retest scanning laser tomographic imaging by same and different operators during same and separate visits. Control subjects had image series spanning at least 3 years. The effect of the positioning of the reference plane on global and regional rim area variability was assessed in intra- and intervisit test-retest images and longitudinal image series and compared with the standard and 320-µm reference planes.
RESULTS. Variability in the experimental reference plane was less in test-retest images and longitudinal data (P < 0.05) and more uniform around the ONH than with other reference planes. Variability in the former was not appreciably affected by testing involving different operators and visits, or by the presence of glaucoma.
CONCLUSIONS. Variability in rim area by the experimental reference plane was significantly less, more uniform around the ONH, not affected by different operators and visits, and less affected by glaucomatous morphology than other reference planes. This difference was pronounced in sequential data and has implications for detecting progression of glaucoma.
Many different reference surfaces have been proposed, but it is not clear which is most useful for assessing progression. Heidelberg Retina Tomograph (HRT; Heidelberg Engineering, Heidelberg, Germany) software previously used a "curved" reference surface with a shape that varied with the height profile of the margin of the ONH.4 Flat reference planes have since been used5 and various landmarks proposed for their positioning, such as peripapillary sclera,6 7 retinal pigment epithelium,8 peripapillary retina,9 10 and retinal height at the margin of the ONH.5 9 However, it is unclear which reference plane is useful in longitudinal analysis, and hence the reproducibility of different reference planes requires study.
We have found that rim area variability differs with reference plane positioning11 and can be explained by fluctuation in height between the ONH and reference plane.12 To optimize reproducibility, we devised an experimental reference plane that maintains an unchanging height relationship with the ONH throughout an imaging series. The reproducibility of this reference plane was evaluated and compared with the HRTs standard reference plane and 320-µm reference plane.
| Methods |
|---|
|
|
|---|
|
LOW5%.
The contour lines lowest region was calculated (LOW5%). Heights 1° apart (360 data points) on the contour lines of single topography images (used to derive the baseline mean topography image) were read from the HRT software. Heights were ranked in each single topography image, from which the mean of the lowest 5% of heights was calculated. The means of the lowest 5% of heights for each of the three single topography images were averaged to arrive at LOW5%.
R.
The reference plane was positioned to ensure that it lay beneath the entire circumference of the contour line. A reference plane lying above part of the contour line underestimates adjacent rim tissue. To determine R, variability was analyzed in longitudinal image series of normal control eyes. R was the distance of the reference plane beneath LOW5%, where rim area variability was least. This was held constant once calculated.
REFdis.
Once positioned at R, the z-axis distance of the reference plane below mean height of the contour line (MHC; which is the mean height of locations on the contour line measured in relation to the reference ring) was calculated for each ONH and called REFdis. MHC was used as a marker of the topographical z-axis position of the ONH. REFdis is unique to each ONH and, once calculated in the baseline image, is kept constant in all the images of an eye. By keeping REFdis constant, we sought to keep the height relationship between the ONH and reference plane constant in image series of the same eye. REFdis can be expressed as:
![]() | (1) |
REFpos.
Position of the reference plane, or REFpos, can thus be expressed as:
![]() | (2) |
![]() | (3) |
Standard Reference Plane.
This plane was set 50 µm posterior to the mean contour line height between 350° and 356° on the contour line (HRT software versions 1.112.01, and HRT II).5 9
The 320-µm Reference Plane.
This plane was offset 320 µm posterior to the mean height of the "reference ring" (HRT software versions 1.091.10). The reference ring is centered on the image frame and located in its periphery, had an outer diameter 94% and a width 3% of the image size, and is the zero-referencing region for pixel height.9
Subjects and Imaging
Variability in rim area was studied in test-retest images of age-matched subjects with glaucoma and normal control subjects, and in longitudinal series of images of normal eyes. Normal subjects and those with glaucoma attended the Ocular Hypertension and Early Glaucoma Research Clinic at Moorfields Eye Hospital and underwent the same protocol of repeat testing. All had previously undergone scanning laser tomography and perimetry at least six times. Table 1 shows subjects demographic information. This study adhered to the tenets of the Declaration of Helsinki and had appropriate Institutional Review Board approval and the subjects informed consent.
|
Test-retest scanning laser tomography (HRT software ver. 2.01; Heidelberg Engineering) of both eyes of 20 normal and 20 age-matched subjects with glaucoma was conducted by experienced operators. Three well-centered 10° single topography images were acquired at each session. Corneal curvature, scan depth, and focus settings were kept constant. Pupils were not dilated. Normal subjects attended two test visits separated by 6 to 8 months. Each visit comprised two imaging sessions separated by at least 1 hour. The same operator performed the scans in both imaging sessions of the first visit and in one session of the second visit. A second operator performed the scans in the other imaging sessions of the second visit. Patients with glaucoma underwent test-retest imaging in only one visit, performed by the same operator. Imaging sessions were randomly scheduled. Test-retest variability was analyzed for testing by: the same operator in the same visit (same operator-same visit) and by different operators in separate visits (different-operator-different-visit). This is relevant because longitudinal imaging in chronic glaucoma is likely to be performed in many different visits over time by different operators.
Nineteen normal subjects had longitudinal image series from at least six imaging visits over at least 3 years. Different operators with various levels of experience had acquired these images because of the turnover of research clinic technical staff. Each series was analyzed separately.
Image Analysis and Statistics
Mean topography images from one randomly selected eye of subjects in the glaucoma and normal groups were analyzed by the three reference planes. Mean images were generated from triplets of single topography images and used if mean pixel standard deviation was less than 50 µm, with pixel height measured by the reference ring. A contour line was outlined corresponding to the inner margin of the scleral ring of Elschnig on a mean topographic image of each subject (all performed by JCHT). Stereoscopic optic disc photographs were referred to if needed. Contour lines were then exported to related test-retest mean topography images. Rim area was evaluated globally and regionally. Patterns of variability around the ONH were examined by assessing regional rim area in 30° sectors around the ONH (0360°).
Variability in test-retest images was represented by agreement14 between pairs of images, expressed in graphs as the width of agreement intervals (95% confidence intervals of differences). For regional variability, analysis was modified so that agreement intervals for rim area sectors were presented in bar graphs by angle (0360°).
Variability in longitudinal image series was represented by standard deviation, which was used to estimate variability so that R could be derived and in longitudinal data. The width of the intervals for the 5th to 95th percentiles and median standard deviations for each sector were plotted in bar graphs. Significance testing in nonparametric data was by the Wilcoxon matched-pairs test (signed rank sum test) for comparing paired data. Statistical analysis was conducted by computer (SPSS ver. 9.0 for Windows; SPSS Inc, Chicago, IL).
| Results |
|---|
|
|
|---|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
No two ONHs are identical, and thus analysis should adapt to suit morphologic variations. The experimental reference plane was devised to lie entirely beneath the margin of the ONH by factoring the height profile of each contour line into calculations. Figure 7 illustrates why this is necessary. LOW5% was the level in each ONH above which the reference plane should not rise. The value for LOW5% was not preselected but was derived from the lowest 5% of contour line heights of each ONH and averaged from multiple topographies to minimize the effect of randomly outlying values. R of 100 µm (100 µm below LOW5%) was relatively deep and compatible with least variability. With this R, the reference plane is more likely to remain entirely below the ONH margin despite topographical variability. The cup is probably steeper here as well. Rim area is affected more when the reference plane shifts on a gradually sloping than a steep cup, as suggested by Figures 1 and 7 . With R deeper than 120 µm, variability tended to increase, possibly because of proximity to major vascular trunks and the lamina cribrosa, especially in shallow cups. REFpos was derived from all heights on the contour line so that the reference planes position would not rely on localized regions or presume selectivity in the pattern of change.
|
Reference planes have usually been positioned by landmarks presumed to be relatively unaffected by glaucoma, although whether this is true has been debated.17 The experimental reference plane differs in being positioned by a landmark that is expected to change in glaucoma. But there are reasons that this position is useful. First, experimental reference-plane data were more reproducible than other reference planes and not appreciably influenced by varied test conditions or glaucomatous morphology. Variability that changes with repeat testing and advancing damage can confound analysis. Second, the experimental reference plane compensated for shifts of the ONH that can affect analysis were the reference plane not to shift in tandem, as illustrated in Figure 8 . Changed rim area should reflect disease rather than such shifts. Third, positioning the reference plane at the same distance below the ONH surface allowed standardized measurement of rim area in imaging series. The ONH may shift in relation to external landmarks, such as the reference ring, which if used to anchor the reference plane can result in measurement artifact (Fig. 8) . Bias in measuring height and volume by the experimental reference plane can be expected, but this is less so with rim area because it is estimated perpendicular to the z-axis. Fourth, detecting progression depends on the systems signal-to-noise ratio. The z-axis position of the ONH fluctuated considerably in image series: 50% had a range of MHC exceeding 120 µm, with this range greater than 200 µm in some. Normal human retinal nerve fiber layer (RNFL) at the ONH margin is 316 to 406 µm thick,18 and so MHC variability is considerable in relation to normal RNFL thickness. MHC variability may even exceed RNFL thickness in glaucoma if significant axons are already lost by the time visual field defects appear. Empiric data suggest that the benefit of significantly reduced variability should outweigh the possibility that MHC is affected by glaucoma. Figure 8 shows that experimental reference plane analysis is sensitive to subtle rim area progression despite markedly reduced MHC.
|
| Footnotes |
|---|
Submitted for publication January 14, 2002; revised April 29, 2002; accepted June 10, 2002.
Commercial relationships policy: N.
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: Roger A. Hitchings, Glaucoma Research Unit, Moorfields Eye Hospital, City Road, London EC1V 2PD, UK; roger.hitchings{at}virgin.net.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. G. Strouthidis, H. Yang, B. Fortune, J. C. Downs, and C. F. Burgoyne Detection of Optic Nerve Head Neural Canal Opening within Histomorphometric and Spectral Domain Optical Coherence Tomography Data Sets Invest. Ophthalmol. Vis. Sci., January 1, 2009; 50(1): 214 - 223. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Poli, N. G. Strouthidis, T. A. Ho, and D. F. Garway-Heath Analysis of HRT Images: Comparison of Reference Planes Invest. Ophthalmol. Vis. Sci., September 1, 2008; 49(9): 3970 - 3975. [Abstract] [Full Text] [PDF] |
||||
![]() |
N G Strouthidis, E T White, V M F Owen, T A Ho, and D F Garway-Heath Improving the repeatability of Heidelberg retina tomograph and Heidelberg retina tomograph II rim area measurements Br. J. Ophthalmol., November 1, 2005; 89(11): 1433 - 1437. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. M. Zangwill, R. N. Weinreb, J. A. Beiser, C. C. Berry, G. A. Cioffi, A. L. Coleman, G. Trick, J. M. Liebmann, J. D. Brandt, J. R. Piltz-Seymour, et al. Baseline Topographic Optic Disc Measurements Are Associated With the Development of Primary Open-Angle Glaucoma: The Confocal Scanning Laser Ophthalmoscopy Ancillary Study to the Ocular Hypertension Treatment Study Arch Ophthalmol, September 1, 2005; 123(9): 1188 - 1197. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. K.-s. Leung, W.-m. Chan, Y.-l. Hui, W.-h. Yung, J. Woo, M.-k. Tsang, and K.-k. Tse Analysis of Retinal Nerve Fiber Layer and Optic Nerve Head in Glaucoma with Different Reference Plane Offsets, Using Optical Coherence Tomography Invest. Ophthalmol. Vis. Sci., March 1, 2005; 46(3): 891 - 899. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. H. Tan, D. Poinoosawmy, and R. A. Hitchings Tomographic Identification of Neuroretinal Rim Loss in High-Pressure, Normal-Pressure, and Suspected Glaucoma Invest. Ophthalmol. Vis. Sci., July 1, 2004; 45(7): 2279 - 2285. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. H. Tan and R. A. Hitchings Optimizing and Validating an Approach for Identifying Glaucomatous Change in Optic Nerve Topography Invest. Ophthalmol. Vis. Sci., May 1, 2004; 45(5): 1396 - 1403. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. C. H. Tan and R. A. Hitchings Approach for Identifying Glaucomatous Optic Nerve Progression by Scanning Laser Tomography Invest. Ophthalmol. Vis. Sci., June 1, 2003; 44(6): 2621 - 2626. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |