Next Article in Journal
Generation of Perfect Vortex Beams with Complete Control over the Ring Radius and Ring Width
Previous Article in Journal
Slot-Waveguide Based All-Optical RF Spectrum Analyzer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination Position and Initial Value of Aspheric Surface for Fisheye Lens Design

1
School of Aeronautical and Mechanical Engineering, Changzhou Institute of Technology, Changzhou 213032, China
2
Department of Precision Mechanical Engineering, Shanghai University, Shanghai 200444, China
3
School of Mechanical Engineering and Rail Transit, Changzhou University, Changzhou 213164, China
4
The Engineering Research Center for CAD/CAM of Fujian Universities, Putian University, Putian 351100, China
*
Authors to whom correspondence should be addressed.
Photonics 2023, 10(12), 1381; https://doi.org/10.3390/photonics10121381
Submission received: 7 November 2023 / Revised: 9 December 2023 / Accepted: 13 December 2023 / Published: 15 December 2023

Abstract

:
The aspheric surface is a commonly used method to improve the imaging quality of the fisheye lens, but it is difficult to determine the position and initial value. Based on the wave aberration theory of the plane-symmetric optical system, a method of using an aspheric surface to design a fisheye lens is proposed, which can quickly determine the appropriate aspheric surface to improve the imaging performance. First, the wave aberration of each optical surface of the fisheye lens is calculated and its aberration characteristics are analyzed. Then, a numerical evaluation function is reported based on the aberration distribution of the fisheye lens on the image plane. According to the functional relationship between the evaluation function and the aspheric coefficient, the position of the aspheric surface and the initial value of the aspheric coefficient can be calculated. Finally, the adaptive and normalized real-coded genetic algorithm is used as the evaluation function to optimize the fisheye lens using an aspheric surface. The proposed method can provide an effective solution for designing a fisheye lens using an aspheric surface.

1. Introduction

In fisheye lens system design, spherical surfaces are the most common because of their simple expressions and easier fabrication. But the parameter optimization of a spherical surface only considers the radius of curvature. The aspheric surface needs to determine the radius of the vertex, the quadratic cone coefficient, and other parameters, which makes the aspheric surface have more freedom for parameter optimization than the spherical surface. Thus, using aspheric surfaces can easily correct the optical system aberration and improve the imaging performance [1,2,3,4,5,6]. However, the processing, inspection, and commissioning of aspheric surfaces are more complicated than spherical surfaces, making the manufacturing cost of optical systems significantly higher [7,8]. To balance the advantages, this paper uses an aspheric surface for a fisheye lens. Reasonable determining position and initial value of the aspheric surface will be significant for improving imaging quality [9].
In the design of the fisheye lens, there are few reports to determine the optimal position and initial value of the aspheric surface. At present, the position of the aspheric surface is usually determined based on experience without theoretical guidance. The fisheye lens usually has a field angle of 180° or even large, and they belong to the plane-symmetric optical systems [10,11]. The Seidel aberration theory cannot be used in aberration analysis. In recent years, Lu et al. have developed an aberration theory of plane-symmetric optical systems by applying the wave aberration method [10], which can analyze the relationship between structural parameters, including aspheric coefficients, and the imaging performance of a fisheye lens. This theory provides a theoretical basis for determining the position and initial value of the aspheric surface. Based on the theory of plane-symmetric optical systems, we report a method to determine the optimal aspheric position in a fisheye lens and solve for its aspheric coefficient.
The paper is organized as follows: Section 2 introduces the wave aberration theory of plane-symmetric optical systems of an ultrawide-angle optical system. Section 3 describes the calculation method of wave aberration distribution of a single optical plane in a fisheye lens. In Section 4, we define an evaluation function based on the aberrations distribution on the image plane, then determine the sensitive aspheric surfaces and corresponding aspheric coefficients in the fisheye lens system. In Section 5, the method is verified by optimizing a fisheye lens, and compared with the original fisheye lens using the MTF (modulation transfer function) curve. Section 6 summarizes the full text.

2. Wave Aberration Theory of the Plane-Symmetric Optical System

2.1. The Chief Ray Transfer Equation

The fisheye lens usually has a very large field angle and no longer approximates a paraxial ray. The tracing of a chief ray should use the exact geometrical relationships of trigonometric tracing formulae of a meridian oblique-incident ray. Figure 1 shows a transmission schematic diagram of an off-axis chief ray in the meridian plane. It shows a chief ray A i O i O i + 1 O i + 2 ¯ with initial field angle ω i 1 to be refracted by one mirror ( k = i ) and then refracted by two optical surfaces ( k = i + 1 , k = i + 2 ). The chief ray passes through optical media with refractive indices n k , and it impinges optical surface k at O k and intersects the optical axis at M k . ω k means the field angle which is included between the chief ray and the optical axis. The sign of the field angle is identical to the chief ray of the normal in the optical surface. α k means the angle of incidence, and β k means the angle of reflection or refraction. The sign is identical to the sign of the chief ray to the normal of the optical surface. Γ k and ρ k are the meridional and sagittal curvature radii of the optical surface k at O k . The distance from the optical surface k to k + 1 is D k D k + 1 ¯ = d k . Their signs are positive on the right side of the optical surface. On the left side, they are negative.
From Figure 1, the parameters of the chief ray at any optical surface can be calculated with its transfer equation [12]:
sin α i + 1 = Γ i + 1 + d i Γ i ρ i + 1 sin ω i + ρ i ρ i + 1 sin β i ,
ω i = ω i 1 + β i α i = ω 0 + i = 1 i ( β i α i ) ,
β i + 1 = sin 1 ( n i n i + 1 sin α i + 1 ) .

2.2. Transformation of Figure Equation Coordinate System

To apply the plane-symmetric aberration theory [13,14] for calculating the wave aberration of an off-axis point object, we first need to take the coordinate system x y z as the coordinate origin O , whose origin O is an intersection of the chief ray and the optical surface, as shown in Figure 2. It shows a chief ray is reflected by a quadric of revolution. The z axis is the normal to the surface and the x axis is along the tangent direction in the meridional plane. θ is the angle between normal at O and the optical axis.
From Figure 2, the expression of the quadrics of revolution in the coordinate system of x y z is expressed as
x 2 + y 2 = a 1 z + a 2 z 2 ,   a 1 = 2 R 0 ,
where R 0 is the curvature radius of the optical surface at point O , R 0 = C 1 O ¯ = Γ 1 . a 2 is a constant that determines the type of quadric of revolution. a 2 > 0 , a 2 = 0 , 1 < a 2 < 0 , a 2 = 1 and a 2 < 1 in Equation (4) represent the quadratic curve as a hyperboloid, a paraboloid, a prolate ellipsoid, a sphere, and an oblate ellipsoid, respectively.
k represents aspheric characteristics in Zemax. The functional relationship between a 2 and k is
a 2 = 1 k .
The transformation relation from the coordinate systems of x y z to x y z is
x = x cos θ + z sin θ + x 0 * , y = y , z = x sin θ + z cos θ + z 0 * . }
Substituting Equation (6) into Equation (4), we can get the figure equation of the optical surface in the coordinate system of x y z ; then, its fourth-order Taylor series expansion is
z = c 2 , 0 x 2 + c 0 , 2 y 2 + c 3 , 0 x 3 + c 1 , 2 x y 2 + c 4 , 0 x 4 + c 0 , 4 y 4 + c 2 , 2 x 2 y 2 ,
where c i , j is the figure coefficient [12],
c 2 , 0 = a 1 2 B 3 , c 2 , 2 = 2 B 7 ( a 1 2 C 2 + 8 A 2 x 0 * 2 ) , c 0 , 2 = 1 B , c 4 , 0 = a 1 2 ( a 1 2 C 2 + 16 A 2 x 0 * 2 ) B 9 , c 0 , 4 = C 2 B 5 , c 3 , 0 = 4 A a 1 2 x 0 * B 6 , c 1 , 2 = 4 A x 0 * B 4 ,
with
A = ( 1 + a 2 ) 4 a 2 x 0 * + a 1 2 , B = ± 4 x 0 * 2 ( 1 + a 2 ) + a 1 2 , C = 4 x 0 * 2 a 2 ( 4 a 2 x 0 * 2 + a 1 2 ) .
From Figure 2, the parameters of Γ 1 and ρ 1 can be derived as
Γ 1 = C 1 O ¯ = z 0 * x 0 * tan θ = a 1 2 + ( 1 + a 2 ) z 0 * , ρ 1 = O C ¯ 1 = 4 x 0 * 2 + ( a 1 + 2 a 2 z 0 * ) 2 2 = 1 2 c 0 , 2 .
The coordinate of O   ( x 0 * , y 0 * ) can be obtained by the geometric relationship of the incident chief ray O M 0 ¯ to the optical surface.
z 0 * = a 1 + 2 S 0 tan 2 ω 0 ± a 1 2 + 4 S 0 ( a 1 + a 2 S 0 ) tan 2 ω 0 2 ( a 2 tan 2 ω 0 ) , x 0 * = tan ω 0 ( z 0 * S 0 ) .
The sign before the root in Equation (11) should be taken as positive for a 1 < 0 ; otherwise, it is negative. The sign of B in Equation (9) should be taken as positive in the case of reflection by a convex mirror or refraction by a concave optical surface, otherwise, it is negative [13].

3. Calculation Method of the Wave Aberration of a Fisheye Lens

3.1. Aperture Ray Wave Aberrations of an Off-Axis Point Object

For an optical system of g elements, the total wave aberration is [10]
W = 1 g i j 4 w i j 0 x i y j , ( i + j 4 ,   j   is   even ) ,
with
w i j 0 = n M i j 0 ( α , r m , r s , 0 ) + n M i j 0 ( β , r m , r s , 0 ) .
In Equation (13), n and n mean the refraction index of the object and image space. M i j 0 means the wave aberration coefficient of object beam pencil, which is calculated by Equations (11) and (13)–(16) in ref. [14]. r m ( i ) and r s ( i ) mean the meridional and sagittal focal distances along the chief ray in the image space of the gth optical surface; they can be obtained with the conditions of the wave aberration coefficients, and w 200 = 0 and w 020 = 0 [10] are applied sequentially to every optical surface as
r m ( i ) = n i cos 2 β i 2 c 2 , 0 ( n i cos α i + n i cos β i ) n i cos 2 α i r m ( i ) , r s ( i ) = n i 2 c 0 , 2 ( n i cos α i + n i cos β i ) n i r s ( i ) .
The distance between object points and surface i + 1 along the chief ray is given by
r m ( i + 1 ) = d ¯ i r m ( i ) , r s ( i + 1 ) = d ¯ i r s ( i ) ,
where d i ¯ means the optical spacing between the ith and (i + 1)th optical surfaces along the chief ray.
d ¯ i = ρ i sin ( ω i 1 α i ) ρ i + 1 sin ( ω i α i + 1 ) sin ω i .
The contributions of the spherical aberration W s p h and the coma aberration W c o m a to any optical surface are
W s p h = w 400 x 4 + w 220 x 2 y 2 + w 040 y 4 , W c o m a = w 300 x 3 + w 120 x y 2 . .

3.2. Chief Ray Wave Aberration of an Off-Axis Point Object

The wave aberrations of the field curvature contributed from the first to the gth optical surface in the meridional and sagittal direction are [9]
W m ( g ) = n g cos 2 β g 2 ( 1 r m ( g ) 1 r c ( g ) ) x g 2 , W s ( g ) = n g 2 ( 1 r s ( g ) 1 r c ( g ) ) y g 2 . .
The wave aberrations of the axial chromatic and transverse chromatic from the first to the gth optical surface are [9]
W C L ( g ) = n g 2 ( 1 r F ( g ) 1 r C ( g ) ) x g _ 0 2 ,
W C T ( g ) = n g x g _ 0 r 0 ( g ) ( y F ( g ) cos ω F ( g ) y C ( g ) cos ω C ( g ) ) .
In Equation (19), r F ( g ) and r C ( g ) are the focal distances in the image space of F and C light at normal incidence, and they can be calculated with Equation (14) in the case of α = 180 ° and β = 0 ° . In Equation (20), the height of the chief ray on the image plane for a light of some color is [9,10]
y g = tan ω g ( ρ g sin β g sin ω g Γ g + r 0 ( g ) ) ,
where r 0 ( g ) represents the distance from the gth optical surface to the Gaussian image plane.
The contributions of these types of wave aberrations from the first to the (g − 1)th optical surface are also calculated with Equations (18)–(20). The wave aberration of field curvature and chromatic contribution from the gth optical surface should be calculated by
W m ( g ) = W m ( g ) W m ( g 1 ) , W s ( g ) = W s ( g ) W s ( g 1 ) , W C L ( g ) = W C L ( g ) W C L ( g 1 ) , W C T ( g ) = W C T ( g ) W C T ( g 1 ) .
In Equations (17)–(20), ( x k , y k ) is the mapping coordinates of the aperture ray of the beam projected on the kth optical surface, and they are obtained by sequential linear approximated from the aperture stop radius [10],
x k = A k x k + 1 , y k = B k y k + 1 ,
where,
A k = r m ( k ) cos α k + 1 r m ( k + 1 ) cos β k , B k = r s ( k ) r s ( k + 1 ) .

4. Determining the Position of Aspheric Surfaces and Aspheric Coefficients

In this section, the optimal aspheric position and initial value of a 2 are solved by the functional curve between the evaluation function and the aspheric coefficient.

4.1. Evaluation Function

