Design method for optical curved surface

Abstract

A design method for an optical curved surface, includes: dividing an optical curved surface into a plurality of rectangular areas that are divided by a lattice; defining an original lattice point on the lattice, a backbone line that crosses the original lattice point on the lattice, and standard lattice points that are the lattice points located on the backbone line except the original lattice point; applying curvatures to all of the lattice points; applying inclinations to the original lattice point and the standard lattice points; applying sag to the original lattice; calculating a sectional shape of the curved surface along the backbone line based on the sag and inclination of the original lattice point and the curvatures of the lattice points on the backbone line; calculating sags of the standard lattice points based on the calculated sectional shape; calculating sectional shapes along orthogonal lines that are orthogonal to the backbone line based on the calculated sags and the applied inclinations of the standard lattice points and the applied curvatures of the lattice points on the orthogonal lines; and representing the rectangular areas as mathematical functions respectively based on the calculated sectional shapes.

Description

BACKGROUND OF THE INVENTION

The present invention relates to a design method for a non-rotationally symmetrical optical curved surface such as a progressive power lens or an fθ lens.

In general, an aspherical surface, particularly a complex aspherical surface, is designed after repeated trial and error. An optical performance of a default shape is evaluated, then at least one parameter of the shape is changed and the performance is evaluated. The change of the parameter and the evaluation are repeated until a satisfactory result is obtained.

An aspherical surface can be represented by a mathematical function. In a first conventional method, the entire area of the aspherical surface is represented by a high-order polynomial. A second conventional method divides the aspherical surface into a plurality of areas and a low-order polynomial is defined for each of the areas. In the first method, since the change of the parameter to improve performance of the specific area influences the performance of the entire area, the first method increases number of times of the trial and error to optimize the entire area with considering the influence.

In the second method, since the change of the parameter of the specific area has no influence with the other area, it eases to improve the local optical performance and it is suited for designing the complex aspherical surface. A lens design method using the second method is disclosed in, for example, Japanese Provisional Patent Publication No. Sho 55-146412, or International Patent Re-publication WO 96/11421.

Sho 55-146412 discloses a design method for an aspherical surface of a progressive power spectacle lens. The disclosed method divides the entire area of the aspherical surface into rectangular areas that are divisions divided by a lattice. Each of the rectangular areas is represented by a bicubic polynomial function and coefficients (parameters) of the function are found to solve simultaneous equations so that the values of the expression itself, the differential and the quadratic differential are continuous, respectively. Namely, since the disclosed method aims to provide distribution of the prismatic powers (inclinations of the lens surface), default distribution of the prismatic powers is applied at the beginning, the shape of the lens surface is calculated based on the default distribution, and then the surface shape is represented by the bicubic polynomial functions for the rectangular areas to evaluate the optical performance.

A design method disclosed in WO 96/11421 sets default values of radius of curvature on principal points, and it divides the entire area of the lens surface into rectangular areas that are divided by lattice. Ray tracing is applied to each of the rectangular areas to evaluate the optical performance of the lens, and then the radii of curvature are corrected according to the evaluation. Bi-cubic expressions are defined for the rectangular areas, and then the optical performance of the surface shape is checked. When there is room for improvement in the optical performance, the radius of curvatures are reset to re-define the bicubic polynomial functions. The cycle of check and reset is repeated until the optimum shape is obtained.

However, Sho 55-146412 only discloses the method for evaluating the optical performance of the lens surface, while it does not disclose how to change the parameters. Further, since the method finds the bicubic polynomial functions based on the differentials and the like after setting the refractive prism powers, the prismatic powers must be the parameters of the design even if the surface shape is changed according to the evaluation. However, since the lens shape represented by the distribution of the prismatic powers is not directly related to refractive power and astigmatism that are important items to evaluate the progressive power lens, it is difficult to determine how to change the prismatic powers based on the evaluation to correct aberrations, which increases number of times of the trial and error for the optimization, requiring huge calculations.

On the other hand, since the optical performance can be evaluated by the computerized ray tracing only after the lens shape is represented by the mathematical function, the method disclosed in WO 96/11421 does not actually evaluate the optical performance. Further, the radius of curvature in an X direction must be different from that in a Y direction in the non-rotationally symmetrical surface such as the progressive power lens, while the disclosed method does not distinguish the difference. Moreover, the publication does not show how to determine the lens shape based on the radii of curvature given for the principal points on the lens surface.

SUMMARY OF THE INVENTION

It is therefore an object of the present invention to provide a design method of an optical curved surface on the condition that areas of a lens surface divided by a lattice are represented by mathematical functions, which is capable of efficiently calculating a shape of the entire lens surface and a shape of the each area of the lens surface based on curvatures (refractive powers) of the lens surface as parameters that are directly related to aberrations.

For the above object, according to the present invention, there is provided a design method of an optical curved surface, which includes:

dividing an optical curved surface into a plurality of rectangular areas that are divided by the lattice;

defining an original lattice point on the lattice, a backbone line that crosses the original lattice point on the lattice, and standard lattice points that are the lattice points located on the backbone line except the original lattice point;

applying curvatures to all of the lattice points;

applying inclinations to the original lattice point and the standard lattice points;

applying sag to the original lattice point;

calculating a sectional shape of the curved surface along the backbone line based on the sag and inclination of the original lattice point and the curvatures of the lattice points on the backbone line;

calculating sags of the standard lattice points based on the calculated sectional shape;

calculating sectional shapes along orthogonal lines that are orthogonal to the backbone line based on the calculated sags and the applied inclinations of the standard lattice points and the applied curvatures of the lattice points on the orthogonal lines; and

representing the rectangular areas as mathematical functions respectively based on the calculated sectional shapes.

With this method, the sectional shape along the backbone line is calculated with integration based on parameters such as curvatures applied to the lattice points at first, and the sectional shapes along the orthogonal lines are calculated with integration. The surface shape in each of the rectangular areas are, for example, represented by a bicubic polynomial function, which enables evaluation of the optical performance. The curvatures applied as parameters can be replaced with sectional surface powers.

When the parameters are changed based on the evaluation result, the curvature is preferable as the parameter because of its direct relationship with the aberration.

The design method may further include:

evaluating optical performance of the optical curved surface represented by the mathematical functions; and

changing at least one of the applied sag, inclinations and curvatures based on the evaluated result.

The calculation of the sectional shape along the backbone line may include:

a first integration of the curvatures with the value corresponding to the inclination of the original lattice point as an integration constant to obtain a distribution of the values corresponding to the inclinations along the backbone line; and

