Klin Monbl Augenheilkd 2003; 220(4): 235-246
DOI: 10.1055/s-2003-38625
Klinische Studie
© Georg Thieme Verlag Stuttgart · New York

Raytracing von Hornhauttopographiedaten zur Ermittlung der optischen Abbildungsqualität des Auges

Assessment of the Optical Image Quality of the Eye Using Raytracing Technique of Corneal Topography DataAchim  Langenbucher1 , Tomas  Sauer2 , Gerhard  J.  van der Heyd3 , Anja  Viestenz1 , Berthold  Seitz1
  • 1Augenklinik mit Poliklinik der Universität Erlangen-Nürnberg, Erlangen (Vorstand: Prof. Dr. med. Dr. h. c. mult. G. O. H. Naumann)
  • 2Mathematisches Institut der Justus-Liebig-Universität Gießen
  • 3Lehrstuhl für Fertigungstechnologie der Universität Erlangen-Nürnberg, Erlangen (Vorstand: Prof. Dr.-Ing. Dr. h. c. mult. M. Geiger)
Die Autoren haben kein finanzielles Interesse an der vorgelegten Studie
Further Information

Publication History

Eingegangen: 13. November 2002

Angenommen: 7. Februar 2003

Publication Date:
15 April 2003 (online)

Zusammenfassung

Hintergrund: Optische Aberrationen im Übertragungskanal beeinträchtigen die Abbildungsqualität des Auges in einer Weise, wie sie durch sphärozylindrische Gläser nicht vollständig ausgeglichen werden können. Die zu erwartende korrigierte Sehschärfe ist somit reduziert. Ziel dieser Studie war, die optische Abbildung einer punktförmigen Lichtquelle sowie eines ausgedehnten Objektes unter Verwendung computertopographischer Messdaten mit exaktem Raytracing abzuleiten. Methoden: Aus den Rohhöhendaten des Hornhauttopographiesystems (TMS-1, Tomey, Erlangen bzw. Orbscan, Orbtec Inc, USA) wurde mittels interpolierendem Subdivisionsschema (modifiziertes Butterfly-Schema) eine konvexe einfach differenzierbare (C1) refraktive Grenzfläche (Hornhautvorderfläche bzw. Vorder- und Rückfläche) berechnet. Die Charakteristika der verbleibenden Grenzflächen wurden dem Augenmodell nach Navarro entnommen. Aus der Berechnung der Punktstreuung (sog. „point-spread function”) mittels exaktem Raytracing (Anwendung des Snellius-Gesetzes) wurde die Fokusebene individuell durch eine Minimierung der Varianz ermittelt. Die Beugung an der Apertur wurde in Form einer Transmissionscharakteristik (radialsymmetrische Bessel-Funktion) an der Eintrittspupille des Systems realisiert. Das Schema wurde in der Programmiersprache C unter LINUX entwickelt und auf einen Normalbefund (Beispiel 1, TMS-1), einen Befund bei Keratokonus im fortgeschrittenem Stadium (Beispiel 2, TMS-1) und einen Befund mit Hornhautnarben (Beispiel 3, Orbscan) angewandt. Ergebnisse: Im Beispiel 1 vergrößerte sich der Fokusabstand (22,5 mm, 22,6 mm und 22,8 mm) geringfügig mit der Pupillenweite (2 mm, 3 mm und 5 mm). Die Varianz der nahezu radialsymmetrischen „point-spread function” in der Fokusebene erreichte einen minimalen Wert bei einer Pupillenweite von 3 mm (0,164, 0,104 und 0,230). Im Beispiel 2 veränderte sich aufgrund der optischen Aberrationen (v. a. Verkippung, Koma und Trifoil) der Fokusabstand uneinheitlich (21,1 mm, 21,0 mm und 21,3 mm) mit der Pupillenweite (2 mm; 3 mm und 5 mm). Die Varianz der deutlich asymmetrischen „point-spread function” in der Fokusebene lag durchwegs über den Vergleichswerten bei Beispiel 1 und erreichte einen minimalen Wert bei einer Pupillenweite von 3 mm (0,255, 0,224 und 0,371). Das optische Auflösungsvermögen (sinusförmig moduliertes Streifenmuster) war aufgrund der Unsymmetrie der „point-spread function” abhängig von der Orientierung (anisotrop). Im Beispiel 3 veränderte sich der Fokusabstand (22,3 mm, 22,3 mm und 22,5 mm) unwesentlich mit der Pupillenweite (2 mm, 3 mm und 5 mm). Die Varianz der „point-spread function” in der Fokusebene lag unter den Vergleichswerten aus Beispiel 2 und änderte sich unwesentlich zwischen Pupillenweiten von 2 und 3 mm (0,231, 0,239 und 0,338). Schlussfolgerung: Raytracing von topographischen Höhendaten auf der Basis eines Augenmodells mit der Option einer Autofokusfunktion bietet die Möglichkeit, das optische Auflösungsvermögen des Auges für punktförmige oder endlich ausgedehnte Lichtquellen abzuleiten. Weitere Studien über das Kontrastempfindlichkeitsverhalten und die Übertragung des optischen Bildes ins Gehirn sind nötig für eine vollständige Modellierung der subjektiven Sehschärfe.

Abstract

Background: Optical aberrations in the optical system may downgrade image quality and cannot be fully compensated by spherocylindrical glasses. The subjectively evaluated visual acuity may be significantly reduced. The purpose of this study was to calculate the image forming properties of the eye using a spotlight source or alternatively extended objects. Methods: A convex and first derivative continuous (C1) surface from the rough height data of the anterior corneal surface (TMS-1, Tomey, Erlangen) or the anterior and posterior corneal surface (Orbscan, Orbtec, USA) was calculated by means of an interpolating subdivision scheme (modified Butterfly algorithm). The characteristics of the residual refractive surfaces were used according to Navarro's eye model. The focal distance was calculated from the exact raytracing calculation (Snellius' law) of the point-spread function by minimising the variance of the point-spread function. The diffraction property of the aperture stop was implemented with a transmission characteristic according to a radially symmetrical Bessel function within the entrance pupil. The algorithm was realised with a C code on the LINUX platform and applied to a normal eye (example 1, TMS-1), an eye with severe keratoconus (example 2, TMS-1) and an eye with corneal scars (example 3, Orbscan). Results: The focal distance in example 1 (22.5 mm, 22.6 mm, and 22.8 mm) increased with the pupil diameter (2 mm, 3 mm, and 5 mm). The variance of the approximately radially symmetrical point-spread function in the focal plane attained a minimum value with a pupil size of 3 mm (0.164, 0.104, and 0.230). In example 2, the focal distance changed inconclusively (21.1 mm, 21.0 mm, and 21.3 mm) with the pupil size (2 mm, 3 mm, and 5 mm). The variance of the markedly asymmetrical point-spread function in the focal plane was systematically higher compared to the values of example 1 and reached a minimum value with a pupil size of 3 mm (0.255, 0.224, and 0.371). The imaging of the sinus-modulated pattern is anisotropic due to the asymmetry of the point-spread function. In example 3, the focal distance (22.3 mm, 22.3 mm, and 22.5 mm) did not change systematically with the pupil size (2 mm, 3 mm, and 5 mm). The variance of the nearly radially symmetrical point-spread function changed only marginally between pupil sizes of 2 mm and 3 mm (0.231, 0.239, and 0.338). Conclusions: Raytracing of corneal topography height data based on refined eye models with the option of auto-focussing has the potential to trace the optical resolution of the eye for arbitrary objects. Further studies on contrast sensitivity and the conversion of the real image to a perceived image by the retina and brain are required for complete modeling of subjective visual acuity.

Literatur

  • 1 Atchinson D A, Smith G, Waterworth M D. Theoretical effect of refractive error and accommodation on longitudinal chromatic aberration of the human eye.  Optom Vis Sci. 1993;  70 716-722
  • 2 Belin M W, Ratliff C D. Evaluating data acquisition and smoothing functions of currently available videokeratoscopes.  J Cataract Refract Surg. 1996;  22 421-426
  • 3 Bronstein I N, Semendjajew K A. Algebra. In: Grosche G, Ziegler V, Ziegler D (Hrsg) Taschenbuch der Mathematik. Frankfurt am Main; Harri Deutsch 1985
  • 4 Camp J J, Maguire L J, Cameron B M, Robb R A. A computer model for the evaluation of the effect of corneal topography on optical performance.  Am J Ophthalmol. 1990;  109 379-386
  • 5 Campbell F W, Green D G. Optical and retinal factors affecting visual resolution.  J Physiol. 1965;  181 576-593
  • 6 Carney L G. Visual loss in keratoconus.  Arch Ophthalmol. 1982;  100 1282-1285
  • 7 Clement R A, Dunne M C, Barnes D A. A method for raytracing through schematic eyes with off-axis components.  Ophthalmic Physiol Opt. 1987;  7 149-152
  • 8 Evans C J, Parks R E, Sullivan P, Taylor J. Visualization of surface figure using Zernike polynomials.  Applied Optics. 1995;  34 7815-7823
  • 9 Ginsburg A P. Spatial filtering and vision: Implications for normal and abnormal vision. In: Proenza LM, Enoch JM, Jampolsky A (eds) Clinical applications of visual psychophysics. Cambridge; Cambridge University Press 1981: 70-106
  • 10 Goodman J W. Introduction to fourier optics. New York; McGraw-Hill 1968: 101-140
  • 11 Gormley D J, Gersten M, Koplin R S, Lubkin V. Corneal modeling.  Cornea. 1988;  7 30
  • 12 Greivenkamp J E, Schwiegerling J, Miller J M, Mellinger M D. Visual acuity modeling using optical raytracing of schematic eyes.  Am J Ophthalmol. 1995;  120 227-240
  • 13 Gullstrand A. Anhang zu Teil 1. In: von Helmholtz H Physiologische Optik, 3. Ausgabe, Band 1. Hamburg; Voss 1909: 350-358
  • 14 von Helmholtz H. Die Dioptrik des Auges: Zerstreuungsbilder auf der Netzhaut (Nachtrag zum ersten Abschnitt). In: Physiologische Optik, 3. Ausgabe, Band 1. Hamburg; Voss 1909: 91-121
  • 15 Kooijman A C. Light distribution on the retina of a wide angle theoretical eye.  J Opt Soc Am. 1983;  73 1544-1550
  • 16 Langenbucher A, Seitz B, Kus M M, Steffen P, Naumann G OH. Fourieranalyse als mathematisches Modell zur Auswertung und Darstellung von postoperativen Hornhauttopographiedaten nach nichtmechanischer perforierender Keratoplastik.  Klin Monatsbl Augenheilkd . 1997;  210 197-206
  • 17 Langenbucher A, Seitz B, Neumann J, Kus M M, Naumann G OH. Modeling of corneal topography surfaces in keratoconic eyes using a subdivision scheme.  Ophthalmic Res. 1999;  31 (Suppl 1) S74
  • 18 Langenbucher A, Seitz B, Kus M M, Naumann G OH. Zernike representation of corneal topography height data after nonmechanical penetrating keratoplasty.  Invest Ophthalmol Vis Sci. 1999;  40 582-591
  • 19 Le Grand Y, El Hage S G. Optics of the eye. In: Le Grand Y, El Hage SG Physiological optics. Berlin, Heidelberg, New York; Springer 1980: 57-69
  • 20 Liang J, Williams D R, Miller D T. Supernormal vision and high-resolution retinal images through adaptive optics.  J Opt Soc Am A. 1997;  14 (11) 2884-2892
  • 21 Lotmar W. Theoretical eye model with aspherics.  J Opt Soc Am. 1971;  61 1522-1529
  • 22 Maloney R K, Bogan S J, Waring G O. Determination of image-forming properties from corneal topography.  Am J Ophthalmol. 1992;  115 31-41
  • 23 Mannis M J, Zadnik K, Johnson C A. The effect of penetrating keratoplasty on contrast sensitivity in keratoconus.  Arch Ophthalmol. 1984;  102 1513
  • 24 Marmor M F. Contrast sensitivity versus visual acuity in retinal disease.  Br J Ophthalmol. 1986;  70 553-555
  • 25 Marmor M F, Gawande A. Effect of visual blur on contrast sensitivity: Clinical implications.  Ophthalmology. 1988;  95 139-142
  • 26 van Meeteren A. Calculation of the optical modulation transfer function of the human eye for white light.  Optica Acta. 1974;  21 395-412
  • 27 Mierdel P, Kämmerer M, Krinke H E, Seiler T. Effects of photorefractive keratectomy and cataract surgery on ocular optical errors of higher order.  Graefe's Arch Clin Exp Ophthalmol. 1999;  237 725-729
  • 28 Mierdel P, Krinke H E, Wiegand W, Kämmerer M, Seiler T. Measuring device for determining monochromatic aberration of the human eye.  Ophthalmologe. 1997;  94 441-445
  • 29 Miller D. Optics and refraction. New York; Gower 1991
  • 30 Naumann G OH, Sautter H. Surgical procedures of the cornea. In: Blodi FC, Mackensen G, Neubauer H (eds) Surgical Ophthalmology 1. Berlin; Springer 1991
  • 31 Naumann G OH. Corneal transplantation in anterior segment diseases. The Bowman Lecture (Number 56) Part II.  Eye. 1995;  9 395-424
  • 32 Navarro R, Santamaria J, Bescos J. Accommodation-dependent model of the human eye with aspherics.  J Opt Soc Am. 1985;  2 1273-1281
  • 33 Roberts C. The accuracy of „power” maps to display curvature data in corneal topography systems.  Invest Ophthalmol Vis Sci. 1994;  35 3525-3532
  • 34 Roberts C. Characterization of the inherent error in a spherically-biased corneal topography system in mapping radially aspheric surface.  J Refract Corneal Surg. 1994;  10 103-116
  • 35 Sami G, El Hage S G, Berny F. Contribution of the crystalline lens to the spherical aberration of the eye.  J Opt Soc Am. 1973;  63 205-211
  • 36 Seiler T, Genth U, Holschbach A, Derse M. Aspheric photorefractive keratectomy with excimer laser.  Refract Corneal Surg. 1993;  9 166-172
  • 37 Seiler T, Reckmann W, Maloney R K. Effective spherical aberration of the cornea as a quantitative descriptor in corneal topography.  J Cataract Refract Surg. 1993;  19 155-165
  • 38 Seitz B, Behrens A, Langenbucher A. Corneal topography.  Curr Opinions Ophthalmol. 1997;  8 (IV) 8-24
  • 39 Smith W J. Modern optical engineering, 2nd ed. New York; McGraw-Hill 1990: 327-363
  • 40 Stiles W S, Crawford B H. The luminous efficiency of rays entering the eye pupil at different points.  Proc R Soc Lond (Biol). 1933;  112 428-450
  • 41 Wang J W, Rice D A, Klyce S D. Analysis of effects of astigmatism and misalignment on corneal surface reconstruction from photokeratoscopic data.  Refract Corneal Surg. 1991;  7 129-140
  • 42 Zadnik K, Mannis M J, Johnson C A. An analysis of contrast sensitivity in identical twins with keratoconus.  Cornea. 1984;  3 99-103

Priv.-Doz. Dr. Achim Langenbucher,Dipl.-Ing. 

Augenklinik mit Poliklinik der Universität Erlangen-Nürnberg

Schwabachanlage 6

91054 Erlangen

Email: achim.langenbucher@augen.imed.uni-erlangen.de

    >