We define an evaluation function that sums types of aberrations at n field angles in its working range:
Q = k = 1 n ε k ( Q x ( k ) 2 + Q y ( k ) 2 + η k Q η ( k ) 2 + μ k Q c ( k ) 2 ) ,
and
Q x ( k ) 2 = 8 π W q L W q / 2 W q / 2 0 1 4 x 2 / W q 2 L / 2 ( x x ¯ ) 2 d x d y , Q y ( k ) 2 = 8 π W q L W q / 2 W q / 2 0 1 4 x 2 / W q 2 L / 2 y 2 d x d y ,
with
x ¯ = 8 π W q L W q / 2 W q / 2 0 1 4 x 2 / W q 2 L / 2 x d x d y
where ε k and μ k mean the corresponding weight factors of the axial chromatic and transverse chromatic; Q x ( k ) and Q y ( k ) are contributed by the aperture aberration of the kth field angle on the image plane, which is obtained by integrating and summing the area covered by the beam [13], respectively.
The evaluation function (Equation (25)) is the sum types of aberrations at n field angles in its working field ranges. This includes the effect of aspherical surfaces on the imaging performance of optical systems.
In Equations (26) and (27), the aperture ray aberrations on the image plane G are calculated by [10]:
x = 1 cos ω g ( d 100 x g + d 200 x g 2 + d 020 y g 2 + d 300 x g 3 + d 120 x g y g 2 ) , y = h 010 y g + h 110 x g y g + h 210 x g 2 y g + h 030 y g 3 ,
where x g and y g mean the ray coordinates on the last optical surface. Since the light beam from an object of an ultrawide-angle impinges on the optical surface with a large incidence angle, the focal line position in the meridian plane and the sagittal plane is severely deviated, and the aperture ray is usually elliptical when transmitted to the last optical surface. Therefore, the integration area in Equations (26) and (27) are the elliptical domain, and W q and L mean the projection length of the aperture ray on the final optical surface along the meridional (x) and sagittal (y) directions, respectively.
Q η ( k ) and Q c ( k ) in Equation (25) mean the numerical estimation of the transverse aberration on the image plane of the axial color aberration and lateral color aberration, respectively.
Q η ( k ) = | 2 x g _ 0 ( r F ( g ) r C ( g ) ) r 0 ( g ) | ,
Q c ( k ) = | y F ( g ) y C ( g ) | .
In the aberration theory of plane-symmetric optical systems, we calculated aberrations with the parameters of the chief ray parameters ω ,   α ,   β ,   r m ,   r s , r m , and r s of each optical surface. These parameters are all related to the figure coefficients c i , j . According to Equations (8) and (9), c i , j is determined by the parameters a 1 and a 2 of the optical surface. a 1 = 2 R 0 represents the radius of curvature of the optical surface at the vertex, and a 2 represents the aspheric coefficient of the optical surface. Therefore, the evaluation function Q in Equation (25) is related to the aspheric coefficient a 2 . We calculate the variation curve of Q with a 2 for each optical surface to determine the optimal position of the aspheric surface and the corresponding initial value of a 2 . This enables the optical system to have the best imaging performance (that is, Q is the minimum).

4.2. Wave Aberration Distribution of the Optical Lens

We will apply the aberration theory of plane-symmetric optical systems [9] to calculate the wave aberration distribution of two lenses, lens I and fisheye lens II, as shown in Figure 3 and Figure 4, and Table 1 and Table 2 list their optical parameters. The characteristic parameters of the systems are the focal lengths of 33.7 mm and 14.6 mm, the diameters of the aperture stop of 5 mm and 6 mm, and the maximum working field angle of view of 20 ° and 80 ° , respectively.
We use Equations (17)–(22) to calculate the wave aberration distribution of each optical surface of lens I at ω 0 = 20 ° and fisheye lens II at ω 0 = 80 ° . Figure 5 and Figure 6 show W m , W s , W C L , W C T , W s p h , and W c o m a , respectively.
From Figure 5 and Figure 6, it can be seen that the total wave aberrations of lens I and fisheye lens II are more serious. The effects of field curvature, lateral color aberration, and coma aberration are particularly significant. The position of the fifth surface is the aperture stop, and its aberration contribution is zero.
We use Equation (25) to calculate the curve of the Q with the aspheric coefficient a 2 of each optical surface of lens I, as shown in Figure 7. Since a 2 in the range of [–3, 1] can cover common easily machinable aspheric surface types [15]. The search range of a 2 is [−3, 2] with a step of 0.01, and the field angles in Equation (25) are set to 0 ° , 5 ° , 10 ° , 15 ° , and 20 ° . Weighting factors ε k , η k , and μ k are set to 1.
Similarly, the curve of the Q with a 2 of each optical surface of fisheye lens II is shown in Figure 8. The search range of a 2 is [−4, 2] with a step of 0.01, and the field angles are set to 0 ° , 10 ° , 20 ° , 30 ° , 40 ° , 50 ° , 60 ° , 70 ° , and 80 ° . ε k ,   η k , and μ k are set to 1.
From Figure 7, the second and third optical surfaces of lens I are more sensitive to their imaging performance (making Q reach a minimum). From Figure 8, the second and fourth optical surfaces of fisheye lens II are more sensitive to their imaging performance (making Q reach a minimum), respectively. The corresponding initial value of a 2 is shown in Table 3.
The wave aberration distributions of lens I and fisheye lens II in the case of ω 0 = 20 ° and ω 0 = 80 ° are recalculated with Equations (17)–(22), respectively. As shown in Figure 9 and Figure 10, only the a 2 of the sensitive optical surface is changed according to Table 3. Figure 9a,b show lens I with a 2 2 = 1.69 and a 2 3 = 0.79 . Figure 10a,b show fisheye lens II with a 2 2 = 2.12 and a 2 4 = 3.42 . Compared with Figure 5 and Figure 6, their total wave aberrations decreased significantly, respectively.
Compared to Figure 9 and Figure 10 with Figure 5 and Figure 6, in lens I and fisheye lens II, only the aspheric coefficient of one sensitive surface is modified, which makes the total wave aberration of these two lenses significantly reduced. It is confirmed that the use of an aspheric surface at a sensitive optical surface has a significant effect on correcting lens aberrations.
Lens I and fisheye lens II have not been designed and optimized. They only contain three and four lenses, respectively, so their aberrations are more serious. Lens I and fisheye lens II are difficult to optimize the lens system with better performance. We will optimize the fisheye lens using an aspheric surface in the next section.

