WO2020186844A1 - Self-adaptive surface absorption spectrum analysis method and system, storage medium, and device - Google Patents

Self-adaptive surface absorption spectrum analysis method and system, storage medium, and device Download PDF

Info

Publication number
WO2020186844A1
WO2020186844A1 PCT/CN2019/125394 CN2019125394W WO2020186844A1 WO 2020186844 A1 WO2020186844 A1 WO 2020186844A1 CN 2019125394 W CN2019125394 W CN 2019125394W WO 2020186844 A1 WO2020186844 A1 WO 2020186844A1
Authority
WO
WIPO (PCT)
Prior art keywords
spectrum
blank
light intensity
absorption
approximation
Prior art date
Application number
PCT/CN2019/125394
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 WO2020186844A1 publication Critical patent/WO2020186844A1/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/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/3103Atomic absorption analysis

Definitions

  • the present invention relates to the field of spectrum analysis, in particular to an adaptive surface absorption spectrum analysis method.
  • Surface absorption spectroscopy detection is a non-destructive spectroscopic analysis technology that uses a wide-band light source to illuminate the surface of the sample, detect the reflected light from the surface of the sample, and determine the content of certain characteristic molecular groups in the sample by the absorption of incident light by the surface material. Because it does not need to destroy the sample morphology and structure, and the detection speed is fast, it is widely used in the field of industrial non-destructive testing. In recent years, it has also shown a good momentum of development in the field of medical in vitro non-destructive diagnosis.
  • the traditional surface absorption spectroscopy analysis method needs to use a sample that does not contain the characteristic substance or characteristic molecular group to be measured as a reference object, measure the surface reflectance spectrum of the reference object as a blank light, and then place the sample in the same relative position to measure its surface reflectance spectrum After dividing by the blank light, it is converted into a reflectance or absorbance curve for qualitative and quantitative analysis.
  • the disadvantage is that the change of the relative position between the sample and the optical probe of the instrument will easily have a greater impact on the measured characteristic peak intensity, and the interference of the absorption band before and after the wavelength position of the characteristic peak cannot be ruled out, making the detection result easy to appear Misjudgment affects the convenience and accuracy of detection.
  • the present invention provides an adaptive surface absorption spectroscopy analysis method.
  • the method uses the characteristics of a wide-band light source with continuous and gentle light intensity bands and slow changes in substrate absorption bands, and avoids potential problems after collecting sample surface reflection spectra.
  • Narrowband absorption peak position using two light intensity bands at a certain distance before and after the potential narrowband absorption peak position to fit and restore a light intensity curve without narrowband absorption as a blank spectrum, and then divide the sample surface reflectance spectrum curve and the blank spectrum curve Logarithm to obtain the surface absorption spectrum curve of the sample.
  • the present invention provides an adaptive surface absorption spectrum analysis method, which includes the following steps:
  • Blank spectrum extraction selecting part of the spectrum data in the reflection spectrum data, and extracting the currently detected blank spectrum by performing data processing on the part of the spectrum data; wherein the part of the spectrum data does not include potential absorption peak data;
  • step S2 further includes the following steps:
  • step S24 is further included:
  • the method further includes:
  • the measured light intensity signal value is a value obtained after smoothing and noise reduction processing.
  • step S1 further includes: irradiating the surface of the sample with a wide-band light source, and obtaining the reflection absorption spectrum of the sample surface with a spectrometer, and the intensity of the reflected light on the sample surface is greater than or equal to 5% of the full range of the spectrometer.
  • An electronic device including: a processor;
  • a memory and a program, wherein the program is stored in the memory and is configured to be executed by a processor, the program including an adaptive surface absorption spectroscopy analysis method.
  • a computer-readable storage medium has a computer program stored thereon, and the computer program is executed by a processor by an adaptive surface absorption spectrum analysis method.
  • the self-adaptive surface absorption spectrum analysis system includes a spectrum data acquisition module, an absorption spectrum calculation module and a blank spectrum extraction module; among them,
  • the acquiring spectrum data module is used to acquire the reflection absorption spectrum of the surface of the sample to obtain reflection spectrum data; wherein the reflection absorption spectrum is a spectrum diagram with wavelength as the abscissa and light intensity signal value as the ordinate;
  • the blank spectrum extraction module is used to select part of the spectrum data in the reflectance spectrum data, and extract the blank spectrum detected once by performing data processing on the part of the spectrum data; wherein the part of the spectrum data does not include potential absorption peak data ;
  • the blank spectrum extraction module includes a selection spectrum collection unit, a calculation unit, and a blank spectrum fitting unit; the selection spectrum collection unit is used to select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak in the reflection spectrum data, Forming a discrete combined spectrum set; the calculation unit is used to calculate the discrete combined spectrum set through the establishment function to obtain an approximation function expression;
  • the absorption spectrum calculation module is used to perform a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value obtained by using the blank spectrum and invert to obtain the absorbance value of the wavelength point.
  • the calculation unit includes a model residual sorting unit, a unit for establishing the optimal approximation condition equation group, and an approximation function solving unit; wherein,
  • the model residual arranging unit is used to establish an approximation function with the wavelength of a discrete point as the abscissa and the light intensity signal value as the ordinate, and the light intensity signal value obtained by the approximation function corresponds to the actual measured light intensity corresponding to the same discrete point
  • the function expression of the difference obtained from the signal value is recorded as the residual expression
  • the unit for establishing the best approximation condition equation group is used to perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the best approximation condition equation group;
  • the approximation function solving unit is used to calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and the coefficient parameters of the approximation function obtain the approximation function expression;
  • the selected spectrum set unit includes a spectrum set processing unit, and the spectrum set processing unit is used to perform smoothing and noise reduction processing on the spectrum set by using an algorithm.
  • the present invention has the following beneficial effects:
  • the self-adaptive surface absorption spectrum analysis method only needs to collect the sample surface reflectance spectrum once during the measurement process. It is not necessary to collect the reflected light of the reference object as a blank spectrum before the measurement, because the blank spectrum used in the calculation is reflected from the sample surface
  • the spectrum curve is extracted synchronously and in real time, so this method can eliminate the influence of the relative position change of the sample and the optical probe of the instrument, and adapt to the change of the distance between the sample and the optical probe within a certain range. Since the fitted blank spectrum curve contains the same substrate absorption band as the sample surface absorption spectrum, this method can eliminate the quantitative error caused by the superposition of substrate absorption.
  • this method adapts to the change of the distance between the sample and the optical probe and the difference of the matrix composition of different samples, and achieves a more accurate quantitative analysis of the narrow-band characteristic absorption peaks in the sample.
  • Figure 1 is a flowchart of the adaptive surface absorption spectroscopy analysis method of the present invention
  • Figure 2 is a specific flow chart of the adaptive surface absorption spectrum analysis method of the present invention.
  • Figure 3(a) is the three reflectance spectra of the surface layer of a gemstone detected by the traditional analysis method
  • Figure 3(b) is the three absorption spectra of the surface layer of a gemstone detected by the traditional analysis method
  • Figure 4 (a) is the three reflection spectra of the surface layer of a gemstone detected by the spectral analysis method of the present invention
  • Figure 4(b) is the three absorption spectra of the surface layer of a certain gemstone detected by the spectral analysis method of the present invention.
  • Figure 5 is a schematic diagram of the adaptive surface absorption spectrum analysis system of the present invention.
  • the adaptive surface absorption spectrum analysis method includes the following steps:
  • the reflectance spectrum data includes wavelength and light intensity signal values
  • a broad-band light source is used to illuminate the sample surface
  • a spectrometer collects the sample surface Reflective absorption spectrum to obtain a sample surface reflection spectrum including potential absorption peaks with wavelength as the abscissa and light intensity AD as the ordinate, where the light intensity AD value is the signal value output by the light intensity through AD conversion, hereinafter referred to as Light intensity AD value.
  • Illuminate the surface of the sample with a broad-band light source collect the reflection and absorption spectrum of the sample surface with a spectrometer, and obtain a sample surface reflection spectrum with wavelength as the abscissa and light intensity AD as the ordinate, that is, the surface reflection spectrum in step S21 As shown in the figure, the reflected light intensity of the sample surface can reach more than 5% of the full range of the spectrometer.
  • the distance between the sample surface and the optical probe of the spectrometer is not limited. Each collection can change the sample surface and the spectrometer optics. The distance between the probes, because the blank spectrum used in the calculation is synchronously and real-time extracted from the sample surface reflectance spectrum curve, which can eliminate the influence of the relative position change of the sample and the optical probe of the instrument, and it can adapt to the sample and the sample within a certain range. The change in distance between optical probes.
  • step S2 further includes the following steps:
  • S21 Select a spectrum set.
  • the potential absorption peak refers to the current target wavelength range for narrowband absorption peak analysis.
  • potential is mainly because there may or may not be a narrowband absorption peak in this wavelength band.
  • the model function is A function expression using wavelength as an independent variable to approximate the aforementioned discrete spectral data points, that is, multiple data points with discrete wavelength as the abscissa and light intensity AD as the ordinate, and the wavelength of each discrete point is substituted into the function expression Subtract the measured light intensity AD value corresponding to the wavelength from the formula, where the measured light intensity AD value is the new value after smoothing, the difference obtained is the residual at this point, and the functional expression of the difference is the residual Difference expression.
  • the model function used in the present invention is one of the commonly used Lagrange polynomials, Hermite polynomials, spline functions, Fourier series, etc., and is not limited to a specific function.
  • step S24 further steps are included:
  • Blank spectrum fitting According to the approximation function expression, bring all wavelengths including potential absorption peaks into the approximation function expression to calculate the corresponding light intensity signal value, and combine all the obtained light intensity signal values to obtain the blank spectrum Graph.
  • the blank spectrum fitting is to calculate the fitted light intensity of all wavelength points through the previously solved best approximation function expression, that is, the wavelength value is substituted into the expression to calculate the corresponding AD value, including the potential peak position That is to say, all wavelength points including the target wavelength band of the current narrowband absorption spectrum peak analysis are calculated, and a new spectrum curve is obtained by fitting. This new spectrum curve is used as a blank spectrum to participate in the subsequent calculation of the absorption spectrum.
  • Absorption spectrum analysis processing performing a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value extracted by using the blank spectrum and inverting the absorbance value of the wavelength point.
  • the measured light intensity AD value of each wavelength point is divided by the fitted blank light intensity AD value of the point, the quotient is logarithmically calculated and then the inverse is multiplied by minus 1 to obtain the absorbance value of the point,
  • the curve formed by wavelength as the abscissa and absorbance as the ordinate is the absorption spectrum.
  • the calculation method of the absorption spectrum in this embodiment is the same as the traditional analysis method, that is, both the measured light intensity and the blank light intensity are punished and logarithmic calculations.
  • the difference lies in the blank light intensity used in the present invention. It is synchronously extracted from the sample reflectance spectrum curve in real time, and its overall intensity changes synchronously with the reflected light from the sample surface, which can eliminate the interference of positional changes and matrix differences, and obtain a more stable and accurate narrow-band characteristic absorption spectrum, providing qualitative and quantitative Accuracy.
  • step S21 further steps are included:
  • the algorithm used for spectral set smoothing and noise reduction processing can be simple multi-point averaging or complex moving window smoothing. It is not limited to a specific algorithm.
  • the spectral set smoothing noise reduction processing reduces noise Interference is that the measured light intensity AD value is more accurate.
  • An electronic device comprising: a processor; a memory; and a program, wherein the program is stored in the memory and is configured to be executed by the processor, and the program includes a method for performing an adaptive surface absorption spectroscopy analysis method .
  • a computer-readable storage medium has a computer program stored thereon, and the computer program is executed by a processor by an adaptive surface absorption spectrum analysis method.
  • the adaptive surface absorption spectrum analysis system includes a spectrum data acquisition module, an absorption spectrum calculation module and a blank spectrum extraction module; among them,
  • the acquiring spectrum data module is used to acquire the reflection absorption spectrum of the surface of the sample to obtain reflection spectrum data; wherein the reflection absorption spectrum is a spectrum diagram with wavelength as the abscissa and light intensity signal value as the ordinate;
  • the blank spectrum extraction module is used to select part of the spectrum data in the reflectance spectrum data, and extract the blank spectrum detected once by performing data processing on the part of the spectrum data; wherein, the part of the spectrum data does not include potential absorption peak data ;
  • the blank spectrum extraction module includes a selection spectrum collection unit, a calculation unit, and a blank spectrum fitting unit; the selection spectrum collection unit is used to select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak in the reflection spectrum data, Forming a discrete combined spectrum set; the calculation unit is used to calculate the discrete combined spectrum set through the establishment function to obtain an approximation function expression;
  • the absorption spectrum calculation module is used to perform a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value obtained by using the blank spectrum and invert to obtain the absorbance value of the wavelength point.
  • the calculation unit includes a model residual sorting unit, a unit for establishing the optimal approximation condition equation group, and an approximation function solving unit;
  • the model residual arranging unit is used to establish an approximation function with the wavelength of a discrete point as the abscissa and the light intensity signal value as the ordinate, and the light intensity signal value obtained by the approximation function corresponds to the actual measured light intensity corresponding to the same discrete point
  • the function expression of the difference obtained from the signal value is recorded as the residual expression
  • the unit for establishing the best approximation condition equation group is used to perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the best approximation condition equation group;
  • the approximation function solving unit is used to calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and the coefficient parameters of the approximation function obtain the approximation function expression;
  • the selected spectrum set unit includes a spectrum set processing unit, and the spectrum set processing unit is used to perform smoothing and noise reduction processing on the spectrum set by using an algorithm.
  • Figure 3(a) and Figure 3(b) are the analysis results of the traditional analysis method.
  • the thicker smooth curve 10 in Figure 3(a) is the blank spectrum measured by using the PTFE whiteboard reference material.
  • the surface absorption spectrum analysis of the reflection spectrum of the gemstone in the second detection uses the gray spectrum as the blank spectrum. After triggering and logarithmic operation, three absorption spectrum curves as shown in Figure 3(b) are obtained.
  • the peak absorbances of the three narrow-band absorption peaks are 0.049, 0.107, and 0.204, respectively, and the peak absorbance values are very different; 3 absorptions
  • the peak-to-valley difference in absorbance of the narrow-band characteristic absorption peak in the spectrum is 0.050, 0.053, and 0.057, respectively. Even if the intensity of the narrow-band characteristic absorption peak is calculated according to the peak-to-valley difference, the average error of the three measurements is greater than 6.8%.
  • Figure 4 (a) and Figure 4 (b) are the analysis results of the analysis method of the present invention.
  • the three thicker smooth curves 201, 202 and 203 in Figure 4 (a) are synchronized in real time from three gem reflection spectra. Fit the extracted 3 blank spectra, use the reflectance spectra detected three times and the blank spectra extracted simultaneously in real time to calculate the absorbance, and obtain the 3 absorption spectrum curves shown in Figure 4(b), which can be seen from the figure
  • the three absorption spectrum curves are highly coincident, and their peak absorbances are 0.064, 0.063, and 0.061, respectively, and the average error of the three measurements is less than 2.7%.
  • this method can eliminate the interference of positional distance changes and matrix differences, and get more Stable and accurate narrow-band characteristic absorption spectrum improves qualitative and quantitative accuracy.
  • the method of the present invention only needs to collect the surface reflectance spectrum of the sample once during the measurement process, and it is not necessary to collect the reflected light of the reference object as the blank spectrum before the measurement, because the blank spectrum used in the calculation is from the surface reflectance curve of the sample Synchronous real-time extraction, so the method of the present invention can eliminate the influence of the relative position change of the sample and the optical probe of the instrument, and adapt to the change of the distance between the sample and the optical probe within a certain range.
  • the method of the present invention can eliminate the quantitative error caused by the superposition of substrate absorption.
  • the method of the present invention adapts to the change of the distance between the sample and the optical probe and the difference of the matrix composition of different samples, and realizes more accurate quantification of the narrow-band characteristic absorption peaks in the sample. analysis.

Abstract

A self-adaptive surface absorption spectrum analysis method, comprising the following steps: obtaining spectrum data (S1), extracting a blank spectrum (S2), and analyzing and processing an absorption spectrum (S3). The present invention relates to a self-adaptive surface absorption spectrum analysis system, an electronic device, and a storage medium. A sample surface reflection spectrum diagram only needs to be collected once in the measurement process, and reflected light of a reference object is not required to be collected before measurement as the blank spectrum; due to the fact that a fitting blank spectrum curve comprises a same substrate absorption band in the sample surface absorption spectrum, the analysis method can eliminate a quantitative error caused by absorption and superposition of a substrate. By synchronously extracting the blank spectrum in real time from a sample reflection spectrum, the analysis method can adapt to a change of the distance between a sample and an optical probe and a difference of different sample matrix components, thereby realizing more accurate quantitative analysis of a narrow-band characteristic absorption peak in the sample.

Description

自适应表面吸收光谱分析方法、系统、存储介质、设备Adaptive surface absorption spectrum analysis method, system, storage medium and equipment 技术领域Technical field
本发明涉及光谱分析领域,尤其涉及自适应表面吸收光谱分析方法。The present invention relates to the field of spectrum analysis, in particular to an adaptive surface absorption spectrum analysis method.
背景技术Background technique
表面吸收光谱检测是一项无损光谱分析技术,利用宽谱带光源照射样品表面、检测样品表层的反射光,通过表层物质对入射光的吸收情况来判断样品中某些特征分子基团的含量。由于无需破坏样品形态结构、检测速度快,其在工业无损检测领域应用较为广泛,近些年在医学体外无损诊断领域也展现出了良好的发展势头。Surface absorption spectroscopy detection is a non-destructive spectroscopic analysis technology that uses a wide-band light source to illuminate the surface of the sample, detect the reflected light from the surface of the sample, and determine the content of certain characteristic molecular groups in the sample by the absorption of incident light by the surface material. Because it does not need to destroy the sample morphology and structure, and the detection speed is fast, it is widely used in the field of industrial non-destructive testing. In recent years, it has also shown a good momentum of development in the field of medical in vitro non-destructive diagnosis.
传统的表面吸收光谱分析方法需要采用一个不含待测特征物质或特征分子基团的样品作为参照物,测量参照物的表面反射光谱作为空白光,然后让样品处于同一相对位置测量其表面反射光谱后与空白光相除转换为反射率或是吸光度曲线来进行定性和定量分析。其缺点是样品与仪器光学探头之间相对位置的变化易对测量所得的特征峰强度产生较大的影响,同时无法排除特征峰所在的波长位置前后存在的吸收带的干扰,使得检测结果容易出现误判,影响着检测的便捷性和准确性。The traditional surface absorption spectroscopy analysis method needs to use a sample that does not contain the characteristic substance or characteristic molecular group to be measured as a reference object, measure the surface reflectance spectrum of the reference object as a blank light, and then place the sample in the same relative position to measure its surface reflectance spectrum After dividing by the blank light, it is converted into a reflectance or absorbance curve for qualitative and quantitative analysis. The disadvantage is that the change of the relative position between the sample and the optical probe of the instrument will easily have a greater impact on the measured characteristic peak intensity, and the interference of the absorption band before and after the wavelength position of the characteristic peak cannot be ruled out, making the detection result easy to appear Misjudgment affects the convenience and accuracy of detection.
发明内容Summary of the invention
为了克服现有技术的不足,本发明提供自适应表面吸收光谱分析方法,该方法利用宽谱带光源光强谱带连续平缓、基底吸收带变化缓慢的特性,采集样品表层反射光谱后避开潜在窄带吸收峰位、采用潜在窄带吸收峰位前后 一定距离外的两段光强谱带拟合还原得到一条无窄带吸收的光强曲线作为空白光谱,然后用样品表层反射光谱曲线与空白光谱曲线相除取对数得到样品的表面吸收光谱曲线图。In order to overcome the shortcomings of the prior art, the present invention provides an adaptive surface absorption spectroscopy analysis method. The method uses the characteristics of a wide-band light source with continuous and gentle light intensity bands and slow changes in substrate absorption bands, and avoids potential problems after collecting sample surface reflection spectra. Narrowband absorption peak position, using two light intensity bands at a certain distance before and after the potential narrowband absorption peak position to fit and restore a light intensity curve without narrowband absorption as a blank spectrum, and then divide the sample surface reflectance spectrum curve and the blank spectrum curve Logarithm to obtain the surface absorption spectrum curve of the sample.
本发明提供自适应表面吸收光谱分析方法,包括如下步骤:The present invention provides an adaptive surface absorption spectrum analysis method, which includes the following steps:
S1、获取光谱数据,获取样品表层的反射光谱数据;其中,所述反射光谱数据包括波长与光强信号值;S1. Obtain spectrum data, and obtain reflectance spectrum data of the sample surface; wherein, the reflectance spectrum data includes wavelength and light intensity signal values;
S2、空白光谱提取,选取反射光谱数据中的部分光谱数据,通过对所述部分光谱数据进行数据处理,提取得到当前检测的空白光谱;其中,所述部分光谱数据不包括潜在吸收峰数据;S2. Blank spectrum extraction, selecting part of the spectrum data in the reflection spectrum data, and extracting the currently detected blank spectrum by performing data processing on the part of the spectrum data; wherein the part of the spectrum data does not include potential absorption peak data;
S3、吸收光谱分析处理,将所述波长点的实测光强信号值与利用空白光谱提取得到的空白光强信号值的商进行对数运算并取反得到所述波长点的吸光度值。S3. Absorption spectrum analysis processing, performing a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value extracted by using the blank spectrum and inverting the absorbance value of the wavelength point.
优选地,在步骤S2中还包括如下步骤:Preferably, step S2 further includes the following steps:
S21、选取光谱集,在所述反射光谱数据中位于潜在吸收峰前后波长间隔处各选取一段谱带,形成离散组合光谱集;S21. Select a spectrum set, and select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak in the reflection spectrum data to form a discrete combined spectrum set;
S22、模型残差整理,建立以离散点波长为横坐标、光强信号值为纵坐标的逼近函数,通过所述逼近函数得出的光强信号值与同一离散点对应的实测光强信号值的差值的函数表达式,记为残差表达式;S22. Organize the residuals of the model, establish an approximation function with the wavelength of the discrete point as the abscissa and the light intensity signal as the ordinate, and the light intensity signal value obtained by the approximation function is the actual measured light intensity signal value corresponding to the same discrete point The function expression of the difference of is recorded as the residual expression;
S23、建立最佳逼近条件方程组,根据最佳逼近评价方法对所述残差表达式进行运算转换得到最佳逼近条件方程组;S23. Establish a best approximation condition equation group, and perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the best approximation condition equation group;
S24、逼近函数求解,根据求解最佳逼近条件方程组计算出逼近函数的全部系数参数,由逼近函数的系数参数得出逼近函数表达式。S24. Solve the approximation function, calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and obtain the approximation function expression from the coefficient parameters of the approximation function.
优选地,在步骤S24之后还包括步骤:Preferably, the step S24 is further included:
S25、空白光谱拟合,根据逼近函数表达式将包括潜在吸收峰在内的全部波长带入逼近函数表达式计算得出相应的光强信号值,将得到的所有光强信号值组合得到空白光谱曲线图。S25. Blank spectrum fitting. According to the approximation function expression, bring all wavelengths including potential absorption peaks into the approximation function expression to calculate the corresponding light intensity signal value, and combine all the obtained light intensity signal values to obtain the blank spectrum Graph.
优选地,在步骤S21之后还包括步骤:Preferably, after step S21, the method further includes:
S211、光谱集处理,利用算法对光谱集进行平滑降噪处理。S211, spectrum set processing, using an algorithm to perform smoothing and noise reduction processing on the spectrum set.
优选地,所述实测光强信号值为经过平滑降噪处理后取得的值。Preferably, the measured light intensity signal value is a value obtained after smoothing and noise reduction processing.
优选地,在步骤S1中还包括,利用宽谱带光源照射样品表面、用光谱仪器获取样品表层的反射吸收光谱,样品表面的反射光的强度大于等于所述光谱仪器的满量程的5%。Preferably, step S1 further includes: irradiating the surface of the sample with a wide-band light source, and obtaining the reflection absorption spectrum of the sample surface with a spectrometer, and the intensity of the reflected light on the sample surface is greater than or equal to 5% of the full range of the spectrometer.
一种电子设备,包括:处理器;An electronic device, including: a processor;
存储器;以及程序,其中所述程序被存储在所述存储器中,并且被配置成由处理器执行,所述程序包括自适应表面吸收光谱分析方法。A memory; and a program, wherein the program is stored in the memory and is configured to be executed by a processor, the program including an adaptive surface absorption spectroscopy analysis method.
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行自适应表面吸收光谱分析方法。A computer-readable storage medium has a computer program stored thereon, and the computer program is executed by a processor by an adaptive surface absorption spectrum analysis method.
自适应表面吸收光谱分析系统,包括获取光谱数据模块、吸收光谱计算模块和空白光谱提取模块;其中,The self-adaptive surface absorption spectrum analysis system includes a spectrum data acquisition module, an absorption spectrum calculation module and a blank spectrum extraction module; among them,
所述获取光谱数据模块用于获取样品表层的反射吸收光谱,得到反射光谱数据;其中,所述反射吸收光谱为以波长为横坐标、光强信号值为纵坐标的光谱图;The acquiring spectrum data module is used to acquire the reflection absorption spectrum of the surface of the sample to obtain reflection spectrum data; wherein the reflection absorption spectrum is a spectrum diagram with wavelength as the abscissa and light intensity signal value as the ordinate;
所述空白光谱提取模块用于选取反射光谱数据中的部分光谱数据,通过对所述部分光谱数据进行数据处理,提取得到一次检测的空白光谱;其中,所述部分光谱数据不包括潜在吸收峰数据;The blank spectrum extraction module is used to select part of the spectrum data in the reflectance spectrum data, and extract the blank spectrum detected once by performing data processing on the part of the spectrum data; wherein the part of the spectrum data does not include potential absorption peak data ;
所述空白光谱提取模块包括选取光谱集单元、计算单元和空白光谱拟合单元;所述选取光谱集单元用于在所述反射光谱数据中位于潜在吸收峰前后波长间隔处各选取一段谱带,形成离散组合光谱集;所述计算单元用于将离散组合光谱集通过建立函数计算后得到逼近函数表达式;The blank spectrum extraction module includes a selection spectrum collection unit, a calculation unit, and a blank spectrum fitting unit; the selection spectrum collection unit is used to select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak in the reflection spectrum data, Forming a discrete combined spectrum set; the calculation unit is used to calculate the discrete combined spectrum set through the establishment function to obtain an approximation function expression;
所述吸收光谱计算模块用于将波长点的实测光强信号值与利用空白光谱得到的空白光强信号值的商进行对数运算并取反得到该波长点的吸光度值。The absorption spectrum calculation module is used to perform a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value obtained by using the blank spectrum and invert to obtain the absorbance value of the wavelength point.
优选地,所述计算单元包括模型残差整理单元、建立最佳逼近条件方程 组单元和逼近函数求解单元;其中,Preferably, the calculation unit includes a model residual sorting unit, a unit for establishing the optimal approximation condition equation group, and an approximation function solving unit; wherein,
所述模型残差整理单元用于建立以离散点波长为横坐标、光强信号值为纵坐标的逼近函数,通过所述逼近函数得出的光强信号值与同一离散点对应的实测光强信号值所得的差值的函数表达式,记为残差表达式;The model residual arranging unit is used to establish an approximation function with the wavelength of a discrete point as the abscissa and the light intensity signal value as the ordinate, and the light intensity signal value obtained by the approximation function corresponds to the actual measured light intensity corresponding to the same discrete point The function expression of the difference obtained from the signal value is recorded as the residual expression;
所述建立最佳逼近条件方程组单元用于根据最佳逼近评价方法对所述残差表达式进行运算转换得到最佳逼近条件方程组;The unit for establishing the best approximation condition equation group is used to perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the best approximation condition equation group;
所述近函数求解单元用于根据求解最佳逼近条件方程组计算出逼近函数的全部系数参数,逼近函数的系数参数得出逼近函数表达式;The approximation function solving unit is used to calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and the coefficient parameters of the approximation function obtain the approximation function expression;
所述选取光谱集单元包括光谱集处理单元,所述光谱集处理单元用于利用算法对光谱集进行平滑降噪处理。The selected spectrum set unit includes a spectrum set processing unit, and the spectrum set processing unit is used to perform smoothing and noise reduction processing on the spectrum set by using an algorithm.
相比现有技术,本发明的有益效果在于:Compared with the prior art, the present invention has the following beneficial effects:
自适应表面吸收光谱分析方法只需在测量的过程中采集一次样品表层反射光谱图即可,无需在测量之前采集参照物的反射光作为空白光谱,由于计算时所使用的空白光谱从样品表层反射光谱曲线中同步实时提取,故本方法可以消除样品与仪器光学探头相对位置变化的影响,在一定的范围内自适应样品与光学探头之间距离远近的变化。由于拟合空白光谱曲线中包含了与样品表面吸收光谱中相同基底吸收带,因而本方法可以消除基底吸收叠加导致的定量误差。通过从样品反射光谱中同步实时提取空白光谱,本方法自适应样品与光学探头之间距离的变化以及不同样品基质成分的差异,实现了对样品中窄带特征吸收峰更为精确的定量分析。The self-adaptive surface absorption spectrum analysis method only needs to collect the sample surface reflectance spectrum once during the measurement process. It is not necessary to collect the reflected light of the reference object as a blank spectrum before the measurement, because the blank spectrum used in the calculation is reflected from the sample surface The spectrum curve is extracted synchronously and in real time, so this method can eliminate the influence of the relative position change of the sample and the optical probe of the instrument, and adapt to the change of the distance between the sample and the optical probe within a certain range. Since the fitted blank spectrum curve contains the same substrate absorption band as the sample surface absorption spectrum, this method can eliminate the quantitative error caused by the superposition of substrate absorption. By synchronously extracting the blank spectrum from the sample reflection spectrum in real time, this method adapts to the change of the distance between the sample and the optical probe and the difference of the matrix composition of different samples, and achieves a more accurate quantitative analysis of the narrow-band characteristic absorption peaks in the sample.
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附图详细说明如后。本发明的具体实施方式由以下实施例及其附图详细给出。The above description is only an overview of the technical solution of the present invention. In order to understand the technical means of the present invention more clearly and implement it in accordance with the content of the description, the preferred embodiments of the present invention are described in detail below with the accompanying drawings. The specific implementation of the present invention is given in detail by the following examples and drawings.
附图说明Description of the drawings
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:The drawings described here are used to provide a further understanding of the present invention and constitute a part of this application. The exemplary embodiments and descriptions of the present invention are used to explain the present invention, and do not constitute an improper limitation of the present invention. In the attached picture:
图1为本发明的自适应表面吸收光谱分析方法的流程图;Figure 1 is a flowchart of the adaptive surface absorption spectroscopy analysis method of the present invention;
图2为本发明的自适应表面吸收光谱分析方法的具体流程图;Figure 2 is a specific flow chart of the adaptive surface absorption spectrum analysis method of the present invention;
图3(a)为传统分析方法某宝石表层3次检测的三条反射光谱图;Figure 3(a) is the three reflectance spectra of the surface layer of a gemstone detected by the traditional analysis method;
图3(b)为传统分析方法某宝石表层3次检测的三条吸收光谱图;Figure 3(b) is the three absorption spectra of the surface layer of a gemstone detected by the traditional analysis method;
图4(a)为利用本发明的光谱分析方法某宝石表层3次检测的三条反射光谱图;Figure 4 (a) is the three reflection spectra of the surface layer of a gemstone detected by the spectral analysis method of the present invention;
图4(b)为利用本发明的光谱分析方法某宝石表层3次检测的三条吸收光谱图;Figure 4(b) is the three absorption spectra of the surface layer of a certain gemstone detected by the spectral analysis method of the present invention;
图5为本发明的自适应表面吸收光谱分析系统示意图;Figure 5 is a schematic diagram of the adaptive surface absorption spectrum analysis system of the present invention;
附图标记:10、采用聚四氟乙烯材质白板参照物实测所得的空白光谱,201、第一次同步实时拟合提取的空白光谱,202、第二次同步实时拟合提取的空白光谱,203、第三次同步实时拟合提取的空白光谱。Reference signs: 10. Blank spectrum obtained by actual measurement using PTFE whiteboard reference material, 201. Blank spectrum extracted by the first synchronization real-time fitting, 202. Blank spectrum extracted by second synchronization real-time fitting, 203 , The third synchronization real-time fitting of the extracted blank spectrum.
具体实施方式detailed description
下面,结合附图以及具体实施方式,对本发明做进一步描述,需要说明的是,在不相冲突的前提下,以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例。Hereinafter, the present invention will be further described with reference to the drawings and specific implementations. It should be noted that, provided that there is no conflict, the following embodiments or technical features can be combined to form new embodiments. .
自适应表面吸收光谱分析方法,如图1、图2所示,包括如下步骤:The adaptive surface absorption spectrum analysis method, as shown in Figure 1 and Figure 2, includes the following steps:
S1、获取光谱数据,获取样品表层的反射光谱数据;其中,所述反射光谱数据包括波长与光强信号值;在一个实施例中,利用宽谱带光源照射样品表面、光谱仪器采集样品表层的反射吸收光谱,得到一幅包括潜在吸收峰的以波长为横坐标、光强AD值为纵坐标的样品表层反射光谱图,其中光强AD值为光强通过AD转换输出的信号值,下面简称光强AD值。用宽谱带光源照 射样品表面、用光谱仪器采集样品表层反射吸收光谱,得到一幅以波长为横坐标、以光强AD值为纵坐标的样品表层反射光谱图即步骤S21中的表面反射光谱图,样品表面的反射光强达到光谱仪满量程的5%以上即可。S1. Acquire spectral data, and acquire reflectance spectrum data of the sample surface; wherein, the reflectance spectrum data includes wavelength and light intensity signal values; in one embodiment, a broad-band light source is used to illuminate the sample surface, and a spectrometer collects the sample surface Reflective absorption spectrum to obtain a sample surface reflection spectrum including potential absorption peaks with wavelength as the abscissa and light intensity AD as the ordinate, where the light intensity AD value is the signal value output by the light intensity through AD conversion, hereinafter referred to as Light intensity AD value. Illuminate the surface of the sample with a broad-band light source, collect the reflection and absorption spectrum of the sample surface with a spectrometer, and obtain a sample surface reflection spectrum with wavelength as the abscissa and light intensity AD as the ordinate, that is, the surface reflection spectrum in step S21 As shown in the figure, the reflected light intensity of the sample surface can reach more than 5% of the full range of the spectrometer.
需要说明的是,一般地,采集样品表面反射光谱需要多次采集,在多次的采集过程中不限定样品表面与光谱仪光学探头之间的距离,每次的采集都可以改变样品表面与光谱仪光学探头之间的距离,因为计算时所使用的空白光谱是从样品表层反射光谱曲线中同步实时提取的,如此可以消除样品与仪器光学探头相对位置变换的影响,在一定的范围内自适应样品与光学探头之间距离远近的变化。It should be noted that in general, collecting the reflectance spectrum of the sample surface requires multiple collections. During the multiple collections, the distance between the sample surface and the optical probe of the spectrometer is not limited. Each collection can change the sample surface and the spectrometer optics. The distance between the probes, because the blank spectrum used in the calculation is synchronously and real-time extracted from the sample surface reflectance spectrum curve, which can eliminate the influence of the relative position change of the sample and the optical probe of the instrument, and it can adapt to the sample and the sample within a certain range. The change in distance between optical probes.
S2、空白光谱提取,选取反射光谱数据中的部分光谱数据,通过对所述部分光谱数据进行数据处理,提取得到当前检测的空白光谱;其中,所述部分光谱数据不包括潜在吸收峰数据;在一个实施例中,在步骤S2中还包括如下步骤:S2. Blank spectrum extraction, selecting part of the spectrum data in the reflection spectrum data, and extracting the currently detected blank spectrum by performing data processing on the part of the spectrum data; wherein, the part of the spectrum data does not include potential absorption peak data; In an embodiment, step S2 further includes the following steps:
S21、选取光谱集,在所述反射光谱数据中位于潜在吸收峰前后波长间隔处各选取一段谱带,形成离散组合光谱集;在本实施例中,通过选取在样品的表面反射光谱图中的潜在窄带吸收峰的前后一定距离外的两段光强谱带,选取的光谱集位于潜在吸收峰前后一定距离即是避开潜在吸收峰。S21. Select a spectrum set. In the reflectance spectrum data, select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak to form a discrete combined spectrum set; in this embodiment, by selecting in the surface reflectance spectrum of the sample For the two light intensity bands at a certain distance before and after the potential narrowband absorption peak, the selected spectrum set is located at a certain distance before and after the potential absorption peak to avoid the potential absorption peak.
需要说明的是,潜在吸收峰指的是当前拟进行窄带吸收谱峰分析的目标波长段,之所以称作潜在主要是因为该波段有可能存在窄带吸收峰,也有可能不存在窄带吸收峰。It should be noted that the potential absorption peak refers to the current target wavelength range for narrowband absorption peak analysis. The reason why it is called potential is mainly because there may or may not be a narrowband absorption peak in this wavelength band.
S22、模型残差整理,建立以离散点波长为横坐标、光强AD值为纵坐标的逼近函数,通过所述逼近函数得出的光强AD值与同一离散点对应的实测光强AD值所得的差值的函数表达式,记为残差表达式;在本实施例中,模型残差整理主要的工作是建立逼近模型函数、计算光谱集各离散点残差,所述的模型函数是一个以波长为自变量、用于逼近前述离散光谱数据点,即以离散波长为横坐标、光强AD值为纵坐标的多个数据点的函数表达式,将各离散点 的波长代入函数表达式再减去该波长对应的实测光强AD值,此处实测光强AD值为平滑处理后的新值,所得的差值即为该点的残差,差值的函数表达式即为残差表达式。本发明采用的模型函数为Lagrange多项式、Hermite多项式、样条函数、Fourier级数等常用中的一种,不局限于某一特定的函数。S22. Organize the residuals of the model, establish an approximation function with the wavelength of the discrete point as the abscissa and the light intensity AD as the ordinate, and the light intensity AD value obtained by the approximation function and the actual measured light intensity AD value corresponding to the same discrete point The function expression of the obtained difference is recorded as the residual expression; in this embodiment, the main work of the model residual sorting is to establish an approximation model function and calculate the residual error of each discrete point in the spectrum set. The model function is A function expression using wavelength as an independent variable to approximate the aforementioned discrete spectral data points, that is, multiple data points with discrete wavelength as the abscissa and light intensity AD as the ordinate, and the wavelength of each discrete point is substituted into the function expression Subtract the measured light intensity AD value corresponding to the wavelength from the formula, where the measured light intensity AD value is the new value after smoothing, the difference obtained is the residual at this point, and the functional expression of the difference is the residual Difference expression. The model function used in the present invention is one of the commonly used Lagrange polynomials, Hermite polynomials, spline functions, Fourier series, etc., and is not limited to a specific function.
S23、建立最佳逼近条件方程组,根据最佳逼近评价方法对所述残差表达式进行运算转换得到最佳逼近条件方程组;在本实施例中,建立最佳逼近条件方程组时先选择一种最佳逼近评价方法,然后根据该方法对残差表达式进行运算转换得到多个条件等式,即为最佳逼近条件方程组。所述的“最佳逼近评价方法”指的是通过残差衡量拟合曲线与原离散数据点的逼近好坏的判别手段,本实施例所采用的最佳逼近评价方法为残差绝对和即绝对值之和最小法或最小二乘法即残差平方和最小法等评价方法中的一种,不局限于某一特定的评价方法。S23. Establish a set of optimal approximation condition equations, and perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the optimal approximation condition equation set; in this embodiment, select first when establishing the optimal approximation condition equation set A best approximation evaluation method, and then based on this method, the residual expression is converted into multiple condition equations, which is the best approximation condition equation group. The "best approximation evaluation method" refers to a method of judging the approximation quality of the fitted curve and the original discrete data point by the residual error. The best approximation evaluation method adopted in this embodiment is the absolute sum of the residuals. One of the evaluation methods such as the minimum sum of absolute value or the least square method, that is, the minimum sum of squares residual method, is not limited to a specific evaluation method.
S24、逼近函数求解,根据求解最佳逼近条件方程组计算出逼近函数的全部系数参数,由逼近函数的系数参数得出逼近函数表达式;在本实施例中,求解逼近函数通过求解最佳逼近条件方程组来计算出逼近函数的全部系数参数,解出系数参数便可得到逼近函数表达式,本实施例中对方程组的求解既可以采用迭代法,也可以采用矩阵分解法或消元法,不局限于某一特定的解法。S24. Solve the approximation function, calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and obtain the approximation function expression from the coefficient parameters of the approximation function; in this embodiment, the solution of the approximation function is obtained by solving the best approximation Condition equations are used to calculate all the coefficient parameters of the approximation function, and the coefficient parameters can be solved to obtain the approximation function expression. In this embodiment, the equations can be solved by iterative method, matrix decomposition method or elimination method , Not limited to a specific solution.
在一个实施例中,在步骤S24之后还包括步骤:In an embodiment, after step S24, further steps are included:
S25、空白光谱拟合,根据逼近函数表达式将包括潜在吸收峰在内的全部波长带入逼近函数表达式计算得出相应的光强信号值,将得到的所有光强信号值组合得到空白光谱曲线图。在本实施例中,空白光谱拟合是通过前述求解出的最佳逼近函数表达式来计算全部波长点的拟合光强,即将波长值代入表达式计算出相应的AD值,包括潜在峰位即当前拟进行窄带吸收谱峰分析的目标波长段在内的全部波长点都进行计算,拟合得到一幅新的光谱曲线图,以此新光谱曲线图作为空白光谱参与后续吸收光谱的计算。S25. Blank spectrum fitting. According to the approximation function expression, bring all wavelengths including potential absorption peaks into the approximation function expression to calculate the corresponding light intensity signal value, and combine all the obtained light intensity signal values to obtain the blank spectrum Graph. In this embodiment, the blank spectrum fitting is to calculate the fitted light intensity of all wavelength points through the previously solved best approximation function expression, that is, the wavelength value is substituted into the expression to calculate the corresponding AD value, including the potential peak position That is to say, all wavelength points including the target wavelength band of the current narrowband absorption spectrum peak analysis are calculated, and a new spectrum curve is obtained by fitting. This new spectrum curve is used as a blank spectrum to participate in the subsequent calculation of the absorption spectrum.
S3、吸收光谱分析处理,将所述波长点的实测光强信号值与利用空白光谱提取得到的空白光强信号值的商进行对数运算并取反得到所述波长点的吸光度值。在一个实施例中,用各波长点的实测光强AD值除以该点的拟合空白光强AD值、对商进行对数运算然后取反即乘负1得到该点的吸光度值,以波长作为横坐标、吸光度作为纵坐标形成的连成的曲线即为吸收光谱图。S3. Absorption spectrum analysis processing, performing a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value extracted by using the blank spectrum and inverting the absorbance value of the wavelength point. In one embodiment, the measured light intensity AD value of each wavelength point is divided by the fitted blank light intensity AD value of the point, the quotient is logarithmically calculated and then the inverse is multiplied by minus 1 to obtain the absorbance value of the point, The curve formed by wavelength as the abscissa and absorbance as the ordinate is the absorption spectrum.
需要说明的是,本实施例中吸收光谱的计算方式与传统的分析方法相同,即都对测量光强和空白光强进行处罚和对数运算,不同之处在于本发明所采用的空白光强是从样品反射光谱曲线中同步实时提取,其整体强度随样品表层反射光同步变化,因而可以排除位置远近变化以及基质差异的干扰,得到更为稳定、准确的窄带特征吸收光谱,提供定性和定量的准确度。It should be noted that the calculation method of the absorption spectrum in this embodiment is the same as the traditional analysis method, that is, both the measured light intensity and the blank light intensity are punished and logarithmic calculations. The difference lies in the blank light intensity used in the present invention. It is synchronously extracted from the sample reflectance spectrum curve in real time, and its overall intensity changes synchronously with the reflected light from the sample surface, which can eliminate the interference of positional changes and matrix differences, and obtain a more stable and accurate narrow-band characteristic absorption spectrum, providing qualitative and quantitative Accuracy.
在一个实施例中,在步骤S21之后还包括步骤:In an embodiment, after step S21, further steps are included:
S211、光谱集处理,利用算法对光谱集进行平滑降噪处理。在本实施例中,光谱集平滑降噪处理采用算法既可以是简单的多点平均,也可以是复杂的移动窗口平滑,不局限于某一特定的算法,光谱集的平滑降噪处理降低噪声干扰,是的实测的光强AD值更加准确。S211, spectrum set processing, using an algorithm to perform smoothing and noise reduction processing on the spectrum set. In this embodiment, the algorithm used for spectral set smoothing and noise reduction processing can be simple multi-point averaging or complex moving window smoothing. It is not limited to a specific algorithm. The spectral set smoothing noise reduction processing reduces noise Interference is that the measured light intensity AD value is more accurate.
一种电子设备,包括:处理器;存储器;以及程序,其中所述程序被存储在所述存储器中,并且被配置成由处理器执行,所述程序包括用于执行自适应表面吸收光谱分析方法。一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行自适应表面吸收光谱分析方法。An electronic device, comprising: a processor; a memory; and a program, wherein the program is stored in the memory and is configured to be executed by the processor, and the program includes a method for performing an adaptive surface absorption spectroscopy analysis method . A computer-readable storage medium has a computer program stored thereon, and the computer program is executed by a processor by an adaptive surface absorption spectrum analysis method.
自适应表面吸收光谱分析系统,如图5所示,包括获取光谱数据模块、吸收光谱计算模块和空白光谱提取模块;其中,The adaptive surface absorption spectrum analysis system, as shown in Figure 5, includes a spectrum data acquisition module, an absorption spectrum calculation module and a blank spectrum extraction module; among them,
所述获取光谱数据模块用于获取样品表层的反射吸收光谱,得到反射光谱数据;其中,所述反射吸收光谱为以波长为横坐标、光强信号值为纵坐标的光谱图;The acquiring spectrum data module is used to acquire the reflection absorption spectrum of the surface of the sample to obtain reflection spectrum data; wherein the reflection absorption spectrum is a spectrum diagram with wavelength as the abscissa and light intensity signal value as the ordinate;
所述空白光谱提取模块用于选取反射光谱数据中的部分光谱数据,通过对所述部分光谱数据进行数据处理,提取得到一次检测的空白光谱;其中, 所述部分光谱数据不包括潜在吸收峰数据;The blank spectrum extraction module is used to select part of the spectrum data in the reflectance spectrum data, and extract the blank spectrum detected once by performing data processing on the part of the spectrum data; wherein, the part of the spectrum data does not include potential absorption peak data ;
所述空白光谱提取模块包括选取光谱集单元、计算单元和空白光谱拟合单元;所述选取光谱集单元用于在所述反射光谱数据中位于潜在吸收峰前后波长间隔处各选取一段谱带,形成离散组合光谱集;所述计算单元用于将离散组合光谱集通过建立函数计算后得到逼近函数表达式;The blank spectrum extraction module includes a selection spectrum collection unit, a calculation unit, and a blank spectrum fitting unit; the selection spectrum collection unit is used to select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak in the reflection spectrum data, Forming a discrete combined spectrum set; the calculation unit is used to calculate the discrete combined spectrum set through the establishment function to obtain an approximation function expression;
所述吸收光谱计算模块用于将波长点的实测光强信号值与利用空白光谱得到的空白光强信号值的商进行对数运算并取反得到该波长点的吸光度值。The absorption spectrum calculation module is used to perform a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value obtained by using the blank spectrum and invert to obtain the absorbance value of the wavelength point.
进一步地,所述计算单元包括模型残差整理单元、建立最佳逼近条件方程组单元和逼近函数求解单元;其中,Further, the calculation unit includes a model residual sorting unit, a unit for establishing the optimal approximation condition equation group, and an approximation function solving unit; wherein,
所述模型残差整理单元用于建立以离散点波长为横坐标、光强信号值为纵坐标的逼近函数,通过所述逼近函数得出的光强信号值与同一离散点对应的实测光强信号值所得的差值的函数表达式,记为残差表达式;The model residual arranging unit is used to establish an approximation function with the wavelength of a discrete point as the abscissa and the light intensity signal value as the ordinate, and the light intensity signal value obtained by the approximation function corresponds to the actual measured light intensity corresponding to the same discrete point The function expression of the difference obtained from the signal value is recorded as the residual expression;
所述建立最佳逼近条件方程组单元用于根据最佳逼近评价方法对所述残差表达式进行运算转换得到最佳逼近条件方程组;The unit for establishing the best approximation condition equation group is used to perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the best approximation condition equation group;
所述近函数求解单元用于根据求解最佳逼近条件方程组计算出逼近函数的全部系数参数,逼近函数的系数参数得出逼近函数表达式;The approximation function solving unit is used to calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and the coefficient parameters of the approximation function obtain the approximation function expression;
所述选取光谱集单元包括光谱集处理单元,所述光谱集处理单元用于利用算法对光谱集进行平滑降噪处理。The selected spectrum set unit includes a spectrum set processing unit, and the spectrum set processing unit is used to perform smoothing and noise reduction processing on the spectrum set by using an algorithm.
下面是运用本发明的自适应表面吸收光谱分析方法与传统分析方法的检测结果效果对比如图3(a)、图3(b)、图4(a)、图4(b)所示,用光谱仪器光学探头对准某一成分结构及其稳定的宝石表面采集其表层反射光谱,逐渐增大仪器光学探头与宝石表面的相对距离得到图3(a)和图4(a)曲线所示的光强逐渐减弱的3幅表层反射光谱图。The following is a comparison of the detection results of the adaptive surface absorption spectrum analysis method of the present invention and the traditional analysis method as shown in Figure 3 (a), Figure 3 (b), Figure 4 (a), and Figure 4 (b). The optical probe of the spectrometer is aligned with a certain component structure and its stable gem surface to collect the surface reflection spectrum, and the relative distance between the optical probe of the instrument and the gem surface is gradually increased to obtain the curve shown in Figure 3 (a) and Figure 4 (a) 3 surface reflectance spectra with decreasing light intensity.
图3(a)、图3(b)为传统分析方法分析效果图,其图3(a)中较粗的平滑曲线10为采用聚四氟乙烯材质白板参照物实测所得的空白光谱,对3次检测宝石反射光谱图的表面吸收光谱分析均以该灰色光谱作为空白光谱,经 触发和对数运算得到如图3(b)所示的三条吸收光谱曲线,从图中可以看出随着宝石样品与光学探头之间相对距离的增大即反射光谱强度降低所得的吸光度曲线逐步向上漂移,三次测量窄带吸收峰的峰值吸光度分别为0.049、0.107、和0.204,峰值吸光度数值相差巨大;3幅吸收光谱图中窄带特征吸收峰的吸光度峰谷差分别为0.050、0.053、和0.057,即便是按照峰谷差来计算窄带特征吸收峰的强度,三次测量均值误差大于6.8%。Figure 3(a) and Figure 3(b) are the analysis results of the traditional analysis method. The thicker smooth curve 10 in Figure 3(a) is the blank spectrum measured by using the PTFE whiteboard reference material. The surface absorption spectrum analysis of the reflection spectrum of the gemstone in the second detection uses the gray spectrum as the blank spectrum. After triggering and logarithmic operation, three absorption spectrum curves as shown in Figure 3(b) are obtained. The increase in the relative distance between the sample and the optical probe, that is, the decrease in the intensity of the reflection spectrum, the absorbance curve gradually drifts upward. The peak absorbances of the three narrow-band absorption peaks are 0.049, 0.107, and 0.204, respectively, and the peak absorbance values are very different; 3 absorptions The peak-to-valley difference in absorbance of the narrow-band characteristic absorption peak in the spectrum is 0.050, 0.053, and 0.057, respectively. Even if the intensity of the narrow-band characteristic absorption peak is calculated according to the peak-to-valley difference, the average error of the three measurements is greater than 6.8%.
图4(a)、图4(b)为本发明分析方法分析效果图,其图4(a)中3条较粗的平滑曲线201、202和203为从3幅宝石反射光谱图中同步实时拟合提取的3幅空白光谱图,用3次检测的反射光谱和各自同步实时提取的空白光谱求取吸光度,得到图4(b)所示的3条吸收光谱曲线,从图中可以看出3条吸收光谱曲线高度重合,其峰值吸光度分别为0.064、0.063、0.061,三次测量均值误差小于2.7%,和图4相比较可得本方法可以排除位置远近变化以及基质差异的干扰,得到更为稳定、准确的窄带特征吸收光谱,提高定性和定量的准确度。本发明中的方法只需在测量的过程中采集一次样品表层反射光谱图即可,无需在测量之前采集参照物的反射光作为空白光谱,由于计算时所使用的空白光谱从样品表层反射光谱曲线中同步实时提取,故本发明中的方法可以消除样品与仪器光学探头相对位置变化的影响,在一定的范围内自适应样品与光学探头之间距离远近的变化。由于拟合空白光谱曲线中包含了与样品表面吸收光谱中相同基底吸收带,因而本发明中的方法可以消除基底吸收叠加导致的定量误差。通过从样品反射光谱中同步实时提取空白光谱,本发明中的方法自适应样品与光学探头之间距离的变化以及不同样品基质成分的差异,实现了对样品中窄带特征吸收峰更为精确的定量分析。Figure 4 (a) and Figure 4 (b) are the analysis results of the analysis method of the present invention. The three thicker smooth curves 201, 202 and 203 in Figure 4 (a) are synchronized in real time from three gem reflection spectra. Fit the extracted 3 blank spectra, use the reflectance spectra detected three times and the blank spectra extracted simultaneously in real time to calculate the absorbance, and obtain the 3 absorption spectrum curves shown in Figure 4(b), which can be seen from the figure The three absorption spectrum curves are highly coincident, and their peak absorbances are 0.064, 0.063, and 0.061, respectively, and the average error of the three measurements is less than 2.7%. Compared with Figure 4, this method can eliminate the interference of positional distance changes and matrix differences, and get more Stable and accurate narrow-band characteristic absorption spectrum improves qualitative and quantitative accuracy. The method of the present invention only needs to collect the surface reflectance spectrum of the sample once during the measurement process, and it is not necessary to collect the reflected light of the reference object as the blank spectrum before the measurement, because the blank spectrum used in the calculation is from the surface reflectance curve of the sample Synchronous real-time extraction, so the method of the present invention can eliminate the influence of the relative position change of the sample and the optical probe of the instrument, and adapt to the change of the distance between the sample and the optical probe within a certain range. Since the fitted blank spectrum curve contains the same substrate absorption band as in the sample surface absorption spectrum, the method of the present invention can eliminate the quantitative error caused by the superposition of substrate absorption. By synchronously extracting the blank spectrum from the sample reflection spectrum in real time, the method of the present invention adapts to the change of the distance between the sample and the optical probe and the difference of the matrix composition of different samples, and realizes more accurate quantification of the narrow-band characteristic absorption peaks in the sample. analysis.
以上,仅为本发明的较佳实施例而已,并非对本发明作任何形式上的限制;凡本行业的普通技术人员均可按说明书附图所示和以上而顺畅地实施本发明;但是,凡熟悉本专业的技术人员在不脱离本发明技术方案范围内,利用以上所揭示的技术内容而做出的些许更动、修饰与演变的等同变化,均为本 发明的等效实施例;同时,凡依据本发明的实质技术对以上实施例所作的任何等同变化的更动、修饰与演变等,均仍属于本发明的技术方案的保护范围之内。The above are only the preferred embodiments of the present invention, and do not limit the present invention in any form; any person of ordinary skill in the industry can smoothly implement the present invention as shown in the drawings and above in the specification; however, wherever Those skilled in the art make use of the technical content disclosed above to make minor changes, modifications, and equivalent changes within the scope of the technical solution of the present invention, all of which are equivalent embodiments of the present invention; meanwhile, Any changes, modifications, and evolutions of any equivalent changes made to the above embodiments based on the essential technology of the present invention still fall within the protection scope of the technical solution of the present invention.

Claims (10)

  1. 自适应表面吸收光谱分析方法,其特征在于,包括如下步骤:The adaptive surface absorption spectrum analysis method is characterized in that it comprises the following steps:
    S1、获取光谱数据,获取样品表层的反射光谱数据;其中,所述反射光谱数据包括波长与光强信号值;S1. Obtain spectrum data, and obtain reflectance spectrum data of the sample surface; wherein, the reflectance spectrum data includes wavelength and light intensity signal values;
    S2、空白光谱提取,选取反射光谱数据中的部分光谱数据,通过对所述部分光谱数据进行数据处理,提取得到当前检测的空白光谱;其中,所述部分光谱数据不包括潜在吸收峰数据;S2. Blank spectrum extraction, selecting part of the spectrum data in the reflection spectrum data, and extracting the currently detected blank spectrum by performing data processing on the part of the spectrum data; wherein the part of the spectrum data does not include potential absorption peak data;
    S3、吸收光谱分析处理,将所述波长点的实测光强信号值与利用空白光谱提取得到的空白光强信号值的商进行对数运算并取反得到所述波长点的吸光度值。S3. Absorption spectrum analysis processing, performing a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value extracted by using the blank spectrum and inverting the absorbance value of the wavelength point.
  2. 如权利要求1所述的自适应表面吸收光谱分析方法,其特征在于,在步骤S2中还包括如下步骤:The adaptive surface absorption spectroscopy analysis method of claim 1, wherein step S2 further comprises the following steps:
    S21、选取光谱集,在所述反射光谱数据中位于潜在吸收峰前后波长间隔处各选取一段谱带,形成离散组合光谱集;S21. Select a spectrum set, and select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak in the reflection spectrum data to form a discrete combined spectrum set;
    S22、模型残差整理,建立以离散点波长为横坐标、光强信号值为纵坐标的逼近函数,通过所述逼近函数得出的光强信号值与同一离散点对应的实测光强信号值的差值的函数表达式,记为残差表达式;S22. Organize the residuals of the model, establish an approximation function with the wavelength of the discrete point as the abscissa and the light intensity signal as the ordinate, and the light intensity signal value obtained by the approximation function is the actual measured light intensity signal value corresponding to the same discrete point The function expression of the difference of is recorded as the residual expression;
    S23、建立最佳逼近条件方程组,根据最佳逼近评价方法对所述残差表达式进行运算转换得到最佳逼近条件方程组;S23. Establish a best approximation condition equation group, and perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the best approximation condition equation group;
    S24、逼近函数求解,根据求解最佳逼近条件方程组计算出逼近函数的全部系数参数,由逼近函数的系数参数得出逼近函数表达式。S24. Solve the approximation function, calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and obtain the approximation function expression from the coefficient parameters of the approximation function.
  3. 如权利要求2所述的自适应表面吸收光谱分析方法,其特征在于,在步骤S24之后还包括步骤:3. The adaptive surface absorption spectroscopy analysis method of claim 2, wherein after step S24, the method further comprises:
    S25、空白光谱拟合,根据逼近函数表达式将包括潜在吸收峰在内的全部波长带入逼近函数表达式计算得出相应的光强信号值,将得到的所有光强信 号值组合得到空白光谱曲线图。S25. Blank spectrum fitting. According to the approximation function expression, bring all wavelengths including potential absorption peaks into the approximation function expression to calculate the corresponding light intensity signal value, and combine all the obtained light intensity signal values to obtain the blank spectrum Graph.
  4. 如权利要求2所述的自适应表面吸收光谱分析方法,其特征在于,在步骤S21之后还包括步骤:The adaptive surface absorption spectroscopy analysis method according to claim 2, wherein after step S21, the method further comprises:
    S211、光谱集处理,利用算法对光谱集进行平滑降噪处理。S211, spectrum set processing, using an algorithm to perform smoothing and noise reduction processing on the spectrum set.
  5. 如权利要求1-4任一项所述的自适应表面吸收光谱分析方法,其特征在于,所述实测光强信号值为经过平滑降噪处理后取得的值。The adaptive surface absorption spectroscopy analysis method according to any one of claims 1 to 4, wherein the measured light intensity signal value is a value obtained after smoothing and noise reduction processing.
  6. 如权利要求1所述的自适应表面吸收光谱分析方法,其特征在于,在步骤S1中还包括,利用宽谱带光源照射样品表面、用光谱仪器获取样品表层的反射吸收光谱,样品表面的反射光的强度大于等于所述光谱仪器的满量程的5%。The self-adaptive surface absorption spectroscopy analysis method of claim 1, wherein step S1 further comprises: irradiating the surface of the sample with a broad-band light source, using a spectrometer to obtain the reflection absorption spectrum of the sample surface, and the reflection of the sample surface The intensity of the light is greater than or equal to 5% of the full scale of the spectrometer.
  7. 一种电子设备,其特征在于包括:处理器;An electronic device, characterized by comprising: a processor;
    存储器;以及程序,其中所述程序被存储在所述存储器中,并且被配置成由处理器执行,所述程序包括用于执行如权利要求1所述的方法。A memory; and a program, wherein the program is stored in the memory and is configured to be executed by a processor, and the program includes a method for performing the method according to claim 1.
  8. 一种计算机可读存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行如权利要求1所述的方法。A computer-readable storage medium having a computer program stored thereon, wherein the computer program is executed by a processor to execute the method according to claim 1.
  9. 自适应表面吸收光谱分析系统,其特征在于,包括获取光谱数据模块、吸收光谱计算模块和空白光谱提取模块;其中,The self-adaptive surface absorption spectrum analysis system is characterized by comprising a spectrum data acquisition module, an absorption spectrum calculation module and a blank spectrum extraction module; among them,
    所述获取光谱数据模块用于获取样品表层的反射吸收光谱,得到反射光谱数据;其中,所述反射吸收光谱为以波长为横坐标、光强信号值为纵坐标的光谱图;The acquiring spectrum data module is used to acquire the reflection absorption spectrum of the surface of the sample to obtain reflection spectrum data; wherein the reflection absorption spectrum is a spectrum diagram with wavelength as the abscissa and light intensity signal value as the ordinate;
    所述空白光谱提取模块用于选取反射光谱数据中的部分光谱数据,通过对所述部分光谱数据进行数据处理,提取得到一次检测的空白光谱;其中,所述部分光谱数据不包括潜在吸收峰数据;The blank spectrum extraction module is used to select part of the spectrum data in the reflectance spectrum data, and extract the blank spectrum detected once by performing data processing on the part of the spectrum data; wherein the part of the spectrum data does not include potential absorption peak data ;
    所述空白光谱提取模块包括选取光谱集单元、计算单元和空白光谱拟合单元;所述选取光谱集单元用于在所述反射光谱数据中位于潜在吸收峰前后波长间隔处各选取一段谱带,形成离散组合光谱集;所述计算单元用于将离 散组合光谱集通过建立函数计算后得到逼近函数表达式;The blank spectrum extraction module includes a selection spectrum collection unit, a calculation unit, and a blank spectrum fitting unit; the selection spectrum collection unit is used to select a section of spectrum bands at the wavelength intervals before and after the potential absorption peak in the reflection spectrum data, Forming a discrete combined spectrum set; the calculation unit is used to calculate the discrete combined spectrum set through the establishment function to obtain an approximation function expression;
    所述吸收光谱计算模块用于将波长点的实测光强信号值与利用空白光谱得到的空白光强信号值的商进行对数运算并取反得到该波长点的吸光度值。The absorption spectrum calculation module is used to perform a logarithmic operation on the quotient of the measured light intensity signal value of the wavelength point and the blank light intensity signal value obtained by using the blank spectrum and invert to obtain the absorbance value of the wavelength point.
  10. 如权利要求9所述的自适应表面吸收光谱分析系统,其特征在于,所述计算单元包括模型残差整理单元、建立最佳逼近条件方程组单元和逼近函数求解单元;其中,The adaptive surface absorption spectrum analysis system according to claim 9, wherein the calculation unit includes a model residual sorting unit, a unit for establishing an optimal approximation condition equation group, and an approximation function solving unit; wherein,
    所述模型残差整理单元用于建立以离散点波长为横坐标、光强信号值为纵坐标的逼近函数,通过所述逼近函数得出的光强信号值与同一离散点对应的实测光强信号值所得的差值的函数表达式,记为残差表达式;The model residual arranging unit is used to establish an approximation function with the wavelength of a discrete point as the abscissa and the light intensity signal value as the ordinate, and the light intensity signal value obtained by the approximation function corresponds to the actual measured light intensity corresponding to the same discrete point The function expression of the difference obtained from the signal value is recorded as the residual expression;
    所述建立最佳逼近条件方程组单元用于根据最佳逼近评价方法对所述残差表达式进行运算转换得到最佳逼近条件方程组;The unit for establishing the best approximation condition equation group is used to perform arithmetic conversion on the residual expression according to the best approximation evaluation method to obtain the best approximation condition equation group;
    所述近函数求解单元用于根据求解最佳逼近条件方程组计算出逼近函数的全部系数参数,逼近函数的系数参数得出逼近函数表达式;The approximation function solving unit is used to calculate all the coefficient parameters of the approximation function according to the optimal approximation condition equation set, and the coefficient parameters of the approximation function obtain the approximation function expression;
    所述选取光谱集单元包括光谱集处理单元,所述光谱集处理单元用于利用算法对光谱集进行平滑降噪处理。The selected spectrum set unit includes a spectrum set processing unit, and the spectrum set processing unit is used to perform smoothing and noise reduction processing on the spectrum set by using an algorithm.
