CN112730290B - Multi-optical path absorption spectrum spectrogram synthesis method - Google Patents
Multi-optical path absorption spectrum spectrogram synthesis method Download PDFInfo
- Publication number
- CN112730290B CN112730290B CN202110095882.4A CN202110095882A CN112730290B CN 112730290 B CN112730290 B CN 112730290B CN 202110095882 A CN202110095882 A CN 202110095882A CN 112730290 B CN112730290 B CN 112730290B
- Authority
- CN
- China
- Prior art keywords
- optical path
- absorption spectrum
- spectrum
- overlapping part
- absorption
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000862 absorption spectrum Methods 0.000 title claims abstract description 93
- 238000001308 synthesis method Methods 0.000 title claims abstract description 9
- 230000003287 optical effect Effects 0.000 claims abstract description 142
- 238000001514 detection method Methods 0.000 claims abstract description 74
- 238000001228 spectrum Methods 0.000 claims abstract description 67
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 15
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 15
- 238000002835 absorbance Methods 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 32
- 238000012937 correction Methods 0.000 claims description 18
- 230000002194 synthesizing effect Effects 0.000 claims description 17
- 230000003595 spectral effect Effects 0.000 claims description 15
- 238000012935 Averaging Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 230000031700 light absorption Effects 0.000 claims 1
- 238000010521 absorption reaction Methods 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 4
- 239000000243 solution Substances 0.000 description 8
- 238000012360 testing method Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004847 absorption spectroscopy Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000000975 dye Substances 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000012447 hatching Effects 0.000 description 1
- IWZKICVEHNUQTL-UHFFFAOYSA-M potassium hydrogen phthalate Chemical compound [K+].OC(=O)C1=CC=CC=C1C([O-])=O IWZKICVEHNUQTL-UHFFFAOYSA-M 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000012086 standard solution Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/314—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/314—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths
- G01N2021/3148—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry with comparison of measurements at specific and non-specific wavelengths using three or more wavelengths
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
Description
技术领域technical field
本发明是关于一种多光程吸收光谱谱图合成方法,涉及光谱检测技术领域。The invention relates to a method for synthesizing a multi-optical path absorption spectrum spectrum, and relates to the technical field of spectrum detection.
背景技术Background technique
进行吸光度(或吸收光谱)检测时,会通过比较多个光程下吸光度的数值,选择处于仪器最佳吸光度检测量程范围内的单个光程下的检测结果,代表样品的性质。对于单一波长的吸光度检测,这样的检测方式能够获取该波长最佳的吸光度检测结果,但是如果涉及到多个波长的吸光度(吸收光谱),这一方法则不能保证存在一个特定光程对不同波长吸光度的检测均能达到最佳。在医药、食品、染料、环境等领域,存在用吸收光谱对样本进行定性分析或化学计量学检测的需求。在这种需求中,吸收光谱形状的高保真度采集尤为重要,因而需要一种多光程吸收光谱谱图合成方法,使得在同一个吸收光谱图中,各个波长下的吸光度检测数据均来源于仪器的最佳检测范围。When performing absorbance (or absorption spectrum) detection, the detection result under a single optical path within the optimal absorbance detection range of the instrument will be selected by comparing the absorbance values under multiple optical paths, representing the properties of the sample. For the absorbance detection of a single wavelength, this detection method can obtain the best absorbance detection result of the wavelength, but if the absorbance (absorption spectrum) of multiple wavelengths is involved, this method cannot guarantee the existence of a specific optical path for the absorbance of different wavelengths. detection can be optimal. In the fields of medicine, food, dyes, environment, etc., there is a need for qualitative analysis or chemometric detection of samples by absorption spectroscopy. In this demand, the high-fidelity acquisition of the shape of the absorption spectrum is particularly important. Therefore, a multi-optical path absorption spectrum synthesis method is required, so that in the same absorption spectrum, the absorbance detection data at each wavelength are derived from The best detection range of the instrument.
现有技术中也有提出变光程检测光路或硬件的设计方法,但是都未涉及多光程吸收谱图合成方法。另外,由于检测设备进行吸收光谱检测时,存在超量程失真、基线漂移、随机误差等方面的干扰,在吸收光谱谱图合成时,需要克服各类检测误差,确保合成谱图的高保真性能。There are also design methods for variable optical path detection optical paths or hardware in the prior art, but none of them involve a multi-optical path absorption spectrum synthesis method. In addition, due to the interference of over-range distortion, baseline drift, random errors, etc., when the detection equipment performs absorption spectrum detection, when synthesizing absorption spectrum spectra, it is necessary to overcome various detection errors to ensure the high-fidelity performance of the synthesized spectrum.
发明内容SUMMARY OF THE INVENTION
针对上述问题,本发明的目的是提供一种能够确保光谱形状保真度的多光程吸收光谱谱图合成方法。In view of the above problems, the purpose of the present invention is to provide a method for synthesizing a multi-optical path absorption spectrum that can ensure the fidelity of the spectral shape.
为实现上述目的,本发明采取以下技术方案:一种多光程吸收光谱谱图合成方法,包括以下内容:To achieve the above object, the present invention adopts the following technical solutions: a method for synthesizing a multi-optical path absorption spectrum, comprising the following contents:
根据多光程吸收光谱检测设备性能选定待测光程组;Select the optical path group to be measured according to the performance of the multi-optical path absorption spectrum detection equipment;
根据选定光程组,通过多光程吸收光谱检测设备进行吸收光谱检测,并对采集的吸收光谱进行单位化处理;According to the selected optical path group, the absorption spectrum is detected by the multi-optical path absorption spectrum detection equipment, and the collected absorption spectrum is unitized;
对基于不同光程采集的吸收光谱进行光谱合成。Spectral synthesis is performed on absorption spectra collected based on different optical lengths.
所述的多光程吸收光谱谱图合成方法,进一步地,根据多光程吸收光谱检测设备性能选定待测光程组包括:The method for synthesizing the multi-optical path absorption spectrum spectrum, further, selecting the optical path group to be measured according to the performance of the multi-optical path absorption spectrum detection equipment includes:
确定多光程吸收光谱检测设备的最佳吸光度检测区间(a,b);Determine the optimal absorbance detection interval (a, b) of the multi-optical path absorption spectrum detection equipment;
将最佳吸光度检测区间的上限与下限的比值b/a作为跨度;Take the ratio b/a of the upper limit and lower limit of the optimal absorbance detection interval as the span;
在多光程吸收光谱检测设备可采集的最大光程和最小光程之间选取数个光程A1…An,使相邻光程比值Ai/Ai+1稍小于跨度b/a。Select several optical paths A1...An between the maximum optical path and the minimum optical path that can be collected by the multi-optical path absorption spectrum detection device, so that the adjacent optical path ratio Ai/Ai+1 is slightly smaller than the span b/a.
所述的多光程吸收光谱谱图合成方法,进一步地,对采集的吸收光谱进行单位化处理包括:In the method for synthesizing the multi-optical path absorption spectrum, further, unitizing the collected absorption spectrum includes:
将光程调节至A1…An并分别进行吸收光谱检测,采集各光程A1…An的吸收光谱,获取系列光程-吸收光谱图;Adjust the optical path to A1...An and perform absorption spectrum detection respectively, collect the absorption spectrum of each optical path A1...An, and obtain a series of optical path-absorption spectrum diagrams;
剔除实测吸收光谱与最佳吸光度检测区间完全不重叠的光程;Eliminate the optical path that does not overlap the measured absorption spectrum and the optimal absorbance detection range;
将有效的光谱除以其光程进行单位化处理。Unitize by dividing the effective spectrum by its optical length.
所述的多光程吸收光谱谱图合成方法,进一步地,根据比尔朗伯定律进行单位化处理,即用最佳吸光度检测区间的光谱除以其光程计算单位光程吸收光谱。In the method for synthesizing the multi-optical path absorption spectrum, further, unitization is performed according to Beer-Lambert's law, that is, the spectrum in the optimal absorbance detection interval is divided by its optical path to calculate the unit optical path absorption spectrum.
所述的多光程吸收光谱谱图合成方法,进一步地,对基于不同光程采集的吸收光谱进行光谱合成包括:In the method for synthesizing multi-optical path absorption spectra, further, performing spectral synthesis on absorption spectra collected based on different optical paths includes:
根据区间交集识别不同光程光谱最佳范围之间的重叠部分;Identify the overlap between the optimal ranges of different optical path spectra according to the interval intersection;
对相邻光程的重叠部分以光程即Ai和Ai+1为权重进行加权平均,获得均值作为重叠部分的最终光谱;The overlapping part of the adjacent optical path is weighted and averaged with the optical path Ai and Ai+1 as the weight, and the average value is obtained as the final spectrum of the overlapping part;
以重叠部分的最终光谱作为基准,计算每个光程的光谱重叠部分在平均化后被修正的平均值;Using the final spectrum of the overlapping part as a reference, calculate the average value of the spectral overlap part of each optical path after averaging;
如果存在最大光程或最小光程的延申部分,则随最大或最小光程的非重叠部分进行调整。If there is an extension of the maximum or minimum path length, adjust with the non-overlapping portion of the maximum or minimum path length.
所述的多光程吸收光谱谱图合成方法,进一步地,对于最大光程和最小光程,以重叠部分的修正值等量调整非重叠部分即可;中间的光程非重叠部分的修正参数为两端修正参数的差分值,从而使重叠部分与非重叠部分在完整光谱中具有较好的连续性。In the method for synthesizing the multi-optical path absorption spectrum, further, for the maximum optical path length and the minimum optical path length, the non-overlapping part can be adjusted by the same amount as the correction value of the overlapping part; the correction parameter of the non-overlapping part of the optical path in the middle The difference value of the parameter is corrected for both ends, so that the overlapping part and the non-overlapping part have better continuity in the complete spectrum.
所述的多光程吸收光谱谱图合成方法,进一步地,如果相邻光程重叠部分的误差偏离常数误差,可以对重叠部分加权平均前后的数值进行线性拟合,以平均前的数值为自变量X,平均后的数值为因变量Y,通过线性拟合获得斜率和截距,用于对非重叠部分的修正;最大光程和最小光程直接应用拟合获得的一次函数对非重叠部分进行修正;中间光程则将两端重叠部分的斜率和截距进行差分运算,获得各个波长吸光度对应的修正系数。In the method for synthesizing the multi-optical path absorption spectrum, further, if the error of the overlapping parts of adjacent optical paths deviates from the constant error, linear fitting can be performed on the values before and after the weighted average of the overlapping parts, and the value before the average is taken as the self-value. Variable X, the averaged value is the dependent variable Y, the slope and intercept are obtained by linear fitting, which are used to correct the non-overlapping part; the maximum optical path and the minimum optical path directly apply the linear function obtained by fitting to the non-overlapping part Correction is performed; for the middle optical path, the slope and intercept of the overlapping parts at both ends are differentially calculated to obtain the correction coefficient corresponding to the absorbance of each wavelength.
本发明由于采取以上技术方案,其具有以下优点:The present invention has the following advantages due to taking the above technical solutions:
1、本发明为同一样本不同光程的吸收光谱合并提供了完整的解决方案,合成后的谱图中可以同时存在吸收系数差异较大的吸收峰,提高了光谱形状的保真度;1. The present invention provides a complete solution for merging absorption spectra of different optical paths of the same sample, and absorption peaks with large differences in absorption coefficients can exist simultaneously in the synthesized spectrum, which improves the fidelity of the spectral shape;
2、本发明对于可进行变光程吸收光谱检测的仪器设备,提供了根据设备性能进行光程选择的准则和方法;2. The present invention provides a criterion and method for selecting an optical path according to the performance of the equipment for the instrument and equipment that can perform variable optical path absorption spectrum detection;
3、本发明通过多光程的合成,相比于单光程吸收光谱扩展了检测量程,尤其对于单光程下量程范围小但可适配光程范围大的设备,扩展量程的效果显著;3. The present invention expands the detection range compared with the single-light-path absorption spectrum through the synthesis of multiple optical paths, especially for equipment with a small range under a single optical path but can be adapted to a large optical path range, the effect of extending the range is remarkable;
4、在进行光谱合成的同时,本发明对于随机误差和基线漂移等具有误差修正的效果;4. While performing spectrum synthesis, the present invention has the effect of error correction for random errors and baseline drift;
综上,本发明对于具有自动变光程吸收光谱检测设备的开发和自动化应用方案提供了支持,可以广泛应用于样品光谱检测中。In conclusion, the present invention provides support for the development and automatic application scheme of an automatic variable optical path absorption spectrum detection device, and can be widely used in sample spectrum detection.
附图说明Description of drawings
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。在整个附图中,用相同的附图标记表示相同的部件。在附图中:Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are for the purpose of illustrating preferred embodiments only and are not to be considered limiting of the invention. The same reference numerals are used to refer to the same parts throughout the drawings. In the attached image:
图1为本发明实施例多光程吸收光谱谱图合成方法流程图;1 is a flowchart of a method for synthesizing a multi-optical path absorption spectrum according to an embodiment of the present invention;
图2为本发明实施例不同光程吸收光谱拼接示意图;2 is a schematic diagram of splicing absorption spectra of different optical paths according to an embodiment of the present invention;
图3为本发明实施例合成前后光谱图,其中(a)为合成前的各光程单位化光谱图,(b)为合成后的光谱图。FIG. 3 is a spectrogram before and after synthesis according to the embodiment of the present invention, wherein (a) is a spectrogram of each optical path unit before synthesis, and (b) is a spectrogram after synthesis.
具体实施方式Detailed ways
下面将参照附图更详细地描述本发明的示例性实施方式。虽然附图中显示了本发明的示例性实施方式,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了能够更透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited by the embodiments set forth herein. Rather, these embodiments are provided so that the present invention will be more thoroughly understood, and will fully convey the scope of the present invention to those skilled in the art.
应理解的是,文中使用的术语仅出于描述特定示例实施方式的目的,而无意于进行限制。除非上下文另外明确地指出,否则如文中使用的单数形式“一”、“一个”以及“所述”也可以表示包括复数形式。术语“包括”、“包含”、“含有”以及“具有”是包含性的,并且因此指明所陈述的特征、步骤、操作、元件和/或部件的存在,但并不排除存在或者添加一个或多个其它特征、步骤、操作、元件、部件、和/或它们的组合。文中描述的方法步骤、过程、以及操作不解释为必须要求它们以所描述或说明的特定顺序执行,除非明确指出执行顺序。还应当理解,可以使用另外或者替代的步骤。It is to be understood that the terminology used herein is for the purpose of describing particular example embodiments only and is not intended to be limiting. As used herein, the singular forms "a," "an," and "the" can also be intended to include the plural forms unless the context clearly dictates otherwise. The terms "comprising", "comprising", "containing" and "having" are inclusive and thus indicate the presence of stated features, steps, operations, elements and/or components, but do not preclude the presence or addition of one or Various other features, steps, operations, elements, components, and/or combinations thereof. Method steps, procedures, and operations described herein are not to be construed as requiring that they be performed in the particular order described or illustrated, unless an order of performance is explicitly indicated. It should also be understood that additional or alternative steps may be used.
为了便于描述,可以在文中使用空间相对关系术语来描述如图中示出的一个元件或者特征相对于另一元件或者特征的关系,这些相对关系术语例如为“内部”、“外部”、“内侧”、“外侧”、“下面”、“上面”等。这种空间相对关系术语意于包括除图中描绘的方位之外的在使用或者操作中装置的不同方位。For ease of description, spatially relative terms may be used herein to describe the relationship of one element or feature to another element or feature as shown in the figures, such as "inner", "outer", "inner" ", "outside", "below", "above", etc. This spatially relative term is intended to include different orientations of the device in use or operation other than the orientation depicted in the figures.
本发明提供一种多光程吸收光谱谱图合成方法,包括根据多光程吸收光谱检测设备性能选定待测光程组;根据选定光程组,通过多光程吸收光谱检测设备进行吸收光谱检测,并对采集的吸收光谱进行单位化处理;对基于不同光程采集的吸收光谱进行光谱合成。本发明为同一样本不同光程的吸收光谱合并提供了完整的解决方案,合成后的谱图中可以同时存在吸收系数差异较大的吸收峰,提高了光谱形状的保真度。The invention provides a method for synthesizing a multi-optical path absorption spectrum spectrum. Spectral detection and unitization of the collected absorption spectra; spectral synthesis of absorption spectra collected based on different optical path lengths. The invention provides a complete solution for merging absorption spectra of different optical paths of the same sample, and absorption peaks with large differences in absorption coefficients can exist simultaneously in the synthesized spectrum, thereby improving the fidelity of the spectral shape.
比尔朗伯定律是进行多光程光谱合并的基本规律,各吸光度与光程成正比,多光程吸收光谱的合成方法需要以比尔朗伯定律为基础,通过流程设计和算法处理,获取高保真度的合成谱图。Beer-Lambert's law is the basic rule for combining multi-optical path spectra. Each absorbance is proportional to the optical path. The synthesis method of multi-path absorption spectra needs to be based on Beer-Lambert's law. Through process design and algorithm processing, high fidelity can be obtained. Synthetic spectrum of degrees.
如图1所示,本发明实施例提供的多光程吸收光谱谱图合成方法,包括以下内容:As shown in FIG. 1 , the method for synthesizing a multi-optical path absorption spectrum provided by an embodiment of the present invention includes the following contents:
S1、对于能够进行多光程吸收光谱检测设备,根据检测设备性能选定待测光程组,具体包括:S1. For multi-optical path absorption spectrum detection equipment, select the optical path group to be measured according to the performance of the detection equipment, specifically including:
S11、通过查阅说明或实际测试确定多光程吸收光谱检测设备的最佳吸光度检测区间。S11. Determine the optimal absorbance detection range of the multi-optical path absorption spectrum detection device by referring to the description or actual testing.
具体地,根据多光程吸收光谱检测设备整体或其主要元器件的原始光强、动态范围、灵敏度、线性范围等参数,获取或估算最佳的吸光度检测区间范围。其中,“获取”即光谱检测设备整体说明书有直接的记录,“估算”是利用元器件的参数对整体的性能进行计算,上述过程均为现有技术,例如可以通过3倍标准偏差进行最佳检测区间的估算,以此为例,不限于此;或者也可利用标准滤光片等吸收光谱已知的器件,对多光程吸收光谱检测设备进行测试,得到对于所有波长均适用的最优吸光度检测区间,定义最佳的吸光度检测下限为a,上限为b,最佳区间为(a,b)。Specifically, the optimal absorbance detection interval range is obtained or estimated according to the original light intensity, dynamic range, sensitivity, linear range and other parameters of the multi-optical path absorption spectrum detection device as a whole or its main components. Among them, "acquisition" means that the overall specification of the spectral detection equipment has a direct record, and "estimation" is to use the parameters of the components to calculate the overall performance. The above processes are all existing technologies. The estimation of the detection interval, for example, is not limited to this; alternatively, a device with a known absorption spectrum, such as a standard filter, can be used to test the multi-optical path absorption spectrum detection equipment to obtain the optimal suitable for all wavelengths. Absorbance detection interval, define the optimal absorbance detection lower limit as a, upper limit as b, and the optimal interval as (a, b).
S12、将最优吸光度检测区间的上限与下限的比值b/a作为跨度。S12, taking the ratio b/a of the upper limit to the lower limit of the optimal absorbance detection interval as the span.
S13、在多光程吸收光谱检测设备可采集的最大光程和最小光程之间选取数个光程A1…An(由大到小),使相邻光程比值Ai/Ai+1稍小于跨度b/a,从而确保相邻两个光程的检测光谱具有部分重叠的最优检测区间,具体取值可以根据实际需要进行选择,在此不做限定。S13. Select several optical paths A1...An (from large to small) between the maximum optical path and the minimum optical path that can be collected by the multi-optical path absorption spectrum detection device, so that the adjacent optical path ratio Ai/Ai+1 is slightly smaller than The span is b/a, so as to ensure that the detection spectra of two adjacent optical paths have a partially overlapping optimal detection range, and the specific value can be selected according to actual needs, which is not limited here.
S2、采集吸收光谱并进行单位化处理,具体包括:S2. Collect absorption spectra and perform unitization processing, including:
S21、将光程调节至A1…An并分别进行吸收光谱检测,采集各光程A1…An的吸收光谱,获取系列光程-吸收光谱图。S21. Adjust the optical path to A1...An and perform absorption spectrum detection respectively, collect the absorption spectrum of each optical path A1...An, and obtain a series of optical path-absorption spectrum diagrams.
S22、剔除实测吸收光谱与最佳吸光度检测区间完全不重叠的光程(只会出现在光程组A1…An的最大端或最小端)。S22. Eliminate the optical path whose measured absorption spectrum and the optimal absorbance detection range do not overlap completely (only appear at the maximum end or the minimum end of the optical path group A1...An).
S23、将有效的光谱(即经过S22剔除后剩余的光程下测得的光谱)除以其光程进行单位化处理。S23: Divide the effective spectrum (that is, the spectrum measured under the remaining optical path after the elimination in S22) by its optical path, and perform unitization processing.
具体地,根据比尔朗伯定律进行单位化处理,即用最佳吸光度检测区间的光谱除以其光程计算单位光程吸收光谱。经过此步骤,获得的光谱在吸光度数轴上如图2所示,光程越大,对应的最佳吸光度检测区间经单位化后的数值越小。由于相邻两个光程的比值小于跨度,且吸收谱图为连续的曲线,因而任意相邻两个光程获得的吸收谱图均存在最佳范围的区间交集的2-7重叠部分,即相邻两个光程最佳检测区间的交集,该波长区间内,相邻光程测得的光谱均具有较高的可靠度。Specifically, the unitization process is performed according to Beer-Lambert's law, that is, the absorption spectrum per unit optical path is calculated by dividing the spectrum in the optimal absorbance detection interval by its optical path. After this step, the obtained spectrum is shown in Figure 2 on the absorbance number axis. The larger the optical path length, the smaller the unitized value of the corresponding optimal absorbance detection interval. Since the ratio of two adjacent optical paths is smaller than the span, and the absorption spectrum is a continuous curve, the absorption spectrum obtained by any two adjacent optical paths has an overlap of 2-7 of the interval intersection of the optimal range, that is, The intersection of the optimal detection interval of two adjacent optical paths, in this wavelength interval, the spectra measured by the adjacent optical paths have high reliability.
S3、由于在单位化时,以光程作为分母,将设备检测中与光程无关的随机误差、基线漂移等均进行了缩放,因而光程越大则该类型误差的引入越小,基于此进行光谱合成的具体过程包括:S3. Since the optical path is used as the denominator during unitization, the random errors and baseline drift that are not related to the optical path in the device detection are scaled. Therefore, the larger the optical path, the smaller the introduction of this type of error. Based on this The specific process of spectral synthesis includes:
S31、根据区间交集识别不同光程光谱最佳范围之间的2-7重叠部分。S31. Identify 2-7 overlapping parts between optimal ranges of different optical path spectra according to the interval intersection.
S32、对相邻光程的重叠部分以光程即Ai和Ai+1为权重进行加权平均,获得均值作为重叠部分的最终光谱,重叠部分仅包含相邻两个光程的区间对应的波长交集,该部分波长的最终光谱由加权平均的方式获得。由于以光程作为权重进行加权平均,大光程下检测的吸收光谱对重叠部分的最终光谱具有的影响,使最终光谱受大光程检测结果影响较大。S32. Perform a weighted average on the overlapping parts of the adjacent optical paths with the optical paths Ai and Ai+1 as the weights to obtain the average value as the final spectrum of the overlapping part, and the overlapping part only includes the wavelength intersection corresponding to the interval of the two adjacent optical paths , the final spectrum of this part of the wavelength is obtained by means of weighted average. Since the optical path is used as the weight to perform the weighted average, the absorption spectrum detected under the large optical path has an influence on the final spectrum of the overlapping part, so that the final spectrum is greatly affected by the detection result of the large optical path.
S33、以重叠部分的最终光谱作为基准,计算每个光程的光谱重叠部分在平均化后被修正的平均值。S33 , using the final spectrum of the overlapping portion as a reference, calculate the average value of the spectral overlap portion of each optical path corrected after averaging.
对于最大光程和最小光程,以重叠部分的修正值等量调整非重叠部分即可;For the maximum optical path length and the minimum optical path length, the non-overlapping part can be adjusted by the same amount as the correction value of the overlapping part;
中间的光程非重叠部分的修正参数为两端修正参数的差分值,从而使重叠部分与非重叠部分在完整光谱中具有较好的连续性。在此过程中,同时实现了对不同光程检测时出现的常数基线漂移的均值修正。The correction parameter of the non-overlapping part of the optical path in the middle is the difference value of the correction parameters at both ends, so that the overlapping part and the non-overlapping part have better continuity in the complete spectrum. In this process, the mean correction of the constant baseline drift that occurs in the detection of different optical lengths is achieved at the same time.
进一步地,如果相邻光程重叠部分的误差偏离常数误差(重叠部分是一段波长区间,其中包含多个波长的吸光度值,如果该区间内,最终光谱和某个光程的光谱的偏差在基本一致,那么就是常数误差,如果区间内从左到右偏差逐步增大或缩小,则是非常数误差),可以对重叠部分加权平均前后的数值进行线性拟合。以平均前的数值为自变量X,平均后的数值为因变量Y,通过线性拟合获得斜率和截距,用于对非重叠部分(重叠部分的最终光谱是加权平均而得,其他部分的数值经过修正衔接到重叠部分上)的修正。最大光程和最小光程直接应用拟合获得的一次函数对非重叠部分进行修正;中间光程则将两端重叠部分的斜率和截距进行差分运算,获得各个波长吸光度对应的修正系数。该方法可以对设备检测基线可能存在的线性漂移进行平均化处理。Further, if the error of the overlapping part of the adjacent optical paths deviates from the constant error (the overlapping part is a wavelength interval, which contains the absorbance values of multiple wavelengths, if within this interval, the deviation between the final spectrum and the spectrum of a certain optical path is basically Consistent, then it is a constant error, if the deviation from left to right in the interval gradually increases or decreases, it is a non-constant error), and the values before and after the weighted average of the overlapping parts can be linearly fitted. Taking the value before averaging as the independent variable X, and the value after averaging as the dependent variable Y, the slope and intercept are obtained by linear fitting, which are used for the non-overlapping part (the final spectrum of the overlapping part is weighted average, and the other parts are The value is corrected to connect to the overlapped part). The maximum optical path length and the minimum optical path directly apply the linear function obtained by fitting to correct the non-overlapping part; for the middle optical path, the slope and intercept of the overlapping parts at both ends are differentially calculated to obtain the correction coefficient corresponding to the absorbance of each wavelength. This method can average the possible linear drift of the device detection baseline.
S34、如果存在最大光程或最小光程的延申部分(吸光度小于最大光程的最佳范围或大于最小光程最佳范围的部分,这些部分不处于任何光程的最佳范围)即最大光程延申范围2-1和2-6最小光程延申范围,则随最大或最小光程的非重叠部分进行调整。S34. If there is an extension of the maximum optical path or the minimum optical path (the absorbance is less than the optimal range of the maximum optical path or greater than the optimal range of the minimum optical path, these parts are not in the optimal range of any optical path), the maximum Optical path extension ranges 2-1 and 2-6 The minimum optical path extension range is adjusted with the non-overlapping portion of the maximum or minimum optical path.
下面通过具体实施例详细说明本发明提供的多光程吸收光谱谱图合成方法的应用。The application of the method for synthesizing the multi-optical path absorption spectrum provided by the present invention will be described in detail below through specific embodiments.
将本实施例的多光程吸收光谱谱图合成方法应用于吸收光谱检测设备的变光程光谱合成,本实施例采用COD标准溶液邻苯二甲酸氢钾作为检测试剂。The method for synthesizing a multi-optical path absorption spectrum spectrum in this embodiment is applied to the variable optical path spectrum synthesis of an absorption spectrum detection device, and this embodiment uses a COD standard solution potassium hydrogen phthalate as a detection reagent.
根据步骤S1,对该吸收光谱检测设备进行光强动态范围测试,其最小检测波长处的最佳吸光度检测范围为(0.07,0.82),其他波长在的最佳吸光度范围均覆盖此区间,因而确定该设备的吸光度最佳范围为(0.07,0.82),计算其跨度为11.7,故可以选择10定为相邻两光程的比值。在该设备的可调光程范围内,可以取3个光程值A1,A2,A3分别为100mm,10mm和1mm。According to step S1, the light intensity dynamic range test is performed on the absorption spectrum detection device. The optimal absorbance detection range at the minimum detection wavelength is (0.07, 0.82), and the optimal absorbance range of other wavelengths covers this range, so it is determined The optimal range of absorbance of the device is (0.07, 0.82), and its span is calculated to be 11.7, so 10 can be selected as the ratio of two adjacent optical paths. Within the adjustable optical path range of the device, 3 optical path values A1, A2 and A3 can be taken as 100mm, 10mm and 1mm respectively.
根据步骤S2,分别检测这3个光程下的吸收光谱。According to step S2, the absorption spectra under the three optical paths are detected respectively.
为方便展示,将200nm到300nm以5mm为间隔记录如表1,各光程的吸收光谱均存在最佳范围内的检测值。各光谱除以其光程得到的单位光程(以m为单位)光谱如图3,并在表1中按照图2用阴影标出各光谱的最佳检测范围。For convenience of presentation, 200nm to 300nm are recorded at 5mm intervals as shown in Table 1, and the absorption spectra of each optical path have detection values within the optimal range. The unit optical path (in m) spectrum obtained by dividing each spectrum by its optical path is shown in Figure 3, and the optimal detection range of each spectrum is marked with hatching in Table 1 according to Figure 2.
表1不同光程的原始光谱、单位化光谱以及合成光谱表Table 1 Original spectrum, normalized spectrum and composite spectrum of different optical paths
根据步骤S3,首先将100mm光程下的光谱与10mm下的光谱合并。According to step S3, the spectrum at 100mm optical path is first combined with the spectrum at 10mm.
对255nm附近的两种光程重叠的最佳吸光度范围进行加权平均,根据步骤S32,100mm光程的权重为100,10mm光程权重为10,可得该区域的合成后的最终光谱,255nm处合成后吸光度为6.28,即(6.10*100+8.14*10)/11。根据步骤S33,须将100mm光程单位化后的数值上调0.18,即6.28-6.10=0.18,10mm光程单位化后的光谱数值下调1.86,即6.28-8.14=-1.86。对10mm光谱与1mm光谱的合并同样采用步骤S32进行2-7重叠区域的加权平均,获得210,215,220nm的最终吸光度分别为87.70,55.18和40.68。最终光谱重叠部分和1mm光谱188.32-87.70=100.62、107.84-55.18=52.66、71.83-40.68=31.15差值不是常数,而是和吸光度大小有关,存在显著的单向减小的趋势,因此发现该段可能存在线性基线误差,因而按照优选方法进行线性拟合。重叠部分均值对于10mm光谱的线性修正函数为Y=1.17X–3.28,对于1mm光谱的线性修正函数为Y=0.41X+11.35。Perform a weighted average of the overlapping optimal absorbance ranges of the two optical paths near 255 nm. According to step S32, the weight of the 100 mm optical path is 100, and the weight of the 10 mm optical path is 10, and the synthesized final spectrum of this region can be obtained, at 255 nm. The absorbance after synthesis was 6.28, ie (6.10*100+8.14*10)/11. According to step S33, the value after unitization of the 100mm optical path must be increased by 0.18, that is, 6.28-6.10=0.18, and the value of the spectrum after the unitization of the 10mm optical path must be decreased by 1.86, that is, 6.28-8.14=-1.86. For the merging of the 10mm spectrum and the 1mm spectrum, step S32 is also used to perform the weighted average of the overlapping areas of 2-7, and the final absorbances at 210, 215, and 220 nm are 87.70, 55.18, and 40.68, respectively. The difference between the final spectral overlap and the 1mm spectrum 188.32-87.70=100.62, 107.84-55.18=52.66, 71.83-40.68=31.15 is not a constant, but is related to the size of the absorbance, and there is a significant unidirectional decreasing trend, so this segment is found There may be a linear baseline error, so a linear fit is performed according to the preferred method. The linear correction function of the mean value of the overlapping portion is Y=1.17X-3.28 for the 10mm spectrum and Y=0.41X+11.35 for the 1mm spectrum.
对于最大光程100mm的光谱数值上调0.18(Y=X+0.18)即可作为最终光谱;对于最小光程1mm的光谱数值进行Y=0.41X+11.35的线性修正即可作为最终光谱;而对于中间光程10mm的光谱则需进行一端为Y=X–1.86,另一端为Y=1.17X–3.28的差分修正(10mm光程与100mm光程重叠部分的修正是Y=X-1.86,10mm光程与1mm光程重叠部分的修正是Y=1.17X-3.28),例如,10mm光程光谱中,从225nm到250nm范围的吸光度需要进行调整,那么对于245nm的吸光度,其调整的斜率为两端的斜率1到1.17之间的5等分点1.034即(1.17-1)/5+1,截距为-1.86和-3.28之间的5等分点-2.144即(-3.28-(-1.86))/5+(-1/86),合成后的光谱如表1最后一列,合成前后的单位化光谱如图3所示。For the spectral value with the maximum optical path length of 100mm, it can be used as the final spectrum by increasing 0.18 (Y=X+0.18). For a spectrum with an optical path of 10mm, a differential correction of Y=X-1.86 at one end and Y=1.17X-3.28 at the other end is required (the correction for the overlapping part of the 10mm optical path and the 100mm optical path is Y=X-1.86, and the 10mm optical path is Y=X-1.86. The correction for the overlapping part with the 1mm optical path is Y=1.17X-3.28), for example, in the 10mm optical path spectrum, the absorbance in the range from 225nm to 250nm needs to be adjusted, then for the absorbance at 245nm, the slope of the adjustment is the slope of the two ends The 5th point between 1 and 1.17 is 1.034 or (1.17-1)/5+1, and the intercept is the 5th point between -1.86 and -3.28 -2.144 which is (-3.28-(-1.86))/ 5+(-1/86), the spectrum after synthesis is shown in the last column of Table 1, and the unitized spectrum before and after synthesis is shown in Figure 3.
最后应说明的是,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, but not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that it can still be Modifications are made to the technical solutions described in the foregoing embodiments, or some technical features thereof are equivalently replaced, and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110095882.4A CN112730290B (en) | 2021-01-25 | 2021-01-25 | Multi-optical path absorption spectrum spectrogram synthesis method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110095882.4A CN112730290B (en) | 2021-01-25 | 2021-01-25 | Multi-optical path absorption spectrum spectrogram synthesis method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112730290A CN112730290A (en) | 2021-04-30 |
CN112730290B true CN112730290B (en) | 2022-04-15 |
Family
ID=75594262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110095882.4A Active CN112730290B (en) | 2021-01-25 | 2021-01-25 | Multi-optical path absorption spectrum spectrogram synthesis method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112730290B (en) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2411413C (en) * | 2000-06-27 | 2011-05-03 | Alberta Research Council Inc. | Multiple pathlength spectrophotometer |
CN105510238B (en) * | 2014-09-28 | 2019-04-05 | 天津先阳科技发展有限公司 | Processing, modeling, prediction technique and the processing unit of multiposition diffusion spectroscopic data |
CN105181645B (en) * | 2015-10-10 | 2017-10-20 | 太原科技大学 | A kind of many optical path devices of screw type for measuring gas concentration |
CN105424185B (en) * | 2015-11-04 | 2017-08-11 | 清华大学 | A kind of computer assisted all band spectrometer wavelength calibration method |
CN205719927U (en) * | 2016-04-11 | 2016-11-23 | 张玉武 | Multichannel many light paths integrated form flow cell group |
CN107328724B (en) * | 2017-07-06 | 2019-08-13 | 中国计量科学研究院 | A kind of high accuracy nucleic acid concentration measuring method based on absorbance |
CN112098365B (en) * | 2020-10-20 | 2022-12-20 | 南京科远智慧科技集团股份有限公司 | Multi-spectral line gas absorption spectrum line type fitting measurement method based on orthogonal decomposition |
-
2021
- 2021-01-25 CN CN202110095882.4A patent/CN112730290B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN112730290A (en) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Terlevich et al. | The colour–magnitude relation for galaxies in the Coma cluster | |
CN109444066B (en) | A Model Transfer Method Based on Spectral Data | |
JP6954358B2 (en) | Fine particle measuring device, information processing device and information processing method | |
Whittle | Virial and jet-induced velocities in Seyfert galaxies. I-A compilation of narrow line region and host galaxy properties | |
WO2015029254A1 (en) | Data processing device for chromatograph, and method | |
JP2009162667A (en) | Method and apparatus for measuring spectroscopic absorbance | |
EP3236240B1 (en) | Spectroscopic analyzer and spectroscopic analysis method | |
US20240011956A1 (en) | System and Method for Background Removal in Spectrometry System | |
CN104990895B (en) | A kind of near infrared spectrum signal standards normal state bearing calibration based on regional area | |
CN106093134B (en) | The compensation method of metal oxide sensor array response drift | |
CN111198165A (en) | A method for determination of water quality parameters based on standardization of spectral data | |
EP2142909A1 (en) | Detector arrangement for a nondispersive infrared gas analyser and method for the detection of a measuring gas component in a gas mixture by means of such a gas analyser | |
US6154708A (en) | Method of processing and correcting spectral data in two-dimensional representation | |
Grdadolnik | Infrared difference spectroscopy: Part I. Interpretation of the difference spectrum | |
JP2018091827A (en) | Gas analyzing device and gas analyzing method | |
Wentzell et al. | Characterization of heteroscedastic measurement noise in the absence of replicates | |
Ryder et al. | Calibration, standardization, and quantitative analysis of multidimensional fluorescence (MDF) measurements on complex mixtures (IUPAC Technical Report) | |
CN112730290B (en) | Multi-optical path absorption spectrum spectrogram synthesis method | |
CN106018331B (en) | The method for estimating stability and pretreatment optimization method of multi-channel spectral system | |
JP2008532052A (en) | Signal processing method and apparatus in spectroscopic measurement using improved apodization function | |
US7002684B2 (en) | Particle size distribution measuring method, device, and program for compensating for particle concentration | |
CN118534030B (en) | A high performance liquid chromatography analysis method for tannic acid component content | |
Ozaki et al. | Spectral analysis in the NIR spectroscopy | |
Norris | Understanding and correcting the factors which affect diffuse transmittance spectra | |
Lei et al. | Window optimisation PMSC–PLS with applications to NIR spectroscopic analyses |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |