CN110081988A - A method of space carrier frequency phase shift algorithm is used for the demodulation of four wave lateral shear interferometer wavefront slopes - Google Patents

A method of space carrier frequency phase shift algorithm is used for the demodulation of four wave lateral shear interferometer wavefront slopes Download PDF

Info

Publication number
CN110081988A
CN110081988A CN201910274361.8A CN201910274361A CN110081988A CN 110081988 A CN110081988 A CN 110081988A CN 201910274361 A CN201910274361 A CN 201910274361A CN 110081988 A CN110081988 A CN 110081988A
Authority
CN
China
Prior art keywords
light
phase
lateral shear
wavefront
wave
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910274361.8A
Other languages
Chinese (zh)
Other versions
CN110081988B (en
Inventor
陈小君
刘峰伟
吴永前
闫峰涛
赵彦
肖向海
徐燕
张娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN201910274361.8A priority Critical patent/CN110081988B/en
Publication of CN110081988A publication Critical patent/CN110081988A/en
Application granted granted Critical
Publication of CN110081988B publication Critical patent/CN110081988B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J9/02Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods
    • G01J9/0215Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by interferometric methods by shearing interferometric methods

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

The present invention relates to a kind of methods that space carrier frequency phase shift algorithm is used for the demodulation of four wave lateral shear interferometer wavefront slopes, belong to field of optical detection.Space carrier frequency phase shift algorithm is used for the demodulation of four wave lateral shear interferometer wavefront slopes by this method, finally to realize unified, high-precision dynamic wavefront measurement.Especially by the relationship between control screen periods and image sensor pixel size, to control the phase-shift phase between carrier fringe adjacent pixel, and the suitable phase shift algorithm of final choice realizes the demodulation of wavefront slope.Fast Fourier Transform (FFT) method (Fast Fourier Transform is utilized compared to traditional four wave lateral shear interferometers, FFT wavefront slope) is demodulated, this method calculating process is simple, high-precision unified dynamic wavefront measurement may be implemented in the interrupted spectra being related to there is no FFT method and edge effect.

Description

It is a kind of that space carrier frequency phase shift algorithm is used for four wave lateral shear interferometer wavefront slopes The method of demodulation
Technical field
The invention belongs to the manufacture of advanced optics and detection fields, and in particular to a kind of that space carrier frequency phase shift algorithm is used for four The method of wave lateral shear interferometer wavefront slope demodulation.
Background technique
Interferometer be using the wavelength of light as a kind of high-precision of measurement unit, highly sensitive metrology and measurement instrument, it It is widely used in optical measurement.Interference fringe is the track of the identical point of optical path difference in interferometer interference field, according to interference Shape, movement, density degree of striped etc. can demodulate to obtain the optical path difference information of tested light.
Four wave lateral shear interferometers are a kind of grating types by propositions such as French scientist Jerome Primo in 2000 Lateral shear interferometer.Light is tested through optical grating diffraction, 90% or more diffraction light energy concentrates on four beam first-order diffraction light.This four Beam first-order diffraction light is along respective diffraction angular direction onwards transmission, and finally mutual dislocation and coherent superposition on the viewing screen, form Interference fringe.By handling interference fringe, demodulation obtains fringe phase, can acquire the wavefront slope of tested light, most Optical path difference information is obtained using restoration algorithm afterwards.
Currently from there are mainly three types of the methods of interference fringe demodulation phase, respectively the Schlieren method, FFT (Fast Fourier Transform) method and phase-shifting method, the wherein precision highest of phase-shifting method, but its interference for needing to acquire several different phase-shift phases Figure, it is more sensitive to extraneous ambient vibration;The Schlieren method and FFT method only need a width interference pattern, influenced by ambient vibration it is smaller, It can be used for dynamically measuring, but processing accuracy is lower.In consideration of it, nineteen ninety Shough et al. proposes space carrier frequency phase shift (Spatial-carrier phase shifting SCPS) method, specific method are the angles by controlling two interfering beams The carrier phase amount for obtaining neighbor pixel is equal to the interference pattern of pi/2, and assumes phase phase of several consecutive points on tested surface Together, then adjacent several pixels can substitute the pixel obtained in time domain by phase shifter phase shift in space, so Phase can be calculated by phase shift formula afterwards.Although space carrier frequency phase shift method only handles single width interference pattern, its processing result The processing accuracy that phase-shifting method can be reached approximately, in combination with FFT method dynamic measurement and the high-precision advantage of phase-shifting method.
The method of four traditional wave lateral shear interferometer demodulated interferential fringe phases is FFT method, but due to side in FFT method The influence of the factors such as edge effect, interrupted spectra and spectral resolution limitation, causes its algorithm operation quantity larger and processing result precision It is affected.
Space carrier frequency phase shifting method is used for four wave lateral shear interferometers by the present invention, passes through control screen periods and image Relationship between sensor pixel size controls the phase-shift phase between carrier fringe neighbor pixel, is realized using phase shift algorithm Interference fringe phase, that is, wavefront slope high-precision is restored, final to realize the measurement of optical path difference Dynamic High-accuracy.
Summary of the invention
The purpose of the present invention is space carrier frequency phase shift method is used for four wave lateral shear interferometer interference fringe phases, i.e. quilt The Dynamic High-accuracy measurement of optical path difference is realized in the demodulation for surveying light wave front slope.This method passes through control screen periods and figure first As the relationship between sensor pixel size, so that the carrier phase amount between adjacent pixel is pi/2, phase shift algorithm is then utilized Solution obtains the wavefront slope of four direction.
The technical solution adopted by the present invention is that: it is a kind of that space carrier frequency phase shift algorithm is used for four wave lateral shear interferometer waves The method of front slope demodulation, comprising:
Step 1: four wave lateral shear interferometer systems of building, comprising:
Two-dimentional complex amplitude grating, for generating four beam first-order diffraction light, the two dimension complex amplitude grating is by amplitude grating and phase Position grating be formed by stacking, the amplitude grating period be p × p, light hole having a size ofPhase grating is by 0 phase and π phase It is alternately arranged, the period is 2p × 2p;
Imaging sensor, for acquiring interference fringe, Pixel Dimensions
Tested light is incident on two-dimentional complex amplitude grating, and through optical grating diffraction, most diffraction light energies concentrate on four beams one Grade diffraction light, this four beams diffraction light form interference fringe, utilize image in position of image sensor mutual dislocation and coherent superposition The interference fringe of sensor acquisition, i.e., the resilient optical path difference information for obtaining tested light;Change screen periods or selection is suitable Imaging sensor so that the carrier phase amount of interference fringe neighbor pixel be pi/2;
Step 2: the light intensity expression of four wave lateral shear interferometers are as follows:
Wherein I (x, y) is interference fringe light intensity, and a (x, y) is interference fringe background light intensity, and b (x, y) is modulation degree, s x The shearing displacement in direction or the direction y, k are wave number, and p is the amplitude grating period;
When assuming that when the corresponding background light intensity of 10 adjacent position points, modulation degree and identical tested light phase, solution point The formula of wavefront slope is as follows at (m, n):
Wherein, I0=I (m, n), I1=I (m-1, n), I2=I (m+1, n), I3=I (m, n-1), I4=I (m, n+1), I5 =I (m-1, n+1), I6=I (m+1, n-1), I7=I (m-1, n-1), I8=I (m+1, n+1), I9=I (m-2, n-2);
Wherein, W is tan before tested light wave-1For arctan function;
The wavefront slope of four direction is found out using formula (2)-(5), final restore obtains tested light optical path difference information.
The advantages of the present invention over the prior art are that:
(1) operation of the present invention is simple, and operand is low, can accelerate to restore speed, realize optical path difference transient measurement.
(2) marginal error of the present invention is small, and recovery accuracy is high, can almost realize unified measurement.
Detailed description of the invention
Fig. 1 is four wave lateral shear interferometer system schematics;
Fig. 2 is four wave lateral shear interferometer two dimension complex amplitude grating schematic diagrams;
Fig. 3 is the light intensity position that wavefront slope needs to use at solution point (m, n);
Fig. 4 is tested wavefront distribution schematic diagram;
Fig. 5 is the wavefront slope in the direction x extracted using the present invention;
Fig. 6 is the wavefront slope in the direction y extracted using the present invention;
Fig. 7 is the wavefront slope in the direction x+y extracted using the present invention;
Fig. 8 is the wavefront slope in the direction x-y extracted using the present invention;
Fig. 9 is that the wavefront restored using the wavefront slope that the present invention extracts is distributed.
Specific embodiment
With reference to the accompanying drawing and specific embodiment further illustrates the present invention.
A kind of method that space carrier frequency phase shift algorithm is used for the demodulation of four wave lateral shear interferometer wavefront slopes of the present invention, If Fig. 1 show four wave lateral shear interferometer system schematics, tested light is incident on two-dimentional complex amplitude grating, spreads out through grating It penetrates, most diffraction light energies concentrate on four beam first-order diffraction light, this four beams diffraction light is mutually wrong in position of image sensor Position and coherent superposition form interference fringe.The interference fringe acquired using imaging sensor, i.e., the resilient light for obtaining tested light Path difference information.
Two-dimentional complex amplitude grating (as shown in Figure 2), for generating four beam first-order diffraction light.The grating is by amplitude grating and phase Position grating be formed by stacking, the amplitude grating period be p × p, light hole having a size ofPhase grating is by 0 phase and π phase It is alternately arranged, the period is 2p × 2p;
Imaging sensor, for acquiring interference fringe, Pixel Dimensions
The light intensity expression of four wave lateral shear interferometers are as follows:
Wherein I (x, y) is interference fringe light intensity, and a (x, y) is interference fringe background light intensity, and b (x, y) is modulation degree, s x The shearing displacement in direction or the direction y, k are wave number, and p is the amplitude grating period.
When assuming that the corresponding background light intensity of 10 adjacent position points shown in Fig. 3, modulation degree and tested light phase are identical When, the formula of wavefront slope is as follows at solution point (m, n):
As shown in figure 3, I0=I (m, n), I1=I (m-1, n), I2=I (m+1, n), I3=I (m, n-1), I4=I (m, n+ 1), I5=I (m-1, n+1), I6=I (m+1, n-1), I7=I (m-1, n-1), I8=I (m+1, n+1), I9=I (m-2, n-2).
Wherein W is tan before tested light wave-1For arctan function.
Shown in the light intensity expression of four wave lateral shear interferometer interference fringes such as claim 2 Chinese style (1), wherein striped is carried on the back Scape light intensity a (x, y), the wavefront slope of modulation degree b (x, y) and four direction It is unknown.Therefore, in order to solve wavefront slope, six equations are at least needed.
For simplicity, it enables,
Then light intensity expression becomes:
After imaging sensor acquires, the light intensity expression of discretization are as follows:
Wherein m, n are the pixel coordinate in transverse and longitudinal direction, due toThereforeLight Strongly expressed formula is rewritten are as follows:
In order to acquire the wavefront slope at point (m, n), the light intensity using as shown in Figure 3 ten consecutive points is needed, and assume Their corresponding background light intensity, modulation degree and phase are identical, that is, assume:
Similarly, it is assumed that a (x, y), b (x, y),Respectively in this ten consecutive points It is worth identical.
As shown in figure 3, this corresponding light intensity expression of ten points is respectively as follows:
Therefore,
According to formula (6), then shown in the solution formula of wavefront slope such as formula (2)~(5).
When wavefront to be measured is defocus aberration shown in (such as Fig. 4), with the direction x acquired of the present invention, the direction y, the direction x+y and Respectively as shown in figures 5-8, the wavefront restored is as shown in Figure 9 for the wavefront slope in the direction x-y.As can be seen that the present invention proposes Method can accurately extract the wavefront slopes of four wave lateral shear interferometers, and finally realize that the high-precision of wavefront is restored.

Claims (1)

1. a kind of method that space carrier frequency phase shift algorithm is used for the demodulation of four wave lateral shear interferometer wavefront slopes, feature exist In, comprising:
Step 1: four wave lateral shear interferometer systems of building, comprising:
Two-dimentional complex amplitude grating, for generating four beam first-order diffraction light, the two dimension complex amplitude grating is by amplitude grating and phase light Grid are formed by stacking, the amplitude grating period be p × p, light hole having a size ofPhase grating is by 0 phase and π phase alternation It arranges, the period is 2p × 2p;
Imaging sensor, for acquiring interference fringe, Pixel Dimensions
Tested light is incident on two-dimentional complex amplitude grating, and through optical grating diffraction, most diffraction light energies concentrate on four beam level-ones and spread out Light is penetrated, this four beams diffraction light forms interference fringe, utilize image sensing in position of image sensor mutual dislocation and coherent superposition The interference fringe of device acquisition, i.e., the resilient optical path difference information for obtaining tested light;Change screen periods or the suitable figure of selection As sensor, so that the carrier phase amount of interference fringe neighbor pixel is pi/2;
Step 2: the light intensity expression of four wave lateral shear interferometers are as follows:
Wherein I (x, y) is interference fringe light intensity, and a (x, y) is interference fringe background light intensity, and b (x, y) is modulation degree, and s is the direction x Or the shearing displacement in the direction y, k are wave number, p is the amplitude grating period;
When assuming that when the corresponding background light intensity of 10 adjacent position points, modulation degree and identical tested light phase, solution point (m, n) The formula for locating wavefront slope is as follows:
Wherein, I0=I (m, n), I1=I (m-1, n), I2=I (m+1, n), I3=I (m, n-1), I4=I (m, n+1), I5=I (m-1, n+1), I6=I (m+1, n-1), I7=I (m-1, n-1), I8=I (m+1, n+1), I9=I (m-2, n-2);
Wherein, W is tan before tested light wave-1For arctan function;
The wavefront slope of four direction is found out using formula (2)-(5), final restore obtains tested light optical path difference information.
CN201910274361.8A 2019-04-08 2019-04-08 Method for applying spatial carrier frequency phase shift algorithm to wavefront slope demodulation of four-wave lateral shearing interferometer Active CN110081988B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910274361.8A CN110081988B (en) 2019-04-08 2019-04-08 Method for applying spatial carrier frequency phase shift algorithm to wavefront slope demodulation of four-wave lateral shearing interferometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910274361.8A CN110081988B (en) 2019-04-08 2019-04-08 Method for applying spatial carrier frequency phase shift algorithm to wavefront slope demodulation of four-wave lateral shearing interferometer

Publications (2)

Publication Number Publication Date
CN110081988A true CN110081988A (en) 2019-08-02
CN110081988B CN110081988B (en) 2021-09-21

Family

ID=67414384

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910274361.8A Active CN110081988B (en) 2019-04-08 2019-04-08 Method for applying spatial carrier frequency phase shift algorithm to wavefront slope demodulation of four-wave lateral shearing interferometer

Country Status (1)

Country Link
CN (1) CN110081988B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111772574A (en) * 2020-06-06 2020-10-16 宁波明星科技发展有限公司 Method and device for eliminating defocusing during wavefront aberration measurement
CN114137713A (en) * 2021-11-17 2022-03-04 华中科技大学 Real-time quantitative phase imaging method and system for unmarked thick sample

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4692027A (en) * 1985-07-03 1987-09-08 Itek Corporation Spatial processing for single or dual shear wavefront sensor
US20090009751A1 (en) * 2007-03-08 2009-01-08 Yukio Taniguchi Aberration measurement apparatus and aberration measurement method
CN101451890A (en) * 2008-12-18 2009-06-10 南京大学 Three light wave transversal shearing interference apparatus and method for extracting differential phase
US20110134438A1 (en) * 2009-12-07 2011-06-09 Canon Kabushiki Kaisha Refractive index distribution measuring method and refractive index distribution measuring apparatus
CN102252765A (en) * 2011-06-24 2011-11-23 北京理工大学 Wavefront aberration detection device and method for beam expansion and collimation system
CN102681365A (en) * 2012-05-18 2012-09-19 中国科学院光电技术研究所 Projection objective wave aberration detection device and method
CN102914374A (en) * 2011-08-03 2013-02-06 佳能株式会社 Wavefront measuring apparatus and wavefront measuring method
CN102967378A (en) * 2012-12-09 2013-03-13 中国科学院光电技术研究所 Miniaturized radial shearing interferometer based on four-step phase-shifting principle
CN103134439A (en) * 2013-01-29 2013-06-05 天津大学 Double-rotation optical wedge space phase shift method used for cutting speckle interference
CN104155011A (en) * 2014-07-22 2014-11-19 北京理工大学 Phase extracting method for two-dimension interferogram
CN204028565U (en) * 2014-08-26 2014-12-17 中国科学院上海光学精密机械研究所 A kind of projection objective wave aberration on-line checkingi interferometer
CN104655291A (en) * 2015-03-16 2015-05-27 中国科学院光电技术研究所 Method for realizing programmable multi-wave lateral shearing interferometer
CN107063477A (en) * 2017-03-16 2017-08-18 中国科学院上海光学精密机械研究所 Grating lateral shearing interference large-numerical aperture wave-front reconstruction method
CN109163816A (en) * 2018-08-28 2019-01-08 中国科学院光电技术研究所 Radial shearing interferometer based on cosine wave band sheet

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4692027A (en) * 1985-07-03 1987-09-08 Itek Corporation Spatial processing for single or dual shear wavefront sensor
US20090009751A1 (en) * 2007-03-08 2009-01-08 Yukio Taniguchi Aberration measurement apparatus and aberration measurement method
CN101451890A (en) * 2008-12-18 2009-06-10 南京大学 Three light wave transversal shearing interference apparatus and method for extracting differential phase
US20110134438A1 (en) * 2009-12-07 2011-06-09 Canon Kabushiki Kaisha Refractive index distribution measuring method and refractive index distribution measuring apparatus
CN102252765A (en) * 2011-06-24 2011-11-23 北京理工大学 Wavefront aberration detection device and method for beam expansion and collimation system
CN102914374A (en) * 2011-08-03 2013-02-06 佳能株式会社 Wavefront measuring apparatus and wavefront measuring method
CN102681365A (en) * 2012-05-18 2012-09-19 中国科学院光电技术研究所 Projection objective wave aberration detection device and method
CN102967378A (en) * 2012-12-09 2013-03-13 中国科学院光电技术研究所 Miniaturized radial shearing interferometer based on four-step phase-shifting principle
CN103134439A (en) * 2013-01-29 2013-06-05 天津大学 Double-rotation optical wedge space phase shift method used for cutting speckle interference
CN104155011A (en) * 2014-07-22 2014-11-19 北京理工大学 Phase extracting method for two-dimension interferogram
CN204028565U (en) * 2014-08-26 2014-12-17 中国科学院上海光学精密机械研究所 A kind of projection objective wave aberration on-line checkingi interferometer
CN104655291A (en) * 2015-03-16 2015-05-27 中国科学院光电技术研究所 Method for realizing programmable multi-wave lateral shearing interferometer
CN107063477A (en) * 2017-03-16 2017-08-18 中国科学院上海光学精密机械研究所 Grating lateral shearing interference large-numerical aperture wave-front reconstruction method
CN109163816A (en) * 2018-08-28 2019-01-08 中国科学院光电技术研究所 Radial shearing interferometer based on cosine wave band sheet

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈小君: "近场光强分布对四波横向剪切干涉仪波前复原的影响研究", 《中国博士学位论文全文数据库工程科技辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111772574A (en) * 2020-06-06 2020-10-16 宁波明星科技发展有限公司 Method and device for eliminating defocusing during wavefront aberration measurement
CN111772574B (en) * 2020-06-06 2023-06-27 宁波明星科技发展有限公司 Method and device for eliminating defocus during wavefront aberration measurement
CN114137713A (en) * 2021-11-17 2022-03-04 华中科技大学 Real-time quantitative phase imaging method and system for unmarked thick sample

Also Published As

Publication number Publication date
CN110081988B (en) 2021-09-21

Similar Documents

Publication Publication Date Title
Leendertz Interferometric displacement measurement on scattering surfaces utilizing speckle effect
Judge et al. Holographic deformation measurements by Fourier transform technique with automatic phase unwrapping
CN101769722B (en) Method for heterodyne temporal series speckle interferometry of object deformation
US20190011247A1 (en) Method for Phase Resolved Heterodyne Shearographic Measurements
CN101915559B (en) Method and system thereof for measuring three-dimensional surface shape of object by electronic speckle phase shift technology
CN103091676A (en) Mining area surface subsidence synthetic aperture radar interferometry monitoring and calculating method
CN102221342A (en) Method for measuring object deformation by time-domain multi-wavelength heterodyne speckle interference
CN104482877A (en) Motion compensation method and system in three-dimensional imaging of dynamic object
CN110081988A (en) A method of space carrier frequency phase shift algorithm is used for the demodulation of four wave lateral shear interferometer wavefront slopes
CN104034277A (en) Measuring method for phase of dual-wavelength micro-nano structure
CN102914256A (en) Synchronous phase shifting interference detection device based on orthogonal double grating and detection method
CN109059787A (en) A kind of sheet thickness distribution measurement method and system based on lateral shearing interference
CN105698702B (en) A kind of diplopore heterodyne ineterferometer based on acousto-optic low frequency differences phase shift
Flores et al. Color deflectometry for phase retrieval using phase-shifting methods
CN102865810A (en) Orthogonal double-grating based detecting device for synchronous phase shift common-light path interference and detecting method therefor
CN105300276A (en) Dual-wavelength single-exposure interference measuring method and system
CN111947600B (en) Robust three-dimensional phase unfolding method based on phase level cost filtering
CN103542870B (en) Ac modulation type low coherence interference demodulating system
CN103267485A (en) Point-diffraction three-dimensional absolute displacement measuring method
CN112212806A (en) Three-dimensional phase unfolding method based on phase information guidance
CN106091974A (en) A kind of object distortion measurement instrument, method and apparatus
CN106643474B (en) A kind of orthonormalization multistep phase-shifting phase measurement method
El-Morsy A novel algorithm based on sub-fringe integration method for direct two-dimensional unwrapping phase reconstruction from the intensity of one-shot two-beam interference fringes
CN201724658U (en) System for measuring three dimensional surface shape of an object by electronic speckle phase shift technology
CN102853761A (en) Space phase shifter

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant