A mathematical model for describing the retinal nerve fiber bundle trajectories in the human eye: average course, variability, and influence of refraction, optic disc size and optic disc position.

Previously we developed a mathematical model for describing the retinal nerve fiber bundle trajectories in the superior-temporal and inferior-temporal regions of the human retina, based on traced trajectories extracted from fundus photographs. Aims of the current study were to (i) validate the existing model, (ii) expand the model to the entire retina and (iii) determine the influence of refraction, optic disc size and optic disc position on the trajectories. A new set of fundus photographs was collected comprising 28 eyes of 28 subjects. From these 28 photographs, 625 trajectories were extracted. Trajectories in the temporal region of the retina were compared to the existing model. In this region, 347 of 399 trajectories (87%) were within the 95% central range of the existing model. The model was extended to the nasal region. With this extension, the model can now be applied to the entire retina that corresponds to the visual field as tested with standard automated perimetry (up to approximately 30° eccentricity). There was an asymmetry between the superior and inferior hemifields and a considerable location-specific inter-subject variability. In the nasal region, we found two "singularities", located roughly at the one and five o'clock positions for the right optic disc. Here, trajectories from relatively widespread areas of the retina converge. Associations between individual deviations from the model and refraction, optic disc size and optic disc position were studied with multiple linear regression. Refraction (P = 0.021) and possibly optic disc inclination (P = 0.09) influenced the trajectories in the superior-temporal region.