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.     

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

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

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.     

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.     

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.     

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 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.     

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.     

Comparison of Hartmann analysis methods

Analysis of Hartmann-Shack wavefront sensors for the eye is traditionally performed by locating and centroiding the sensor spots. These centroids provide the gradient, which is integrated to yield the ocular aberration. Fourier methods can replace the centroid stage, and Fourier integration can replace the direct integration. The two--demodulation and integration--can be combined to directly retrieve the wavefront, all in the Fourier domain. Now we applied this full Fourier analysis to circular apertures and real images. We performed a comparison between it and previous methods of convolution, interpolation, and Fourier demodulation. We also compared it with a centroid method, which yields the Zernike coefficients of the wavefront. The best performance was achieved for ocular pupils with a small boundary slope or far from the boundary and acceptable results for images missing part of the pupil. The other Fourier analysis methods had much higher tolerance to noncentrosymmetric apertures.     

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.     

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.     

Mathematical construction and perturbation analysis of Zernike discrete orthogonal points

Zernike functions are orthogonal within the unit circle, but they are not over the discrete points such as CCD arrays or finite element grids. This will result in reconstruction errors for loss of orthogonality. By using roots of Legendre polynomials, a set of points within the unit circle can be constructed so that Zernike functions over the set are discretely orthogonal. Besides that, the location tolerances of the points are studied by perturbation analysis, and the requirements of the positioning precision are not very strict. Computer simulations show that this approach provides a very accurate wavefront reconstruction with the proposed sampling set.     

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.     

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.     

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.     

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.     

Three-dimensional relationship between high-order root-mean-square wavefront error, pupil diameter, and aging

We report root-mean-square (RMS) wavefront error (WFE) for individual aberrations and cumulative high-order (HO) RMS WFE for the normal human eye as a function of age by decade and pupil diameter in 1 mm steps from 3 to 7 mm and determine the relationship among HO RMS WFE, mean age for each decade of life, and luminance for physiologic pupil diameters. Subjects included 146 healthy individuals from 20 to 80 years of age. Ocular aberration was measured on the preferred eye of each subject (for a total of 146 eyes through dilated pupils; computed for 3, 4, 5, 6, and 7 mm pupils; and described with a tenth-radial-order normalized Zernike expansion. We found that HO RMS WFE increases faster with increasing pupil diameter for any given age and pupil diameter than it does with increasing age alone. A planar function accounts for 99% of the variance in the 3-D space defined by mean log HO RMS WFE, mean age for each decade of life, and pupil diameter. When physiologic pupil diameters are used to estimate HO RMS WFE as a function of luminance and age, at low luminance (9 cd/m2) HO RMS WFE decreases with increasing age. This normative data set details (1) the 3-D relationship between HO RMS WFE and age for fixed pupil diameters and (2) the 3-D relationship among HO RMS WFE, age, and luminance for physiologic pupil diameters.     

Holography-based wavefront sensing

We describe a modal wavefront sensing technique of using multiplexed holographic optical elements (HOEs). The phase pattern of a set of aberrations is angle multiplexed in a HOE, and the correlated information is obtained with a position sensing detector. The recorded aberration pattern is based on an orthogonal basis set, the Zernike polynomials, and a spherical reference wave. We show that only two recorded holographic patterns for any particular aberration type are sufficient to allow interpolated readout of aberrations to lambda/50. In this paper, we demonstrate the capability of detecting errors between +/-2lambda PV for each orthogonal set at rates limited only by the speeds of the detection electronics, which could be up to 1 MHz. We show how we take advantage of the unavoidable intermodal and intramodal cross talks in determining the type, amplitude, and orientation of the wavefront aberrations.     

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 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.