CN115290655B - Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking - Google Patents

Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking Download PDF

Info

Publication number
CN115290655B
CN115290655B CN202210823537.2A CN202210823537A CN115290655B CN 115290655 B CN115290655 B CN 115290655B CN 202210823537 A CN202210823537 A CN 202210823537A CN 115290655 B CN115290655 B CN 115290655B
Authority
CN
China
Prior art keywords
laser
image
thermal
output end
defect
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.)
Active
Application number
CN202210823537.2A
Other languages
Chinese (zh)
Other versions
CN115290655A (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.)
Harbin Institute of Technology
Wuhu Robot Technology Research Institute of Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Wuhu Robot Technology Research Institute of Harbin Institute of 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 Harbin Institute of Technology, Wuhu Robot Technology Research Institute of Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202210823537.2A priority Critical patent/CN115290655B/en
Publication of CN115290655A publication Critical patent/CN115290655A/en
Application granted granted Critical
Publication of CN115290655B publication Critical patent/CN115290655B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • G01N21/8851Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00Investigating or analyzing materials by the use of thermal means
    • G01N25/72Investigating presence of flaws
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/246Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/84Systems specially adapted for particular applications
    • G01N21/88Investigating the presence of flaws or contamination
    • G01N21/8851Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges
    • G01N2021/8887Scan or image signal processing specially adapted therefor, e.g. for scan signal adjustment, for detecting different kinds of defects, for compensating for structures, markings, edges based on image processing techniques
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Quality & Reliability (AREA)
  • Signal Processing (AREA)
  • Multimedia (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The invention provides a defect accurate detection photo-thermal fusion imaging device and method based on heat flow diffusion tracking. The invention provides a discrete cosine-Rucaskarad method based on thermal wave time-frequency characteristics, which is used for accurately tracking the time-frequency domain change process of the thermal wave characteristics and constructing the plane size of a defect by adopting a reverse reconstruction method. Finally, the defect size detection error of the composite material, the metal material and the like is less than 2.5 percent.

Description

Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking
Technical Field
The invention belongs to the technical fields of photothermal science and detection technology, signal analysis and feature extraction, image processing and identification, and particularly relates to a device and a method for precisely detecting photo-thermal fusion imaging based on heat flow diffusion tracking. The imaging device and the imaging method are suitable for the fields of precise nondestructive detection and evaluation of defects/damages of materials such as aerospace, microelectronics, micro-nano structures and the like.
Background
As an emerging nondestructive testing technology, the infrared thermal wave imaging testing technology has the advantages of non-contact, intuitiveness, large detection area and the like, and is suitable for nondestructive testing, quality evaluation and the like of thin metal components, composite material components and high polymer materials. The infrared thermal wave imaging detection technology adopts an external excitation source to carry out active thermal excitation loading on a detected test piece, heat flow permeates from the surface of the detected test piece to the inside, and when defects exist in the detected test piece, the defects are obviously different from the detected test piece material in the thermal physical parameters, so that the heat flow diffusion process is accelerated or hindered. At the moment, the blocking effect of the defects is diffused to the surface of the detected test piece through heat conduction, and the infrared radiation difference of the detected test piece can be identified by utilizing the thermal infrared imager, so that the effective identification of the defects is realized. The propagation of the heat flow completely depends on the temperature gradient in the detected test piece, and the heat flow diffuses from a high-temperature area to a low-temperature area, so that the amplification detection effect on the defects exists, and the transverse heat diffusion of the heat flow inevitably leads to the defect that the detected size is larger than the actual size. Although the thermal wave imaging detection method can realize the amplification detection of small-scale defects by using a non-uniform diffusion detection means of heat flow, the defect size detection deviation is large, so that the accurate defect detection based on the infrared thermal wave imaging detection technology faces a great challenge, and if the transverse thermal diffusion effect of the heat flow is effectively inhibited in the defect detection process, the defect detection method becomes a key technical bottleneck for realizing the accurate defect detection.
In the prior art, wei Jiacheng, guizhou university, et al, proposed a similar optical flow method (Likeness optical flow, LOF) that can estimate the thermal flow gradient diffusion field in two-dimensional corresponding thermal signature images and utilize the method to reduce the impact of lateral thermal diffusion on detection accuracy. Simulation and experiments prove that the LOF method plays a good role in inhibiting the transverse heat spreading effect of the heat flow of the infrared image sequence, and meanwhile, the method has a good characteristic enhancement effect for a specific characteristic extraction method, but in the LOF iterative operation process, strict criteria are not formulated for the selection of an initial image and a termination image, so that the accuracy of an obtained result is not significant when the defect scale is unknown. The LOF method is optimized by Harbin university of industry Wang Fei and the like, a global smooth assumption is added on the basis of an optical flow basic constraint equation based on a Horn-Schink algorithm with gradient change, the completeness of the equation is realized, defect characteristics are extracted through a characteristic extraction algorithm, defect depth information is obtained through theoretical calculation, the depth information is fed back to a characteristic image iteration model, and finally effective suppression of the defects on transverse heat diffusion is realized. However, since the Horn-Schuk algorithm is applied on the premise that the image moves slowly (i.e. the heat flow spreads slowly), the defect equivalent diameter accuracy still has a large deviation from the actual defect diameter. In the method, in order to realize accurate measurement of the debonding defect of the thermal barrier coating structure, the femtosecond laser is adopted to carry out active thermal excitation loading on a test piece, so that the transverse thermal diffusion time of heat flow is shortened, and the accurate measurement of the defect size is realized, but the method also brings the problem of low signal-to-noise ratio of defect detection, and can not be recognized for smaller debonding defects. The university of capital and university of Bejing and Feng Lichun and the like (patent number ZL200510077750. X) of New technology Co., ltd.) propose a method for reconstructing an infrared thermal wave detection tomographic image, in which three-dimensional reconstruction of defects is realized only according to a theoretical formula of diffusion characteristics of heat flow in the depth direction, and the planar diffusion amplification effect of the heat flow is not considered, so that the obtained tomographic image is pseudo-tomographic image only in the depth direction, and the obtained equivalent diameter of the defects has larger deviation from the actual diameter. In summary, the current research on infrared thermal wave imaging detection methods mainly focuses on how to three-dimensionally locate defects, and because of the specificity of heat flow gradient diffusion, there is a great challenge to precisely quantitatively characterize the sizes of defects, so that the research on the aspect is generally less.
In order to fully consider the heat flow gradient diffusion characteristic and further accurately quantify and characterize the defect size, the invention relates to a defect accurate detection photo-thermal fusion imaging device and method based on heat flow diffusion tracking.
Disclosure of Invention
The invention aims to solve the problems in the prior art and provides a photo-thermal fusion imaging device and a photo-thermal fusion imaging method for accurately detecting defects based on heat flow diffusion tracking. The invention is suitable for nondestructive testing of composite materials, metal materials and the like in the fields of aerospace, microelectronics, micro-nano structures and the like.
The invention is realized by the following technical scheme that the invention provides a defect accurate detection photo-thermal fusion imaging device based on heat flow diffusion tracking, which comprises a thermal infrared imager, a first BNC control line, a second BNC control line, a first Ethernet line, a conducting optical fiber, an air-cooling refrigerator, a semiconductor refrigerator, a laser, a control line, a laser power supply, a third BNC control line, a function generator, a first USB communication line, a fourth BNC control line, a lock-in amplifier, a fifth BNC control line, a data acquisition card, a second Ethernet line, an optical camera, a computer, a second USB communication line, a motion controller, a motion control line, a two-dimensional mobile station, a cylindrical lens, a first polarizer, an engineering diffuser, a sleeve, a collimating mirror, a clamp, a test piece and a second polarizer;
the computer is provided with four signal ends, one signal input end of the computer is connected with the signal output end of the thermal infrared imager through a first Ethernet wire, the second signal input end of the computer is connected with the signal output end of the optical camera through a second Ethernet wire, the third signal output end of the computer is connected with the input end of the function generator through a first USB communication wire, the fourth signal output end of the computer is connected with the input end of the motion controller through a second USB communication wire, one signal input end of the thermal infrared imager is connected with the output end of the function generator through a first BNC control wire, the second signal output end of the thermal infrared imager is connected with the input end of the data acquisition card through a second BNC control wire, the air cooling refrigerator, the semiconductor refrigerator and the laser are mechanically connected with the collimating mirror through a conducting optical fiber, the signal output end of the laser is connected with the collimating mirror through a cylindrical lens, the first polarizer, the engineering body, the sleeve and the collimating mirror through threads are connected on the sleeve, the signal output end of the laser power supply is connected with the input end of the laser through a control wire, the signal output end of the laser power supply is connected with the input end of the motion controller through the BNC control wire, the signal output end of the laser power supply is connected with the output end of the thermal infrared imager through the second BNC control wire, the signal output end of the thermal infrared imager is connected with the optical fiber through the output end of the motion controller through the second BNC control wire, the signal output end of the thermal infrared imager is connected with the optical fiber is connected with the output end of the optical fiber through the optical fiber.
The invention provides an imaging method for precisely detecting defects of a photo-thermal fusion imaging device based on heat flow diffusion tracking, which specifically comprises the following steps:
step one: defining a test piece to be detected, and placing the test piece on a two-dimensional mobile station;
step two: starting a defect accurate detection photo-thermal fusion imaging device based on heat flow diffusion tracking;
step three: after the laser power supply is turned on for a plurality of minutes, the laser, the air cooling refrigerator and the semiconductor refrigerator already control the temperature of the laser to be at a preset temperature, and the laser is turned on at the moment;
step four: starting the thermal infrared imager and the optical camera, imaging the test sample in real time, irradiating the test sample by using the laser, controlling the motion controller by the computer, and further adjusting the two-dimensional moving table to ensure that the laser beam irradiation positions are all in the fields of view of the thermal infrared imager and the optical camera;
step five: the computer controls the function generator to generate a pulse trigger signal, one path of the pulse trigger signal controls the laser to generate constant pulse laser with fixed power, the other path of the pulse trigger signal controls the infrared thermal imager to acquire a thermal radiation image of a test piece according to fixed frequency, meanwhile, the computer controls the two-dimensional moving table and the optical camera to work, the moving direction of the two-dimensional moving table is perpendicular to the direction of the line laser, the positions of the infrared thermal imager, the laser and the optical camera are relatively static, the test piece moves at a uniform speed along with a certain direction of the two-dimensional moving table, the position of the line scanning laser acting on a test field is fixed, and image sequence reconstruction is carried out after scanning is completed;
step six: the reconstructed image sequence is equivalent to the thermal image response stimulated by the pulse heat source, and the reconstructed image is subjected to discrete cosine change to obtain the characteristic response image sequence;
step seven: characteristic image sequences of different frequency components can be obtained through discrete cosine transformation, and a corresponding relation curve of phase characteristic differences of the defect position and the defect-free position of the test piece and different frequencies is selected at the moment;
step eight: based on two characteristic images, the motion velocity v of the image along the x and y directions is obtained by adopting the method of Rucaskaner x ,v y
Step nine: acquiring a final characteristic image;
step ten: registering and fusing the acquired characteristic image and the image acquired by the optical camera, and acquiring a characteristic value between pixel gaps of the original infrared characteristic image by adopting a fitting method;
step eleven: based on the obtained photo-thermal fusion imaging detection result, performing binarization threshold segmentation and connected domain marking to realize accurate defect detection;
step twelve: after the test is finished, 5 minutes later, the laser power supply, the laser, the function generator, the data acquisition card, the lock-in amplifier, the optical camera and the thermal infrared imager are turned off.
Further, in step two, the computer, the laser power supply, the data acquisition card, the function generator and the lock-in amplifier are turned on.
Further, in step three, the laser power was turned on for 2 minutes and the laser temperature was controlled at 15 ℃.
Further, in the fourth step, the polarization direction of the polarizer arranged at the front end of the laser is consistent with that of the polarizer in front of the thermal infrared imager, so that effective filtering of signals is achieved, the engineering diffuser arranged at the front end of the laser converts a Gaussian light source into uniform light, and the cylindrical lens converts the uniform light source into a linear laser light source.
Further, in the fifth step, the Pixel of Δpixel_x column of the line laser heat source in the field of view perpendicular to the moving direction thereof is taken, and defined as the number of collected pixels for reconstruction, and during the scanning process, the Pixel of Δpixel_x column will sweep the surface of the detected object along with the heat source and collect and record according to a certain frequency, and in order to enable the Pixel of Δpixel_x column to record the heat wave signal of the whole sweeping range, the following relation needs to be satisfied:
Figure SMS_1
wherein v is the moving speed of the two-dimensional mobile station, PL is the actual physical position corresponding to a single pixel, and f is the acquisition frequency of the thermal infrared imager;
after the line laser completely scans the tested test piece, the image sequence is reconstructed by using the following formula,
Figure SMS_2
wherein T is seqi The method is characterized in that an ith frame image acquired by a thermal infrared imager is T seq1 1 The reconstructed first image.
Further, in the sixth step, the discrete cosine transform after reconstructing the image is specifically:
Figure SMS_3
Figure SMS_4
wherein C is u [T seq (t n )]The discrete cosine transform after reconstructing the image, a (u) is an orthogonal transform coefficient, and N is the number of reconstructed image frames.
Further, in step seven, the phase characteristic difference is defined as,
DiPH=|Ph Nondefect {C u [T seq (t n )]}-Ph Defect {C u [T seq (t n )]}| (4)
wherein DipH is the difference in phase characteristics, ph Nondefect Mean value of phase characteristics of defect-free position, ph Defect The phase characteristic mean value of the defect position; selecting the phase characteristic images of the frequencies corresponding to the DipH peak position and the 1/2 position as P respectively peak And P peak/2
Further, in step eight, the rate of change of the two directions is obtained:
Figure SMS_5
Figure SMS_6
further, the finally acquired feature image is:
Ph=Ph peak/2 -Ph peak ·v x -Ph peak ·v y (6)。
the beneficial effects of the invention are as follows:
(1) The invention discloses a heat flow diffusion tracking-based photo-thermal fusion imaging method for accurately detecting defects, and provides a discrete cosine-Rucaskard method based on heat wave time-frequency characteristics. Finally, the defect size detection error is less than 2.5%.
(2) The invention provides a multi-source multi-sensing data fusion imaging detection method based on an infrared thermal imager and a monocular optical camera, which is used for carrying out high-precision fusion on a characteristic image and an optical image obtained based on infrared thermal wave imaging, so that the defect detection resolution of an infrared thermal wave imaging detection technology can be improved, and finally, the defect positioning precision error is less than 0.1%.
(3) The invention discloses a defect accurate detection photo-thermal fusion imaging device based on heat flow diffusion tracking, which adopts a lock-in amplifier to monitor real-time synchronization of heat flow modulation signals and infrared thermal imager acquisition signals, and feeds back synchronization delay information to a function generator to realize accurate control of the signals.
Drawings
FIG. 1 is a line laser scanning schematic;
FIG. 2 is a photo-thermal fusion imaging device for accurately detecting defects based on heat flow diffusion tracking;
FIG. 3 is a schematic diagram of the detection result.
In the figure: the infrared imaging system comprises a 1-infrared thermal imager, a 2-first BNC control wire, a 3-second BNC control wire, a 4-first Ethernet wire, a 5-conducting optical fiber, a 6-air cooling refrigerator, a 7-semiconductor refrigerator, an 8-laser, a 9-control wire, a 10-laser power supply, a 11-third BNC control wire, a 12-function generator, a 13-first USB communication wire, a 14-fourth BNC control wire, a 15-lock-in amplifier, a 16-fifth BNC control wire, a 17-data acquisition card, a 18-second Ethernet wire, a 19-optical camera, a 20-computer, a 21-second USB communication wire, a 22-motion controller, a 23-motion control wire, a 24-two-dimensional moving table, a 25-cylindrical lens, a 26-first polarizer, a 27-engineering diffuser, a 28-sleeve, a 29-collimating mirror, a 30-clamp, a 31-test piece and a 32-second polarizer.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention aims to provide a defect accurate detection photo-thermal fusion imaging method based on heat flow diffusion tracking, which aims to solve the problem of low defect detection precision caused by heat flow transverse heat diffusion in the current infrared heat wave imaging detection process. Finally, the defect size detection error of the composite material, the metal material and the like is less than 2.5 percent.
The invention provides a multi-source multi-sensor fusion imaging detection method based on a thermal infrared imager and a monocular optical camera, which aims to solve the problem of lower defect position positioning precision in the detection process by simply relying on the thermal infrared imager at present. At present, the pixels of the thermal infrared imager are generally lower (640 x 512 pixels are the common thermal infrared imager), the resolution of the optical camera is generally higher, and the characteristic image and the optical image acquired based on the thermal infrared imaging are fused with high precision, so that the defect detection resolution of the thermal infrared imaging detection technology can be improved, and finally, the defect positioning precision error of a composite material, a metal material and the like can be less than 0.1%.
Referring to fig. 1-3, the invention provides a defect accurate detection photo-thermal fusion imaging device based on heat flow diffusion tracking, which comprises a thermal infrared imager, a first BNC control line, a second BNC control line, a first ethernet line, a conducting optical fiber, an air-cooled refrigerator, a semiconductor refrigerator, a laser, a control line, a laser power supply, a third BNC control line, a function generator, a first USB communication line, a fourth BNC control line, a lock-in amplifier, a fifth BNC control line, a data acquisition card, a second ethernet line, an optical camera, a computer, a second USB communication line, a motion controller, a motion control line, a two-dimensional mobile station, a cylindrical lens, a first polarizer, an engineering diffuser, a sleeve, a collimating mirror, a clamp, a test piece and a second polarizer;
the computer is provided with four signal ends, one signal input end of the computer is connected with the signal output end of the thermal infrared imager through a first Ethernet wire, the second signal input end of the computer is connected with the signal output end of the optical camera through a second Ethernet wire, the third signal output end of the computer is connected with the input end of the function generator through a first USB communication wire, the fourth signal output end of the computer is connected with the input end of the motion controller through a second USB communication wire, one signal input end of the thermal infrared imager is connected with the output end of the function generator through a first BNC control wire, the second signal output end of the thermal infrared imager is connected with the input end of the data acquisition card through a second BNC control wire, the air cooling refrigerator, the semiconductor refrigerator and the laser are mechanically connected with the collimating mirror through a conducting optical fiber, the signal output end of the laser is connected with the collimating mirror through a cylindrical lens, the first polarizer, the engineering body, the sleeve and the collimating mirror through threads are connected on the sleeve, the signal output end of the laser power supply is connected with the input end of the laser through a control wire, the signal output end of the laser power supply is connected with the input end of the motion controller through the BNC control wire, the signal output end of the laser power supply is connected with the output end of the thermal infrared imager through the second BNC control wire, the signal output end of the thermal infrared imager is connected with the optical fiber through the output end of the motion controller through the second BNC control wire, the signal output end of the thermal infrared imager is connected with the optical fiber is connected with the output end of the optical fiber through the optical fiber.
The invention provides an imaging method for precisely detecting defects of a photo-thermal fusion imaging device based on heat flow diffusion tracking, which specifically comprises the following steps:
step one: defining a test piece to be detected, and placing the test piece on a two-dimensional mobile station;
step two: starting a defect accurate detection photo-thermal fusion imaging device based on heat flow diffusion tracking; in the second step, the computer, the laser power supply, the data acquisition card, the function generator and the lock-in amplifier are started.
Step three: after the laser power supply is turned on for a plurality of minutes, the laser, the air cooling refrigerator and the semiconductor refrigerator already control the temperature of the laser to be at a preset temperature, and the laser is turned on at the moment; in the third step, the laser power is turned on for 2 minutes, and the laser temperature is controlled at 15 ℃.
Step four: starting the thermal infrared imager and the optical camera, imaging the test sample in real time, irradiating the test sample by using the laser, controlling the motion controller by the computer, and further adjusting the two-dimensional moving table to ensure that the laser beam irradiation positions are all in the fields of view of the thermal infrared imager and the optical camera; in the fourth step, the polarization direction of the polaroid arranged at the front end of the laser is consistent with that of the polaroid in front of the thermal infrared imager, so that effective filtering of signals is realized, the engineering diffuser arranged at the front end of the laser converts a Gaussian light source into uniform light, and meanwhile, the cylindrical lens converts the uniform light source into a linear laser light source.
Step five: the computer controls the function generator to generate a pulse trigger signal, one path of the pulse trigger signal controls the laser to generate constant pulse laser with fixed power, the other path of the pulse trigger signal controls the infrared thermal imager to acquire a thermal radiation image of a test piece according to fixed frequency, meanwhile, the computer controls the two-dimensional moving table and the optical camera to work, the moving direction of the two-dimensional moving table is perpendicular to the direction of the line laser, the positions of the infrared thermal imager, the laser and the optical camera are relatively static, the test piece moves at a uniform speed along with a certain direction of the two-dimensional moving table, the position of the line scanning laser acting on a test field is fixed, and image sequence reconstruction is carried out after scanning is completed as shown in fig. 1;
in the fifth step, a Δpixel_x column Pixel perpendicular to the moving direction of the line-of-view laser heat source is taken, and is defined as the number of collected pixels for reconstruction, and in the scanning process, the Δpixel_x column Pixel will sweep the surface of the detected object along with the heat source and collect and record according to a certain frequency, and in order to enable the Δpixel_x column Pixel to record the heat wave signal of the whole sweeping range, the following relation needs to be satisfied:
Figure SMS_7
wherein v is the moving (scanning) speed of the two-dimensional mobile station, PL is the actual physical position corresponding to a single pixel, and f is the acquisition frequency of the thermal infrared imager;
after the line laser completely scans the tested test piece, the image sequence is reconstructed by using the following formula,
Figure SMS_8
wherein T is seqi The method is characterized in that an ith frame image acquired by a thermal infrared imager is T seq1 1 The reconstructed first image.
Step six: the reconstructed image sequence is equivalent to the thermal image response stimulated by the pulse heat source, and the reconstructed image is subjected to discrete cosine change to obtain the characteristic response image sequence;
in the sixth step, the discrete cosine transform after the image is reconstructed is specifically:
Figure SMS_9
Figure SMS_10
wherein C is u [T seq (t n )]The discrete cosine transform after reconstructing the image, a (u) is an orthogonal transform coefficient, and N is the number of reconstructed image frames.
Step seven: characteristic image sequences of different frequency components can be obtained through discrete cosine transformation, and a corresponding relation curve of phase characteristic differences of the defect position and the defect-free position of the test piece and different frequencies is selected at the moment;
in step seven, the phase characteristic difference is defined as,
DiPH=|Ph Nondefect {C u [T seq (t n )]}-Ph Defect {C u [T seq (t n )]}| (4)
wherein DipH is the difference in phase characteristics, ph Nondefect Mean value of phase characteristics of defect-free position, ph Defect The phase characteristic mean value of the defect position; selecting the phase characteristic images of the frequencies corresponding to the DipH peak position and the 1/2 position as P respectively peak And P peak/2
Step eight: based on two characteristic images, the motion velocity v of the image along the x and y directions is obtained by adopting the method of Rucaskaner x ,v y
In step eight, the rate of change in both directions is obtained:
Figure SMS_11
Figure SMS_12
step nine: acquiring a final characteristic image;
the finally acquired characteristic image is as follows:
Ph=Ph peak/2 -Ph peak ·v x -Ph peak ·v y (6)。
step ten: registering and fusing the acquired characteristic images with the images acquired by the optical camera, and acquiring characteristic values among pixel gaps of the original infrared characteristic images by adopting a fitting method to realize pixel expansion and improvement of image resolution;
step eleven: based on the obtained photo-thermal fusion imaging detection result, performing binarization threshold segmentation and connected domain marking to realize accurate defect detection;
step twelve: after the test is finished, 5 minutes later, the laser power supply, the laser, the function generator, the data acquisition card, the lock-in amplifier, the optical camera and the thermal infrared imager are turned off.
Examples
The computer 20 of the invention is provided with four signal ends, one signal input end of the computer 20 is connected with a signal output end of the thermal infrared imager 1 through a first Ethernet line 4, a second signal input end of the computer 20 is connected with a signal output end of an optical camera 19 through a second Ethernet line 18, a third signal output end of the computer 20 is connected with an input end of the function generator 12 through a first USB communication line 13, a fourth signal output end of the computer 20 is connected with an input end of the motion controller 22 through a second USB communication line 21, one signal input end of the thermal infrared imager 1 is connected with an output end of the function generator 12 through a first BNC control line 2, a second signal output end of the thermal infrared imager 1 is connected with an input end of the optical camera 19 through a second BNC control line 3, the air-cooled refrigerator 6, the semiconductor refrigerator 7 and the laser 8 are mechanically connected through bolts, the output end of the laser 8 is connected with the collimating lens 29 through the conducting optical fiber 5, the cylindrical lens 25, the first polaroid 26, the engineering diffuser 27, the sleeve 28 and the collimating lens 29 are connected on the sleeve 28 through threads, the signal output end of the laser power supply 10 is connected with the input end of the laser 8 through the control line 9, the signal input end of the laser power supply 10 is connected with the signal output end of the function generator 12 through the third BNC control line 11, the signal output end of the function generator 12 is connected with the signal input end of the lock-in amplifier 15 through the fourth BNC control line 14, the output end of the data acquisition card 17 is connected with the signal input end of the lock-in amplifier 15 through the fifth BNC control line 16, the output end of the motion controller 22 is connected with the signal input end of the two-dimensional moving table 24 through the motion control line 23, the clamp 30 is fixed on the two-dimensional moving table 24 through threads, the clamp 30 is used for clamping a test piece 31, and the second polaroid 32 is connected with the thermal infrared imager 1 through threaded connection.
According to the precise defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking shown in fig. 2, in the embodiment, the thermal infrared imager 1 is of the model of FLIR SC6520, the response wavelength is 3.6-5.2 μm, the pixel size is 320×256, the maximum frame frequency is 100Hz, the optical camera pixel is 1920×1080 pixels, and the test piece 31 is a carbon fiber composite material with a prefabricated flat bottom hole for simulating a debonding defect.
The invention provides an imaging method for precisely detecting defects of a photo-thermal fusion imaging device based on heat flow diffusion tracking, which specifically comprises the following steps:
step one: defining a test piece 31 to be detected, and placing the test piece 31 on a two-dimensional mobile station;
step two: starting a defect accurate detection photo-thermal fusion imaging device based on heat flow diffusion tracking, wherein the method comprises the steps of starting a computer 20, a laser power supply 10, a data acquisition card 17, a function generator 12, a lock-in amplifier 15 and other devices;
step three: after the laser power supply 10 is turned on for about 2 minutes, the laser air-cooled refrigerator 6 and the semiconductor refrigerator 7 have controlled the temperature of the laser 8 to 15 ℃, and at this time, the laser 8 is turned on;
step four: the thermal infrared imager 1 and the optical camera 19 are started, real-time imaging is carried out on a test sample, the laser 8 is adopted to irradiate the test sample, the computer 20 controls the motion controller 22, the two-dimensional moving table 24 is further adjusted to ensure that the laser beam irradiation positions are all in the visual field of the thermal infrared imager 1 and the optical camera 19, the polarization direction of the polarization plate 26 arranged at the front end of the laser is consistent with that of the polarization plate 32 arranged in front of the thermal infrared imager 1, effective filtering of signals is achieved, the Gaussian light source is converted into uniform light by the engineering diffuser 27 arranged at the front end of the laser, and meanwhile, the uniform light source is converted into a linear laser light source by the cylindrical lens 25.
Step five: the computer 20 controls the function generator 12 to generate a pulse trigger signal, one path of the pulse trigger signal controls the laser 8 to generate constant pulse laser with fixed power, the other path of the pulse trigger signal controls the thermal imaging system 1 to acquire thermal radiation images of the test piece 31 according to fixed frequency, meanwhile, the computer 20 controls the two-dimensional moving table 24 and the optical camera 19 to work, wherein the moving direction of the two-dimensional moving table 24 is vertical to the line laser direction, the positions of the thermal imaging system 1, the laser 8 and the optical camera 19 are relatively static, the test piece 31 moves at a constant speed along with a certain direction of the two-dimensional moving table 24, and the position of the line scanning laser acting on a test field is fixed.
Taking the delta pixel_x column Pixel of the line-of-view laser heat source perpendicular to the moving direction of the line-of-view laser heat source, defining the Pixel as the collection Pixel number for reconstruction, wherein the delta pixel_x column Pixel can sweep the surface of a detected object along with the heat source in the scanning process and collect and record according to a certain frequency, and the following relational expression needs to be satisfied in order to enable the delta pixel_x column Pixel to record the heat wave signal of the whole sweeping range:
Figure SMS_13
where v is the moving (scanning) speed of the two-dimensional mobile station, PL is the actual physical position corresponding to a single pixel, and f is the acquisition frequency of the thermal infrared imager.
After the line laser completely scans the tested test piece, the image sequence is reconstructed by using the following formula,
Figure SMS_14
wherein T is seqi For the ith frame image acquired by the thermal infrared imager 1, T seq1 1 The reconstructed first image.
Step six: the reconstructed image sequence is equivalent to the thermal image response after being excited by the pulse heat source, the reconstructed image is utilized to obtain the characteristic response image sequence by discrete cosine change,
Figure SMS_15
Figure SMS_16
wherein C is u [T seq (t n )]The discrete cosine transform after reconstructing the image, a (u) is an orthogonal transform coefficient, and N is the number of reconstructed image frames.
Step seven: characteristic image sequences of different frequency components can be obtained through discrete cosine transformation, at the moment, phase characteristic differences of the defect position and the defect-free position of the test piece are selected to correspond to curves of different frequencies, wherein the phase characteristic differences are defined as,
DiPH=|Ph Nondefect {C u [T seq (t n )]}-Ph Defect {C u [T seq (t n )]}| (10)
wherein DipH is the difference in phase characteristics, ph Nondefect Mean value of phase characteristics of defect-free position, ph Defect Is the defect position phase characteristic mean value. Selecting the phase characteristic images of the frequencies corresponding to the DipH peak position and the 1/2 position as P respectively peak And P peak/2
Step eight: based on two characteristic images, the motion velocity v of the image along the x and y directions is obtained by adopting the method of Rucaskaner x ,v y Wherein, the method can describe the transformation process of two characteristic images, thereby obtaining the change rates of two directions,
Figure SMS_17
Figure SMS_18
step nine: the characteristic image which is finally acquired is that,
Ph=Ph peak/2 -Ph peak ·v x -Ph peak ·v y (12)
step ten: and registering and fusing the acquired characteristic image with the image acquired by the optical camera 19, and acquiring the characteristic value between the pixel gaps of the original infrared characteristic image by adopting a fitting method to realize pixel expansion and image resolution improvement.
Step eleven: and performing binarization threshold segmentation and connected domain marking based on the obtained photo-thermal fusion imaging detection result to realize accurate defect detection.
Step twelve: after the test is finished, 5 minutes after the test is finished, the laser power supply 10, the laser 8, the function generator 12, the data acquisition card 17, the lock-in amplifier 15, the optical camera 19, the thermal infrared imager 1 and other devices are turned off.
The test results for the test piece are shown in fig. 3, in which the error of detection is less than 2% for the three simulated flat bottom hole defect diameters.
The invention provides a defect accurate detection photo-thermal fusion imaging device and method based on heat flow diffusion tracking, which are described in detail, wherein specific examples are applied to illustrate the principle and implementation of the invention, and the description of the above examples is only used for helping to understand the method and core ideas of the invention; meanwhile, as those skilled in the art will have variations in the specific embodiments and application scope in accordance with the ideas of the present invention, the present description should not be construed as limiting the present invention in view of the above.

Claims (3)

1. The imaging method of the defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking is characterized in that the imaging device comprises a thermal infrared imager, a first BNC control line, a second BNC control line, a first Ethernet line, a conducting optical fiber, an air-cooled refrigerator, a semiconductor refrigerator, a laser, a control line, a laser power supply, a third BNC control line, a function generator, a first USB communication line, a fourth BNC control line, a lock-in amplifier, a fifth BNC control line, a data acquisition card, a second Ethernet line, an optical camera, a computer, a second USB communication line, a motion controller, a motion control line, a two-dimensional mobile station, a cylindrical lens, a first polarizer, an engineering diffuser, a sleeve, a collimating mirror, a clamp, a test piece and a second polarizer;
the computer is provided with four signal ends, one signal input end of the computer is connected with the signal output end of the thermal infrared imager through a first Ethernet wire, the second signal input end of the computer is connected with the signal output end of the optical camera through a second Ethernet wire, the third signal output end of the computer is connected with the input end of the function generator through a first USB communication wire, the fourth signal output end of the computer is connected with the input end of the motion controller through a second USB communication wire, one signal input end of the thermal infrared imager is connected with the output end of the function generator through a first BNC control wire, the second signal output end of the thermal infrared imager is connected with the input end of the data acquisition card through a second BNC control wire, the air cooling refrigerator, the semiconductor refrigerator and the laser are mechanically connected with the collimating mirror through a conducting optical fiber, the cylindrical lens, the first polarizer, the engineering body and the collimating mirror are connected on the sleeve through threads, the signal output end of the laser power supply is connected with the input end of the laser through a control wire, the signal output end of the laser power supply is connected with the input end of the laser through a second USB communication wire, the signal output end of the thermal infrared imager is connected with the output end of the laser through a second BNC control wire, the signal output end of the thermal infrared imager is connected with the output end of the optical fiber through a second BNC control wire, the optical fiber is connected with the output end of the optical fiber through a two-dimensional amplifier, and the optical fiber is connected with the output end of the optical fiber through the optical fiber;
the imaging method specifically comprises the following steps:
step one: defining a test piece to be detected, and placing the test piece on a two-dimensional mobile station;
step two: starting a defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking;
step three: after the laser power supply is turned on for a plurality of minutes, the laser, the air cooling refrigerator and the semiconductor refrigerator already control the temperature of the laser to be at a preset temperature, and the laser is turned on at the moment;
step four: starting the thermal infrared imager and the optical camera, imaging a test piece to be detected in real time, irradiating the test piece by adopting a laser, controlling the motion controller by a computer, and further adjusting the two-dimensional moving table to ensure that the irradiation positions of the laser beams are all in the fields of view of the thermal infrared imager and the optical camera;
step five: the computer controls the function generator to generate a pulse trigger signal, one path of the pulse trigger signal controls the laser to generate constant pulse laser with fixed power, the other path of the pulse trigger signal controls the infrared thermal imager to acquire a thermal radiation image of the test piece according to fixed frequency, meanwhile, the computer controls the two-dimensional moving table and the optical camera to work, the moving direction of the two-dimensional moving table is perpendicular to the direction of the line laser, the positions of the infrared thermal imager, the laser and the optical camera are relatively static, the test piece moves at a uniform speed along with a certain direction of the two-dimensional moving table, the position of the line scanning laser acting on a visual field is fixed, and image sequence reconstruction is carried out after scanning is completed;
step six: the reconstructed image sequence is equivalent to the thermal image response stimulated by the pulse heat source, and the reconstructed image is subjected to discrete cosine change to obtain the characteristic response image sequence;
step seven: characteristic response image sequences of different frequency components can be obtained through discrete cosine transformation, and a corresponding relation curve of phase characteristic differences of the defect position and the defect-free position of the test piece and different frequencies is selected at the moment;
step eight: based on the two characteristic response images, the motion velocity v of the image along the x and y directions is obtained by adopting the method of Rucaskaner x ,v y
Step nine: acquiring a final characteristic response image;
step ten: registering and fusing the acquired characteristic response image and the image acquired by the optical camera, and acquiring a characteristic value between pixel gaps of the original infrared characteristic image by adopting a fitting method;
step eleven: based on the obtained photo-thermal fusion imaging detection result, performing binarization threshold segmentation and connected domain marking to realize defect detection;
step twelve: after the test is finished, 5 minutes later, the power supply of the laser, the function generator, the data acquisition card, the lock-in amplifier, the optical camera and the thermal infrared imager are turned off;
in the fourth step, the polarization direction of a polaroid arranged at the front end of the laser is consistent with that of a polaroid in front of the thermal infrared imager, so that effective filtering of signals is realized, an engineering diffuser arranged at the front end of the laser converts a Gaussian light source into uniform light, and a cylindrical lens converts the uniform light source into a linear laser light source;
in step five, the DeltaPixel x column pixels, which are perpendicular to the movement direction of the laser heat source in the line of view, are defined as being used for
The number of reconstructed acquisition pixels, the Δpixel_x column pixels sweep the surface of the detected object along with the heat source in the scanning process and acquire and record according to a certain frequency, and in order to enable the Δpixel_x column pixels to record the heat wave signals of the whole sweeping range, the following relational expression needs to be satisfied:
Figure FDA0004203587600000022
wherein v is the moving speed of the two-dimensional mobile station, PL is the actual physical position corresponding to a single pixel, and f is the acquisition frequency of the thermal infrared imager;
after the line laser completely scans the tested test piece, the image sequence is reconstructed by using the following formula,
Figure FDA0004203587600000031
wherein T is seqi The method is characterized in that an ith frame image acquired by a thermal infrared imager is T seq1 1 A reconstructed first image;
in the sixth step, the discrete cosine transform after the image is reconstructed is specifically:
Figure FDA0004203587600000032
Figure FDA0004203587600000033
wherein C is u [T seq (t n )]Discrete cosine transform after reconstructing an image, a (u) is an orthogonal transform coefficient, and N is the number of frames of the reconstructed image;
in step seven, the phase characteristic difference is defined as,
DiPH=|Ph Nondefect {C u [T seq (t n )]}-Ph Defect {C u [T seq (t n )]}| (4)
wherein DipH is the difference in phase characteristics, ph Nondefect Mean value of phase characteristics of defect-free position, ph Defect The phase characteristic mean value of the defect position; selecting the phase characteristic response images of the frequencies corresponding to the DipH peak position and the 1/2 position as P respectively peak And P peak/2
In step eight, the rate of change in both directions is obtained:
Figure FDA0004203587600000034
Figure FDA0004203587600000035
the finally obtained characteristic response image is as follows:
Ph=Ph peak/2 -Ph peak ·v x -Ph peak ·v y (6)。
2. the imaging method of claim 1, wherein in step two, the computer, the laser power supply, the data acquisition card, the function generator, and the lock-in amplifier are turned on.
3. The imaging method of claim 2, wherein in step three, the laser is powered on for 2 minutes and the laser temperature is controlled at 15 ℃.
CN202210823537.2A 2022-07-13 2022-07-13 Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking Active CN115290655B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210823537.2A CN115290655B (en) 2022-07-13 2022-07-13 Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210823537.2A CN115290655B (en) 2022-07-13 2022-07-13 Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking

Publications (2)

Publication Number Publication Date
CN115290655A CN115290655A (en) 2022-11-04
CN115290655B true CN115290655B (en) 2023-07-11

Family

ID=83821525

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210823537.2A Active CN115290655B (en) 2022-07-13 2022-07-13 Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking

Country Status (1)

Country Link
CN (1) CN115290655B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117030733B (en) * 2023-07-07 2024-08-16 哈尔滨工业大学 Multimode laser excitation wide spectrum multi-source fusion imaging detection system and method based on FPGA
CN117782327B (en) * 2023-08-01 2024-08-13 哈尔滨工业大学 Composite photo-excitation pearson related infrared photothermal imaging method for helicopter rotor blade impact damage
CN117388316B (en) * 2023-12-11 2024-03-15 瑞金市金联机械有限公司 Raw material defect detection and marking equipment on aluminum film balloon machine

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234953B (en) * 2013-04-16 2015-03-11 南京诺威尔光电系统有限公司 Laser scanning thermal wave tomography system and method
CN106404835B (en) * 2016-09-13 2019-04-09 哈尔滨工业大学 A kind of infrared relevant thermal wave imaging system and the detection method based on the system
CN106768367A (en) * 2017-02-20 2017-05-31 哈尔滨工业大学 A kind of three-dimensional infrared thermal wave image-forming detecting system and method
CN106940332A (en) * 2017-03-31 2017-07-11 哈尔滨工业大学 A kind of Infrared Non-destructive Testing system and method based on phase code modulation system
CN109613063A (en) * 2019-01-17 2019-04-12 四川大学 A kind of device and method based on face battle array pulse laser excitation detection thermal barrier coating face
DE102019124344B4 (en) * 2019-09-11 2021-11-11 Edevis Gmbh Method and device for checking the internal condition of a test object by means of an impressed heat flow

Also Published As

Publication number Publication date
CN115290655A (en) 2022-11-04

Similar Documents

Publication Publication Date Title
CN115290655B (en) Imaging method of defect detection photo-thermal fusion imaging device based on heat flow diffusion tracking
CN103926274B (en) Infrared thermal wave radar imaging nondestructive testing method for defects of carbon fiber reinforced plastic (CFRP) plywood
CN108426839B (en) Additive manufacturing component detection method based on mechanical arm scanning laser ultrasonic signal correlation analysis
CN106768367A (en) A kind of three-dimensional infrared thermal wave image-forming detecting system and method
CN104040327A (en) Method for measuring thermal conductivity
CN108760546A (en) A kind of fatigue crack growth rate measurement method based on Infrared Thermography Technology
Welch et al. Remote measurement of in‐plane diffusivity components in plates
CN109990985B (en) Triangular linear array infrared detector modulation transfer function test method
CN111692991B (en) Point cloud data acquisition method for measuring batten bonding surface based on white light interference
CN108444921B (en) Additive manufacturing component online detection method based on signal correlation analysis
CN103644854A (en) Film thickness detection method based on laser scanning thermal wave imaging technology
CN109164136B (en) High-flux measurement system and method for thermoelectric transport parameters
CN113295092B (en) Line laser three-dimensional measurement system and method for thin-wall part
CN115112636B (en) Photo-thermal three-dimensional characteristic reconstruction system and reconstruction method based on phase control array laser beam concentration excitation
CN105444878A (en) High-precision mass measurement device and high-precision mass measurement method of chemical oxygen iodine laser far-field beam
CN113960032B (en) Online laser cleaning effect detection method and three-dimensional detection device
CN105092585A (en) Subsurface measuring device and method based on total internal reflection and optical coherence tomography
CN103673904A (en) Laser-scanning thermal wave imaging film thickness measuring method
CN210220974U (en) Automatic detection and positioning device for laser spot focus
CN113884538B (en) Infrared thermal image detection method for micro defects in large wind turbine blade
CN102155931A (en) Sub-surface damage detection method based on temperature field finite element analysis and simulation
CN113447527A (en) Dual-mode laser infrared thermal imaging detection system and method
Wang et al. Parameters impact analysis of CFRP defect detection system based on line laser scanning thermography
CN110987357B (en) Two-dimensional focusing laser differential interferometer and flat boundary layer density pulsation measurement method
CN113281002A (en) Flow field measuring system and measuring method of transient flow field under flow fluctuation working condition

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