Abstract
Purpose:
To characterize the effect of sampling window size on maps of foveal cone density derived from adaptive optics scanning light ophthalmoscope (AOSLO) images of the cone mosaic.
Methods:
Forty-four AOSLO-derived montages of the foveal cone mosaic (300 x 300µm) were used for this study (from 44 individuals with normal vision). Cone photoreceptor coordinates were semi-automatically identified by one experienced grader. From these coordinates, cone density matrices across each foveal montage were derived using 10 different sampling window sizes containing 5, 10, 15, 20, 40, 60, 80, 100, 150, or 200 cones. For all 440 density matrices, we extracted the location and value of peak cone density (PCD), the cone density centroid (CDC) location, and cone density at the CDC.
Results:
Across all window sizes, PCD values were larger than those extracted at the CDC location, though the difference between these density values decreased as the sampling window size increased (p<0.0001). Overall, both PCD (r=-0.8099, p=0.0045) and density at the CDC (r=-0.7596, p=0.0108) decreased with increasing sampling window size. This reduction was more pronounced for PCD, with a 27.8% lower PCD value on average when using the 200-cone versus the 5-cone window (compared to only a 3.5% reduction for density at the CDC between these same window sizes). While the PCD and CDC locations did not occur at the same location within a given montage, there was no significant relationship between this PCD-CDC offset and sampling window size (p=0.8919). The CDC location was less variable across sampling windows, with an average per-participant 95% confidence ellipse area across the 10 window sizes of 47.56µm² (compared to 844.10µm² for the PCD location, p<0.0001).
Conclusion:
CDC metrics appear more stable across varying sampling window sizes than PCD metrics. Understanding how density values change according to the method used to sample the cone mosaic may facilitate comparing cone density data across different studies.
1 Introduction
Adaptive optics scanning light ophthalmoscopy (AOSLO) enables direct observation of the human photoreceptor mosaic with single-cell resolution (1–3). Of particular interest are the cone photoreceptors at the fovea – a highly specialized region of the retina that supports our high-acuity vision (4, 5). Relative to the parafoveal cone mosaic, foveal cone density is higher and the cones are contiguously packed (due in part to the absence of rod photoreceptors) (6). Given the vital importance of foveal cone photoreceptors for our vision, and the involvement of the fovea in many pathologies, there is significant interest in advancing quantitative biomarkers to assess the foveal cone mosaic. Numerous metrics exist for describing geometric properties of the cone mosaic, though cone density and cone spacing are used most often (7–10). Despite convergence on these metrics to characterize the foveal cone mosaic, there remain differences in the methods used to derive them – limiting the ability to compare results from different studies and ultimately limiting progress toward the development of robust clinical biomarkers.
A frequent inconsistency between studies relates to the size and shape of the sampling window used to derive metrics of the cone mosaic. The steep cone density gradient in the central fovea requires special consideration of how large of an area to use when deriving a local estimate of cone density (6). Larger sampling windows will result in lower values of peak density due to inclusion of lower density neighboring regions, while smaller sampling windows result in more noisy density maps that can make derivation of the global peak cone density (PCD) location challenging. Despite this, fixed-area sampling windows have been used in many studies (6, 8, 9, 11–20). This issue is partially overcome by using sampling windows that include a fixed number of cones, regardless of the retinal location. While this approach has been used by some investigators (10, 21–24), there is inconsistency in the number of cones sampled. Here we examined the relationship between the number of cones included in the sampling window and cone mosaic metrics (location and value of PCD, location of the cone density centroid (CDC), and density at the CDC). We also report the effect of variable sampling windows on the overall cone density topography for individual participants.
2 Methods
2.1 Participants
This study followed the tenets of the Declaration of Helsinki and was approved by the Institutional Review Board at the Medical College of Wisconsin (PRO 30741). Included in this study were right eyes of 44 participants (age range: 12-61 years; 17 males and 27 females) with self-reported normal vision and previously imaged as part of Cava et al. (2020) and Wynne et al. (2022) (23, 24). Prior to AOSLO imaging, axial length (IOL Master, Carl Zeiss Meditec, Dublin CA, USA) and autorefraction (KR-800S, Topcon Corporation, Tokyo, Japan) measurements were taken. The study eye was then dilated and accommodation suspended with one drop of 2.5% phenylephrine hydrochloride (Akorn, Lake Forest, IL, USA) and one drop of 1% tropicamide (Akorn, Lake Forest, IL, USA) for participants aged 18 years and older. For participants under age 18, one drop of 1% cyclopentolate hydrochloride (Cyclomydril, Alcon Laboratories Inc., Fort Worth, TX, USA) was used for dilation and accommodation suspension.
2.2 AOSLO imaging, processing, and montaging
A dental impression bite bar was used for head stabilization during AOSLO imaging. Per our imaging protocol, confocal AOSLO videos consisting of 150-200 frames each were acquired at various locations within the fovea. Videos were taken at 1.5° field of view (FOV) with approximately 1° overlap between video locations or at 1.0° FOV with approximately 0.5° overlap. A 775nm or 790nm super luminescent diode (SLD) was used for imaging. Several techniques were used to improve the resolution of foveal cones during imaging. These techniques included imaging with 0.75° and 0.5° FOVs, using a 680nm SLD (incident power: 32.5µW), or using a sub-airy disk pinhole (0.5-0.7 Airy disk diameter) (24).
Each video was processed to create high-resolution .tif images as previously described (24). From each video, a minimally distorted reference frame was automatically chosen (25) and used to register and average the remaining frames in the video with a strip-based registration algorithm (25, 26). Frames were then repaired to remove additional distortion with a de-warping repair script (https://github.com/OCVL/Eye-Motion-Repair) based on a previous method (27). This repair script works by estimating the bias of random eye motion distortions throughout the reference frames, based on the median translation observed in each row of the registered images. The median translation fixes the distortion of the frame in a same-magnitude but opposite-direction approach (24, 28). These processing techniques create a .tif image with high signal-to-noise ratio from each acquired video. Additional improvement in the resolution of the central most foveal cones was possible by averaging multiple processed images from videos acquired at the same retinal location – either at different focal planes or simply at different time points. This approach reduces the between cell variation in reflectivity, making discrimination of cones in the mosaic easier (24, 29).
The scale of each montage (µm/pixel) was determined from the AOSLO system scale as well as the individual’s axial length (2) at the time of image acquisition. This calculation was previously described (19) using the following equation:
where T represents the periodicity of a Ronchi ruling (µm/cycle), f1 represents the focal length of the model eye in our system (µm), Ts represents the sampling period of the lines in the model eye image of the Ronchi ruling (pixels/cycle), RMF represents the assumed retinal magnification factor (291 µm/degree) of an eye with a 24.0mm axial length (30), and AL is the participant’s axial length (mm).
A custom MATLAB script (https://github.com/BrainardLab/AOAutomontaging) was used to automatically align the processed .tif images (31) into a single montage. Alignment of overlapping images was manually corrected using Adobe Photoshop CS6 (Adobe Systems, Inc., San Jose, CA, USA). Individual frames were blended to create a flattened, seamless image of the area containing the subjective peak cone density as well as its adjacent frames. Finally, a 300 x 300µm region of interest (ROI) centered at the subjective location of approximate highest cone density (locus of smallest cones/tightest packing) was cropped for analysis.
2.3 Cone counting and density mapping
Cones were identified using a semi-automated cone counting software (Mosaic Analytics; Translation Imaging Innovations, Hickory, NC, USA). One experienced human grader (JC) manually corrected the cone markings automatically identified by the algorithm. For each foveal montage, a cone coordinate file was generated containing the (x,y) locations of the individual cone markings. For this study we analyzed the effect of sampling window size when determining bound cone density for individual foveal montages. We used a square sampling window and specified the number of cells to be included in each window using a custom density mapping script (https://github.com/AOIPLab/Metricks/releases/tag/Warr_et_al_2024). The criterion sampling window size for generating the density maps was adjusted to include 5, 10, 15, 20, 40, 60, 80, 100, 150, or 200 cells. The sampling window expands at each coordinate until it includes at least the specified number of bound Voronoi cells. Due to the irregularity of Voronoi cells at the fovea, the sampling window is often unable to be sized such that it exactly matches the specified number of bound Voronoi cells. In these instances, we set the sampling window to the smallest size that still exceeds the specified number of cells, then randomly remove Voronoi cells around the edge of the sampling window until the exact number of specified cells is achieved. For edge cells in the montage, the sampling window expands on two or three sides toward the center of the montage instead of expanding proportionally from the target coordinate (this way, density is only calculated within the bounds of the montage). The density at each coordinate is calculated by taking the number of bound cells divided by the sum of their Voronoi areas. The window moves to each coordinate in the coordinate file and a pixelwise linear interpolated density map is created for each montage (Figure 1).
Figure 1
The interpolated pixel in each foveal montage containing the absolute maximum density is the location of the PCD and the subsequent density value at that point is taken as the PCD. The CDC was determined by generating an 80% isodensity contour of the region that includes pixels containing the top 20th percentile of density values within the foveal montage. We then fit an ellipse around this 80% isodensity region and define the center of this best-fit ellipse as the CDC location (see Supplementary Figure 1). This method for determining the CDC is based off the method described by Reiniger et al. (2021) (32). We determined PCD and CDC metrics from each sampling window condition across all 440 density maps (44 participants with 10 density maps each). The 10 density maps for one participant are shown in Figure 2. To assess intra-individual reproducibility of the PCD and CDC locations, we computed the 95% confidence ellipse areas for both PCD and CDC locations from all 10 window sizes using a custom MATLAB script (https://github.com/AOIPLab/Metricks/releases/tag/Warr_et_al_2024). Example PCD and CDC 95% confidence ellipses for a single participant are shown in Figure 2K.
Figure 2
2.4 Assessing overall cone mosaic topography
To assess how the sampling window affects overall aspects of foveal cone topography (not just PCD and CDC), we created two separate composite maps (mean density and density standard deviation) for each participant using the data from their 10 individual density maps. From these two maps, we extracted a horizontal and vertical cross section through the average CDC location (derived from that participant’s 10 individual density maps). For most participants, the average CDC (x,y) coordinates were not whole numbers. In these cases, we used a weighted average of the values in the rows and/or columns neighboring the average CDC point (for example, if the average CDC x-coordinate was 150.3, the 150th column of values were weighted at 70% while the 151st column of values were weighted at 30%). Shown in Figure 3 are topographical plots from both cross-sections illustrating the average trend from all 44 participants. All MATLAB scripts utilized above can be found at: https://github.com/AOIPLab/Metricks/releases/tag/Warr_et_al_2024.
Figure 3
2.5 Statistics
Linear units were used to represent PCD and CDC density (cones/mm²) as well as PCD and CDC locations (µm). We also repeated the comparisons using angular units to represent density (cones/deg²). See Supplementary Table 1 for participant level data. A Shapiro-Wilk normality test was used to determine the use of parametric or nonparametric approaches. All statistics were calculated using GraphPad Prism (Prism 9.0.0; GraphPad Software, San Diego, CA, USA).
3 Results
Across all window sizes, average PCD values were larger than those extracted at the CDC location (p=0.0012, paired t test; Figure 4A). On average, PCD and CDC both decreased with increasing sampling window size, though the effect was larger for PCD (r=-0.8099, p=0.0045) than CDC (r=-0.7596, p=0.0108). As the density values at the CDC with the smallest window sizes behaved in an inconsistent fashion, we also examined the relationship between density and window size using only window sizes of 20 cones and larger (Figure 4B). We again see a strong relationship between PCD and window size (r=-0.9192, p=0.0034) and also for CDC density values (r=-0.8595, p=0.0132). As summarized in Supplementary Table 2, our density values are comparable to previous estimates in individuals with normal vision (10, 13, 21, 23, 24, 33–36). For example, with a window size of 100 cones, the average PCD from all 44 participants was 185,501 cones/mm² (range=128,007-244,921 cones/mm²) and the average density at the CDC was 174,696 cones/mm² (range=121,169-242,614 cones/mm²). At a window size of 150 cones, the average PCD from all 44 participants was 182,184 cones/mm² (range=126,323-239,362 cones/mm²) and the average density at the CDC was 174,295 cones/mm² (range=121,891-235,663 cones/mm²). The intra-participant differences in density estimates at the PCD and CDC locations decreased significantly as window size increased (p<0.0001, ANOVA, Figure 4C). This was driven largely by the increased PCD values obtained with smaller sampling windows (Tukey’s multiple comparisons test, see Supplementary Table 3). The same trends were observed when reporting density in angular units (cones/deg2), see Supplementary Figure 2. Across, the 44 participants, as the number of cones in the sampling window increase, we observed a proportional increase in window area (y=9.21x + 1.7146, R² = 1; where x is the window size in number of cones and y is the average bound area window size across all participants). For example, the average area for the 200-cone sampling window was 1,841.14µm² while the average area for a 5-cone sampling window was 45.98µm² (see Supplementary Table 4).
Figure 4
In addition to density values, we observed variation in PCD and CDC location as a function of window size. Across all 440 ROIs (44 participants with 10 window sizes each), there were no instances where the reported PCD value occurred at more than one point. The average PCD confidence ellipse area across all participants was 844.10µm2, which was significantly greater than the average CDC confidence ellipse area (47.56µm2, p<0.0001, Wilcoxon matched-pairs signed-rank test). PCD 95% confidence ellipse areas were greater than CDC 95% confidence ellipse areas for all participants with one exception (JC_11830). Across the 44 participants, PCD 95% confidence ellipse area ranged from 105.28 to 8,389.90µm2, while CDC 95% confidence ellipse area ranged from 2.50 to 481.73µm2. Example 95% confidence ellipses are shown for four participants in Figure 5. There was variation between participants in the relationship between the magnitude of the offset between PCD and CDC location and sampling window size. On average, the 44 participants demonstrate relative stability in PCD-CDC location offset as sampling window size increases (Figure 6A). However, when assessed individually, some participants showed convergence of the PCD and CDC location with increasing sampling window, some showed divergence, while others showed inconsistent trends (Figure 6B). This heterogeneity in how window size impacts PCD-CDC offset is likely due to the underlying global topography of each cone mosaic.
Figure 5
Consistent with prior studies of cone topography (6, 22, 35), we observed a decline in average cone density as a function of eccentricity in both the horizontal and vertical meridians (Figure 3). Additionally, we observed greater within-participant standard deviation near the CDC location compared to locations toward the edge of the foveal ROI – indicating that the variability induced by using different sampling windows affects density more in areas of rapidly changing density (foveal center) compared to areas where density is more uniform (edge of the foveal ROI used here).
Figure 6
4 Discussion
Consistent with previous studies of the parafoveal cone mosaic (16, 37–39), our results highlight the importance of defining the sampling window used in deriving maps of foveal cone density. Importantly, we find that the value and location of PCD (versus CDC) is more significantly affected by using a variable sampling window size. This is demonstrated by the steeper decline in PCD as window size increases as compared to density at the CDC (Figure 4A). Additionally, the smaller average CDC 95% confidence ellipse area compared to PCD 95% confidence ellipse area demonstrates greater stability in CDC location compared to PCD location as sampling window size changes. These findings are expected based on how we define PCD versus density at the CDC. PCD is the interpolated point of absolute highest density in a given montage whereas density at the CDC describes a relative maximum based on a pooled region of highly dense locations. Thus, when the sampling window is small, we can resolve several high-density regions in the foveal mosaic, visible as small “islands” of 80% isodensity contour (Figures 2A–D). To derive the CDC, we generate a best fit ellipse of the 80% isodensity contour and define the center of that ellipse as the CDC location. This method identifies a point whose density is an average of a region of high cone density. Therefore, despite changes in sampling window size, the CDC location of an individual ROI remains more stable than the PCD location. As well, previous studies have demonstrated superior intersession repeatability for CDC location compared to PCD location (32) as well as superior intergrader reproducibility in CDC location relative to PCD location (23). Taken together, the emerging picture is that CDC represents a more robust metric to describe the “center” of the foveal cone mosaic.
As reported by Wang et al. (2019), use of smaller sampling windows may produce unreliable estimates of PCD and PCD location (34), though it is important to note they used a fixed dimension sampling window, which may contain variable numbers of cones across participants. Additionally, their sampling window of 10 arcmin contains over 200 cones, which would be larger than the largest sampling window used in our study. Our data support the notion that smaller window sizes impact PCD and we show that this also applies to the CDC (albeit to a lesser degree). When excluding the smallest sampling window sizes, we observe a strong relationship between sampling window size and both PCD and density at CDC (Figure 4B). This observation may allow comparison of density estimates between studies that used different sampling windows to construct their respective cone density maps. However, while this approach may work on average, there is variation between individual cone mosaics in how density estimates and locations vary as a function of sampling window size (Figure 6B). This is likely due to more global features of cone topography, including individual variation in kurtosis of the cone density versus eccentricity function. Further work is needed to understand if and how best to combine data across studies employing different methodology, as large and robust normative databases are needed to advance the use of AOSLO imaging.
Limitations of our study include reliance on a single observer’s cone markings – errors or biases in their coordinates may exist (23). As such, generalization of these trends should be applied carefully to other datasets where coordinates were derived from different observers. Additional factors related to sampling the cone mosaic could be considered in future studies – such as inclusion of “edge” cells within the window, shape of the window (square, arcuate, circular), and interpolation methods. Additionally, our sample of 44 participants was almost 2:1 female (27F, 17M) and almost exclusively white (approximately 70% of participants self-reported race as white). As it is unknown whether cone mosaic topography varies with race and ethnicity, our analysis should be repeated in a more heterogenous population.
Ultimately, the “cost” of variation in any metric (including those assessed here: PCD/CDC value and PCD/CDC location) depends on the subsequent application. For this reason, we do not identify one sampling window as the ideal size when extracting density metrics. Correlation of cone density with measures of function (e.g., acuity or sensitivity) may wish to utilize sampling windows with sizes comparable to some feature of the visual system (such as fixational stability) (34). Moreover, correlation with measures of function would be best made using cone density at the retinal location used for the specific task (40–42), which may not occur at the PCD or the CDC (32, 43, 44). Thus, in these scenarios, the PCD or CDC location may not be overly relevant. Likewise, correlation of mosaic metrics to other structural features in other imaging modalities would require co-registration of images to ensure alignment of the cone mosaic measures with these other retinal features (45, 46). One application that depends on having a reliable “anchor” within the cone mosaic is selecting parafoveal and perifoveal regions of interest from a larger montage for subsequent analysis. Whether PCD or CDC, this anchor serves as the (0,0) reference point from which retinal eccentricity of a given region of interest is computed (47). Though we observe an offset in PCD and CDC location at all window sizes, if choosing PCD or CDC as an anchor point this offset would remain a fixed value at all retinal eccentricities. However, as cone (and rod) density changes with eccentricity (6, 35), this fixed offset in the location assigned to a given region of interest could result in variably misleading conclusions about the health of the underlying mosaic. Our findings support the use of CDC as a more reliable anchor for such applications, even where resolution of the most densely packed cones is possible and PCD measures are desired for different purposes.
Statements
Data availability statement
The raw data (foveal montages, cone coordinate files, density matrices) supporting the conclusions of this article will be made available upon request to the corresponding author (jcarroll@mcw.edu), without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical College of Wisconsin IRB. The studies were conducted in accordance with local legislation and institutional requirements. The participants (or when appropriate, their parent’s or guardian’s) provided their written informed consent to participate in this study.
Author contributions
EW: Conceptualization, Formal analysis, Visualization, Writing – original draft, Writing – review & editing. JG: Software, Writing – review & editing. RFC: Funding acquisition, Software, Writing – review & editing. JC: Conceptualization, Formal analysis, Funding acquisition, Supervision, Visualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The project described was supported by National Center for Advancing Translational Sciences, National Institutes of Health (NIH), award number UL1TR001436. This study was supported by the National Eye Institute of the NIH under award numbers R01EY017607, R01EY033580, R44EY031278. Funds also came from the Foundation Fighting Blindness through FFB-BR-CL-0720-0784-MCW and FFB-CC-CL-0620-0785-MRQ and The Gene and Ruth Posner Foundation. This study was conducted in part in a facility constructed with support from Research Facilities Improvement Program, grant number C06RR016511 from the National Center for Research Resources, NIH.
Acknowledgments
The authors would like to thank Drs. Heather Heitkotter and Niamh Wynne for their assistance with this study.
Conflict of interest
Authors RFC and JC both have equity interest in company Translational Imaging Innovations.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor WH declared a past co-authorship with one of the authors JC.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Author disclaimer
The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fopht.2024.1348950/full#supplementary-material
References
1
RoordaARomero-BorjaFDonnellyWJIIIQueenerHHebertTCampbellM. Adaptive optics scanning laser ophthalmoscopy. Opt Express. (2002) 10:405–12. doi: 10.1364/oe.10.000405
2
DubraASulaiYNorrisJLCooperRFDubisAMWilliamsDRet al. Noninvasive imaging of the human rod photoreceptor mosaic using a confocal adaptive optics scanning ophthalmoscope. BioMed Opt Express. (2011) 2:1864–76. doi: 10.1364/BOE.2.001864
3
ScolesDSulaiYNLangloCSFishmanGACurcioCACarrollJet al. In vivo imaging of human cone photoreceptor inner segments. Invest Ophthalmol Vis Sci. (2014) 55:4244–51. doi: 10.1167/iovs.14-14542
4
ProvisJMHendricksonAE. The foveal avascular region of developing human retina. Arch Ophthalmol. (2008) 126:507–11. doi: 10.1001/archopht.126.4.507
5
TutenWSHarmeningWM. Foveal vision. Curr Biol. (2021) 31:R701–R3. doi: 10.1016/j.cub.2021.03.097
6
CurcioCASloanKRKalinaREHendricksonAE. Human photoreceptor topography. J Comp Neurol. (1990) 292:497–523. doi: 10.1002/cne.902920402
7
LittsKMCooperRFDuncanJLCarrollJ. Photoreceptor-based biomarkers in aoslo retinal imaging. Invest Ophthalmol Vis Sci. (2017) 58:255–67. doi: 10.1167/iovs.17-21868
8
MirhajianmoghadamHJnawaliAMusialGQueenerHMPatelNBOstrinLAet al. In vivo assessment of foveal geometry and cone photoreceptor density and spacing in children. Sci Rep. (2020) 10:8942. doi: 10.1038/s41598-020-65645-2
9
SongHChuiTYZhongZElsnerAEBurnsSA. Variation of cone photoreceptor packing density with retinal eccentricity and age. Invest Ophthalmol Vis Sci. (2011) 52:7376–84. doi: 10.1167/iovs.11-7199
10
CooperRFWilkMATarimaSCarrollJ. Evaluating descriptive metrics of the human cone mosaic. Invest Ophthalmol Vis Sci. (2016) 57:2992–3001. doi: 10.1167/iovs.16-19072
11
DeesEWDubraABaraasRC. Variability in parafoveal cone mosaic in normal trichromatic individuals. BioMed Opt Express. (2011) 2:1351–8. doi: 10.1364/BOE.2.001351
12
ChuiTYSongHBurnsSA. Individual variations in human cone photoreceptor packing density: variations with refractive error. Invest Ophthalmol Vis Sci. (2008) 49:4679–87. doi: 10.1167/iovs.08-2135
13
PutnamNMHoferHJDobleNChenLCarrollJWilliamsDR. The locus of fixation and the foveal cone mosaic. J Vis. (2005) 5:632–9. doi: 10.1167/5.7.3
14
ChuiTYPSongHXBurnsSA. Adaptive-optics imaging of human cone photoreceptor distribution. J Opt Soc Am A. (2008) 25:3021–9. doi: 10.1364/JOSAA.25.003021
15
LombardoMSerraoSDucoliPLombardoG. Adaptive optics photoreceptor imaging. Ophthalmology. (2012) 119:1498–.e2. doi: 10.1016/j.ophtha.2012.03.019
16
LombardoMSerraoSDucoliPLombardoG. Influence of sampling window size and orientation on parafoveal cone packing density. BioMed Opt Express. (2013) 4:1318–31. doi: 10.1364/BOE.4.001318
17
ParkSPChungJKGreensteinVTsangSHChangS. A study of factors affecting the human cone photoreceptor density measured by adaptive optics scanning laser ophthalmoscope. Exp Eye Res. (2013) 108:1–9. doi: 10.1016/j.exer.2012.12.011
18
LombardoMSerraoSLombardoG. Technical factors influencing cone packing density estimates in adaptive optics flood illuminated retinal images. PloS One. (2014) 9:e107402. doi: 10.1371/journal.pone.0107402
19
LittsKMGeorgiouMLangloCSPattersonEJMasteyRRKalitzeosAet al. Interocular symmetry of foveal cone topography in congenital achromatopsia. Curr Eye Res. (2020) 45:1257–64. doi: 10.1080/02713683.2020.1737138
20
LangloCSPattersonEJHigginsBPSummerfeltPRazeenMMErkerLRet al. Residual foveal cone structure in cngb3-associated achromatopsia. Invest Ophthalmol Vis Sci. (2016) 57:3984–95. doi: 10.1167/iovs.16-19313
21
LiKYTiruveedhulaPRoordaA. Intersubject variability of foveal cone photoreceptor density in relation to eye length. Invest Ophthalmol Vis Sci. (2010) 51:6858–67. doi: 10.1167/iovs.10-5499
22
CurcioCASloanKR. Packing geometry of human cone photoreceptors: variation with eccentricity and evidence for local anisotropy. Vis Neurosci. (1992) 9:169–80. doi: 10.1017/S0952523800009639
23
WynneNCavaJAGaffneyMHeitkotterHScheidtAReinigerJLet al. Intergrader agreement of photoreceptor topography at the foveal center on adaptive optics scanning light ophthalmoscopy. BioMed Opt Express. (2022) 13:4445–54. doi: 10.1364/BOE.460821
24
CavaJAAllphinMTMasteyRRGaffneyMLindermanRECooperRFet al. Assessing interocular symmetry of the foveal cone mosaic. Invest Ophthalmol Vis Sci. (2020) 61:23. doi: 10.1167/iovs.61.14.23
25
SalmonAECooperRFLangloCSBaghaieADubraACarrollJ. An automated reference frame selection (Arfs) algorithm for cone imaging with adaptive optics scanning light ophthalmoscopy. Transl Vis Sci Tech. (2017) 6:9. doi: 10.1167/tvst.6.2.9
26
DubraAHarveyZ. Registration of 2d images from fast scanning ophthalmic instruments. In: FischerBDawantBLorenzC, editors. Biomedical image registration. Lecture notes in computer science, 1. Springer-Verlag, Berlin (2010). p. 60–71.
27
BedggoodPMethaA. De-warping of images and improved eye tracking for the scanning laser ophthalmoscope. PloS One. (2017) 12:e0174617. doi: 10.1371/journal.pone.0174617
28
ChenMCooperRFGeeJCBrainardDHMorganJIW. Automatic longitudinal montaging of adaptive optics retinal images using constellation matching. BioMed Opt Express. (2019) 10:6476–96. doi: 10.1364/boe.10.006476
29
CooperRFDubisAMPavaskarARhaJDubraACarrollJ. Spatial and temporal variation of rod photoreceptor reflectance in the human retina. BioMed Opt Express. (2011) 2:2577–89. doi: 10.1364/BOE.2.002577
30
EmsleyH. Visual optics. 5th ed. London, UK: Butterworth & Co (1953), p. 1–444.
31
ChenMCooperRFHanGKGeeJBrainardDHMorganJI. Multi-modal automatic montaging of adaptive optics retinal images. BioMed Opt Express. (2016) 7:4899–918. doi: 10.1364/BOE.7.004899
32
ReinigerJLDomdeiNHolzFGHarmeningW. Human gaze is systematically offset from the center of cone topography. Curr Biol. (2021) 31:1–6. doi: 10.1016/j.cub.2021.07.005
33
ZhangTGodaraPBlancoERGriffinRLWangXCurcioCAet al. Variability in human cone topography assessed by adaptive optics scanning laser ophthalmoscopy. Am J Ophthalmol. (2015) 160:290–300. doi: 10.1016/j.ajo.2015.04.034
34
WangYBensaidNTiruveedhulaPMaJRavikumarSRoordaA. Human foveal cone photoreceptor topography and its dependence on eye length. eLife. (2019) 8:e47148. doi: 10.1101/589135
35
Wells-GrayEMChoiSSBriesADobleN. Variation in rod and cone density from the fovea to the mid-periphery in healthy human retinas using adaptive optics scanning laser ophthalmoscopy. Eye. (2016) 30:1135–43. doi: 10.1038/eye.2016.107
36
WilkMAWilkBMLangloCSCooperRFCarrollJ. Evaluating outer segment length as a surrogate measure of peak foveal cone density. Vision Res. (2017) 130:57–66. doi: 10.1016/j.visres.2016.10.012
37
FengSGaleMJFayJDFaridiATitusHEGargAKet al. Assessment of different sampling methods for measuring and representing macular cone density using flood-illuminated adaptive optics. Invest Ophthalmol Vis Sci. (2015) 56:5751–63. doi: 10.1167/iovs.15-16954
38
GarriochRLangloCDubisAMCooperRFDubraACarrollJ. Repeatability of in vivo parafoveal cone density and spacing measurements. Optom Vis Sci. (2012) 89:632–43. doi: 10.1097/OPX.0b013e3182540562
39
GianniniDLombardoGMariottiLDevaneyNSerraoSLombardoM. Reliability and agreement between metrics of cone spacing in adaptive optics images of the human retinal photoreceptor mosaic. Invest Ophthalmol Vis Sci. (2017) 58:3127–37. doi: 10.1167/iovs.16-20890
40
RatnamKDomdeiNHarmeningWMRoordaA. Benefits of retinal image motion at the limits of spatial vision. J Vis. (2017) 17:30. doi: 10.1167/17.1.30
41
RossiEARoordaA. The relationship between visual resolution and cone spacing in the human fovea. Nat Neurosci. (2010) 13:156–7. doi: 10.1038/nn.2465
42
ClarkAMIntoyJRucciMPolettiM. Eye drift during fixation predicts visual acuity. Proc Natl Acad Sci USA. (2022) 119:e2200256119. doi: 10.1073/pnas.2200256119
43
WilkMADubisAMCooperRFSummerfeltPDubraACarrollJ. Assessing the spatial relationship between fixation and foveal specializations. Vision Res. (2017) 132:53–61. doi: 10.1016/j.visres.2016.05.001
44
DomdeiNReinigerJLHolzFGHarmeningW. The relationship between visual sensitivity and eccentricity, cone density and outer segment length in the human foveola. Invest Ophthalmol Vis Sci. (2021) 62:31. doi: 10.1167/iovs.62.9.31
45
FooteKGRinellaNTangJBensaidNZhouHZhangQet al. Cone structure persists beyond margins of short-wavelength autofluorescence in choroideremia. Invest Ophthalmol Vis Sci. (2019) 60:4931–42. doi: 10.1167/iovs.19-27979
46
SunLWJohnsonRDLangloCSCooperRFRazeenMMRussilloMCet al. Assessing photoreceptor structure in retinitis pigmentosa and usher syndrome. Invest Ophthalmol Vis Sci. (2016) 57:2428–42. doi: 10.1167/iovs.15-18246
47
RoshandelDSampsonDMMackeyDAChenFK. Impact of reference center choice on adaptive optics imaging cone mosaic analysis. Invest Ophthalmol Vis Sci. (2022) 63:12. doi: 10.1167/iovs.63.4.12
Summary
Keywords
fovea, cone density, cone spacing, adaptive optics, retina
Citation
Warr E, Grieshop J, Cooper RF and Carroll J (2024) The effect of sampling window size on topographical maps of foveal cone density. Front. Ophthalmol. 4:1348950. doi: 10.3389/fopht.2024.1348950
Received
03 December 2023
Accepted
13 March 2024
Published
09 April 2024
Volume
4 - 2024
Edited by
Wolf Harmening, University Hospital Bonn, Germany
Reviewed by
Rigmor C. Baraas, University of South-Eastern Norway (USN), Norway
Julius Ameln, University of Bonn, Germany
Updates
Copyright
© 2024 Warr, Grieshop, Cooper and Carroll.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joseph Carroll, jcarroll@mcw.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.