5. Numerical Validation

To verify the performance, we optimize fisheye lens III using an aspheric surface in the following steps: the sensitive optical surface in fisheye lens III and the corresponding initial value of a 2 are solved; then, based on the MTF calculated by the aberration theory as the evaluation function, a self-adaptive and normalized real-coded genetic algorithm is used to optimize the optical parameters of fisheye lens III. Note that the optimization parameters do not include the refractive index of the lens’s material.

5.1. Determining the Position and Initial Value of Aspheric Surface in the Fisheye Lens

Fisheye lens III [16] is composed of 18 standard optical spheres, as shown in Figure 11, where the 11th surface is the aperture stop. Table 4 lists optical parameters, with a focal length of 9 mm, F/# = 5, 2 ω 0 = 160 ° .
Taking the aspheric coefficient a 2 of each optical surface as a variable in sequence, the evaluation function Q of fisheye lens III is calculated by Equation (25). The minimum value of each curve and the corresponding initial value of a 2 are calculated, and the curve is shown in Figure 12. We set the field angles in Equation (25) to be 0 ° , 10 ° , 20 ° , 30 ° , 40 ° , 50 ° , 60 ° , 70 ° , and 80 ° . The weighting factors of ε k , η k , and μ k are set to 1.
In Figure 12, the sixth optical surface using an aspheric surface is most sensitive to the imaging performance of fisheye lens III. It is a more significant effect than the other 17 optical surfaces. Therefore, we determine the sixth optical surface as the best aspheric surface position in fisheye lens III, where a 2 ( 6 )  is  1.04 and Q ( 6 )  is  0.0385 .

5.2. Optimizing Fisheye Lens III Using an Aspheric Surface

Based on the fourth-order wave aberration theory of the plane-symmetric optical systems, we use the MTF as the evaluation function and apply a self-adaptive and normalized real-coded genetic algorithm to optimize the optical parameters of fisheye lens III. The evaluation function is defined as [12]
Q = η Q m + ξ Q s ,
and
Q m = 1 k ε i MTF m ( i ) , Q s = 1 k ε i MTF s ( i ) ,
where Q m and Q s mean the evaluation function; and M T F m ( i ) and M T F s ( i ) mean modulation transfer function in the meridional and sagittal directions, respectively. ε i is the corresponding weight factor at the ith field angle. In this optimization, the field angle and the corresponding ε i are listed in Table 5. We set the weighting factors of the meridional and sagittal directions as η = 4 , ξ = 1 , respectively.
The optimization of fisheye lens III is performed in two rounds. In the first round, we search for the optimal values of the upper limit is set to 1.7 times, and the lower limit is set to 0.5 times the original radius of fisheye lens III. The upper limit is set to 1.3 times and the lower limit is set to 0.7 times the original optical distance of fisheye lens III; the search range for a 2 is set to [−2, 0]. The optimal position of the image plane is set to a range of ± 0.5 mm around r 0 .
The number of generations in the genetic algorithm is taken to be 600. We run 20 times for the original fisheye lens III and choose two designs of good performance as the starting designs for the second round of optimization.
In the second round of optimization, we take a relatively small search range of parameters, setting the radius to ± 10 % and the optical distances to ± 5 % of the starting designs. Around ± 0.2 of the starting design, the search range of a 2 is set. The optimal position of the image plane is in a small range of ± 0.2 mm around r 0 . The number of generations in the genetic algorithm is taken to be 200, and each starting design is run 10 times. The optimization result for the aspheric coefficient is a 2 ( 6 ) = 1.065 . Consequently, we obtained the best design of the fisheye lens from the optimization, and Table 4 lists its optical parameters. We apply the FFT method of Zemax [17,18] to calculate the MTF curves of the original and optimized fisheye lens III, respectively. We apply Zemax to calculate its MTF curves for the spatial frequencies of 10 lp/mm and 30 lp/mm in the case of F/# = 5 as shown in Figure 13; (a) shows the MTF curve of the original fisheye lens III; (b) shows the MTF curve of the optimized fisheye lens III. By comparing the MTF curves, the proposed approach in the paper makes the imaging performance of fisheye lenses superior. Our next step will be actual production and further testing and research.
From Table 4, compared with the original optical parameters of fisheye lens III, the optimized parameters have little change. The aspheric coefficient a 2 ( 6 ) = 1.065 of the sixth optical surface is close to the standard spherical surface. From Figure 13, we obtain the best design of fisheye lens III from the optimization.

6. Conclusions

Based on the wave aberration theory of plane-symmetric optical systems, this paper develops a method by which to determine the position and initial value of the aspheric surface, which provides an effective means for the design and optimization of fisheye lens systems. Through aberration calculation and example verification, the conclusions are as follows:
(1)
In fisheye lenses, the optical surface with a small curvature radius generally adopts an aspheric surface, which is more advantageous for improving its imaging performance. The aspheric coefficient is also closer to the spherical surface, making it easy to manufacture.
(2)
In the case of a moderate acceptance aperture of fisheye lenses, the contribution of the aspheric surface to the balance of field curvature aberration is more significant.
The method in this paper is also suitable for the optimal design of paraxial optical imaging systems and other wide angle lens systems. We hope that this work may be helpful to optical designers in creating high-quality fisheye lens systems.

Author Contributions

Conceptualization, L.F. and L.L.; methodology, L.F.; software, L.F. and K.Y.; validation, G.Q. and L.F.; formal analysis, S.G. and H.Z.; investigation, L.F.; resources, L.F.; data curation, K.Y.; writing—original draft preparation, L.F.; writing—review and editing, L.F.; visualization, G.Q.; supervision, L.L. and K.Y.; project administration, K.Y.; funding acquisition, K.Y. All authors have read and agreed to the published version of the manuscript.

Funding

National Natural Science Foundation of China (NSFC) (61975111); National Natural Science Foundation of China (NSFC) (52205552); Open Project of The Engineering Research Center for CAD/CAM of Fujian Universities (K202206); Basic Research Program of Jiangsu Province (Natural Science Foundation)—Youth Fund Project (BK20200983).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Conflicts of Interest

The author declares that there is no conflict of interest related to this article.

References

  1. Shannon, R.R. The Art and Science of Optical Design; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  2. Li, X.; Cen, Z.; Optics, G. Aberrations and Optical Design; Zhejiang University: Hangzhou, China, 2014. [Google Scholar]
  3. Yabe, A. Optimal selection of aspheric surfaces in optical design. Opt. Express 2005, 13, 7233–7242. [Google Scholar] [CrossRef] [PubMed]
  4. Yabe, A. Global Optimization with Traveling Aspheric-Aspheric Surface Number as Continuous Variable. In Proceedings of the International Optical Design Conference and Optical Fabrication and Testing, OSA Technical Digest (CD), Jackson Hole, WY, USA, 13–17 June 2010; Optical Society of America: Washington, DC, USA, 2010. [Google Scholar]
  5. Feng, Y.; Cheng, H.; Wang, T.; Dong, Z.; Tam, H. Optimal strategy for fabrication of large aperture aspheric surfaces. Appl. Opt. 2014, 53, 147–155. [Google Scholar] [CrossRef] [PubMed]
  6. Zhu, W.; Duan, F.; Tatsumi, K.; Beaucamp, A. Monolithic topological honeycomb lens for achromatic focusing and imaging. Optica 2022, 9, 100–107. [Google Scholar] [CrossRef]
  7. Pan, J. Design, Manufacture and Test of the Aspherical Optical Surfaces; Suzhou University: Suzhou, China, 2004. [Google Scholar]
  8. Jing, H.; Wang, S.; Yan, L.; Yan, L.; Shi, C.; Cha, J. Manufacturing of an aluminum alloy off-axis parabolic mirror. Def. Manuf. Technol. 2019, 4, 31–35. [Google Scholar]
  9. Fan, L.; Lu, L. Calculation of the wave aberration of field curvature and color aberrations of an ultrawide-angle optical system. Optics Commun. 2021, 479, 126414. [Google Scholar] [CrossRef]
  10. Lu, L. Aberration theory of plane-symmetric grating systems. J. Synchrotron Radiat. 2008, 15, 399–410. [Google Scholar] [CrossRef] [PubMed]
  11. Chrisp, M.P. Aberrations of holographic toroidal grating systems. Appl. Opt. 1983, 22, 1508–1518. [Google Scholar] [CrossRef]
  12. Fan, L.; Lu, L. Design of a simple fisheye lens. Appl. Opt. 2019, 58, 5311–5319. [Google Scholar] [CrossRef]
  13. Lu, L.-J.; Hu, X.-Y.; Sheng, C.-Y. Optimization method for ultra-wide-angle and panoramic optical systems. Appl. Opt. 2012, 51, 3776–3786. [Google Scholar] [CrossRef] [PubMed]
  14. Lu, L.; Deng, Z. Geometric characteristics of aberrations of plane-symmetric optical system. Appl. Opt. 2009, 48, 6946–6960. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, B.; Li, D.; Zhang, S.; Wang, J. Design of aspheric fisheye lens and study of distortion correction algorithm. Acta Optica Sinica 2014, 34, 222–227. [Google Scholar]
  16. The Inst. of Optics, Handbook of Lenses, No. 6; (National Defense Technology, 1983).
  17. Zemax User’s Guide; Zemax Development Corporation: Bellevue, WA, USA, 2007.
  18. Lin, X. ZEMAX Optical Design Super Learning Manual; Posts and Telecommunications Press: Beijing, China, 2014. [Google Scholar]
