WO2023221282A1 - Model correction method, spectroscopy device, computer device, and storage medium - Google Patents

Model correction method, spectroscopy device, computer device, and storage medium Download PDF

Info

Publication number
WO2023221282A1
WO2023221282A1 PCT/CN2022/106217 CN2022106217W WO2023221282A1 WO 2023221282 A1 WO2023221282 A1 WO 2023221282A1 CN 2022106217 W CN2022106217 W CN 2022106217W WO 2023221282 A1 WO2023221282 A1 WO 2023221282A1
Authority
WO
WIPO (PCT)
Prior art keywords
spectral data
model
predicted
correction
original
Prior art date
Application number
PCT/CN2022/106217
Other languages
French (fr)
Chinese (zh)
Inventor
潘从元
张兵
贾军伟
Original Assignee
合肥金星智控科技股份有限公司
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 合肥金星智控科技股份有限公司 filed Critical 合肥金星智控科技股份有限公司
Publication of WO2023221282A1 publication Critical patent/WO2023221282A1/en

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/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/71Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Definitions

  • the present application relates to the technical field of laser-induced breakdown spectroscopy analysis, and in particular to a model correction method, spectral equipment, computer equipment and storage media.
  • Embodiments of the present application provide a model correction method, spectral equipment, computer equipment and storage media.
  • Predictive model corrections for multiple focus distances are based on errors between the predicted spectral data and the original spectral data.
  • the model correction method of the embodiment of the present application can realize the migration correction of the prediction model of laser-induced breakdown spectroscopy equipment under different focusing distances, and improve the accuracy and stability of the prediction model.
  • Laser-induced breakdown spectroscopy equipment can quickly adapt to changes in on-site test conditions, solving the problem of prediction model failure caused by changes in focus distance due to changes in the position of detected materials at the application site. There is no need to recalibrate, saving manpower, material resources and time.
  • predictive model correction for multiple focus distances based on an error between the predicted spectral data and the original spectral data includes:
  • the prediction model of multiple focus distances is corrected using the determined value of the correction function.
  • establishing a correction function between original spectral data and predicted spectral data includes:
  • the correction function is established based on the original concentration and the predicted concentration.
  • the correction function adopts the following relationship:
  • c is the predicted concentration of the element at the focusing distance l of the prediction model
  • C is the original concentration of the element
  • Z(l) is the corresponding correction function at the focusing distance l.
  • determining a value of the correction function based on the relationship includes:
  • the minimum value is used to determine the value of the correction function.
  • the loss function adopts the following relationship:
  • J is the loss function
  • i is the sample number
  • j is the focusing distance number
  • M is the number of samples
  • N is the number of focusing distances
  • C i is the original concentration of the element of the i-th sample
  • Z(l j ) is the j-th
  • c j i is the predicted concentration of the element at the jth focusing distance l j of the i-th sample of the prediction model.
  • the model correction method includes:
  • the corrected prediction model is used to predict the elements in the sample at the corresponding focusing distance.
  • the spectroscopic equipment in the embodiment of the present application includes a sample stage, a laser and a spectrometer.
  • the sample stage is used to carry the sample
  • the laser is used to emit laser light to the sample
  • the spectrometer is used to receive the laser light reflected by the sample.
  • the spectrometer includes a processor, and the processor is used to implement the model correction method in any of the above embodiments.
  • the spectrometer of the spectroscopic equipment in the embodiment of the present application has a processor that can implement the model correction method.
  • the model correction method can be used to collect effective spectra without recalibration, saving manpower, material resources and time.
  • the computer device in the embodiment of the present application includes a memory and a processor.
  • the memory stores a computer program.
  • the processor executes the computer program, it implements the model correction method in any of the above embodiments.
  • the non-volatile computer-readable storage medium of computer-executable instructions when the computer-executable instructions are executed by one or more processors, cause the processor to execute any one of the above-mentioned embodiments.
  • the model calibration method when the computer-executable instructions are executed by one or more processors, cause the processor to execute any one of the above-mentioned embodiments.
  • Figure 1 is a schematic flow chart of the model correction method in the embodiment of the present application.
  • Figure 2 is a schematic structural diagram of the spectroscopic equipment in the embodiment of the present application.
  • Figure 3 is a schematic diagram of the model correction method in the embodiment of the present application to obtain spectral data at multiple focusing distances;
  • Figure 4 is a diagram showing the results of the prediction model established at the focus distance of 152.5cm by the model correction method in the embodiment of the present application, respectively predicting different focus distance data before correction;
  • Figure 5 is a graph showing the relationship between the focusing distance and error of one of the samples in Figure 4.
  • Figure 6 is a schematic flow chart of the model correction method in the embodiment of the present application.
  • Figure 7 is a schematic flow chart of the model correction method in the embodiment of the present application.
  • Figure 8 is a schematic flow chart of the model correction method in the embodiment of the present application.
  • Figure 9 is a diagram showing the results of predicting different focus distance data after correcting the prediction model established by the model correction method in the embodiment of the present application at a focus distance of 152.5 cm.
  • Sample stage 10 Sample stage 10, laser 20, spectrometer 30, beam expander 40, reflector 50, focusing lens 60, collection lens 70, optical fiber 80, processor 31.
  • first and second are used for descriptive purposes only and cannot be understood as indicating or implying relative importance or implicitly indicating the quantity of indicated technical features.
  • features defined as “first” and “second” may explicitly or implicitly include one or more of the described features.
  • “plurality” means two or more than two, unless otherwise explicitly and specifically limited.
  • the term “above” or “below” a first feature on a second feature may include direct contact between the first and second features, or may also include the first and second features. Not in direct contact but through additional characteristic contact between them.
  • the terms “above”, “above” and “above” a first feature on a second feature include the first feature being directly above and diagonally above the second feature, or simply mean that the first feature is higher in level than the second feature.
  • “Below”, “under” and “under” the first feature is the second feature includes the first feature being directly below and diagonally below the second feature, or simply means that the first feature is less horizontally than the second feature.
  • the model correction method in the embodiment of this application includes:
  • Step S10 Based on laser-induced breakdown spectroscopy technology, obtain original spectral data at multiple focusing distances;
  • Step S20 Establish a prediction model using the original spectral data at one of the focusing distances
  • Step S30 Use the prediction model to predict the predicted spectral data of the remaining focusing distances
  • Step S40 Calibrate the prediction model of multiple focus distances based on the error between the predicted spectral data and the original spectral data.
  • the spectroscopic equipment 100 in the embodiment of the present application includes a sample stage 10, a laser 20 and a spectrometer 30.
  • the sample stage 10 is used to carry the sample
  • the laser 20 is used to emit laser light to the sample
  • the spectrometer 30 is used to receive the laser light reflected by the sample.
  • the spectrometer 30 includes a processor 31, which is used to implement the above model correction method.
  • the processor 31 is used to obtain original spectral data of multiple focusing distances based on laser-induced breakdown spectroscopy technology; and is used to establish a prediction model using the original spectral data of one focusing distance; and is used to use the prediction model to predict Predicted spectral data for the remaining focus distances; and prediction model correction for multiple focus distances based on errors between the predicted spectral data and the original spectral data.
  • the spectrometer 30 of the spectroscopic device 100 in the embodiment of the present application has a processor 31 that can implement the model correction method.
  • the model correction method can be used to collect effective spectra without recalibration, saving manpower and Material resources and time.
  • the spectroscopy device 100 can be a device for spectrum analysis such as a laser-induced breakdown spectroscopy device.
  • the laser 20 of the spectrometry device 100 can emit laser light.
  • the laser can be a high-energy pulse laser.
  • the laser light emitted by the laser 20 can pass through a beam expander. 40 and then hit the reflective mirror 50, which can reflect the laser to the sample stage 10.
  • a focusing lens 60 can be provided between the sample stage 10 and the reflector 50.
  • the laser light reflected to the sample stage 10 can be scattered after contacting the sample on the sample stage 10.
  • the laser light scattered by the sample stage 10 can be The collection lens 70 collects, and then the collected laser light can be transmitted to the spectrometer 30 through the optical fiber 80 and analyzed by the processor 31, so that the spectroscopic device 100 can collect the original spectral data of the sample.
  • the prediction model migration correction of the laser-induced breakdown spectroscopy device under different focusing distances can be realized, thereby improving the accuracy and stability of the prediction model.
  • Laser-induced breakdown spectroscopy equipment can quickly adapt to changes in on-site test conditions, solving the problem of prediction model failure caused by changes in focus distance due to changes in the position of detected materials at the application site. There is no need to recalibrate, saving manpower, material resources and time.
  • the laser-induced breakdown spectroscopy (LIBS) technology based on the model calibration method is a composition analysis technology based on the one-to-one correspondence between the wavelengths of atomic spectra and ion spectra and specific elements, and the spectral signal intensity is related to the corresponding element.
  • the content also has a certain quantitative relationship.
  • Plasma is formed by focusing the surface of the sample with a high-energy pulse laser.
  • a spectrometer is used to record the spectral information emitted by the plasma. By analyzing the plasma spectrum, the elements in the sample are analyzed based on the position of the characteristic wavelength and the spectral intensity. Qualitative and quantitative analysis.
  • the model calibration method takes step S10, which can use LIBS technology equipment to acquire the original spectral data of the sample at multiple focusing distances.
  • the sample may be a metal sample, such as steel, alloy, etc., the number of samples may be multiple, and multiple samples may have multiple samples of the same element type.
  • Multiple focusing distances may be the distance at which the optical system of the LIBS technology device focuses on the sample. For example, as shown in Figure 3, equipment using LIBS technology acquires raw spectral data from the sample at two focusing distances.
  • the focusing distance from the sample to the optical system of the LIBS technology equipment is l 1
  • the original spectral data is obtained from the sample at position a with a focusing distance of l 1
  • the sample can be placed at position b
  • the original spectral data is obtained from the sample at position b with a focusing distance of l 2 .
  • step S20 the model correction method takes step S20 to establish a prediction model using the original spectral data of one of the multiple focusing distances obtained in step S10.
  • the prediction model established in step S20 can predict the spectral data of the remaining focusing distances among the plurality of focusing distances based on the original spectral data of one of the focusing distances.
  • the model correction method takes step S40 to correct the prediction model established in step S20 based on the error between the predicted spectral data of multiple focusing distances in step S30 and the original spectral data of multiple focusing distances in step S10.
  • the sample is carbon steel, the sample quantity is 10 pieces, and the sample number is 0 to 9. There are 11 focusing distances between 152.5cm and 376.5cm.
  • LIBS technology equipment can be used to obtain raw spectral data of 11 carbon steel samples at 11 focusing distances between 152.5cm and 376.5cm.
  • step S20 the original spectral data of 10 carbon steel samples collected at a focusing distance of 152.5 cm can be used to establish a prediction model that can predict the spectral data.
  • step S30 the predicted spectral data of 10 pieces of carbon steel with 11 focusing distances between 152.5cm and 376.5cm can be predicted according to the prediction model, and further the quantitative analysis data of the concentration of manganese (Mn) element in the carbon steel can be obtained.
  • the results are shown in Figure 4.
  • step S40 calculate the error between the predicted value and the original value of the 11 focusing distances of each piece of carbon steel at the focusing distance of 152.5cm to 376.5cm.
  • one piece of carbon steel is at the focusing distance of 152.5cm.
  • the prediction model established on the original spectral data of 10 carbon steel samples collected at a focusing distance of 152.5cm can be corrected based on the calculated error.
  • the prediction model correction (step S40) for multiple focusing distances includes:
  • Step S41 Establish a correction function between the original spectral data and the predicted spectral data
  • Step S42 Based on the correction function, establish the relationship between the loss function and the correction function;
  • Step S43 Based on the relationship, determine the value of the correction function
  • Step S44 Use the determined values of the correction function to correct the prediction models of multiple focus distances.
  • the processor 31 is configured to establish a correction function of the original spectral data and the predicted spectral data; and to establish a relationship between the loss function and the correction function based on the correction function; and to determine the correction based on the relationship. a numerical value of the function; and a predictive model correction for a plurality of focus distances using the determined numerical value of the correction function.
  • the functional relationship is used to establish the correction function, the loss function, and the relationship between the correction function and the loss function to correct the prediction model and improve the accuracy and stability of the prediction model.
  • the prediction model for multiple focusing distances can be corrected.
  • the model correction method can first take step S41, using the error between the original spectrum data and the predicted spectrum data to establish a correction function of the original spectrum data and the predicted spectrum data, and then take step S42, based on the correction function established in step S41. , establish the relationship between the loss function and the correction function.
  • the loss function can be a fitting function of the focusing distance and error established using the principle function and construct a loss function.
  • the error is the error data between the original spectral data and the predicted spectral data.
  • step S43 can be taken, and the specific value of the correction function can be obtained based on the relationship between the loss function and the correction function established in step S42.
  • the specific value of the correction function can be obtained by solving the functional relationship between the loss function and the correction function.
  • step S44 is taken, and the specific values of the correction function obtained in step S43 are brought into corresponding prediction models of multiple focus distances to correct the prediction model.
  • establishing a correction function between original spectral data and predicted spectral data includes:
  • Step S411 Determine the original concentration of the element based on the original spectral data
  • Step S412 Determine the predicted concentration of the element based on the predicted spectrum data
  • Step S413 Establish a correction function based on the original concentration and the predicted concentration.
  • the processor 31 is configured to determine an original concentration of the element based on the original spectral data; and to determine a predicted concentration of the element based on the predicted spectral data; and to establish a correction function based on the original concentration and the predicted concentration.
  • the correction function established based on the original spectral data and predicted spectral data is converted into a correction function established based on the original concentration and predicted concentration, thereby improving the accuracy and stability of the prediction model.
  • the model calibration method may first take step S411 to determine the element in the sample detected at the corresponding focus distance based on the relationship between the original spectral data obtained in step S10 and a certain element of the sample detected at the corresponding focus distance. The original concentration of an element.
  • Step S412 may then be taken to determine the predicted concentration of the element in the sample detected at the corresponding focus distance based on the relationship between the predicted spectrum data predicted in step S30 and a certain element in the sample detected at the corresponding focus distance.
  • step S413 may be taken to establish a correction function corresponding to the focus distance based on the original concentration in step S411 and the predicted concentration in step S412.
  • the correction function takes the following relationship:
  • c is the predicted concentration of the element at the focusing distance l of the prediction model
  • C is the original concentration of the element
  • Z(l) is the corresponding correction function at the focusing distance l.
  • the predicted concentration of an element is based on the relationship between the original concentration of the element and the correction function.
  • the correction function can predict the predicted concentration of elements at other focusing distances through the original concentration of the element, realizing the migration correction of the prediction model at different focusing distances. .
  • the relational expression of the correction function determines the relationship between the source, spectral data and element concentration.
  • LIBS technology is based on the one-to-one correspondence between the wavelengths of atomic spectra and ion spectra and specific elements, and mainly involves spectral information, high-energy pulse lasers , plasma, characteristic wavelength and other data.
  • PD is the laser power density
  • E is the pulse energy
  • w is the pulse width
  • d is the focused spot diameter
  • w is the pulse width
  • d is the focused spot diameter
  • i represents the i energy level
  • j represents the j energy level
  • a ij is the transition probability
  • g i is the high energy level degeneracy
  • ⁇ ij is the radiation wavelength
  • U is the current temperature.
  • the matching function corresponding to the ion is the following, E i is the high energy level energy, K B is the Boltzmann constant, where F is the scaling factor related to system parameters, plasma temperature, element characteristics, etc., C is the original concentration of the element, I ij is the element characteristic;
  • I is the characteristic of the element
  • C is the original concentration of the element
  • c is the predicted concentration of the element at the focusing distance l of the prediction model
  • C is the original concentration of the element
  • Z(l) is the corresponding correction function at the focusing distance l.
  • determining the value of the correction function includes:
  • Step S431 Use the least squares method to determine the minimum value of the loss function
  • Step S432 Use the minimum value to determine the value of the correction function.
  • the processor 31 is configured to determine a minimum value of the loss function using a least squares method; and is configured to determine a value of the correction function using the minimum value.
  • the least squares method can be used to solve the minimum value of the loss function, thereby determining the value of the correction function.
  • the model correction method can first take step S431 based on the loss function constructed in step S42, bring the loss function into the data and use the least squares method to obtain the minimum value of the loss function, and then take the step S431.
  • Step S432 Since there is a relationship between the loss function and the correction function, the specific value of the correction function can be further obtained by solving the minimum value of the loss function through the least squares method.
  • the loss function adopts the following relationship:
  • J is the loss function
  • i is the sample number
  • j is the focusing distance number
  • M is the number of samples
  • N is the number of focusing distances
  • C i is the original concentration of the element of the i-th sample
  • Z(l j ) is the j-th
  • c j i is the predicted concentration of the element at the jth focusing distance l j of the i-th sample of the prediction model.
  • the loss function relationship includes the original concentration and predicted element concentration of multiple groups of samples at multiple focusing distances, so that the loss function relationship can achieve data fitting.
  • the relational expression of the loss function can be a functional relational expression constructed from the square of the difference between the predicted concentration of the element obtained according to the correlation correction function of the prediction model and the original concentration of the element.
  • the loss function can be solved by using the least squares method. By making the loss function reach the minimum value, the loss function can be solved.
  • the loss The function can take formula, and then substitute the original concentration of the element, the predicted concentration of the element and other relevant data of the sample corresponding to the sample number in the loss function at the corresponding focusing distance number into the loss function, and solve the minimum value of the loss function through the least squares method.
  • the specific values of the corresponding correction functions are obtained, and then the prediction models of the multiple focus distances can be corrected through the determined values of the correction functions of the multiple focus distances.
  • the model correction method includes:
  • Step S50 Use the corrected prediction model to predict the elements in the sample at the corresponding focusing distance.
  • the processor 31 is configured to use the corrected prediction model to predict elements in the sample at the corresponding focusing distance.
  • the prediction model corrected by the model calibration method is more accurate in predicting the spectral data and element concentration of the sample than the prediction model before correction.
  • the model correction method can use step S50 to correct the prediction model.
  • the prediction model corresponding to the focus distance established in step S20 is corrected based on the specific value of the correction function obtained in step S40.
  • the corrected prediction model can be By predicting the elements in the sample at multiple focusing distances, more accurate element measurement data of the sample can be obtained.
  • the sample is carbon steel, the sample quantity is 10 pieces, and the sample number is 0 to 9. There are 11 focusing distances between 152.5cm and 376.5cm.
  • step S30 the prediction model predicts the predicted spectrum data of 10 pieces of carbon steel at 11 focusing distances between 152.5cm and 376.5cm, and obtains the quantitative analysis data of the concentration of manganese (Mn) element in the carbon steel.
  • the result is: Figure 4.
  • step S50 the elements in the sample at the corresponding focusing distance are predicted according to the corrected prediction model, and 10 pieces of carbon steel at 11 focusing distances between 152.5cm and 376.5cm are predicted according to the corrected prediction model.
  • the predicted spectral data were corrected to obtain the corrected quantitative analysis data of the concentration of manganese (Mn) element in carbon steel. The results are shown in Figure 9.
  • the computer device in the embodiment of the present application includes a memory and a processor.
  • the memory stores a computer program.
  • the processor executes the computer program, it implements the model correction method in any of the above embodiments.
  • the non-volatile computer-readable storage medium of computer-executable instructions in the embodiment of the present application when the computer-executable instructions are executed by one or more processors, causes the processor to execute the model correction method in any of the above-mentioned embodiments.
  • the computer device in the embodiment of the present application can be a calculator, a programmable controller, a desktop computer, a laptop computer, a tablet computer, a server, and other devices.
  • the computer device can include a processor, a memory, and a computer connected through a system bus. Communication interfaces, displays and input devices.
  • the processor of the computer device may be a central processing unit (Central Processing Unit, CPU).
  • the processor can also be other general-purpose processors, Digital Signal Processor (DSP), Application Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA) or other Chips such as programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or combinations of these types of chips.
  • DSP Digital Signal Processor
  • ASIC Application Specific Integrated Circuit
  • FPGA Field-Programmable Gate Array
  • Chips such as programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or combinations of these types of chips.
  • the computer program can be stored in the memory.
  • the memory as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs, non-transitory computer executable programs and modules, as described in the method in the above method embodiment.
  • the processor executes various functional applications and data processing of the processor by running non-transient software programs, instructions and modules stored in the memory, that is, implementing the method in the above method embodiment.
  • a "computer-readable medium” may be any device that can contain, store, communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
  • Non-exhaustive list of computer readable media include the following: electrical connections with one or more wires (electronic device), portable computer disk cartridges (magnetic device), random access memory (RAM), Read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and portable compact disc read-only memory (CDROM).
  • the computer-readable medium may even be paper or other suitable medium on which the program may be printed, as the paper or other medium may be optically scanned, for example, and subsequently edited, interpreted, or otherwise suitable as necessary. process to obtain the program electronically and then store it in computer memory.
  • the processor can be a Central Processing Unit (CPU), other general-purpose processors, Digital Signal Processor (DSP), Application Specific Integrated Circuit (ASIC), or off-the-shelf programmable processors. Gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.
  • a general-purpose processor may be a microprocessor or the processor may be any conventional processor, etc.
  • various parts of the embodiments of the present application can be implemented in hardware, software, firmware, or a combination thereof.
  • various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system.
  • a logic gate circuit with a logic gate circuit for implementing a logic function on a data signal.
  • Discrete logic circuits application specific integrated circuits with suitable combinational logic gates, programmable gate arrays (PGA), field programmable gate arrays (FPGA), etc.
  • the program can be stored in a computer-readable storage medium.
  • the program can be stored in a computer-readable storage medium.
  • each functional unit in various embodiments of the present application can be integrated into a processing module, each unit can exist physically alone, or two or more units can be integrated into one module.
  • the above integrated modules can be implemented in the form of hardware or software function modules. If the integrated module is implemented in the form of a software function module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.
  • the storage media mentioned above can be read-only memory, magnetic disks or optical disks, etc.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

A model correction method, a spectroscopy device, a computer device, and a storage medium. The model correction method comprises: acquiring original spectral data of a plurality of focusing distances on the basis of a laser-induced breakdown spectroscopy technology (S10); establishing a prediction model by using the original spectral data of one of the focusing distances (S20); predicting predicted spectral data of the remaining focusing distances by using the prediction model (S30); and correcting the prediction model for the plurality of focusing distances on the basis of an error between the predicted spectral data and the original spectral data (S40).

Description

模型校正方法、光谱设备、计算机设备和存储介质Model calibration method, spectroscopic equipment, computer equipment and storage media
优先权信息priority information
本申请请求2022年05月18日向中国国家知识产权局提交的、专利申请号为202210538088.7的专利申请的优先权和权益,并且通过参照将其全文并入此处。This application requests the priority and rights of the patent application with patent application number 202210538088.7, which was submitted to the State Intellectual Property Office of China on May 18, 2022, and its full text is incorporated herein by reference.
技术领域Technical field
本申请涉及激光诱导击穿光谱分析技术领域,尤其涉及一种模型校正方法、光谱设备、计算机设备和存储介质。The present application relates to the technical field of laser-induced breakdown spectroscopy analysis, and in particular to a model correction method, spectral equipment, computer equipment and storage media.
背景技术Background technique
当前LIBS设备应用现场受周期性设备检修或堆料更换,待测物料的液面和料面发生较大变化以及设备移动后的测量距离变化,都会导致无法采集到有效光谱,需调整聚焦点,当聚焦点调整后原始模型的精度下降,如果重新定标需耗费大量的人力、物力以及时间。Current LIBS equipment application sites are subject to periodic equipment maintenance or stockpile replacement. Large changes in the liquid and material levels of the materials to be measured, as well as changes in the measurement distance after the equipment is moved, will result in the inability to collect effective spectra, and the focus needs to be adjusted. When the focus point is adjusted, the accuracy of the original model decreases, and recalibration requires a lot of manpower, material resources, and time.
发明内容Contents of the invention
本申请实施方式提供一种模型校正方法、光谱设备、计算机设备和存储介质。Embodiments of the present application provide a model correction method, spectral equipment, computer equipment and storage media.
本申请实施方式的模型校正方法包括:The model correction method in the embodiment of this application includes:
基于激光诱导击穿光谱技术,获取多个聚焦距离的原始光谱数据;Based on laser-induced breakdown spectroscopy technology, original spectral data at multiple focusing distances are obtained;
利用其中一个聚焦距离的原始光谱数据建立预测模型;Build a predictive model using raw spectral data at one of the focusing distances;
利用所述预测模型,预测其余聚焦距离的预测光谱数据;Using the prediction model, predict the predicted spectral data for the remaining focus distances;
基于所述预测光谱数据和所述原始光谱数据之间的误差,对多个聚焦距离的预测模型校正。Predictive model corrections for multiple focus distances are based on errors between the predicted spectral data and the original spectral data.
本申请实施方式的模型校正方法可以实现对激光诱导击穿光谱设备在不同聚焦距离下的预测模型迁移校正,提高预测模型的准确性与稳定性。激光诱导击穿光谱设备能够快速适应现场测试条件变化,解决了应用现场因检测物料位置变化后改变聚焦距离导致的预测模型失效问题,无需重新标定,节省人力、物力和时间。The model correction method of the embodiment of the present application can realize the migration correction of the prediction model of laser-induced breakdown spectroscopy equipment under different focusing distances, and improve the accuracy and stability of the prediction model. Laser-induced breakdown spectroscopy equipment can quickly adapt to changes in on-site test conditions, solving the problem of prediction model failure caused by changes in focus distance due to changes in the position of detected materials at the application site. There is no need to recalibrate, saving manpower, material resources and time.
在某些实施方式中,基于所述预测光谱数据和所述原始光谱数据之间的误差,对多个聚焦距离的预测模型校正,包括:In some embodiments, predictive model correction for multiple focus distances based on an error between the predicted spectral data and the original spectral data includes:
建立原始光谱数据与预测光谱数据的校正函数;Establish a correction function between original spectral data and predicted spectral data;
基于所述校正函数,建立损失函数与所述校正函数之间的关系;Based on the correction function, establish a relationship between the loss function and the correction function;
基于所述关系,确定所述校正函数的数值;Based on the relationship, determining a value of the correction function;
利用确定的所述校正函数的数值对多个聚焦距离的预测模型校正。The prediction model of multiple focus distances is corrected using the determined value of the correction function.
在某些实施方式中,建立原始光谱数据与预测光谱数据的校正函数,包括:In some embodiments, establishing a correction function between original spectral data and predicted spectral data includes:
基于所述原始光谱数据确定元素的原始浓度;determining an original concentration of an element based on the original spectral data;
基于所述预测光谱数据确定元素的预测浓度;determining a predicted concentration of an element based on the predicted spectral data;
基于所述原始浓度和所述预测浓度建立所述校正函数。The correction function is established based on the original concentration and the predicted concentration.
在某些实施方式中,所述校正函数采用以下关系式:In some embodiments, the correction function adopts the following relationship:
C=Z(l)cC=Z(l)c
其中,c为预测模型在聚焦距离l处的元素的预测浓度,C为元素的原始浓度,Z(l)为聚焦距离l处对应的校正函数。Among them, c is the predicted concentration of the element at the focusing distance l of the prediction model, C is the original concentration of the element, and Z(l) is the corresponding correction function at the focusing distance l.
在某些实施方式中,基于所述关系,确定所述校正函数的数值,包括:In some embodiments, determining a value of the correction function based on the relationship includes:
利用最小二乘法确定所述损失函数的最小值;Determine the minimum value of the loss function using the least squares method;
利用所述最小值确定所述校正函数的数值。The minimum value is used to determine the value of the correction function.
在某些实施方式中,所述损失函数采用以下关系式:In some implementations, the loss function adopts the following relationship:
Figure PCTCN2022106217-appb-000001
Figure PCTCN2022106217-appb-000001
其中,J为损失函数,i为样品编号,j为聚焦距离编号,M为样品数量,N为聚焦距离数量,C i为第i个样品的元素的原始浓度,Z(l j)为第j个聚焦距离l j处对应的校正函数,c j i为预测模型的第i个样品在第j个聚焦距离l j处的元素的预测浓度。 Among them, J is the loss function, i is the sample number, j is the focusing distance number, M is the number of samples, N is the number of focusing distances, C i is the original concentration of the element of the i-th sample, Z(l j ) is the j-th The correction function corresponding to the focusing distance l j , c j i is the predicted concentration of the element at the jth focusing distance l j of the i-th sample of the prediction model.
在某些实施方式中,所述模型校正方法包括:In some embodiments, the model correction method includes:
利用校正后的预测模型对相应聚焦距离上的样品中的元素进行预测。The corrected prediction model is used to predict the elements in the sample at the corresponding focusing distance.
本申请实施方式的光谱设备包括样品台、激光器和光谱仪,所述样品台用于承载样品,所述激光器用于向所述样品发出激光,所述光谱仪用于接收所述样品反射的激光,所述光谱仪包括处理器,所述处理器用于实现上述任一项实施方式中所述模型校正方法。The spectroscopic equipment in the embodiment of the present application includes a sample stage, a laser and a spectrometer. The sample stage is used to carry the sample, the laser is used to emit laser light to the sample, and the spectrometer is used to receive the laser light reflected by the sample. The spectrometer includes a processor, and the processor is used to implement the model correction method in any of the above embodiments.
本申请实施方式的光谱设备的光谱仪具有能够实现模型校正方法的处理器,光谱设备在不同测量距离测量时,采用模型校正方法无需重新定标便可采集到有效光谱,节省人力、物力和时间。The spectrometer of the spectroscopic equipment in the embodiment of the present application has a processor that can implement the model correction method. When the spectroscopic equipment is measured at different measurement distances, the model correction method can be used to collect effective spectra without recalibration, saving manpower, material resources and time.
本申请实施方式的计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述任一项实施方式中所述模型校正方法。The computer device in the embodiment of the present application includes a memory and a processor. The memory stores a computer program. When the processor executes the computer program, it implements the model correction method in any of the above embodiments.
本申请实施方式的计算机可执行指令的非易失性计算机可读存储介质,当所述计算机 可执行指令被一个或多个处理器执行时,使得所述处理器执行上述任一项实施方式中所述模型校正方法。The non-volatile computer-readable storage medium of computer-executable instructions according to the embodiment of the present application, when the computer-executable instructions are executed by one or more processors, cause the processor to execute any one of the above-mentioned embodiments. The model calibration method.
本申请的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本申请的实践了解到。Additional aspects and advantages of the application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the application.
附图说明Description of the drawings
本申请的上述和/或附加的方面和优点从结合下面附图对实施方式的描述中将变得明显和容易理解,其中:The above and/or additional aspects and advantages of the present application will become apparent and readily understood from the description of the embodiments in conjunction with the following drawings, in which:
图1是本申请实施方式中的模型校正方法的流程示意图;Figure 1 is a schematic flow chart of the model correction method in the embodiment of the present application;
图2是本申请实施方式中的光谱设备的结构示意图;Figure 2 is a schematic structural diagram of the spectroscopic equipment in the embodiment of the present application;
图3是本申请实施方式中的模型校正方法获取多个聚焦距离的光谱数据的示意图;Figure 3 is a schematic diagram of the model correction method in the embodiment of the present application to obtain spectral data at multiple focusing distances;
图4是本申请实施方式中的模型校正方法在聚焦距离152.5cm下建立的预测模型在校正前分别预测不同聚焦距离数据的结果图;Figure 4 is a diagram showing the results of the prediction model established at the focus distance of 152.5cm by the model correction method in the embodiment of the present application, respectively predicting different focus distance data before correction;
图5是图4的其中一块样品的聚焦距离与误差的关系图;Figure 5 is a graph showing the relationship between the focusing distance and error of one of the samples in Figure 4;
图6是本申请实施方式中的模型校正方法的流程示意图;Figure 6 is a schematic flow chart of the model correction method in the embodiment of the present application;
图7是本申请实施方式中的模型校正方法的流程示意图;Figure 7 is a schematic flow chart of the model correction method in the embodiment of the present application;
图8是本申请实施方式中的模型校正方法的流程示意图;Figure 8 is a schematic flow chart of the model correction method in the embodiment of the present application;
图9是本申请实施方式中的模型校正方法在聚焦距离152.5cm下建立的预测模型校正后分别预测不同聚焦距离数据的结果图。Figure 9 is a diagram showing the results of predicting different focus distance data after correcting the prediction model established by the model correction method in the embodiment of the present application at a focus distance of 152.5 cm.
主要元件符号附图说明:Description of main component symbols and drawings:
光谱仪器100; Spectroscopic Instruments 100;
样品台10、激光器20、光谱仪30、扩束镜40、反射镜50、聚焦透镜60、收集透镜70、光纤80、处理器31。 Sample stage 10, laser 20, spectrometer 30, beam expander 40, reflector 50, focusing lens 60, collection lens 70, optical fiber 80, processor 31.
具体实施方式Detailed ways
下面详细描述本申请的实施方式,所述实施方式的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施方式是示例性的,仅用于解释本申请,而不能理解为对本申请的限制。Embodiments of the present application are described in detail below, examples of which are illustrated in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements with the same or similar functions. The embodiments described below with reference to the drawings are exemplary and are only used to explain the present application and cannot be understood as limiting the present application.
在本申请的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或 元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个所述特征。在本申请的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。In the description of this application, it needs to be understood that the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", " The directions or positional relationships indicated by "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", etc. are based on the directions shown in the accompanying drawings or positional relationship is only for the convenience of describing the present application and simplifying the description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation of the present application. In addition, the terms “first” and “second” are used for descriptive purposes only and cannot be understood as indicating or implying relative importance or implicitly indicating the quantity of indicated technical features. Thus, features defined as “first” and “second” may explicitly or implicitly include one or more of the described features. In the description of this application, "plurality" means two or more than two, unless otherwise explicitly and specifically limited.
在本申请中,除非另有明确的规定和限定,第一特征在第二特征之“上”或之“下”可以包括第一和第二特征直接接触,也可以包括第一和第二特征不是直接接触而是通过它们之间的另外的特征接触。而且,第一特征在第二特征“之上”、“上方”和“上面”包括第一特征在第二特征正上方和斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”包括第一特征在第二特征正下方和斜下方,或仅仅表示第一特征水平高度小于第二特征。In this application, unless otherwise explicitly stated and limited, the term "above" or "below" a first feature on a second feature may include direct contact between the first and second features, or may also include the first and second features. Not in direct contact but through additional characteristic contact between them. Furthermore, the terms "above", "above" and "above" a first feature on a second feature include the first feature being directly above and diagonally above the second feature, or simply mean that the first feature is higher in level than the second feature. “Below”, “under” and “under” the first feature is the second feature includes the first feature being directly below and diagonally below the second feature, or simply means that the first feature is less horizontally than the second feature.
下文的公开提供了许多不同的实施方式或例子用来实现本申请的不同结构。为了简化本申请的公开,下文中对特定例子的部件和设置进行描述。当然,它们仅仅为示例,并且目的不在于限制本申请。此外,本申请可以在不同例子中重复参考数字和/或参考字母,这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施方式和/或设置之间的关系。此外,本申请提供了的各种特定的工艺和材料的例子,但是本领域普通技术人员可以意识到其他工艺的应用和/或其他材料的使用。The following disclosure provides many different embodiments or examples for implementing the various structures of the present application. To simplify the disclosure of the present application, the components and arrangements of specific examples are described below. Of course, they are merely examples and are not intended to limit the application. Furthermore, this application may repeat reference numbers and/or reference letters in different examples, such repetition being for the purposes of simplicity and clarity and does not by itself indicate a relationship between the various embodiments and/or arrangements discussed. In addition, this application provides examples of various specific processes and materials, but one of ordinary skill in the art will recognize the application of other processes and/or the use of other materials.
请参阅图1,本申请实施方式的模型校正方法包括:Please refer to Figure 1. The model correction method in the embodiment of this application includes:
步骤S10:基于激光诱导击穿光谱技术,获取多个聚焦距离的原始光谱数据;Step S10: Based on laser-induced breakdown spectroscopy technology, obtain original spectral data at multiple focusing distances;
步骤S20:利用其中一个聚焦距离的原始光谱数据建立预测模型;Step S20: Establish a prediction model using the original spectral data at one of the focusing distances;
步骤S30:利用预测模型,预测其余聚焦距离的预测光谱数据;Step S30: Use the prediction model to predict the predicted spectral data of the remaining focusing distances;
步骤S40:基于预测光谱数据和原始光谱数据之间的误差,对多个聚焦距离的预测模型校正。Step S40: Calibrate the prediction model of multiple focus distances based on the error between the predicted spectral data and the original spectral data.
请参阅图2,本申请实施方式的光谱设备100包括样品台10、激光器20和光谱仪30,样品台10用于承载样品,激光器20用于向样品发出激光,光谱仪30用于接收样品反射的激光,光谱仪30包括处理器31,处理器31用于实现上述模型校正方法。或者说,处理器31用于基于激光诱导击穿光谱技术,获取多个聚焦距离的原始光谱数据;及用于利用其中一个聚焦距离的原始光谱数据建立预测模型;及用于利用预测模型,预测其余聚焦距离的预测光谱数据;以及用于基于预测光谱数据和原始光谱数据之间的误差,对多个聚焦距离的预测模型校正。Please refer to Figure 2. The spectroscopic equipment 100 in the embodiment of the present application includes a sample stage 10, a laser 20 and a spectrometer 30. The sample stage 10 is used to carry the sample, the laser 20 is used to emit laser light to the sample, and the spectrometer 30 is used to receive the laser light reflected by the sample. , the spectrometer 30 includes a processor 31, which is used to implement the above model correction method. In other words, the processor 31 is used to obtain original spectral data of multiple focusing distances based on laser-induced breakdown spectroscopy technology; and is used to establish a prediction model using the original spectral data of one focusing distance; and is used to use the prediction model to predict Predicted spectral data for the remaining focus distances; and prediction model correction for multiple focus distances based on errors between the predicted spectral data and the original spectral data.
本申请实施方式的光谱设备100的光谱仪30具有能够实现模型校正方法的处理器31,光谱设备100在不同测量距离测量时,采用模型校正方法无需重新定标便可采集到有效光 谱,节省人力、物力和时间。The spectrometer 30 of the spectroscopic device 100 in the embodiment of the present application has a processor 31 that can implement the model correction method. When the spectroscopic device 100 is measured at different measurement distances, the model correction method can be used to collect effective spectra without recalibration, saving manpower and Material resources and time.
具体地,光谱设备100可以是激光诱导击穿光谱设备等用于光谱分析的设备,光谱设备100的激光器20可发射出激光,激光可以是高能脉冲激光,激光器20射出的激光可经过扩束镜40后打在反射镜50上,反射镜50可将激光反射至样品台10。样品台10与反射镜50之间可设有聚焦透镜60,反射至样品台10的激光接触到样品台10上的样品后可发生散射,经过样品台10散射的激光可被样品台10旁的收集透镜70收集,然后收集到的激光可经过光纤80传输至光谱仪30经过处理器31进行分析,由此,光谱设备100可以采集到样品的原始光谱数据。Specifically, the spectroscopy device 100 can be a device for spectrum analysis such as a laser-induced breakdown spectroscopy device. The laser 20 of the spectrometry device 100 can emit laser light. The laser can be a high-energy pulse laser. The laser light emitted by the laser 20 can pass through a beam expander. 40 and then hit the reflective mirror 50, which can reflect the laser to the sample stage 10. A focusing lens 60 can be provided between the sample stage 10 and the reflector 50. The laser light reflected to the sample stage 10 can be scattered after contacting the sample on the sample stage 10. The laser light scattered by the sample stage 10 can be The collection lens 70 collects, and then the collected laser light can be transmitted to the spectrometer 30 through the optical fiber 80 and analyzed by the processor 31, so that the spectroscopic device 100 can collect the original spectral data of the sample.
本申请实施方式的模型校正方法和光谱设备100中,可以实现对激光诱导击穿光谱设备在不同聚焦距离下的预测模型迁移校正,提高预测模型的准确性与稳定性。激光诱导击穿光谱设备能够快速适应现场测试条件变化,解决了应用现场因检测物料位置变化后改变聚焦距离导致的预测模型失效问题,无需重新标定,节省人力、物力和时间。In the model correction method and the spectroscopic device 100 of the embodiment of the present application, the prediction model migration correction of the laser-induced breakdown spectroscopy device under different focusing distances can be realized, thereby improving the accuracy and stability of the prediction model. Laser-induced breakdown spectroscopy equipment can quickly adapt to changes in on-site test conditions, solving the problem of prediction model failure caused by changes in focus distance due to changes in the position of detected materials at the application site. There is no need to recalibrate, saving manpower, material resources and time.
具体地,模型校正方法基于的激光诱导击穿光谱(LIBS)技术是一种成分分析技术,基于原子光谱和离子光谱的波长与特定的元素一一对应的关系,而且光谱信号强度与对应元素的含量也具有一定的量化关系,通过高能脉冲激光聚焦样品表面形成等离子体,用光谱仪记录等离子体发射的光谱信息,通过解析等离子光谱,并由特征波长的位置和光谱强度来对样品中的元素进行定性及定量分析。Specifically, the laser-induced breakdown spectroscopy (LIBS) technology based on the model calibration method is a composition analysis technology based on the one-to-one correspondence between the wavelengths of atomic spectra and ion spectra and specific elements, and the spectral signal intensity is related to the corresponding element. The content also has a certain quantitative relationship. Plasma is formed by focusing the surface of the sample with a high-energy pulse laser. A spectrometer is used to record the spectral information emitted by the plasma. By analyzing the plasma spectrum, the elements in the sample are analyzed based on the position of the characteristic wavelength and the spectral intensity. Qualitative and quantitative analysis.
首先模型校正方法采取步骤S10可利用LIBS技术设备对样品在多个聚焦距离的原始光谱数据进行获取。First, the model calibration method takes step S10, which can use LIBS technology equipment to acquire the original spectral data of the sample at multiple focusing distances.
样品可以是金属样品,例如钢材、合金等,样品的数量可以是多个,多个样品可以具有同一元素类型的多个样品。多个聚焦距离可以是LIBS技术设备的光学系统聚焦到样品的距离。示例性地,如图3所示,采用LIBS技术的设备对样品在两个聚焦距离获取原始光谱数据,当样品位于a处时,样品至LIBS技术设备的光学系统的聚焦距离为l 1,对处于a处的聚焦距离为l 1的样品获取原始光谱数据,然后可将样品置于b处,再对处于b处的聚焦距离为l 2的样品获取原始光谱数据。 The sample may be a metal sample, such as steel, alloy, etc., the number of samples may be multiple, and multiple samples may have multiple samples of the same element type. Multiple focusing distances may be the distance at which the optical system of the LIBS technology device focuses on the sample. For example, as shown in Figure 3, equipment using LIBS technology acquires raw spectral data from the sample at two focusing distances. When the sample is located at a, the focusing distance from the sample to the optical system of the LIBS technology equipment is l 1 , and for The original spectral data is obtained from the sample at position a with a focusing distance of l 1 , and then the sample can be placed at position b, and then the original spectral data is obtained from the sample at position b with a focusing distance of l 2 .
然后,模型校正方法采取步骤S20可利用步骤S10获取到的多个聚焦距离其中的一个聚焦距离的原始光谱数据建立预测模型。Then, the model correction method takes step S20 to establish a prediction model using the original spectral data of one of the multiple focusing distances obtained in step S10.
然后,模型校正方法采取步骤S30通过步骤S20建立的预测模型可基于其中一个聚焦距离的原始光谱数据来预测多个聚焦距离中其余的聚焦距离的光谱数据。Then, the model correction method takes step S30. The prediction model established in step S20 can predict the spectral data of the remaining focusing distances among the plurality of focusing distances based on the original spectral data of one of the focusing distances.
再然后,模型校正方法采取步骤S40基于步骤S30中多个聚焦距离的预测光谱数据和步骤S10中多个聚焦距离的原始光谱数据之间的误差,对步骤S20建立的预测模型校正。Then, the model correction method takes step S40 to correct the prediction model established in step S20 based on the error between the predicted spectral data of multiple focusing distances in step S30 and the original spectral data of multiple focusing distances in step S10.
综上所述,用LIBS技术设备采集一定聚焦距离的样品的光谱数据的实施例来进一步说 明:In summary, the embodiment of using LIBS technology equipment to collect spectral data of samples at a certain focusing distance will be further explained:
样品为碳钢,样品数量为10块,样品编号为0到9。聚焦距离为152.5cm至376.5cm之间的11个聚焦距离。The sample is carbon steel, the sample quantity is 10 pieces, and the sample number is 0 to 9. There are 11 focusing distances between 152.5cm and 376.5cm.
首先如步骤S10,可采用LIBS技术设备分别获取10块碳钢样品在聚焦距离152.5cm至376.5cm间的11个聚焦距离的原始光谱数据。First, as in step S10, LIBS technology equipment can be used to obtain raw spectral data of 11 carbon steel samples at 11 focusing distances between 152.5cm and 376.5cm.
然后如步骤S20,可利用在聚焦距离152.5cm处采集的10块碳钢样品的原始光谱数据,建立可以预测光谱数据的预测模型。Then in step S20, the original spectral data of 10 carbon steel samples collected at a focusing distance of 152.5 cm can be used to establish a prediction model that can predict the spectral data.
然后如步骤S30,可根据预测模型预测聚焦距离152.5cm至376.5cm之间11个聚焦距离的10块碳钢的预测光谱数据,进一步可以得到碳钢中锰(Mn)元素的浓度定量分析数据,结果如图4所示。Then in step S30, the predicted spectral data of 10 pieces of carbon steel with 11 focusing distances between 152.5cm and 376.5cm can be predicted according to the prediction model, and further the quantitative analysis data of the concentration of manganese (Mn) element in the carbon steel can be obtained. The results are shown in Figure 4.
然后如步骤S40,计算每块碳钢在聚焦距离152.5cm至376.5cm的11个聚焦距离的预测值和原始值之间的误差,如图5所示是其中1块碳钢在聚焦距离152.5cm至376.5cm的11个聚焦距离的Mn元素浓度预测值和原始值的误差。可根据计算的误差对聚焦距离152.5cm处采集的10块碳钢样品的原始光谱数据所建立的预测模型校正。Then in step S40, calculate the error between the predicted value and the original value of the 11 focusing distances of each piece of carbon steel at the focusing distance of 152.5cm to 376.5cm. As shown in Figure 5, one piece of carbon steel is at the focusing distance of 152.5cm. The error between the predicted value and the original value of Mn element concentration for 11 focusing distances to 376.5cm. The prediction model established on the original spectral data of 10 carbon steel samples collected at a focusing distance of 152.5cm can be corrected based on the calculated error.
请参阅图6,在某些实施方式中,基于预测光谱数据和原始光谱数据之间的误差,对多个聚焦距离的预测模型校正(步骤S40),包括:Referring to Figure 6, in some embodiments, based on the error between the predicted spectral data and the original spectral data, the prediction model correction (step S40) for multiple focusing distances includes:
步骤S41:建立原始光谱数据与预测光谱数据的校正函数;Step S41: Establish a correction function between the original spectral data and the predicted spectral data;
步骤S42:基于校正函数,建立损失函数与校正函数之间的关系;Step S42: Based on the correction function, establish the relationship between the loss function and the correction function;
步骤S43:基于关系,确定校正函数的数值;Step S43: Based on the relationship, determine the value of the correction function;
步骤S44:利用确定的校正函数的数值对多个聚焦距离的预测模型校正。Step S44: Use the determined values of the correction function to correct the prediction models of multiple focus distances.
在某些实施方式中,处理器31用于建立原始光谱数据与预测光谱数据的校正函数;及用于基于校正函数,建立损失函数与校正函数之间的关系;及用于基于关系,确定校正函数的数值;以及用于利用确定的校正函数的数值对多个聚焦距离的预测模型校正。In some embodiments, the processor 31 is configured to establish a correction function of the original spectral data and the predicted spectral data; and to establish a relationship between the loss function and the correction function based on the correction function; and to determine the correction based on the relationship. a numerical value of the function; and a predictive model correction for a plurality of focus distances using the determined numerical value of the correction function.
如此,利用函数关系建立校正函数、损失函数以及校正函数与损失函数之间的关系对预测模型校正,提高预测模型的准确性与稳定性。In this way, the functional relationship is used to establish the correction function, the loss function, and the relationship between the correction function and the loss function to correct the prediction model and improve the accuracy and stability of the prediction model.
具体地,模型校正方法中得到预测光谱数据和原始光谱数据之间的误差数据后可对多个聚焦距离的预测模型校正。Specifically, after obtaining the error data between the predicted spectral data and the original spectral data in the model correction method, the prediction model for multiple focusing distances can be corrected.
模型校正方法为实现步骤S40可先采取步骤S41,利用原始光谱数据与预测光谱数据之间的误差建立原始光谱数据与预测光谱数据的校正函数,然后可采取步骤S42,基于步骤S41建立的校正函数,建立损失函数与校正函数之间的关系,损失函数可以是利用原理函数建立的聚焦距离与误差的拟合函数并构建损失函数,误差为原始光谱数据与预测光谱数据之间的误差数据。In order to implement step S40, the model correction method can first take step S41, using the error between the original spectrum data and the predicted spectrum data to establish a correction function of the original spectrum data and the predicted spectrum data, and then take step S42, based on the correction function established in step S41. , establish the relationship between the loss function and the correction function. The loss function can be a fitting function of the focusing distance and error established using the principle function and construct a loss function. The error is the error data between the original spectral data and the predicted spectral data.
然后可采取步骤S43,基于步骤S42建立的损失函数与校正函数之间的关系可得出校正函数的具体数值,校正函数的具体数值可以是损失函数与校正函数之间的函数关系求解得到。然后再采取步骤S44,根据步骤S43中得到的校正函数的具体数值带入多个聚焦距离的对应预测模型中对预测模型进行校正。Then step S43 can be taken, and the specific value of the correction function can be obtained based on the relationship between the loss function and the correction function established in step S42. The specific value of the correction function can be obtained by solving the functional relationship between the loss function and the correction function. Then step S44 is taken, and the specific values of the correction function obtained in step S43 are brought into corresponding prediction models of multiple focus distances to correct the prediction model.
请参阅图7,在某些实施方式中,建立原始光谱数据与预测光谱数据的校正函数(步骤S41),包括:Referring to Figure 7, in some embodiments, establishing a correction function between original spectral data and predicted spectral data (step S41) includes:
步骤S411:基于原始光谱数据确定元素的原始浓度;Step S411: Determine the original concentration of the element based on the original spectral data;
步骤S412:基于预测光谱数据确定元素的预测浓度;Step S412: Determine the predicted concentration of the element based on the predicted spectrum data;
步骤S413:基于原始浓度和预测浓度建立校正函数。Step S413: Establish a correction function based on the original concentration and the predicted concentration.
在某些实施方式中,处理器31用于基于原始光谱数据确定元素的原始浓度;及用于基于预测光谱数据确定元素的预测浓度;及用于基于原始浓度和预测浓度建立校正函数。In certain embodiments, the processor 31 is configured to determine an original concentration of the element based on the original spectral data; and to determine a predicted concentration of the element based on the predicted spectral data; and to establish a correction function based on the original concentration and the predicted concentration.
如此,利用光谱数据与元素浓度的关系,将根据原始光谱数据和预测光谱数据建立校正函数转化为根据原始浓度与预测浓度建立校正函数,提高预测模型的准确性与稳定性。In this way, using the relationship between spectral data and element concentration, the correction function established based on the original spectral data and predicted spectral data is converted into a correction function established based on the original concentration and predicted concentration, thereby improving the accuracy and stability of the prediction model.
具体地,模型校正方法为实现步骤S41可先采取步骤S411,基于步骤S10获取的原始光谱数据与所对应聚焦距离检测的样品的某种元素之间的关系来确定对应聚焦距离检测的样品中该元素的原始浓度。然后可采取步骤S412,基于步骤S30预测的预测光谱数据与所对应聚焦距离检测的样品的某种元素之间的关系来确定对应聚焦距离检测的样品中该元素的预测浓度。再然后可采取步骤S413,基于步骤S411的原始浓度和步骤S412的预测浓度建立对应聚焦距离的校正函数。Specifically, in order to implement step S41, the model calibration method may first take step S411 to determine the element in the sample detected at the corresponding focus distance based on the relationship between the original spectral data obtained in step S10 and a certain element of the sample detected at the corresponding focus distance. The original concentration of an element. Step S412 may then be taken to determine the predicted concentration of the element in the sample detected at the corresponding focus distance based on the relationship between the predicted spectrum data predicted in step S30 and a certain element in the sample detected at the corresponding focus distance. Then step S413 may be taken to establish a correction function corresponding to the focus distance based on the original concentration in step S411 and the predicted concentration in step S412.
在某些实施方式中,校正函数采用以下关系式:In some embodiments, the correction function takes the following relationship:
C=Z(l)cC=Z(l)c
其中,c为预测模型在聚焦距离l处的元素的预测浓度,C为元素的原始浓度,Z(l)为聚焦距离l处对应的校正函数。Among them, c is the predicted concentration of the element at the focusing distance l of the prediction model, C is the original concentration of the element, and Z(l) is the corresponding correction function at the focusing distance l.
如此,元素的预测浓度建立在元素的原始浓度与校正函数的关系上,校正函数可以通过元素的原始浓度预测处于其他聚焦距离的元素的预测浓度,实现了不同聚焦距离下的预测模型的迁移校正。In this way, the predicted concentration of an element is based on the relationship between the original concentration of the element and the correction function. The correction function can predict the predicted concentration of elements at other focusing distances through the original concentration of the element, realizing the migration correction of the prediction model at different focusing distances. .
具体地,校正函数的关系式确定可来源与光谱数据和元素浓度之间的关系,LIBS技术基于原子光谱和离子光谱的波长与特定的元素一一对应的关系,主要涉及光谱信息、高能脉冲激光、等离子体、特征波长等数据。Specifically, the relational expression of the correction function determines the relationship between the source, spectral data and element concentration. LIBS technology is based on the one-to-one correspondence between the wavelengths of atomic spectra and ion spectra and specific elements, and mainly involves spectral information, high-energy pulse lasers , plasma, characteristic wavelength and other data.
校正函数采用的关系式可由具体公式推导得到:The relational expression used in the correction function can be derived from the specific formula:
由于激光脉冲功率密度公式为Since the laser pulse power density formula is
Figure PCTCN2022106217-appb-000002
Figure PCTCN2022106217-appb-000002
式中,PD为激光功率密度,E为脉冲能量,w为脉冲宽度,d为聚焦光斑直径;In the formula, PD is the laser power density, E is the pulse energy, w is the pulse width, and d is the focused spot diameter;
而等离子体温度T与脉冲功率密度PD有关,则等离子体温度可表示为T=f(E,w,d),其中脉冲能量为E=g(l),带入等离子体绝对强度公式The plasma temperature T is related to the pulse power density PD, then the plasma temperature can be expressed as T = f (E, w, d), where the pulse energy is E = g (l), which is brought into the plasma absolute intensity formula
Figure PCTCN2022106217-appb-000003
中,则
Figure PCTCN2022106217-appb-000003
in the middle, then
Figure PCTCN2022106217-appb-000004
Figure PCTCN2022106217-appb-000004
式中w为脉冲宽度,d为聚焦光斑直径,i代表i能级,j代表j能级,A ij为跃迁几率,g i为高能级简并度,λ ij为辐射波长,U为当前温度下该离子对应的匹配函数,E i为高能级能量,K B为玻尔兹曼常数,其中F为与系统参数、等离子体温度、元素特性等相关的比例因子,C为元素的原始浓度,I ij为元素特征; In the formula, w is the pulse width, d is the focused spot diameter, i represents the i energy level, j represents the j energy level, A ij is the transition probability, g i is the high energy level degeneracy, λ ij is the radiation wavelength, and U is the current temperature. The matching function corresponding to the ion is the following, E i is the high energy level energy, K B is the Boltzmann constant, where F is the scaling factor related to system parameters, plasma temperature, element characteristics, etc., C is the original concentration of the element, I ij is the element characteristic;
当前测量参数一致,在忽略自吸收影响时,上述公式中c前的各系数可记为常数a,因此公式可改写为:The current measurement parameters are consistent. When the influence of self-absorption is ignored, each coefficient before c in the above formula can be recorded as a constant a, so the formula can be rewritten as:
I=aCI=aC
C=I/aC=I/a
其中I为元素特征,C为元素的原始浓度;Among them, I is the characteristic of the element, and C is the original concentration of the element;
但是当前实验条件存在聚焦距离的变化,因此公式改写为However, the current experimental conditions have changes in focusing distance, so the formula is rewritten as
C=Z(l)cC=Z(l)c
其中,c为预测模型在聚焦距离l处的元素的预测浓度,C为元素的原始浓度,Z(l)为聚焦距离l处对应的校正函数。Among them, c is the predicted concentration of the element at the focusing distance l of the prediction model, C is the original concentration of the element, and Z(l) is the corresponding correction function at the focusing distance l.
请参阅图8,在某些实施方式中,基于关系,确定校正函数的数值(步骤S43),包括:Referring to Figure 8, in some embodiments, based on the relationship, determining the value of the correction function (step S43) includes:
步骤S431:利用最小二乘法确定损失函数的最小值;Step S431: Use the least squares method to determine the minimum value of the loss function;
步骤S432:利用最小值确定校正函数的数值。Step S432: Use the minimum value to determine the value of the correction function.
在某些实施方式中,处理器31用于利用最小二乘法确定损失函数的最小值;以及用于利用最小值确定校正函数的数值。In some embodiments, the processor 31 is configured to determine a minimum value of the loss function using a least squares method; and is configured to determine a value of the correction function using the minimum value.
如此,确定损失函数的最小值时采用最小二乘法的方式能够便于对损失函数的最小值 求解,从而确定校正函数的数值。In this way, when determining the minimum value of the loss function, the least squares method can be used to solve the minimum value of the loss function, thereby determining the value of the correction function.
具体地,模型校正方法为实现步骤S43可先依据步骤S42中构建得到的损失函数采取步骤S431,将损失函数通过带入数据并采用最小二乘法的求解方法得到损失函数的最小值,然后可采取步骤S432,由于损失函数与校正函数之间存在关系式,通过最小二乘法求解损失函数的最小值可进一步得到校正函数的具体数值。Specifically, in order to implement step S43, the model correction method can first take step S431 based on the loss function constructed in step S42, bring the loss function into the data and use the least squares method to obtain the minimum value of the loss function, and then take the step S431. Step S432: Since there is a relationship between the loss function and the correction function, the specific value of the correction function can be further obtained by solving the minimum value of the loss function through the least squares method.
在某些实施方式中,损失函数采用以下关系式:In some implementations, the loss function adopts the following relationship:
Figure PCTCN2022106217-appb-000005
Figure PCTCN2022106217-appb-000005
其中,J为损失函数,i为样品编号,j为聚焦距离编号,M为样品数量,N为聚焦距离数量,C i为第i个样品的元素的原始浓度,Z(l j)为第j个聚焦距离l j处对应的校正函数,c j i为预测模型的第i个样品在第j个聚焦距离l j处的元素的预测浓度。 Among them, J is the loss function, i is the sample number, j is the focusing distance number, M is the number of samples, N is the number of focusing distances, C i is the original concentration of the element of the i-th sample, Z(l j ) is the j-th The correction function corresponding to the focusing distance l j , c j i is the predicted concentration of the element at the jth focusing distance l j of the i-th sample of the prediction model.
如此,损失函数的关系式中包括了多组样品在多个聚焦距离的元素原始浓度和元素预测浓度,使得损失函数关系式可以实现数据的拟合。In this way, the loss function relationship includes the original concentration and predicted element concentration of multiple groups of samples at multiple focusing distances, so that the loss function relationship can achieve data fitting.
具体地,可以理解,损失函数的关系式可以是样品的根据预测模型的关联校正函数得到的元素的预测浓度,与元素的原始浓度之间的差值平方所构建的函数关系式。损失函数的求解可采用最小二乘法的原理,通过使损失函数达到最小值,可对损失函数求解。Specifically, it can be understood that the relational expression of the loss function can be a functional relational expression constructed from the square of the difference between the predicted concentration of the element obtained according to the correlation correction function of the prediction model and the original concentration of the element. The loss function can be solved by using the least squares method. By making the loss function reach the minimum value, the loss function can be solved.
示例性地,当样品数量为10个,LIBS技术设备对样品获取的聚焦距离的数量为11个时,可将样品编号从1编号至10,不同聚焦距离由1编号至11,由此,损失函数可以采用
Figure PCTCN2022106217-appb-000006
Figure PCTCN2022106217-appb-000007
式,再将损失函数中的对应样品编号的样品在对应的聚焦距离编号上获取的元素的原始浓度、元素的预测浓度等相关数据代入损失函数中,通过最小二乘法对损失函数求解最小值可得到对应的校正函数的具体数值,然后可通过确定的多个聚焦距离的校正函数的数值对多个聚焦距离的预测模型进行校正。
For example, when the number of samples is 10 and the number of focusing distances obtained by the LIBS technology equipment for the samples is 11, the samples can be numbered from 1 to 10, and the different focusing distances can be numbered from 1 to 11. Therefore, the loss The function can take
Figure PCTCN2022106217-appb-000006
Figure PCTCN2022106217-appb-000007
formula, and then substitute the original concentration of the element, the predicted concentration of the element and other relevant data of the sample corresponding to the sample number in the loss function at the corresponding focusing distance number into the loss function, and solve the minimum value of the loss function through the least squares method. The specific values of the corresponding correction functions are obtained, and then the prediction models of the multiple focus distances can be corrected through the determined values of the correction functions of the multiple focus distances.
请参阅图1,在某些实施方式中,模型校正方法包括:Referring to Figure 1, in some embodiments, the model correction method includes:
步骤S50:利用校正后的预测模型对相应聚焦距离上的样品中的元素进行预测。Step S50: Use the corrected prediction model to predict the elements in the sample at the corresponding focusing distance.
在某些实施方式中,处理器31用于利用校正后的预测模型对相应聚焦距离上的样品中的元素进行预测。In some embodiments, the processor 31 is configured to use the corrected prediction model to predict elements in the sample at the corresponding focusing distance.
如此,模型校正方法经校正后的预测模型相比校正前的预测模型对样品的光谱数据和元素浓度的预测更精确。In this way, the prediction model corrected by the model calibration method is more accurate in predicting the spectral data and element concentration of the sample than the prediction model before correction.
具体地,模型校正方法为实现对预测模型的校正可采用步骤S50,结合步骤S40中得到的校正函数的具体数值对步骤S20中建立的对应聚焦距离的预测模型进行校正,校正后的预测模型可再对多个聚焦距离的样品中的元素进行预测,可以得到更为精确的样品的元素测量数据。Specifically, the model correction method can use step S50 to correct the prediction model. The prediction model corresponding to the focus distance established in step S20 is corrected based on the specific value of the correction function obtained in step S40. The corrected prediction model can be By predicting the elements in the sample at multiple focusing distances, more accurate element measurement data of the sample can be obtained.
仍然采用上述实施例中LIBS技术设备采集一定聚焦距离的样品的光谱数据的来进一步说明:The LIBS technology equipment in the above embodiment is still used to collect spectral data of samples at a certain focusing distance to further explain:
样品为碳钢,样品数量为10块,样品编号为0到9。聚焦距离为152.5cm至376.5cm之间的11个聚焦距离。The sample is carbon steel, the sample quantity is 10 pieces, and the sample number is 0 to 9. There are 11 focusing distances between 152.5cm and 376.5cm.
在步骤S30中,根据预测模型预测聚焦距离152.5cm至376.5cm之间11个聚焦距离的10块碳钢的预测光谱数据,得到的碳钢中锰(Mn)元素的浓度定量分析数据,结果为图4。In step S30, the prediction model predicts the predicted spectrum data of 10 pieces of carbon steel at 11 focusing distances between 152.5cm and 376.5cm, and obtains the quantitative analysis data of the concentration of manganese (Mn) element in the carbon steel. The result is: Figure 4.
在步骤S50中,根据校正后的预测模型对相应聚焦距离上的样品中的元素进行预测,根据校正后的预测模型对预测聚焦距离152.5cm至376.5cm之间11个聚焦距离的10块碳钢的预测光谱数据进行校正,得到的碳钢中锰(Mn)元素的校正后的浓度定量分析数据,结果为图9。In step S50, the elements in the sample at the corresponding focusing distance are predicted according to the corrected prediction model, and 10 pieces of carbon steel at 11 focusing distances between 152.5cm and 376.5cm are predicted according to the corrected prediction model. The predicted spectral data were corrected to obtain the corrected quantitative analysis data of the concentration of manganese (Mn) element in carbon steel. The results are shown in Figure 9.
可以理解,相比图4,经校正后的图9中预测模型预测的各编号的样品的在不同聚焦距离处预测得到的浓度差值大幅减小,因此,模型校正方法可以实现了不同聚焦距离下对预测模型校正。It can be understood that compared with Figure 4, the predicted concentration difference of each numbered sample predicted by the prediction model in Figure 9 at different focusing distances is greatly reduced after correction. Therefore, the model correction method can achieve different focusing distances. Calibrate the prediction model below.
本申请实施方式的计算机设备包括存储器和处理器,存储器存储有计算机程序,处理器执行计算机程序时实现上述任一项实施方式中模型校正方法。The computer device in the embodiment of the present application includes a memory and a processor. The memory stores a computer program. When the processor executes the computer program, it implements the model correction method in any of the above embodiments.
本申请实施方式的计算机可执行指令的非易失性计算机可读存储介质,当计算机可执行指令被一个或多个处理器执行时,使得处理器执行上述任一项实施方式中模型校正方法。The non-volatile computer-readable storage medium of computer-executable instructions in the embodiment of the present application, when the computer-executable instructions are executed by one or more processors, causes the processor to execute the model correction method in any of the above-mentioned embodiments.
具体地,本申请实施方式中的计算机设备可以是计算器、可编程控制器、台式电脑、膝上型电脑、平板电脑、服务器等设备,计算机设备可以包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。Specifically, the computer device in the embodiment of the present application can be a calculator, a programmable controller, a desktop computer, a laptop computer, a tablet computer, a server, and other devices. The computer device can include a processor, a memory, and a computer connected through a system bus. Communication interfaces, displays and input devices.
计算机设备的处理器可以为中央处理器(Central Processing Unit,CPU)。处理器还可以为其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、 分立硬件组件等芯片,或者上述各类芯片的组合。The processor of the computer device may be a central processing unit (Central Processing Unit, CPU). The processor can also be other general-purpose processors, Digital Signal Processor (DSP), Application Specific Integrated Circuit (ASIC), Field-Programmable Gate Array (FPGA) or other Chips such as programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or combinations of these types of chips.
计算机程序可以被存储在存储器中,存储器作为一种非暂态计算机可读存储介质,可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块,如上述方法实施例中的方法所对应的程序指令/模块。处理器通过运行存储在存储器中的非暂态软件程序、指令以及模块,从而执行处理器的各种功能应用以及数据处理,即实现上述方法实施例中的方法。The computer program can be stored in the memory. The memory, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs, non-transitory computer executable programs and modules, as described in the method in the above method embodiment. Corresponding program instructions/modules. The processor executes various functional applications and data processing of the processor by running non-transient software programs, instructions and modules stored in the memory, that is, implementing the method in the above method embodiment.
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本申请的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本申请的实施例所属技术领域的技术人员所理解。Any process or method descriptions in flowcharts or otherwise described herein may be understood to represent modules, segments, or portions of code that include one or more executable instructions for implementing the specified logical functions or steps of the process. , and the scope of the preferred embodiments of the present application includes additional implementations in which functions may be performed out of the order shown or discussed, including in a substantially simultaneous manner or in the reverse order, depending on the functionality involved, which shall It should be understood by those skilled in the technical field to which the embodiments of this application belong.
在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理模块的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或多个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。The logic and/or steps represented in the flowcharts or otherwise described herein, for example, may be considered a sequenced list of executable instructions for implementing the logical functions, and may be embodied in any computer-readable medium, For use by, or in combination with, instruction execution systems, devices or equipment (such as computer-based systems, systems including processing modules, or other systems that can fetch instructions from and execute instructions from the instruction execution system, device or equipment) or equipment. For the purposes of this specification, a "computer-readable medium" may be any device that can contain, store, communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. More specific examples (non-exhaustive list) of computer readable media include the following: electrical connections with one or more wires (electronic device), portable computer disk cartridges (magnetic device), random access memory (RAM), Read-only memory (ROM), erasable and programmable read-only memory (EPROM or flash memory), fiber optic devices, and portable compact disc read-only memory (CDROM). Additionally, the computer-readable medium may even be paper or other suitable medium on which the program may be printed, as the paper or other medium may be optically scanned, for example, and subsequently edited, interpreted, or otherwise suitable as necessary. process to obtain the program electronically and then store it in computer memory.
处理器可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。The processor can be a Central Processing Unit (CPU), other general-purpose processors, Digital Signal Processor (DSP), Application Specific Integrated Circuit (ASIC), or off-the-shelf programmable processors. Gate array (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or the processor may be any conventional processor, etc.
应当理解,本申请的实施方式的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如,如果用硬件来实现,和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列 (PGA),现场可编程门阵列(FPGA)等。It should be understood that various parts of the embodiments of the present application can be implemented in hardware, software, firmware, or a combination thereof. In the above embodiments, various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if it is implemented in hardware, as in another embodiment, it can be implemented by any one or a combination of the following technologies known in the art: a logic gate circuit with a logic gate circuit for implementing a logic function on a data signal. Discrete logic circuits, application specific integrated circuits with suitable combinational logic gates, programmable gate arrays (PGA), field programmable gate arrays (FPGA), etc.
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。Those of ordinary skill in the art can understand that all or part of the steps involved in implementing the methods of the above embodiments can be completed by instructing relevant hardware through a program. The program can be stored in a computer-readable storage medium. The program can be stored in a computer-readable storage medium. When executed, one of the steps of the method embodiment or a combination thereof is included.
此外,在本申请的各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。In addition, each functional unit in various embodiments of the present application can be integrated into a processing module, each unit can exist physically alone, or two or more units can be integrated into one module. The above integrated modules can be implemented in the form of hardware or software function modules. If the integrated module is implemented in the form of a software function module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.
上述提到的存储介质可以是只读存储器,磁盘或光盘等。The storage media mentioned above can be read-only memory, magnetic disks or optical disks, etc.
尽管上面已经示出和描述了本申请的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本申请的限制,本领域的普通技术人员在本申请的范围内可以对上述实施方式进行变化、修改、替换和变型。Although the embodiments of the present application have been shown and described above, it can be understood that the above-mentioned embodiments are illustrative and cannot be understood as limitations of the present application. Those of ordinary skill in the art can make modifications to the above-mentioned embodiments within the scope of the present application. The embodiments are subject to changes, modifications, substitutions and variations.

Claims (10)

  1. 一种模型校正方法,其特征在于,包括:A model correction method, characterized by including:
    基于激光诱导击穿光谱技术,获取多个聚焦距离的原始光谱数据;Based on laser-induced breakdown spectroscopy technology, original spectral data at multiple focusing distances are obtained;
    利用其中一个聚焦距离的原始光谱数据建立预测模型;Build a predictive model using raw spectral data at one of the focusing distances;
    利用所述预测模型,预测其余聚焦距离的预测光谱数据;Using the prediction model, predict the predicted spectral data for the remaining focus distances;
    基于所述预测光谱数据和所述原始光谱数据之间的误差,对多个聚焦距离的预测模型校正。Predictive model corrections for multiple focus distances are based on errors between the predicted spectral data and the original spectral data.
  2. 根据权利要求1所述的模型校正方法,其特征在于,基于所述预测光谱数据和所述原始光谱数据之间的误差,对多个聚焦距离的预测模型校正,包括:The model correction method according to claim 1, characterized in that, based on the error between the predicted spectral data and the original spectral data, the prediction model correction for multiple focusing distances includes:
    建立原始光谱数据与预测光谱数据的校正函数;Establish a correction function between original spectral data and predicted spectral data;
    基于所述校正函数,建立损失函数与所述校正函数之间的关系;Based on the correction function, establish a relationship between the loss function and the correction function;
    基于所述关系,确定所述校正函数的数值;Based on the relationship, determining a value of the correction function;
    利用确定的所述校正函数的数值对多个聚焦距离的预测模型校正。The prediction model of multiple focus distances is corrected using the determined value of the correction function.
  3. 根据权利要求2所述的模型校正方法,其特征在于,建立原始光谱数据与预测光谱数据的校正函数,包括:The model correction method according to claim 2, characterized in that establishing a correction function between original spectral data and predicted spectral data includes:
    基于所述原始光谱数据确定元素的原始浓度;determining an original concentration of an element based on the original spectral data;
    基于所述预测光谱数据确定元素的预测浓度;determining a predicted concentration of an element based on the predicted spectral data;
    基于所述原始浓度和所述预测浓度建立所述校正函数。The correction function is established based on the original concentration and the predicted concentration.
  4. 根据权利要求3所述的模型校正方法,其特征在于,所述校正函数采用以下关系式:The model correction method according to claim 3, characterized in that the correction function adopts the following relationship:
    C=Z(l)cC=Z(l)c
    其中,c为预测模型在聚焦距离l处的元素的预测浓度,C为元素的原始浓度,Z(l)为聚焦距离l处对应的校正函数。Among them, c is the predicted concentration of the element at the focusing distance l of the prediction model, C is the original concentration of the element, and Z(l) is the corresponding correction function at the focusing distance l.
  5. 根据权利要求2所述的模型校正方法,其特征在于,基于所述关系,确定所述校正函数的数值,包括:The model correction method according to claim 2, characterized in that, based on the relationship, determining the value of the correction function includes:
    利用最小二乘法确定所述损失函数的最小值;Determine the minimum value of the loss function using the least squares method;
    利用所述最小值确定所述校正函数的数值。The minimum value is used to determine the value of the correction function.
  6. 根据权利要求2所述的模型校正方法,其特征在于,所述损失函数采用以下关系式:The model correction method according to claim 2, characterized in that the loss function adopts the following relationship:
    Figure PCTCN2022106217-appb-100001
    Figure PCTCN2022106217-appb-100001
    其中,J为损失函数,i为样品编号,j为聚焦距离编号,M为样品数量,N为聚焦距离数量,C i为第i个样品的元素的原始浓度,Z(l j)为第j个聚焦距离l j处对应的校正函数,c j i为预测模型的第i个样品在第j个聚焦距离l j处的元素的预测浓度。 Among them, J is the loss function, i is the sample number, j is the focusing distance number, M is the number of samples, N is the number of focusing distances, C i is the original concentration of the element of the i-th sample, Z(l j ) is the j-th The correction function corresponding to the focusing distance l j , c j i is the predicted concentration of the element at the jth focusing distance l j of the i-th sample of the prediction model.
  7. 根据权利要求1所述的模型校正方法,其特征在于,所述模型校正方法包括:The model correction method according to claim 1, characterized in that the model correction method includes:
    利用校正后的预测模型对相应聚焦距离上的样品中的元素进行预测。The corrected prediction model is used to predict the elements in the sample at the corresponding focusing distance.
  8. 一种光谱设备,其特征在于,所述光谱设备包括:A kind of spectroscopic equipment, characterized in that the spectroscopic equipment includes:
    样品台,用于承载样品;Sample stage, used to carry samples;
    激光器,用于向所述样品发出激光;a laser for emitting laser light to the sample;
    光谱仪,所述光谱仪用于接收所述样品反射的激光,所述光谱仪包括处理器,所述处理器用于实现如权利要求1-7中任一项所述模型校正方法。A spectrometer, the spectrometer is used to receive the laser light reflected by the sample, the spectrometer includes a processor, the processor is used to implement the model correction method according to any one of claims 1-7.
  9. 一种计算机设备,其特征在于,所述计算机设备包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现如权利要求1-7中任一项所述模型校正方法。A computer device, characterized in that the computer device includes a memory and a processor, the memory stores a computer program, and when the processor executes the computer program, it implements the method described in any one of claims 1-7 Model calibration method.
  10. 一种计算机可执行指令的非易失性计算机可读存储介质,其特征在于,当所述计算机可执行指令被一个或多个处理器执行时,使得所述处理器执行权利要求1-7中任一项所述模型校正方法。A non-volatile computer-readable storage medium of computer-executable instructions, characterized in that when the computer-executable instructions are executed by one or more processors, the processor is caused to execute claims 1-7 Any one of the model correction methods.
PCT/CN2022/106217 2022-05-18 2022-07-18 Model correction method, spectroscopy device, computer device, and storage medium WO2023221282A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210538088.7A CN114636688B (en) 2022-05-18 2022-05-18 Model correction method, spectroscopic apparatus, computer apparatus, and storage medium
CN202210538088.7 2022-05-18

Publications (1)

Publication Number Publication Date
WO2023221282A1 true WO2023221282A1 (en) 2023-11-23

Family

ID=81952767

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/106217 WO2023221282A1 (en) 2022-05-18 2022-07-18 Model correction method, spectroscopy device, computer device, and storage medium

Country Status (2)

Country Link
CN (1) CN114636688B (en)
WO (1) WO2023221282A1 (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150346103A1 (en) * 2014-05-29 2015-12-03 Bwt Property, Inc. Laser Induced Breakdown Spectroscopy (LIBS) Apparatus and Method for Performing Spectral Imaging of a Sample Surface
US20170219494A1 (en) * 2016-02-01 2017-08-03 Bwt Property, Inc. Laser Induced Breakdown Spectroscopy (LIBS) Apparatus with Automatic Wavelength Calibration
CN107817239A (en) * 2016-09-13 2018-03-20 中国科学院沈阳自动化研究所 A kind of LIBS spectral correction methods based on plasma position information
CN109358036A (en) * 2018-12-20 2019-02-19 中国科学院光电研究院 Laser induced breakdown spectroscopy signal errors corrects system and method
CN111272735A (en) * 2020-03-26 2020-06-12 中国科学院空天信息创新研究院 Detection method of laser-induced breakdown spectroscopy
CN111289496A (en) * 2020-03-26 2020-06-16 中国科学院空天信息创新研究院 Detection method and device for long-distance zoom laser-induced breakdown spectroscopy
CN114324302A (en) * 2022-03-08 2022-04-12 合肥金星智控科技股份有限公司 Automatic positioning method and system for laser focusing position
CN114509425A (en) * 2022-01-28 2022-05-17 山东大学 Material composition surface analysis method based on variable defocus amount laser-induced breakdown spectroscopy

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101271063A (en) * 2008-05-07 2008-09-24 陕西恒智科技发展有限公司 Process of vector machine emendation model method supported by gas infrared spectrum analysis
EA201390408A1 (en) * 2010-10-01 2013-10-30 Текнолоджикал Ресорсиз Пти. Лимитед LASER-SPARK EMISSION SPECTROSCOPIC ANALYZER
CA2950163A1 (en) * 2016-11-30 2018-05-30 National Research Council Of Canada Method and system for characterizing oil sands by using laser-induced breakdown spectroscopy
CN106546576B (en) * 2016-12-06 2019-02-26 山西大学 A kind of laser induced breakdown spectroscopy bearing calibration based on homogeneous substance
CN206945540U (en) * 2017-07-14 2018-01-30 中国科学技术大学 A kind of efficient detection device based on Laser-induced plasma spectroscopy
CN111044504B (en) * 2019-12-16 2021-03-30 华南理工大学 Coal quality analysis method considering uncertainty of laser-induced breakdown spectroscopy
CN111488926B (en) * 2020-04-07 2024-03-29 中国科学院南京土壤研究所 Soil organic matter determination method based on optimization model
CN113204879B (en) * 2021-05-07 2022-11-01 西华大学 Improved Hankel matrix prediction model modeling method based on fluorescent oil film and application

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150346103A1 (en) * 2014-05-29 2015-12-03 Bwt Property, Inc. Laser Induced Breakdown Spectroscopy (LIBS) Apparatus and Method for Performing Spectral Imaging of a Sample Surface
US20170219494A1 (en) * 2016-02-01 2017-08-03 Bwt Property, Inc. Laser Induced Breakdown Spectroscopy (LIBS) Apparatus with Automatic Wavelength Calibration
CN107817239A (en) * 2016-09-13 2018-03-20 中国科学院沈阳自动化研究所 A kind of LIBS spectral correction methods based on plasma position information
CN109358036A (en) * 2018-12-20 2019-02-19 中国科学院光电研究院 Laser induced breakdown spectroscopy signal errors corrects system and method
CN111272735A (en) * 2020-03-26 2020-06-12 中国科学院空天信息创新研究院 Detection method of laser-induced breakdown spectroscopy
CN111289496A (en) * 2020-03-26 2020-06-16 中国科学院空天信息创新研究院 Detection method and device for long-distance zoom laser-induced breakdown spectroscopy
CN114509425A (en) * 2022-01-28 2022-05-17 山东大学 Material composition surface analysis method based on variable defocus amount laser-induced breakdown spectroscopy
CN114324302A (en) * 2022-03-08 2022-04-12 合肥金星智控科技股份有限公司 Automatic positioning method and system for laser focusing position

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
CHAOYANG LI, AN NING; PAN CONGYUAN: "Application of Laser-Induced Breakdown Spectroscopy to", ANHUI SCIENCE & TECHNOLOGY, vol. 2, 20 February 2021 (2021-02-20), pages 46 - 48, XP093108154 *

Also Published As

Publication number Publication date
CN114636688A (en) 2022-06-17
CN114636688B (en) 2022-08-12

Similar Documents

Publication Publication Date Title
JP4247559B2 (en) X-ray fluorescence analyzer and program used therefor
US10551325B2 (en) Systems for parsing material properties from within SHG signals
EP2710353B1 (en) SPECTROSCOPIC APPARATUS AND METHOD of DETERMINING COMPONENTS PRESENT IN A SAMPLE
CN109738881B (en) Calibration method and device of time-of-flight depth module and readable storage medium
JP2023505380A (en) Systems and methods for analyzing unknown sample composition using predictive models based on emission spectra
JP2006058294A (en) Inspection method of pattern formed on substrate, and inspection apparatus for implementing same
KR102589702B1 (en) Method and apparatus for evaluating statistically distributed measured values in the inspection of elements of a photolithographic process
WO2023221282A1 (en) Model correction method, spectroscopy device, computer device, and storage medium
JP7482782B2 (en) Sequence-based protein structure and properties determination
Liu et al. Long-term repeatability improvement using beam intensity distribution for laser-induced breakdown spectroscopy
TWI805876B (en) Loosely coupled inspection and metrology system for high-volume production process monitoring
JP4634413B2 (en) measuring device
US20160131526A1 (en) Spectroscopic Analysis System and Method
CN110031416B (en) Gas concentration detection device and method
CN117216558A (en) Coal quality analysis model training method and device, electronic equipment and storage medium
CN115753740A (en) Underwater ore in-situ detection method and system based on combination of LIBS and pulse Raman
CN114660048B (en) Predictive value correction method, spectroscopic apparatus, computer apparatus, and storage medium
JP5248249B2 (en) Particle size measuring apparatus, particle size measuring method, and computer program
WO2019113955A1 (en) Sub-film defect detection method and sub-film defect detection device
CN114688964A (en) Critical dimension measurement correction method, system and computer readable storage medium
CN111965166A (en) Rapid measurement method for biomass briquette characteristic index
JP2006275901A (en) Device and method for crystal evaluation
US20170212050A1 (en) Method for calibrating and using calibration of microplates
JP7190751B2 (en) X-ray fluorescence analyzer
JP2007322405A (en) Inspection apparatus

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22942313

Country of ref document: EP

Kind code of ref document: A1