a second integration of the distribution of the values corresponding to the inclinations with the sags of the original lattice point as an integration constant to obtain the sectional shape along the backbone line.

The calculation of the sectional shape along the orthogonal line may include:

a first integration of the curvatures with the value corresponding to the inclination of the standard lattice point on the orthogonal lattice point as an integration constant to obtain a distribution of the values corresponding to the inclinations along the orthogonal line; and

a second integration of the distribution of the values corresponding to the inclinations with the sags of the standard lattice point on the orthogonal line as an integration constant to obtain the sectional shape along the orthogonal line.

When the optical curved surface has a distribution of the variation of curvature, it is preferable that the rectangular areas are relatively small where the variation of curvature is relatively large and the rectangular areas are relatively large where the variation of curvature is relatively small.

Further, the optical curved surface may be divided by a first lattice whose pitch is relatively large and by a second lattice whose pitch is relatively small. In such a case, the parameters are applied to the lattice points of the first lattice, and the design method further includes:

calculating the parameters for the lattice points of the second lattice by interpolating the parameters applied to the lattice points of the first lattice,

calculating a sectional shape of the curved surface along the backbone line based on the parameters of the lattice points of the second lattice, and

calculating a sectional shape of the curved surface along orthogonal lines that are orthogonal to the backbone line based on the parameters of the lattice points of the second lattice.

DESCRIPTION OF THE ACCOMPANYING DRAWINGS

FIG. 1 is a flowchart that shows a design method according to a first embodiment;

FIG. 2 shows parameters that are applied to the design method of FIG. 1 for designing an fθ lens;

FIG. 3 is a table showing concrete numerical values of the parameters in FIG. 2 ;

FIGS. 4A-4E show calculation steps of the design method in FIG. 1 ;

FIG. 5 is a flowchart that shows a design method according to a second embodiment;

FIG. 6 shows parameters of a lattice points of a first lattice that are applied to the design method of FIG. 5 for designing a progressive power lens;

FIG. 7 is a table showing concrete numerical values of the parameters of the lattice points of the first lattice in FIG. 6 ;

FIG. 8 shows parameters of a second lattice that are applied to the design method of FIG. 5 for designing the progressive power lens;

FIG. 9 is a table showing concrete numerical values of the parameters of the lattice points of the second lattice in FIG. 8 where 3≦m≦13;

FIG. 10 is a table showing concrete numerical values of the parameters of the second lattice in FIG. 8 where 14≦m≦30;

FIG. 11 is a table showing concrete numerical values of the parameters of the second lattice in FIG. 8 where 31≦m≦47;

FIGS. 12A-12E show calculation steps of the design method in FIG. 5 ;

FIG. 13A shows distribution of surface astigmatism of the progressive power lens designed by the method of the second embodiment;

FIG. 13B shows distribution of average refractive power of the progressive power lens designed by the method of the second embodiment;

FIG. 14 is a graph showing a target value of refractive power in optimization;

FIGS. 15A-15E are graphs showing progress of the optimization by the design method of FIG. 5 ;

FIGS. 16A-16E are graphs showing progress of the optimization by a conventional SSD (Shape Spline Differential) method;

FIGS. 17A-17E are graphs showing progress of the optimization by a conventional PLY (Polynomial) method; and

FIG. 18 is a graph showing the optimizing speeds of the design method of FIG. 5 , the SSD method and the PLY method.

DESCRIPTION OF THE EMBODIMENTS

First Embodiment

A first embodiment shows the design method to design a non-rotationally asymmetrical fθ lens. FIG. 1 is a flowchart showing a design method according to a first embodiment, FIG. 2 shows parameters that are applied to the design method of FIG. 1 , FIG. 3 is a table showing concrete numerical values of the parameters and FIGS. 4A-4E show calculation steps of the design method.

Settings of coordinate axes and rectangular areas are described with reference to FIG. 2. A Z-axis is defined as an axis that is perpendicular to the sheet of FIG. 2 , an XY-plane is perpendicular to the Z-axis. In the XY plane, a curved surface of the fθ lens is divided into a plurality of rectangular areas that are divided by lattice. The lattice is shown by I×J matrix, particularly by 4×8 matrix in FIG. 2 . Intersections of boundary lines shown by dotted lines are defined as lattice points whose XY coordinates are represented by (Xi, Yj) [0≦i≦I, 0≦j≦J, i, j are integers]. An original lattice point is defined as (XK, YL). The original lattice point is located in the center of an effective area of the curved surface in this embodiment, that is, K=2 and L=4. Since K and L are constants and i and j are discrete variables, XK is a constant to indicate the X-coordinate of the original lattice point, and Xi is a discrete variable to indicate the X-coordinate of one of the lattice points where 0≦i≦I, X is a continuous variable to indicate the X-coordinate of an arbitrary point.

In a step S 1 of FIG. 1 , various parameters, which are a sag in the Z-axis direction ZKL, an inclination in the X-axis direction BXKL, an inclination in the Y-axis direction BYKL, a curvature in the X-axis direction CXKL and a curvature in the Y-axis direction CYKL, are applied to the original lattice point (XK, YL).

In a step S 2 , various parameters, which are the inclinations in the X-axis direction BXKj, the curvatures in the X-axis direction CXKj and the curvatures in the Y-axis direction CYKj, are applied to each of lattice points (XK, Yj) [j≠L]. In the description, the straight line “X=XK” that is parallel to the Y-axis and crosses the original lattice point (XK, YL) is referred as the “backbone line” and the lattice points that are arranged on the backbone line except the original lattice point (XK, Yj) [j≠L] are referred as standard lattice points. The term “lattice points on the backbone line” includes the original lattice point and the standard lattice points.

In a step S 3 , the curvature in the X-axis direction CXij is applied to each of the lattice points (Xi, Yj) [i≠K] that are all of lattice points except the lattice points on the backbone line.

The above described three steps set the parameters. In FIG. 2 , a symbol “·” represents Zij, “−” represents BXij, “|” represents BYij, “Å” represents CXij and “<” represents CYij. Namely, the curvatures in the X-axis direction CXij is applied to the every lattice point (Xi, Yj), the inclination in the X-axis direction BXKj and the curvature in the Y-axis direction CYKj are applied to the lattice points (XK, Yj) on the backbone line, and the sag ZKL and the inclination BYXL is only applied to the original lattice point (XK,YL).

FIG. 3 is a table showing concrete numerical values of the parameters in FIG. 2 . Symbols ΔCYKj and ΔCXij show differences of the curvatures CYKj and CXij from the values of CYKL and CXKL at the original lattice point. The curved surface of the fθ lens has width of ±60.00 mm in the Y-axis direction and width of ±6.00 mm in the X-axis direction. The curved surface is divided into 8 areas in the Y-axis direction and into 4 areas in the X-axis direction. The values ΔCYKj and BXKj are applied to the lattice points (XK, Yj) on the backbone line, as shown in the left side of the table. The values ΔCXij are applied to the other lattice points that are not arranged on the backbone line as shown in the right side of the table. The parameters applied to the original lattice point (XK, YL) are shown as follows.

CXKL = −8.28 × 10                                                                        −3 BYKL = 0.00
CYKL = −3.73 × 10                                                                        −3 ZKL = 0.00

In the following steps, mathematical functions that represent the surface shape are calculated based upon the applied parameters described above Sag Z, inclination BX and BY, curvatures CX and CY at the arbitrary lattice point (X, Y) are defined as the following functions.

z=f(X, Y)

BX=g(X, Y)

BY=h(X, Y)

CX=u(X, Y)

CY=v(X, Y)

The values of the parameters at each of the lattice points are defined as follows.

Zij=f(Xi, Yj)

BXij=g(Xi, Yj)

BYij=h(Xi, Yj)

CXij=u(Xi, Yj)

CYij=v(Xi, Yj)

The functions u and v, which represent distributions of the curvatures, are integrated to obtain the functions that represent distribution of values corresponding to the inclinations. Since the integration of the functions u and v do not directly get the functions g and h that represent distribution of the inclination, intermediate functions p and q are defined as follows. The function f that represents distribution of the sags is found by integrating the functions g and h. $$ p ⁡ ( X , Y ) ≡ g 1 + g 2 $$ q ⁡ ( X , Y ) ≡ h 1 + h 2

The above described definitions result the following relationships. $$ g ⁡ ( X , Y ) ≡ ∂ f ∂ X $$ h ⁡ ( X , Y ) ≡ ∂ f ∂ Y $$ u ⁡ ( X , Y ) ≡ ∂ p ∂ X $$ v ⁡ ( X , Y ) ≡ ∂ q ∂ Y

Further, the differential w of the function f(X, Y) by X and Y is defined as follows. $$ w ⁡ ( X , Y ) ≡ ∂ 2 ⁢ f ∂ X ⁢ ∂ Y

In a step S 4 , the function v(XK, Y) that defines the continuous variation of the curvatures in the Y-axis direction CYKj along the backbone line is found, and the function f(XK, Y) that defines the continuous variation of the sags ZKj along the backbone line is calculated from the function v(XK, Y) by the process including two-step integration. The function v(XK, Y) defines a curve to smoothly connect the discrete curvatures of the lattice points on the backbone line on a graph showing relationship between the curvature and the Y-coordinate. In the same manner, the function f(XK, Y) defines a curve to smoothly connect the discrete sags of the lattice points on the backbone line on a graph showing relationship between the sag and the Y-coordinate. The step S 4 corresponds to the calculation shown in FIG. 4 A. It is noted that the curved surface is divided into the areas of 2×4 matrix in order to simplify the drawings in FIGS. 4A-4E , while the process in each step is common in spite of the size of the matrix.

The process of the step S 4 includes two integral steps. In the first integral step, the function v(XK, Y) is integrated with respect to Y with the integration constant q(XK, YL) that corresponds to the inclination in the Y-axis direction BYKL at the original lattice point (XK, YL). This process results in the function q(XK, Y) that defines the continuous variation of the inclination in the Y-axis direction along the backbone line as follows. The function q(XK, Y) continuously connects the discrete distribution of the inclination in the Y-axis direction BYKj (=q(XK, Yj)) of the lattice points on the backbone line. $$ q ⁡ ( X K , Y ) = ∫ Y L Y ⁢ v ⁡ ( X K , Y ) ⁢   ⁢ ⅆ Y + q ⁡ ( X K , Y L )

Where the integration constant q(XK,YL) is expressed by the following equation. $$ q ⁡ ( XK , YL ) = h ⁡ ( XK , YL ) 1 + h ⁡ ( XK , YL ) 2 = BYKL 1 + BYKL 2

Further, the function h (XK, Y) that defines the continuous variation of the inclinations in the Y-axis direction is found based on the function q(XK, Y) according to the above described relationship between the function q(X, Y) and the function h as shown in follows. $$ h ⁡ ( XK , Y ) = q ⁡ ( XK , Y ) 1 - q ⁡ ( XK , Y ) 2

In the second integral step at the step S 4 , the function h(XK, Y) is integrated with respect to Y with the sag ZKL of the original lattice point to calculate the function f(XK, Y) as shown by the following equation. $$ f ⁡ ( XK , Y ) = ∫ Y L Y ⁢ h ⁡ ( X K , Y ) ⁢   ⁢ ⅆ Y + f ⁡ ( X K , Y L ) = ∫ Y L Y ⁢ h ⁡ ( X K , Y ) ⁢   ⁢ ⅆ Y + Z KL

The process of the step S 4 determines BYKj=h(XK,Yj) and ZKj=f(XK,Yj), which finds the sectional shape of the curved surface along the backbone line.

In a step S 5 , the functions u(X, Yj) are calculated first. The functions u(X, Yj) define the continuous variations of the curvatures in the X-axis direction CXij of the lattice points (Xi, Yj) arranged on orthogonal lines Y=Yj. The orthogonal line cross the backbone line at right angles. Further, the functions f(X, Yj) that define the continuous variations of the sags ZKj of the lattice points on the orthogonal lines Y=Yj are calculated from the functions u(X, Yj) by the process including two integration steps. The process of the step S 5 corresponds to the calculation of FIG. 4 B. This process finds the sectional shapes of the curved surface along the orthogonal lines Y=Yj.

The process of the step S 5 includes two integral steps for the every orthogonal line. In the first integral step, as shown in the following equation, the functions u(X, Yj) are integrated with respect to X with the integration constants p (XK, Yj) that correspond the inclinations in the X-axis direction on the lattice points (XK, Yj) on the backbone line to obtain the functions p(X, Yj) that define the continuous variations of the inclinations in the X-axis direction at the lattice points (XK, Yj) arranged on the orthogonal lines. The functions p(X, Yj) continuously connect the discrete distributions BXij (=p(Xi, Yj)) of the inclinations in the X-axis direction. $$ p ⁡ ( X , Y j ) = ∫ X K X ⁢ u ⁡ ( X , Y j ) ⁢   ⁢ ⅆ X + p ⁡ ( X K , Y j )

Where the integration constanst p(XK,Yj) are expressed by the following equation. $$ p ⁡ ( XK , Yj ) = g ⁡ ( XK , Yj ) 1 + g ⁡ ( XK , Yj ) 2 = BXKj 1 + BXKj 2

Further, the function g (x, Yj) that defines the continuous variation of the inclinations in the X-axis direction at the lattice points arranged on the orthogonal line Y=Yj is found based on the function p(X, Yj) according to the above described relationship between the function p(X, Y) and the function g as shown in follows. $$ g ⁡ ( X , Yj ) = p ⁡ ( X , Yj ) 1 - p ⁡ ( X , Yj ) 2

In the second integral step at the step S 5 , the function g(X, Yj) is integrated with respect to X with the sag ZKj of the lattice point on the backbone line as the integration constant to calculate the function f(X, Yj) for the every orthogonal line as shown by the following equation. $$ f ⁡ ( X , Yj ) = ∫ X K X ⁢ g ⁡ ( X , Yj ) ⁢   ⁢ ⅆ X + f ⁡ ( X K , Yj ) = ∫ X K X ⁢ g ⁡ ( X , Yj ) ⁢   ⁢ ⅆ X + Z Kj

The process in the step S 5 determines BXij=g(Xi,Yj) and Zij=f(Xi,Yj), which finds the sectional shape of the entire area of the curved surface.

Steps S 6 and S 7 are processes to calculate the inclinations in the Y-axis direction BYij. The inclinations in the X-axis direction BXij that are calculated as rows in the X-axis direction are regarded as rows in the Y-axis direction to calculate the function g(Xi, Y) that defines the continuous variation of the inclinations in the X-axis direction along the Y-axis direction. Then, the differential of the function g(Xi, Y) by Y is integrated with respect to X to calculate the inclinations in the Y-aixs direction BYij.

In the step S 6 , the function g(Xi, Y) that continuously connects the discrete distribution of the inclinations BXij is differentiated by Y to find differential Wij at each lattice point (Xi, Yj) as shown in the following equations. $$ w ⁡ ( Xi , Y ) = ∂ g ⁡ ( Xi , Y ) ∂ Y $$ Wij = w ⁡ ( Xi , Yj )

The process in the step S 6 corresponds to the calculation of FIG. 4 c . Symbol “&Asteriskpseud;” represents the differential Wij of each lattice point.

In the step S 7 , the functions w(X, Yj), which defines continuous variations of the differentials Wij of the lattice points (Xi, Yj) on the orthogonal lines Y=Yj, is integrated with respect to X with the integration constants h(XK, Yj) that are the inclination at the lattice points on the backbone line to calculate the functions h(X, Yj) as shown in the following equation. $$ h ⁡ ( X , Yj ) = ∫ X K X ⁢ w ⁡ ( X , Yj ) ⁢ ⅆ X + h ⁡ ( X K , Yj ) = ∫ X K X ⁢ w ⁡ ( X , Yj ) ⁢ ⅆ X + BY Kj

The functions h(X, Yj) define continuous variations of the inclination in the Y axis direction along the orthogonal lines Y=Yj. In accordance with the function h(X, Yj), the inclinations in the Y-axis direction BYij (=h(Xi, Yj)) at each lattice point (Xi, Yj) is calculated. The process of the step S 7 corresponds to the calculation of FIG. 4 D.

The processes from the steps S 1 -S 7 determine four kinds of parameters Zij, BXij, BYij and Wij that show the shape of the curved surface for all of the lattice points (Xi, Yj) [0≦i≦I, 0≦j≦J]. In a step S 8 , a rectangular area Rij={ (X, Y) |X i−1 ≦X≦Xi, Y j−1 ≦Y≦Yj} is represented by the following bicubic polynomial function with respect to X and Y. $$ fij ⁢ ( X , Y ) = ∑ a = 0 3 ⁢ ∑ b = 0 3 ⁢ ( γ ab ( ij ) ⁢ ( X - X i - 1 ) ) a · ( Y - Y j - 1 ) b

Sixteen coefficients included in the bicubic polynomial function are defined by simultaneous equations using sixteen parameters that are the four parameters for each of four lattice points surrounding the rectangular area. The coefficients γab (ij) for the rectangular area Rij whose order is ij where the lattice points (X i−1 , Y j−1 ) and (Xi, Yj) are opposite corners are calculated from the following matrix. $$ [   ⁢ γ 00 γ 01 γ 02 γ 03 γ 10 γ 11 γ 12 γ 13 γ 20 γ 21 γ 22 γ 23 γ 30 γ 31 γ 32 γ 33 ] =   ⁢ ⁢ [   ⁢ 1 0 0 0 0 1 0 0 - 3 / Δ ⁢   ⁢ X i - 1 2 - 2 / Δ ⁢   ⁢ X i - 1 3 / Δ ⁢   ⁢ X i - 1 2 - 1 / Δ ⁢   ⁢ X i - 1 2 / Δ ⁢   ⁢ X i - 1 3 1 / Δ ⁢   ⁢ X i - 1 2 - 2 / Δ ⁢   ⁢ X i - 1 3 1 / Δ ⁢   ⁢ X i - 1 2 ⁢   ] ⁢   · ⁢   ⁢ &AutoLeftMatch; [   ⁢ Z i - 1 ⁢   ⁢ j - 1 BY i - 1 ⁢ j - 1 Z i - 1 ⁢ j BY i - 1 ⁢ j BX i - 1 ⁢ j - 1 W i - 1 ⁢ j - 1 BX i - 1 ⁢ j W i - 1 ⁢ j Z ij - 1 BY ij - 1 Z ij BY ij BX ij - 1 W ij - 1 BX ij W ij ] · [   ⁢ 1 0 - 3 / Δ ⁢   ⁢ Y j - 1 2 2 / Δ ⁢   ⁢ Y j - 1 3 0 1 - 2 / Δ ⁢   ⁢ Y j - 1 1 / Δ ⁢   ⁢ Y j - 1 2 0 0 3 / Δ ⁢   ⁢ Y j - 1 2 - 2 / Δ ⁢   ⁢ Y j - 1 3 0 0 - 1 / Δ ⁢   ⁢ Y j - 1 1 / Δ ⁢   ⁢ Y j - 1 2 ⁢   ]

Where ΔX i−1 =Xi−X i−1 and ΔY j−1 =Yj−Y j−1 .

The determination of the coefficients γab (ij) enables the calculation of the sags Z (=fij (X, Y)) at a arbitrary point in the rectangular area Rij and the differentials thereof by X and Y. Therefore, the coefficients γab (ij) are determined for every rectangular area, which determines the expression Z=f(X, Y) that represents the entire area of the curved surface as a set of the expressions fij (X, Y) [0≦i≦I, 0≦j≦J]. The expression Z=f(X, Y) guarantees that the values of the expression itself, the differential and the quadratic differential are continuous in the entire areas. The process of the step S 8 corresponds to the calculation of FIG. 4 E.

In a step S 9 , the optical performance is simulated based on the determined shape of the curved surface. And then, the optical performance is evaluated based on the result of the simulation. If the evaluation result shows satisfactory optical performance, the design process is terminated. On the other hand, if the evaluation result is no good, the curvatures CXij and CYKj are changed to optimize the performance in a step S 11 and the process returns to the step S 4 . The parameter to be changed in the step S 11 is not only limited to the curvatures, the inclinations in the X-axis direction BXKj of the lattice points on the backbone line can be adjusted.

In the above described embodiment, since the curved surface of the fθ lens is divided into a plurality of rectangular areas and each of the areas is represented by the respective mathematical function, the change of the parameter only affects the areas in the vicinity of the target area whose parameters are changed, which enables local optimization of the performance, providing the appropriate method for designing complex aspherical surfaces.

Since the curvatures CXij and CYKj that have direct relationship with aberrations are used as parameters, the effect of the change of the parameter can be directly grasped, which provides good design prospect. Further, in the fθ lens, the curvatures CYKj control linearity of scanning and field curvature in a principal scanning direction, the curvatures CXKj control field curvature in an auxiliary scanning direction and the inclinations BXKj control curvature of a scanning line (Bow), respectively. Therefore, the operator (or the computer software) can foresee the kind of the parameter and the changing amount thereof in order to reduce the aberration to improve the optical performance.

Second Embodiment

A second embodiment will be described hereinafter. In the second embodiment, the design method of the present invention is applied to design of a progressive power lens. An optical curved surface of the progressive lens is divided by a first lattice whose pitch is relatively large and by a second lattice whose pitch is relatively small in the second embodiment. Values are applied to the parameters of the lattice points of the first lattice. The parameters of the lattice point of the second lattice are calculated by interpolating the parameters applied to the lattice points of the first lattice. The surface shape is determined based on the parameters of the lattice point of the second lattice.

FIG. 5 is a flowchart that shows a design method according to the second embodiment, FIG. 6 shows parameters of the lattice points of the first lattice, FIG. 7 is a table showing concrete numerical values of the parameters of the first lattice in FIG. 6 , FIG. 8 shows parameters of the lattice points of the second lattice, FIGS. 9-11 are tables showing concrete numerical values of the parameters of the lattice points of the second lattice and FIGS. 12A-12E show calculation steps of the design method of FIG. 5 .

In the following description, a coordinate of a lattice point of the first lattice is defined as (Xi, Yj), and that of the second lattice is defined as (X˜m, Y˜n). Further, an original lattice point (XK, YL) of the first lattice is coincident with an original lattice point (X˜P, Y˜Q) of the second lattice, and thus the backbone lines of the lattices are coincident with each other.

The entire rectangular area of 100 mm×100 mm that includes the lens surface is divided into a plurality of rectangular areas by the first lattice that has 22×22 matrix, and it is divided by the second lattice that has 50×50 matrix. It will be noted K=11, L=12, I=J=22, P=Q=25, M=N=50. The size of the rectangular area divided by the first lattice is relatively small where the variation of curvature is relatively large and it is relatively large where the variation of curvature is relatively small. Since the progressive lens has an intermediate corridor where the variation of curvature is relatively large in the lower portion of the lens, the pitch of the first lattice is smaller in the intermediate corridor than in the other portion. The pitch of the second lattice is constant.

In step S 21 of FIG. 5 , various parameters, which are a sag in the Z-axis direction ZKL, an inclination in the X-axis direction BXKL, an inclination in the Y-axis direction BYKL, a curvature in the X-axis direction CXKL and a curvature in the Y-axis direction CYKL, are applied to the original lattice point (XK, YL) of the first lattice.

In step S 22 , various parameters, which are the inclinations in the X-axis direction BXKj, the curvatures in the X-axis direction CXKj and the curvatures in the Y-axis direction CYKJ, are applied to each of the standard lattice points (XK, Yj) [j≠L] of the first lattice arranged on the backbone line X=XK.

In step S 23 , the curvatures in the X-axis direction CXij is applied to each of the lattice points (Xi, Yj) [i≠K] of the first lattice except the lattice points on the backbone line. FIG. 6 shows the applied parameters for the lattice points of the first lattice. Meanings of the symbols are identical to that of FIG. 2 .

The above described three steps set the parameters. FIG. 7 is a table showing concrete numerical values of the parameters in FIG. 6 . In the table of FIG. 7 , sectional refractive powers Dxij and DYKj are applied to the lattice points in place of the curvatures CXij and CYKj. The sectional refractive power is directly related to the lens power, and thus it is preferable to be used as a parameter. The curvatures can be calculated by the following expressions.

CXij=Dxij/(In−1)

CYKj=DYKj/(In−1)

Where In is a refractive index of the lens.

The values DYKj and BXKj are applied to the lattice points (XK, Yj) on the backbone line as shown in the left side of the table. The values DXij are applied to the other lattice points that are not arranged on the backbone line as shown in the right side of the table. The parameters applied to the original lattice point (XK, YL) are shown as follows.

BYKL=0.00 ZKL=0.00

Since a spectacle lens has a circular outline in an uncut fashion, it is not necessary to give the parameters DXij for all of the lattice points (Xi, Yj) [0≦i≦22, 0≦j≦22]. In this embodiment, the parameters are applied to the lattice points covered by a circle whose diameter is 90 mm. The circle is set to cover the outline of an effective diameter 80 mm of the progressive lens. This reduces the number of parameters, reducing calculation. FIG. 7 shows the value ij min that is the minimum number of i on a j-th line being parallel to the X-axis and the value ij max that is the maximum number of i on the j-th line. In the tables, the parameters are not defined in the lattice points whose cells show “#####”.

In the second embodiment, the function f(X, Y), which defines the variation of the sectional shape of the curved surface, can be determined by the parameters applied to the lattice points of the first lattice as well as the first embodiment. However, since the maximum pitch of the first lattice is larger than 5 mm, the variation of the sectional shape will become like a broken line (a polygonal line), the change of the sectional shape cannot be smooth enough for a spectacle lens. In the second embodiment, therefore, the second lattice having a smaller pitch is used, and the function f(X, Y) is determined by the parameters of the lattice points of the second lattice.

In step 24 , the lens surface is divided by the second lattice having a smaller pitch than the first lattice. The parameters Z˜PQ, BY˜PQ, BX˜Pn, CY˜Pn and CX˜mn of the lattice points of the second lattice are calculated by interpolating the parameters ZKL, BYKL, BXKJ, CYKj and CXij applied to the lattice points of the first lattice. The following description is directed to the calculation in the step S 24 .

Since the second lattice is set so that the original lattice points are coincident with each other (X˜P=XK and Y˜Q=YL), the parameters on the original lattice point are coincident with each other (Z˜PQ=ZKL and BY˜PQ=BYKL). The curvatures in the Y-axis direction CY˜Pn and the inclinations BX˜Pn of the lattice points (X˜P, Y˜n) of the second lattice on the backbone line are calculated by interpolating the curvatures CYKj and the inclinations BXKj of the lattice points of the first lattice using a one-dimensional spline function along the backbone line X=XK.

Further, the curvatures CX˜mj on the intersections between the lines X=X˜m crossing the lattice points of the second lattice and the lines Y=Yj are calculated by interpolating the curvatures CXij of the lattice points of the first lattice using the one-dimensional spline function along the lines Y=Yj. Then curvatures CX˜mn of the lattice points (X˜m, Y˜n) of the second lattice are calculated by interpolating the curvatures CX˜mj using the one-dimensional spline function along the lines X=X˜m crossing the lattice points of the second lattice. These two steps for calculating the curvatures CX˜mn can be replaced by the single step using a two dimensional spline function. The one or two dimensional spline function should be selected so that at least first order differentials are continuous.

FIG. 8 shows the calculated parameters of the lattice point of the second lattice in the step S 24 Meanings of the symbols are identical to FIG. 2 . FIGS. 9-11 show the concrete values of the parameters for the lattice points of the second lattice except peripheral areas (0≦n≦2, 48≦n≦50, 0≦m≦2, 48≦m≦50) where the parameters are not established. FIG. 9 shows the parameters in the range of 3≦m≦13, FIG. 10 covers the range of 14≦m≦30 and FIG. 11 covers the range of 31≦m≦47. The tables shown in FIGS. 9-11 are originally formed as one table, where FIG. 9 is left side, FIG. 11 is right side of FIG. 10 , respectively.

The processes in steps S 25 through S 29 are similar to the steps S 4 through S 8 in FIG. 1 , while the rectangular areas divided by the second lattice are represented by the bicubic polynomial functions calculated from the parameters of the lattice points (X˜m, Y˜n) of the second lattice. FIGS. 12A-12E show the processes to find the bicubic polynomial functions that correspond to the steps S 25 -S 29 . In FIGS. 12A-12E , the second lattice is shown by 6×6 matrix for the convenience of description.

In step 30 , the optical performance is simulated based on the determined shape of the curved surface. And then, in step 31 , the optical performance is evaluated based on the result of the simulation. If the evaluation result shows satisfactory optical performance, the design process is terminated. On the other hand, if the evaluation result is no good, the inclinations BXKj, the curvatures CXij and CYKj are changed to optimize the performance in step 32 and the process returns to the step S 24 .

As described above, the use of the second lattice to represent the curved surface causes smooth variation of the curvatures as compared with the process that only uses the first lattice. Further, the parameters are applied to the lattice points of the coarse first lattice and the parameters of the lattice points of the fine second lattice are automatically calculated from the parameters of the first lattice, which reduces the number of the parameters to be set by an operator as compared with the case where the parameters are directly set to the lattice points of the fine second lattice, reducing the burden of the parameter setting.

In the second embodiment, since the progressive surface is divided into a plurality of rectangular areas and each of the areas is represented by the respective mathematical function, the change of the parameter only affects the areas in the vicinity of the target area whose parameters are changed, which enables local optimization of the performance.

The refractive powers DYKj and DXKj control astigmatism and focal power along a main meridian that extends along the backbone line, the refractive powers DXij control focal power of each area. The inclinations BXKj control inset that defines displacement of a near portion towards the nose, and the inclination BYKL controls a balance of the sags along the Y-axis direction. The operator (or the computer software) can foresee the kind of the parameter and the changing amount thereof in order to reduce the aberration to improve the optical performance.

FIG. 13A shows distribution of surface astigmatism of the progressive power lens designed by the method of the second embodiment, and FIG. 13B shows distribution of average surface refractive power thereof.

Next, the optimizing ability of the method for designing a progressive power lens according to the second embodiment of the present invention is now described. Default values are distant from optimum values in order to compare the ability of the design method shown in FIG. 5 with the conventional design method. It will be noted that the following descriptions are only directed for the optimization of the sectional shape along the backbone line to ease the description.

In the steps S 21 -S 23 in FIG. 5 , the default values are applied to the parameters ZKL, BYKL and BXKj of the lattice points of the first lattice. A lens shape of the opposite side of the progressive surface, a center thickness of the lens, refractive index of the lens or the like may be included in the parameters. In the steps S 24 -S 29 , the sectional shapes of the rectangular areas of the second lattice are represented by the bicubic polynomial functions.

In step 30 , the optical performance and the shape performance are simulated considering wearing conditions of glasses. The simulation of the optical performance includes analyzing the sectional refractive power by differential geometry, calculations of aberrations caused by the ray transmitted by the lens using ray tracing. The simulation of the shape performance includes a calculation of the edge thickness of the lens.

In the step S 31 , the performance of the lens is evaluated by comparing the simulated performance with the predetermined target performance. If the evaluation shows a poor result, the parameters are changed to improve the performance. The parameter is changed with damped least square method.

FIG. 14 is a graph showing a target value of transmitting refractive power along the backbone line in the optimization, and FIGS. 15A-15E are graphs showing progress of the optimization by the design method according to the second embodiment using the damped least square method for changing the parameters. In the following descriptions, the design method of the second embodiment is referred as Power Spline Integral (PSI) method.

The default parameters are, as shown in FIG. 15A , the surface refractive powers DYKj at the lattice points of the first lattice whose pitch is 3 mm. The applied default values represent a spherical surface, surface refractive power is 2 Diopter at the every lattice point. That is, the radius of curvature is 300 mm and the curvature is 3.33 m −1 when the refractive index equals 1.60.

On the basis of the discrete distribution of the surface refractive powers shown in FIG. 15A , the continuous variation of the sectional shape (FIG. 15 B), the continuous variation of the surface refractive power ( FIG. 15C ) and the continuous variation of the transmitting refractive power ( FIG. 15D ) are calculated. FIG. 15E shows the difference of the variation of the transmitting refractive power shown in FIG. 15D from the target variation thereof shown in FIG. 14 . In FIG. 15A , a symbol “∘” represents the default value, “Δ” represents the values after first optimization and “&Circlesolid;” represents the values after second optimization. Further, in FIGS. 15B-15E , a long dotted line represents the default value, a short dotted line represents the value after first optimization and a solid line represents the value after second optimization.

As shown in FIG. 15E , the PSI method reaches the target value after two optimization iterations. The difference between the object distance of the upper side of the lens and that of the lower side causes the difference of the performance at the default condition.

Two conventional methods, Shape Spline Differential (SSD) and Polynomial (PLY) method, are described hereinafter to verify the abilities of the PSI method. The SSD method calculates the continuous variation of the shape based on the discrete distribution of the shape using spline interpolation. The PLY method simulates the optical performance using a polynomial expression to represent the surface shape. In the following examples, the SSD method and the PLY method also use the damped least squares method for changing the parameters.

FIGS. 16A-16E are graphs showing progress of the optimization by the SSD method. The distribution of the sags ZKj at the lattice points whose pitch is 3 mm are applied as default values in the SSD method. The distribution of the default values in FIG. 16A is equivalent to that in FIG. 15A , which indicates the spherical curve whose radius of curvature is 300 mm, while the difference of the forms of the expression of the default values causes the difference of the graphs. That is, the horizontal axis of FIG. 16A shows the sag, while the horizontal axis of FIG. 15A shows the surface refractive power.

FIG. 16A shows the discrete distributions of the sags as the default values, values after first optimization and values after tenth optimization. FIGS. 16B , 16 C, 16 D and 16 E show the continuous variation of the sectional shape, the continuous variation of the surface refractive power, the continuous variation of the transmitting refractive power and the difference from the target values, respectively.

As shown in FIG. 16E , the SSD method does not reach the target value even after ten times of the optimization.

FIGS. 17A-17E are graphs showing progress of the optimization by the PLY method. In the PLY method of this example, the progressive lens surface is represented by the following polynomial expression. Sag Z is a 26th order function of a variable Y and coefficients Ai except zero and first orders are parameters to be changed. $$ Z ⁢ ( Y ) = ∑ i = 2 26 ⁢ Ai · Y i

The default shape is a spherical surface whose radius of curvature is 300 mm as well as the example of the PSI method. FIG. 17A shows aspherical components at the point Y=30. The aspherical components are represented by Zi=Ai·30 I and they correspond to the values of Ai. FIGS. 17B , 17 C, 17 D and 17 E show the continuous variation of the sectional shape, the continuous variation of the surface refractive power, the continuous variation of the transmitting refractive power and the difference from the target values, respectively.

As shown in FIG. 17E , a little difference from the target value remains after ten times of the optimization in the PLY method.

FIG. 18 is a graph showing a relationship between the number of the optimizations and the reduction of the merit function's value to compare the optimizing speeds of the PSI method, the SDD method and the PLY method. As compared with the PSI method, the SSD method has slower optimizing speed because of difficulty changing the sags by evaluating the transmitting refractive power, and the PLY method has also slower optimizing speed because of difficulty changing the local performance using the single polynomial expression. On the other hand, the PSI method reduces the merit function to about zero level after first optimization, which shows clearly superior effect as compared with the conventional methods.

The present disclosure relates to the subject matter contained in Japanese Patent Application No. HEI 10-216244, filed on Jul. 30, 1998, which is expressly incorporated herein by reference in its entirety.

Claims

1. A design method of an optical curved surface, comprising:dividing an optical curved surface into a plurality of rectangular areas that are divided by a lattice; defining an original lattice point on said lattice, a backbone line that crosses said original lattice point on said lattice, and standard lattice points that are the lattice points located on said backbone line except said original lattice point; applying curvatures to all of the lattice points; applying inclinations to said original lattice point and said standard lattice points; applying sag to said original lattice point; calculating a sectional shape of said curved surface along said backbone line based on said sag and inclination of said original lattice point and said curvatures of the lattice points on said backbone line; calculating sags of said standard lattice points based on the calculated sectional shape; calculating sectional shapes along orthogonal lines that are orthogonal to said backbone line based on the calculated sags and the applied inclinations of said standard lattice points and the applied curvatures of the lattice points on said orthogonal lines; and representing said rectangular areas as mathematical functions respectively based on the calculated sectional shapes.

2. The design method according to claim 1, further comprising:evaluating optical performance of said optical curved surface represented by said mathematical functions; and changing at least one of the applied sag, inclinations and curvatures based on the evaluated result.

3. The design method according to claim 1, wherein the calculation of the sectional shape along said backbone line includes:a first integration of said curvatures with the value corresponding to said inclination of said original lattice point as an integration constant to obtain a distribution of the values corresponding to said inclinations along said backbone line; and a second integration of said distribution of the values corresponding to said inclinations with said sag of said original lattice point as an integration constant to obtain said sectional shape along said backbone line.

4. The design method according to claim 1, wherein the calculation of the sectional shape along said orthogonal line includes:a first integration of said curvatures with the value corresponding to said inclination of the standard lattice point on said orthogonal line as an integration constant to obtain a distribution of the values corresponding to said inclinations along said orthogonal line; and a second integration of said distribution of the values corresponding to said inclinations with said sag of said standard lattice point on said orthogonal line as an integration constant to obtain said sectional shape along said orthogonal line.

5. The design method according to claim 1, wherein said rectangular areas are relatively small where the variation of curvature is relatively large and said rectangular areas are relatively large where the variation of curvature is relatively small.

6. The design method according to claim 1, wherein said optical curved surface is divided by a first lattice whose pitch is relatively large and by a second lattice whose pitch is relatively small, andwherein the parameters are applied to the lattice points of said first lattice, said design method further including: calculating the parameters for the lattice points of said second lattice by interpolating the parameters applied to said lattice points of the first lattice, calculating a sectional shape of said curved surface along said backbone line based on the parameters of the lattice points of said second lattice, and calculating a sectional shape of said curved surface along orthogonal lines that are orthogonal to said backbone line based on the parameters of the lattice points of said second lattice.

7. The design method according to claim 6, further comprising;evaluating optical performance of said optical curved surface; and changing at least one parameter based on the evaluated result.

