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.     

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.     

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

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

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.     

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.     

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.     

Ray-tracing and Aberration Formulae for a General Optical System

Abstract

All aspects of the ray-tracing and of the calculation of the monochromatic and chromatic aberrations for a completely general optical system are treated. These include the specification of the system, the ray-tracing formulae for refraction, reflection and transfer, the introduction of pseudo-paraxial variables in the object and image spaces, the determination of the pupil domain and the best image plane, plus calculation of aberrations along rays. Hopkins' canonical coordinates are employed in the object and image space to facilitate the calculation of the aberration polynomial coefficients, for use in the evaluation of diffraction-based image criteria, such as the point spread function. All the formulae involved are always determinate and of good accuracy, no matter whether the object, image or either pupil, of the system is at finite distance or at infinity.     

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.     

Effect of a Low Efficient Shutter on Balanced Aberrations

The effect of a 50 per cent efficient between-the-lens shutter which may be considered as one of the typical low efficient shutters is studied from the viewpoint of balanced aberrations. A low efficient between-the-lens shutter which operates circularly produces an effect similar to a shaded aperture in which the amplitude transmittance decreases from the centre of the aperture to its rim. The results of balanced aberrations calculated from the Strehl definition of an exposure image of the point source are compared with both the clear aperture and the shaded aperture. It is shown that the balanced aberrations are affected by a 50 per cent efficient shutter. When the aberration has only primary spherical aberration or primary coma, the shift of best focal setting caused by the aberration is reduced by using a 50 per cent efficient shutter.     

Maréchal Intensity Criteria for Apertures with Polynomial Non-uniform Transmission. Evaluation of the Diffraction Focus and the Strehl Ratio

Abstract

General formulae for the wave aberration weighted variance are obtained for rotationally symmetric systems with a non-uniform transmission pupil expressed by polynomials. They are valid for any combination of residual aberrations and any polynomial transmission function on the pupil. Expressions to obtain the position of the diffraction focus and the Strehl ratio are given. The accuracy in the evaluation of the diffraction focus is studied for optical systems with different residual aberrations.     

复合静电-磁浸没透镜的宽束曲轴聚焦性质及像差理论

研究了新型的复合静电 磁浸没透镜的曲光轴高斯电子光学性质及像差特性。应用宽束曲轴理论 ,导出了三维局部正交坐标系中的电子运动的中心轨迹方程和曲光轴的近轴轨迹方程 ;利用数学软件Mathematica推导出了复合静电 磁浸没透镜全部曲轴二级像差系数。作为实例 ,文中针对轴上磁场和电场具有某种解析表达式的旋转对称的静电 磁浸没透镜系数 ,计算了它的曲轴高斯聚焦特性和二级像差 ,并给出了二级像差分布图形     

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.     

On the computation of the nijboer-zernike aberration integrals at arbitrary defocus

Abstract

We present a new computation scheme for the integral expressions describing the contributions of single aberrations to the diffraction integral in the context of an extended Nijboer-Zernike approach. Such a scheme, in the form of a power series involving the defocus parameter with coefficients given explicitly in terms of Bessel functions and binomial coefficients, was presented recently by the authors with satisfactory results for small-to-medium-large defocus values. The new scheme amounts to systemizing the procedure proposed by Nijboer in which the appropriate linearization of products of Zernike polynomials is achieved by using certain results of the modern theory of orthogonal polynomials. It can be used to compute point-spread functions of general optical systems in the presence of arbitrary lens transmission and lens aberration functions and the scheme provides accurate data for any, small or large, defocus value and at any spatial point in one and the same format. The cases with high numerical aperture, requiring a vectorial approach, are equally well handled. The resulting infinite series expressions for these point-spread functions, involving products of Bessel functions, can be shown to be practically immune to loss of digits. In this respect, because of its virtually unlimited defocus range, the scheme is particularly valuable in replacing numerical Fourier transform methods when the defocused pupil functions require intolerably high sampling densities.     

Effect of sampling on real ocular aberration measurements

The minimum number of samples necessary to fully characterize the aberration pattern of the eye is a question under debate in the clinical as well as the scientific community. We performed repeated measurements of ocular aberrations in 12 healthy nonsurgical human eyes and in 3 artificial eyes, using different sampling patterns (hexagonal, circular, and rectangular with 19 to 177 samples, and 3 radial patterns with 49 sample coordinates corresponding to zeros of the Albrecht, Jacobi, and Legendre functions). For each measurement set we computed two different metrics based on the root-mean-square (RMS) of difference maps (RMS_Diff) and the proportional change in the wavefront (W%). These metrics are used to compare wavefront estimates as well as to summarize results across eyes. We used computer simulations to extend our results to "abnormal eyes" (keratoconic, post-LASIK, and post-radial keratotomy eyes). We found that the spatial distribution of the samples can be more important than the number of samples for both our measured as well as our simulated "abnormal" eyes. Experimentally, we did not find large differences across patterns except, as expected, for undersampled patterns.     

Ocular aberrations with ray tracing and Shack-Hartmann wave-front sensors: does polarization play a role?

Ocular aberrations were measured in 71 eyes by using two reflectometric aberrometers, employing laser ray tracing (LRT) (60 eyes) and a Shack-Hartmann wave-front sensor (S-H) (11 eyes). In both techniques a point source is imaged on the retina (through different pupil positions in the LRT or a single position in the S-H). The aberrations are estimated by measuring the deviations of the retinal spot from the reference as the pupil is sampled (in LRT) or the deviations of a wave front as it emerges from the eye by means of a lenslet array (in the S-H). In this paper we studied the effect of different polarization configurations in the aberration measurements, including linearly polarized light and circularly polarized light in the illuminating channel and sampling light in the crossed or parallel orientations. In addition, completely depolarized light in the imaging channel was obtained from retinal lipofuscin autofluorescence. The intensity distribution of the retinal spots as a function of entry (for LRT) or exit pupil (for S-H) depends on the polarization configuration. These intensity patterns show bright corners and a dark area at the pupil center for crossed polarization, an approximately Gaussian distribution for parallel polarization and a homogeneous distribution for the autofluorescence case. However, the measured aberrations are independent of the polarization states. These results indicate that the differences in retardation across the pupil imposed by corneal birefringence do not produce significant phase delays compared with those produced by aberrations, at least within the accuracy of these techniques. In addition, differences in the recorded aerial images due to changes in polarization do not affect the aberration measurements in these reflectometric aberrometers.     

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.     

Influences of reference plane and direction of measurement on eye aberration measurement

We explored effects of measurement conditions on wave aberration estimates for uncorrected, axially myopic model eyes. Wave aberrations were initially referenced to either the anterior corneal pole or the natural entrance pupil of symmetrical eye models, with rays traced into the eye from infinity (into the eye) to simulate normal vision, into the eye from infinity and then back out of the eye from the retinal intercepts (into/out of the eye), or out of the eye from the retinal fovea (out of the eye). The into-the-eye and out-of-the-eye ray traces gave increases in spherical aberration as myopia increased, but the into/out-of-the-eye ray trace showed little variation in spherical aberration. Reference plane choice also affected spherical aberration. Corresponding residual aberrations were calculated after the models had been optically corrected, either by placing the object or image plane at the paraxial far point or by modifying corneas to simulate laser ablation corrections. Correcting aberrations by ablation was more complete if the original aberrations were referenced to the cornea rather than to the entrance pupil. For eyes corrected by spectacle lenses, failure to allow for effects of pupil magnification on apparent entrance pupil diameter produced larger changes in measured aberrations. The general findings regarding choice of reference plane and direction of measurement were found to be equally applicable to eyes that lacked rotational symmetry.