Figure 1. The diagram shows a chief ray A i O i O i + 1 O i + 2 ¯ passing through a reflecting and two refracting optical surfaces. The chief ray passes through the media of refractive indices n i 1 , n i , n i + 1 , and n i + 2 in turn. It impinges optical surface k at O k and intersects the optical axis at M k . ω k means the field angle. α k and β k mean the angle of incidence and reflection or refraction of optical surface k . The optical surfaces are assumed to be quadrics of revolution. The Γ k and ρ k are the meridional and sagittal curvature radius of the optical surface k at O k . The spacing between the optical surface k and k + 1 is D k D k + 1 ¯ = d k .
Figure 1. The diagram shows a chief ray A i O i O i + 1 O i + 2 ¯ passing through a reflecting and two refracting optical surfaces. The chief ray passes through the media of refractive indices n i 1 , n i , n i + 1 , and n i + 2 in turn. It impinges optical surface k at O k and intersects the optical axis at M k . ω k means the field angle. α k and β k mean the angle of incidence and reflection or refraction of optical surface k . The optical surfaces are assumed to be quadrics of revolution. The Γ k and ρ k are the meridional and sagittal curvature radius of the optical surface k at O k . The spacing between the optical surface k and k + 1 is D k D k + 1 ¯ = d k .
Photonics 10 01381 g001
Figure 2. The diagram shows a chief ray reflected by a quadric of revolution. The coordinate system of x y z . The coordinate system x y z with O as the origin is transformed into a coordinate system x y z , whose origin is at the intersection point of the chief ray to the optical surface O . θ is the angle between the normal at O and the optical axis.
Figure 2. The diagram shows a chief ray reflected by a quadric of revolution. The coordinate system of x y z . The coordinate system x y z with O as the origin is transformed into a coordinate system x y z , whose origin is at the intersection point of the chief ray to the optical surface O . θ is the angle between the normal at O and the optical axis.
Photonics 10 01381 g002
Figure 3. The optical system of lens I is an objective lens. Table 1 lists its optical parameters. The numbers in the figure represent the optical surfaces, where the fifth surface is the aperture stop.
Figure 3. The optical system of lens I is an objective lens. Table 1 lists its optical parameters. The numbers in the figure represent the optical surfaces, where the fifth surface is the aperture stop.
Photonics 10 01381 g003
Figure 4. The optical system of lens II is a fisheye lens. Table 2 lists its optical parameters. The numbers in the figure represent the optical surfaces, where the fifth surface is the aperture stop.
Figure 4. The optical system of lens II is a fisheye lens. Table 2 lists its optical parameters. The numbers in the figure represent the optical surfaces, where the fifth surface is the aperture stop.
Photonics 10 01381 g004
Figure 5. The diagram shows the wave aberration distribution of lens I at ω 0 = 20 ° . Table 1 lists its optical parameters.
Figure 5. The diagram shows the wave aberration distribution of lens I at ω 0 = 20 ° . Table 1 lists its optical parameters.
Photonics 10 01381 g005
Figure 6. The diagram shows the wave aberration distribution of fisheye lens II at ω 0 = 80 ° . Table 2 lists its optical parameters.
Figure 6. The diagram shows the wave aberration distribution of fisheye lens II at ω 0 = 80 ° . Table 2 lists its optical parameters.
Photonics 10 01381 g006
Figure 7. The abscissa indicates the aspheric coefficient a 2 ; the ordinate indicates the evaluation function Q. The curve shows the relationship between Q and a 2 of lens I.
Figure 7. The abscissa indicates the aspheric coefficient a 2 ; the ordinate indicates the evaluation function Q. The curve shows the relationship between Q and a 2 of lens I.
Photonics 10 01381 g007
Figure 8. The abscissa indicates the aspheric coefficient a 2 ; the ordinate indicates the evaluation function Q. The curve shows the relationship between Q and a 2 of fisheye lens II.
Figure 8. The abscissa indicates the aspheric coefficient a 2 ; the ordinate indicates the evaluation function Q. The curve shows the relationship between Q and a 2 of fisheye lens II.
Photonics 10 01381 g008
Figure 9. The diagram shows the wave aberration distributions of lens I at ω 0 = 20 ° . They are recalculated with Equations (17)–(22), and only the a 2 of the sensitive optical surface is changed. (a) shows a 2 2 = 1.69 and (b) shows a 2 3 = 0.79 , respectively.
Figure 9. The diagram shows the wave aberration distributions of lens I at ω 0 = 20 ° . They are recalculated with Equations (17)–(22), and only the a 2 of the sensitive optical surface is changed. (a) shows a 2 2 = 1.69 and (b) shows a 2 3 = 0.79 , respectively.
Photonics 10 01381 g009
Figure 10. The diagram shows the wave aberration distributions of fisheye lens II at ω 0 = 80 ° . (a) shows a 2 2 = 2.12 and (b) shows a 2 4 = 3.42 , respectively.
Figure 10. The diagram shows the wave aberration distributions of fisheye lens II at ω 0 = 80 ° . (a) shows a 2 2 = 2.12 and (b) shows a 2 4 = 3.42 , respectively.
Photonics 10 01381 g010aPhotonics 10 01381 g010b
Figure 11. The diagram shows the layout of fisheye lens III. Table 4 lists its optical parameters.
Figure 11. The diagram shows the layout of fisheye lens III. Table 4 lists its optical parameters.
Photonics 10 01381 g011
Figure 12. The abscissa indicates the aspheric coefficient a 2 ; the ordinate indicates the evaluation function Q. The curve shows the relationship between Q and a 2 of fisheye lens III.
Figure 12. The abscissa indicates the aspheric coefficient a 2 ; the ordinate indicates the evaluation function Q. The curve shows the relationship between Q and a 2 of fisheye lens III.
Photonics 10 01381 g012
Figure 13. The diagram shows the MTF curves of fisheye lens III. (a) shows the MTF curve of the original fisheye lens III; (b) shows the MTF curve of the optimized fisheye lens III, with a 2 ( 6 ) = 1.065 . Their optical parameters are listed in Table 5.
Figure 13. The diagram shows the MTF curves of fisheye lens III. (a) shows the MTF curve of the original fisheye lens III; (b) shows the MTF curve of the optimized fisheye lens III, with a 2 ( 6 ) = 1.065 . Their optical parameters are listed in Table 5.
Photonics 10 01381 g013
Table 1. Optical parameters of lens I.
Table 1. Optical parameters of lens I.
Surface iRadius/mmSpacing/mmIndexGlass
ObjectInfinite2000.0
135.05.01.84666N-SF57
215.01.0
315.012.01.45600N-FK58
4−25.05.0
5 (STO)Infinite10.0
6−130.03.01.68893P-SF8
7−30.015.0
Table 2. Optical parameters of fisheye lens II.
Table 2. Optical parameters of fisheye lens II.
Surface iRadius/mmSpacing/mmIndexGlass
ObjectInfinite2000.0
1141.4453.8081.71300N-LAK8
223.05917.392
387.8894.0991.71300N-LAK8
425.68541.513
5 (STO)Infinite2.0
630.04.01.71300N-LAK8
7−80.015.0
8−24.04.01.45600N-FK58
9−50.031.101--
Table 3. Optimal optical aspheric position in lens I and fisheye lens II, and the corresponding initial value of a 2 .
Table 3. Optimal optical aspheric position in lens I and fisheye lens II, and the corresponding initial value of a 2 .
Lens IFisheye Lens II
Surface2324
a 2 −1.690.79−2.12−3.42
Table 4. Optical parameters of the original and optimized fisheye lens III.
Table 4. Optical parameters of the original and optimized fisheye lens III.
Surface iRadius/mmSpacing/mmIndexGlass
OriginalOptimized OriginalOptimized
1162.236148.8527.6765.9461.5163BK7HT
269.16835.58228.97529.988--
398.31460.4434.1863.7161.6127SK4
442.47435.56013.05410.810--
564.48441.5204.2103.5211.6148SSK3
612.62313.64413.53613.741--
7−105.027−170.4401.2971.2191.4875N-FK5
818.95314.9574.5355.0501.7847SF56A
9−232.713−273.5251.1751.1181.7555P-LAF37
1075.79174.98913.86212.265--
11
(STO)
2.5282.843--
1257511.353039.161.3641.3961.7847SF56A
1315.42513.3478.5867.2031.7433N-LAF35
14−46.059−63.2380.1420.117--
15−259.655−316.8904.7514.0431.7555P-LAF37
16−37.371−36.43818.49316.484--
1738.373134.7782.5242.1691.7847SF56A
1824.43819.82919.82023.8871.6204N-SK16
19−53.609−43.55812.96312.006--
Table 5. Weight factors at different field angles of fisheye lens III.
Table 5. Weight factors at different field angles of fisheye lens III.
ω i ( ) 01020304050607080
Frist round ε i 111122222
Second round ε i 112244667
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fan, L.; Yan, K.; Qiao, G.; Lu, L.; Gao, S.; Zheng, H. Determination Position and Initial Value of Aspheric Surface for Fisheye Lens Design. Photonics 2023, 10, 1381. https://doi.org/10.3390/photonics10121381

AMA Style

Fan L, Yan K, Qiao G, Lu L, Gao S, Zheng H. Determination Position and Initial Value of Aspheric Surface for Fisheye Lens Design. Photonics. 2023; 10(12):1381. https://doi.org/10.3390/photonics10121381

Chicago/Turabian Style

Fan, Lirong, Ketao Yan, Guodong Qiao, Lijun Lu, Shuyuan Gao, and Huadong Zheng. 2023. "Determination Position and Initial Value of Aspheric Surface for Fisheye Lens Design" Photonics 10, no. 12: 1381. https://doi.org/10.3390/photonics10121381

APA Style

Fan, L., Yan, K., Qiao, G., Lu, L., Gao, S., & Zheng, H. (2023). Determination Position and Initial Value of Aspheric Surface for Fisheye Lens Design. Photonics, 10(12), 1381. https://doi.org/10.3390/photonics10121381

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop