Method for high-accuracy wavefront measurement base on grating shearing interferometry

Abstract

A method for high-accuracy wavefront measurement based on grating shearing interferometry, which adopts a grating shearing interferometer system comprising an illuminating system, an optical imaging system under test, an object plane diffraction grating plate, an image plane diffraction grating plate, a two-dimensional photoelectric sensor, and a calculation processing unit. The object plane diffraction grating plate and the image plane diffraction grating plate are respectively arranged on the object plane and the image plane of the optical imaging system under test. The shearing phase of 1.sup.st-order diffracted beam and −1.sup.st-order diffracted beam is exactly extracted through phase shifting method, and the original wavefront is obtained by carrying out reconstruction algorithm according to a shear ratio of 2s, such that the accuracy of wavefront measurement of the optical imaging system under test is improved, wherein s is the shear ratio of the grating shearing interferometer.

Claims

1. A method for high-accuracy wavefront measurement based on grating shearing interferometry, comprising (i) placing an optical imaging system under test (3) in a grating shearing interferometer system, positioning an illuminating system (8) on an object space of the optical imaging system under test (3) and an image plane diffraction grating plate (4) on an image space of the optical imaging system under test (3), adjusting a first three-dimensional stage (2) to position an object plane diffraction grating plate (1) on an object plane of the optical imaging system under test (3), and adjusting a second three-dimensional stage (5) to position the image plane diffraction grating plate (4) on an image plane of the optical imaging system under test (3); (ii) moving the first three-dimensional stage (2) to enable grating line on the object plane diffraction grating plate (1) along a first grating (101) in a y axis direction to move into an object space field-of-view point position of the optical imaging system under test (3), moving the second three-dimensional stage (5) to enable a checkerboard grating or a one-dimensional grating in a corresponding direction on the image plane diffraction grating plate (4) to move into the image space field-of-view point position of the imaging system under test (3), wherein an included angle between a diagonal line direction of the checkerboard grating and an x axis or the y axis is 45 degrees, obtaining a series of shearing interferograms in the x axis direction by object plane or image plane grating phase-shifting, and solving and unwrapping to obtain a shearing phase φ.sub.x in the x axis direction: $\begin{array}{cc}{\phi }_x=\mathrm{unwrap}\left(\arctan \left[\frac{2a_1\left(\sin {\phi }_1+\sin {\phi }_{-1}\right)}{2a_1\left(\cos {\phi }_1+\cos {\phi }_{-1}\right)}\right]\right),& \end{array}$ wherein, at the moment, φ.sub.−1=ϕ(x, y)−ϕ(x−Δ, y) represents a shearing phase of 0.sup.th-order beam and −1.sup.st-order diffracted beam in the x axis direction, and φ.sub.1=ϕ(x−Δ, y)−ϕ(x, y) represents a shearing phase of 1.sup.st-order diffracted beam and 0.sup.th-order beam in the x axis direction, and A is an offset of 1.sup.st-order diffracted beam relative to 0.sup.th-order beam received on a two-dimensional photoelectric detector (6), and function unwrap (x) represents unwrapping x; (iii) moving the first three-dimensional stage (2) to enable a second grating (102) of the grating line on the object plane diffraction grating plate (1) along the x axis direction to move into the object space field-of-view point position of the optical imaging system under test (3); obtaining a series of shearing interferograms in the y axis direction by object plane or image plane grating phase-shifting, and solving and unwrapping to obtain a shearing phase ϕ.sub.y in the y axis direction: ${\phi }_y=\mathrm{unwrap}\left(\arctan \left[\frac{2a_1\left(\sin {\phi }_1+\sin {\phi }_{-1}\right)}{2a_1\left(\cos {\phi }_1+\cos {\phi }_{-1}\right)}\right]\right),$ wherein, at the moment, φ.sub.−1=ϕ(x, y)−ϕ(x−Δ, y) represents the shearing phase of 0.sup.th-order beam and −1.sup.st-order diffracted beam in they axis direction, and φ.sub.1=ϕ(x+Δ, y)−ϕ(x, y) represents the shearing phase of 1.sup.st-order diffracted beam and 0.sup.th-order beam in the y axis direction; (iv) defining a shear ratio s between +1.sup.st-order diffracted beam and 0.sup.th-order beam; multiplying the shearing phase φ.sub.x in the x axis direction by 2 to obtain a differential wavefront with a shear ratio of 2s in the x axis direction ΔW.sub.x=2φ.sub.x=ϕ(x+Δ, y)−ϕ(x−Δ, y); multiplying the shearing phase φ.sub.y in the y axis direction by 2 to obtain a differential wavefront with a shear ratio of 2s in the y axis direction ΔW.sub.y=2φ.sub.y=ϕ(x, y+Δ)−ϕ(x, y−Δ); and obtaining a wavefront aberration W of the optical imaging system under test (3) by using wavefront reconstruction algorithm of shearing interferometer on ΔW.sub.x and ΔW.sub.y according to a shear ratio of 2s, wherein the grating shearing interferometer system comprise the illuminating system (8), the object plane diffraction grating plate (1), the first three-dimensional stage (2), the image plane diffraction grating plate (4), the second three-dimensional stage (5), the two-dimensional photoelectric sensor (6), and a calculation processing unit (7); the illuminating system (8) outputs spatially incoherent light; the object plane diffraction grating plate (1) is fixed on the first three-dimensional stage (2); the image plane diffraction grating plate (4) is fixed on the second three-dimensional stage (5); the object plane diffraction grating plate (1) comprises two sets of one-dimensional gratings with vertical grating line directions; the image plane diffraction grating plate (4) comprises one set of checkerboard grating or two groups of one-dimensional grating with vertical grating line directions; an output end of the two-dimensional photoelectric sensor (6) is connected with the calculation processing unit (7), and a xyz coordinate system is established where z axis direction is along the optical axis direction of the grating shearing interferometer, the x axis is along a grating line direction of the second grating (102) on the object plane diffraction grating plate (1), the y axis is along the grating line direction of the first grating (101) on the object plane diffraction grating plate (1), and establishing that movement axes of the first three-dimensional stage (2) and the second three-dimensional stage (5) are respectively the x axis, the y axis, and the z axis.

