|
|
||||||||
1From the Department of Visual Science, Institute of Ophthalmology, London, United Kingdom; 2Glaucoma Research Unit, Moorfields Eye Hospital, London, United Kingdom; 3St. Bartholomews Hospital, London, United Kingdom.
| Abstract |
|---|
|
|
|---|
METHODS. A computer-generated line stimulus was presented with displacements randomized between 0 and 40 minutes of arc (min arc). Displacement thresholds (50% seen) were compared for stimuli of equal area but different edge length (orthogonal to the direction of motion) at four retinal locations. Also, MDT thresholds were recorded at five values of Michelson contrast (25%84%) for each of five line lengths (11128 min arc) at a single nasal location (27,3). Frequency-of-seeing (FOS) curves were generated and displacement thresholds and interquartile ranges (IQR, 25%75% seen) determined by probit analysis.
RESULTS. Equivalent displacement thresholds were found for stimuli of equal area but half the edge length. Elevations of thresholds and IQR were demonstrated as line length and contrast were reduced. Equivalent displacement thresholds were also found for stimuli of equivalent energy (stimulus area x [stimulus luminance background luminance]), in accordance with Riccos law. There was a linear relationship (slope 0.5) between log MDT threshold and log stimulus energy.
CONCLUSIONS. Stimulus area, rather than edge length, determined displacement thresholds within the experimental conditions tested. MDT thresholds are linearly related to the square root of the total energy of the stimulus. A new law, the threshold energy-displacement (TED) law, is proposed to apply to MDT summation properties, giving the relationship T = K
E where, T is the MDT threshold, Kis the constant, and E is the stimulus energy.
A fundamental requisite for understanding the normal psychophysiological response to a stimulus is to select appropriate stimulus parameters to obtain good test sensitivity and specificity that will distinguish normal and diseased states. Specifically, in this study we investigated the spatial summation characteristics of the MDT stimulus.
We have previously demonstrated equivalent MDT thresholds for stimuli of equivalent line length, whether the stimulus comprises a single line or three shorter lines (independent of line separation up to 128 minutes arc [min arc]).10 The present study extends our earlier work and specifically explores summation properties of MDT stimuli in relation to line edge length, stimulus area, luminance, and contrast.
Scobey and Johnson11 12 investigated the effect of changing line length and luminance using a unidirectional stimulus. Elevation of the displacement threshold was found as line length and luminance were decreased, with an increase in threshold mean and SD with eccentricity. Similar findings were reported by MacVeigh et al.13 for an oscillating stimulus of low temporal frequency (2 Hz) in the presence of stationary references. The Moorfields MDT stimulus differs from these examples in that it undergoes three oscillations per presentation, which results in lower displacement thresholds when compared with a single oscillatory displacement.6 It is presented at 5 Hz in the absence of stationary references. Direct comparison of the methodologies is therefore difficult.
It has been suggested that vernier thresholds are determined by the stimulus edge length.14 15 16 17 18 19 20 In this study, we tested the hypothesis that it is the MDT stimulus area, rather than the edge length, that determines the threshold response, by comparing MDT thresholds for constant stimulus area but changing line edge length.
Previous studies have been performed to investigate vernier acuity as a function of contrast or line length.11 12 21 In the present study we took a novel approach and calculated displacement thresholds as a function of the stimulus energy (stimulus area x luminance) to explore the application of Riccos Law to MDT stimuli. Standard automated perimetry (SAP) measures differential light sensitivity (DLS) and, where the stimulus lies within Riccos critical area, the product of the area and intensity of the stimulus at threshold is constant.22 There are no publications establishing whether Riccos Law is applicable to an MDT-like stimulus, responses to which historically have been assumed to be mediated by the magnocellular (M) pathway.23
| Methods |
|---|
|
|
|---|
3%. Peripheral refractive error was estimated by retinoscopy in patients under mydriasis (tropicamide 1%). The subject fixated a central target, and the required corrective lenses were centered on a peripheral line stimulus, with an addition of +3.0 DSph to allow for the fixation distance of 30 cm. The line stimulus was programmed to pass through five presentations of each of seven randomized displacements (024 min arc). In each case, the vector of motion was orthogonal to the long edge of the line, which was aligned with the horizontal meridian. Each stimulus presentation consisted of three oscillations at 200 ms per cycle.7 The subject task was to press the computer mouse each time a displacement was seen. Each test paradigm was applied in a randomized order and tested three times, so that each displacement was presented 15 times for each stimulus condition. The duration of each test was approximately 2 minutes 30 seconds. The number of tests was limited to eight per session, with rest intervals of 5 minutes between tests.
Experiment 1: Variation of Stimulus Line Edge Length
The effect of halving the MDT stimulus line length, while maintaining equivalent area, was investigated for seven different configuration pairs of equal area.
Four retinal locations were tested using the specifications summarized in Table 1 . Two subjects were tested for each configuration, with five perimetrically experienced subjects taking part in the investigation (age: mean 37 years, range 2851).
|
|
|
|
|
|
|
|
Experiment 2: Additional Data Analysis.
The results of subjects 1 and 2 were averaged. The mean MDT threshold was plotted as a function of contrast for each line length on both linear and log coordinates. The interquartile range (IQR; difference between 75% and 25% of the probit fitted curve) was calculated to give an indication of the psychometric function slope.24 25 The relation between mean log MDT threshold and log IQR was explored by linear regression analysis. Stimulus energy was calculated relative to background luminance as follows: relative stimulus energy = stimulus area x (stimulus luminance background luminance). The background luminance was deducted from the stimulus luminance to account for the fact that, when the MDT stimulus is displaced, it moves to an area that is not black, but of constant background luminance (10 cd/m2). The relation between mean log MDT threshold and log relative stimulus energy was explored by linear regression analysis.
| Results |
|---|
|
|
|---|
Experiment 2
The FOS curves of each subject illustrate a similar pattern of shallowing of the slope and widening of the 95% CI and IQR, as contrast and line length is reduced, as shown by the selected examples (Fig. 2) . Note the steep curve of the longest line length (128 min arc) at the highest contrast (85%), with a low IQR and narrow confidence intervals (Fig. 2A) . Differences in psychometric function slopes remained when plotted on a logarithmic scale.
No responses were recorded for the smallest range of displacements (0 24 min arc) when testing the smaller line lengths at low contrast (subjects 1 and 2: 11, 22, and 43 min arc line lengths at 25% contrast and the 11 min arc line length at 40% contrast; and subject 1: 22 min arc line length at 40% contrast). The increase in the upper limit of displacement to 40 min arc enabled only one further FOS curve to be measured for each subject (subject 1: 22 min arc line length at 40% contrast; subject 2: 43 min arc line length at 25% contrast).
Analysis of MDT threshold as a function of contrast showed a curvilinear relationship for each line length, with an increase in threshold and widening of 95% CI, as contrast was reduced (Fig. 3) . This relationship became linear when plotted in loglog coordinates (Fig. 4 ; Table 2 ). There was a linear relationship (slope: 1.3) between log MDT threshold and log IQR (Fig. 5) . Log MDT thresholds were also linearly related to log relative stimulus energy, with a slope of approximately 0.5, indicating that the MDT threshold is related to the square root of the stimulus energy (Fig. 6) . The slope for each subject was subject 1, 0.46 (95% CI 0.51 to 0.41; R2 = 0.94, P < 0.001) and subject 2, 0.43 (95% CI 0.36 to 0.49; R2 = 0.88, P < 0.001).
|
| Discussion |
|---|
|
|
|---|
Displacement acuity is conventionally considered to be a function of edge jump, with summation calculations being made in relation to line edge length.14 19 The results of the first experiment (Fig. 1) suggest that it is the total area of the stimulus, as opposed to its edge length, that is the determining factor (for the range of parameters tested in the experiment). The wide CIs for the shorter but wider stimulus at retinal location 9,15 (Fig. 1C) may be explained by the stimulus movment across a blood vessel (Fig. 7) .
The second experiment shows a clear trend of decreasing MDT threshold with increasing line length and luminance (Fig. 3) , which is consistent with previously reported findings of summation experiments for line displacement stimuli.12 13 A re-analysis of 14 studies, which examined the relationship between contrast and vernier threshold as a power function, found loglog slopes ranging between 0.4 and 1.0.21 In this study, we found slopes ranging between 0.8 and 1.2 (Table 2) .
The IQR scales directly with the MDT threshold (Fig. 5) , with a slope of 1.26 (95% CI 0.98 to 1.54) on loglog coordinates. As line length and contrast are reduced, there is a shallowing of the FOS slope with elevation of IQR (Fig. 2) which is of an appearance similar to that seen in glaucoma,25 where retinal ganglion cell (RGC) function is impaired. The increased variability seen with reduction in sensitivity is similar to that found by Henson when assessing SAP in normal and diseased eyes, and which he attributed to a difference in functional ganglion cell density.29 An alternative explanation is that the precision of the central detector is directly related to its spatial tuning, with greater precision associated with higher spatial frequency tuning.
There remains a widespread assumption that tests of motion are directed to the M pathway, although this may be an oversimplification. Clear segregation of the retinogeniculate pathways exists up to the lateral geniculate nucleus (LGN), but within the cortex multiple interconnections occur, making the pathways less independently driven than was once thought.30 31
The conventional view is that the M pathway acts as a motion-detection system which shows greater sensitivity, relative to the P pathway, to stimuli of low-contrast, high temporal frequency, and low spatial frequency.32 33 34 The MDT conditions of this study are at a moderate temporal frequency (5 Hz), for stimuli of high spatial frequency, over a wide range of contrast (25%85%). It is therefore possible that responses from both M and P contributions are invoked. However, this view is contradicted by a recent study, which suggests that the M pathway is primarily responsible for vernier acuity thresholds, with higher neural mechanisms being responsible for spatiotemporal integration.35 Oscillatory displacement thresholds, measured in the presence of stationary references, have been shown to saturate at contrasts below 15%, consistent with the M pathway contrast response function.36
The energy plot (Fig. 6) shows a constant stimulusresponse relationship over the range of contrasts assessed, which could apply to either a single or multiple pathway contribution. There is no evidence of a transitional shift in the threshold or slope, which would be a possible finding if different pathways were to predominate under different stimulus conditions within the experimental design. All stimuli were presented at contrasts equal to or greater than 25%, which is above the level at which the M pathway responses are generally considered to saturate. However, the results of this study suggest that it is not the contrast per se that modulates the response, but the stimulus energy (Fig. 6) so that saturation for higher contrast stimuli occurs only if the stimulus is also large. Saturation is generally measured with grating stimuli and so the "total energy" within a "receptive field" may be much larger than with the single-line displacement stimuli used in this experiment. It is therefore not possible to hypothesize mediation of the response by a particular cell type for these experimental conditions.
Table 3 shows the estimated number of RGCs stimulated by lines of various lengths. Values were approximated using estimates of human ganglion cell density (GCD)37 and coverage values38 (see the Appendix). The data in the table may underestimate the actual number of receptive fields stimulated, in that no allowance is given for the optical effect of line spread.39
|
The anatomic extent of the dendritic tree corresponds to the electrophysiological central receptive field, but may be smaller than the psychophysical receptive field estimate.43 In humans, dendritic tree measurements approximate to 60-min arc diameter for M cells and 20-min arc diameter for P cells, at 30° eccentricity in the nasal field (Fig. 8) .38 Human dendritic trees are larger than their counterparts in the macaque, predicting larger receptive fields.38 The estimated psychophysical receptive field size in the macaque, at 30° eccentricity, for drifting sinusoidal gratings approximates 42 min arc.32 44 45
|
The threshold in SAP is determined by the relationship of the area (A) and intensity (I) of the stimulus, given by IAn = constant. The summation coefficient (n) varies for different conditions: Where the stimulus lies within Riccos critical area, n = 1. (AI = constant), and Riccos law is said to hold. For larger stimuli, summation is incomplete, and the threshold is determined by probability summation, which is achieved by recruitment of adjacent receptive areas.46 Pipers law holds when n = 0.5.22
The results of this study suggest that the same stimulusthreshold relationship may hold for the MDT stimulus. Longer lines extend over many ganglion cells, yet result in the same threshold as shorter, brighter lines of equivalent energy (area x luminance; Fig. 6 ). This is consistent with linear summation of RGC responses for the line lengths used in these experiments. Riccos law holds for the stimulus parameters tested in this study (Fig. 6) , so that there is a constant MDT threshold where the product of line length and luminance is constant.
Riccos area is psychophysically defined22 and estimates of the size of Riccos critical area vary according to the type of the stimulus,47 indicating different patterns of spatial filtering.44 Riccos area changes with eccentricity and this is thought to reflect changes in the underlying ganglion cell type, density, and neural convergence.47 48 In humans, Riccos area measures larger than the estimated anatomic dendritic field size, suggesting summation across RGCs and possible higher neural input.49
By contrast with DLS, electrophysiological studies of spot light displacement in primates indicate that all RGCs respond to movement of <10% the diameter of their receptive field. It is thought that intrareceptive field directional sensitivity is achieved by pooling of the excitatory and inhibitory responses of overlapping on-off receptive fields.27 50 The displacement threshold is thus determined by the sensitivity gradient between the two stimulus positions.12 51 It is possible that, as the line stimulus oscillates from one position to another, different parts of the concentric on/off ganglion cell receptive fields are stimulated as an alternating decrementalincremental response.
A route to understanding the relationship between stimulus energy and the MDT threshold lies in ideal detector modeling, whereby a hypothetical model of optimal performance is created for known variables. Geisler52 predicted a slope of 0.5 between log threshold and log stimulus intensity for a hyperacuity stimulus calculated on the basis of preneural factors (line-spread function, pupil size, retinal receptor anatomy, and Poisson distribution of photon absorption) and an ideal detector. This conclusion is consistent with the findings of the present study. The optimal receptive field was built by simple linear summation and inhibition. In a later paper, Geisler and Diehl53 advise consideration of both preneural and hierarchal neural factors, with receptive fields regarded as free parameters.
The quantal absorption of photons by the retinal detectors is a source of noise described by the Poisson distribution. Noise is a constant feature in all sensory pathways, and signal detection is dependent on the average number of quantal absorptions necessary to elicit one extra impulse. This level determines gain control, and the signal-to-noise ratio grows in proportion to the square root of the stimulus area.54 55
Previous studies have suggested that vernier acuity is predicted from changes in the relative luminance of the retinal image.56 57 58 The results of the present study indicate that MDT sensitivity is achieved by detection of the total energy difference as the stimulus is displaced with respect to its environment. The MDT stimulus is displaced in position over time and is therefore probably best modeled using a three-dimensional space-space-time filter input.59 Adelson and Bergen60 suggested a two-stage hierarchal model of local motion signals, in which the first stage represents pairs of linear filters that are oriented in space and time and tuned in spatial frequency. The response of each pair is squared and summed. The second stage computes "motion energy" on the assumption of neural opponence. The MDT is an oscillating "A-B-A" square wave stimulus of apparent motion. The displacement threshold may therefore represent the summed difference of the on-off responses and also opposing directional filters, as the stimulus is displaced from position "A" to second position "B" and back to starting position "A."
The novel finding of this study is that the MDT threshold is related to the square root of the stimulus energy for the stimulus conditions described. Thus, a new principle, TED (threshold energy displacement), is proposed to apply to MDT summation properties, giving the relationship: T = K
E (where T is the MDT threshold, K is the constant, and E is the stimulus energy: stimulus area x stimulus luminance background luminance). TED offers a numerical basis for modeling the structuralfunctional relationship of MDT thresholds and will be applied to the new multilocation MDT under development. It is anticipated that there is the potential to quantify neural losses from glaucoma with MDT using the principles described for SAP.44 61 62 Future work will explore the relationship of MDT thresholds with stimulus energy at other retinal locations.
| Appendix 1 |
|---|
|
|
|---|
Eccentricity =
(x2) + (y2).
Therefore, eccentricity at x = 27 and y = 3 is
(27 x 27) + (3 x 3) = 27.16.
Magnification = 0.286 [0.000014(ecc2)].
Therefore, magnification at 27.16 is 0.286 [0.000014(27.16 x 27.16)] = 0.276.
Number of ganglion cells stimulated = line length (in millimeters) x linear estimate of ganglion cell numbers.
Estimated ganglion cell density at 27.16° = 877 ganglion cells per mm2.37 62
Linear estimate =
GCD = 30 ganglion cells per millimeter.
Example of Calculation of the Number of RGCs Stimulated at 27,3 by a 128-min arc Length Line
Line length is 128 min arc = 2.13° = (0.276 x 2.13) mm = 0.588 mm.
Line width is 5 min arc = 0.08° = 0.02208 mm.
Therefore stimulus length stimulates 0.588 x 30 = 17.6
18 ganglion cells (all types).
Width stimulates: 0.022 x 30 = 0.66 ganglion cells (all types).
Calculation of the Number of P and M Cells and Coverage Values for each "On" and "Off" Mosaic
The number of P (midget) and M (parasol) cells was estimated from percentage estimates of each cell type from primate histologic data. At 30° eccentricity, this approximated to 20% parasol and 45% midget.40
Coverage is defined as the average number of cells with receptive-field centers that include the same point in visual space. It is identical with overlap (the average number of cells with receptive fields that cover that of a given cell).63 The anatomic estimate of coverage is given by cell density (in square millimeters) x dendritic field area (in square millimeters).40 The calculations are adjusted for coverage of each "on" and "off" mosaic.38 42 The estimated coverage is 3:1 for the parasol cells and 1:1 for the midget cells.40 Calculations are further adjusted for the higher cell density of "off" P cells, located in the outer inner plexiform layer (IPL), compared with "on" P cells, which are located in the inner IPL (outer: inner density ratio of 1.7:1).42 M cells are assumed to be divided equally between "off" and "on" mosaics.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication November 21, 2005; revised April 18 and May 22, 2006; accepted September 5, 2006.
Disclosure: G.M. Verdon-Roe, Pfizer (F); M.C. Westcott, Pfizer (F); A.C.Viswanathan, Pfizer (F); F. W. Fitzke, Pfizer (F); D.F. Garway-Heath, Carl Zeiss Meditec and Pfizer (F)
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: Gay Mary Verdon-Roe, Glaucoma Research Unit, Moorfields Eye Hospital, London, UK; g.m.verdon-roe{at}ucl.ac.uk.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |