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 aberration polynomials for anamorphic optical imaging systems with circular pupils

In a recent paper, we considered the classical aberrations of an anamorphic optical imaging system with a rectangular pupil, representing the terms of a power series expansion of its aberration function. These aberrations are inherently separable in the Cartesian coordinates (x,y) of a point on the pupil. Accordingly, there is x-defocus and x-coma, y-defocus and y-coma, and so on. We showed that the aberration polynomials orthonormal over the pupil and representing balanced aberrations for such a system are represented by the products of two Legendre polynomials, one for each of the two Cartesian coordinates of the pupil point; for example, L(l)(x)L(m)(y), where l and m are positive integers (including zero) and L(l)(x), for example, represents an orthonormal Legendre polynomial of degree l in x. The compound two-dimensional (2D) Legendre polynomials, like the classical aberrations, are thus also inherently separable in the Cartesian coordinates of the pupil point. Moreover, for every orthonormal polynomial L(l)(x)L(m)(y), there is a corresponding orthonormal polynomial L(l)(y)L(m)(x) obtained by interchanging x and y. These polynomials are different from the corresponding orthogonal polynomials for a system with rotational symmetry but a rectangular pupil. In this paper, we show that the orthonormal aberration polynomials for an anamorphic system with a circular pupil, obtained by the Gram-Schmidt orthogonalization of the 2D Legendre polynomials, are not separable in the two coordinates. Moreover, for a given polynomial in x and y, there is no corresponding polynomial obtained by interchanging x and y. For example, there are polynomials representing x-defocus, balanced x-coma, and balanced x-spherical aberration, but no corresponding y-aberration polynomials. The missing y-aberration terms are contained in other polynomials. We emphasize that the Zernike circle polynomials, although orthogonal over a circular pupil, are not suitable for an anamorphic system as they do not represent balanced aberrations for such a system.     

Zernike-gauss polynomials and optical aberrations of systems with gaussian pupils

In the first two Notes of this series,(l,2) we discussed Zernike circle and annular polynomials that represent optimally balanced classical aberrations of systems with uniform circular or annular pupils, respectively. Here we discuss Zernike-Gauss polynomials which are the corresponding polynomials for systems with Gaussian circular or annular pupils.(3-5) Such pupils, called apodized pupils, are used in optical imaging to reduce the secondary rings of the pointspread functions of uniform pupils.(6) Propagation of Gaussian laser beams also involves such pupils.     

Zernike circle polynomials and optical aberrations of systems with circular pupils

Zernike circle polynomials, their numbering scheme, and relationship to balanced optical aberrations of systems with circular pupils are discussed.     

Transformation of Zernike coefficients: scaled, translated, and rotated wavefronts with circular and elliptical pupils

Zernike polynomials and their associated coefficients are commonly used to quantify the wavefront aberrations of the eye. When the aberrations of different eyes, pupil sizes, or corrections are compared or averaged, it is important that the Zernike coefficients have been calculated for the correct size, position, orientation, and shape of the pupil. We present the first complete theory to transform Zernike coefficients analytically with regard to concentric scaling, translation of pupil center, and rotation. The transformations are described both for circular and elliptical pupils. The algorithm has been implemented in MATLAB, for which the code is given in an appendix.     

Orthonormal aberration polynomials for anamorphic optical imaging systems with rectangular pupils

The classical aberrations of an anamorphic optical imaging system, representing the terms of a power-series expansion of its aberration function, are separable in the Cartesian coordinates of a point on its pupil. We discuss the balancing of a classical aberration of a certain order with one or more such aberrations of lower order to minimize its variance across a rectangular pupil of such a system. We show that the balanced aberrations are the products of two Legendre polynomials, one for each of the two Cartesian coordinates of the pupil point. The compound Legendre polynomials are orthogonal across a rectangular pupil and, like the classical aberrations, are inherently separable in the Cartesian coordinates of the pupil point. They are different from the balanced aberrations and the corresponding orthogonal polynomials for a system with rotational symmetry but a rectangular pupil.     

Comparison of annular wavefront interpretation with Zernike circle polynomials and annular polynomials

A general wavefront fitting procedure with Zernike annular polynomials for circular and annular pupils is proposed. For interferometric data of typical annular wavefronts with smaller and larger obscuration ratios, the results fitted with Zernike annular polynomials are compared with those of Zernike circle polynomials. Data are provided demonstrating that the annular wavefront expressed with Zernike annular polynomials is more accurate and meaningful for the decomposition of aberrations, the calculation of Seidel aberrations, and the removal of misalignments in interferometry. The primary limitations of current interferogram reduction software with Zernike circle polynomials in analyzing wavefronts of annular pupils are further illustrated, and some reasonable explanations are provided. It is suggested that the use of orthogonal basis functions on the pupils of the wavefronts analyzed is more appropriate.     

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.     

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 annular polynomials and optical aberrations of systems with annular pupils

Zernike annular polynomials that represent orthogonal andbalanced aberrations suitable for systems with annular pupilsare described. Their numbering scheme is the same asfor Zernike circle polynomials. Expressions for standard deviationof primary and balanced primary aberrations are given.     

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 expansion of separable functions of cartesian coordinates

A Zernike expansion over a circle is given for an arbitrary function of a single linear spatial coordinate. The example of a half-plane mask (Hilbert filter) is considered. The expansion can also be applied to cylindrical aberrations over a circular pupil. A product of two such series can thus be used to expand an arbitrary separable function of two Cartesian coordinates.     

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.     

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.     

Zernike polynomials: a guide

In this paper we review a special set of orthonormal functions, namely Zernike polynomials which are widely used in representing the aberrations of optical systems. We give the recurrence relations, relationship to other special functions, as well as scaling and other properties of these important polynomials. Mathematica code for certain operations are given in the Appendix.     

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.     

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.     

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.