8. The design method according to claim 6, wherein said rectangular areas divided by said first lattice are relatively small where the variation of curvature is relatively large and said rectangular areas divided by said first lattice are relatively large where the variation of curvature is relatively small.

9. The design method according to claim 6, wherein said rectangular areas divided by said second lattice have constant size in the entire area of said optical curved surface.

10. A design method of an optical curved surface, comprising:dividing an optical curved surface into a plurality of rectangular areas that are divided by a lattice; defining an original lattice point on said lattice, a backbone line that crosses said original lattice point on said lattice, and standard lattice points that are the lattice points located on said backbone line except said original lattice point; applying sectional surface powers to all of the lattice points; applying inclinations to said original lattice point and said standard lattice points; applying sag to said original lattice point; calculating a sectional shape of said curved surface along said backbone line based on said sag and inclination of said original lattice point and said sectional surface powers of the lattice points on said backbone line; calculating sags of said standard lattice points based on the calculated sectional shape; calculating sectional shapes along orthogonal lines that are orthogonal to said backbone line based on the calculated sags and the applied inclinations of said standard lattice points and the applied sectional surface power of the lattice points on said orthogonal lines; and representing said rectangular areas as mathematical functions respectively based on the calculated sectional shapes.

11. The design method according to claim 10, further comprising:evaluating optical performance of said optical curved surface represented by said mathematical functions; and changing at least one of the applied sag, inclinations and sectional surface powers based on the evaluated result.

12. A design method of an optical curved surface, comprising:dividing an optical curved surface into a plurality of rectangular areas by a first lattice whose pitch is relatively large; dividing said optical curved surface into a plurality of rectangular areas by a second lattice whose pitch is relatively small; applying parameters to the lattice points of said first lattice; calculating the parameters for the lattice points of said second lattice by interpolating the parameters applied to said lattice points of the first lattice; representing said rectangular areas divided by said second lattice as mathematical functions respectively based on the calculated parameters of the lattice points of said second lattice; evaluating optical performance of said optical curved surface represented by said mathematical functions; and changing at least one parameter applied to the lattice points of said first lattice based on the evaluated result.

13. A design method of an optical curved surface, comprising:dividing an optical curved surface in an XY plane, which is perpendicular to a Z axis showing a reference direction, into a plurality of rectangular areas that are divided by a lattice that is shown by an I×J matrix, XY coordinate of each lattice point being represented by (Xi, Yj) [0≦i≦I, 0≦j≦J, i, j are integers]; applying a sag in the Z-axis direction ZKL, an inclination in the X-axis direction BXKL, an inclination in the Y-axis direction BYKL, a curvature in the X-axis direction CXKL and a curvature in the Y-axis direction CYKL to an original lattice point (XK, YL); applying inclinations in the X-axis direction BXKj, curvatures in the X-axis direction CXKj and curvatures in the Y-axis direction CYKj to standard lattice points (XK, Yj) [j≠L] that are located on a predetermined backbone line X=XK that is parallel to the Y-axis and crosses said original lattice point (XK, YL); applying curvatures in the X-axis direction CXij to all of the lattice points (Xi, Yj) [i≠K] except the lattice points on said backbone line; integrating a function v(XK, Y) that defines the continuous variation of CYKj of the standard lattice points (XK, Yj) with respect to Y with the integration constant that is calculated from BYKL, said integration of the function v(XK, Y) resulting in a function h(XK, Y) that defines the continuous variation of the inclination in the Y-axis direction BYKj of the lattice points on said backbone line; integrating said function h(XK, Y) with respect to Y with the integration constant ZKL to find a function f(XK, Y) that defines the continuous variation of the sags ZKj of the lattice points (XK, Yj) on said backbone line; integrating functions u(X, Yj), which defines continuous variations of CXij of the lattice points (Xi, Yj) arranged on orthogonal lines Y=Yj that are orthogonal to said backbone line, with respect to X with the integration constants that are calculated from BXKj, said integration of said functions u(X, Yj) resulting functions g(X, Yj) that define continuous variations of the inclinations in the X-axis direction BXij of the lattice points on said orthogonal lines Y=Yj; integrating said functions g(X, Yj) with respect to X with integration constants ZKj that are the sags of the lattice points on said backbone line to find functions f(X, Yj) that define continuous variations of the sag Zij of the lattice points on said orthogonal lines Y=Yj; calculating BXij of the lattice points (Xi, Yj) based on said functions g(X, Yj); differentiating said functions g(Xi, Y), which define the continuous variations of BXij on lines X=Xi that is parallel to the Y-aixs, by Y to find differentials Wij on the lattice points (Xi, Yj); integrating functions w(X, Yj), which define the continuous variations of Wij of the lattice points (Xi, Yj) on said orthogonal lines Y=Yj, with respect to X to find the inclinations in the Y-axis direction BYij of the lattice points (Xi, Yi); representing the curved surfaces in said rectangular areas as bicubic polynomial functions based on four parameters Zij, BXij, BYij and Wij for each of four lattice points surrounding each rectangular area; evaluating optical performance of said optical curved surface represented by a set of said bicubic polynomial functions; and changing at least one of CXij, CYKJ and BXKj based on the evaluated result.

14. The design method according to claim 13, wherein said optical curved surface is divided by a first lattice that has an I×j matrix and by a second lattice that has a M×N matrix having a pitch smaller than said first lattice, XY coordinate of each lattice point of said second lattice being represented by (X˜m, Y˜n) [0≦m≦M, 0≦n≦N, m, m are integers], an original lattice point (X˜P, Y˜Q) of said second lattice corresponding to said original lattice point (XK, YL) of said first lattice; andwherein the sag ZKL, the inclinations in the X-axis direction BXKj, the inclination in the Y-axis direction BYKL, the curvatures in the X-axis direction CXij and the curvatures in the Y-axis direction CYKj are applied to the lattice points (Xi, Yj) of said first lattice, said design method further including: calculating the sag Z˜PQ, the inclinations in the X-axis direction BX˜Pn, the inclination in the Y-axis direction BY˜PQ, the curvatures in the X-axis direction CX˜mn and the curvatures in the Y-axis direction CX˜Pn for the lattice points (X˜m, Y˜n) of said second lattice by interpolating the corresponding parameters applied to said lattice points of the first lattice; calculating a sectional shape of said curved surface based on the parameters of the lattice points of said second lattice; representing said rectangular areas divided by said second lattice as mathematical functions respectively based on the calculated parameters of the lattice points of said second lattice; evaluating optical performance of said optical curved surface represented by said mathematical functions; and changing at least one parameter applied to the lattice point of said first lattice based on the evaluated result.