PCT/CN2019/125394 2019-03-19 2019-12-13 Self-adaptive surface absorption spectrum analysis method and system, storage medium, and device WO2020186844A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201910206770.4A CN109991181B (en) 2019-03-19 2019-03-19 Adaptive surface absorption spectrum analysis method, system, storage medium and device
CN201910206770.4 2019-03-19

Publications (1)

Publication Number Publication Date
WO2020186844A1 true WO2020186844A1 (en) 2020-09-24

Family

ID=67129729

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/125394 WO2020186844A1 (en) 2019-03-19 2019-12-13 Self-adaptive surface absorption spectrum analysis method and system, storage medium, and device

Country Status (2)

Country Link
CN (1) CN109991181B (en)
WO (1) WO2020186844A1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109991181B (en) * 2019-03-19 2020-08-18 谱诉光电科技(苏州)有限公司 Adaptive surface absorption spectrum analysis method, system, storage medium and device
CN112782109B (en) * 2021-01-27 2023-07-04 中国计量科学研究院 Spectrum parameter acquisition method and device, electronic equipment and storage medium
CN113889196B (en) * 2021-12-07 2022-03-22 谱诉光电科技(苏州)有限公司 Spectrogram analysis method and device for component detection and electronic equipment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5910838A (en) * 1996-04-19 1999-06-08 Bodenseewerk Perkin-Elmer Gmbh Method of determining concentrations by means of continuous-spectrum atomic absorption spectroscopy
CN102938151A (en) * 2012-11-22 2013-02-20 中国人民解放军电子工程学院 Method for detecting anomaly of hyperspectral image
CN103076296A (en) * 2010-03-16 2013-05-01 新疆富科思生物技术发展有限公司 Background interference elimination method for UV visible absorption spectrum of optical fiber in situ drug dissolution/releasing rate tester
CN107167430A (en) * 2017-05-15 2017-09-15 浙江全世科技有限公司 A kind of method of emission spectrum background correction
CN109991181A (en) * 2019-03-19 2019-07-09 谱诉光电科技(苏州)有限公司 Adaptive surface absorption spectrum analyzing method, system, storage medium, equipment

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000071980A1 (en) * 1999-05-21 2000-11-30 Bookham Technology Plc Method and system for adaptive interpretation of spectrometric data combined with continual re-calibration
US8370243B1 (en) * 2008-09-01 2013-02-05 Prospercuity, LLC Financial planner and portfolio simulator
CN104237158B (en) * 2014-09-04 2016-07-27 浙江科技学院 A kind of Qualitative Analysis of Near Infrared Spectroscopy method with universality
CN104931518B (en) * 2015-06-09 2017-07-28 东南大学 A kind of method for X-ray fluorescence spectra background rejection
CN106770024B (en) * 2017-01-19 2019-04-30 东南大学 Based on the gas concentration measuring method for exempting from the direct absorption spectrum of baseline length scanning

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5910838A (en) * 1996-04-19 1999-06-08 Bodenseewerk Perkin-Elmer Gmbh Method of determining concentrations by means of continuous-spectrum atomic absorption spectroscopy
CN103076296A (en) * 2010-03-16 2013-05-01 新疆富科思生物技术发展有限公司 Background interference elimination method for UV visible absorption spectrum of optical fiber in situ drug dissolution/releasing rate tester
CN102938151A (en) * 2012-11-22 2013-02-20 中国人民解放军电子工程学院 Method for detecting anomaly of hyperspectral image
CN107167430A (en) * 2017-05-15 2017-09-15 浙江全世科技有限公司 A kind of method of emission spectrum background correction
CN109991181A (en) * 2019-03-19 2019-07-09 谱诉光电科技(苏州)有限公司 Adaptive surface absorption spectrum analyzing method, system, storage medium, equipment

Also Published As

Publication number Publication date
CN109991181A (en) 2019-07-09
CN109991181B (en) 2020-08-18

Similar Documents

Publication Publication Date Title
WO2020186844A1 (en) Self-adaptive surface absorption spectrum analysis method and system, storage medium, and device
CN101105446B (en) Differential optical absorption spectroscopy air quality detection system
CN108956505B (en) Method and device for detecting concentration of chlorophyll a in small water body based on Sentinel-2 image
CN109001207B (en) Method and system for detecting surface and internal defects of transparent material
CN105548070B (en) A kind of apple soluble solid near infrared detection position compensation method and system
CN110632021A (en) Spectrum detection method and system based on portable near-infrared spectrometer
CN105784672A (en) Drug detector standardization method based on dual-tree complex wavelet algorithm
CN111999258A (en) Spectral baseline correction-oriented weighting modeling local optimization method
US10548517B2 (en) Spectroscopic apparatus for biological material and spectroscopic method using the same
CN109001182A (en) The Raman spectrum non-destructive determination method of alcohol content in closed container
CN110987899B (en) Method for measuring content of total acid total esters in white spirit and identifying quality grade of white spirit
JP2007155630A (en) Thickness measuring method and thickness measuring apparatus of laminated thin film
JPH02271211A (en) Evaluating method of painting sharpness
CN116026780A (en) Method and system for online detection of coating moisture absorption rate based on series strategy wavelength selection
JPH1164217A (en) Component quantity detecting device for spectral analyzer
CN109459408A (en) A kind of Near-Infrared Quantitative Analysis method based on sparse regression LAR algorithm
CN112729555B (en) Method for synchronously diagnosing plasma temperature field by standard temperature method and relative spectral line method
CN110836878B (en) Convolution interpolation coupling Gaussian mixture model rapid three-dimensional fluorescence peak searching method
Harasim et al. Accuracy improvement of Fiber Bragg Grating peak wavelength demodulation using wavelet transform and various center wavelength detection algorithms
JP4233423B2 (en) Quantitative method and spectrum measuring apparatus
JPH02290537A (en) Method for estimating eating taste value by near infrared ray
CN111398212A (en) Method for establishing pepper detection model based on portable near-infrared spectrometer
CN117454098B (en) Dust concentration measuring method and system based on laser scattering
CN103592258A (en) Detection method of tea polyphenol contents in tea
CN112161708B (en) Spectrum preprocessing method for calculating color tristimulus values

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: 19920185

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19920185

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 19920185

Country of ref document: EP

Kind code of ref document: A1

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 01/04/2022)

122 Ep: pct application non-entry in european phase

Ref document number: 19920185

Country of ref document: EP

Kind code of ref document: A1