Several attempts have been made to model the spectral reflection of the fundus of the eye by assuming reflection at different retinal layers interspaced with ocular absorbers. In 1986, Van Norren and Tiemeijer1 measured the reflectance of the fovea, a peripheral site, and the optic disc at 14 wavelengths. Measurements were rather cumbersome, as each wavelength needed a separate session. Their instrument optimized the directional reflection from the receptors (the optical Stiles-Crawford effect, or OSCE) by using entry and exit pupils that were small and close together. Their model had two reflecting layers, one at the retinal pigment epithelium (RPE) and one at the sclera. In addition, four absorbers were distinguished; lens, macular pigment, melanin, and blood. In that same year, Van Blokland and Van Norren2 presented reflectance measurements with an instrument that enabled detailed information about the shape of the directional reflection, but only four laser wavelengths were available. In a descriptive model, the authors suggested several trajectories that light takes in the receptor layer.
A few years later, Delori and Pflibsen3 build a reflectometer based on a grating spectrometer, allowing all wavelengths to be measured simultaneously. The optical design used the exit pupil of a Zeiss fundus camera and thus was not optimized for obtaining the directional reflection from the receptors. Their model added a graded reflection from the choroidal space to the model of Van Norren and Tiemeijer.1 They achieved a better estimate of absorption in melanin and blood, but the receptor layer was ignored. Van Norren and Van de Kraats presented an instrument that could measure 15 wavelengths quasi-simultaneously.4 The instrument was again optimized for the directional reflection, and it used low light levels to enable retinal reflection spectra with the nearly full complement of visual pigments. A new model was derived by Van de Kraats with better funded trajectories of light in the receptor layer.5 Spectral reflection was measured at two pupillary positions, one at the maximum of the directional reflection and one at a position with the entry and exit pupils shifted horizontally by . A reflection from the stack of cone receptor discs was incorporated in the model, based on measurements with high and low levels of visual pigment and with high and low directional reflection. The model permitted the isolation of pre-receptor, receptor, and post-receptor signals, thereby improving the accuracy of estimating the pre-receptor parameters and lens and macular pigment density.
A new instrument proposed by Van de Kraats and elaborated by Zagers 6 combined the directional resolution earlier presented by Van Blokland and Van Norren,2 with the spectral resolution of Delori and Pflibsen.3 In 1 second, a two-dimensional (2-D) image could be measured from a foveal spot illuminated with intense white light, containing in one dimension the spectral aspects and in the other the directional aspects. The fast measurement allows large samples. In their data analysis, Zagers, 6 as a first approach, used the directional component of the signal to estimate the pre-receptor absorbers and parameters for the optical quality of the receptor layer.
The purpose of the present study was to propose and test a model that fully explains both the spectral and the directional aspects of the reflection of the human fovea, i.e., the directional cone reflection and the nondirectional background reflection. The model generates quantitative information about reflectors and absorbers. In addition, the dependence on age of these parameters was studied. A validation of the estimates of macular pigment density was obtained with three other methods in an identical subgroup of 53 subjects.
Optical Stiles-Crawford Effect in the Pupil Plane
The angular reflection, , from an illuminated spot at the fovea shows up in the pupil plane as a Gaussian-shaped spatial intensity distribution (the directional component with amplitude ) on a pedestal of background light (the nondirectional component ):2, 7, 8the directionality, and the horizontal location in mm in the pupil plane relative to the center of the pupil. The peak of the OSCE is generally not exactly at the center of the pupil; the difference is . Measured curves are cut off by the edges of the (generally dilated) pupil (Fig. 1 ).
Spectral Behavior of Directionalitya waveguide component inherent to the receptors that is constant over wavelength, a scalar preceding a scatter component that varies with wavelength due to interference effects between receptors, and the wavelength in nm (Fig. 2 ; the data are described in Sec. 5). The equation stems from Marcos, 9 but it was slightly modified by the normalization at to enable easy comparison of the contributions of and .
Origin of the Directional Component
The origin of the directional reflection has been located on the discs of the cone receptors5 or at the end of the outer segments.12 Experiments with the visual pigments in a bleached condition as used throughout this paper cannot discriminate between these two options. In line with our earlier interpretation5 based on measurements in both the dark and the bleached states of the visual pigments, we opted for choosing the discs as the source of the directional reflection. The reflection from the layers posterior to the receptor layer is considered nondirectional. The principle is illustrated in Fig. 3 , where ray tracing through a lens and through a more complex refractive element, a receptor, results in the same reflection value as from the bare surface. In the retina, finally all rays through the receptor layer, guided or nonguided, hit the layers posterior to it, including the highly scattering choroid. Reciprocity theory states that the reflected rays take the same paths on their way back, preserving the wide-angle scattering nature, at least over the full width of the pupil.
To prove this, Van de Kraats 5 described an experiment with an annulus around the measuring area to illuminate the receptor layer only from the choroidal side. As expected, the reflected light from the central, nonilluminated part showed heavy absorption because of blood and melanin, but no directionality. Similarly, in experiments from the Burns group using the fluorescence of lipofuscin in the pigment epithelium as a new source of light located behind the receptors, no directional light was found. 12, 13, 14, 15 The implication of both types of experiment is not that light avoids the pathway through the photoreceptors on its way back,12 but rather that the Gaussian-shaped directional light guided through (recaptured by) the photoreceptor outer segments is perfectly complemented with inverted Gaussian-shaped light through the pathways passing the outer segment wall. These pathways are in reversal to that of the (nonperpendicular) light escaping the outer segments on entering the cones from the direction of the pupil (Fig. 3). In conclusion, any directional light must originate from reflections in the cones themselves (including the outer segment tips), and not from light that also traveled the layers behind it.
Cone Capture Area
From (optical) antenna theory,16 it follows that the capture area of a cone is not purely determined by its physical dimensions, but can be larger.5, 17 Estimations using waveguide theory also give a significantly higher effective capture area compared to the physical frontal area of a cone.18 We therefore assumed that all perpendicular light enters the cones and that none reaches the interspaces. For oblique angles, light escapes from the cone outer segment into the interspaces, according to the generally accepted explanation of the Stiles-Crawford effect.19 Energy distributions calculated in rods showed almost identical behavior for wavelengths of 475, 505, and .20 On this basis, we also assumed that, again for the perpendicular case, the light captured by the foveal cones is constant over wavelength. In the bleached condition, the light passes the cone outer segments into the deeper layers. Due to the reciprocity principle, a large part of the light reflected from the deeper layers must be recaptured by the cones.
An earlier version of the model with an extensive discussion was published by Van de Kraats 5 The main deviations in the present model concern the regaining of scattered light of the nondirectional component, a wavelength-dependent reflection from the cones, the use of a tapered blood layer thickness, and new templates for the eye media and the macular pigment. A schematic view of the model is given in Fig. 4 . Briefly, with the directional reflection originating from the cone receptors, only absorbers anterior to the cones, i.e., lens and macular pigment, leave their spectral fingerprint.5, 21 The nondirectional amplitude is in addition influenced by absorption in blood and melanin.3 In the anterior eye, light is absorbed by the media. The first (nondirectional) reflection occurs in front of the receptor layer, at the internal limiting membrane (ILM). Although the reflection at the ILM is specular and therefore directional, it does illuminate the pupil quasi-uniformly and can therefore be treated as nondirectional because it emanates from an image of the source located very close to the retina.22 Next, the only directional reflection takes place at the cone outer segments. A further nondirectional reflection occurs at the retinal pigment epithelium. Last, light is reflected at the choroid, with absorption in blood and melanin. To avoid complex formulas, losses at reflecting layers of a few percent were ignored in the calculation of posterior layers. Visual pigments were supposed to be fully bleached.
Absorbers in Both Directional and Nondirectional Components
The transmission of the various absorbing layers was described by their optical density (Fig. 5 ). In general,23 is a spectral template for losses in Rayleigh scatter in the cornea, humors, and lens; is a template for losses in the young lens; and is for losses in the aging lens. The parameters , , and , derived from fitting the measured data to the model, were the age-dependent densities at . A term for a neutral density of the media to provide for some losses due to scatter by relatively large scattering particles was ignored. is the density of of water; we used24 a fixed total thickness of water of .
For the case of nondirectional reflected light, we defined , almost similar to the preceding case, but with set to zero. This assumption was made after fitting the data from the isolated nondirectional light in the pilot analysis, where showed a strong tendency toward zero. The rationale is that the light that is scattered somewhere anterior to the cones changes angle of direction. It is therefore easily lost for the directional reflection, while it is not easily lost for the nondirectional component consisting of light scattered perhaps a second time in deeper layers. This assumption of regaining scattered light in the nondirectional condition was tested with nearly noiseless mean data from a young and old group (Sec. 5.3).
The transmission of the macular pigment was described by25: to 1. Thus, parameter is the density of macular pigment at .
The OSCE was described by a Gaussian shape in the pupil plane:the directionality, and the horizontal position relative to the entrance beam in the pupil plane of the instrument. Due to the alignment procedure (see Sec. 4.3), the position of the maximum of the OSCE at was always close to zero. Note first that we did not use a reference position relative to the center of the pupil as in the general form in Eq. 1. Note, too, that Eq. 1 also includes the nondirectional component.
, the directionally reflected light originating from the cone discs, was assumed to decrease with wavelength asthe reflection from the outer segments. The wavelength dependency deviates from an earlier description5 for reasons addressed in Sec. 6.5.
At the level of the receptor layer, the directional reflection becomes
The nondirectional component is easiest considered from the choroidal space. The backscattering of light in the choroidal tissues was taken as a neutral reflection . Light in the choroidal space is absorbed by blood and melanin.
We defined the density of a layer of 95% oxygenated blood with a thickness of as .26 To account for the variety of path lengths through the center and periphery of small and large blood vessels, we assumed for the template a wedge-shaped blood layer, with path lengths from zero to . These fixed path lengths were taken to avoid complex path length corrections during the fitting process. To obtain the final blood thickness parameter as a scalar in , we divided the wedge results by 100. The template for a layer of blood of with correction for different path lengths can now be calculated asin is
For the transmission through a uniform layer of melanin, we took the spectral data of Gabel, 27 approximated by
The term 2.45 served to normalize the function to 1 at so that is the parameter for the density of melanin at . The density of melanin becomesthe parameter for a neutral reflection at the RPE.
The complete description of the nondirectional reflection isrepresents a parameter for a neutral reflection at the ILM, and is the special case for the density of the media without Rayleigh scatter losses, as explained in Sec. 3.2.
Complete Model1 ). In a single measurement, there are about 4500 free parameters (90 at 5-nm intervals in the spectral range from 400 to at 0.1-mm intervals in the pupil range; see Sec. 4.1).
Mean parameters for the complete group of subjects.
|Parameter (reference wavelength)||Value (standard deviation)||Correlation with age||p -value for age effect||Repro (%)2|
The instrument described by Zagers 6 was redesigned in a desktop version that could be aligned with a joystick. In addition, the spectral range was extended to 400 to , and the switching mirrors were replaced by beamsplitters, enabling continued observation of both pupil and retina during measurements. A halogen lamp L (12-V, 30-W Wotan 64260, Osram, Munich, Germany) illuminated a spot of 1.8-deg diameter on the fovea (Fig. 6 ). The light was spectrally filtered by F ( BG26 Schott AG, Mainz, Germany; 1-mm Schott UG3; Unaxis TL60, Linos, Goettingen, Germany) for the comfort of the subject and to prevent overloading the sensitive CCD camera (KX85 Apogee Instruments, Inc., Auburn, Massachusetts) that served as the detector. The intensity of the spot was Td; calculations showed that it could be viewed safely for .28 The filament of the halogen lamp was imaged in the pupil plane of the eye, where it defined the illumination pupil. With a separation of below this illumination pupil, a slit-shaped detection pupil S of the instrument of formed the input for a prism-based imaging spectrometer. The 15-mm slit length allows ample room for aligning an eye with a dilated pupil and an eccentric OSCE maximum (see Sec. 4.3). The 2-D image on the CCD detector had in one dimension the directional information from the intensity distribution over the slit-shaped exit pupil. In the other dimension, it contained the spectral information, as light from each point of the slit was decomposed by the prism. Only light from the central 1.5-deg foveal spot was used for detection. Because of the prism, the image at the CCD with pixels had an alinear spectral axis. It was transformed by software, using the dispersion calculation, to an intensity-normalized and linear reflection image of 250 pixels at 0.1-mm intervals at the pupil axis, and 90 pixels at 5-nm intervals at the spectral axis (400 to ). Refraction errors were compensated by adjusting a Badal-type front lens system (L1 and L2).
A mercury lamp (90-W, Type 93136, Philips, Eindhoven, The Netherlands) was used for spectral calibration. Before each measurement session, the instrument was calibrated for sensitivity by mounting two different calibration tubes on the front lens. An instrument stray-light image was assessed with a dark reference tube, containing a light trap. All subsequent images were corrected with this image. A white reference tube contained at the end a surface painted with white paint (Kodak 6080 White Reflectance Coating, Eastman Kodak Company, Rochester, New York). Because the distance of the white surface was at from the instrument’s pupil plane (10 times the focal distance of the eye), the measured white reference image represented a 1% reflection from the eye.
A group of 102 subjects, mean age 49 years (range 18 to 75), took part in the measurements. The data of 53 subjects were published earlier using another analysis technique.29 The tenets of the Declaration of Helsinki were followed, and the local Medical Ethics Committee approved the protocol. Before the experiment started, the nature of the experiment was explained to the subjects, and written informed consent was obtained. The pupil of one eye was dilated with one or two drops of tropicamide 0.5%. Subjects were aligned to the instrument, with the entrance beam clearly visible in the anterior eye. Initially, the entrance beam was placed somewhere above the center of the pupil, to allow space for the invisible detecting slit below the entrance beam. A chin rest and temple pads helped maintain the position of the head. The subjects were instructed to fixate the center of the illuminated spot. The front lens was adjusted for a sharp image of the retina at the monitor. Next, the subject leaned back, and the dark reference image was taken. Measurements at a rate of about 2 per second were displayed at the computer screen in the form of two cross sections of the 2-D image, the spectral shape at pupil profile position zero, and the pupil profile at . The alignment was optimized by searching in the pupil plane for the peak of the OSCE. This process was facilitated by simultaneously displaying (in a different color) the highest profile from the start of the measurements. The optimal horizontal position in the pupil plane was easily found by moving the instrument until the peak position of the Gaussian-shaped directional profile fell symmetrically around profile position zero. Profile position zero corresponds to the horizontal center of the entry beam. The optimal vertical position in the pupil plane was found by trial and error using the highest profile trace. During optimizing, the spectral view was watched to keep unwanted reflections from the anterior eye low. These reflections were recognized by their relatively high reflectances near . The optimization procedure took about , enough to bleach the visual pigments to almost 100%. At the optimal position, five measurements with an integration time of each were obtained. If necessary, head position was readjusted when the head drifted away from the optimal position, recognized by a drop in amplitude.
The reflection profile, as seen in the entrance slit of the instrument, is cut off by the pupil of the eye. Left and right limits in these pupil plane positions were determined as follows. First, the data in the spectral range from 500 to (where the signal has low noise) were binned to , yielding 20 profiles. Next, left and right from the central pupil position, those data points where the signal dropped to 20% of the central value were taken. From those 20 left and right limits, we took the lowest. In the final step, the left and right limits were moved to the center and applied to the complete measurement to avoid any interference from the pupil edge.
Subselection of Young and Old Groups
To develop and evaluate the model, two subselections were made from the group of subjects: a group years (mean 24), and a group years (mean 58). From the five measurements per subject (see Sec. 4.3), a further selection was made by allowing measurements only over at least , with pupil limits outside the range , (nasal to temporal), and a minimum of three measurements per subject. This gave 45 measurements from 10 subjects for the young group and 81 measurements from 18 subjects for the old group. To find a template for as a function of wavelength, the pupil profiles at every wavelength (binned to ) for the young group were analyzed for with the Gaussian model from Eq. 1. The means of the standard deviations within a subject were calculated to provide an indication for the error in the data points. Next, a “young measurement” was calculated by averaging the data points in the young group within the pupil range they had in common in their original measurements; similarly, an “old measurement” was calculated.
Fitting the Data
For the pilot explorations with the young and old group data, the Solver in Microsoft Excel 2003 was used to fit the data with the model by minimizing chi-squares. All spectral templates were convolved with the bandwidth of the spectrometer. For the final analysis of individual measurements with the 2-D model, the Levenberg-Marquardt routine30 was applied. The noise in the data points, used for weighting, was calculated from the square root of the photons counted. For analyzing the statistics, SPSS 12.0.1 for Windows (SPSS, Inc., Chicago) was used.
To illustrate the spatial profile and the spectral behavior of the directional and the nondirectional components, the data from a single measurement of a 20-year-old subject are presented in Fig. 7 . The amplitudes of the directional component and the nondirectional component were estimated by fitting the elevated Gaussian curve from Eq. 1, at each wavelength to the pupil profile. The result at is presented in Fig. 7A. The amplitudes versus wavelength are presented in Fig. 7B. The total reflection at the peak of the profile was found by adding the amplitudes of the directional component and the nondirectional component. Going from left to right, both curves show a sharp increase beyond , where the lens absorption ends (young eye), and a further increase near , where the absorption from the macular pigment (MP) ends. The effect of macular pigment is stronger in the directional curve than in the nondirectional one. The nondirectional curve displays an increase near , where blood no longer absorbs, together with a gradual increase due to increasingly less absorption in melanin. The directional curve lacks these features. Beyond , the curves show the influence of absorption by water, leading to a strong decline at .
The more complex shape of the nondirectional component due to the addition of light from the choroid, visible above , is evident. The dilution of the macular pigment fingerprint near , due to preretinal reflectors, is also clearly visible.
When the directional component becomes low in the absolute sense, or low relative to the nondirectional component, the estimation of becomes increasingly noisy, visible below 500 and above . Through interaction in parameters, noise in then also causes noise in the separation of the directional and nondirectional amplitude. To resolve this, a smooth template for as function of wavelength was derived (Sec. 5.2). The single parameter for , instead of a value at each wavelength, also keeps the total number of parameters in the model limited to 13 instead of more than 100.
Template. Below , the result was very noisy due to the high density of the aged lens. Also, at the longer wavelengths above with the increasing amplitude of the nondirectional component, the more irregular shape of the reflection in the pupil plane in the older group appeared as noise in the directional component. The error bars in Fig. 2 are the means of the standard deviations within a subject. These values were used for weighing the data in the fitting process. Equation 18 normalized to 1 at was used as a template for as a function of wavelength.
Test of on Mean Directional and Nondirectional Spectra for the Young and Old Groups
The almost noiseless mean data from the young group and the mean data from the old group were used to test the general concept of the spectral model from Sec. 3, and in particular the assumption of for the nondirectional part. With the template for the spectrum of [Eq. 18], data were fit using Eq. 1 to derive the directional and nondirectional amplitudes. Results of these analyses are shown as the data points in Fig. 8 . Model curves were drawn using the equations from Secs. 3.3, 3.4. The directional data were fitted first to Eq. 10 (with and fixed to the calculated values according to the mean age of the groups), yielding values for , , and . The nondirectional data were fitted next to Eq. 16 (with , , and fixed to the values obtained in the directional fit, and fixed as above and set to zero), yielding values for , , and , , and .
An almost perfect fit to the data could be obtained (chi-square ; chi-square ). When was not set to zero for the fitting of the nondirectional spectra, but also taken as a fixed parameter according to the mean age of the groups [by replacing in Eq. 16 by ], chi-square values were more than an order of magnitude higher (chi-square ; chi-square ). The two-step strategy of first fitting the media and macular pigment parameters on the directional component, and so avoiding interaction with deeper layer parameters, was applied in this section only to strengthen the assumption of setting to zero for the nondirectional component. In the rest of this paper, fitting to the data of individual measurements, all free parameters were fitted simultaneously.
Next, the model was applied to the individual measurements, allowing the study of age effects in the parameters. In total, 510 measurements from 102 subjects were analyzed. At first, all 13 parameters were allowed to vary, except water (fixed to ). As expected, the media term showed a high correlation with age (0.843, ). The media parameter , however, showed an unexpected negative trend with age.23 Inspection of the original data showed that with the increased density of the lens at higher ages together with the attenuation by , the measured reflection near reached the noise floor of the instrument. This limited the necessary dynamic range for proper detection of . Other parameters have a much wider effective spectral range of their templates and were less sensitive to this problem. In the receptor layer, only showed a significant correlation with age ( , ). In the deeper layers, reached the value . Last, a very low correlation was found for (0.209, ).
The media scatter term showed a low (but not significant, ) correlation with age, as found in an earlier study.23 To avoid problems due to the low signal at , both media parameters and were set to the calculated values for a mean age of 40,23 leaving as a free parameter for the media. This left 10 free parameters for the final fitting of the data. The mean results are summarized in Table 1. Parameters with a significant age trend were plotted versus age squared in Fig. 9 . Age squared was used because it consistently yielded smaller chi-squares than a linear fit.23 Reproducibility (the repro column in Table 1) was defined as the average coefficient of variation ( standard deviation in the five measurements per subject/mean parameter value).
Parameter Cross Influences
The mean covariances between parameters of all individual measurements were calculated from covariances provided by the fit routine. Only three appeared higher than 0.6: the value 0.694 between and the reflection from the RPE, 0.681 between the ILM reflection and the density of macular pigment, and 0.620 between melanin and the reflection from the choroid.
The Pearson bivariate correlations between parameters after correction for age were also calculated. Neglecting correlations below 0.266, corresponding to 7.1% of the variance explained, only two remained with . The first is between and (0.545), the second between melanin and blood .
We developed a model for fundus reflection that for the first time simultaneously describes the spectra of the directional and nondirectional reflection of the human fovea at a high spectral and spatial resolution. The excellent fits support the assumption that the deeper layers generate no directional component. Van de Kraats 5 already gave experimental and theoretical evidence, and this was also confirmed by Burns 14 using the autofluorescence of lipofuscin. As the fluorescence is by nature omnidirectional, it formed a new source of light behind the receptor layer, undisturbed by preretinal reflections. We found that a large part of the light from deeper layers must be recaptured by the cones due to the reciprocity principle (Sec. 2.4). At first glance, a surprising aspect of the model is that Rayleigh scatter losses in the eye media and in preretinal layers have to be included only in the directional pathway. This will be discussed in the next section.
In the initial analysis of the complete group of 102 subjects, the lens young component and the Rayleigh scatter component were free parameters. The mean value of was , and was . In a pilot analysis, the reflection near appeared so low due to the combined attenuation caused by and that it reaches the noise level of the instrument, limiting the dynamic range for proper detection. In a second analysis, therefore, and were fixed to mean values corresponding to an age of 40 years (1.26 and 0.5).23 This prevented less meaningful results, especially at older ages when increases substantially. The decrease with age of the lens young component as shown by Zagers probably was contaminated by this effect.21 Going from 20 to 60 years, the density for the total media at in the Zagers study increased by units. Our data showed an increase of at . This suggests that the Zagers increase in age of the lens old component was compensated by a too large decrease of the lens young component.
The significant improvement of the fits by using the value for of 0.5 in the directional reflection, but setting it to zero in the nondirectional one, was surprising at first. At second thought, it was explained by part of the light that is scattered anterior but near to the cones. Light scattered in that layer is changing the angle of direction, and therefore is easily lost for the directional reflection, while it is not easily lost in the case of the diffusely scattered nondirectional component. Van den Berg 31 presented optical scatter losses in the cornea, but in fact they included the humors and preretinal layers as well. Their magnitude of the Rayleigh scatter losses for a 1-deg field can be calculated as 0.34. Estimates of losses of light scattered outside a retinal field of 1-deg by the cornea and lens for the relatively young eye from Vos 32 result in a density of only 0.03. For our larger illumination field of , it should be even less. According to this, a large part of our Rayleigh scatter value of 0.5 originates from additional scatter losses in preretinal layers anterior to the cone photoreceptors. The model is however unable to discriminate between these losses and scatter losses in the eye media.
The aging slope of for the media attributed here solely to ( was fixed) is comparable to the sum of the aging values and from a literature study.23 Apparently because of fixing the value of to 0.5, some of its (small) aging aspect has moved into .
The macular pigment optical density in a subgroup of 53 subjects was analyzed before with two previous models and compared with densities obtained with heterochromatic flickerphotometry (HFP), scanning laser ophthalmoscope (SLO), reflectance difference at 488 and , and SLO autofluorescence.29 Analysis of the data with the 1996 version of the fundus reflection model5 resulted in 0.44, compared with the current 0.37. A simple model containing only the directional component21 yielded 0.53, being too high because of the earlier mentioned parameter interaction with the Zagers lens templates.21
Apart from the current wavelength-dependent reflection from the cones, some of the differences can be due to the different absorption data for macular pigment used. The Handelman data25 we used show a somewhat different density near (see Sec. 3.2) than previously used data. Comparing the absolute numbers with results from other instruments (Table 2 ) is also more difficult because of the strong effect of the retinal area under test (here, ). Unexpected at first sight, the models, while using a larger field, yield the highest densities. This is explained by the fact that the models compensate for all disturbing factors, like the diluting effect from light reflected at the ILM. The necessity for a good model is seen in Fig. 7B, where the macular pigment fingerprint is unmistakably higher in the directional component compared to that of the nondirectional component. The latter clearly shows the diluting influence of the ILM reflection. The SLO reflection method was lacking such compensation. More difficult to understand is the lower value for the SLO autofluorescence technique, because here dilution by preretinal autofluorescence is not expected.
Macular pigment density using different methods.
|Method||Field size (deg)||dMP||Correlation with current|
We repeated the calculation of the Pearson bivariate correlations on macular pigment between the current model results, previous models, and the other techniques (Table 2). All correlations were highly significant . The parameter showed no significant change with age, as in most other studies,29, 33, 34 with the notable exception of a recent study including over 800 subjects reported by Nolan using HFP that showed an unexpected decrease of 0.17 from 20 to 60 years.35
For the complete group of 102 subjects, despite the covariance of macular pigment with the reflection from the ILM of 0.681, the model seemed to succeed well in separating them, as shown by the low and not significant correlation after correction for age (0.113, ).
For the absorption in blood, we used a template derived from a wedge-shaped layer of blood. This is a more realistic model of the blood layer, because it contains a large range of path lengths through the blood layer instead of always the same length as with a homogeneous layer thickness. The result is that the large dynamic range in the transmissions over the whole spectrum in the case of a homogeneous layer thickness is compressed, yet without adding more parameters to the model.3, 36 declines significantly with age, probably due to the accumulation of basal laminar deposits on the choriocapillaries and Bruchs membrane,37, 38, 39 screening the choriocapillaris.
Measurements of choroidal blood volume with a laser Doppler flow technique also shows a decrease of about 45% from 20 to 70 years.40 Another study using morphometric analysis showed a 45% decrease in the density and a 34% decrease in the diameter of the lumen of the choriocapillaris from the first to the tenth decade.41 It is difficult to say how these phenomena impact on , and what fraction of is from the choriocapillaris.
The correlation between blood and melanin of is likely to have a true physiological background. Their templates are so different that the interaction of parameters in the model is not expected. This is expressed by the low covariance between these parameters (0.138). One obvious explanation is that blood and melanin in some degree crowd each other out.
Melanin was not found to be age dependent. This is a minor contradiction with determinations from young and old donor material based on only the retinal pigment epithelium, decreasing about 37% in the macular area.42 In analyzing the results of the macular area of only 16 subjects in another study, an age relation could not be found.43 Measured in optical density units, the fraction of the melanin in the RPE is about one-third of the total melanin.27, 44 If two-thirds of the melanin fraction in the choroidal space is constant with age, age dependency is expected to be much lower. While pure melanin decreases with age, more complex granules of the type melanolipofuscin and melanolysosomes increase by 72%.42 With their absorption spectrum not very different from pure melanin,45 discrimination with our reflection technique is not possible—another argument for an almost zero age dependency of melanin. More important is that in the case of using reflection techniques to derive the total melanin content, it better reflects the functionality of melanin for the stray-light reduction in the retina. Other techniques based on pure melanin only are perhaps less important in this respect.
Despite the covariance of melanin with the reflection from the choroid of 0.62, the model seemed to succeed well in separating them, as shown by the low correlation after correction for age (0.252, ).
The first, small reflection (approximately 0.1%) occurs at the ILM. A reflection at the outer limiting membrane was not incorporated because it cannot be discriminated from the reflection (0.56%) located at the receptor layer/retinal pigment layer interface. Although a reflection at Bruchs membrane seems physiologically probable, it would add at least two parameters ( and ) to the deeper layers, as in Delori 3 Without these extra parameters, we obtained very good fits. The price to be paid is that the absolute values of and have limited physiological meaning; trends are relevant, however. Reflection of the choroid increases somewhat with age (Fig. 9), possibly due to the same deposits as discussed for the increases of blood layer thickness with age.
In previous papers, we assumed that the origin of the cone directional reflection lay in spectrally neutral Fresnel reflections by small differences in the refractive index of the discs and the interstitional fluids.5, 21, 46 Although this was supported by measurements of Zagers on pseudophakic subjects,21 Choi 47 argued that the wavelength range used was too limited to support the conclusion. Other authors48, 49 assumed the integral volume of the Gaussian shape to show spectrally neutral behavior, as is commonly found in diffraction-limited processes.16 Consequently, the peak amplitude should drop with . With the current, much-extended wavelength range, the peak amplitude was clearly seen to decrease at the longer wavelengths, where other absorbers play no role [Figs. 7B, 8]. We therefore also opted for behavior in the new model. Following Marcos 9, 10 the of the waveguiding part of the OSCE is constant with wavelength. Another scatter component from interference of light from neighboring cones introduces behavior, as seen in Fig. 2. Summarizing, we have shown that the light intake by the cones on perpendicular entrance is neutral because of the large capture area (Sec. 2.4). The summated reflection from the discs itself still remains a neutral Fresnel reflection. The light escaping the receptor in the backward direction forms a Gaussian intensity distribution at the pupil plane, with the properties shown in Figs. 1 and 2. The volume, however, is taken as neutral.
Choi 47 stated, on the basis of three wavelengths (550, 650, and ), that the guided fraction under the 2-D OSCE is neutral with wavelength. Our extensive wavelength data (Fig. 10 ) follow the model curve and clearly show a strong wavelength dependency (note that absorption in the media is cancelled in the calculation); a maximum is at around . Such wavelength dependence is obvious because below , the directional cone signal is reduced by macular pigment, while the ILM part of the nondirectional reflection is not. At wavelengths above , the nondirectional reflection from the deeper layers increases because blood and melanin become more and more transparent.
The correlation between and of 0.545 could in principle be caused by a non-Gaussian behavior in the skirts of the directional reflection. If the theoretical function is low in that region compared to the real directional profile, it uses to fill the gap. Physiological interaction is less likely with the relatively high covariance of 0.694 in mind, an indication that the model likes to interact with these parameters. A practical solution in further applying the model could be fixing and using as the indicator of healthy cones.
For the present purpose with bleached visual pigments, might as well be taken as the outer-segment end reflection. That the discs are the origin could only be made plausible by experiments with dark adapting and bleaching of the visual pigments.5 Figure 9 shows a large decrease in the disc reflection as a function of age. That the large reflection from the choroid (8%) and that from the retinal pigment epithelium (0.56%) have no such a decrease is a strong argument for the receptors as origin. Because the width of the Gaussian as given by has no age dependency and there is neither evidence for a large decrease in the number of cones nor for a large decrease in length of the outer segments,17, 50 we might speculate that the indexes of refraction of the discs, and the interspaces between them, slowly drift toward each other at older ages. In diseases like age-related macular degeneration (AMD), might decrease at a faster rate because here the number of cones or the lengths of the outer segments are possibly affected.
A new model gives fundamental insight in to how the directional and nondirectional spectral reflection from the fovea are related to each other. Previous attempts were less successful because the transmission in the eye media was taken as equal in both components. Setting the Rayleigh scatter to zero in the nondirectional reflection component caused a dramatic improvement of the model fit. This part is, however, presently not fully understood; we could not find support from other sources in the literature for additional scatter losses in preretinal layers. Further work on this issue will be needed.
The model generates quantitative estimates of several important optical absorbers in the eye as macular pigment, lens density, blood, and melanin, as well as parameters related to the optical quality of the cones, the cone disc reflection, and . In addition, the model estimates the reflection from the retinal pigment epithelium and the choroid. A validation of the macular pigment optical density was obtained by the significant correlations obtained with independent estimates.
The authors are indebted to Dr. Tos Berendschot, Martijn Kanis, and Robert Wisse for allowing us to use their data.