The present invention relates to an eccentricity measuring method, a lens manufacturing method, and an eccentricity measuring apparatus.
An aspheric lens has recently been used for an optical apparatus, such as a camera and a semiconductor exposure apparatus, in order to reduce the size of its optical system. In addition, highly accurately manufacturing the aspheric lens is required for the high definition of an image obtained by the optical apparatus. Thus, a technique is demanded which measures the relative eccentricity (inter-surface eccentricity) of both surfaces of the aspheric lens with high accuracy.
Japanese Patent No. (“JP”) 4767255 discloses a method for obtaining an inter-surface eccentricity by holding an aspheric lens with a holder having a through-hole. This method initially scans a first surface of the aspheric lens and the through-hole and finds a positional relationship between them. Next, this method inverts the aspheric lens together with the holder, scans a second surface of the aspheric lens and the through-hole, and finds a positional relationship between them. Thereafter, the positional relationship between the first surface and the second surface or the inter-surface eccentricity (or decentering) is obtained with reference to the position of the through-hole.
However, the method of JP 4767255 needs to scan both surfaces of the aspheric lens by inverting the aspheric lens, and thus requires a long time to measure the inter-surface eccentricity.
The present invention provides an eccentricity measuring method, a lens manufacturing method, and an eccentricity measuring apparatus, each of which can measure inter-surface eccentricity of a test lens at high speed and with high accuracy.
An eccentricity measuring method according to one aspect of the present invention includes a forming step of dividing reflected light from a first surface and a second surface of a test lens by a plurality of optical elements and of forming a first spot group and a second spot group, and an eccentricity calculating step of calculating an eccentricity amount of the first surface relative to the second surface based on the first spot group and the second spot group. A lens manufacturing method according to another aspect of the present invention includes a producing step for producing a test lens, and an eccentricity measuring step for measuring an eccentricity amount of the test lens using the eccentricity measuring method.
An eccentricity measuring apparatus according to another aspect of the present invention includes a plurality of optical elements configured to divide reflected light from a first surface and a second surface of a test lens and to form a first spot group and a second spot group; and a calculator configured to calculate an eccentricity amount of the first surface relative to the second surface based on the first spot group and the second spot group.
Further features of the present invention will become apparent from the following description of exemplary embodiments with reference to the attached drawings.
Referring now to the accompanying drawings, a description will be given of embodiments according to the present invention. Corresponding elements in respective figures will be designated by the same reference numerals, and a duplicate description thereof will be omitted.
When the xyz orthogonal coordinate system is defined as illustrated in
Light from the light source 1 is emitted from the optical fiber connector 1b via the single mode optical fiber 1a, and travels as a spherical wave along a measurement optical axis 1c. The light source 1 uses a monochromatic laser in this embodiment, but may use a light emitting diode or the like.
The stage 7 moves the test lens 12 based on a command from the stage controller 7a. In the stage 7, similar to the test lens 12, the xyz orthogonal coordinate system illustrated in
The lens 4 converts the light from the optical fiber connector 1a that has transmitted through the half-mirror 8, into convergent light, and illuminates the first surface 12a of the test lens 12. The half-mirror 8 is disposed at a slope of 45° relative to the measurement optical axis 1c. Thereby, a traveling direction of the reflected light of the test lens 12 is converted from a direction parallel to the measurement optical axis 1c to a direction parallel to a measurement optical axis 1d orthogonal to the measurement optical axis 1c. The measurement optical axes 1c and 1d intersect each other on the half-mirror 8. The lens 5, together with the half-mirror 8 and the lens 4, constitutes an optical system that focuses the light has just been reflected by the first surface 12a, on the detector 9 with an imaging magnification m. In other words, the lens 5, the half-mirror 8, and the lens 4 serve as a imaging lens 14, and the first surface 12a and the detector 9 are conjugate with each other via the imaging lens 14.
The focal length and effective diameter of each of the lenses 4 and 5 are determined by the effective diameter and curvature radius of the first surface 12a and the size of the detection surface of the detector 9. The distance between the lens 4 and the test lens 12 is set so that the light that has transmitted through the lens 4 converges near the center of curvature of the first surface 12a. Thereby, the light that has transmitted through the lens 4 enters the first surface 12a almost perpendicularly over its entire surface, and the reflected light travels through the same optical path as the incident light up to the half-mirror 8. The imaging lens 14 may include a movable mechanism or an exchange mechanism so that the above condition is always satisfied for the test lens 12 having a different design shape. However, an angle of the reflected light of the first surface 12a depends on the aspheric amount (deviation from the spherical surface) and the shape error of the first surface 12a. When the aspherical amount of the first surface 12a is large, the angle of the reflected light on the first surface 12a shifts from the angle of the incident light.
In this embodiment, MLA 2 is made to coincide with the conjugate surface of the first surface 12a. However, other surfaces in the detector 9, such as the light receiving surface 3a, may coincide with the conjugate surface of the first surface 12a.
The position of the optical axis of ML 2a (ξ0,j,k, η0,j,k) is previously acquired, for example, by a method disclosed in Applied Optics Vol. 44, no. 30, p 6419. The detector 9 is assembled so that the distance 1 is equal to the focal length f, but the assembly has a finite error. The distance 1 is also calibrated to acquire an accurate value in advance. In this embodiment, the distance between all the MLs 2a and the light receiving surface 3a is treated as being uniform with the surface interval 1, but the distances 1j, k relative to the light receiving surface 3a are obtained for each ML 2a and reflected in a calculation expression for the back surface inclination described later.
The processor 10 includes a CPU 10a having a calculation function and serves as a calculator. The processor 10 can input an output signal of the image sensor 3 via an unillustrated interface, and output position control information of the test lens 12 to the stage controller 7a. Using these functions, the processor 10 performs measurement processing in accordance with the measurement procedure of the inter-surface eccentricity of the lens 12 to be described later. The measurement processing needs a program for the measurement processing, information on the configuration of the detector 9, information on the shapes and arrangement of the lenses 4 and 5, the half-mirror 8 and the test lens 12. These data are stored in a predetermined area of a memory such as a ROM 10b and a RAM 10c provided in the processor 10, for example.
The processor 10 includes a communicator 10d includes a communicator 10d, such as a network interface of IEEE802.3 standard, and is connected to a display unit 10e. The CPU 10a displays, for example, the measurement result of the inter-surface eccentricity of the test lens 12 or the evaluation result of the test lens 12 based on the measurement result on the display unit 10e, or transmits via the communicator 10d to another equipment in a lens manufacturing plant in which the eccentricity measuring apparatus 100 is installed.
When the test lens 12 is installed in the holder 7b, part of the light converged by the lens 4 is reflected by the first surface 12a. The wavefront of the light that has just been reflected, reflects the aspherical shape of the first surface 12a. The light reflected by the first surface 12a (first light) passes through the lens 4, is reflected by the half-mirror 8, is approximately collimated by the lens 5, and enters the detector 9. Part of the light converged by the lens 4 passes through the first surface 12a and is reflected by the second surface 12b. The light reflected by the second surface 12b (second light) once converges near the first surface 12a and then becomes divergent light, is approximately collimated by the lens 4, and enters the lens 5. However, since the second light entering the lens 5 has a light beam diameter that greatly exceeds the aperture diameter of the lens 5, most of the second light is blocked by the holder 5a that holds the lens 5. Only the second light near the measurement optical axis 1d that is not blocked by the holder 5a is converged by the lens 5 and enters the central portion of the detector 9.
Assume that (ξu, j, k, ηu, j, k) is a position of the spot formed by the first light L1 passing through “ML on the k-th row and j-th column,” and (ξb, j, k, ηb, j, k) is a position of the spot formed by the second light L2 passing through it. When the “ML on the k-th row and the j-th column” is near the measurement optical axis 1d, these spot positions are approximated by the following expressions (1) to (4).
ξu,j,k=ξu,0,0+jqu (1)
ηu,j,k=ηu,0,0+kqu (2)
ξb,j,k=ξb,0,0+jqb (3)
ηb,j,k=ηb,0,0+kqb (4)
qu is the spot interval of the spot group SP1, and qb is the spot interval of the spot group SP2. (ξu, 0, 0, ηu, 0, 0) is the position of the spot formed by the first light L1 that has transmitted through “ML on 0-th row and 0-th column,” and the shift amount of the spot group SP1 (first shift amount). (ξb, 0, 0, ηb, 0, 0) is the position of the spot formed by the second light L2 that has passed through “ML on the 0-th row and 0-th column,” and represents the shift amount (second shift amount) of the spot group SP2.
In the step S101, the processor 10 calculates a relevant parameter group (spot parameter group) from a device set value for a spot that may be imaged by the image sensor 3.
Herein, ρu and ρb represent the curvatures of the wavefronts W1 and W2 incident on the detector 9, respectively. For p, l, ξ0, 0, 0, and η0, 0, 0, previously acquired values may be substituted. ρu and ρb may be analytically calculated from the designed value of the curvature of the light wavefront incident on the test lens 12, the designed value of the imaging magnification of an imaging lens 14, and the designed shape and refractive index of the test lens 12. Alternatively, the ray tracing may be performed based on these parameters.
In the step S101b, the processor 10 calculates a parameter group W (=(wu, wb, Iu, Ib)). wu and wb represent the spot radii constituting the spot groups SP1 and SP2, respectively. Iu and Ib represent the peak intensities of the spot groups SP1 and SP2, respectively. Both parameters relate to the spot shape. The spot radii wu and wb are calculated using a light propagation expression such as the Fresnel diffraction or Fraunhofer diffraction from the designed value of ML 2a and the surface interval 1 between MLA 2 and the light receiving surface 3a. Instead of using the light propagation formula, the light propagation may be calculated by the angular spectrum method, the FDTD method, or the like, and the shape of the wavefront W1 of the first light L1 or the wavefront W2 of the second light L2 may be considered in the calculation. The peak intensities Iu and Ib are calculated from the irradiation densities of the first light L1 and the second light L2 entering the detector 9, the spot radii wu and wb, and the light receiving sensitivity of the image sensor 3. The irradiation densities of the first light L1 and the second light L2 may be calculated analytically based on the output of the light source 1, the imaging magnification of the imaging lens 14, or the like, or may be simulated by ray tracing.
In the step S101c, the processor 10 calculates a parameter group R (=ru, rb, rap, ra)). As illustrated in
The processing in the step S101 thus ends. The parameter groups X, W, and R calculated in the step S101 are stored in the ROM 10b (or RAM 10c).
The processing in the step S101 may be performed in advance to the measurement of the inter-surface eccentricity. The parameter group is calculated based on the device design value in this embodiment, but may be actually measured on the eccentricity measuring apparatus 100.
In the step S102, the test lens 12 is placed on the holder 7b. At that time, the aspherical axis 12c is made approximately to coincide with the measurement optical axis 1c. Hence, a mark may be previously provided on the holder 7b, and the test lens 12 may be installed in accordance with the mark, the holder 7b may include an unillustrated positioning pin which the test lens 12 contacts. The position of the stage 7 in the θz direction is set to 0°. As a result of installing the test lens 12 on the holder 7b, the spot groups SP1 and SP2 are formed on the light receiving surface 3a as illustrated in
In the step S103a, the alignment of the test lens 12 is made in the x direction, the y direction, the θx direction, and the θy direction.
In the step S103b, the alignment of the test lens 12 is made in the z direction.
In the step S104, the processor 10 causes the image sensor 3 to capture a spot image based on the spot groups SP1 and SP2, and acquires as a spot image the light intensity distribution I(ξ, η) on the light receiving surface 3a from the image sensor 3. The processing in the step S104 corresponds to the imaging step.
The processing after the step S105 below corresponds to the eccentricity calculation step.
In the step S105, the processor 10 calculates a provisional back surface inclination (Δθl, x, Δθl, y) from the spot image acquired in the step S104.
Among the three steps constituting the step S105, the second step S105b for fitting the spot image will be described first. This fitting calculates a parameter group X representing the spot position that minimizes the difference between the acquired spot image and the spot image calculated by the following expression (11).
The calculation may use a nonlinear optimization program such as a solver function of Excel software or an fminsearch function of Matlab software. The spot image to be fitted is limited to the analysis area AR1 having the radius ra calculated in the step S101c.
A description will be given of the expression (11) used for the spot image modeling in the fitting. In the expression (11), the intensity distribution of each spot is expressed by a Gaussian function, and the spot image is expressed as a superposition thereof. As long as the intensity distribution of the spot is appropriately expressed, it is not limited to the Gaussian function, and a Bessel function, a sinc function, or the like may be used. It is assumed that the spots are arranged at regular intervals in both the spot groups SP1 and SP2 and follow the expressions (1) to (4). This premise is based on the fact that the analysis area AR1 is limited to the central portion of the light receiving surface 3a. As a result, the first term on the right side representing the spot group SP1 is a two-dimensional periodic function having a period qu and a phase (2πξu, 0, 0/qu, 2πηu, 0, 0/qu). The second term representing the spot group SP2 is a two-dimensional periodic function having a period qb and a phase (2πξb, 0, 0/qb, 2πηb, 0, 0/qb). In other words, the entire right side is the sum of periodic functions. The shift amount of the spot group is proportional to the phase of the periodic function and corresponds one-to-one.
A description will now be given of a relationship between the back surface inclination of the test lens 12 and the expression (11). As described above, the back surface inclination is defined as the difference between the inclination of the first surface 12a and the inclination of the second surface 12b. As the first surface 12a inclines, the wavefront W1 of the first light L1 is inclined and the spot group SP1 is shifted. As the second surface 12b inclines, the wavefront W2 of the second light L2 is inclined and the spot group SP2 is shifted. As a result, the back surface inclination is reflected on a shift between the shift amount (ξu, 0, 0, ηu, 0, 0) of the spot group SP1 and the shift amount (ξb, 0, 0, ηb, 0, 0) of the spot group SP2 or the phase shift between the first term on the right side and the second term on the right side in the expression (11). The parameter group X includes a parameter representing these shift amounts (phases), which is equivalent to information on the back surface inclination. This embodiment addresses this point, and this is the reason why the parameter group X is calculated in the step S105b.
In order to correctly fit the spot image using the expression (11) in the step S105b, it is necessary to input the highly accurate parameter group W. The step S101b calculates the parameter group W, which is calculated from the device design value, and does not reflect the manufacturing error of the eccentricity measuring apparatus 100. Hence, the parameter group W calculated in the step S101b is insufficient in accuracy to be used in the step S105b.
In fitting a spot image in a wide area with a nonlinear optimization program, high accuracy is required for initial values of the spot intervals qu and qb constituting the parameter group X. If the initial values of the spot intervals qu and qb, are deviated from the actual spot interval, the fitting is started from a state in which the positions of the spots on the outer periphery are largely deviated, and the convergence lowers. In the step S101a, the spot intervals qu and qb are calculated, but these are calculated from the device design values, and the manufacturing error of the eccentricity measuring apparatus 100 is not reflected. Hence, the accuracies of the spot intervals qu and qb calculated in the step S101a are insufficient to be used in the step S105b.
Accordingly, this embodiment performs the processing of the step S105a for pre-fitting the spot image before the processing of the step S105b. In the step S105a, pre-fitting is performed for the spot image in the analysis area AR2 smaller than the analysis area AR1 using the expression (11), and the parameter groups X and W are calculated. The values calculated in the step S101 are used as initial values of the parameter groups X and W in the pre-fitting. Since the pre-fitting analysis area AR2 is sufficiently small, a high convergence can be obtained even if the parameter groups X and W based on the device design values are set as initial values. In the fitting in the step S105b, the parameter group W is fixed to the value calculated in the step S105a, and the parameter group X calculated in the step S105a is set as an initial value. By introducing such pre-fitting processing in advance, highly accurate fitting can be realized in the step S105b.
The pre-fitting in the step S105a uses the expression (11), and the fitting initial values of the spot shift amounts (ξu, 0, 0, ηu, 0, 0) and (ξb, 0, 0, ηb, 0, 0) use the results of the expressions (7) to (10). As described above, the expressions (7) to (10) are based on the assumption that the optical axis of the wavefront W1 coincides with the measurement optical axis 1d. In other words, the spot image to be pre-fitted needs to be acquired in a state where the optical axis of the wavefront W1 coincides with the measurement optical axis 1d. Thus, the aspherical axis 12c needs to coincide with the measurement optical axis 1c. If these axes do not coincide with each other, the wavefronts of the light reflected by the first surface 12a and the second surface 12b significantly incline, and the wavefronts W1 and W2 are also significantly incline. As a result, the shift amounts of the spot groups SP1 and SP2 significantly shift, and stable convergence cannot be expected even when fitting is performed using the shift amount calculated by the expression (11) as an initial value. The step S102 makes the aspherical axis 12c substantially coincide with the measurement optical axis 1c using a contact mechanism or the like, but this mechanism alone cannot provide the installation accuracy necessary for pre-fitting.
Accordingly, in this embodiment, before the spot image is captured in the step S104, the step S103 makes an alignment of the test lens 12 in the x direction, the y direction, the θx direction, and the θy direction. More specifically, similar to the step S104, a spot image is captured first. As described above, the first light L1 and the second light L2 have entered the processor 10, and the light beam diameter of the first light L1 is larger. Thus, the light receiving surface 3a has a donut-shaped area where only the spot group SP1 formed by the first light L1 is formed. Since the spot group SP2 does not exist in this area, the position of the spot group SP1 can be detected according to the conventional simple method. Accordingly, this embodiment detects the position of the spot group SP1 in the donut-shaped area through the processor 10 according to the method described in Japanese Patent Laid-Open No. 2016-38300, and calculates an inclination component and a coma aberration component of the wavefront W1 of the first light L1 detected therefrom. Moreover, following the method described in Japanese Patent Laid-Open No. 2015-75396, etc., this embodiment calculates positional shifts of the aspherical axis 12c in the x direction, the y direction, the θx direction, and the θy direction relative to the measurement optical axis 1c from the inclination component and the coma aberration component of the wavefront W1. In other words, the aspherical axis 12c is calculated from the spot group SP1. Then, the processor 10 issues a command to the stage controller 7a to drive the stage 7 so as to cancel the calculated positional shift amount. This embodiment makes an alignment by referring to the inclination component of the wavefront W1. Since the inclination component of the wavefront W1 is obtained from the shift amount of the spot group SP1, the alignment is made based on the shift amount of the spot group SP1 in the step S103a.
The step S103a makes the alignment of the test lens 12, and the aspherical axis 12c coincides with the measurement optical axis 1c. Thereby, the shift amount of the spot group SP1 to be imaged substantially matches the value calculated by the expressions (7) to (8), and the shift amount of the spot group SP2 approaches to the value calculated by the expressions (9) to (10). As a result, the convergence of pre-fitting in the step S105a is improved.
In the spot S105c, the processor 10 calculates a provisional back surface inclination. Now assume that the aspherical axis 12c coincides with the measurement optical axis 1c, the lenses 4 and 5 are accurately axisymmetric relative to the measurement optical axes 1c and 1d, respectively, and the back surface inclination has a finite value. At this time, since the first surface 12a is not inclined to the measurement optical axis 1c, the inclination (φξ, b, φη, b) of the wavefront W2 incident on the processor 10 is expressed by the following expressions (12) and (13) using the back surface inclination (Δθl, x, Δθl, Y), the imaging magnification m of the imaging lens 14, and the refractive index n of the test lens 12.
On the other hand, the shift amount (ξb, 0, 0, ηb, 0, 0) of the spot group SP2 has the relationship with the inclination (φξ, b, φη, b) of the wavefront W2 as expressed in the following expressions (14) and (15).
From the following expressions (12) to (15), the back surface inclination of the test lens 12 is calculated by the following expressions (16) and (17).
The step S105c calculates the back surface inclination using the expressions (16) and (17) and the shift amount (ξb, 0, 0, ηb, 0, 0) calculated in the step S105b and the spot interval qb, and sets it to the provisional back surface (Δθ′l, x(θz=0), Δθ′l, y (θz=0)).
This embodiment detects the shift amount of the spot group SP2 by coinciding the shift amount of the spot group SP1 with the position specified by the expressions (7) to (10) through the alignment, and calculates the back surface inclination by substituting the result for the expression (16) and (17). This calculation corresponds to calculating the back surface inclination based on the shift amounts of both the spot groups SP1 and SP2 or calculating the back surface inclination based on a shift between the phase of the periodic function formed by the spot group SP1 and the phase of the periodic function formed by the spot group SP2.
The position of the detector 9 may be adjusted so as to satisfy ξ0, 0, 0=η0, 0, 0=0 before the measurement of the inter-surface eccentricity is performed. At this time, from the simplified expressions (16) and (17) (Δθl, x=mηb, 0, 0/(2nl) and Δθl, y=−mξb, 0, 0/(2nl)), the back surface inclination can be calculated without using the spot interval qu and qb. In addition, when an extremely accurate measurement value of the inter-surface eccentricity is unnecessary, the provisional back surface inclination may be calculated from the parameter group X acquired by the pre-fitting in the step S105a without performing the step S105b.
The imaging magnification of the imaging lens 14 is m, which is a designed value, when the first surface 12a is installed on the conjugate surface of the detector 9. However, when the first surface 12a shifts from the conjugate surface, the imaging magnification shifts from m, and an error occurs in the back surface inclination calculated by the expressions (16) and (17). Accordingly, in order to coincide the imaging magnification of the imaging lens 14 with m, the step S103b is provided which makes an alignment of the test lens 12 in the z direction and coincides the first surface 12a with the conjugate plane before the step S104 for capturing the spot image is processed. The step S103b moves the stage 7 in the z direction so as to cancel the positional shift while monitoring the z position shift of the test lens 12 through the length measuring unit 15. At this time, the z position measurement result by the length measuring unit 15 is input into the processor 10, and the stage driving amount is transmitted from the processor 10 to the stage controller 7a based on the result. The length measuring unit 15 is provided with an unillustrated movable mechanism and usually installed away from the optical path, but is inserted between the lens 4 and the test lens 12 in the step S103b. The length measuring unit 15 may use a length measuring unit by a triangulation, a laser length measuring unit, or a white interferometer.
The expressions (12) and (13) assume that the lenses 4 and 5 are accurately axisymmetric with respect to the measurement optical axes 1c and 1d, respectively, and the expressions (16) and (17) derived from them are based on the same assumption. However, the lenses 4 and 5 actually have manufacturing errors and installation errors, so they are not strictly axially symmetric. Hence, a systematic error derived from the non-axisymmetric property of the lenses 4 and 5 occurs in the back surface inclination calculated by the expressions (16) and (17). This embodiment extracts and corrects a systematic error derived from the non-axisymmetric property of the lenses 4 and 5 from the measured provisional back surface inclination. As described above, since the test lens 12 can rotate around the optical axis of the test lens 12, this embodiment corrects the systematic error by calculating the eccentricity data when the test lens 12 is located at a first rotating position and a second rotating position.
In the step S106, the processor 10 determines whether the position of the stage 7 in the θz direction is 180°. When the position of the stage 7 in the θz direction is 180°, the flow proceeds to the step S108, and when the position of the stage 7 in the θz direction is not 180°, the flow proceeds to step S107.
In the step S107, the processor 10 rotates the stage 7 by 180° in the θz direction. Thereafter, the processing of the steps S103a, S103b, S104, and S105 are performed again, and a provisional back surface inclination (Δθ′l, x (θz=π), Δθ′l, y (θz=π)) is calculated.
In the step S108, the processor 10 calculates the systematic error (Δθs, x, Δθs, y) included in the provisional back surface inclination using the following expressions (18) and (19).
The expressions (18) and (19) are derived from the fact that the provisional back surface inclination is expressed by the following expressions (20) to (23).
Δθ′l,x(0)=Δθl,x+Δθs,x (20)
Δθ′l,y(0)=Δθl,y+Δθs,y (21)
Δθ′l,x(π)=−Δθl,x+Δθs,x (22)
Δθ′l,y(π)=−Δθl,y+Δθs,y (23)
In the step S109, the processor 10 calculates the back surface inclination by removing the systematic error (Δθx, x, Δθx, y) from the provisional back surface inclination (Δθ′l, x(0), Δθ′l, y(0)). In other words, the systematic error included in the provisional back surface inclination (eccentric data) is corrected. Calculating the systematic error and acquiring the provisional back surface inclination (Δθ′(π)l, x, Δθ′(π)l, y) necessary for the calculation are performed only by the initial measurement after the assembly of the eccentricity measuring apparatus 100, and in the step S109 for the second and subsequent times, the provisional back surface inclination may be corrected with the systematic error acquired for the first time. The back surface inclination may be directly calculated by deleting Δθs, x and Δθs, y from the expressions (20) to (23) without performing the step S108.
In the step S110, the processor 10 multiplies the back surface inclination acquired in the step S109 by a designed value R of the curvature radius of the second surface 12b, and converts it into the optical axis shift amounts Δl, x and Δl, y of the second surface 12b. Moreover, The processor 10 outputs the optical axis shift amounts Δl, x and Δl, y as the eccentric data. In outputting the eccentricity data, it may be displayed on the display unit 10e or may be transmitted to another device via the communicator 10d. Moreover, the back surface inclination acquired in the step S109 may be directly output as eccentricity data. Whether or not the test lens 12 can be shipped may be determined based on the output eccentric data, or the test lens 12 may be further processed to suppress the eccentricity.
As described above, the configuration of the processor 10 is similar to that of a general SHS. SHS is known as a wavefront sensor, and measures a light beam having a gentle single wavefront (equiphase surface). On the other hand, in this embodiment, the measurement target of the processor 10 is a light beam having two wavefronts.
On the other hand, as an optical measurement device configured to measure a light flux not composed of a single wavefront, the eccentricity measurement using an autocollimator is known. This method collects the light reflected on the first surface and the light reflected on the second surface of the test lens, through a common lens, and detects the position of each one spot (totally two spots) formed on the image sensor. Thereafter, the inclination difference (eccentricity) between the first surface and the second surface is obtained from these spot positions. However, this method cannot detect the aspherical axis of the aspheric lens, or cannot measure the back surface inclination based on this result.
The SHS has a huge image sensor and a plurality of microlenses, and can be regarded as “a plurality of autocollimators arranged in parallel.” The SHS can detect an aspherical axis that cannot be detected by a single autocollimator. This embodiment considers this point, and adds a function of analyzing two spot groups to the SHS conventionally known as a (single) wavefront sensor and applies it to the measurement of a light flux that does not form a single wavefront. Thereby, this method can measure the back surface inclination based on the aspherical axis that cannot be measured by the conventional autocollimator.
In general, the light receiving surface of the image sensor has pixels divided into about 10 μm. On the other hand, as described above, the processor 10 can be replaced with general SHS, and the general SHS spot diameter is about 40 μm. In other words, the pixel size is an about quarter of the spot size and is not sufficiently small. Under such circumstances, it is not obvious that the shift amount between the two spot groups can be calculated with high accuracy by the fitting in the step S105b.
A simulation is performed on the relationship between the pixel size of the image sensor and the calculation error of the shift amount of the spot group.
In the simulation, the one-dimensional spot group SP1 and the spot group SP2 are synthesized and discretized with a predetermined pixel size to calculate a spot image.
Since the accurate shift amount of each spot group can be acquired by fitting, this embodiment can acquire highly accurate inter-surface eccentricity. Since it is unnecessary to scan the measurement probe or the test lens, the inter-surface eccentricity can be quickly acquired.
This embodiment is the same as the first embodiment in using the eccentricity measuring apparatus 100 in
The first embodiment calculates the provisional back surface inclination using the expressions (16) and (17). The expressions (16) and (17) assume that the aspherical axis 12c completely coincides with the measurement optical axis 1c in consideration of the alignment of the test lens 12 in the step S103a. However, there is a finite error in the wavefront measurement result serving as the reference of the alignment, and there is also the finite installation error in the test lens 12 after the alignment. This embodiment calculates the provisional back surface inclination in consideration of the installation error of the test lens 12.
When the aspherical axis 12c inclines to the measurement optical axis 1c by θx, u in the θx direction and θy, u in the θy direction, the inclination of the wavefront W1 of the first light L1 entering the detector 9 (φξ, u, φη, u) is expressed by the following expressions (24) and (25).
When the second surface 12b inclines in the θx direction by θx, b, and θy, b in the θy direction, the inclination (φξ, b, φη, b) of the wavefront W2 of the second light L2 entering the detector 9 is expressed by the following expressions (26) and (27).
The back surface inclination (Δθl, x, Δθl, y), which is the difference in inclination between the first surface 12a and the second surface 12b, is calculated by the following expressions (28) and (29).
The step S105c calculates the provisional back surface inclination using the expressions (28) and (29).
This embodiment acquires the inter-surface eccentricity in consideration of the alignment error generated in the step S103a, and can realize a measurement with higher accuracy in comparison with the first embodiment.
This embodiment is the same as the first embodiment in using the eccentricity measuring apparatus 100 in
Since the processing of the steps S301 to S305 is the same as the processing of the steps S101 to S105 in
In this embodiment, in comparison with the first embodiment, the number of times of calculating the provisional back surface inclination is changed from 2 times to N times (N≥3).
In the step S306, the processor 10 determines whether or not the position of the stage 7 in the θz direction is 2π(N−1)/N. When the position of the stage 7 in the θz direction is 2π(N−1)/N, the flow proceeds to the step S308, and when the position of the stage 7 in the θz direction is not 2π(N−1)/N, the flow proceeds to the step S307.
In the step S307, the processor 10 rotates the stage 7 by 2π/N in the θz direction. Thereby, N sets of provisional back surface inclinations (Δθ′l, x (θz=2πi/N), Δθ′l, y (θz=2πi/N)) (i=0, 1, . . . , (N−1)) is obtained. In this embodiment, the stage 7 is rotated in the θz direction by 2π/N each time, but if the rotation angle is known, it is not always necessary to have the same angle every time.
In the step S308, the processor 10 calculates the systematic error of the back surface inclination, but the calculation method is different from that according to the first embodiment. This embodiment fits the obtained provisional back surface inclination value by the following expressions (30) and (31).
Δθ′l,x(θz)=Δθs,x+Δθl cos(θz−θz,0) (30)
Δθ′l,y(θz)=Δθs,y+Δθl sin(θz−θz,0) (31)
The systematic error (Δθs, x, Δθx, y) is calculated by the expressions (30) and (31).
This embodiment calculates the systematic error based on a lot of provisional back surface inclination data in comparison with the first embodiment. Hence, the systematic error can be calculated with higher accuracy and the provisional back surface inclination can be calculated with higher accuracy.
An experiment that measures the back surface inclination of the aspheric lens whose second surface is a flat surface based on this embodiment can show an error of 0.1 minutes or less. This experiment also suggests the effectiveness of this embodiment.
This embodiment is the same as the first embodiment in using the eccentricity measuring apparatus 100 in
The step S105b in the first embodiment assumes that MLs 2a are arranged in a square matrix in the ξ and η directions. However, in reality, there are a manufacturing error of MLA 2 and an installation error in the rotational direction around the optical axis of MLA 2, and MLs 2a are not strictly arranged in the square matrix. Therefore, the step S105b in this embodiment fits the spot image using, for example, the following the expression (32) in consideration of these errors.
As described above, the position (ξ0, j, k, η0, j, k) of each ML 2a has been calibrated, and this calibration value is substituted for the expression (32) in the fitting.
In this embodiment, the spot image is fitted in consideration of the manufacturing error and arrangement error of MLA 2. Thus, the fitting accuracy is more improved than the first embodiment, and a more accurate back surface inclination value can be obtained.
The measurement method according to this embodiment assumes that the MLs are designed to be arranged in the square matrix and intends to cope with the manufacturing error and the arrangement error, but is applicable to randomly arranged MLs.
This embodiment is the same as the first embodiment in using the eccentricity measuring apparatus 100 in
More specifically, the following expression (33) is used instead of the expression (11) when the spot image is pre-fitted in the step S105a constituting the step S105 or the spot image is fitted in the step S105b.
Thereby, the inclination (φξ, u, φη, u) of the wavefront W1 of the first light L1 and the inclination (φξ, b, φη, b) of the wavefront W2 of the second light L2 entering the detector 9 are acquired. In the step S105c, the acquired wavefront inclination is substituted for the following expressions (34) and (35) to calculate the provisional back surface inclination.
The expressions (34) and (35) do not explicitly include the shift amounts of the spot groups SP1 and SP2 when the provisional back surface inclination (eccentric data) is calculated in this embodiment. However, they include the inclination (φξ, u, φη, u) of the wavefront W1 and the inclination (φξ, b, φη, b) of the wavefront W2, and these are in one-to-one correspondence with the shift amounts of the spot groups SP1 and SP2. In other words, this embodiment also corresponds to the calculation of the inter-surface eccentricity based on the shift amounts of the spot groups SP1 and SP2.
In the first embodiment, when the first light L1 and the second light L2 enter the processor 10, the wavefronts W1 and W2 can be approximated by spherical waves only near the measurement optical axis 1d. Under this premise, for a model of the expressions (1) to (4) in which the spot groups SP1 and SP2 are both arranged at regular intervals, a spot image at the center is fitted. However, when the aspheric amount of the first surface 12a is large or when it is necessary to improve the accuracy by enlarging the fitting area, it is necessary to treat the wavefronts W1 and W2 as aspherical waves and the model of the first embodiment is not satisfied. This embodiment will describe a method for measuring the eccentricity of the lens 12 with high accuracy even in this case.
When the aspherical amount of the first surface 12a increases, the wavefronts W1 and W2 have axially symmetrical shapes ΔWu(ξ, η) and ΔWb(ξ, η) represented, for example, by the following expressions (36) and (37), respectively, H) is added.
cu, i and cb, i represent Zernike coefficients, and Zi(ξ, η) represents a Zernike function. The Zernike function Zi(ξ, η) is defined by the following expressions (38) and (39).
Z
9(ξ,η)=(1−6ρ2+6ρ4) (38)
Z
16(ξ,η)=−1+12ρ2−30 ρ4+20ρ6 (39)
Herein, ρ2=ξ2+η2. This embodiment sets the maximum number of terms to 16 when the wavefront is developed by the Zernike function, but may add a higher order axially symmetric component. The spot image at this time is represented by the following expression (40) instead of the expression (11).
In this embodiment, the step S105b uses the Zernike coefficient and the parameter group X as variation parameters to fit the spot image with the calculated value calculated from the expression (40). Other procedures follow the first embodiment.
According to this embodiment, the spot images are modeled with the wavefronts W1 and W2 as aspherical waves and thus can be fitted with high accuracy even when the aspheric amount of the first surface 12a is large. As a result, this embodiment can obtain a highly accurate inter-surface eccentricity as compared with the first embodiment.
This embodiment is the same as the first embodiment in using the eccentricity measuring apparatus 100 of
When a spot image is taken as the sum of periodic functions, it can also be expressed by the following expression (41).
Herein, ku and kb are wave numbers of a periodic function and correspond to the spot intervals qu and qb. (ψu, ξ, ψb, ξ) and (ψu, η, ψb, η) represent phases of the periodic function, and correspond to shift amounts (ξu, 0, 0, ηu, 0, 0), (ξb, 0, 0, ηb, 0, 0) of the spot group. When a two-dimensional complex Fourier transform is performed for this spot image, a Fourier transform image of the following expression (42) is obtained.
I(kξ,kη)∝[ΣmIu,m exp(−imφu,ξ)Δ(kξ−mku)][ΣnIu,n exp(−inφu,η)Δ(kη−nku)]+[ΣmIb,m exp(−imφb,ξ)Δ(kξ−mkb)][ΣnIb,n exp(−inφb,η)Δ(kη−nkb)] (42)
A peak is periodically formed in an absolute value map of the complex Fourier transform image. Among them, for the peak of m=n=1, information on the wave number (spot interval) is included at that position, and information on the phase (shift amount of the spot group) is included in an argument at that position. In other words, the information on the parameter group X required to calculate the back surface inclination is included.
Accordingly, in this embodiment, the step S105 calculates the provisional back surface inclination using a two-dimensional complex Fourier transform.
In the step S701, the processor 10 performs a two-dimensional complex Fourier transform for the spot image.
In the step S702, the processor 10 calculates the parameter group X from the Fourier transform image. In selecting the peak of m=n=1, the spot intervals qu and qb calculated from the device design value in the step S101a may be referred to.
In the step S703, the processor 10 calculates the provisional back surface inclination.
This embodiment requires no fitting work, and can more quickly and more stably calculate the back surface inclination than the first embodiment.
The present invention is applicable not only to an aspheric lens but also to a spherical lens. This embodiment will describe a test lens as the aspheric lens and the measurement of the inter-surface eccentricity between the surfaces.
This embodiment is the same as the first embodiment in using the eccentricity measuring apparatus 100 in
In the step S102, the test lens 13 is installed so that the reference surface 13c is parallel to the xy plane. In this embodiment, the holder 7b has an unillustrated reference surface for receiving the reference surface 13c, and the tilt of the stage 7 is adjusted so that the reference surface is parallel to the xy plane.
In the step S103a for making an alignment of the test lens 13, the stage 7 is adjusted in the x direction and the y direction so as to suppress the tilt component of the wavefront W1 of the first light L1. Unlike the first embodiment, the aspherical axis 12c is not made to coincide with the measurement optical axis 1c, and the wavefront W1 is not adjusted in the θx direction and the θy direction while the wavefront W1 is monitored.
Since there is no step of coinciding the aspherical axis with the measurement optical axis, this embodiment can measure the inter-surface eccentricity between the spherical lens.
This embodiment measures the inter-surface eccentricity according to the method described in the first embodiment, and then calculates the shape of the first surface 12a. This processing corresponds to the shape calculating step.
In calculating the shape of the first surface 12a, the stage 7 is set to θz=0, and the spot image I(ξ, η) and the parameter group X acquired in the steps S104 and S105b, respectively, are substituted for the following expression (43). Thereby, a spot image I′(ξ, η) is calculated in which only the spot group SP1 is extracted and the spot group SP2 caused by the second light L2 is eliminated.
Thereafter, the shape of the first surface 12a is calculated from the spot image I′(ξ, η) according to the method described in Japanese Patent Laid-Open No. 2013-186024.
This embodiment is the same as the first embodiment in using the eccentricity measuring apparatus 100 illustrated in
In the step S400, the processor 10 acquires a first temporary spot image. Details of the first temporary spot image and its acquisition procedure will be described later.
The steps S401 to S404 are performed, and a spot image is acquired. Since the processing of the steps S401 to S404 is the same as the processing of the steps S101 to S104 in
In the step S405, a provisional back surface inclination is calculated from the spot image and the first temporary spot image. Details will be described later.
The processing of the steps S406 to S410 is executed, and eccentricity data is calculated. Since the processing of the steps S406 to S410 is the same as that of the steps S306 to S310 in
An unillustrated spare lens 16 manufactured based on the same designed shape as that of the test lens 12 is used to obtain the first temporary spot image. The spare lens 16 has a third surface and a fourth surface based on the same designed shape as that of each of the first surface 12a and the second surface 12b of the test lens. The fourth surface of the spare lens 16 has received an antireflection treatment. Examples of the antireflection treatment include a formation of an antireflection film and roughening by sandblasting. Alternatively, a gel material may be applied which has a refractive index close to that of the test lens.
The step S405f subtracts the light intensity distribution (simply referred to as a first temporary spot image hereinafter) constituting a first temporary spot image ISP3(ξ, η) (shown as one temporary spot image) from the light intensity distribution (simply referred to as a spot image hereinafter) I(ξ, η) constituting the spot image acquired in the step S404. The subtraction result is a first differential spot image ΔI1(ξ, η) (=I−ISP3).
The spot image I(ξ, η) is a sum (I=ISP1+ISP2) of the light intensity distribution ISP1(ξ, η) constituting the spot group SP1 formed by the first light and the intensity distribution ISP2(ξ, η) constituting the spot group SP2 formed by the second light. On the other hand, the first temporary spot image ISP3 is constituted by the spot group SP3 formed by the third light of the spare lens 16, and is similar to the light intensity distribution ISP1. Thus, in the first differential spot image ΔI1 (=I−ISP3) that is the difference between the spot image I(ξ, η) and the one temporary spot image ISP3, the light intensity distribution ISP2 derived from the spot group SP2 remains as it is, and the light intensity distribution ISP1 derived from the spot group SP1 is reduced.
Herein, processing for further reducing the light intensity distribution derived from the spot group SP1 included in the first differential spot image ΔI1 may be added. For example, in the first differential spot image ΔI1, a signal smaller than a predetermined threshold or a negative signal may be converted into zero.
In the step S405g, the first differential spot image ΔI1 is analyzed, and the position (ξb, j, k, ηb, j, k) of each spot constituting the spot group SP2 is detected. Since the optical signal intensity derived from the spot group SP1 is reduced in the first differential spot image ΔI1, the position of each spot constituting the spot group SP2 can be detected according to the same method as that in the step S103a in the first embodiment.
In the step S405h, based on the spot position (ξb, j, k, ηb, j, k) detected in the step S405g, the light intensity distribution I2(ξ, η) (simply referred to as a second temporary spot image hereinafter) constituting the second temporary spot image is calculated. More specifically, for example, the second temporary spot image I2 is calculated by substituting the spot position (ξb, j, k, ηb, j, k) for the expression (44). Thereby, it is possible to obtain the second temporary spot image I2 including only the light intensity distribution derived from the spot group SP2 and not including the light intensity distribution derived from the spot group SP1.
The step S405i subtracts the second temporary spot image I2 (˜ISP2) from the spot image I(ξ, η) (=ISP1+ISP2) and calculates a second differential spot image ΔI2(ξ, η) (=I−I2˜ISP1). In the second differential spot image ΔI2, the light intensity distribution ISP2 derived from the spot group SP2 is reduced, and the light intensity distribution ISP1 derived from the spot group SP1 remains as it is. Similar to the step S405f, processing for further reducing the light intensity signal derived from the spot group SP2 may be inserted.
In the step S405j, the second differential spot image ΔI2 is analyzed, and the position (ξu, j, k, ηu, j, k) of the spot constituting the spot group SP1 is detected. As described above, since the light intensity signal derived from the spot group SP2 is reduced in the second differential spot image ΔI2, the spot group SP1 is analyzed by analyzing it according to the same method as in the step S103a of the first embodiment, and the position of each spot to be configured can be detected.
The above processing acquires the position (ξu, j, k, ηu, j, k) of each spot constituting the spot group SP1 formed by the first light of the test lens 12, and the position (ξb, j, k, ηb, j, k) of each spot constituting the spot group SP2 formed by the second light. In the step S405k, the shift amounts of the spot groups SP1 and SP2 are calculated using the acquired spot positions. More specifically, for example, the position of each spot (ξu, j, k, ηu, j, k, ξb, j, k, ηb, j, k) and the calibration value of MLA 2 (ξ0, j, k, η0, j, k) and the designed values ρu and ρb of the curvatures of the wavefronts W1 and W2 are substituted for the expressions (45) to (48), and the shift amounts (ξu, 0, 0, ηu, 0, 0) and (ξb, 0, 0, ηb, 0, 0) are obtained.
Assume that Nspot is the number of spots used for the calculation. The second term in the symbol Σ in the expressions (45) to (48) shows the spot position where the test lens 12 and the imaging lens 14 are formed as designed and axisymmetric about the measurement optical axes 1c and 1d.
The processing using the expressions (45) to (48) does not necessarily have to be performed for all spots constituting the spot group, but may be performed by extracting a spot group near the measurement optical axis 1d, a specific spot group, or the like.
In the step S405l, the provisional back surface inclination is calculated. The provisional back surface inclination is set to the back surface inclinations Δθl, x, Δθl, y calculated by substituting the shift amount (ξu, 0, 0, ηu, 0, 0) of the spot group SP1 obtained in the step S405k and the shift amount (ξb, 0, 0, ηb, 0, 0) of the spot group SP2 for the expressions (28) and (29) in the second embodiment
The above calculation method uses for the spare lens 16 a lens formed based on the same designed shape as the test lens 12, but as long as the third surface and the first surface 12a may have the same designed shape, the shape of the fourth surface, the thickness of the lens, etc. are not questioned. The designed shape is not necessarily the same as that of the third surface. For example, the shape may be similar only in the area to be analyzed, or only the curvature component may be matched.
In calculating the shift amount in the step S405k, the expressions (16) and (17) of the first embodiment may be used. At this time, since the shift amount (ξu, 0, 0, ηu, 0, 0) of the spot group SP1 of the first light is unnecessary, the calculation of the second temporary spot image I2 is unnecessary. Accordingly, the steps S405h to S405j are omitted.
In addition, before the second temporary spot image I2 is calculated in the step S405h, the intensity Ib, j, k and the spot radius wb,j,k of each spot may be acquired from the first differential spot image ΔI1. For example, the intensity Ib, j, k and the radius wb,j,k of each spot are acquired by cutting the spots one by one from the spot image I(ξ, η) and by fitting them with the expression (49).
This embodiment does not need to fit a spot image including a plurality of spots with a predetermined function, and improves the robustness of the analysis in comparison with the first embodiment.
The first differential spot image ΔI1 includes a noise ΔISP1 (=ISP1−ISP3) caused by the difference between the spot group SP1 and the spot group SP3. Thus, in the step S405g, in detecting the position of each spot constituting the spot group SP2, a detection error occurs due to this noise ΔISP1. However, in calculating the shift amount of the spot group, the influence of this detection error is insignificant. This point will be described below.
The noise ΔISP1 is generated in accordance with the difference between the spot groups SP1 and SP3 and periodically appears similarly to the spot groups SP1 and SP3, and its period is about the same as the period (=spot interval) qu of the spot group SP1. On the other hand, as shown by the expressions (5) and (6), the first surface 12a and the second surface 12b of the test lens 12 have different curvatures from each other, so the period of the spot group SP2 is qb, which is different from qu of the period of the spot group SP1. Thus, the noise appearing in the first differential spot image ΔI1 shows a period different from that of the spot group SP2. Hence, the positional relationship between each spot constituting the spot group SP2 and the adjacent noise differs for each spot, and the position detecting error caused by the noise in the step S405g shows a different magnitude and direction for each spot.
As shown in the expressions (45) to (48), in calculating the shift amount for the entire spot group, the shift amounts of the Nspot spots are averaged. Hence, the shift amount calculating error for the entire spot group is reduced down to about 1/√{square root over ( )}Nspot of the position detecting error of each spot due to the averaging effect. For example, if Nspot=100, it is reduced down to about 1/10. If the position detecting error of each spot is distributed in the range of −100 to 100 nm, the error of the shift amount of the spot group is about 10 nm. If the distance between the image sensor 3 and the MLA 2 is 5 mm, the calculation error of the surface inclination is assumed to be about 2 μrad, which is sufficiently small compared to the eccentric accuracy required for a general lens.
As described above, the first differential spot image ΔI1 includes the noise ΔISP1 with the period qu. In calculating the second differential spot image ΔI2 from the first differential spot image ΔI1, this embodiment converts the period from qu to qb by introducing the step of calculating the second temporary spot image I2 through the steps S405g and S405h. Therefore, the noise position is separated from the spot in the second differential spot image ΔI2 where the spots are arranged at the interval qu. As a result, in detecting the position of the spot included in the second differential spot image ΔI2 in the step S405j, the detection error is suppressed, and in the subsequent step S405k, the shift amount of the spot group SP1 is accurately calculated.
Assume that the first differential spot image ΔI1 is directly set to the second temporary spot image I2 without performing the step S405h. In this case, since the second temporary spot image I2 includes the noise ΔISP1 (=ISP1−ISP3) in addition to the light intensity distribution ISP2 caused by the spot group SP2, I2=ISP2+ISP1−ISP3 is established. The step S405i calculates the second differential spot image ΔI2 as the difference (=I−I2) between the spot image I (=ISP1+ISP2) and the second temporary spot image I2 (=ISP2+ISP1−ISP3). As a result, the second differential spot image ΔI2 coincides with the first temporary spot image ISP3, and the steps S405j and S405k calculate the shift amount of the spot group SP3 instead of the shift amount of the spot group SP1. In other words, the difference between the shift amounts of the spot groups SP1 and SP3 is added as an error to the shift amount of the spot group SP1 calculated in the step S405k.
Hence, this embodiment performs the steps S405g and S405h to calculate the second temporary spot image I2.
This embodiment is the same as the tenth embodiment in using the eccentricity measuring apparatus 100 illustrated in
The step S400 in the tenth embodiment acquires the first temporary spot image ISP3 using the spare lens 16, whereas the step S400 according to this embodiment provides the calculation based on the designed value of the test lens 12 and the calibration value of the MLA 2. More specifically, for example, the first temporary spot image ISP3 is simulated using the expression (50).
The expression (50) is derived from the expressions (5) to (11) described above, and refers to the designed value of the test lens 12 and the calibration value of MLA 2. As illustrated in the first embodiment, the spot image I may be fitted by the expression (11), and the light intensity distribution calculated by the second term on the right side may be set to the first temporary spot image ISP3. The initial values of the parameter groups X, W, and R required for this fitting are calculated based on the designed value of the test lens 12 and the calibration value of the MLA 2 in the same procedure as the step S101 in the first embodiment.
The steps S400 and S405 according to this embodiment simulate the first temporary spot image ISP3 and form the second temporary spot image I2 based on the result, but may simulate the second temporary spot image I2 and form the first temporary spot image ISP3 based on the result. Further, both the first temporary spot image ISP3 and the second temporary spot image I2 may be obtained by simulation, and the first differential spot image ΔI1 and the second differential spot image ΔI2 may be obtained by subtracting each from the spot image I.
Since no spare lens is necessary, this embodiment can reduce the measurement cost as compared with the tenth embodiment.
This embodiment is the same as the tenth embodiment in using the eccentricity measuring apparatus 100 illustrated in
The step S400 uses an illustrated spare lens 17. The spare lens 17 has a fifth surface and a sixth surface, and is formed based on the designed shapes of the first surface 12a and the second surface 12b of the test lens 12. An antireflection film is formed on the fifth surface, and most of the light incident on the fifth surface transmits through it. In the step S400, similar to the tenth embodiment, the second temporary spot image I2 is captured according to the flowchart in
In the step S405 in this embodiment, a provisional back surface inclination is calculated according to a flowchart in
Thus, the steps S505f, S505g, S505k, and S505l are performed to calculate the provisional back surface inclination. Since the processing in the steps S505f, S505g, S505k, and S505l are the same as that in the steps S405f, S405g, S405k, and S405l in
Similar to the case of the tenth embodiment, since it is unnecessary to fit a spot image including a plurality of spots, this embodiment improves the robustness of the analysis as compared with the first embodiment.
The step S400 may include both the step of capturing the first temporary spot image ISP3 using the spare lens 16 and the step of capturing the second temporary spot image I2 using the spare lens 17. In this case, in the step S405, the step S505m may be omitted from the flowchart in
Embodiment(s) of the present invention can also be realized by a computer of a system or apparatus that reads out and executes computer executable instructions (e.g., one or more programs) recorded on a storage medium (which may also be referred to more fully as a ‘non-transitory computer-readable storage medium’) to perform the functions of one or more of the above-described embodiment(s) and/or that includes one or more circuits (e.g., application specific integrated circuit (ASIC)) for performing the functions of one or more of the above-described embodiment(s), and by a method performed by the computer of the system or apparatus by, for example, reading out and executing the computer executable instructions from the storage medium to perform the functions of one or more of the above-described embodiment(s) and/or controlling the one or more circuits to perform the functions of one or more of the above-described embodiment(s). The computer may comprise one or more processors (e.g., central processor (CPU), microprocessor (MPU)) and may include a network of separate computers or separate processors to read out and execute the computer executable instructions. The computer executable instructions may be provided to the computer, for example, from a network or the storage medium. The storage medium may include, for example, one or more of a hard disk, a random-access memory (RAM), a read only memory (ROM), a storage of distributed computing systems, an optical disk (such as a compact disc (CD), digital versatile disc (DVD), or Blu-ray Disc (BD)™), a flash memory device, a memory card, and the like.
While the present invention has been described with reference to exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.
This application claims the benefit of Japanese Patent Application No. 2018-207797, filed on Nov. 2, 2018, which is hereby incorporated by reference herein in its entirety.
Number | Date | Country | Kind |
---|---|---|---|
2018-207797 | Nov 2018 | JP | national |