Using matrices to understand geometrical optics
José B . Almeida2000, SPIE Proceedings
Sign up for access to the world's latest research
checkGet notified about relevant papers
checkSave papers to use in your research
checkJoin the discussion with peers
checkTrack your impact
Abstract
The optimization of an optical system benefits greatly from a study of its aberrations and an identification of each of its elements' contribution to the overall aberration figures. The matrix formalism developed by one of the authors was the object of a previous paper and allows the expression of image-space coordinates as high-order polynomials of object-space coordinates. In this paper we approach the question of aberrations, both through the evaluation of the wavefront evolution along the system and its departure from the ideal spherical shape and the use of ray density plots. Using seventh-order matrix modeling, we can calculate the optical path between any two points of a ray as it travels along the optical system and we define the wavefront as the locus of the points with any given optical path; the results are presented on the form of traces of the wavefront on the tangential plane, although the formalism would also permit sagital plane plots. Ray density plots are obtained by actual derivation of the seventh-order polynomials.
Related papers
Wavefront and ray-density plots using seventh-order matrices
Optik - International Journal for Light and Electron Optics, 2005
The optimization of an optical system benefits greatly from a study of its aberrations and an identification of each of its elements' contribution to the overall aberration figures. The matrix formalism developed by the author was the object of a previous paper and allows the expression of image-space coordinates as high-order polynomials of objectspace coordinates. In this paper we approach the question of aberrations, both through the evaluation of the wavefront evolution along the system and its departure from the ideal spherical shape and the use of ray density plots. Using seventh-order matrix modelling, we can calculate the optical path between any two points of a ray as it travels along the optical system and we define the wavefront as the locus of the points with common optical path length; the results are presented on the form of traces of the wavefront on the tangential plane, although the formalism would also permit sagittal plane plots. Ray density plots are obtained by actual derivation of the seventh-order polynomials.
Analytical Study of Optical Wavefront Aberrations Using Maple
1997
This paper describes a package for analytical ray tracing of relatively simple optical systems. AESOP (An Extensible Symbolic Optics Package) enables analysis of the effects of small optical element misalignments or other perturbations. (It is possible to include two or more simultaneous independent perturbations.) Wavefront aberrations and optical path variations can be studied as functions of the perturbation parameters. The power of this approach lies in the fact that the results can be manipulated algebraically, allowing determination of misalignment tolerances as well as developing physical intuition, especially in the picometer regime of optical path length variations.
Programming matrix optics into Mathematica
Optik - International Journal for Light and Electron Optics, 2005
The various non-linear transformations incurred by the rays in an optical system can be modelled by matrix products up to any desired order of approximation. Mathematica software has been used to find the appropriate matrix coefficients for the straight path transformation and for the transformations induced by conical surfaces, both direction change and position offset. The same software package was programmed to model optical systems in seventh-order. A Petzval lens was used to exemplify the modelling power of the program.
Notes on Matrix Method of Geometrical Optics - Basic ideas of Matrix Methods in Geometrical Optics: Lecture notes (17 pages)
DESCRIPTION Lecture notes (17 pages) : The matrix methods in geometrical optics can be developed in terms of three basic operators (matrices) for reflection, refraction and translation and can be applied with ease even for an optical system consisting of a large number of optical elements. An elementary but thorough discussion of this method is presented here.
Transformations of aberrations in optical systems
JOSA A, 2011
The general transformation properties of aberrations in cylindrically symmetric optical systems are described using matrix analysis. The aberrations in an optical system can change with any change in the position of the aperture stop, or when the system operates at new conjugates even if the relative positions of the refractive surfaces and their powers are fixed. Expanding the wavefront aberration function in terms of aberration coefficients allows the new aberration coefficients to be written as linear combinations of the old aberration coefficients for every order. A pattern is established by which higher-order aberration transformations can be calculated.
Visual optics under the wavefront perspective
Arquivos Brasileiros de Oftalmologia, 2014
Some intriguing concepts of visual optics cannot be explained by ray tracing. However, they can be clarified using wavefront formalism. Its main advantage is in the use of the concept of vergence, which is very helpful in interpreting the optical phenomena involved in the neutralization of the ametropias. In this line of thinking, the major role of a lens is in the creation of a new light source (the image point) that orientates the refracted waves. Once the nature and position of this source is known, one can easily predict the behavior of the wavefronts. The formalism also allows for an easier understanding on how wavefronts relate to light rays and on how algebraic signs are assigned to optical distances.
Expression for Aberration Coefficients Using Nonlinear Transforms--Extension to Spherocylindrical Surfaces
2000
Previously we have used a matrix method to do nonlinear calculations and presented exact expressions for the third and fifth order aberration coefficients for a general conic surface as functions of the latus rectum and eccentricity. In this paper we show how the method can be extended to include the case of a general spherocylindrical surface. (Optom Vis Sci 2000;77:156-162)
Calculation of the aberration spot: improved numerical algorithm
Journal of Modern Optics, 1998
In the paper a numerical algorithm for calculation of the point spread function from the complex light amplitude known on the output pupil of the imaging system is described. First step of this algorithm includes a global polynomial approximation of the input data given in a number of points distributed arbitrarily over the entire output pupil. In the next step the numerical calculation of the diffraction integral is performed. After transition to the polar co-ordinates and dividing the region of integration into a number of infinitesimal subregions of square shape the exponent in the integrand is approximated locally by a polynomial and then the whole integrand is approximated by a polynomial multiplied by the exponential function with linear exponent. Finally the respective integration is performed analytically.
Paraxial Optics of Astigmatic Systems: Relations Between the Wavefront and the Ray Picture Approaches
Optometry and Vision Science, 2007
Purpose. The paraxial propagation of astigmatic wavefronts through astigmatic optical systems can be described by the augmented step-along method (ASAM). Its equivalence to the linear ray optics approach is considered in detail. Methods. The ASAM is exploited to derive paraxial ray paths through a general coaxial astigmatic system. Results. Starting from the information inherent in the ASAM all 2 ϫ 2 submatrices rendering the general 4 ϫ 4 transference of linear optics can be generated. This proves the complete equivalence of both approaches. Additionally, we show that the symplectic relations are automatically obeyed in the ASAM. Conclusions. The ASAM offers a complete alternative to describe the paraxial optics of astigmatic optical systems. According to the ASAM, an optical system is fully characterized by the back vertex vergence and the angular magnification matrix. Hence, a complete description of the paraxial optics of an eye should not only report the state of refraction but the angular magnification matrix as well, although it is not yet very common. The magnification matrix might be important in cases of anisometropia or the design of progressive addition lenses. Yet, a simple clinical procedure to determine the angular magnification matrix is missing. (Optom Vis Sci 2007;84:E72-E78)
Characteristic functions and the aberrations of symmetric optical systems. I. Transverse aberrations when the eikonal is given
Journal of the Optical Society of America A, 1988
We describe the transverse aberrations of a symmetric optical system with a given eikonal (angle characteristic function). The aberrations are expressed in the coordinates of the image field and the pupil parameters according to Schwarzschild. Conditions for freedom from aberrations and for aberrations in the focal plane are considered. We show how the transverse aberrations up to the ninth order can be calculated from eikonal coefficients; we give the conditions under which the eikonal coefficients are identical to the aberration coefficients.
Related topics
Related papers
Comparison of matrix method and ray tracing in the study of complex optical systems
SPIE Proceedings, 2000
In the context of the classical study of optical systems within the geometrical Gauss approximation, the cardinal elements are efficiently obtained with the aid of the transfer matrix between the input and output planes of the system. In order to take into account the geometrical aberrations, a ray tracing approach, using the Snell-Descartes laws, has been implemented in an interactive software. Both methods are applied for measuring the correction to be done to a human eye suffering from ametropia. This software may be used by optometrists and ophthalmologists for solving the problems encountered when considering this pathology. The ray tracing approach gives a significant improvement and could be very helpful for a better understanding of an eventual surgical act.
Lie-Hamilton geometric aberration optics: symmetries and symbolic computation
1997
Recibido el 28 de noviembre de 1996;aceptado el 13 de febrero de 1997
Simulations of Four Types of Optical Aberrations using Zernik Polynomials
Iraqi Journal of Science, 2017
In this paper, a computer simulation is implemented to generate of an optical aberration by means of Zernike polynomials. Defocus, astigmatism, coma, and spherical Zernike aberrations were simulated in a subroutine using MATLAB function and applied as a phase error in the aperture function of an imaging system. The studying demonstrated that the Point Spread Function (PSF) and Modulation Transfer Function (MTF) have been affected by these optical aberrations. Areas under MTF for different radii of the aperture of imaging system have been computed to assess the quality and efficiency of optical imaging systems. Phase conjugation of these types aberration has been utilized in order to correct a distorted wavefront. The results showed that the largest effect on the PSF and MTF is due to the contribution of the third type coma aberrated wavefront.
Derivation of the propagation equations for higher order aberrations of local wavefronts
Journal of the Optical Society of America A, 2011
From the literature the analytical calculation of local power and astigmatism of a wavefront after refraction and propagation is well known; it is, e.g., performed by the Coddington equation for refraction and the classical vertex correction formula for propagation. Recently the authors succeeded in extending the Coddington equation to higher order aberrations (HOA). However, equivalent analytical propagation equations for HOA do not exist. Since HOA play an increasingly important role in many fields of optics, e.g., ophthalmic optics, it is the purpose of this study to extend the propagation equations of power and astigmatism to the case of HOA (e.g., coma and spherical aberration). This is achieved by local power series expansions. In summary, with the results presented here, it is now possible to calculate analytically the aberrations of a propagated wavefront directly from the aberrations of the original wavefront containing both low-order and high-order aberrations.
Aberration Matrices of an Axial Bundle and an Investigation of Their Elements
Applied Optics, 1966
Two aberration matrices are presented. The first one is the aberrations contribution matrix formulated as a difference of two Herzberger refraction matrices of an optical surface. The second one is called aberration matrix and describes the aberrations of an axial bundle in the Gaussian image plane of an optical system. The interpretation of the elements of the aberration matrix is given, and the connection of these elements with the spherical aberration and Abbe's sine condition is shown.
Calculational technique for optical systems in the region of third-order aberrations
Journal of Optical Technology, 1999
This paper proposes a computerized calculational technique for third-order aberrations of optical systems ͑including zoom systems͒ that consist of thin components. This technique makes it possible to correct all third-order monochromatic aberrations, as well as first-order chromatic aberration. In essence, it consists of solving a system of aberration inequalities. It is shown that solving a system of aberration inequalities is a problem in nonlinear programming, in particular quadratic programming. An analytical expression is presented for an estimating function that minimizes the deviation of the aberrations from given values.
ABCD matrix for calculating third-order aberrations of gradient index optical elements
2010
Third-order aberrations need the paraxial raytracing of a marginal as well as principal ray through the entire system to be determined. When the optical system has GRIN media, even the paraxial raytrace could be cumbersome. This work presents the use of the ABCD matrix for obtaining the raytracing of those rays and, thus, it could ease the calculation of the third-order aberration coefficients, i.e. contribution of GRIN to surface refraction and transfer, for each individual surface of GRIN optical elements. A procedure to be implemented in any modern commercial optical software has been presented and evaluated for simple systems.
MOTO: a Matlab Object-oriented programming Toolbox for Optics
Education and Training in Optics and Photonics, 2007
The ray optics is the branch of optics in which all the wave effects are neglected: the light is considered as travelling along rays which can only change their direction by refraction or reflection. On one hand, a further simplifying approximation can be made if attention is restricted to rays travelling close to the optical axis and at small angles: the well-known linear or paraxial approximation introduced by Gauss. On the other hand, in order to take into account the geometrical aberrations, it is sometimes necessary to pay attention to marginal rays with the aid of a ray tracing procedure. This contribution describes a toolbox for the study of optical systems which implements both approaches. It has been developed in the framework of an educational project, but it is general enough to be useful in most of the cases.
Lie algebraic treatment of dioptric power and optical aberrations
Journal of the Optical Society of America A, 1998
The dioptric power of an optical system can be expressed as a four-component dioptric power matrix. We generalize and reformulate the standard matrix approach by utilizing the methods of Lie algebra. This generalization helps one deal with nonlinear problems (such as aberrations) and further extends the standard matrix formulation. Explicit formulas giving the relationship between the incident and the emergent rays are presented. Examples include the general case of thick and thin lenses. The treatment of a graded-index medium is outlined. © 1998 Optical Society of America [S0740-3232(98)03909-X]
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.