CN113155279A - Single-point measurement method for laser beam quality - Google Patents

Single-point measurement method for laser beam quality Download PDF

Info

Publication number
CN113155279A
CN113155279A CN202110625100.3A CN202110625100A CN113155279A CN 113155279 A CN113155279 A CN 113155279A CN 202110625100 A CN202110625100 A CN 202110625100A CN 113155279 A CN113155279 A CN 113155279A
Authority
CN
China
Prior art keywords
laser beam
formula
quality
laser
light
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.)
Pending
Application number
CN202110625100.3A
Other languages
Chinese (zh)
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.)
Changchun University of Science and Technology
Original Assignee
Changchun University of Science and Technology
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 Changchun University of Science and Technology filed Critical Changchun University of Science and Technology
Publication of CN113155279A publication Critical patent/CN113155279A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00Photometry, e.g. photographic exposure meter

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Testing Of Optical Devices Or Fibers (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention discloses a single-point measurement method for laser beam quality, which obtains x/y direction differential phase and light intensity information by transforming a four-wave transverse shearing interference pattern, reconstructs the phase information to be measured by adopting a differential zernike algorithm, reconstructs complex amplitude according to the phase information and the light intensity information, obtains the complex amplitude and the light spot size omega (z) of different axial positions by respectively adopting an angular spectrum algorithm and a second-order moment method, and further fits a laser transmission curve equation, thereby obtaining a far field divergence angle theta (theta) of a laser beam0Beam waist spot size omega0And a beam quality factor M2. The measuring method can realize real-time monitoring of the quality of the laser beam, is also suitable for measuring the quality of the pulse laser beam, and has the advantages of simplicity, high measuring efficiency and the like.

Description

Single-point measurement method for laser beam quality
Technical Field
The invention relates to a single-point measurement method for the beam quality of a laser, and belongs to the technical field of laser.
Background
The laser has been widely used in the fields of national economy, national defense safety and the like due to the characteristics of good monochromaticity, good coherence, good directivity, high brightness and the like. The quality of the laser beam has been a focus of attention of laser development, production and application workers. A multi-spot method beam quality measuring instrument represented by Spiricon, Thorlabs, etc., in the united states, which samples and measures spot sizes at 10 or more different positions in the laser beam axial direction mainly by moving an image detector (such as a CCD camera, a CMOS camera) in the beam axial direction, and requires that at least half or more of the positions are located in the rayleigh range, and then obtains the beam quality factor of the laser beam by hyperbolic fitting, has been widely used in practice. The more axial position samples, the higher the fitting accuracy, and the higher the measurement accuracy. Because the image detector needs to be moved along the axial direction of the light beam during each sampling, the light beam needs to be ensured to be vertically incident to the detection surface of the image detector all the time in the whole measuring process before measurement, the measuring period is long, the light path is difficult to adjust, and the introduced measuring error is large, so that the real-time monitoring of the quality of the laser beam cannot be met, and the method is not suitable for measuring the quality of the pulse laser beam.
Disclosure of Invention
Aiming at the problems of large measurement error, long measurement period and difficult light path adjustment of the quality of laser beams by a multipoint method and the problem that the real-time quality of the laser beams cannot be metThe invention aims to provide a single-point measurement method for measuring the beam quality of a laser, which is used for acquiring a far-field divergence angle theta by a four-wave transverse shear interference method0Beam waist spot omega0Beam quality factor M2And the like.
As shown in fig. 1, the present invention provides a single-point measurement method for laser beam quality, which comprises the following steps:
acquiring a four-wave transverse shearing interference image with pixel number width W and pixel number length H; the image has no background noise and saturated pixels;
the four-wave transverse shear interferogram in the first step is acquired by: splitting a light beam of a laser to be detected by adopting a two-dimensional phase grating, blocking light beams of other diffraction orders by using an aperture diaphragm, allowing only four sub-light beams of the (plus or minus 1) orders to pass through, and focusing the four sub-light beams on a detection surface of an image sensor through a lens, thereby obtaining a four-wave transverse shearing interference pattern on the image sensor;
performing Fourier transform on the four-wave transverse shearing interference graph to obtain a frequency spectrum, wherein the Fourier transform formula is as formula (1), the four-wave transverse shearing interference graph in the frequency domain comprises a central zero frequency and four pairs of fundamental frequencies, a midpoint O of the central zero frequency is taken as a coordinate origin, a horizontal right direction is taken as an x-axis positive direction, a vertical upward direction is taken as a y-axis positive direction, and a rectangular coordinate system xOy is established;
Figure BDA0003100693450000011
step three, respectively screening out a positive level frequency spectrum and a zero level frequency spectrum in the x direction and the y direction by adopting a frequency domain filtering method, and performing inverse Fourier transform by using a formula (2) to obtain phase distribution in the x direction and the y direction so as to obtain differential phase information delta phi in the x direction and the y directiontAnd light intensity information It
Figure BDA0003100693450000012
Step four, utilizing the differential phase information delta phi in the x direction and the y direction obtained in the step threetPerforming least square fitting to calculate zernike coefficient, and reconstructing phase information phi of the measured laser beam by using the zernike coefficient and zernike polynomialt
Step five, according to the obtained light intensity information ItAnd phase information
Figure BDA0003100693450000021
Respectively reconstructing the complex amplitude E of the laser beam at the position z of 0 by using the formula (3) to the formula (5)0(x, y) complex amplitude conjugate E0(x, y) and intensity distribution I0(x,y);
Figure BDA0003100693450000022
Figure BDA0003100693450000023
I(x,y)=E(x,y)*E*(x,y) (5)
Sixthly, complex amplitude information E of the axial coordinates z at different positions of the laser beam is obtained according to a formula (3) and an angular spectrum algorithm formula (6)z(x,y);
Figure BDA0003100693450000024
Where λ is the wavelength of the laser, fxAnd fyIs a fourier transform pair in the frequency domain;
step seven, solving the light intensity distribution I of more than 10 different axial positions in the Rayleigh range according to the formulas (4) to (6)z(x,y);
Step eight, solving the spot size omega of the different axial positions of the light beams according to the formula (7)x(z) and ωy(z);
Figure BDA0003100693450000025
Step nine, the size omega of the light spotx(z)、ωy(z) substituting the axial position coordinate z into equation (8) and fitting the coefficients A, B and C by means of the least square method to obtain a spatial transmission curve equation of the X/Y-direction laser beam;
ωx(z)2=AxZ2+BxZ+Cx,ωy(z)2=AyZ2+ByZ+Cy (8)
Figure BDA0003100693450000026
Figure BDA0003100693450000027
Figure BDA0003100693450000028
step ten, respectively obtaining the far field divergence angle theta of the laser beam in the X/Y direction according to the formulas (9) to (11)0Beam waist spot size omega0And a beam quality factor M2
Has the advantages that: according to the laser beam quality single-point measurement method provided by the invention, a laser transmission characteristic curve equation can be fitted by detecting the vertical axis light spot at any position of the measured laser beam along the axial direction, so that the beam waist light spot size, the far field divergence angle and the light beam quality factor of the measured laser beam are solved, the real-time monitoring of the laser beam quality can be realized, the method is also suitable for measuring the pulse laser beam quality, and the method has the advantages of simplicity in measurement method, high measurement efficiency and the like.
Drawings
FIG. 1 is a flow chart of a single point measurement method for laser beam quality.
Detailed Description
Embodiment 1 a laser beam quality single point measurement method.
As shown in fig. 1, the present invention provides a single-point measurement method for laser beam quality, which comprises the following steps:
acquiring a four-wave transverse shearing interference image with pixel number width W and pixel number length H; the image has no background noise and saturated pixels;
the four-wave transverse shear interferogram in the first step is acquired by: splitting a light beam of a laser to be detected by adopting a two-dimensional phase grating, blocking light beams of other diffraction orders by using an aperture diaphragm, allowing only four sub-light beams of the (plus or minus 1) orders to pass through, and focusing the four sub-light beams on a detection surface of an image sensor through a lens, thereby obtaining a four-wave transverse shearing interference pattern on the image sensor;
performing Fourier transform on the four-wave transverse shearing interference graph to obtain a frequency spectrum, wherein the Fourier transform formula is as formula (1), the four-wave transverse shearing interference graph in the frequency domain comprises a central zero frequency and four pairs of fundamental frequencies, a midpoint O of the central zero frequency is taken as a coordinate origin, a horizontal right direction is taken as an x-axis positive direction, a vertical upward direction is taken as a y-axis positive direction, and a rectangular coordinate system xOy is established;
Figure BDA0003100693450000031
step three, respectively screening out a positive level frequency spectrum and a zero level frequency spectrum in the x direction and the y direction by adopting a frequency domain filtering method, and performing inverse Fourier transform by using a formula (2) to obtain phase distribution in the x direction and the y direction so as to obtain differential phase information delta phi in the x direction and the y directiontAnd light intensity information It
Figure BDA0003100693450000032
Step four, utilizing the differential phase information delta phi in the x direction and the y direction obtained in the step threetPerforming least square fitting to calculate Zernike coefficient, and reconstructing phase information phi of the measured laser beam by using the Zernike coefficient and Zernike polynomialt
Step five, according to the obtained light intensity information ItAnd phase information
Figure BDA0003100693450000033
Respectively reconstructing the complex amplitude E of the laser beam at the position z of 0 by using the formula (3) to the formula (5)0(x, y) complex amplitude conjugate E0(x, y) and intensity distribution I0(x,y);
Figure BDA0003100693450000034
Figure BDA0003100693450000035
I(x,y)=E(x,y)*E*(x,y) (5)
Sixthly, complex amplitude information E of the axial coordinates z at different positions of the laser beam is obtained according to a formula (3) and an angular spectrum algorithm formula (6)z(x,y);
Figure BDA0003100693450000036
Where λ is the wavelength of the laser, fxAnd fyIs a fourier transform pair in the frequency domain;
step seven, solving the light intensity distribution I of more than 10 different axial positions in the Rayleigh range according to the formulas (4) to (6)z(x,y);
Step eight, solving the spot size omega of the different axial positions of the light beams according to the formula (7)x(z) and ωy(z);
Figure BDA0003100693450000041
Step nine, the size omega of the light spotx(z)、ωy(z) substituting the axial position coordinate z into equation (8) and fitting the coefficients A, B and C by means of the least square method to obtain a spatial transmission curve equation of the X/Y-direction laser beam;
ωx(z)2=AxZ2+BxZ+Cx,ωy(z)2=AyZ2+ByZ+Cy (8)
Figure BDA0003100693450000042
Figure BDA0003100693450000043
Figure BDA0003100693450000044
step ten, respectively obtaining the far field divergence angle theta of the laser beam in the X/Y direction according to the formulas (9) to (11)0Beam waist spot size omega0And a beam quality factor M2

Claims (1)

1. A single-point measurement method for the quality of a laser beam is characterized by comprising the following steps:
acquiring a four-wave transverse shearing interference image with pixel number width W and pixel number length H; the image has no background noise and saturated pixels;
the four-wave transverse shear interferogram in the first step is acquired by: splitting a light beam of a laser to be detected by adopting a two-dimensional phase grating, blocking light beams of other diffraction orders by using an aperture diaphragm, allowing only four sub-light beams of the (plus or minus 1) orders to pass through, and focusing the four sub-light beams on a detection surface of an image sensor through a lens, thereby obtaining a four-wave transverse shearing interference pattern on the image sensor;
performing Fourier transform on the four-wave transverse shearing interference graph to obtain a frequency spectrum, wherein the Fourier transform formula is as formula (1), the four-wave transverse shearing interference graph in the frequency domain comprises a central zero frequency and four pairs of fundamental frequencies, a midpoint O of the central zero frequency is taken as a coordinate origin, a horizontal right direction is taken as an x-axis positive direction, a vertical upward direction is taken as a y-axis positive direction, and a rectangular coordinate system xOy is established;
Figure FDA0003100693440000011
step three, respectively screening out a positive level frequency spectrum and a zero level frequency spectrum in the x direction and the y direction by adopting a frequency domain filtering method, and performing inverse Fourier transform by using a formula (2) to obtain phase distribution in the x direction and the y direction so as to obtain differential phase information in the x direction and the y direction
Figure FDA0003100693440000017
And light intensity information It
Figure FDA0003100693440000012
Step four, utilizing the differential phase information in the x direction and the y direction obtained in the step three
Figure FDA0003100693440000018
Performing least square fitting to calculate zernike coefficient, and reconstructing phase information of the measured laser beam by using the zernike coefficient and zernike polynomial
Figure FDA0003100693440000019
Step five, according to the obtained light intensity information ItAnd phase information
Figure FDA00031006934400000110
Respectively reconstructing the complex amplitude E of the laser beam at the position z of 0 by using the formula (3) to the formula (5)0(x, y) complex amplitude conjugate E0(x, y) and intensity distribution I0(x,y);
Figure FDA0003100693440000013
Figure FDA0003100693440000014
I(x,y)=E(x,y)*E*(x,y) (5)
Sixthly, complex amplitude information E of the axial coordinates z at different positions of the laser beam is obtained according to a formula (3) and an angular spectrum algorithm formula (6)z(x,y);
Figure FDA0003100693440000015
Where λ is the wavelength of the laser, fxAnd fyIs a fourier transform pair in the frequency domain;
step seven, solving the light intensity distribution I of more than 10 different axial positions in the Rayleigh range according to the formulas (4) to (6)z(x,y);
Step eight, solving the spot size omega of the different axial positions of the light beams according to the formula (7)x(z) and ωy(z);
Figure FDA0003100693440000016
Step nine, the size omega of the light spotx(z)、ωy(z) substituting the axial position coordinate z into equation (8) and fitting the coefficients A, B and C by means of the least square method to obtain a spatial transmission curve equation of the X/Y-direction laser beam;
Figure FDA0003100693440000021
Figure FDA0003100693440000022
Figure FDA0003100693440000023
Figure FDA0003100693440000024
step ten, respectively obtaining the far field divergence angle theta of the laser beam in the X/Y direction according to the formulas (9) to (11)0Beam waist spot size omega0And a beam quality factor M2
CN202110625100.3A 2020-10-30 2021-06-04 Single-point measurement method for laser beam quality Pending CN113155279A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202011190725.3A CN112326021A (en) 2020-10-30 2020-10-30 Single-point measurement method for laser beam performance parameters
CN2020111907253 2020-10-30

Publications (1)

Publication Number Publication Date
CN113155279A true CN113155279A (en) 2021-07-23

Family

ID=74296311

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202011190725.3A Pending CN112326021A (en) 2020-10-30 2020-10-30 Single-point measurement method for laser beam performance parameters
CN202110625100.3A Pending CN113155279A (en) 2020-10-30 2021-06-04 Single-point measurement method for laser beam quality

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN202011190725.3A Pending CN112326021A (en) 2020-10-30 2020-10-30 Single-point measurement method for laser beam performance parameters

Country Status (1)

Country Link
CN (2) CN112326021A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113639967A (en) * 2021-08-09 2021-11-12 南京理工大学 Laser beam quality real-time measuring method based on special-shaped prism

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113639967A (en) * 2021-08-09 2021-11-12 南京理工大学 Laser beam quality real-time measuring method based on special-shaped prism
CN113639967B (en) * 2021-08-09 2024-06-07 南京理工大学 Laser beam quality real-time measurement method based on special prism

Also Published As

Publication number Publication date
CN112326021A (en) 2021-02-05

Similar Documents

Publication Publication Date Title
CN104482875A (en) Single-slit spatial carrier shearing speckle interferometry measuring system and measuring method
CN105928455B (en) The coaxial striking rope type synchronous phase shift interferometer of space light splitting and its measurement method
US6229913B1 (en) Apparatus and methods for determining the three-dimensional shape of an object using active illumination and relative blurring in two-images due to defocus
JP6622515B2 (en) Lens array, wavefront sensor, wavefront measuring device, shape measuring device, aberration measuring device, optical element manufacturing method, and optical device manufacturing method
CN104006765A (en) Phase extraction method and detecting device for single width carrier frequency interference fringes
CN104320598A (en) Aberration-free image reconstruction method based on CCD array pixel response function frequency domain calibration
US9927371B2 (en) Confocal line inspection optical system
CN111121675A (en) Visual field expansion method for microsphere surface microscopic interferometry
EP3183535A1 (en) Interferometer
CN109186956A (en) A kind of transient state wavefront distortion measurement method for the online situation of beam shaping element
CN113155279A (en) Single-point measurement method for laser beam quality
Podanchuk et al. Holographic wavefront sensor based on the Talbot effect
CN116907804B (en) Single-pulse ultrafast laser space-time coupling distortion measuring instrument
CN110927116B (en) Method, device and system for measuring mark structure
WO2018146456A1 (en) Compact interferometer
CN103698022A (en) Wavefront measurement method of lateral shear interferometer
CN108107020B (en) Device and method for measuring nonlinear refractive index coefficient of material
CN109780993B (en) Differential phase contrast microscopy system and method
CN107589542B (en) The selection method of central wavelength in broadband phase difference image reconstruction
CN110608677B (en) Three-dimensional displacement measurement method based on circular grating radial shearing interferometer
CN113218519B (en) Radial shear wavefront measurement system based on double-layer sub-wave slot structure
CN115560848A (en) Intensity-wave front-wavelength measuring instrument based on super-surface radial shear
US20160320485A1 (en) Image processing method
CN114758019A (en) Infrared detector simulation image generation method and device
JP2011232243A (en) Measuring device

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