2. The method according to claim 1, wherein a ratio of a period of the one-dimensional grating on the object plane diffraction grating plate (1) to a period of the checkerboard grating or the one-dimensional grating on the image plane diffraction grating plate (4) is equal to a magnification of the optical imaging system under test (3).

3. The method according to claim 1, wherein the duty-cycle of the object plane diffraction grating and the image plane diffraction grating is 50%.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 shows the wavefront measurement device of grating shearing interferometer in the present invention.

(2) FIG. 2 shows the object plane diffraction grating plate in the present invention.

(3) FIG. 3 shows the checkerboard grating of the image plane diffraction grating plate in the present invention.

(4) FIG. 4 is a diagram showing the relationship between the shear ratio s and the numerical aperture N.A. of the grating shearing interferometer system in the present invention.

(5) Reference numbers used in the figures refer to the following structure: 1—object plane diffraction grating plate; 2—first three-dimensional stage; 3—optical imaging system under test; 4—image plane diffraction grating plate; 5—second three-dimensional stage; 6—two-dimensional photoelectric sensor; 7—calculation processing unit; 8—illuminating system.

DETAILED DESCRIPTION OF THE INVENTION

(6) The present invention is further described in detail embodiments in connection with the drawings, however, the scope of the present invention should not be limited by the embodiments.

(7) The present invention discloses a method for high-accuracy wavefront measurement based on grating shearing interferometry that adopts a grating shearing interferometer as shown in FIG. 1. The device of the present invention comprises an illuminating system 8, an object plane diffraction grating plate 1, a first three-dimensional stage 2, an image plane diffraction grating plate 4, a second three-dimensional stage 5, a two-dimensional photoelectric sensor 6 and a calculation processing unit 7, wherein the illuminating system 8 outputs spatially incoherent light, the object plane diffraction grating plate 1 is fixed on the first three-dimensional stage 2, the image plane diffraction grating plate 4 is fixed on the second three-dimensional stage 5, and the output end of the two-dimensional photoelectric sensor 6 is connected with the calculation processing unit 7.

(8) In the present invention, an xyz coordinate system is established, wherein the z axis direction is along the optical axis direction of the shearing interferometer, the x axis is along the grating line direction of the second grating 102 on the object plane diffraction grating plate 1, the y axis is along the grating line direction of the first grating 101 on the object plane diffraction grating plate 1, and the motion axes of the first three-dimensional stage 2 and the second three-dimensional stage 5 are respectively an x axis, a y axis and a z axis; the first three-dimensional stage 2 is used for moving the first grating 101 and the second grating 102 in the object plane diffraction grating plate 1 to the center of the object plane field of the optical imaging system under test 3; the second three-dimensional stage 5 is used for moving the checkerboard grating in the image plane diffraction grating plate 4 to the center of the image plane field of the optical imaging system under test 3, and carrying out specific periodic movement in the x axis direction and the y axis direction on the image plane diffraction grating plate 4; the two-dimensional photoelectric sensor 6 can be a charge coupled device (CCD) or a CMOS image sensor, and the test plane receives the shearing interference interferogram generated by the diffraction of a checkerboard grating; and the calculation processing unit 7 is used for acquiring and storing shearing interferograms, then processing and analyzing the shearing interferograms.

(9) As shown in FIG. 2, the object plane diffraction grating plate 1 comprises two one-dimensional diffraction gratings, which are a first grating 101 with a grating line along the y axis direction and a second grating 102 with a grating line along the x axis direction respectively, wherein the period of the one-dimensional diffraction grating is P1 and the duty ratio is 50%. The first grating 101 and the second grating 102 can be phase gratings or amplitude gratings.

(10) As shown in FIG. 3, in the chessboard grating of the image plane diffraction grating plate 4, the chessboard grating is with a period of P2 and a duty ratio of 50%; the chessboard grating consists of square grids, and the diagonal direction of the square is along the x axis direction or the y axis direction; the chessboard grating can be a phase grating or an amplitude grating; the period P1 of the one-dimensional grating and the period P2 of the two-dimensional grating satisfy:
P1=M.Math.P2  (1)
wherein M is the magnification of the optical imaging system under test 3.

(11) FIG. 4 shows the relationship between the shear ratio s of the shearing interferometer system and the numerical aperture NA of the optical imaging system under test 3 and the period P2 of the image plane grating. The shear ratio is defined as the ratio of the diffraction angle to the full aperture angle (normalized shear ratio between 1.sup.st-order diffracted beam and 0.sup.th-order beam):

(12) $\begin{array}{cc}s=\frac{\beta }{\alpha }=\frac{\arcsin \left(\frac{\lambda }{P2}\right)}{2ars\sin \left(\mathrm{NA}\right)}& \left(2\right)\end{array}$

(13) wherein, β is the diffraction angle of 1.sup.st-order diffracted beam, and α is the full aperture angle of the beam.

(14) Based on the grating shearing interferometer system, only high-order diffraction beams with odd-order are generated after passing through the chessboard grating of the image plane diffraction grating plate 4, and the spatial coherence of the optical field is modulated by combining with the linear gratings 101 or 102 on the object plane diffraction grating plate 1, so that only the interference between 0.sup.th-order beam and other high-order diffracted beams exists in the optical field, and there is no interference between these high-order diffracted beams. Taking the shearing direction along the x direction as an example, the interference field received by the two-dimensional photoelectric sensor 6 can be written as:

(15) $\begin{array}{cc}\begin{array}{c}I\left(x,y\right)=A_0^2+\underset{m=-\left(2n+1\right)}{\overset{m=2n+1}{.Math.}}{A}_n^2+\mathrm{.Math.}+\\ 2\underset{m=1}{\overset{m=2n+1}{.Math.}}{A}_0{A}_m{\gamma }_m\cos \left[\varphi \left(x+m\Delta ,y\right)-\varphi \left(x,y\right)-{\alpha }_m\right]+\\ 2\underset{m=1}{\overset{m=2n+1}{.Math.}}{A}_0{A}_m{\gamma }_m\cos \left[\varphi \left(x,y\right)-\varphi \left(x-m\Delta ,y\right)-{\alpha }_m\right]\end{array}(n=0,1,2,\mathrm{.Math.})& \left(3\right)\end{array}$

(16) In the equation, A.sub.0 is the amplitude of the 0.sup.th-order beam, A.sub.m is the amplitude of the m.sup.th-order diffracted beam, is the optical path difference between the m.sup.th-order diffracted beam and the 0.sup.th-order beam, ϕ(x, y) is the wavefront to be measured, γ.sub.m is the degree of coherence between the m.sup.th-order diffracted beam and 0.sup.th-order beam, in is the diffraction order, and A is the offset of 1.sup.st-order diffracted beam relative to the 0.sup.th-order beam received by the two-dimensional photoelectric detector 6. Assuming A.sub.0 is 1, the A.sub.m and γ.sub.m coefficients satisfy the following relationship:

(17) $\begin{array}{cc}{A}_m={\gamma }_m=\frac{2}{m\pi }& \left(4\right)\end{array}$

(18) When introducing a phase-shifting, the above equation can be written as:

(19) $\begin{array}{cc}I=I0+\underset{m=1}{\overset{m=2n-1}{.Math.}}{a}_m\left(\cos \left({\phi }_{-m}+m\delta \right)+\cos \left({\phi }_m+m\delta \right)\right)& \left(5\right)\end{array}$

(20) Wherein φ.sub.−m=ϕ(x, y)−ϕ(x−mΔ, y) represents the shearing phase of 0.sup.th-order beam and −m.sup.th-order diffracted beam; φ.sub.m=ϕ(x+mΔ, y)−ϕ(x, y) indicates the shearing phase of m.sup.th-order diffracted beam and 0.sup.th-order beam.

(21) The method for high-accuracy wavefront measurement using the above described grating shearing interferometer of the present invention comprises the following steps:

(22) (1) placing the optical imaging system under test 3 in the grating shearing interferometer to enable the illuminating system 8 to be positioned on the object space of the optical imaging system under test 3, and the image plane diffraction grating plate 4 to be positioned on the image space of the optical imaging system under test 3, adjusting the first three-dimensional stage 2 to enable the object plane diffraction grating plate 1 to be positioned on the object plane of the optical imaging system under test 3, and adjusting the second three-dimensional stage 5 to enable the image plane diffraction grating plate 4 to be positioned on the image plane of the optical imaging system under test 3;

(23) (2) determining the phase shifting amounts according to the shear ratio s of the grating shearing interferometer: firstly, determining the maximum diffraction order in the interference field to be

(24) $m=\mathrm{ceil}\left(\frac{1}{s}\right)-1,$
with the diffraction orders of the grating shearing interferometer system being sequentially ±1, ±3, . . . , ±(2n−1), wherein n is the total number of positive high-order diffracted beam or the total number of negative high-order diffracted beams in the shearing interferometer system,

(25) $n=\mathrm{fix}\left(\frac{m+1}{2}\right),$
function ceil (X) returns a minimum integer greater than or equal to X, and function fix (X) returns a maximum integer less than or equal to X; then according to n, determining the movement of the image plane diffraction grating plant 4 when an shearing interferogram acquired being 0, ¼, ½, ¾ and

(26) $\frac{{\alpha }_i}{2\pi },\frac{\pi -{\alpha }_i}{2\pi },\frac{2\pi -{\alpha }_i}{2\pi }$
times the chessboard grating period (where i=2, 3 . . . n, α.sub.i<π and

(27) 0${\alpha }_i\ne \frac{\pi }{2}),$
the total number of steps being 3n+1;

(28) (3) moving the first three-dimensional stage 2 to enable the first grating 101 of the grating line on the object plane diffraction grating plate 1 along the y axis direction to move into the object space field-of-view point position of the optical imaging system under test 3, and moving the second three-dimensional stage 5 to enable the chessboard grating or a one-dimensional grating in a corresponding direction on the image plane diffraction grating plate 4 to move into an image space field-of-view point position of the imaging system under test 3, wherein an included angle between a diagonal line of the chessboard grating and the x axis (or the y axis) is 45 degrees;

(29) (4) moving the second three-dimensional stage 5 along the x axis direction according to the period 0, ¼, ½, ¾ and

(30) $\frac{{\alpha }_i}{2\pi },\frac{\pi -{\alpha }_i}{2\pi },\frac{2\pi -{\alpha }_i}{2\pi }$
times the chessboard grating period (wherein i=2, 3 . . . n, α.sub.i<π and

(31) ${\alpha }_i\ne \frac{\pi }{2}),$
wherein the two-dimensional photoelectric sensor 6 acquires one shearing interferogram I.sub.α.sub.i after each movement and transmits the shearing interferograms to the data processing unit 7 such that a total of 3n+1 shearing interferograms are obtained, an equation set consisting of n linear equation in the x direction is obtained, then the shearing phase is calculated according to the following method:

(32) $\begin{array}{cc}[\begin{array}{cccc}1& 1& \mathrm{.Math.}& 1\\ \cos {\alpha }_2& \cos 3{\alpha }_2& \mathrm{.Math.}& \cos \left(2n-1\right){\alpha }_2\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \cos {\alpha }_n& \cos 3{\alpha }_n& \mathrm{.Math.}& \cos \left(2n-1\right){\alpha }_n\end{array}].\hspace{1em}[\begin{array}{c}2a_1\left(\cos {\phi }_1+\cos {\phi }_{-1}\right)\\ 2a_3\left(\cos {\phi }_3+\cos {\phi }_{-3}\right)\\ \mathrm{.Math.}\\ 2a_{2\mathrm{n­}1}\left(\cos {\phi }_{2n­1}+\cos {\phi }_{-(2n­1)}\right)\end{array}]=\hspace{1em}\hspace{1em}\left[\begin{array}{c}{I}_0-I_{\pi }\\ I_{{\alpha }_2}-I_{{\mathrm{\pi ­\alpha }}_2}\\ \mathrm{.Math.}\\ I_{{\alpha }_n}-I_{{\mathrm{\pi ­\alpha }}_n}\end{array}\right]& \left(7\right)\\ \left[\begin{array}{cccc}\sin \frac{\pi }{2}& \sin \frac{3\pi }{2}& \mathrm{.Math.}& \sin \frac{(2n-1)\pi }{2}\\ \sin {\alpha }_2& \sin 3{\alpha }_2& \mathrm{.Math.}& \sin \left(2n-1\right){\alpha }_2\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ \sin {\alpha }_n& \sin 3{\alpha }_n& \mathrm{.Math.}& \sin \left(2n-1\right){\alpha }_n\end{array}\right].\hspace{1em}\left[\begin{array}{c}2a_1\left(\sin {\phi }_1+\sin {\phi }_{-1}\right)\\ 2a_3\left(\sin {\phi }_3+\sin {\phi }_{-3}\right)\\ \mathrm{.Math.}\\ 2a_{2n-1}\left(\sin {\phi }_{2n-1}+\sin {\phi }_{-\left(2n-1\right)}\right)\end{array}\right]=\hspace{1em}\hspace{1em}\left[\begin{array}{c}{I}_{\frac{\pi }{2}}-I_{\frac{3\pi }{2}}\\ I_{{\alpha }_2}-I_{2{\mathrm{\pi ­\alpha }}_2}\\ \mathrm{.Math.}\\ I_{{\alpha }_n}-I_{2{\mathrm{\pi ­\alpha }}_n}\end{array}\right]& \left(8\right)\end{array}$
2a.sub.1(cos φ.sub.1+ cos φ.sub.−1) and 2a.sub.1(sin φ.sub.1+ sin φ.sub.−1) and in the x direction are obtained by solving the linear equation sets (7) and (8) respectively, and the shearing phase φ.sub.x in the x axis direction is obtained by solving and unwrapping:

(33) $\begin{array}{cc}{\phi }_x=\mathrm{unwrap}\left(\arctan \left[\frac{2a_1\left(\sin {\phi }_1+\sin {\phi }_{-1}\right)}{2a_1\left(\cos {\phi }_1+\cos {\phi }_{-1}\right)}\right]\right),& \left(9\right)\end{array}$
at the moment, φ.sub.−1=ϕ(x, y)−ϕ(x−Δ, y) represents the shearing phase of 0.sup.th-order beam and −1.sup.st-order diffracted beam in the x axis direction, and φ.sub.1=ϕ(x+Δ, y)−ϕ(x, y) represents the shearing phase of 1.sup.st-order diffracted beam and 0.sup.th-order beam in the x axis direction; and the function unwrap(x) represents unwrapping x;

(34) (5) moving the first three-dimensional stage 2 to enable the second grating 102 of the grating line on the object plane diffraction grating plate 1 along the x axis direction to move into the object space field-of-view point position of the optical imaging system under test 3; performing the same movement on the second three-dimensional stage 5 along the y axis direction, acquiring one shearing interferogram I.sub.α.sub.i by the two-dimensional photoelectric sensor 6 after each movement and transmitting the shearing interferogram to the data processing unit 7, obtaining a total of 3n+1 shearing interference interferograms, obtaining an equation set consists of n linear equations in the y direction, calculating 2a.sub.1(cos φ.sub.1+ cos φ.sub.−1) and 2a.sub.1 (sin φ.sub.1+ sin φ.sub.−1) in the y direction according to the equation (7) and the equation (8) as well, and solving and unwrapping to obtain a shearing phase φ.sub.y in the y axis direction:

(35) $\begin{array}{cc}{\phi }_y=\mathrm{unwrap}\left(\arctan \left[\frac{2a_1\left(\sin {\phi }_1+\sin {\phi }_{-1}\right)}{2a_1\left(\cos {\phi }_1+\cos {\phi }_{-1}\right)}\right]\right),& \left(10\right)\end{array}$
at the moment, φ.sub.−1=ϕ(x, y)−ϕ(x−Δ, y) representing the shearing phase of 0.sup.th-order beam and −1.sup.st-order diffracted beam in the y axis direction, and φ.sub.1=ϕ(x+Δ, y)−ϕ(x, y) representing the shearing phase of 1.sup.st-order diffracted beam and 0.sup.th-order beam in the y axis direction; and

(36) (6) defining a shear ratio s between +1.sup.st-order diffracted beam and 0.sup.th-order beam;

(37) multiplying the shearing phase φ.sub.x in the x axis direction by 2 to obtain a differential wavefront with a shear ratio of 2s in the x axis direction ΔW.sub.x=2φ.sub.x=ϕ(x+Δ, y)−ϕ(x−Δ, y); multiplying the shearing phase φ.sub.y in the y axis direction by 2 to obtain a differential wavefront with a shear ratio of 2s in the y axis direction ΔW.sub.y=2φ.sub.y=ϕ(x, y+Δ)−ϕ(x, y−Δ); and obtaining the wavefront aberration W of the optical imaging system under test 3 by wavefront reconstruction algorithm on ΔW.sub.x and ΔW.sub.y according to the shear ratio of 2s.

(38) Additional description of step 6 is as follows:

(39) deducing the shearing phase in the x axis direction is as follows:

(40) $\begin{array}{c}{\phi }_x=\arctan \left[\frac{2a_1\left(\sin {\phi }_1+\sin {\phi }_{-1}\right)}{2a_1\left(\cos {\phi }_1+\cos {\phi }_{-1}\right)}\right]\\ =\arctan \left[\frac{2a_{1}2\sin \frac{{\phi }_1+{\phi }_{-1}}{2}\cos \frac{{\phi }_1-{\phi }_{-1}}{2}}{2a_{1}2\cos \frac{{\phi }_1+{\phi }_{-1}}{2}\cos \frac{{\phi }_1-{\phi }_{-1}}{2}}\right]\\ =\frac{{\phi }_1+{\phi }_{-1}}{2}\\ =\frac{1}{2}\{[\varphi \left(x+\Delta ,y\right)-\varphi \left(\left(x,y\right)\right]+\left[\varphi \left(x,y\right)-\varphi \left(\left(x-\Delta ,y\right)\right]\right\}\\ =\frac{\varphi \left(x+\Delta ,y\right)-\varphi \left(x-\Delta ,y\right)}{2}\end{array}$
wherein it can be seen that the shearing phase φ.sub.x in x direction is equivalent to half of the shearing phase of 1.sup.st-order diffracted beam and −1.sup.st-order diffracted beam in the physical sense, and the shear ratio is 2s;

(41) in the same way, the shearing phase in the y axis direction can be represented as

(42) ${\phi }_y=\frac{\varphi \left(x,y+\Delta \right)-\varphi \left(x,y-\Delta \right)}{2},$
which is also equivalent to half of the shearing phase of 1.sup.st-order diffracted beam and −1.sup.st-order diffracted beam, with a shear ratio of 2s, where s is defined as the shear ratio between 1.sup.st-order diffracted beam and 0.sup.th-order beam.

(43) In Chinese Invention Patent Applications 201910183242.1 and 201910183243.6 to Lu, Yunjun et al., as, approximate processing φ.sub.x≈φ.sub.1≈φ.sub.−1 is performed in the shearing phase extraction, the shearing phase φ.sub.x is approximately equal to the shearing phase of 1.sup.st-order diffracted beam and 0.sup.th-order beam in the physical sense, and the shear ratio is s; in the same way, the shearing phase φ.sub.y is approximately equal to the shearing phase of 1.sup.st-order diffracted beam and 0.sup.th-order beam, and the shear ratio is s. The differential phase retrieval error resulting from the approximation is

(44) $\frac{{\phi }_1-{\phi }_{-1}}{2},$
and the inaccuracy in the shear ratio s during the wavefront reconstruction also results in wavefront error. In general, the larger the shear ratio s, the larger the error resulting from the φ.sub.1≈φ.sub.−1 approximation processing, and the larger the final measured wavefront aberration error.

(45) The method of the present invention has the advantages of a large range of the measurable numerical aperture and the adjustable shear ratio of the grating interferometer, as well as the capability of accurately measuring the wavefront aberration of the optical imaging system under test, and the like.

(46) Compared with Chinese Invention Patent Applications 201910183242.1 and 201910183243.6 to Lu, Yunjun et al, the present invention adopts the differential wavefront completely equivalent to 1.sup.st-order diffracted beam and −1.sup.st-order diffracted beam, performs wavefront reconstruction strictly according to the shear ratio of 2s, obtains the original wavefront aberration from the shearing phase retrieval to the wavefront reconstruction, obtains the wavefront aberration accurately in the whole process without any approximation, as a result, improves the accuracy of wavefront aberration measurement of the optical imaging system under test. The method is optimized just by optimization of shearing phase retrieval algorithm with no need of extra mechanical structure, the interference of high-order diffracted beams are eliminated, as a result, the accuracy of wavefront aberration measurement is improved.