In recent years, large field of view (FOV) angle optical imaging system, such as fish-eye lens and catadioptric panoramic imaging system, have attracted more and more attention. It has been widely used in robot or vehicle navigation, positioning, scene monitoring, external space detection, micro intelligent system and so on. In daily life, the fish-eye lenses are also used by many photography enthusiasts. In April 2015, Ricoh released a 360° panoramic camera, which used two fish-eye lens systems with a field of view angle more than 180° to realize the imaging of the three-dimensional angle of the whole space, causing people’s attention all over the world.

The field of view angle of fish-eye lens system reach or even exceed 180°. The beam focus, wave array surface and aberration of the optical system in the meridian and sagittal plane are totally different, which has the imaging characteristics of the plane symmetric optical system. At the same time, people often pursue the large aperture imaging (e.g. F#<3). Therefore, this kind of system has three characteristics: plane symmetry, super large field of view, and large aperture imaging, which make the design of fish-eye lens optical system very complex. Aberration theory is an important means to study optical system. Because of the plane symmetry imaging characteristics for this kind of lens system. Seidel primary aberration and high-order aberration theories that develop from axisymmetric optical system are no longer suitable for aberration analysis and design of this super large field of view optical system.

The conventional lens design of different structures and their aberration characteristics are quite clear^[1]. From the book, fish-eye lens optics written by Wang Yongzhong^[2], it is also pointed out that there is no a law in the structure and parameter design for fish-eye lens system. Most people use patents as the initial structure, and then use some commercial software to improve and optimize the initial structure repeatedly. Due to the lack of the aberration theory guidance, the selection of the initial structure of the fish-eye lens optical system is blind. We do not know what kind of aberration characteristics it has, and whether the initial structure is reasonable enough to determine the final performance of the fish-eye lens. Secondly, the commercial software application to optimize system parameters depends on algorithm to a great extent, and its design process is essentially "black box operation". If we can use reasonably the precise wave aberration theory for the optical system with super large field of view, we can analyze the function and contribution and all kinds of aberration of the different optics lens in the fish-eye lens system, we can analysis how to design the optical system structure and parameters and achieve the ideal aberration balance.

For the optical system of large aperture imaging, the effect of higher-order aberrations becomes significant^[3]. The aberration accuracy using fourth-order wave aberration theory is not enough. We use the modulation transfer function based on the fourth-order wave aberration theory as the evaluation function to optimize the design of Nikon 16 mm/F2.8 fish-eye lens. We take the original design of the optical lens parameters as the initial value, and search in a large range for its optimization. The results show that the performance of the optimization results is better than the original design in the case of small aperture (F#<5), but worse in the case of large aperture. Therefore, only further development of the higher-order aberrations optical system theory with super large field of view, can we objectively research the imaging law of such fish-eye lens system.

It is a classical method to solve optical system aberration with wave aberration theory, and it is more suitable for the aberration analysis and calculation of multi-element optical system. For a multi-element optical system, the total system wave aberration is the sum of the wave aberrations for each elements. However, the derivation model of the wave aberration expression with a single element generally assume that the object point is an ideal geometric object point. According to the first-order optical surface parameters, the correlation ray path is expanded into a series of pupil coordinates and field angle, and then every order wave aberration expression (called intrinsic wave aberrations) is obtained. Except for the first optical plane, the other optical plane in a multi-element optical system, due to the aberration influence by its front optical system, the assumption of the ideal geometry point is not tenable in a strict sense. Therefore, the intrinsic wave aberration only reflect the single system element aberration in a real sense. The intrinsic wave aberration ahead of system will produce a correction amount of wave aberration (called derivative wave aberration), which is related to the intrinsic wave aberration of the optical plane and the optical system ahead of it. Therefore, the aberration with multi-element optical system is composed of the intrinsic aberration and the derivative aberration. For large aperture imaging optical system, the contribution of derivative aberration cannot be ignored. The evolution process of higher-order intrinsic aberration is basically the same as that of lower-order intrinsic aberration, but the evolution ideas of derivative aberrations and intrinsic aberrations are totally different. The research on derivative aberrations for the plane symmetric optical system cannot directly apply the derivative aberrations methods by axisymmetric optical system. Secondly, for the paraxial optical system aberration, the light pupil coordinates transmission on the adjacent optical surface is generally treated as in a linear way, the nonlinear term contribution is generally small and its influence on the aberration can be ignored. However, for a large field of view optical system, such as a fish-eye lens, the light beams from object with large field of view hit the optical surface. If the pupil coordinates transmission between adjacent optical surfaces is still treated as a linear way, there will be more calculation error with the increase of aperture. Through the research on many fish-eye lens systems, it is shown that the nonlinear term in the pupil coordinate transfer relationship has a significant effect on the object point wave aberration in a large field of view optical system, even more than the derivative wave aberration.

The basic structure of the fish-eye lens is composed of the former group with several large negative meniscus optical power lenses and the latter group with several common objective lenses. The former group is mainly used to compress the field angle. Most of the aberration analysis involved in the fish-eye lens research reports is still based on Seidel’s primary aberration theory. Dai Jianning, in his doctoral dissertation on fish-eye lens research^[4], the aberration of several field rays and aperture rays is calculated by ray tracing to measure the high-order spherical aberration and high-order coma approximately, and it is included in the evaluation index of optimization design. The aberration theory can play an important guiding role in the initial structure and parameter design for the fish-eye lens optical system^[5].

In some professional optical design software, the sampling optical light path difference is mostly used to simulate the wave aberration distribution at the exit pupil of optical system, and then calculate the optical system transfer function to evaluate its imaging and its accuracy, that all depend on the number of sampling light. The evaluation function based on the wave aberration theory includes only one or several aberrations contribution. In addition, it is different from the calculation method of light sampling.

The sixth-order wave aberration theory mainly include the sixth-order intrinsic wave aberration, the fifth order aberration, and the derivative wave aberration and the influence of the second order precision of pupil coordinates on the wave aberration of plane symmetric optical system.

According to the fourth-order wave aberration theory of plane symmetric optical system developed by Lu Lijun, the sixth-order accuracy expressions of meridian and sagittal focal line with respect to pupil coordinates are obtained by using the torus as the reference wave surface^[6]. Then, by using ray equation and polynomial fitting method, the mapping relationship between the light coordinates of the incident/outgoing wave array and the optical element surface is obtained. According to the wave aberration definition of the plane symmetric optical system, the intrinsic wave aberration coefficients from the first-order to the sixth-order can be obtained by expanding the optical path of the relevant ray path with respect to the pupil coordinates according to the sixth-order number. The aberration coefficients from the first-order to the fourth-order are the same as those in the original fourth-order aberration theory^[6]. The fifth-order and the sixth order aberration coefficients are w₁₄₀, w₃₂₀, w₅₀₀, w₆₀₀, w₄₂₀, w₂₄₀ and w₀₆₀respectively.

The wave aberration calculation is based on the first-order optical system parameters. The parameters of the main light ray in the optical path are calculated using the main light ray transfer equation. Based on the first-order and the second-order wave aberration coefficients (monochromatic aberration is calculated for D light), the object and the image distance of each optical plane in the meridian and sagittal plane are determined. These basic optical path parameters can be used to calculate the each order wave aberration coefficients of the optical system^[7].

The solution of the fifth-order aberration is the same as the original fourth-order wave aberration theory, but the correlation expression is required to keep the sixth-order accuracy. The first-order to the third-order aberration coefficients are consistent with the original fourth-order wave aberration theory. The fourth-order and fifth-order meridional and sagittal aberration coefficients are d₄₀₀, d₂₂₀, d₀₄₀, d₅₀₀, d₃₂₀, d₁₄₀, h₃₁₀, h₁₃₀, h₄₁₀, h₂₃₀, h₀₅₀ respectively.

As shown in Fig.1, it shows the relevant coordinate systems used in the main ray transmission and aberration calculation. In the process of calculation, all the relevant coordinate systems adopt the right hand rule, and the positive direction of y-axis is from the inside to the outside, x₁O₁z₁ and x₂O₂z₂ are the optical surface coordinate systems of M₁ and M₂ respectively. P₀P₁P₂ is any aperture ray, and A₁O₁O₂ is the main ray. $ {\alpha }_{1} $, $ \;{\beta }_{1} $, $ {\alpha }_{2} $ and $ \;{\beta }_{2} $ are the incident angle and refraction angle of the optical plane M₁ and M₂ respectively, which are determined by the main ray transfer equation. But in the correlation processing of other aberrations, the incident angle of each optical plane should be increased by 180°, because the incident or refraction angle of the ray is defined by reference to the z axis. In the calculation of wave aberration and other relevant aberration, the ray coordinates on the final optical plane are set as the pupil reference coordinates.

Figure 2 shows two element systems consisting of M₁ and M₂ optical surfaces. P₁P₂B is any aperture ray, A₁O₁O₂ is the main ray, x′₀₁O₁y′₀₁ is the outgoing wave array coordinate system of M₁, x₀₂O₂y₀₂ and x′₀₂O₂y′₀₂ are the incoming and outgoing wave array coordinate systems of M₂, and $ \overline {{d}_{1}} $ is the distance between M₁ and M₂ along the main light ray.

The coordinate (x₂, y₂) of the light on M₂ are the pupil reference optical coordinate system. The reverse optical path is used to calculate the wave aberration of optical surface M₂ on M₁, and to correct the pupil coordinates of M₁. In Fig.2, $ {x}_{2}^{*}{O}_{2}{y}_{2}^{*} $ represents the image coordinate system at the distance from M₂ is $ {r}_{0}^{*} $ in the calculation of the reverse light path. Applying the fifth-order aberration expression (only including the wave aberration contribution of M₂), and setting $ {r}_{0}^{*}=\overline {{d}_{1}} $, the light aberration ( $ {\Delta }{x}_{2}^{*},\;\;{\Delta }{y}_{2}^{*} $) of the image at O₁ in the coordinate system $ {x}_{2}^{*}{O}_{2}{y}_{2}^{*} $ is obtained, that is $ \overline {{B}_{1}{B}_{1}^{{'}}} $, and B₁ is only determined by the calculation of aperture ray according to the second-order accuracy transfer relationship, and B₁′ also includes the influence of derivative aberration. Convert $ {\Delta }{x}_{2}^{*} \;{\rm{and}}\;{\Delta }{y}_{2}^{*} $ to coordinate system x′₀₁O₁y′₀₁.

$ \Delta {x}_{01}^{{{'}}}=- {\Delta }{x}_{2}^{*} $  (1)

$ \Delta {y}_{01}^{{{'}}}= {\Delta }{y}_{2}^{*} $  (2)

Therefore

$ \begin{split} \Delta {x}_{01}^{{{'}}}=&-({d}_{200\left(2\right)}^{{{'}}}{x}_{2}^{2}+{d}_{020\left(2\right)}^{{{'}}}{y}_{2}^{2}+{d}_{300\left(2\right)}^{{{'}}}{x}_{2}^{3}+{d}_{120\left(2\right)}^{{{'}}}{{x}_{2}y}_{2}^{2}+\\ & {d}_{400\left(2\right)}^{{{'}}}{x}_{2}^{4}+{d}_{220\left(2\right)}^{{{'}}}{x}_{2}^{2}{y}_{2}^{2}+{d}_{040\left(2\right)}^{{{'}}}{y}_{2}^{4}+{d}_{500\left(2\right)}^{{{'}}}{x}_{2}^{5}+\\ & {d}_{320\left(2\right)}^{{{'}}}{x}_{2}^{3}{y}_{2}^{2}+{d}_{140\left(2\right)}^{{{'}}}{x}_{2}{y}_{2}^{4}) \end{split}$ (3)

$ \begin{split} \Delta {y}_{01}^{{{'}}}=&\left({h}_{110\left(2\right)}^{{{'}}}{x}_{2}{y}_{2}+{h}_{210\left(2\right)}^{{{'}}}{x}_{2}^{2}{y}_{2}+{h}_{030\left(2\right)}^{{{'}}}{y}_{2}^{3}+{h}_{310\left(2\right)}^{{{'}}}{x}_{2}^{3}{y}_{2}+\right.\\ & \left.{h}_{130\left(2\right)}^{{{'}}}{x}_{2}{y}_{2}^{3}+{h}_{410\left(2\right)}^{{{'}}}{x}_{2}^{4}{y}_{2}+{h}_{230\left(2\right)}^{{{'}}}{x}_{2}^{2}{y}_{2}^{3}+{h}_{050\left(2\right)}^{{{'}}}{y}_{2}^{5}\right) \end{split} $ (4)

In the above formula, $ {d}_{{{ijk}}\left(2\right)}^{{{'}}} $ and $ {h}_{{{ijk}}\left(2\right)}^{{{'}}} $ are the aberration coefficients with only M₂ wave aberration contribution. Due to the application of the reverse light path to calculate the light aberration, in the expression of $ {d}_{{{ijk}}\left(2\right)}^{{{'}}} $ and $ {h}_{{{ijk}}\left(2\right)}^{{{'}}} $, it is necessary to replace the relevant parameters in the aberration coefficients calculated by the forward optical path: $ {\;\beta }_{2}\to {\alpha }_{2} $, $ {r}_{{\rm{m}}2}^{{{'}}}\to {r}_{{\rm{m}}2} $, $ {r}_{{\rm{s}}2}^{{{'}}}\to {r}_{{\rm{s}}2} $.

Since the aberration Δx₀₁′ and Δy₀₁′ are generally very small, they are transformed into the surface coordinate system of M₁ by linear approximation.

$ {\Delta x}_{1}=\frac{{\Delta x}_{01}^{{{'}}}}{{\cos}{\beta }_{1}} $ (5)

$ {\Delta y}_{1}={\Delta y}_{01}^{{{'}}} $ (6)

The total wave aberration of multi-element optical system is the linear sum of each element wave aberration. When calculating the wave aberration with aberration theory, it is necessary to convert the light pupil coordinates on each element into the pupil reference coordinates of the final optical element. As shown in Fig.1, by using the spatial ray equation and polynomial fitting method (excluding the influence of wave aberration), the second-order accuracy relation between the coordinates P₁(x₁,y₁) and P₂(x₂,y₂) at the intersection of ray and optical surface M₁ and M₂ is obtained. As follows:

$ \begin{split} &{x_1} = {A_1}{x_2} + {A_1}{\varGamma _{R\left( 1 \right)}}x_2^2 + {A_1}{\varGamma _{\rho \left( 1 \right)}}y_2^2\\ & {y_1} = {B_1}{y_2} + {B_1}{\varGamma _{12\left( 1 \right)}}{x_2}{y_2} \end{split} $  (7)

In the formula above,

$ {\varGamma _{R\left( k \right)}} = {A_k}{c_{2,0\left( k \right)}}\tan {\beta _k} - {c_{2,0\left( {k + 1} \right)}}\tan {\alpha _{k + 1}} + \frac{{\sin \left( {{\alpha _{k + 1}} - {\beta _k}} \right)}}{{{r_{m\left( {k + 1} \right)}}\cos {\beta _k}}} $  (8)

$ {\varGamma _{\rho \left( k \right)}} = \frac{{B_k^2{c_{0,2\left( k \right)}}\tan {\beta _k}}}{{{A_k}}} - {c_{0,2\left( {k + 1} \right)}}\tan {\alpha _{k + 1}} $  (9)

$ {\varGamma _{12\left( k \right)}} = \frac{{\sin {\alpha _{k + 1}}}}{{{r_{s\left( {k + 1} \right)}}}} - \frac{{{A_k}\sin {\beta _k}}}{{{{r'}_{s\left( k \right)}}}} $  (10)

$ {A}_{k}={r}_{m\left(k\right)}^{{'}}\cos{\alpha }_{k+1}/\left({r}_{m\left(k+1\right)}\cos{\beta }_{k}\right) $  (11)

$ {B}_{k}=-{r}_{s\left(k\right)}^{{'}}/{r}_{s\left(k+1\right)} $ (12)

$ {c}_{20\left(k\right)}=1/\left(2{R}_{k}\right) $, $ {c}_{02\left(k\right)}=1/\left(2{\rho }_{k}\right) $, $ {R}_{k} $ and $\; {\rho }_{k} $ are the meridional and sagittal radii of curvature for the k-th optical plane. $ {\alpha }_{k} $ and $ {\;\beta }_{k} $ are the incident angle and refraction angle for the k-th optical plane. $ {r}_{m\left(k\right)} $, $ {r}_{m\left(k\right)}^{{'}} $, $ {r}^{{'}}_{s\left(k\right)} $ are the distance between object and image in meridian and sagittal plane for the k-th optical plane respectively. For the two elements system, the parameter subscript k=1 in the above expression. In equation (7), the first term on the right side of equal sign represents the linear approximation transfer relation of pupil coordinates between two adjacent optical surfaces. According to the above discussion, considering the correction of the derivative aberration and the second-order accuracy of the pupil coordinate to the pupil coordinate of the optical surface M₁, its actual wave aberration is obtained. As follows:

$\begin{split} {W_1} =& {w_{120\left( 1 \right)}}\left( {{x_1} + \Delta {x_1}} \right){\left( {{y_1} + \Delta {y_1}} \right)^2} + {w_{300\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^3}\; + \\ & {w_{400\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^4} + {w_{220\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^2}{\left( {{y_1} + \Delta {y_1}} \right)^2}\; + \\ & {w_{040\left( 1 \right)}}{\left( {{y_1} + \Delta {y_1}} \right)^4} + {w_{500\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^5} + \\ & {w_{140\left( 1 \right)}}\left( {{x_1} + \Delta {x_1}} \right){\left( {{y_1} + \Delta {y_1}} \right)^4} + \\ & {w_{320\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^3}{\left( {{y_1} + \Delta {y_1}} \right)^2}\; + \\ & {w_{600\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^6} + {w_{240\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^2}{\left( {{y_1} + \Delta {y_1}} \right)^4} + \\ & {w_{420\left( 1 \right)}}{\left( {{x_1} + \Delta {x_1}} \right)^4}{\left( {{y_1} + \Delta {y_1}} \right)^2} + {w_{060\left( 1 \right)}}{\left( {{y_1} + \Delta {y_1}} \right)^6} \end{split} $  (13)

By substituting the equations (5), (6) and (7) into (13) and expanding them in series with respect to (x₂, y₂), the wave aberration correction expressions caused by the derivative wave aberration and the second-order transmission accuracy of pupil coordinates can be gotten. In the following application of aberration theory to the imaging simulation of fish-eye lens, the former ignores items related to $ {{\varGamma }}_{R} $, $ {{\varGamma }}_{\rho } $, $ {{\varGamma }}_{12} $ and items $ {w}_{ij0\left(1\right)}{w}_{{i}^{{'}}{j}^{{'}}0\left(2\right)}^{2} $. The derived wave aberration coefficient has a linear combination form for $ {w}_{ij0\left(1\right)}{w}_{{i}^{{'}}{j}^{{'}}0\left(2\right)}^{2} $ correlation. The latter only consider the primary term of $ {{\varGamma }}_{R} $, $ {{\varGamma }}_{\rho } $, $ {{\varGamma }}_{12} $.

Referring to the solution method of the derivative aberration of the axisymmetric optical system^[8-9]. For an optical system composed of k+1 optical planes, when we study the derivative aberrations of the k+1 optical plane, we can regard the first to the k-th optical plane as the system1 (M₁), while the (k+1)-th optical plane can be regarded as system 2 (M₂). The physical basis of such processing is that the pupil aberration does not change the shape of the wave, but only change the light position coordinates of the wave. Therefore, it is necessary to replace the corresponding intrinsic aberrations of system 1 and system 2 in the above discussion by intrinsic aberration of the first to the k-th optical plane and the (k+1)-th optical plane.

For multi-element system, pupil coordinates are transferred between adjacent elements according to the equation (7). After induction, the relationship expression between the pupil coordinates and the pupil reference coordinates on each optical plane can be obtained. Combined with the above process, the wave aberration correction expression by the second-order accuracy of pupil coordinates of each optical surface is obtained. As follows:

$ w_{400}^{\left( S \right)} = 3\overline {w_{300}^{\varPsi 1}} $  (14)

$ w_{220}^{\left( S \right)} = 3\overline {w_{300}^{\varPsi 2}} + \overline {w_{120}^{\varPsi 1}} + 2\overline {w_{120}^{\varPsi 3}} $  (15)

$ w_{040}^{\left( S \right)} = \overline {w_{120}^{\varPsi 2}} $  (16)

$ w_{500}^{\left( S \right)} = 4\overline {w_{400}^{\varPsi 1}} $  (17)

$ w_{140}^{\left( S \right)} = 2\overline {w_{220}^{\varPsi 2}} + 4\overline {w_{040}^{\varPsi 3}} $  (18)

$ w_{320}^{\left( S \right)} = 2\overline {w_{220}^{\varPsi 1}} + 2\overline {w_{220}^{\varPsi 3}} + 4\overline {w_{400}^{\varPsi 2}} $  (19)

$ w_{600}^{\left( S \right)} = 5\overline {w_{500}^{\varPsi 1}} $  (20)

$ w_{240}^{\left( S \right)} = \overline {w_{140}^{\varPsi 1}} + 4\overline {w_{140}^{\varPsi 3}} + 3\overline {w_{320}^{\varPsi 2}} $  (21)

$ w_{420}^{\left( S \right)} = 5\overline {w_{500}^{\varPsi 2}} + 3\overline {w_{320}^{\varPsi 1}} + 2\overline {w_{320}^{\varPsi 3}} $  (22)

$ w_{060}^{\left( S \right)} = \overline {w_{140}^{\varPsi 2}} $  (23)

The expressions of $ \overline {{w}_{ij0}^{{{\varPsi }}{n}}} $, $ {{\varPsi }}_{1\left(i\right)} $, $ {{\varPsi }}_{2\left(i\right)} $ and $ {{\varPsi }}_{3\left(i\right)} $ are as follows:

$ \overline {w_{ij0}^{\varPsi n}} = \sum\limits_{l = 1}^{g - 1} {A_{l|g}^i} B_{l|g}^j{\varPsi _{n\left( l \right)}}{w_{ij0\left( l \right)}}\;\;\left( {n = 1{,}\;2{,}\;3} \right) $  (24)

$ {\varPsi _{1\left( i \right)}} = \sum\limits_{l = i}^{g - 1} {{A_{l + 1|g}}{\varGamma _{R\left( l \right)}}} $  (25)

$ {\varPsi _{2\left( i \right)}} = \sum\limits_{l = i}^{g - 1} {\frac{{B_{l + 1|g}^2}}{{{A_{l + 1|g}}}}} {\varGamma _{\rho \left( l \right)}} $  (26)

$ {\varPsi _{3\left( i \right)}} = \sum\limits_{l = i}^{g - 1} {{A_{l + 1|g}}} {\varGamma _{12\left( l \right)}} $  (27)

where

$ {A_{k|g}} = \frac{{{x_k}}}{{{x_g}}} = \frac{{{{r'}_{m\left( k \right)}}{{r'}_{m\left( {k + 1} \right)}} \cdots {{r'}_{m\left( {g - 1} \right)}}\cos {\alpha _{k + 1}}\cos {\alpha _{k + 2}} \cdots \cos {\alpha _g}}}{{{r_{m\left( {k + 1} \right)}}{r_{m\left( {k + 2} \right)}} \cdots {r_{m\left( g \right)}}\cos {\beta _k}\cos {\beta _{k + 1}} \cdots \cos {\beta _{g - 1}}}} $  (28)

$ {B_{k|g}} = \frac{{{y_k}}}{{{y_g}}} = {\left( { - 1} \right)^{g - k}} \cdot \frac{{{{r'}_{s\left( k \right)}}{{r'}_{s\left( {k + 1} \right)}} \cdots {{r'}_{s\left( {g - 1} \right)}}}}{{{r_{s\left( {k + 1} \right)}}{r_{s\left( {k + 2} \right)}} \cdots {r_{s\left( g \right)}}}} $  (29)

The total wave aberration of an optical system should be the sum of three items through the system wave aberration correction, that are the intrinsic wave aberration, the derivative wave aberration and the second-order accuracy of the pupil coordinates:

$ {w_{ij0\left( T \right)}}{\rm{ = }}{\bar w_{ij0}} + w_{ij0}^{\left( E \right)} + w_{ij0}^{\left( S \right)} $  (30)

The fish-eye lens system designed is a lens in the visible-light band with the F number that no more than 3.5, with the focal length no less than 5 mm when the system FOV angle is 180°. The modulation-transfer-function (MTF) requirement is that, when the space frequency is 60 lp/mm, the modulus of the optical transfer function (OTF) is no less than 0.5. The size of the image plane is approximately 18 mm, and the aspect ratio of the image plane is 4:3. The design specifications of the fish-eye lens system are listed in Tab.1.

Generally, fish-eye lens consists of several negative meniscus lenses for compressed FOV in the former group and several objective lenses in the latter group. The design idea is firstly to determine the parameters of the former group optical system according to the design constraints, and then select the structure and parameters of the latter group objective lenses according to the space field angle and aberration distribution of the former group optical system. The reason for adopting the above strategy is that the parameters of the former group optical system should be used to meet the relevant constraints. For example, the small radius of the negative meniscus lens should meet the requirements of the lateral and longitudinal size limits of the fish-eye lens. At the same time, some aberrations (mainly field curvature) of the lens should be controlled. According to the structure of the lens, many parameters can be used to optimize the system aberration.

The wave aberrations of the field curvature, axial chromatic aberration and vertical chromatic aberration of the former optical group are calculated by using the wave aberration theory of the plane symmetric optical system. According to the wave aberration balance condition of the former and the latter optical groups^[10], the optical power and optical spacing of each lens of the latter optical group objective lens (i.e. the first-order optical parameters) are calculated. The flow chart of fish-eye lens system design using the sixth-order wave aberration theory is shown in Fig.3.

By using the sixth-order wave aberration theory, a large aperture imaging simulation of a fish-eye lens optical system is carried out, and the simulation results are compared with the numerical results of ZEMAX ray tracing. The optical parameters of the fish-eye lens system are shown in Tab.2. The monochromatic aberration distribution of visible light at five field angles of 0°, 27°, 45°, 63° and 90° is calculated respectively. It is assumed that the F-number of the fish-eye lens system is 3.2. The distance between the object point and the first optical plane is infinite at different field angles. The fish-eye lens optical system with maximum field of view angle of 180° is designed. The Optical path diagram of the optical system is shown in Fig.4.

In order to improve the accuracy of aberration calculation, it is necessary to consider the sixth-order wave aberration and more accurate pupil coordinate transfer processing. Through the calculation and analysis of this fish-eye lens systems similar to the above cases, the results show that the wave aberrations correction including higher-order aberrations and the second-order accuracy of pupil coordinates significantly improves the accuracy of aberrations calculation^[11], which is related to the specific design of the system.

It is also an auxiliary method to verify the wave aberration theory to analyze whether the simplified results of aberration expression in processing paraxial optical system are consistent with the previous results. In the sixth-order wave aberration theory, if $ \alpha ={180}^{ \circ } $, $\; \beta ={0}^{ \circ }, $ and the optical surface is a sphere, it represents the aberration of the object point on the axis of the axisymmetric optical system. As a result, the fourth-order intrinsic wave aberration coefficient is the same as that of Seidel, the third-order and fifth-order are zero, the sixth-order intrinsic wave aberration coefficient is the same as that of w₀₆₀ reported in reference^[13]. Moreover, the expression of the aberration in the direction of meridian and sagittal also becomes the same, which conforms to the geometric symmetry with axial symmetry optical system.

The modulation transfer function (MTF) curve is not only related to the optical systems aberrations, but also to the diffraction effects of optical systems. It has the advantages of objectivity and reliability to evaluate the imaging quality of optical systems. It can be applied to optical systems with small aberrations and optical systems with large aberrations at the same time.

Figure 5 shows the fish-eye lens system with the focal length of 5.989 mm and 180° FOV angle. The modulus of the OTF is no less than 0.561 in full FOV when the spatial frequency is 60 lp/mm.

This fish-eye lens system can satisfy the requirement of minimum the modulus of the OTF 0.561 when the spatial frequency is 60 lp/mm, and the MTF curve is smooth, and the imaging quality of the edge and central FOV is better guaranteed. The imaging quality can be maintained stably, and the imaging effect of millions of pixels of resolution can be obtained^[14].

Figure 6 shows the field curvature and distortion curves of the fish-eye system at focal lengths of 5.989 mm. It can be seen that the absolute field curvature of all FOV does not exceed 0.037 mm, and the maximum distortion of all the field of view is 4.34%.

For a panoramic lens, the distortion itself is very large. Such a large distortion can produce special imaging effects in practical application. Large distortion is in high demand in some specialized fields. Such a large distortion can be further corrected^{[5, 15]}.

Figure 7 shows the relative illumination curve of the fish-eye lens system at focal lengths of 5.989 mm. Comprehensive analysis shows that the relative illumination of the edge is no less than 0.9, which is mainly due to the decrease of the edge illumination due to the setting of the vignetting coefficient. The uniformity of the image surface illumination is very good. The system guarantees the uniformity of illumination of the image plane and the problem of illumination attenuation caused by a large FOV has been solved.

The method by developing the derivative wave aberration of plane symmetric optical system is given, that is different from the method by developing the derivative wave aberration of paraxial optical system, and its calculation expression is given. The large aperture imaging for such system is put forward, whose pupil coordinates are in the adjacent optical plane. In order to meet the accuracy requirement of wave aberration calculation, the second-order polynomial relation is used to deal with the transfer between them. The results show that the sixth-order wave aberration can objectively reflect the pupil aberration of large aperture imaging in this kind of fish-eye lens system.

The aberration expression is not a function of the field angle, it reflects the pupil aberration of each field object point, but it cannot be used to calculate the field aberration directly, such as the field curvature aberration. We use the main ray transfer equation to calculate the parameters of the main ray, and then calculate the field curvature aberration on the image surface from the geometric relationship, and then transmit it to the corresponding wave aberration. The disadvantage is that we can't deal with the field curvature aberration of any single optical surface in the optical system. In addition, there are similar limitations in the treatment of axial chromatic aberration and transverse chromatic aberration. Therefore, the calculation methods of field curvature and chromatic aberration to overcome the above limitations need to be further explored.

The sixth-order wave aberration theory should be further explored to study the laws of large aperture imaging for such fish-eye lens optical systems, including the balance between derived aberration and intrinsic aberration, and the relationship between optical system structure and parameters. The design laws of such fish-eye lens optical systems, including the balance between derived aberration and systems should be summarized, and the new fish-eye lens system should be developed.