CN106885631B - A kind of spectral method of detection based on digital projection - Google Patents
A kind of spectral method of detection based on digital projection Download PDFInfo
- Publication number
- CN106885631B CN106885631B CN201710087279.5A CN201710087279A CN106885631B CN 106885631 B CN106885631 B CN 106885631B CN 201710087279 A CN201710087279 A CN 201710087279A CN 106885631 B CN106885631 B CN 106885631B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- projection
- spectrometer
- resolution
- point detector
- 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.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 28
- 230000003595 spectral effect Effects 0.000 title claims description 39
- 238000000034 method Methods 0.000 title claims description 5
- 238000001228 spectrum Methods 0.000 claims abstract description 95
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 9
- 230000003287 optical effect Effects 0.000 claims description 6
- 239000013307 optical fiber Substances 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims 1
- 229910003460 diamond Inorganic materials 0.000 claims 1
- 239000010432 diamond Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 238000013461 design Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000004907 flux Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- QNRATNLHPGXHMA-XZHTYLCXSA-N (r)-(6-ethoxyquinolin-4-yl)-[(2s,4s,5r)-5-ethyl-1-azabicyclo[2.2.2]octan-2-yl]methanol;hydrochloride Chemical compound Cl.C([C@H]([C@H](C1)CC)C2)CN1[C@@H]2[C@H](O)C1=CC=NC2=CC=C(OCC)C=C21 QNRATNLHPGXHMA-XZHTYLCXSA-N 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2803—Investigating the spectrum using photoelectric array detector
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
- G01J3/0208—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using focussing or collimating elements, e.g. lenses or mirrors; performing aberration correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/02—Details
- G01J3/0205—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows
- G01J3/0216—Optical elements not provided otherwise, e.g. optical manifolds, diffusers, windows using light concentrators or collectors or condensers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2803—Investigating the spectrum using photoelectric array detector
- G01J2003/282—Modified CCD or like
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J2003/283—Investigating the spectrum computer-interfaced
- G01J2003/2843—Processing for eliminating interfering spectra
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种基于数字投影的光谱检测方法,包括:依据光谱仪光电检测器或单点探测器的像素尺寸大小,设计一种宽度呈梯度变化的入口部分,使入射光斑最终呈梯度投影在光谱仪光电检测器或单点探测器上;根据横向投影像素点数由少到多,提取出分辨率由高到低的一组光谱图;将不同分辨率的光谱图输入数字投影算法,进行迭代计算,直到结果收敛;收敛结果可还原为分辨率最佳的光谱,即投影全部集中在光谱仪光电检测器或单点探测器上一个像素时的光谱。本发明通过在光谱检测系统前添加宽度呈梯度变化的入口以及后期配套的数据处理手段,克服传统光谱仪所获取的光谱无法同时具有高分辨率和高信噪比的问题。
The invention discloses a spectrum detection method based on digital projection, comprising: according to the pixel size of a photoelectric detector or a single-point detector of a spectrometer, designing an entrance part whose width changes in a gradient, so that the incident light spot is finally projected in a gradient On the photoelectric detector or single-point detector of the spectrometer; according to the number of horizontal projection pixels from few to many, a set of spectrograms with resolutions from high to low is extracted; spectrograms with different resolutions are input into the digital projection algorithm for iterative calculation , until the result converges; the converged result can be restored to the spectrum with the best resolution, that is, the spectrum when all the projections are concentrated on one pixel on the photodetector or single-point detector of the spectrometer. The invention overcomes the problem that the spectrum obtained by the traditional spectrometer cannot have high resolution and high signal-to-noise ratio at the same time by adding an entrance with a gradient width before the spectrum detection system and supporting data processing means in the later stage.
Description
技术领域technical field
本发明涉及光谱检测技术领域,尤其涉及一种基于数字投影的光谱检测方法。The invention relates to the technical field of spectrum detection, in particular to a digital projection-based spectrum detection method.
背景技术Background technique
光谱是复色光经过色散系统分光后,被色散开的单色光按波长(或频率)大小而依次排列的图案。其中光波由原子内部运动的电子产生,各种物质的原子内部电子的运动情况不同,所以它们发射的光波也不同。因此研究不同物质的发光和吸收光的情况,有重要的理论和实际意义。The spectrum is a pattern in which the dispersed monochromatic light is arranged sequentially according to the wavelength (or frequency) after the polychromatic light is split by the dispersion system. The light waves are generated by the electrons moving inside the atoms, and the electrons in the atoms of various substances move differently, so the light waves they emit are also different. Therefore, it is of great theoretical and practical significance to study the luminescence and light absorption of different substances.
目前光谱检测技术比较单一,其结构主要应用矩形狭缝,通过调整矩形狭缝的宽窄来控制光谱的分辨率与信噪比。狭缝的形状可以最终投影在光谱仪的光电检测器(CCD)上,投影全部集中在CCD上一个像素时,光谱分辨率最佳,投影像素数越多光谱分辨率越差。因此,当矩形狭缝较窄时,光谱可以有着较高的分辨率,但是由于光通量变小,会直接导致光谱信噪比下降,影响光谱质量;相对应的,当矩形狭缝变宽时,虽然光通量变大了可以相对提高信噪比,但是又会因为CCD上光谱投影像素过多而使得分辨率降低。也就是说,传统的光谱检测技术无法同时具有较高的分辨率与信噪比。At present, the spectral detection technology is relatively simple, and its structure mainly uses rectangular slits, and the resolution and signal-to-noise ratio of the spectrum are controlled by adjusting the width of the rectangular slits. The shape of the slit can be finally projected on the photodetector (CCD) of the spectrometer. When the projections are all concentrated on one pixel on the CCD, the spectral resolution is the best. The more the number of projected pixels, the worse the spectral resolution. Therefore, when the rectangular slit is narrow, the spectrum can have a higher resolution, but because the luminous flux becomes smaller, it will directly lead to a decrease in the spectral signal-to-noise ratio and affect the spectral quality; correspondingly, when the rectangular slit becomes wider, Although the increased luminous flux can relatively improve the signal-to-noise ratio, but the resolution will be reduced due to too many spectral projection pixels on the CCD. That is to say, traditional spectral detection technology cannot have high resolution and signal-to-noise ratio at the same time.
发明内容Contents of the invention
本发明提供了一种基于数字投影的光谱检测方法,本发明通过在光谱检测系统前添加宽度呈梯度变化的入口以及后期配套的数据处理手段,克服传统光谱仪所获取的光谱无法同时具有高分辨率和高信噪比的问题,详见下文描述:The invention provides a spectrum detection method based on digital projection. The invention overcomes the fact that the spectra acquired by traditional spectrometers cannot have high resolution at the same time by adding an entrance with a gradient width before the spectrum detection system and supporting data processing means in the later stage. and high signal-to-noise ratio, see the description below:
一种基于数字投影的光谱检测方法,所述光谱检测方法包括以下步骤:A kind of spectral detection method based on digital projection, described spectral detection method comprises the following steps:
依据光谱仪光电检测器或单点探测器的像素尺寸大小,设计一种宽度呈梯度变化的入口部分,使入射光斑最终呈梯度投影在光谱仪光电检测器或单点探测器上;According to the pixel size of the spectrometer photodetector or single-point detector, design an entrance part with a gradient width, so that the incident light spot is finally projected on the spectrometer photoelectric detector or single-point detector in a gradient;
其中,所述宽度呈梯度变化的入口部分具体为:光纤式与狭缝式,狭缝式分为固定式与可调式;Wherein, the entrance part whose width changes in a gradient is specifically: an optical fiber type and a slit type, and the slit type is divided into a fixed type and an adjustable type;
根据横向投影像素点数由少到多,提取出分辨率由高到低的一组光谱图;According to the number of horizontal projection pixels from less to more, a set of spectral images with resolution from high to low is extracted;
将不同分辨率的光谱图输入数字投影算法,进行迭代计算,直到结果收敛;Input the spectrograms with different resolutions into the digital projection algorithm, and perform iterative calculations until the results converge;
收敛结果可还原为分辨率最佳的光谱,即投影全部集中在光谱仪光电检测器或单点探测器上一个像素时的光谱。The converged result is reverted to the best-resolution spectrum, that is, the spectrum when the projections are all centered on a single pixel on the spectrometer's photodetector or single-point detector.
进一步地,所述固定式为菱形狭缝。Further, the fixed type is a diamond-shaped slit.
其中,所述使入射光斑最终呈梯度投影在光谱仪光电检测器或单点探测器上的步骤具体为:Wherein, the step of making the incident light spot finally gradiently projected on the spectrometer photodetector or single-point detector is specifically:
设计一光路结构,该光路结构包括:入口部分、准直镜、光栅、汇聚镜以及光谱仪光电检测器或单点探测器,Design an optical path structure, the optical path structure includes: entrance part, collimating mirror, grating, converging mirror and spectrometer photodetector or single point detector,
光从宽度呈梯度变化的入口部分入射,通过准直镜发散为平行光射到光栅上,光栅进行分光后,被色散开的单色光按波长或频率大小而依次排列;The light is incident from the entrance part with a gradient width, diverges into parallel light through the collimator and hits the grating. After the grating is split, the dispersed monochromatic light is arranged in sequence according to the wavelength or frequency;
再经物镜汇聚,投影在光谱仪光电检测器或单点探测器,当投影全部集中在光谱仪光电检测器或单点探测器上一个像素时光谱分辨率最佳,投影像素数越多光谱分辨率越差。After converging through the objective lens, it is projected on the photoelectric detector or single-point detector of the spectrometer. When the projection is all concentrated on one pixel of the photoelectric detector or single-point detector of the spectrometer, the spectral resolution is the best. The more projection pixels, the better the spectral resolution. Difference.
所述将不同分辨率的光谱图输入数字投影算法,进行迭代计算,直到结果收敛得步骤具体为:The specific steps of inputting the spectrograms with different resolutions into the digital projection algorithm and performing iterative calculation until the results converge are as follows:
根据数字投影算法,输入光谱组等于投影矩阵与输出光谱的乘积;According to the digital projection algorithm, the input spectrum group is equal to the product of the projection matrix and the output spectrum;
根据输入光谱组做预估计,得到一个估计光谱,对估计光谱做投影变换;将投影变换结果与输入光谱组做比,并用比值对估计光谱进行修正;Pre-estimate according to the input spectrum group, obtain an estimated spectrum, and perform projection transformation on the estimated spectrum; compare the projection transformation result with the input spectrum group, and use the ratio to correct the estimated spectrum;
反复迭代直到比值小于阈值,此时估计光谱符合预计要求,将此时的估计光谱作为结果输出。Iterate repeatedly until the ratio is less than the threshold, at this time, the estimated spectrum meets the expected requirements, and the estimated spectrum at this time is output as the result.
其中,所述将投影变换结果与输入光谱组做比,并用比值对估计光谱进行修正的步骤具体为:Wherein, the step of comparing the projection transformation result with the input spectrum group and using the ratio to correct the estimated spectrum is specifically:
其中,fdiv(n)为光谱偏差数组;ain为投影矩阵的列向量;m为投影矩阵的行数;gdiv(n)为投影偏差数组;f(n)为本次估计光谱数组;f′(n)为上一次估计光谱数组;hsum(n)为投影比重数组。Among them, f div (n) is the spectral deviation array; a in is the column vector of the projection matrix; m is the number of rows of the projection matrix; g div (n) is the projection deviation array; f(n) is the estimated spectral array; f'(n) is the last estimated spectral array; h sum (n) is the projection specific gravity array.
本发明提供的技术方案的有益效果是:本发明通过在光谱检测系统前添加宽度呈梯度变化的入口来获取分辨率不同的一组光谱图,经过数字投影可以得到该光谱仪极限分辨率的光谱。该光谱检测方法不仅能够简单、快速地克服传统光谱仪所获取的光谱无法同时具有高分辨率和高信噪比的问题,同时还可以降低激光平均强度,实现吸热材料光谱的测量,防止激光能量过高破坏样品;同时,由于入口的形状的改变,允许通过大半径的光斑,可以降低激光平均强度,实现吸热材料光谱的测量。The beneficial effect of the technical solution provided by the present invention is: the present invention acquires a group of spectrograms with different resolutions by adding an entrance with a gradient width before the spectral detection system, and the spectrum with the limit resolution of the spectrometer can be obtained through digital projection. This spectral detection method can not only simply and quickly overcome the problem that the spectra obtained by traditional spectrometers cannot have high resolution and high signal-to-noise ratio at the same time, but also can reduce the average laser intensity, realize the measurement of the spectrum of heat-absorbing materials, and prevent laser energy from Excessive damage to the sample; at the same time, due to the change in the shape of the entrance, it is allowed to pass through a large-radius spot, which can reduce the average intensity of the laser and realize the measurement of the spectrum of the endothermic material.
附图说明Description of drawings
图1是本发明提供的一种基于数字投影的光谱检测方法的流程图;Fig. 1 is the flowchart of a kind of spectral detection method based on digital projection provided by the present invention;
图2是本发明提供的仪器光路结构简图;Fig. 2 is a schematic diagram of the optical path structure of the instrument provided by the present invention;
其中,1为入口部分;2为准直镜;3为光栅;4为汇聚镜;5为光电检测器,下面以CCD为例进行说明。Among them, 1 is the entrance part; 2 is the collimating mirror; 3 is the grating; 4 is the converging mirror;
图3是本发明提供的三种可行入口;Fig. 3 is three kinds of possible entrances provided by the present invention;
其中,(a)为固定狭缝(菱形狭缝);(b)为可调狭缝;(c)为呈梯度排列的光纤。Among them, (a) is a fixed slit (rhombic slit); (b) is an adjustable slit; (c) is an optical fiber arranged in a gradient.
图4是本发明提供的普通矩形狭缝与梯度入口的区别图,以及通过四种入口后光斑在CCD上的投影图;Fig. 4 is a difference diagram between a common rectangular slit and a gradient inlet provided by the present invention, and the projection diagram of the light spot on the CCD after passing through four kinds of inlets;
图5是本发明提供的效果图;Figure 5 is an effect diagram provided by the present invention;
其中,(a)为光斑投影覆盖CCD上1个像素点时所产生的光谱;(b)为光斑投影覆盖CCD上95个像素点时所产生的光谱;(c)为数字投影算法的迭代结果;(d)为(c)和(a)之前的差异值。Among them, (a) is the spectrum produced when the spot projection covers 1 pixel on the CCD; (b) is the spectrum produced when the spot projection covers 95 pixels on the CCD; (c) is the iterative result of the digital projection algorithm ; (d) is the difference value before (c) and (a).
图6是数字投影算法示意图。Fig. 6 is a schematic diagram of a digital projection algorithm.
其中,(a)和(c)为两种原始信号;(b)和(d)为错位叠加后的信号。Among them, (a) and (c) are two kinds of original signals; (b) and (d) are signals after dislocation superimposition.
具体实施方式Detailed ways
为使本发明的目的、技术方案和优点更加清楚,下面对本发明实施方式作进一步地详细描述。In order to make the purpose, technical solution and advantages of the present invention clearer, the implementation manners of the present invention will be further described in detail below.
数字投影算法是一种数据恢复算法,信号由于错位叠加导致分辨率降低,可以通过数字投影算法进行还原,通过多次迭代可以近似还原出原始信号。The digital projection algorithm is a data recovery algorithm. The resolution of the signal is reduced due to dislocation and superposition, and it can be restored by the digital projection algorithm. The original signal can be approximated by multiple iterations.
为了使光谱同时具有高分辨率与高信噪比,克服背景技术中存在的问题,本发明实施例提出光谱检测技术与数字投影算法结合的新思路,通过改变通光入口的形状,使光通量较大的情况下依旧可以借助后期处理,得到高分辨率的光谱。In order to make the spectrum have high resolution and high signal-to-noise ratio at the same time, and overcome the problems existing in the background technology, the embodiment of the present invention proposes a new idea of combining spectrum detection technology and digital projection algorithm. By changing the shape of the light entrance, the luminous flux is relatively small In large cases, post-processing can still be used to obtain high-resolution spectra.
实施例1Example 1
一种基于数字投影技术的光谱检测方法,参见图1,该光谱检测方法主要包括以下步骤:A kind of spectral detection method based on digital projection technology, referring to Fig. 1, this spectral detection method mainly comprises the following steps:
101:依据光谱仪光电检测器(CCD)或单点探测器(下文以CCD为例)的像素尺寸大小,设计一种宽度呈梯度变化的入口,使入射光斑最终呈梯度投影在CCD上;101: According to the pixel size of the spectrometer photodetector (CCD) or single-point detector (hereinafter, CCD is taken as an example), design an entrance with a gradient width, so that the incident light spot is finally projected on the CCD in a gradient;
102:根据横向投影像素点数由少到多,可以从中提取出分辨率由高到低的一组光谱图;102: According to the number of horizontal projection pixels from less to more, a set of spectral images with resolution from high to low can be extracted from it;
其中,当投影全部集中在CCD上一个像素时光谱分辨率最佳,投影像素数越多光谱分辨率越差。Among them, the spectral resolution is the best when the projections are all concentrated on one pixel on the CCD, and the spectral resolution is poorer as the number of projected pixels increases.
103:将不同分辨率的光谱图输入数字投影算法,进行迭代计算,直到结果收敛;103: Input the spectrograms with different resolutions into the digital projection algorithm, and perform iterative calculations until the results converge;
其中,步骤103具体为:Wherein, step 103 is specifically:
根据数字投影算法,输入光谱组等于投影矩阵与输出光谱的乘积;According to the digital projection algorithm, the input spectrum group is equal to the product of the projection matrix and the output spectrum;
根据输入光谱组做预估计,得到一个估计光谱,对估计光谱做投影变换;将投影变换结果与输入光谱组做比,并用比值对估计光谱进行修正;Pre-estimate according to the input spectrum group, obtain an estimated spectrum, and perform projection transformation on the estimated spectrum; compare the projection transformation result with the input spectrum group, and use the ratio to correct the estimated spectrum;
反复迭代直到比值小于阈值,此时估计光谱符合预计要求,将此时的估计光谱作为结果输出。Iterate repeatedly until the ratio is less than the threshold, at this time, the estimated spectrum meets the expected requirements, and the estimated spectrum at this time is output as the result.
104:收敛结果即可还原为分辨率最佳的光谱,也就是投影全部集中在CCD上一个像素时的光谱。104: The convergence result can be restored to the spectrum with the best resolution, that is, the spectrum when all the projections are concentrated on one pixel on the CCD.
根据之前所述的数字投影原理,测量光谱是由分辨率最佳的光谱经过投影产生的,数字投影算法模拟这个投影过程,以测量光谱作为输入,通过多次迭代进行逆运算,还原投影过程,得到的结果即为分辨率最佳的光谱。According to the digital projection principle mentioned above, the measured spectrum is produced by projection of the spectrum with the best resolution. The digital projection algorithm simulates this projection process, takes the measured spectrum as input, and performs inverse calculation through multiple iterations to restore the projection process. The result is the best-resolution spectrum.
综上所述,本发明实施例通过上述步骤101-步骤104实现了光谱同时具有高分辨率与高信噪比的要求;同时,由于入口形状的改变允许通过大半径的光斑,可以降低激光平均强度,实现吸热材料光谱的测量。In summary, the embodiment of the present invention realizes the requirements of high resolution and high signal-to-noise ratio of the spectrum through the above steps 101 to 104; at the same time, since the change of the entrance shape allows the light spot with a large radius to pass, the laser average can be reduced. Intensity, to achieve the measurement of the spectrum of endothermic materials.
实施例2Example 2
下面结合图2至图6、以及具体实例对实施例1中的方案进行进一步地介绍,详见下文描述:Below in conjunction with Fig. 2 to Fig. 6, and specific examples, the scheme in embodiment 1 is further introduced, see the following description for details:
本发明实施例提出了一种基于数字投影的光谱检测方法,设计了一种宽度呈梯度变化的入口采集具有高信息量的光谱图,并配合后期数据处理手段,还原出同时具有高分辨率和高信噪比的光谱。该光谱检测方法可以分为入口和检测两部分,入口部分1为本发明实施例设计的重点。The embodiment of the present invention proposes a spectral detection method based on digital projection, and designs an entrance with a gradient width to collect spectral images with high information content, and cooperates with later data processing methods to restore spectral images with high resolution and Spectra with high signal-to-noise ratio. The spectral detection method can be divided into two parts: entrance and detection. The entrance part 1 is the focus of the design of the embodiment of the present invention.
入口部分1可以有两种选择,分别是光纤式与狭缝式,而狭缝式又可以分为固定式与可调式。检测部分可以匹配现有其他光谱检测仪器的检测部分,一般采用光栅分光,利用光电检测器(CCD)或者数字微镜器件(DMD)与单点探测器配合进行光谱采集。对于检测部分,本发明实施例均以光栅2分光后用CCD5进行采集为例。There are two options for the entrance part 1 , which are optical fiber type and slit type, and the slit type can be divided into fixed type and adjustable type. The detection part can match the detection part of other existing spectrum detection instruments. Generally, the grating is used to split the spectrum, and the photoelectric detector (CCD) or digital micromirror device (DMD) is used to cooperate with the single-point detector to collect the spectrum. As for the detection part, the embodiments of the present invention take the example of collecting light with CCD5 after grating 2 splitting light.
图2是本发明实施例提供的仪器光路简图,入口部分1为可调式的狭缝,可以由图3所示入口部分替代。图3是三种入口部分的示意图,包括但不仅限于如图三种,只要是宽度呈梯度变化的入口部分都满足。FIG. 2 is a schematic diagram of the optical path of the instrument provided by the embodiment of the present invention. The entrance part 1 is an adjustable slit, which can be replaced by the entrance part shown in FIG. 3 . Fig. 3 is a schematic diagram of three kinds of entrance parts, including but not limited to the three kinds shown in Fig. 3, as long as the width of the entrance part changes gradually.
如图2所示,光从狭缝入射,通过准直镜2发散为平行光射到光栅3上,光栅3进行分光后,被色散开的单色光按波长(或频率)大小而依次排列,再经物镜4汇聚,投影在光电检测器5(本发明实施例以CCD为例进行说明)上。当投影全部集中在CCD5上一个像素时光谱分辨率最佳,投影像素数越多光谱分辨率越差。As shown in Figure 2, the light is incident from the slit, diverges into parallel light through the collimating mirror 2 and hits the grating 3, after the grating 3 splits the light, the dispersed monochromatic light is in order according to the wavelength (or frequency) Arranged, converged by the objective lens 4, and projected on the photodetector 5 (the embodiment of the present invention uses a CCD as an example for illustration). The spectral resolution is the best when the projections are all concentrated on one pixel on the CCD5, and the spectral resolution is worse when the number of projected pixels is more.
本发明实施例提出了一种宽度呈梯度变化的入口部分1,如菱形狭缝(下文均以菱形狭缝为例进行说明),如图4所示,因为菱形狭缝会以1比1的比例投影在CCD5上,所以传统矩形狭缝如果想得到极限分辨率(即在CCD上的投影为一个像素),狭缝宽度就必须小于CCD5一个像素的宽度,这样会导致光通量过小,信噪比过低。菱形狭缝不同于常用的矩形狭缝,拥有着更大的光通量与更好的信噪比。光斑从菱形狭缝入射,最后在CCD5上的投影也相对应变为菱形,可以从中提取出分辨率由高到低的一组光谱图。The embodiment of the present invention proposes an inlet portion 1 whose width changes gradually, such as a diamond-shaped slit (hereinafter, the rhombus-shaped slit is used as an example), as shown in FIG. The ratio is projected on the CCD5, so if the traditional rectangular slit wants to obtain the ultimate resolution (that is, the projection on the CCD is one pixel), the width of the slit must be smaller than the width of one pixel of the CCD5, which will cause the light flux to be too small and the signal-to-noise ratio too low. The diamond-shaped slit is different from the commonly used rectangular slit, and has a larger luminous flux and a better signal-to-noise ratio. The light spot is incident from a rhombus-shaped slit, and the projection on the CCD5 is relatively rhombus-shaped, from which a set of spectral images with resolutions ranging from high to low can be extracted.
光谱分辨率降低的原理如图5所示,当光斑投影在CCD5上多个像素点时,会出现不同波长光谱的重叠,从而导致分辨率的降低,因此(a)的分辨率远低于(b)。如图5所示,选取第75组至第95组数据(即光斑投影覆盖CCD5上75个像素点至95个像素点时所产生的光谱)作为输入,其中(a)为光斑投影覆盖CCD5上1个像素点时所产生的光谱,有着最好的分辨率,(b)为光斑投影覆盖CCD5上95个像素点时所产生的光谱,(c)为迭代后得到的结果,极大的提高了光谱分辨率,(d)为(a)和(c)的差别。由此可以使光谱同时具有高分辨率与高信噪比;同时,由于狭缝形状的改变允许通过大半径的光斑,可以降低激光平均强度,从而实现吸热材料光谱的测量。The principle of spectral resolution reduction is shown in Figure 5. When the light spot is projected on multiple pixels on the CCD5, there will be overlapping of different wavelength spectra, which will lead to a reduction in resolution. Therefore, the resolution of (a) is much lower than that of ( b). As shown in Figure 5, select the 75th to 95th group of data (that is, the spectrum generated when the spot projection covers 75 pixels to 95 pixels on the CCD5) as input, where (a) is the spot projection covering the CCD5 The spectrum generated by 1 pixel has the best resolution, (b) is the spectrum generated when the spot projection covers 95 pixels on the CCD5, (c) is the result obtained after iteration, which greatly improves Spectral resolution, (d) is the difference between (a) and (c). Therefore, the spectrum can have high resolution and high signal-to-noise ratio at the same time; at the same time, because the change of the shape of the slit allows a large-radius spot to pass, the average intensity of the laser can be reduced, thereby realizing the measurement of the spectrum of the endothermic material.
数字投影算法是一种数据恢复算法,信号由于错位叠加导致分辨率降低,如图6所示,可以通过数字投影算法进行还原,通过多次迭代可以近似还原出原始信号。通过反演迭代得到极限分辨率光谱,也就是投影全部集中在CCD5上一个像素时的光谱图。The digital projection algorithm is a data recovery algorithm. The resolution of the signal is reduced due to dislocation and superposition. As shown in Figure 6, it can be restored by the digital projection algorithm, and the original signal can be approximated by multiple iterations. The limit resolution spectrum is obtained through inversion iteration, that is, the spectrum map when the projections are all concentrated on one pixel on the CCD5.
其中,实施例1中步骤103的操作具体为:Wherein, the operation of step 103 in embodiment 1 is specifically:
设输入光谱组为g0;投影矩阵为H=Hm*n(aij)m*n;H为投影矩阵;aij为投影矩阵的元素;m为投影矩阵的行数;n为投影矩阵的列数;估计光谱为f;输出光谱为f0。Let the input spectrum group be g 0 ; the projection matrix is H=H m*n (a ij ) m*n ; H is the projection matrix; a ij is the element of the projection matrix; m is the number of rows of the projection matrix; n is the projection matrix The number of columns; the estimated spectrum is f; the output spectrum is f 0 .
根据数字投影算法,有g0=H*f0;According to the digital projection algorithm, g 0 =H*f 0 ;
其中,ain为投影矩阵的列向量;hsum(n)为投影比重数组。in, a in is the column vector of the projection matrix; h sum (n) is the projection proportion array.
根据输入光谱组g0做预估计,得到一个估计光谱f=HT*g0;对估计光谱做投影变换,g=H*f;T为转置。Perform pre-estimation according to the input spectrum group g 0 to obtain an estimated spectrum f=H T *g 0 ; perform projection transformation on the estimated spectrum, g=H*f; T is transpose.
结果g与输入光谱组g0做比,gdiv=g0/g,并用比值gdiv对估计光谱进行修正。The result g is compared with the input spectrum group g 0 , g div =g 0 /g, and the estimated spectrum is corrected by the ratio g div .
fdiv(n)=(ain)m*1*gdiv(n);gdiv(n)为投影偏差数组;fdiv(n)为光谱偏差数组。f div (n)=(a in ) m*1 *g div (n); g div (n) is an array of projection deviations; f div (n) is an array of spectral deviations.
f(n)为本次估计光谱数组;f′(n)为上一次估计光谱数组。 f(n) is the estimated spectrum array for this time; f'(n) is the last estimated spectrum array.
反复迭代直到估计光谱符合预计要求,即gdiv小于阈值D,其中阈值D的设定根据实际应用中的需要进行设定,本发明实施例对此不做限制。The iterations are repeated until the estimated spectrum meets the expected requirement, that is, g div is smaller than the threshold D, where the threshold D is set according to the actual application requirements, which is not limited in the embodiment of the present invention.
最后的估计光谱作为结果输出,即f0=f。The final estimated spectrum is output as a result, ie f 0 =f.
综上所述,本发明实施例通过上述操作实现了光谱同时具有高分辨率与高信噪比的要求;同时,由于入口形状的改变允许通过大半径的光斑,可以降低激光平均强度,实现吸热材料光谱的测量。To sum up, the embodiment of the present invention realizes the requirements of high resolution and high signal-to-noise ratio of the spectrum through the above operations; at the same time, because the change of the shape of the entrance allows the light spot with a large radius to pass, the average intensity of the laser can be reduced, and the absorption can be achieved. Measurement of thermal material spectra.
本发明实施例对各器件的型号除做特殊说明的以外,其他器件的型号不做限制,只要能完成上述功能的器件均可。In the embodiments of the present invention, unless otherwise specified, the models of the devices are not limited, as long as they can complete the above functions.
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。Those skilled in the art can understand that the accompanying drawing is only a schematic diagram of a preferred embodiment, and the serial numbers of the above-mentioned embodiments of the present invention are for description only, and do not represent the advantages and disadvantages of the embodiments.
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the protection of the present invention. within range.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710087279.5A CN106885631B (en) | 2017-02-17 | 2017-02-17 | A kind of spectral method of detection based on digital projection |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710087279.5A CN106885631B (en) | 2017-02-17 | 2017-02-17 | A kind of spectral method of detection based on digital projection |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106885631A CN106885631A (en) | 2017-06-23 |
CN106885631B true CN106885631B (en) | 2018-11-30 |
Family
ID=59179379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710087279.5A Expired - Fee Related CN106885631B (en) | 2017-02-17 | 2017-02-17 | A kind of spectral method of detection based on digital projection |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106885631B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103116162A (en) * | 2012-12-14 | 2013-05-22 | 西安电子科技大学 | High-resolution sonar location method based on sparsity of objective space |
CN103163548A (en) * | 2013-03-07 | 2013-06-19 | 北京辛耕普华医疗科技有限公司 | Radioactive substance detection method, device and system based on gamma camera |
CN104266755A (en) * | 2014-09-29 | 2015-01-07 | 上海交通大学 | Spectrum measurement method and system for improving luminous flux |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS55134323A (en) * | 1979-04-06 | 1980-10-20 | Hitachi Ltd | Spectroscope |
JP6167985B2 (en) * | 2014-05-09 | 2017-07-26 | 株式会社島津製作所 | Spectroscopic sensor |
-
2017
- 2017-02-17 CN CN201710087279.5A patent/CN106885631B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103116162A (en) * | 2012-12-14 | 2013-05-22 | 西安电子科技大学 | High-resolution sonar location method based on sparsity of objective space |
CN103163548A (en) * | 2013-03-07 | 2013-06-19 | 北京辛耕普华医疗科技有限公司 | Radioactive substance detection method, device and system based on gamma camera |
CN104266755A (en) * | 2014-09-29 | 2015-01-07 | 上海交通大学 | Spectrum measurement method and system for improving luminous flux |
Also Published As
Publication number | Publication date |
---|---|
CN106885631A (en) | 2017-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jeon et al. | Compact snapshot hyperspectral imaging with diffracted rotation | |
Chapin et al. | SCUBA-2: iterative map-making with the Sub-Millimetre User Reduction Facility | |
Wang et al. | Dual-camera design for coded aperture snapshot spectral imaging | |
CN110326022B (en) | Apparatus and method for recovering hyperspectral data from image | |
Jang et al. | Experimental demonstration of adaptive infrared multispectral imaging using plasmonic filter array | |
CN104101430B (en) | A kind ofly push away the quick absolute radiation calibration method of sweeping color dispersion-type imaging spectrometer | |
Saglia et al. | The photometric classification server for Pan-STARRS1 | |
JP2013535931A (en) | Reduced image acquisition time for compression imaging devices | |
García-González et al. | A mid-infrared statistical investigation of clumpy torus model predictions | |
US10298813B2 (en) | High resolution turbulence free imaging | |
CN106815878A (en) | A kind of image rebuilding method based on the multispectral tomography of various visual angles | |
Noda et al. | Spatial deconvolution of spectropolarimetric data: an application to quiet Sun magnetic elements | |
CN107421640B (en) | Multispectral light field imaging system and method based on the principle of chromatic aberration amplification | |
Brown et al. | The dearth of UV-bright stars in M32: Implications for stellar evolution theory | |
JP7594717B2 (en) | Enhanced Resolution in Semiconductor Manufacturing Data Acquisition Instruments Using Machine Learning | |
Guaita et al. | Magellan/MMIRS near-infrared multi-object spectroscopy of nebular emission from star-forming galaxies at 2< z< 3 | |
EP3700197B1 (en) | Imaging device and method, and image processing device and method | |
Fabricius et al. | Regrowth of stellar discs in mature galaxies: the two-component nature of NGC 7217 revisited with VIRUS-W | |
CN106885631B (en) | A kind of spectral method of detection based on digital projection | |
Kobayashi et al. | Space-resolved visible spectroscopy for two-dimensional measurement of hydrogen and impurity emission spectra and of plasma flow in the edge stochastic layer of LHD | |
Tang et al. | AzTEC survey of the central molecular zone: data reduction, analysis, and preliminary results | |
CN103558160B (en) | A kind of method and system improving light spectrum image-forming spatial resolution | |
Zhang et al. | Sub-pixel dispersion model for coded aperture snapshot spectral imaging | |
CN115147287A (en) | Spectrum recovery method and spectrum recovery device for spectrum chip and electronic equipment | |
Sun et al. | Image reconstruction based on deep learning for the SPIDER Optical Interferometric System |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181130 Termination date: 20200217 |
|
CF01 | Termination of patent right due to non-payment of annual fee |