|
|
||||||||
From the Departments of 1 Ophthalmology, 2 Physiology and Biophysics, and 3 Mathematics and Statistics, Dalhousie University, Halifax, Nova Scotia, Canada.
| Abstract |
|---|
|
|
|---|
METHODS. The 256 x 256-pixel array of topographic height values obtained with each image from the Heidelberg Retina Tomograph (Heidelberg Engineering, Heidelberg, Germany) was divided into an array of 64 x 64 superpixels, where each superpixel contained 16 (i.e., 4 x 4) pixels. An analysis of variance technique was developed to analyze each superpixel with three baseline and three follow-up images. The performance of the technique was tested with and without adjustment for spatial correlation of topographic values using computer simulations and with real data from a normal control subject and a patient with progressive glaucomatous disc change.
RESULTS. Computer simulation with fixed population means and variance, and varying spatial correlation showed a monotonically increasing number of superpixels with significant test results (false positives), with 20% false-positives when the spatial correlation was 0.8 (the approximate median value in real patient data). The number of false-positive results was similar (17%) in serial images of a normal subject. When corrected for spatial correlation, the number of false-positives was independent of the level of spatial correlation and remained at the expected value of less than 5% in both simulations and real data. Although the number of significant test results in the patient with progressive glaucoma decreased after correction for spatial correlation, the change was readily apparent. Statistical power to detect mean differences in topographic values ranging from 0.5 to 4.0 SDs in computer simulation showed low power for changes of 1 SD or less, but increased dramatically with larger changes.
CONCLUSIONS. This technique has a high level of sensitivity to detect changes in the optic disc while maintaining a high level of specificity.
| Introduction |
|---|
|
|
|---|
Confocal scanning laser tomography has been introduced as an adjunct or alternative method for the clinical evaluation of the optic disc. The technique has been described elsewhere,6 7 but briefly, confocal sections of the optic disc, where the focal plane of the laser and the detector plane are optically conjugate, are obtained. The focal plane of the laser is changed incrementally to obtain the sections. The optical setup ensures that the information contained in a given image section is derived largely from the focal plane of the laser. After the confocal sections are aligned and processed, topographic heights of discrete locations in the scanned area are estimated.
Confocal scanning laser tomography allows reproducible estimates of the topographic measurements at individual measured locations, or pixels.8 9 10 The most important advantages of the technique are ease of operation, rapid image acquisition and processing times, and, unlike conventional photography, the ability to obtain images with natural pupils in most patients.11
The ability of detecting change due to disease depends largely on the testretest variability of measurements. When the variability of measurements is high, there is little statistical confidence in detecting small changes over time. If, however, the variability is low, small changes can be detected with confidence. Previous studies with confocal scanning laser tomography have shown that local variability measurements depend critically on topographic gradients and whether the measurements are made on blood vessels.10 12 Because these anatomic features are unique to each eye, local variability estimates should be made for each eye to gauge whether corresponding local topographic differences between two sets of images separated by time are statistically significant.
The purpose of this study was to describe a statistical technique for detecting changes in the optic disc and the peripapillary retina using a commercially available device, the Heidelberg Retina Tomograph (Heidelberg Engineering, Heidelberg, Germany). The robustness of the technique was tested with computer simulation as well as with real follow-up data from a glaucoma patient with progressive disc change and a healthy control subject.
| Materials and Methods |
|---|
|
|
|---|
The 32 sections in each scan are aligned for horizontal and vertical shifts to compensate for any eye movements during image acquisition. The reflectivity profile for each aligned pixel is determined by plotting reflectivity as a function of distance along the z-axis. For a given pixel, the area under the reflectivity profile is the sum reflectivity, and the position of maximum reflectivity along the z-axis is assumed to be the topographic height value. After the calculations have been made for all pixels, reflectivity and topographic images of the scanned area are determined. Therefore, the result of each scan is a topographic image representing the topographic height of each pixel from the focal plane of the eye.
Typically, in each session multiple scans are obtained in a subject (usually three9 ) from which a mean topographic and reflectivity image are computed after horizontal, vertical, and rotational alignments are made with a correlation procedure using the reflectivity images. Alignment for depth and tilt are made using the topographic images.
Statistical Analysis
The 256 x 256-pixel array of each topographic image was
divided into a 64 x 64-array of superpixels with each superpixel
containing 16 (i.e., 4 x 4) topographic values. Pooling over a
larger area allows more reliable estimates of testretest variability
using three baseline scans.10
In this case, there are 48
measurements (3 x 4 x 4) per superpixel. Because the size
of a superpixel in a 10° x 10°-scan is approximately 47 x 47
µm, variability estimates based on the 48 measurements will be
influenced by the topography of the imaged structure in the superpixel.
Therefore, although the pooling has minimal effect on a superpixel
situated in an area with flat topography, in an area with steep
contours such as the optic cup edge, the variability estimates are
increased. To remove the topographic component in variability
estimates, we subtracted the respective topographic measurement in each
pixel from its respective mean across the three images to determine the
adjusted value. After this process, we computed three corrected images
from which estimates of testretest variability were made. These
estimates are expressed as variability maps and have been described in
detail elsewhere.10
Assuming that there are I baseline and follow-up images and
that the topographic measurements from a superpixel of 4 x 4
pixels is a vector of length 16 indexed by l, an analysis of
variance model for each superpixel in vector form is:
![]() | (1) |
e2.
With respect to equation 1
, determining the significance of temporal
changes in topographic values within each superpixel corresponds to
testing for the main effect of time and the interaction of time and
location simultaneously leading to the test statistic
![]() | (2) |
; where
SS{T} is the sum of squares associated with time;
SS{TL} is the sum of squares associated with time by
location interaction; SS{I(T)} is the sum of squares
associated with image within time; SS{e} is the sum
of squares associated with the residual or error; and degrees of
freedom
= 2(I - 1)16.
If Tt,
TLtl, and
I(T)it are assumed
to be fixed effects, then the statistic in equation 2
has an
F distribution with 16 and
degrees of freedom. (See
Appendix A for a detailed description of the test statistic).
Given that topographic values in neighboring pixels are likely to
be correlated, accounting has to be made for this spatial dependence.
If I(T)it is a
random variable with a mean of zero and variance
I(T)2 that
are independent of etli, then the model
(equation 1) allows for spatial correlation between topographic values
within a superpixel. The correlation between
htli and
htmi (two locations within the same
superpixel) is
![]() |
measures the degree of spatial dependence
between topographic values within a superpixel.
The effect of spatial correlation is to decrease the amount of
available information. The Satterthwaite correction13
appropriately reduces the degrees of freedom of the approximating
F distribution for the test statistic to
![]() | (3) |
![]() |
, and are 16 and
for
= 0. The effect
of the reduced degrees of freedom is to increase P (reduce
significance).
Computer Simulation and Case Reports
The specificity and sensitivity of the procedure described above
was tested using computer simulation and illustrated in case reports of
follow-up images of a healthy subject and a patient with progressing
glaucoma. The study protocol for the subjects was approved by the
Research Ethics Committee of the Queen Elizabeth II Health Sciences
Center. Informed consent was obtained, and the tenets of the
Declaration of Helsinki were followed.
Specificity estimates were made by varying the degrees of spatial
dependence in baseline and follow-up images in computer simulations.
Simulated baseline and follow-up images were randomly generated by
keeping the population mean and variance of topographic values within a
corresponding aligned superpixel identical at baseline and follow-up.
Ten thousand simulations were performed for each level of spatial
correlation ranging from 0 (no spatial correlation) to 0.99 (almost
perfect spatial correlation), and the number of superpixels showing
significant change (false positives) at
= 0.05 was recorded.
We compared the results with and without adjustment for spatial
correlation.
Sensitivity estimates were made by randomly generating baseline and
follow-up images where the mean change between corresponding aligned
superpixels varied from 0.5 to 4 SDs, whereas the variance was kept
identical. Ten thousand simulations were performed, for the same levels
of spatial correlation used for the specificity estimates. The number
of significant differences that were detected at
= 0.05 were
recorded.
| Results |
|---|
|
|
|---|
|
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
In this study we have described a new method for the serial analysis of topographic images based on an empiric probabilistic approach. Variability estimates for small discrete areas (superpixels), each comprising a square of 4 x 4 pixels, are estimated. An analysis of variance is conducted on aligned superpixels based on a calculation of serial pointwise differences. The degrees of freedom value used to calculate the F statistic is adjusted to account for the spatial correlation between topographic values of pixels within a superpixel.
Our study showed that the correlation between the topographic values in neighboring pixels is considerable with a median intraclass correlation coefficient of 0.8, which, if not corrected for, would yield erroneous results in the follow-up of patients and healthy control subjects. Computer simulations using baseline and follow-up images with the same population mean and variance showed that as the spatial correlation increased, the number of false-positive test results also increased, with more than 20% false-positive results when the spatial correlation approached 0.8. Analyzing serial images of a normal subject, we found a similar number of false-positive results. When we applied the Satterthwaite correction to adjust the test statistic, the number of false-positive results was invariant to increases in spatial correlation in computer simulation and reduced the number of false-positive results in the same normal subject to levels expected by chance alone.
Computer simulation also allowed us to determine the power of our proposed method to detect changes as a function of change in topographic values between serial images and spatial correlation. For changes of 2 SDs or more, there was adequate power to detect changes in superpixels where the spatial correlation is 0.6 or less; however, increasingly larger changes are required in superpixels that have a high degree of spatial correlation. The correction for spatial correlation is likely not to inhibit the ability to detect clinically meaningful glaucomatous disc change.
Because scanning laser tomography can be performed rapidly and without pupil dilation in most patients, it is feasible to obtain optic disc images at each clinical visit. The availability of a series of serial images would enhance the ability of the clinician to assess optic disc progression. A form of serial change analysis in which each follow-up image is compared to a baseline image is clearly possible. Using these data, it is possible to determine locations where repeatably significant change from baseline has occurred and assign probabilities to these changes, both in terms of size of change (number of clustered superpixels) and time (temporally overlapping significant superpixels).
We have described a statistical technique for the serial analysis of optic disc topography with scanning laser tomography with computer simulation and case examples. The validity and merit of the technique in a clinical situation, of course, require data from a larger number of subjects. Additionally, the guidelines for the clinical use of this and other techniques for determining change require some convergence between statistical significance and clinical significance. The performance of this technique in a prospective longitudinal cohort of glaucoma patients and normal control subjects and comparing results with the topographic indices, conventional optic disc photography, and visual field progression are beyond the scope of this article and will the subject of a future publication.
| Appendix A |
|---|
|
|
|---|
The test statistic is a ratio of two terms MS{NUM} and MS{DEN}.MS{NUM} is a scaled measure of the total squared deviation between the estimated mean heights at the two times, whereas MS{DEN} is a scaled measure of the total squared deviation of individual heights around their means. Under the null hypothesis of no mean differences among the images at each time and of no mean difference between the two timesthat is, no T, TL, or I(T) effectsthe F ratio is approximately equal to one. Large values of F give evidence against the null hypothesis of no difference and suggest that the differences are real.
Assuming no differences between images at time t, the mean
heights are µtl = µ +
Tt + Ll +
TLtl. These means are estimated by the
average heights at each location l
![]() |
![]() |
![]() |
2/I, and
SS{NUM}/(2
2/I) has a
2 distribution with 16 degrees of freedom. The
numerator of the test statistic is
![]() |
Under the same assumptions, the deviation
htli -
tl
estimates the error term etli and their
sum of squares
![]() |

2, where
=
2(I - 1)16, so an estimate of
2 is given by
![]() |
2 has
a
2 distribution with
degrees of freedom,
independent of SS{NUM}. It can be shown that
![]() |
It follows that, under the assumptions, the test statistic (equation 2) has an F distribution with 16 and
degrees of freedom. If
I(T)it is assumed to be a
random effect, however, then the Satterthwaite approximation is used to
approximate the distribution of the test statistic as an F
with fN and
fD degrees of freedom, as in equation 3
.
| Appendix B |
|---|
|
|
|---|
![]() |
![]() |
where
is the degree of spatial correlation and
=
2(I - 1)16. If there is no spatial dependence
(
= 0), then the expected number of degrees of freedom are
(16,
), which agrees with the classic result. As
increases, the
number of degrees of freedom in both the numerator and denominator
decrease in a monotonic fashion. In the limit, as
approaches 1,
E[fN] approaches 1 and
E[fD] approaches
2(I - 1).
| Footnotes |
|---|
Submitted for publication July 7, 1999; revised September 8, 1999; accepted October 26, 1999.
Commercial relationships policy: N.
Corresponding author: Balwantray C. Chauhan, Department of Ophthalmology, Dalhousie University, 2nd Floor Centennial Building, Queen Elizabeth II Health Sciences Centre, Halifax, NS, Canada B3H 2Y9. bal{at}is.dal.ca
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. C. Chauhan, D. M. Hutchison, P. H. Artes, J. Caprioli, J. B. Jonas, R. P. LeBlanc, and M. T. Nicolela Optic Disc Progression in Glaucoma: Comparison of Confocal Scanning Laser Tomography to Optic Disc Photographs in a Prospective Study Invest. Ophthalmol. Vis. Sci., April 1, 2009; 50(4): 1682 - 1691. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Bowd, M. Balasubramanian, R. N. Weinreb, G. Vizzeri, L. M. Alencar, N. O'Leary, P. A. Sample, and L. M. Zangwill Performance of Confocal Scanning Laser Tomograph Topographic Change Analysis (TCA) for Assessing Glaucomatous Progression Invest. Ophthalmol. Vis. Sci., February 1, 2009; 50(2): 691 - 701. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Fortune, G. A. Cull, and C. F. Burgoyne Relative Course of Retinal Nerve Fiber Layer Birefringence and Thickness and Retinal Function Changes after Optic Nerve Transection Invest. Ophthalmol. Vis. Sci., October 1, 2008; 49(10): 4444 - 4452. [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] |
||||
![]() |
H. Yang, J. C. Downs, A. Bellezza, H. Thompson, and C. F. Burgoyne 3-D Histomorphometry of the Normal and Early Glaucomatous Monkey Optic Nerve Head: Prelaminar Neural Tissues and Cupping Invest. Ophthalmol. Vis. Sci., November 1, 2007; 48(11): 5068 - 5084. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. DeLeon Ortega, L. M. Sakata, B. Kakati, G. McGwin Jr, B. E. Monheit, S. N. Arthur, and C. A. Girkin Effect of Glaucomatous Damage on Repeatability of Confocal Scanning Laser Ophthalmoscope, Scanning Laser Polarimetry, and Optical Coherence Tomography Invest. Ophthalmol. Vis. Sci., March 1, 2007; 48(3): 1156 - 1163. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. M. F. Owen, N. G. Strouthidis, D. F. Garway-Heath, and D. P. Crabb Measurement Variability in Heidelberg Retina Tomograph Imaging of Neuroretinal Rim Area Invest. Ophthalmol. Vis. Sci., December 1, 2006; 47(12): 5322 - 5330. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Patterson, D. F. Garway-Heath, and D. P. Crabb Improving the repeatability of topographic height measurements in confocal scanning laser imaging using maximum-likelihood deconvolution. Invest. Ophthalmol. Vis. Sci., October 1, 2006; 47(10): 4415 - 4421. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. G. Strouthidis, A. Scott, N. M. Peter, and D. F. Garway-Heath Optic disc and visual field progression in ocular hypertensive subjects: detection rates, specificity, and agreement. Invest. Ophthalmol. Vis. Sci., July 1, 2006; 47(7): 2904 - 2910. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Shunmugam and A. Azuara-Blanco The quality of reporting of diagnostic accuracy studies in glaucoma using the heidelberg retina tomograph. Invest. Ophthalmol. Vis. Sci., June 1, 2006; 47(6): 2317 - 2323. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. T. Nicolela, A. S. Soares, M. M. Carrillo, B. C. Chauhan, R. P. LeBlanc, and P. H. Artes Effect of moderate intraocular pressure changes on topographic measurements with confocal scanning laser tomography in patients with glaucoma. Arch Ophthalmol, May 1, 2006; 124(5): 633 - 640. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Patterson, D. F. Garway-Heath, N. G. Strouthidis, and D. P. Crabb A New Statistical Approach for Quantifying Change in Series of Retinal and Optic Nerve Head Topography Images Invest. Ophthalmol. Vis. Sci., May 1, 2005; 46(5): 1659 - 1667. [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] |
||||
![]() |
J. C. H. Tan and R. A. Hitchings Reference Plane Definition and Reproducibility in Optic Nerve Head Images Invest. Ophthalmol. Vis. Sci., March 1, 2003; 44(3): 1132 - 1137. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. H. Morgan, B. C. Chauhan, D.-Y. Yu, S. J. Cringle, V. A. Alder, and P. H. House Optic Disc Movement with Variations in Intraocular and Cerebrospinal Fluid Pressure Invest. Ophthalmol. Vis. Sci., October 1, 2002; 43(10): 3236 - 3242. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. C. Chauhan, T. A. McCormick, M. T. Nicolela, and R. P. LeBlanc Optic Disc and Visual Field Changes in a Prospective Longitudinal Study of Patients With Glaucoma: Comparison of Scanning Laser Tomography With Conventional Perimetry and Optic Disc Photography Arch Ophthalmol, October 1, 2001; 119(10): 1492 - 1499. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |