Following data collection, the OCT images were exported from the instrument and analyzed using custom written software. Each of the two sets of scans for each subject were analyzed to segment the outer surface of the RPE, and the inner surface of the chorioscleral interface (CSI), to determine choroidal thickness across the 30° width of each scan. The RPE initially was segmented using an automated method based on graph theory.
30 An experienced masked observer then manually segmented the CSI using a method that has been described in detail previously,
23 and involves the observer manually selecting a series of points along the CSI. The software then automatically fits a smooth function (spline fit) to these points to define the boundary. Additionally, the observer also checked the integrity of the automated segmentation of the RPE and manually corrected any segmentation errors. The center of the fovea, defined as the deepest point in the central foveal pit, also was marked in each scan by the observer.
Following the segmentation of the OCT scans, the choroidal thickness data from each scan had the transverse scaling corrected to account for ocular magnification effects, based upon each participant's biometric and refraction measurements. The length in millimeters of each 30° scan was determined based upon the distance from the retina to the eye's second nodal point, assuming a spherical retina.
31 The position of the second nodal point for each individual subject was determined with the “step along” method,
32 using each subject's measured ocular refraction and ocular biometry measures, applying the methods outlined by Bennett
33 to determine the equivalent power of the eye and crystalline lens.
The thickness data subsequently was analyzed to determine the subfoveal choroidal thickness, and the average choroidal thickness across a series of concentric zones around the fovea, including the central foveal zone (central 1-mm diameter), the inner macula zone (from an inner diameter of 1 mm to an outer diameter of 3 mm), and the outer macula zone (inner diameter of 3 mm and outer diameter of 6 mm,
Fig. 1b). This analysis provided the average thickness at 8 locations (temporal, superior-temporal, superior, superior-nasal, nasal, inferior-nasal, inferior, and inferior-temporal), across each of the three zones (central fovea, inner macula, and outer macula).
To assess the reliability and repeatability of the thickness data derived from the manual segmentation of the OCT images, the scans from 20 randomly selected subjects were analyzed twice by the observer, who was masked to the choroidal thickness results from the initial analysis. The subfoveal and parafoveal (central fovea, inner macula, and outer macula zones) thickness data were analyzed using the methods described by Bland and Altman.
34 Since there appeared to be a trend for the differences between the two analyses to be slightly greater for thicker choroids, the mean differences and limits of agreement were determined based upon the ratios of the differences between the two analyses.
34
The influence of refractive error and sex upon subfoveal choroidal thickness was examined with a 2-way ANOVA. Additionally, a stepwise multiple regression was performed to examine the influence of demographic (age and sex) and biometric ocular factors (spherical equivalent refractive error, axial length, central corneal thickness, anterior chamber depth, and lens thickness) upon the subfoveal choroidal thickness measures. To examine the topographic distribution of choroidal thickness across the posterior pole, a repeated measures ANOVA was performed with two within-subjects factors, including choroidal location (temporal, superior-temporal, superior, superior-nasal, nasal, inferior-nasal, inferior, or inferior-temporal) and choroidal zone (central foveal, inner macula, and outer macula), and one between-subjects factor (refractive error group). The five children who exhibited poor fixation during OCT imaging were excluded from the topographic analysis of choroidal thickness. All results presented represent the mean ± SD.