Scaling Zernike expansion coefficients to different pupil sizes

Recent developments in technologies to correct aberrations in the eye have fostered extensive research in wave-front sensing of the eye, resulting in many reports of Zernike expansions of wave-front errors of the eye. For different reports of Zernike expansions, to be compared, the same pupil diameter is required. Since no standard pupil size has been established for reporting these results, a technique for converting Zernike expansion coefficients from one pupil size to another is needed. This investigation derives relationships between the Zernike expansion coefficients for two different pupil sizes.     

Statistical variation of aberration structure and image quality in a normal population of healthy eyes

A Shack-Hartmann aberrometer was used to measure the monochromatic aberration structure along the primary line of sight of 200 cyclopleged, normal, healthy eyes from 100 individuals. Sphero-cylindrical refractive errors were corrected with ophthalmic spectacle lenses based on the results of a subjective refraction performed immediately prior to experimentation. Zernike expansions of the experimental wave-front aberration functions were used to determine aberration coefficients for a series of pupil diameters. The residual Zernike coefficients for defocus were not zero but varied systematically with pupil diameter and with the Zernike coefficient for spherical aberration in a way that maximizes visual acuity. We infer from these results that subjective best focus occurs when the area of the central, aberration-free region of the pupil is maximized. We found that the population averages of Zernike coefficients were nearly zero for all of the higher-order modes except spherical aberration. This result indicates that a hypothetical average eye representing the central tendency of the population is nearly free of aberrations, suggesting the possible influence of an emmetropization process or evolutionary pressure. However, for any individual eye the aberration coefficients were rarely zero for any Zernike mode. To first approximation, wave-front error fell exponentially with Zernike order and increased linearly with pupil area. On average, the total wave-front variance produced by higher-order aberrations was less than the wave-front variance of residual defocus and astigmatism. For example, the average amount of higher-order aberrations present for a 7.5-mm pupil was equivalent to the wave-front error produced by less than 1/4 diopter (D) of defocus. The largest pupil for which an eye may be considered diffraction-limited was 1.22 mm on average. Correlation of aberrations from the left and right eyes indicated the presence of significant bilateral symmetry. No evidence was found of a universal anatomical feature responsible for third-order optical aberrations. Using the Marechal criterion, we conclude that correction of the 12 largest principal components, or 14 largest Zernike modes, would be required to achieve diffraction-limited performance on average for a 6-mm pupil. Different methods of computing population averages provided upper and lower limits to the mean optical transfer function and mean point-spread function for our population of eyes.     

Orthonormal polynomials in wavefront analysis: error analysis

Zernike circle polynomials are in widespread use for wavefront analysis because of their orthogonality over a circular pupil and their representation of balanced classical aberrations. However, they are not appropriate for noncircular pupils, such as annular, hexagonal, elliptical, rectangular, and square pupils, due to their lack of orthogonality over such pupils. We emphasize the use of orthonormal polynomials for such pupils, but we show how to obtain the Zernike coefficients correctly. We illustrate that the wavefront fitting with a set of orthonormal polynomials is identical to the fitting with a corresponding set of Zernike polynomials. This is a consequence of the fact that each orthonormal polynomial is a linear combination of the Zernike polynomials. However, since the Zernike polynomials do not represent balanced aberrations for a noncircular pupil, the Zernike coefficients lack the physical significance that the orthonormal coefficients provide. We also analyze the error that arises if Zernike polynomials are used for noncircular pupils by treating them as circular pupils and illustrate it with numerical examples.     

Orthonormal polynomials in wavefront analysis: analytical solution

Zernike circle polynomials are in widespread use for wavefront analysis because of their orthogonality over a circular pupil and their representation of balanced classical aberrations. In recent papers, we derived closed-form polynomials that are orthonormal over a hexagonal pupil, such as the hexagonal segments of a large mirror. We extend our work to elliptical, rectangular, and square pupils. Using the circle polynomials as the basis functions for their orthogonalization over such pupils, we derive closed-form polynomials that are orthonormal over them. These polynomials are unique in that they are not only orthogonal across such pupils, but also represent balanced classical aberrations, just as the Zernike circle polynomials are unique in these respects for circular pupils. The polynomials are given in terms of the circle polynomials as well as in polar and Cartesian coordinates. Relationships between the orthonormal coefficients and the corresponding Zernike coefficients for a given pupil are also obtained. The orthonormal polynomials for a one-dimensional slit pupil are obtained as a limiting case of a rectangular pupil.     

矩形孔径上折射率均匀性数据拟合研究

泽尼克圆多项式在圆形光瞳的正交性和能够代表经典像差而被广泛应用到波前分析中,用泽尼克圆多项式作为矩形光瞳基底函数,通过推导得到在矩形光瞳上正交的多项式.这个在矩形光瞳上正交的多项式不仅是唯一的,而且也能够表示经典像差,就像泽尼克圆多项式在表示圆形光瞳时具有这样的特性一样.矩形光瞳上正交多项式像泽尼克圆多项式一样即可以用...     

Hartmann-Shack technique and refraction across the horizontal visual field

We compared refractions across the horizontal visual field, based on different analyses of wave aberration obtained with a Hartmann-Shack instrument. The wave aberrations had been determined for 6-mm-diameter pupils up to at least the sixth Zernike order in five normal subjects [J. Opt. Soc. Am. A 19, 2180 (2002)]. The polynomials were converted into refractions based on 6-mm pupils and second-order Zernike aberrations (6 mm/2nd order), 3-mm pupils and second-order aberrations (3 mm/2nd order), 1-mm pupils and second-order aberrations (1 mm/2nd order), and 6-mm pupils with both second- and fourth-order aberrations (6 mm/4th order). The 3-mm/2nd-order and 6-mm/2nd-order refractions differed by as much as 0.9 D in mean sphere on axis, but the differences reduced markedly toward the edges of the visual field. The cylindrical differences between these two analyses were small at the center of the visual field (<0.3 D) but increased into the periphery to be greater than 1.0 D for some subjects. Much smaller differences in mean sphere and cylinder were found when 3-mm/2nd-order refractions and either the 1-mm/2nd-order refractions or the 6-mm/4th-order refractions were compared. The results suggest that, for determining refractions based on wave aberration data with large pupils, similar results occur by either restricting the analysis to second-order Zernike aberrations with a smaller pupil such as 3 mm or using both second- and fourth-order Zernike aberrations. Since subjective refraction is largely independent of the pupil size under photopic conditions, objective refractions based on either of these analyses may be the most useful.     

Zernike representation and Strehl ratio of optical systems with variable numerical aperture

We consider optical systems with variable numerical aperture (NA) on the level of the Zernike coefficients of the correspondingly scalable pupil function. We thus present formulas for the Zernike coefficients and their first two derivatives as a function of the scaling factor ε ≤ 1, and we apply this to the Strehl ratio and its derivatives of NA-reduced optical systems. The formulas for the Zernike coefficients of NA-reduced optical systems are also useful for the forward calculation of point-spread functions and aberration retrieval within the Extended Nijboer–Zernike (ENZ) formalism for optical systems with reduced NA or systems that have a central obstruction. Thus, we retrieve a Gaussian, comatic pupil function on an annular set from the intensity point-spread function in the focal region under high-NA conditions.     

Zernike polynomials and optical aberrations

The use of Zernike polynomials to calculate the standard deviation of a primary aberration across a circular, annular, or a Gaussian pupil is described. The standard deviation of secondary aberrations is also discussed briefly.     

General method to derive the relationship between two sets of Zernike coefficients corresponding to different aperture sizes

Zernike polynomials have been widely used to describe the aberrations in wavefront sensing of the eye. The Zernike coefficients are often computed under different aperture sizes. For the sake of comparison, the same aperture diameter is required. Since no standard aperture size is available for reporting the results, it is important to develop a technique for converting the Zernike coefficients obtained from one aperture size to another size. By investigating the properties of Zernike polynomials, we propose a general method for establishing the relationship between two sets of Zernike coefficients computed with different aperture sizes.     

高级像差对人眼成像质量和视觉的影响

定量地研究高级像差对人眼成像质量和视觉的影响对人眼像差矫正具有重要的实验和临床意义.利用Hartmann-Shack波前传感器人眼像差仪测量了正常人眼6mm瞳孔的波前像差,由波前像差计算出人眼光学系统的光学调制传递函数MTF和Strehl比率,并由MTF和视网膜空间像调制度AIM曲线计算出人眼视锐度和对比敏感度函数CSF.根据MTF和Strehl比率分析了高级像差对人眼成像质量的影响,根据视锐度和对比敏感度函数CSF分析了高级波像差对视觉的影响.研究表明Zernik前6级像差对人眼成像质量和视觉的影响是不可忽略的,更高级的像差对人眼成像质量和视觉的影响较小,甚至可以忽略.对Zemik前6级像差进行矫正,可以得到相当好的视觉.     

Extended Nijboer-Zernike approach to aberration and birefringence retrieval in a high-numerical-aperture optical system

The judgment of the imaging quality of an optical system can be carried out by examining its through-focus intensity distribution. It has been shown in a previous paper that a scalar-wave analysis of the imaging process according to the extended Nijboer-Zernike theory allows the retrieval of the complex pupil function of the imaging system, including aberrations as well as transmission variations. However, the applicability of the scalar analysis is limited to systems with a numerical aperture (NA) value of the order of 0.60 or less; beyond these values polarization effects become significant. In this scalar retrieval method, the complex pupil function is represented by means of the coefficients of its expansion in a series involving the Zernike polynomials. This representation is highly efficient, in terms of number and magnitude of the required coefficients, and lends itself quite well to matching procedures in the focal region. This distinguishes the method from the retrieval schemes in the literature, which are normally not based on Zernike-type expansions, and rather rely on point-by-point matching procedures. In a previous paper [J. Opt. Soc. Am. A 20, 2281 (2003)] we have incorporated the extended Nijboer-Zernike approach into the Ignatowsky-Richards/Wolf formalism for the vectorial treatment of optical systems with high NA. In the present paper we further develop this approach by defining an appropriate set of functions that describe the energy density distribution in the focal region. Using this more refined analysis, we establish the set of equations that allow the retrieval of aberrations and birefringence from the intensity point-spread function in the focal volume for high-NA systems. It is shown that one needs four analyses of the intensity distribution in the image volume with different states of polarization in the entrance pupil. Only in this way will it be possible to retrieve the "vectorial" pupil function that includes the effects of birefringence induced by the imaging system. A first numerical test example is presented that illustrates the importance of using the vectorial approach and the correct NA value in the aberration retrieval scheme.     

Direct transformation of Zernike eye aberration coefficients between scaled, rotated, and/or displaced pupils

In eye aberrometry it is often necessary to transform the aberration coefficients in order to express them in a scaled, rotated, and/or displaced pupil. This is usually done by applying to the original coefficients vector a set of matrices accounting for each elementary transformation. We describe an equivalent algebraic approach that allows us to perform this conversion in a single step and in a straightforward way. This approach can be applied to any particular definition, normalization, and ordering of the Zernike polynomials, and can handle a wide range of pupil transformations, including, but not restricted to, anisotropic scalings. It may also be used to transform the aberration coefficients between different polynomial basis sets.     

Zernike annular polynomials and atmospheric turbulence

Imaging through atmospheric turbulence by systems with annular pupils is discussed using the Zernike annular polynomials. Fourier transforms of these polynomials are derived analytically to facilitate the calculation of variance and covariance of the aberration coefficients. Zernike annular shape functions are derived and used to calculate the Strehl ratio and the residual phase structure and mutual coherence functions when a certain number of modes are corrected using, say, a deformable mirror. Special cases of long- and short-exposure images are also considered. The results for systems with a circular pupil are obtained as a special case of the annular pupil.     

Zernike coefficients of a scaled pupil

By expressing a scaled Zernike radial polynomial as a linear combination of the unscaled radial polynomials, we give a simple derivation for determining the Zernike coefficients of an aberration function of a scaled pupil in terms of their values for a corresponding unscaled pupil.     

Numerical parametric lens for shifting, magnification, and complete aberration compensation in digital holographic microscopy

The concept of numerical parametric lenses (NPL) is introduced to achieve wavefront reconstruction in digital holography. It is shown that operations usually performed by optical components and described in ray geometrical optics, such as image shifting, magnification, and especially complete aberration compensation (phase aberrations and image distortion), can be mimicked by numerical computation of a NPL. Furthermore, we demonstrate that automatic one-dimensional or two-dimensional fitting procedures allow adjustment of the NPL parameters as expressed in terms of standard or Zernike polynomial coefficients. These coefficients can provide a quantitative evaluation of the aberrations generated by the specimen. Demonstration is given of the reconstruction of the topology of a microlens.     

Dynamics of the eye's wave aberration

It is well known that the eye's optics exhibit temporal instability in the form of microfluctuations in focus; however, almost nothing is known of the temporal properties of the eye's other aberrations. We constructed a real-time Hartmann-Shack (HS) wave-front sensor to measure these dynamics at frequencies as high as 60 Hz. To reduce spatial inhomogeneities in the short-exposure HS images, we used a low-coherence source and a scanning system. HS images were collected on three normal subjects with natural and paralyzed accommodation. Average temporal power spectra were computed for the wave-front rms, the Seidel aberrations, and each of 32 Zernike coefficients. The results indicate the presence of fluctuations in all of the eye's aberration, not just defocus. Fluctuations in higher-order aberrations share similar spectra and bandwidths both within and between subjects, dropping at a rate of approximately 4 dB per octave in temporal frequency. The spectrum shape for higher-order aberrations is generally different from that for microfluctuations of accommodation. The origin of these measured fluctuations is not known, and both corneal/lenticular and retinal causes are considered. Under the assumption that they are purely corneal or lenticular, calculations suggest that a perfect adaptive optics system with a closed-loop bandwidth of 1-2 Hz could correct these aberrations well enough to achieve diffraction-limited imaging over a dilated pupil.     

Systematic comparison of the use of annular and Zernike circle polynomials for annular wavefronts

The theory of wavefront analysis of a noncircular wavefront is given and applied for a systematic comparison of the use of annular and Zernike circle polynomials for the analysis of an annular wavefront. It is shown that, unlike the annular coefficients, the circle coefficients generally change as the number of polynomials used in the expansion changes. Although the wavefront fit with a certain number of circle polynomials is identically the same as that with the corresponding annular polynomials, the piston circle coefficient does not represent the mean value of the aberration function, and the sum of the squares of the other coefficients does not yield its variance. The interferometer setting errors of tip, tilt, and defocus from a four-circle-polynomial expansion are the same as those from the annular-polynomial expansion. However, if these errors are obtained from, say, an 11-circle-polynomial expansion, and are removed from the aberration function, wrong polishing will result by zeroing out the residual aberration function. If the common practice of defining the center of an interferogram and drawing a circle around it is followed, then the circle coefficients of a noncircular interferogram do not yield a correct representation of the aberration function. Moreover, in this case, some of the higher-order coefficients of aberrations that are nonexistent in the aberration function are also nonzero. Finally, the circle coefficients, however obtained, do not represent coefficients of the balanced aberrations for an annular pupil. The various results are illustrated analytically and numerically by considering an annular Seidel aberration function.     

Measuring ocular aberrations in the peripheral visual field using Hartmann-Shack aberrometry

We have used the Hartmann-Shack technique previously to measure ocular aberrations along the horizontal meridian of the visual field. This requires considerable modifications from the technique for measuring the aberrations corresponding to the fovea. We now further develop the technique so that it can be used for any meridian of the visual field. Allowance is made for any auxiliary optics placed in front of the eye to compensate for the limited range of the Hartmann-Shack technique and for the case where aberrations are estimated at a wavelength other than the measuring wavelength. Zernike wave aberrations are converted to peripheral refractions. Examples are presented showing the developments, and we discuss change in wave aberrations when converting from a circular to an elliptical pupil.