Multi-source heterogeneous data fusion analysis method for target laser ablation effect experiment

Abstract

The invention particularly relates to a multi-source heterogeneous data fusion analysis method for a target laser ablation effect experiment. The technical problems that an existing laser ablation effect experiment is low in data utilization rate, insufficient in information mining and high in time cost, economic cost and labor cost are solved. According to the target material laser ablation effect experiment multi-source heterogeneous data fusion analysis method, target material front surface light spot distribution, temperature field and ablation hole data can be fused and analyzed, and information in target material laser ablation effect experiment data is fully excavated; the corresponding relation between the light intensity and the temperature and between the light intensity and the linear ablation rate on the outline of the ablation hole and the result of time evolution of statistical information of the size and the area of the ablation hole are obtained, and a large number of experiment rules are obtained through one target material laser ablation effect experiment.

Description

Technical Field

[0001] The present invention specifically relates to a method for fusion analysis of multi-source heterogeneous data of target material laser ablation effect experiments, which is suitable for fusion analysis of experimental data such as scattered light signals, spot online monitoring data, and thermal map sequences of target material laser ablation effect. Background Art

[0002] In applications such as laser processing and laser defense, the ablation perforation time of materials, the linear ablation rate and the change law of the incident laser intensity are important effect data that attract attention. Such effect data are usually obtained using laser ablation effect experimental methods.

[0003] Among them, the experimental data obtained from the laser ablation effect experiment usually includes the voltage curve output by the photoelectric detector detecting the scattered light on the front surface of the target and the rear protective plate surface, the thermal image sequence file output by the thermal imager monitoring the temperature field on the front surface of the target, and the spot measurement equipment monitoring the target surface and other optical paths. When processing this type of multi-source heterogeneous data, it is usually processed separately. Except for the statistical information such as perforation time, maximum light intensity / average light intensity, and maximum temperature, other information is not associated between the multi-source heterogeneous data.

[0004] In summary, in order to obtain sufficient information to understand the laws of ablation effect, a large number of laser ablation effect experiments need to be carried out. However, on the one hand, the data utilization efficiency of this experiment is low and the information mining is insufficient. On the other hand, the time, economic and labor costs are high. Summary of the invention

[0005] The purpose of the present invention is to solve the technical problems of low data utilization, insufficient information mining and high time, economic and labor costs in existing laser ablation effect experiments, and to provide a multi-source heterogeneous data fusion analysis method for target material laser ablation effect experiments.

[0006] In order to achieve the above object, the present invention adopts the following technical solution:

[0007] A method for fusion analysis of multi-source heterogeneous data in target laser ablation effect experiment, which is special in that it includes the following steps:

[0008] 1] Define the data acquisition timing signal V 时序 , taking the laser emission time as zero, adjust the data acquisition timing signal V 时序 The laser emission time t = 0;

[0009] 2] Using the thermal map sequence files of the target material front surface and the target material rear protection plate, the target material front surface temperature field, and the spot data files at the equal optical path of the target material front surface, calculate the laser light emission time, the target material perforation time, the thermal map frame number at the target material perforation time, the thermal map frame number at the laser light stop time, the spot frame number at the laser light emission time, and the spot frame number at the laser light stop time; then align the time axes corresponding to the scattered light on the target material front surface and the scattered light on the target material rear protection plate surface, the time axis of the thermal map sequence files of the target material front surface temperature field, and the time axis of the spot data files at the equal optical path of the target material front surface with the laser light emission time t=0, and set the thermal map frame number list {n} that needs to be processed;

[0010] 3] Calculate the coordinates of the inner point of the ablation hole, the contour coordinates, the centroid coordinates of the ablation hole, the size and area of ​​the ablation hole, and the spot frame number and the adjacent frame number interval in the spot data file corresponding to the current time [n] according to the thermal map sequence file, the spot data file and the thermal map frame number list {n} t1,S ,n t2,S ); read the adjacent frame number interval [n t1,S ,n t2,S ) corresponding to the image frame data I nt1 ,I nt2 , and interpolate it to get the spot matrix at the current moment, and calculate the centroid coordinate X of the spot according to the spot matrix C and the peak coordinate X P ;

[0011] 4] Set the centroid coordinate X C Or peak coordinate X P Overlap with the centroid coordinates of the ablation hole to obtain the row-column correspondence between the spot matrix and the temperature field matrix; obtain the corresponding contour of the ablation hole in the spot according to the row-column correspondence and the contour coordinates of the ablation hole, and calculate the light intensity on the corresponding contour of the ablation hole in the spot; then obtain the corresponding relationship between the light intensity and the temperature on the contour according to the temperature at the contour, and finally calculate the corresponding relationship between the light intensity and the linear ablation rate;

[0012] 5] Determine whether the thermal map frame number list {n} has reached the end point. If so, output the time series of the relationship between light intensity and temperature, light intensity and line ablation rate on the contour, and the results of the evolution of the size and area of ​​the ablation hole over time to complete the multi-source heterogeneous data fusion analysis of the target laser ablation effect experiment; otherwise, return to step 3].

[0013] Further, step 2] is specifically as follows:

[0014] 2.1. Calculate the laser emission time t using the target front surface and the target rear protection plate 激光 and target perforation time t 穿孔, and the time axis corresponding to the scattered light on the front surface of the target and the scattered light on the surface of the protective plate behind the target is connected to the data acquisition timing signal V 时序 Zero-time alignment of

[0015] 2.2. Processing the thermal map sequence file F of the target front surface temperature field 热图 , get the zero-pre-start time t of the heat map sequence 热图，零前起始 , thermal image frame number n at the time of target perforation 穿孔,T and the thermal image frame number n at the laser stop time 停光,T , the heat map sequence file F 热图 The time axis of is aligned with the zero time of the data acquisition timing signal, and a list of heat map frame numbers to be processed {n} is set; wherein the subscript T represents the heat map;

[0016] 2.3. Process the spot data file F at the equal optical path on the front surface of the target 光斑 , get the zero-front starting time t of the spot sequence 光斑，零前起始 , the spot frame number n at the time of laser emission 出光,S and the spot frame number n at the moment when the laser stops 停光,S , the spot data file F 光斑 The time axis of is aligned with the zero time of the data acquisition timing signal; wherein the subscript S represents the light spot.

[0017] Furthermore, step 2.1 is specifically as follows:

[0018] 2.1.1. Draw the scattered light output voltage curve V on the front surface of the target laser，f ; Draw the scattered light output voltage curve V on the surface of the protective plate behind the target laser，b ; Voltage curve V laser，f And the voltage curve V laser，b The data are collected simultaneously and have the same time axis; wherein f represents the front surface of the target material and b represents the surface of the protective plate;

[0019] 2.1.2. Read the voltage curve V laser，f The rising edge time t 上升沿，f and falling edge time t 下降沿，f , the laser emission time t is calculated by the following formula 激光 ; Read the voltage curve V laser，b The rising turning point t 上升转折，b , the target perforation time t is calculated by the following formula 穿孔 ;

[0020] t 激光 =t 下降沿，f -t 上升沿，f

[0021] t 穿孔 =t 上升转折，b -t上升沿，f ;

[0022] 2.1.3. Subtract t from the time axis corresponding to the scattered light on the front surface of the target and the scattered light on the surface of the protective plate behind the target. 上升沿，f , so that the origin of the time axis is adjusted to the laser emission moment and the data acquisition timing signal V 时序 Zero time alignment.

[0023] Furthermore, step 2.2 is specifically as follows:

[0024] 2.2.1. According to the heat map sequence file F 热图 The image frame determines the thermal image frame number n corresponding to the laser emission time t 出光,T , the zero-pre-start time t of the heat map sequence is calculated by the following formula 热图，零前起始 , and according to the zero-preceding start time t of the heat map sequence 热图，零前起始 The heat map sequence file F 热图 The time axis of is aligned with the zero time of the data acquisition timing signal;

[0025] t 热图，零前起始 =-n 出光,T / f 热图

[0026] Where, subscript T represents heat map, f 热图 is the frame rate of the heat map;

[0027] 2.2.2. Calculate the thermal image frame number n at the time of target perforation by the following formula 穿孔,T Heat map frame number n with laser stop moment 停光,T :

[0028]

[0029] In the formula, int() is the rounding function;

[0030] 2.2.3. Set the list of heat map frame numbers {n} that need to be fused and sort the list from small to large, where the i-th heat map frame number satisfies n 穿孔，T ≤n i ≤n 停光，T .

[0031] Furthermore, step 2.3 is specifically as follows:

[0032] 2.3.1. According to the spot data file F 光斑 The image frame determines the laser light emission time t corresponding to the image spot frame number n 出光,S The zero-front start time t of the spot sequence is calculated by the following formula: 光斑，零前起始 , and according to the zero-front starting time t of the spot sequence 光斑，零前起始 The spot data file F 光斑The time axis of is aligned with the zero time of the data acquisition timing signal;

[0033] t 光斑，零前起始 =-n 出光,S / f 光斑

[0034] Where, subscript S represents the light spot, f 光斑 F is the spot data file 光斑 The frame rate;

[0035] 2.3.2. Calculate the spot frame number n at the laser stop moment by the following formula 停光,S :

[0036] n 停光,S = int((t 激光 -t 光斑，零前起始 )×f 光斑 ).

[0037] Further, step 3] is specifically as follows:

[0038] 3.0. Take out the heat map frame number n in the order of the heat map frame number list {n} 当前,T , according to the thermal image frame number n at the time of target perforation 穿孔,T , the thermal image frame number n when the laser stops 停光,T and heatmap sequence file F 热图 The frame rate f 热图 , calculate the heat map frame number n 当前,T The corresponding current time t 当前 ;

[0039] 3.1. Read the heat map frame number n obtained in step 2.4 当前,T Corresponding heat map data, get the temperature field matrix At the same time, the temperature field matrix is ​​converted Convert to grayscale image and calculate the coordinates X of the point inside the ablation hole 孔 , contour coordinate X 孔轮廓 ; Read the temperature field matrix The temperature at the corresponding coordinates is the temperature T at the contour 孔轮廓 , calculate the centroid coordinate X of the ablation hole on the grayscale image C，孔 and the size d and area S of the ablated hole;

[0040] 3.2. According to the spot frame number n at the time of laser emission 出光,S , the spot frame number n when the laser stops 停光,S and spot data file F 光斑 The frame rate f 光斑 , calculate the current time t 当前 Corresponding spot data file F 光斑 Middle spot frame number n 当前,Sand the adjacent frame number interval [n t1,S ,n t2,S );

[0041] 3.3、Read the adjacent frame number interval [n t1,S ,n t2,S ) corresponding to the image frame data And interpolate it to get the current time t 当前 The spot matrix I t，当前 , according to the spot matrix I t，当前 Calculate the center of mass coordinate X of the light spot C and the peak coordinate X P .

[0042] Furthermore, in step 2.4, the heat map frame number n 当前,T The corresponding current time t 当前 for:

[0043]

[0044] In step 3.1, the centroid coordinate X of the ablation hole on the grayscale image C，孔 And the size d and area S of the ablation hole are:

[0045]

[0046] S=NΔxΔy

[0047]

[0048] Where Δx is the physical distance between adjacent points in the x direction, Δy is the physical distance between adjacent points in the y direction, and X 孔,i is the centroid coordinate of the ablation hole at point i, and N is the number of data points included in the ablation hole;

[0049] In step 3.2, the spot frame number n 当前,S and the adjacent frame number interval [n t1,S ,n t2,S ) are:

[0050] n 当前，S =t 当前 f 光斑 +n 出光，S

[0051] n 出光，S ≤n t1，S ≤n 当前，S <n t2，S ≤n 停光，S

[0052] n t2，S =n t1，S +1;

[0053] In step 3.3, the spot matrix I t，当前 , centroid coordinate X C and the peak coordinate X P They are:

[0054]

[0055] In the formula, subscript C represents the centroid, P represents the peak value, and X i is the coordinate of point i in the spot matrix.

[0056] Further, step 4] is specifically as follows:

[0057] 4.1. Set the current time t 当前 Corresponding spot matrix I t，当前 The centroid coordinate X C Or peak coordinate X P and the centroid coordinate X of the ablation hole C,孔 Overlap to get the spot matrix I t，当前 and the temperature field matrix The row and column correspondence of ;

[0058] 4.2. According to the row-column correspondence and the contour coordinate X of the ablation hole 孔轮廓 Get the corresponding contour of the ablation hole in the light spot, and calculate the light intensity I on the corresponding contour of the ablation hole in the light spot 孔轮廓 ;

[0059] 4.3. According to the temperature T at the contour 孔轮廓 Get the light intensity I on the contour 孔轮廓 Corresponding relationship with temperature;

[0060] 4.4. The light intensity I on the contour is calculated by the following formula: 孔轮廓 With the line ablation rate v ablate The corresponding relationship;

[0061] Line ablation rate v ablate for:

[0062]

[0063] Where h is the initial thickness of the target material at the contour of the ablation hole.

[0064] Furthermore, the heat map frame number n 出光,T It is obtained by using a neural network method, or by the rising edge of the highest temperature curve of the spot irradiation area of ​​the heat map, or by image grayscale conversion, image binarization, and image contour extraction;

[0065] Heatmap frame number n 穿孔,TIt is obtained by using a neural network method, or by combining the falling edge of the average temperature curve of the thermal map ablation hole area with the thermal map image, or by image grayscale conversion, image binarization, and image contour extraction.

[0066] Furthermore, the thermal map sequence file F of the target front surface temperature field 热图 To use a thermal imager to monitor the temperature field on the front surface of the target, output the thermal image sequence file F 热图 ;

[0067] The spot data file F of the target material front surface at the equal optical path 光斑 In order to use the spot measurement equipment to monitor the equal optical path of the front surface of the target, output the online monitoring spot data file F 光斑 .

[0068] Beneficial effects of the present invention:

[0069] 1. The present invention provides a multi-source heterogeneous data fusion and analysis method for target material laser ablation effect experiments, which can fuse and analyze the spot distribution, temperature field, and ablation hole data on the front surface of the target material, fully mine the information in the target material laser ablation effect experimental data, obtain the corresponding relationship between light intensity and temperature, light intensity and linear ablation rate on the ablation hole contour, and the size and area statistical information of the ablation hole evolving over time, and obtain a large number of experimental laws through a single target material laser ablation effect experiment.

[0070] 2. The present invention provides a method for fusion and analysis of multi-source heterogeneous data in target material laser ablation effect experiments, which is convenient for batch processing of multi-source heterogeneous experimental data obtained on the same target material laser ablation effect experimental platform.

[0071] 3. The multi-source heterogeneous data fusion analysis method for target laser ablation effect experiment of the present invention reduces the experimental time, economic and labor costs. BRIEF DESCRIPTION OF THE DRAWINGS

[0072] Figure 1 It is a flow chart of an embodiment of a method for fusion and analysis of multi-source heterogeneous data in a target laser ablation effect experiment of the present invention;

[0073] Figure 2 1 is a schematic diagram of the time axis in an embodiment of a method for fusion analysis of multi-source heterogeneous data in a target laser ablation effect experiment of the present invention, wherein 1 is a data acquisition timing signal, 21 is a time axis of scattered light on the front surface of the target, 22 is a time axis of scattered light on the surface of the protective plate behind the target, and 3 is a thermal map sequence file F 热图 The time axis, 4 is the spot data file F 光斑 5 is the start time of detecting scattered light from the front surface of the target and the surface of the protective plate behind the target, and 6 is the zero-pre-start time t when the spot measurement device starts to collect the spot 光斑，零前起始 , 7 is the starting time t before the thermal imager starts to collect thermal images热图，零前起始 , 81 is the laser emission time t, 82 is the target perforation time t 穿孔 ;

[0074] Figure 3 In the embodiment of the target laser ablation effect experiment multi-source heterogeneous data fusion analysis method of the present invention, the voltage curve V of completing the laser zero time alignment is shown in FIG. laser，f With voltage curve V laser，b Schematic diagram;

[0075] FIG4( a ) is a schematic diagram of a heat map when laser is not irradiated in an embodiment of a method for fusion and analysis of multi-source heterogeneous data in a target laser ablation effect experiment of the present invention;

[0076] FIG4( b ) is a heat map schematic diagram of the laser light emission time in an embodiment of a method for fusion and analysis of multi-source heterogeneous data in a target laser ablation effect experiment of the present invention;

[0077] FIG4( c ) is a schematic diagram of a heat map of the target material perforation moment in an embodiment of a method for fusion and analysis of multi-source heterogeneous data in a target material laser ablation effect experiment of the present invention;

[0078] Figure 5 This is a schematic diagram of the fusion analysis results in an embodiment of a multi-source heterogeneous data fusion analysis method for a target laser ablation effect experiment of the present invention, wherein 1 is the x-axis and y-axis after the light spot and the thermal map coordinates are aligned, 2 is the light spot intensity distribution on the x-axis, and 3 is the temperature distribution on the x-axis;

[0079] Figure 6 It is a schematic diagram of the corresponding relationship between the light intensity at the ablation hole contour on the x-axis and the line ablation rate obtained by fusion analysis in an embodiment of a multi-source heterogeneous data fusion analysis method for a target laser ablation effect experiment of the present invention. DETAILED DESCRIPTION

[0080] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0081] like Figure 1 As shown, a method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment includes the following steps:

[0082] 1] Define the data acquisition timing signal V 时序 , taking the laser emission time as zero, adjust the data acquisition timing signal V 时序 The laser emission time t = 0;

[0083] In this embodiment, the data acquisition timing signal V 时序 A unified time axis voltage can be generated by the data acquisition system, and the photoelectric detector, thermal imager, and spot measurement device of the scattered light acquisition system can be triggered by the time axis voltage to perform acquisition.

[0084] 2] If Figure 2 , Figure 3 As shown in FIG4 , the laser emission time, target perforation time, thermal map frame number at the target perforation moment, thermal map frame number at the laser stop moment, spot frame number at the laser emission moment, and spot frame number at the laser stop moment are calculated by using the thermal map sequence file of the target front surface and the target rear protection plate, the temperature field of the target front surface, and the spot data file at the equal optical path of the target front surface; then the time axes corresponding to the scattered light on the target front surface and the scattered light on the target rear protection plate surface, the time axis of the thermal map sequence file of the target front surface temperature field, and the time axis of the spot data file at the equal optical path of the target front surface are aligned with the laser emission moment t=0 respectively; according to the multi-source heterogeneous data fusion analysis requirements of the target laser ablation effect experiment, the thermal map frame number list {n} to be processed is set as follows:

[0085] 2.1. Calculate the laser emission time t using the target front surface and the target rear protection plate 激光 and target perforation time t 穿孔 ; and the time axis corresponding to the scattered light on the front surface of the target and the scattered light on the surface of the protective plate behind the target and the data acquisition timing signal V 时序 Zero-time alignment of

[0086] 2.1.1. Draw the scattered light output voltage curve V on the front surface of the target laser，f ; Draw the scattered light output voltage curve V on the surface of the protective plate behind the target laser，b ; The voltage curve V laser，f And the voltage curve V laser，b The data are collected simultaneously and have the same time axis; wherein f represents the front surface of the target material and b represents the surface of the protective plate;

[0087] 2.1.2. Read the voltage curve V obtained in step 2.1.1 laser，f The rising edge time t 上升沿，f and falling edge time t 下降沿，f , the laser emission time t is calculated by the following formula 激光 ; Read the voltage curve V obtained in step 2.1.1 laser，b The rising turning point t 上升转折，b , the target perforation time t is calculated by the following formula 穿孔 ;

[0088] t 激光 =t 下降沿，f-t 上升沿，f

[0089] t 穿孔 =t 上升转折，b -t 上升沿，f ;

[0090] 2.1.3. Subtract the time axis corresponding to the scattered light on the front surface of the target and the scattered light on the surface of the protective plate behind the target from the time axis obtained in step 2.1.2. 上升沿，f , so that the origin of the time axis is adjusted to the laser emission moment and the data acquisition timing signal V 时序 Zero time alignment.

[0091] 2.2. Processing the thermal map sequence file F of the target front surface temperature field 热图 , get the zero-pre-start time t of the heat map sequence 热图，零前起始 , thermal image frame number n at the time of target perforation 穿孔,T and the thermal image frame number n at the laser stop time 停光,T , the heat map sequence file F 热图 The time axis of is aligned with the zero time of the data acquisition timing signal, and a list of heat map frame numbers to be processed {n} is set; wherein the subscript T represents the heat map;

[0092] 2.2.1. Use a thermal imager to monitor the temperature field on the front surface of the target and output the thermal image sequence file F 热图 ; According to the heat map sequence file F 热图 The image frame determines the thermal image frame number n corresponding to the laser emission time t 出光,T , the zero-pre-start time t of the heat map sequence is calculated by the following formula 热图，零前起始 , and according to the zero-preceding start time t of the heat map sequence 热图，零前起始 The heat map sequence file F 热图 The time axis of is aligned with the zero time of the data acquisition timing signal;

[0093] t 热图，零前起始 =-n 出光,T / f 热图

[0094] Where, subscript T represents heat map, f 热图 is the frame rate of the heat map;

[0095] In other embodiments, the thermal image frame number n corresponding to the laser emission time t 出光,T It can also be obtained through any of the following methods:

[0096] (1) Through the neural network method, a large number of thermal image images at the time of laser non-irradiation and laser emission are annotated and used as training sets to obtain the laser emission prediction model. The prediction model is then used to train the thermal image sequence file F 热图Perform frame-by-frame analysis and take the thermal image frame number corresponding to the first frame of laser light judgment result as n 出光,T ;

[0097] (2) The thermal image frame number corresponding to the laser light emission moment is obtained through the rising edge of the highest temperature curve of the spot irradiation area of ​​the thermal image as n 出光,T ;

[0098] (3) The thermal image is converted into grayscale, binarized, and image hole detected to obtain the thermal image frame number at the time of laser emission as n 出光,T .

[0099] 2.2.2. According to the zero-pre-start time t of the heat map sequence obtained in step 2.2.1 热图，零前起始 The thermal image frame number n at the time of target perforation is calculated by the following formula 穿孔,T Heat map frame number n with laser stop moment 停光,T :

[0100]

[0101] In the formula, int() is the rounding function;

[0102] In other embodiments, the thermal image frame number n corresponding to the laser emission time t 穿孔,T It can also be obtained through any of the following methods:

[0103] (1) A large number of heat maps without perforation and heat maps with perforation are annotated by a neural network method and used as training sets to obtain a perforation judgment model. The heat map sequence file F is trained according to the perforation judgment model. 热图 The frame number of the heat map corresponding to the first heat map of the perforation is taken as n. 穿 Kong, T.;

[0104] (2) The falling edge of the average temperature curve of the ablation hole area of ​​the thermal map is combined with the thermal map image to obtain the thermal map frame number corresponding to the target material perforation time as n 穿孔,T ;

[0105] (3) The thermal image is converted into grayscale, binarized, and image hole detected to obtain the thermal image frame number n at the time of target perforation. 穿孔,T ;

[0106] (4) Determine n by viewing the thermal image frame by frame 穿孔,T .

[0107] 2.2.3. According to the analysis needs, set the list of heat map frame numbers {n} that need to be fused and analyzed, and arrange the list in ascending order, where the i-th heat map frame number satisfies n 穿孔，T ≤n i ≤n 停光，T .

[0108] 2.3. Process the spot data file F at the equal optical path on the front surface of the target 光斑 , get the zero-front starting time t of the spot sequence 光斑，零前起始 , the spot frame number n at the time of laser emission 出光,S and the spot frame number n at the moment when the laser stops 停光,S , the spot data file F 光斑 The time axis of is aligned with the zero time of the data acquisition timing signal; wherein the subscript S represents the light spot;

[0109] 2.3.1. Use the spot measurement equipment to monitor the equal optical path of the front surface of the target and output the online monitoring spot data file F 光斑 ; According to the spot data file F 光斑 The image frame determines the laser light emission time t corresponding to the image spot frame number n 出光,S The zero-front start time t of the spot sequence is calculated by the following formula: 光斑，零前起始 , and according to the zero-front starting time t of the spot sequence 光斑，零前起始 The spot data file F 光斑 The time axis of is aligned with the zero time of the data acquisition timing signal;

[0110] t 光斑，零前起始 =-n 出光,S / f 光斑

[0111] Where, subscript S represents the light spot, f 光斑 F is the spot data file 光斑 Frame rate;

[0112] 2.3.2. According to the zero-pre-start time t of the spot sequence obtained in step 2.3.1 光斑，零前起始 The spot frame number n at the moment when the laser stops is calculated by the following formula 停光,S :

[0113] n 停光,S = int((t 激光 -t 光斑，零前起始 )×f 光斑 );

[0114] 3] Calculate the coordinates of the inner point of the ablation hole, the contour coordinates, the centroid coordinates of the ablation hole, the size and area of ​​the ablation hole, and the spot frame number and the adjacent frame number interval in the spot data file corresponding to the current time [n] according to the thermal map sequence file, the spot data file and the thermal map frame number list {n} t1,S ,n t2,S ); read the adjacent frame number interval [n t1,S ,n t2,S ) corresponding to the image frame data I nt1 ,I nt2, and interpolate it to get the spot matrix at the current moment, and calculate the centroid coordinate X of the spot according to the spot matrix C and the peak coordinate X P ;

[0115] 3.0. Take out the heat map frame number n in the order of the heat map frame number list {n} 当前,T , according to the thermal image frame number n of the target perforation moment obtained in step 2.2.2 穿孔,T , the thermal image frame number n when the laser stops 停光,T and heatmap sequence file F 热图 The frame rate f 热图 , calculate the heat map frame number n 当前,T The corresponding current time t 当前 :

[0116]

[0117] 3.1. Read the heat map frame number n obtained in step 2.4 当前,T Corresponding heat map data, get the temperature field matrix At the same time, the temperature field matrix is ​​converted to The grayscale image is converted to a grayscale image, and the coordinates of the inner point of the ablation hole are calculated by image binarization, image hole detection, and image contour extraction. 孔 , contour coordinate X 孔轮廓 ; Read the temperature field matrix The temperature at the corresponding coordinates is the temperature T at the contour 孔轮廓 , calculate the centroid coordinate X of the ablation hole on the grayscale image C，孔 and the size d and area S of the ablated hole;

[0118]

[0119] S=NΔxΔy

[0120]

[0121] Where Δx is the physical distance between adjacent points in the x direction, Δy is the physical distance between adjacent points in the y direction, and X 孔,i is the centroid coordinate of the ablation hole at point i, and N is the number of data points included in the ablation hole;

[0122] 3.2. The spot frame number n at the laser emission time obtained in step 2.3.1 出光,S , the spot frame number n at the laser stop moment obtained in step 2.3.2 停光,S and spot data file F 光斑 The frame rate f 光斑 , calculate the current time t 当前 Corresponding spot data file F 光斑Middle spot frame number n 当前,S and the adjacent frame number interval [n t1,S ,n t2,S );

[0123] n 当前，S =t 当前 f 光斑 +n 出光，S

[0124] n 出光，S ≤n t1，S ≤n 当前，S <n t2，S ≤n 停光，S

[0125] n t2，S =n t1，S +1

[0126] 3.3、Read the adjacent frame number interval [n t1,S ,n t2,S ) corresponding to the image frame data And interpolate it to get the current time t 当前 The spot matrix I t，当前 , according to the spot matrix I t，当前 Calculate the center of mass coordinate X of the light spot C and the peak coordinate X P :

[0127]

[0128] In the formula, subscript C represents the centroid, P represents the peak value, and X i is the coordinate of point i in the spot matrix.

[0129] 4] If Figure 5 , Figure 6 As shown, the centroid coordinate X C Or peak coordinate X P and the centroid coordinate X of the ablation hole C，孔 Overlap is performed to obtain the row-column correspondence between the spot matrix and the temperature field matrix; the corresponding contour of the ablation hole in the spot is obtained according to the row-column correspondence and the contour coordinates of the ablation hole, and the light intensity on the corresponding contour of the ablation hole in the spot is calculated; then the correspondence between the light intensity and the temperature on the contour is obtained according to the temperature at the contour, and finally the correspondence between the light intensity and the linear ablation rate is calculated, which is specifically:

[0130] 4.1. Set the current time t 当前 Corresponding spot matrix I t，当前 The centroid coordinate X C Or peak coordinate X P and the centroid coordinate X of the ablation hole C，孔 Overlap to get the spot matrix I t，当前and the temperature field matrix The row and column correspondence of ;

[0131] 4.2. According to the row-column correspondence and the contour coordinate X of the ablation hole 孔轮廓 Get the corresponding contour of the ablation hole in the light spot, and calculate the light intensity I on the corresponding contour of the ablation hole in the light spot 孔轮廓 ;

[0132] 4.3. According to the temperature T at the contour 孔轮廓 Get the light intensity I on the contour 孔轮廓 Corresponding relationship with temperature;

[0133] 4.4. The light intensity I on the contour is calculated by the following formula: 孔轮廓 With the line ablation rate v ablate The corresponding relationship: the line ablation rate v ablate for:

[0134]

[0135] Where h is the initial thickness of the target at the ablation hole contour, t 当前 is the time corresponding to the current contour.

[0136] 5] Determine whether the heat map frame number list {n} reaches the end point. If so, output the light intensity I on the contour 孔轮廓 The time series of the relationship between temperature, light intensity and linear ablation rate, as well as the results of the evolution of the size and area of ​​the ablation hole over time, complete the multi-source heterogeneous data fusion analysis of the target laser ablation effect experiment; otherwise, return to step 3.0.

[0137] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions within the technical scope disclosed by the present invention should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims

1. A method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiments. It is characterized in that The following steps are involved: 1] Define the data acquisition timing signal V 时序 , taking the laser emission time as zero, adjust the data acquisition timing signal V 时序 The laser emission time t = 0; 2] Using the thermal map sequence files of the target material front surface and the target material rear protection plate, the target material front surface temperature field, and the spot data files at the equal optical path of the target material front surface, calculate the laser light emission time, the target material perforation time, the thermal map frame number at the target material perforation time, the thermal map frame number at the laser light stop time, the spot frame number at the laser light emission time, and the spot frame number at the laser light stop time; then align the time axes corresponding to the scattered light on the target material front surface and the scattered light on the target material rear protection plate surface, the time axis of the thermal map sequence files of the target material front surface temperature field, and the time axis of the spot data files at the equal optical path of the target material front surface with the laser light emission time t=0, and set the thermal map frame number list {n} that needs to be processed; 3] Calculate the coordinates of the inner point of the ablation hole, the contour coordinates, the centroid coordinates of the ablation hole, the size and area of ​​the ablation hole, and the spot frame number and the adjacent frame number interval in the spot data file corresponding to the current time [n] according to the thermal map sequence file, the spot data file and the thermal map frame number list {n} t1,S ,n t2,S ); read the adjacent frame number interval [n t1,S ,n t2,S ) corresponding to the image frame data I nt1 ,I nt2 , and interpolate it to get the spot matrix at the current moment, and calculate the centroid coordinate X of the spot according to the spot matrix C and the peak coordinate X P ; 4] Set the centroid coordinate X C Or peak coordinate X P Overlap with the centroid coordinates of the ablation hole to obtain the row-column correspondence between the spot matrix and the temperature field matrix; obtain the corresponding contour of the ablation hole in the spot according to the row-column correspondence and the contour coordinates of the ablation hole, and calculate the light intensity on the corresponding contour of the ablation hole in the spot; then obtain the corresponding relationship between the light intensity and the temperature on the contour according to the temperature at the contour, and finally calculate the corresponding relationship between the light intensity and the linear ablation rate; 5] Determine whether the thermal map frame number list {n} has reached the end point. If so, output the time series of the relationship between light intensity and temperature, light intensity and line ablation rate on the contour, and the results of the evolution of the size and area of ​​the ablation hole over time to complete the multi-source heterogeneous data fusion analysis of the target laser ablation effect experiment; otherwise, return to step 3].

2. According to the method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment of claim 1, it is characterized in that step 2] specifically comprises: 2.

1. Calculate the laser emission time t using the target front surface and the target rear protection plate 激光 and target perforation time t 穿孔 , and the time axis corresponding to the scattered light on the front surface of the target and the scattered light on the surface of the protective plate behind the target is connected to the data acquisition timing signal V 时序 Zero-time alignment of 2.

2. Processing the thermal map sequence file F of the target front surface temperature field 热图 , get the zero-pre-start time t of the heat map sequence 热图，零前起始 , thermal image frame number n at the time of target perforation 穿孔,T and the thermal image frame number n at the laser stop time 停光,T , the heat map sequence file F 热图 The time axis of is aligned with the zero time of the data acquisition timing signal, and a list of heat map frame numbers to be processed {n} is set; in, Subscript T indicates heat map; 2.

3. Process the spot data file F at the equal optical path on the front surface of the target 光斑 , get the zero-front starting time t of the spot sequence 光斑，零前起始 , the spot frame number n at the time of laser emission 出光,S and the spot frame number n at the moment when the laser stops 停光,S , the spot data file F 光斑 The time axis of is aligned with the zero time of the data acquisition timing signal; wherein the subscript S represents the light spot.

3. According to the method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment of claim 2, It is characterized in that Step 2.1 is as follows: 2.1.

1. Draw the scattered light output voltage curve V on the front surface of the target laser，f ; Draw the scattered light output voltage curve V on the surface of the protective plate behind the target laser，b ; The voltage curve V laser，f And the voltage curve V laser，b Simultaneous acquisition, which has the same time axis; in, f represents the front surface of the target, and b represents the surface of the protective plate; 2.1.

2. Read the voltage curve V laser，f The rising edge time t 上升沿，f and falling edge time t 下降沿，f , the laser emission time t is calculated by the following formula 激光 ; Read the voltage curve V laser，b The rising turning point t 上升转折，b , the target perforation time t is calculated by the following formula 穿孔 ; t 激光 =t 下降沿，f -t 上升沿，f t 穿孔 =t 上升转折，b -t 上升沿，f ; 2.1.

3. Subtract t from the time axis corresponding to the scattered light on the front surface of the target and the scattered light on the surface of the protective plate behind the target. 上升沿，f , so that the origin of the time axis is adjusted to the laser emission moment and the data acquisition timing signal V 时序 Zero time alignment.

4. According to claim 3, a method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment, It is characterized in that Step 2.2 is as follows: 2.2.

1. According to the heat map sequence file F 热图 The image frame determines the thermal image frame number n corresponding to the laser emission time t 出光,T , the zero-pre-start time t of the heat map sequence is calculated by the following formula 热图，零前起始 , and according to the zero-preceding start time t of the heat map sequence 热图，零前起始 The heat map sequence file F 热图 The time axis of is aligned with the zero time of the data acquisition timing signal; t 热图，零前起始 =-n 出光,T / f 热图 Where, subscript T represents heat map, f 热图 is the frame rate of the heat map; 2.2.

2. Calculate the thermal image frame number n at the time of target perforation by the following formula 穿孔,T Heat map frame number n with laser stop moment 停光,T : In the formula, int() is the rounding function; 2.2.

3. Set the list of heat map frame numbers {n} that need to be fused and sort the list from small to large, where the i-th heat map frame number satisfies n 穿孔，T ≤n i ≤n 停光，T .

5. The method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment according to claim 4, characterized in that step 2.3 specifically comprises: 2.3.

1. According to the spot data file F 光斑 The image frame determines the laser light emission time t corresponding to the image spot frame number n 出光,S The zero-front start time t of the spot sequence is calculated by the following formula: 光斑，零前起始 , and according to the zero-front starting time t of the spot sequence 光斑，零前起始 The spot data file F 光斑 The time axis of is aligned with the zero time of the data acquisition timing signal; t 光斑，零前起始 =-n 出光,S / f 光斑 Where, subscript S represents the light spot, f 光斑 F is the spot data file 光斑 The frame rate; 2.3.

2. Calculate the spot frame number n at the laser stop moment by the following formula 停光,S : n 停光,S =int((t 激光 -t 光斑，零前起始 )×f 光斑 )。 6. The method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment according to claim 5, characterized in that step 3] specifically comprises: 3.

0. Take out the heat map frame number n in the order of the heat map frame number list {n} 当前,T , according to the thermal image frame number n at the time of target perforation 穿孔,T , the thermal image frame number n when the laser stops 停光,T and heatmap sequence file F 热图 The frame rate f 热图 , calculate the heat map frame number n 当前,T The corresponding current time t 当前 ; 3.

1. Read the heat map frame number n obtained in step 2.4 当前,T Corresponding heat map data, get the temperature field matrix At the same time, the temperature field matrix is ​​converted to Convert to grayscale image and calculate the coordinates X of the point inside the ablation hole 孔 , contour coordinate X 孔轮廓 ; Read the temperature field matrix The temperature at the corresponding coordinates is the temperature T at the contour 孔轮廓 , calculate the centroid coordinate X of the ablation hole on the grayscale image C，孔 and the size d and area S of the ablated hole; 3.

2. According to the spot frame number n at the time of laser emission 出光,S , the spot frame number n when the laser stops 停光,S and spot data file F 光斑 The frame rate f 光斑 , calculate the current time t 当前 Corresponding spot data file F 光斑 Middle spot frame number n 当前,S and the adjacent frame number interval [n t1,S ,n t2,S ); 3.3、Read the adjacent frame number interval [n t1,S ,n t2,S ) corresponding to the image frame data And interpolate it to get the current time t 当前 The spot matrix I t，当前 , according to the spot matrix I t，当前 Calculate the center of mass coordinate X of the light spot C and the peak coordinate X P .

7. The method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment according to claim 6, Features: In step 2.4, the heat map frame number n 当前,T The corresponding current time t 当前 for: In step 3.1, the centroid coordinate X of the ablation hole on the grayscale image C，孔 And the size d and area S of the ablation hole are: S=NΔxΔy Where Δx is the physical distance between adjacent points in the x direction, Δy is the physical distance between adjacent points in the y direction, and X 孔,i is the centroid coordinate of the ablation hole at point i, and N is the number of data points included in the ablation hole; In step 3.2, the light spot frame number n 当前,S and the adjacent frame number interval [n t1,S ,n t2,S ) are: n 当前，S =t 当前 f 光斑 +n 出光，S n 出光，S ≤n t1，S ≤n 当前，S <n t2，S ≤n 停光，S n t2，S =n t1，S +1; In step 3.3, the light spot matrix I t，当前 , centroid coordinate X C and the peak coordinate X P They are: In the formula, subscript C represents the centroid, P represents the peak value, and X i is the coordinate of point i in the spot matrix.

8. The method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment according to claim 7, Its characteristics are, Step 4] Specifically: 4.

1. Set the current time t 当前 Corresponding spot matrix I t，当前 The centroid coordinate X C Or peak coordinate X P and the centroid coordinate X of the ablation hole C,孔 Overlap to get the spot matrix I t，当前 and the temperature field matrix The row and column correspondence of ; 4.

2. According to the row-column correspondence and the contour coordinate X of the ablation hole 孔轮廓 Get the corresponding contour of the ablation hole in the light spot, and calculate the light intensity I on the corresponding contour of the ablation hole in the light spot 孔轮廓 ; 4.

3. According to the temperature T at the contour 孔轮廓 Get the light intensity I on the contour 孔轮廓 Corresponding relationship with temperature; 4.

4. The light intensity I on the contour is calculated by the following formula: 孔轮廓 With the line ablation rate v ablate The corresponding relationship; The line ablation rate v ablate for: Where h is the initial thickness of the target material at the contour of the ablation hole.

9. The method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment according to claim 8, Its characteristics are: The heat map frame number n 出光,T It is obtained by using a neural network method, or by the rising edge of the highest temperature curve of the spot irradiation area of ​​the heat map, or by image grayscale conversion, image binarization, and image contour extraction; The heat map frame number n 穿孔,T It is obtained by using a neural network method, or by combining the falling edge of the average temperature curve of the thermal map ablation hole area with the thermal map image, or by image grayscale conversion, image binarization, and image contour extraction.

10. The method for fusion and analysis of multi-source heterogeneous data in target laser ablation effect experiment according to claim 9, Its characteristics are: The thermal map sequence file F of the target front surface temperature field 热图 To use a thermal imager to monitor the temperature field on the front surface of the target, output the thermal image sequence file F 热图 ; The target material front surface equal optical path spot data file F 光斑 In order to use the spot measurement equipment to monitor the equal optical path of the front surface of the target, output the online monitoring spot data file F 光斑 .

Images