CN114543998A - Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging - Google Patents
Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging Download PDFInfo
- Publication number
- CN114543998A CN114543998A CN202210147293.0A CN202210147293A CN114543998A CN 114543998 A CN114543998 A CN 114543998A CN 202210147293 A CN202210147293 A CN 202210147293A CN 114543998 A CN114543998 A CN 114543998A
- Authority
- CN
- China
- Prior art keywords
- spectral
- signal
- temperature
- measurement
- spectral imaging
- 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.)
- Pending
Links
- 238000000701 chemical imaging Methods 0.000 title claims abstract description 77
- 238000005259 measurement Methods 0.000 title claims abstract description 61
- 230000003595 spectral effect Effects 0.000 claims abstract description 94
- 238000005070 sampling Methods 0.000 claims abstract description 16
- 230000003287 optical effect Effects 0.000 claims abstract description 11
- 238000012634 optical imaging Methods 0.000 claims abstract description 11
- 230000005855 radiation Effects 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 18
- 238000003384 imaging method Methods 0.000 claims description 16
- 108091006146 Channels Proteins 0.000 claims description 15
- 238000002834 transmittance Methods 0.000 claims description 15
- 238000007906 compression Methods 0.000 claims description 12
- 230000006835 compression Effects 0.000 claims description 11
- 239000004038 photonic crystal Substances 0.000 claims description 11
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 8
- 229910052710 silicon Inorganic materials 0.000 claims description 8
- 239000010703 silicon Substances 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 4
- 238000000691 measurement method Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000035945 sensitivity Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000009529 body temperature measurement Methods 0.000 abstract description 28
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000004861 thermometry Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
Images
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
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- 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
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J2005/0077—Imaging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明属于温度场时空分布测量装置技术领域,具体涉及一种基于凝视型快照式光谱成像的温度场时空分布测量装置,包括光成像收集单元、光谱成像单元、信号采样单元和电脑PC端,所述光成像收集单元的光路方向上设置有光谱成像单元,所述光谱成像单元电性连接有信号采样单元,所述信号采样单元电性连接有电脑PC端。本发明采用光谱成像技术实现,光谱成像能够同时获取待测目标的图像信息和光谱信息,图像信息能够实现温度场空间分布测量,光谱信息能够实现每个像元对应目标的温度准确测量。
The invention belongs to the technical field of temperature field spatiotemporal distribution measuring devices, and in particular relates to a temperature field spatiotemporal distribution measuring device based on staring-type snapshot spectral imaging, comprising an optical imaging collection unit, a spectral imaging unit, a signal sampling unit and a computer PC terminal. The optical imaging collection unit is provided with a spectral imaging unit in the direction of the optical path, the spectral imaging unit is electrically connected with a signal sampling unit, and the signal sampling unit is electrically connected with a computer PC terminal. The invention adopts the spectral imaging technology to achieve, the spectral imaging can simultaneously obtain the image information and spectral information of the target to be measured, the image information can realize the spatial distribution measurement of the temperature field, and the spectral information can realize the accurate temperature measurement of each pixel corresponding to the target.
Description
技术领域technical field
本发明属于温度场时空分布测量装置技术领域,具体涉及一种基于凝视型快照式光谱成像的温度场时空分布测量装置。The invention belongs to the technical field of temperature field spatiotemporal distribution measuring devices, and in particular relates to a temperature field spatiotemporal distribution measuring device based on staring-type snapshot spectral imaging.
背景技术Background technique
温度场测量不仅能够获得燃料燃烧、爆炸信息及效能评估,也可以获得航空、航天发动机健康状态和隐身性能评估,在民生安全、国防装备领域具有非常重要的作用。然而,温度测量,特别是温度场的时间分布、空间分布测量一直是没有解决好的技术难题。目前,温度的测量方式大致可以分为两大类:接触式和非接触式。接触式测温利用热电偶等传感器与被测对象良好接触,利用温度引起温敏材料性质的变化来实现温度测量,传感器的加入几乎不改变对象的温度,但要求被测温度不超过传感器能承受的上限温度,热电偶测温上限一般不超过2000℃,此外,接触测温通常限制了腐蚀、冲击等恶劣环境的温度测量。非接触式有红外测温、比色测温和多光谱辐射测温。红外测温采用热像仪成像方式实现温度测量,能够实现温度场的分布测量,以美国FLIRE公司的产品为代表,测温上限不超过3500℃,成像时间在ms量级,限制了该技术在超高速、超高温场的测量应用,并且红外测温要求知道被测对象光谱发射率,而被测对象光谱发射率不仅与被测对象的材质有关,而且还与温度有关,因此在实际应用中红外测温精度有限。比色测温测量两个波长的光谱辐射来实现温度测量,利用被测物体相邻波长发射率近似相等的关系,通过相邻波长辐射能量比值消除了发射率的影响,在一定温度范围内实现了较高精度的温度测量,但该技术不适用于大动态范围温度测量。The temperature field measurement can not only obtain fuel combustion, explosion information and performance evaluation, but also obtain the health status and stealth performance evaluation of aviation and aerospace engines, which plays a very important role in the fields of people's livelihood security and national defense equipment. However, the measurement of temperature, especially the measurement of time distribution and spatial distribution of temperature field, has always been an unsolved technical problem. At present, temperature measurement methods can be roughly divided into two categories: contact and non-contact. Contact temperature measurement uses sensors such as thermocouples to make good contact with the measured object, and uses temperature-induced changes in the properties of temperature-sensitive materials to achieve temperature measurement. The addition of sensors hardly changes the temperature of the object, but the measured temperature is required not to exceed the sensor can withstand The upper limit of temperature, the upper limit of thermocouple temperature measurement generally does not exceed 2000 ℃, in addition, contact temperature measurement usually limits the temperature measurement in harsh environments such as corrosion and shock. Non-contact infrared temperature measurement, colorimetric temperature measurement and multi-spectral radiation temperature measurement. Infrared temperature measurement adopts the imaging method of thermal imager to realize temperature measurement, which can realize the distribution measurement of temperature field. It is represented by the products of FLIRE Company in the United States. Ultra-high-speed and ultra-high temperature field measurement applications, and infrared temperature measurement requires knowing the spectral emissivity of the measured object, and the spectral emissivity of the measured object is not only related to the material of the measured object, but also related to temperature, so in practical applications Infrared temperature measurement has limited accuracy. Colorimetric thermometry measures the spectral radiation of two wavelengths to achieve temperature measurement. Using the relationship that the emissivity of adjacent wavelengths of the measured object is approximately equal, the influence of the emissivity is eliminated through the ratio of the radiation energy of adjacent wavelengths, and the realization is achieved within a certain temperature range. A higher accuracy temperature measurement can be achieved, but this technique is not suitable for large dynamic range temperature measurement.
多光谱测温同时实现多个光谱辐射测量,根据测量光谱信息结合普朗克辐射定律实现待测目标温度和光谱发射率的同时测量,具有测温准确,测温动态范围宽等优势,是目前非接触测温方式中极具竞争力的测温方法。但目前多光谱测量技术多采用光栅光谱仪,傅里叶变换光谱仪,不仅测量速度有限,而且无法实现成像测量。因此,多光谱测温在温度场分布测量应用受到限制。Multi-spectral temperature measurement can realize multiple spectral radiation measurements at the same time. According to the measured spectral information combined with Planck's radiation law, the simultaneous measurement of the temperature of the target to be measured and the spectral emissivity can be achieved. It has the advantages of accurate temperature measurement and wide dynamic range of temperature measurement. The most competitive temperature measurement method in the non-contact temperature measurement method. However, the current multispectral measurement technology mostly adopts grating spectrometer and Fourier transform spectrometer, which not only has limited measurement speed, but also cannot realize imaging measurement. Therefore, the application of multispectral thermometry in the measurement of temperature field distribution is limited.
发明内容SUMMARY OF THE INVENTION
针对上述的技术问题,本发明提供了一种测量准确、适用范围广、安全性高的基于凝视型快照式光谱成像的温度场时空分布测量装置。In view of the above technical problems, the present invention provides a temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging with accurate measurement, wide application range and high safety.
为了解决上述技术问题,本发明采用的技术方案为:In order to solve the above-mentioned technical problems, the technical scheme adopted in the present invention is:
基于凝视型快照式光谱成像的温度场时空分布测量装置,包括光成像收集单元、光谱成像单元、信号采样单元和电脑PC端,所述光成像收集单元的光路方向上设置有光谱成像单元,所述光谱成像单元电性连接有信号采样单元,所述信号采样单元电性连接有电脑PC端。A temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging includes an optical imaging collection unit, a spectral imaging unit, a signal sampling unit and a computer PC terminal. The optical imaging collection unit is provided with a spectral imaging unit in the direction of the optical path. The spectral imaging unit is electrically connected with a signal sampling unit, and the signal sampling unit is electrically connected with a computer PC terminal.
所述光谱成像单元包括光谱成像编码模块和探测器阵列,所述探测器阵列设置在光谱成像编码模块的光路方向上。The spectral imaging unit includes a spectral imaging coding module and a detector array, and the detector array is arranged in the optical path direction of the spectral imaging coding module.
所述光成像收集单元采用反射式望远镜结构、卡塞格林望远镜或者牛顿望远镜。The light imaging collection unit adopts a reflective telescope structure, a Cassegrain telescope or a Newtonian telescope.
所述光谱成像编码模块采用硅基平板光子晶体,所述硅基平板光子晶体通过十字孔、圆孔或三角孔平板光子晶体空气孔隙来进行目标辐射光谱伪随机透过率编码调制。The spectral imaging coding module adopts a silicon-based flat photonic crystal, and the silicon-based flat photonic crystal performs pseudo-random transmittance coding modulation of the target radiation spectrum through a cross hole, a circular hole or a triangular hole flat photonic crystal air hole.
所述探测器阵列采用CMOS探测器阵列,所述探测器阵列的波长范围为400-900nm,所述探测器阵列的像元尺寸为5.6μm,所述探测器阵列的像素不小于2048×1024,所述探测器阵列的帧频不低于100fps。The detector array adopts a CMOS detector array, the wavelength range of the detector array is 400-900 nm, the pixel size of the detector array is 5.6 μm, and the pixels of the detector array are not less than 2048×1024, The frame rate of the detector array is not lower than 100fps.
所述光谱成像单元光谱分辨率大于5nm,所述光谱成像单元的光谱通道数为N≥100,所述光谱成像单元实际压缩测量次数M≤25,所述光谱成像单元的压缩比为N:M≥4:1,所述光谱成像单元的每个光谱测量像元采用5×5的硅基平板光子晶体制备。The spectral resolution of the spectral imaging unit is greater than 5 nm, the number of spectral channels of the spectral imaging unit is N≧100, the number of actual compression measurements of the spectral imaging unit M≦25, and the compression ratio of the spectral imaging unit is N:M ≥4:1, each spectral measurement pixel of the spectral imaging unit is prepared with a 5×5 silicon-based flat photonic crystal.
一种基于凝视型快照式光谱成像的温度场时空分布测量装置的测量方法,包括下列步骤:A measurement method of a temperature field spatiotemporal distribution measurement device based on staring-type snapshot spectral imaging, comprising the following steps:
S1、待测目标辐射光经光成像收集单元汇聚成像在光谱成像单元上;S1, the radiation light of the target to be measured is converged and imaged on the spectral imaging unit through the optical imaging collection unit;
S2、光谱成像单元将光信号转变为电信号,探测器阵列对经光谱成像编码模块伪随机透过率编码的成像光信号进行伪随机成像探测;S2. The spectral imaging unit converts the optical signal into an electrical signal, and the detector array performs pseudo-random imaging detection on the imaging optical signal encoded by the pseudo-random transmittance of the spectral imaging encoding module;
S3、电信号经信号采样单元完成采集传输至电脑PC端;S3. The electrical signal is collected and transmitted to the computer PC through the signal sampling unit;
S4、在电脑PC端中完成光谱数据和图像数据处理,然后进一步反演获得温度场的空间分布和时间分布信息。S4. Complete the spectral data and image data processing in the computer PC terminal, and then further invert to obtain the spatial distribution and time distribution information of the temperature field.
所述光谱成像单元获得的压缩测量信号在电脑PC端复原获得待测目标光谱信号,依据每个光谱测量像元上获得的光谱信号,结合普朗克辐射定律,依次获得不少于100个与光谱发射率和温度关联方程,利用线性最小二乘法对测量数据进行线性拟合计算获得待测目标光谱发射率和温度。The compressed measurement signal obtained by the spectral imaging unit is restored on the computer PC to obtain the spectral signal of the target to be measured. According to the spectral signal obtained on each spectral measurement pixel, combined with Planck's radiation law, no less than 100 and The correlation equation between spectral emissivity and temperature, and the linear least squares method is used to linearly fit the measured data to obtain the spectral emissivity and temperature of the target to be measured.
每个所述光谱测量像元获得的温度,按光谱测量像元位置进行二维排列获得温度场空间分布测量,将每个时刻获得的温度场在时间维度进行复原便掌握待测目标温度变化规律,获得温度场的时间分布测量。The temperature obtained by each of the spectral measurement pixels is two-dimensionally arranged according to the position of the spectral measurement pixels to obtain the spatial distribution measurement of the temperature field, and the temperature field obtained at each moment is restored in the time dimension to grasp the temperature change rule of the target to be measured. , to obtain a time distribution measurement of the temperature field.
所述依次获得不少于100个与光谱发射率和温度关联方程的方法为:待测目标光强V(λ)在被探测器阵列探测到时,由光谱成像编码模块的透过率调控矩阵Qi(λ)进行了调控,第i次的调控光谱信号被探测器阵列探测到,构成压缩测量信号si为The method for sequentially obtaining no less than 100 equations related to spectral emissivity and temperature is: when the light intensity V (λ) of the target to be measured is detected by the detector array, the transmittance control matrix of the spectral imaging coding module The Q i(λ) is regulated, and the i-th regulation spectral signal is detected by the detector array, which constitutes the compressed measurement signal si as
si=∫V(λ)Qi(λ)η(λ)dλ i=1,2,…,Ms i =∫V (λ) Q i(λ) η (λ) dλ i=1,2,…,M
所述η(λ)表示探测器的光谱响应转换系数,所述η(λ)由探测器阵列量子效率和探测器阵列外围电路积分系数决定,通过对探测器阵列进行波长辐射定标得到;The η (λ) represents the spectral response conversion coefficient of the detector, and the η (λ) is determined by the quantum efficiency of the detector array and the integral coefficient of the peripheral circuit of the detector array, obtained by performing wavelength radiation calibration on the detector array;
所述光谱成像单元实现的压缩测量信号si被信号采样单元采样为离散数字量Si,将待测目标光强信号离散为N维向量表示,对应光谱通道数设为N,且满足多次压缩测量信号Si满足Si∈RM×1,所述所述代表透过率调控压缩测量矩阵;The compressed measurement signal si realized by the spectral imaging unit is sampled as a discrete digital quantity Si by the signal sampling unit, the light intensity signal of the target to be measured is discretely represented by an N-dimensional vector, and the number of corresponding spectral channels is set to N, and Satisfy The multiple compression measurement signal S i satisfies S i ∈ R M×1 , the said represents the transmittance regulation compression measurement matrix;
光谱信号稀疏表示为I=Ψα,所述α为光谱信号稀疏后的K-稀疏向量,所述K-稀疏向量包含K个非零元素,所述K<<N;所述Ψ是稀疏矩阵;根据压缩感知的理论框架,将改写为S=QV=QΨα=Aα,所述A=QΨ代表传感矩阵,是由压缩测量矩阵Q和稀疏矩阵Ψ共同决定的;然后对S=QV=QΨα=Aα的逆问题进行求解:The spectral signal is sparsely represented as I=Ψα, the α is the K-sparse vector after the spectral signal is sparsed, the K-sparse vector contains K non-zero elements, and the K<<N; the Ψ is a sparse matrix; According to the theoretical framework of compressed sensing, the Rewritten as S=QV=QΨα=Aα, the A=QΨ represents the sensing matrix, which is jointly determined by the compressed measurement matrix Q and the sparse matrix Ψ; then the inverse problem of S=QV=QΨα=Aα is solved:
通过寻找最优光谱稀疏信号然后将待测光谱信号从稀疏信号中正确地恢复出来依据每个光谱测量像元上获得的光谱信号基于普朗克辐射定律进行温度反演,根据普朗克辐射定律,绝对温度T的光谱辐射温度的公式为:By finding the optimal spectrally sparse signal Then correctly recover the spectral signal to be measured from the sparse signal According to the spectral signal obtained on each spectral measurement pixel, the temperature is inverted based on Planck's radiation law. According to Planck's radiation law, the formula for the spectral radiation temperature of the absolute temperature T is:
所述L(λ,T)为物体的辐射亮度,所述λ为波长,所述T为绝对温度;所述ε(λ,T)为物体的光谱发射率;所述C1=3.7415×108W·μm4·m-2;所述C2=1.43879×104μm4·K;The L(λ, T) is the radiance of the object, the λ is the wavelength, and the T is the absolute temperature; the ε(λ, T) is the spectral emissivity of the object; the C 1 =3.7415×10 8 W·μm 4 ·m −2 ; the C 2 =1.43879×10 4 μm 4 ·K;
在短波段时C2/λT>>1,此时用维恩公式来近似代替普朗克公式:When C 2 /λT>>1 in the short waveband, the Wien formula is used to approximate the Planck formula:
探测器阵列第i个通道的光谱辐射强度输出信号为:The spectral radiation intensity output signal of the i-th channel of the detector array is:
Vi=τ(λi)S(λi)L(λi,T),所述τ(λi)为光谱透过率,S(λi)为探测器灵敏度;V i =τ(λ i )S(λ i )L(λ i ,T), the τ(λ i ) is the spectral transmittance, and S(λ i ) is the detector sensitivity;
则第i个通道的输出信号表示为:Then the output signal of the i-th channel is expressed as:
在短波段第i个通道的输出信号表示为:The output signal of the ith channel in the short band is expressed as:
而目标真温的求解需要构建辐亮度温度Ti与目标绝对温度T之间的关系,由公式表示为: The solution of the target true temperature needs to construct the relationship between the radiance temperature T i and the target absolute temperature T, which is expressed by the formula:
建立100波长光谱发射率模型,表示为:ε(λ,T)=exp(α0+α1λ+α2λ2+…+αmλm),所述α0,α1,α2,…,αm为常数,此时一共获得100个方程式,利用线性最小二乘法对测量数据进行线性拟合计算,获得待测目标光谱发射率和温度。A 100-wavelength spectral emissivity model is established, expressed as: ε(λ, T)=exp(α 0 +α 1 λ+α 2 λ 2 +...+α m λ m ), the α 0 ,α 1 ,α 2 ,…,α m is a constant, a total of 100 equations are obtained at this time, and the linear least squares method is used to linearly fit the measured data to obtain the spectral emissivity and temperature of the target to be measured.
本发明与现有技术相比,具有的有益效果是:Compared with the prior art, the present invention has the following beneficial effects:
1、本发明采用光谱成像技术实现,光谱成像能够同时获取待测目标的图像信息和光谱信息,图像信息能够实现温度场空间分布测量,光谱信息能够实现每个像元对应目标的温度准确测量。1. The present invention is realized by using spectral imaging technology. Spectral imaging can simultaneously obtain image information and spectral information of the target to be measured, the image information can realize the measurement of the spatial distribution of the temperature field, and the spectral information can realize the accurate measurement of the temperature of each pixel corresponding to the target.
2、本发明采用多光谱辐射测温方式实现温度测量,能够同时实现待测目标光谱发射率和温度测量,不需要实现知道光谱发射率,温度测量准确,较宽范围的光谱信息获取能够实现较大动态范围温度测量。2. The present invention adopts the multi-spectral radiation temperature measurement method to realize temperature measurement, which can realize the spectral emissivity and temperature measurement of the target to be measured at the same time, and does not need to know the spectral emissivity, the temperature measurement is accurate, and the acquisition of a wide range of spectral information can achieve a relatively wide range of spectral information acquisition. Large dynamic range temperature measurement.
3、本发明采用凝视型快照式光谱成像方式,光谱图像信息获取速度快,可以实现待测目标实时连续的温度测量,在温度场空间分布测量基础上获得每个温度场随时间变化规律。3. The present invention adopts the staring snapshot type spectral imaging method, the spectral image information acquisition speed is fast, the real-time continuous temperature measurement of the target to be measured can be realized, and the time-varying law of each temperature field can be obtained on the basis of the measurement of the spatial distribution of the temperature field.
4、本发明采用凝视型快照式光谱成像来实现温度时空分布测量,为非接触光谱成像测量方法,本发明采用望远成像镜头将温度场成像在光谱成像采集模块上,为非接触遥测方法,能够规避爆炸、燃烧场冲击、腐蚀、振动等恶劣环境影响。4. The present invention adopts staring-type snapshot spectral imaging to realize temperature distribution measurement, which is a non-contact spectral imaging measurement method. The present invention uses a telephoto imaging lens to image the temperature field on the spectral imaging acquisition module, which is a non-contact telemetry method. It can avoid harsh environmental influences such as explosion, combustion field impact, corrosion, vibration and so on.
附图说明Description of drawings
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引申获得其它的实施附图。In order to illustrate the embodiments of the present invention or the technical solutions in the prior art more clearly, the following briefly introduces the accompanying drawings that are required to be used in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only exemplary, and for those of ordinary skill in the art, other implementation drawings can also be derived from the provided drawings without any creative effort.
本说明书所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容能涵盖的范围内。The structures, proportions, sizes, etc. shown in this specification are only used to cooperate with the contents disclosed in the specification, so as to be understood and read by those who are familiar with the technology, and are not used to limit the conditions for the implementation of the present invention, so there is no technical The substantive meaning, any modification of the structure, the change of the proportional relationship or the adjustment of the size, without affecting the effect that the present invention can produce and the purpose that can be achieved, should still fall within the technical content disclosed in the present invention. within the scope of coverage.
图1为本发明的结构示意图;Fig. 1 is the structural representation of the present invention;
图2为本发明光谱成像单元示意图。FIG. 2 is a schematic diagram of the spectral imaging unit of the present invention.
其中:1为光成像收集单元,2为光谱成像单元,3为信号采样单元,4为电脑PC端,2-1为光谱成像编码模块,2-2为探测器阵列。Among them: 1 is an optical imaging collection unit, 2 is a spectral imaging unit, 3 is a signal sampling unit, 4 is a computer PC terminal, 2-1 is a spectral imaging coding module, and 2-2 is a detector array.
具体实施方式Detailed ways
为使本发明实施例的目的、技术方案和优点更加清楚,下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例,这些描述只是为进一步说明本发明的特征和优点,而不是对本发明权利要求的限制;基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。In order to make the purposes, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be described clearly and completely below. Obviously, the described embodiments are only part of the embodiments of the present application, not All embodiments, these descriptions are only to further illustrate the features and advantages of the present invention, rather than to limit the claims of the present invention; based on the embodiments in this application, those of ordinary skill in the art can obtain without creative work. All other embodiments belong to the scope of protection of the present application.
下面结合附图和实施例,对本发明的具体实施方式做进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。The specific embodiments of the present invention will be described in further detail below with reference to the accompanying drawings and embodiments. The following examples are intended to illustrate the present invention, but not to limit the scope of the present invention.
在本申请的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体的连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本申请中的具体含义。In the description of this application, it should be noted that, unless otherwise expressly specified and limited, the terms "installed", "connected" and "connected" should be understood in a broad sense, for example, it may be a fixed connection or a detachable connection Connection, or integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium, and it can be the internal communication of two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present application can be understood in specific situations.
如图1所示,待测目标辐射光经光成像收集单元1汇聚成像在光谱成像单元2上,光谱成像单元2将光信号转变为电信号,电信号经信号采样单元3完成采集传输至电脑PC端4,在电脑中完成数据处理。As shown in FIG. 1 , the radiation light of the target to be measured is collected and imaged on the spectral imaging unit 2 by the optical imaging collection unit 1 . The spectral imaging unit 2 converts the optical signal into an electrical signal, and the electrical signal is collected and transmitted to the computer through the
在本实施例中,光成像收集单元1工作在可见/近红外波段,采用反射式望远镜结构实现非接触遥测,光成像收集单元1选用卡塞格林望远镜或者牛顿望远镜。将待测目标辐射光收集聚焦成像到光谱成像单元2上。In this embodiment, the optical imaging collection unit 1 works in the visible/near-infrared band, and uses a reflective telescope structure to realize non-contact telemetry. The optical imaging collection unit 1 selects a Cassegrain telescope or a Newtonian telescope. The radiant light of the target to be measured is collected and imaged on the spectral imaging unit 2 .
如图2所示光谱成像单元示意图,成像单元包括光谱成像编码模块2-1和探测器阵列2-2。光谱编码模块2-1,采用硅基平板光子晶体设计,采用十字孔、圆孔、三角孔等平板光子晶体空气孔隙来进行目标辐射光谱伪随机透过率编码调制;探测器阵列2-2选用CMOS探测器阵列,对经光谱编码模块2-1伪随机透过率编码的成像光信号进行伪随机成像探测。As shown in FIG. 2, a schematic diagram of a spectral imaging unit is shown, and the imaging unit includes a spectral imaging coding module 2-1 and a detector array 2-2. Spectral coding module 2-1 is designed with a silicon-based flat photonic crystal, and uses cross holes, circular holes, triangular holes and other flat photonic crystal air holes to perform pseudo-random transmittance coding modulation of the target radiation spectrum; the detector array 2-2 is selected The CMOS detector array performs pseudo-random imaging detection on the imaging optical signal encoded by the pseudo-random transmittance of the spectral encoding module 2-1.
如图2所示光谱成像是光谱成像单元2将待测目标具有二维空间信息(x,y)和一维光谱信息。目标每个成像单元入射光在整个光谱范围内通过光谱成像单元2的光谱编码模块2-1入射到探测器阵列2-2。待测目标光强V(λ)在被探测器探测到时,由编码板的透过率调控矩阵Qi(λ)进行了调控,第i次的调控光谱信号被探测器探测到,构成压缩测量信号si可描述为As shown in FIG. 2, spectral imaging means that the spectral imaging unit 2 provides the target to be measured with two-dimensional spatial information (x, y) and one-dimensional spectral information. The incident light of each imaging unit of the target is incident on the detector array 2-2 through the spectral encoding module 2-1 of the spectral imaging unit 2 in the entire spectral range. When the light intensity V (λ) of the target to be measured is detected by the detector, it is regulated by the transmittance regulation matrix Q i(λ) of the coding plate, and the i-th regulation spectral signal is detected by the detector, which constitutes a compression The measurement signal si can be described as
si=∫V(λ)Qi(λ)η(λ)dλ i=1,2,…,Ms i =∫V (λ) Q i(λ) η (λ) dλ i=1,2,…,M
其中,η(λ)表示探测器的光谱响应转换系数,主要由探测器量子效率和探测器外围电路积分系数决定,通常通过对探测器进行波长辐射定标得到。Among them, η (λ) represents the spectral response conversion coefficient of the detector, which is mainly determined by the quantum efficiency of the detector and the integral coefficient of the peripheral circuit of the detector, which is usually obtained by calibrating the wavelength radiation of the detector.
光谱成像单元2实现的压缩测量信号,被信号采样单元3采样为离散数字量Si。将待测目标光强信号离散为N维向量表示,对应光谱通道数设为N,且满足多次压缩测量光谱信号Si满足Si∈RM×1。运用这些矩阵,测量过程可以描述为:The compressed measurement signal realized by the spectral imaging unit 2 is sampled by the
代表透过率调控压缩测量矩阵。本实施实例中,探测器阵列2-2选用CMOS探测器阵列有效工作波长选择工作在400-900nm,像元尺寸为5.6μm,像素不小于2048×1024,帧频不低于100fps。整个光谱成像单元2光谱分辨率优于5nm,光谱通道数N设置为N≥100,压缩比N:M≥4:1,其中光谱编码模块2-1每个光谱测量像元采用5×5的硅基平板光子晶体制备。光谱信号稀疏表示为I=Ψα,其中α为光谱信号稀疏后的K-稀疏向量(包含K个非零元素,其中K<<N),Ψ是稀疏矩阵。根据压缩感知的理论框架,上式可以改写为 Represents the transmittance modulation compression measurement matrix. In this embodiment, the detector array 2-2 selects a CMOS detector array with an effective working wavelength of 400-900 nm, a pixel size of 5.6 μm, a pixel size of not less than 2048×1024, and a frame rate of not less than 100fps. The spectral resolution of the entire spectral imaging unit 2 is better than 5 nm, the number of spectral channels N is set to N≥100, and the compression ratio is N:M≥4:1. Each spectral measurement pixel of the spectral coding module 2-1 uses a 5×5 Preparation of silicon-based flat-plate photonic crystals. The spectral signal is sparsely expressed as I=Ψα, where α is a K-sparse vector (containing K non-zero elements, where K<<N) after the spectral signal is sparsed, and Ψ is a sparse matrix. According to the theoretical framework of compressed sensing, the above formula can be rewritten as
S=QV=QΨα=AαS=QV=QΨα=Aα
上式中A=QΨ代表传感矩阵,是由压缩测量矩阵Q和稀疏矩阵Ψ共同决定的。然后对上式的逆问题进行求解。In the above formula, A=QΨ represents the sensing matrix, which is jointly determined by the compressed measurement matrix Q and the sparse matrix Ψ. Then solve the inverse problem of the above equation.
通过寻找最优光谱稀疏信号然后将待测光谱信号从稀疏信号中正确地恢复出来 By finding the optimal spectrally sparse signal Then correctly recover the spectral signal to be measured from the sparse signal
依据每个光谱测量像元上获得的光谱信号基于普朗克辐射定律进行温度反演。普朗克定律描述了物体在不同温度下的辐射亮度随波长变化的分布规律。According to the spectral signal obtained on each spectral measurement pixel, temperature inversion is performed based on Planck's radiation law. Planck's law describes the distribution of the radiance of an object at different temperatures as a function of wavelength.
根据普朗克辐射定律,绝对温度T的光谱辐射温度可以由公式表示:According to Planck's law of radiation, the spectral radiation temperature of the absolute temperature T can be expressed by the formula:
式中L(λ,T)为物体的辐射亮度(W·m-2·μm-1sr-1),λ为波长(μm),T为绝对温度(K);ε(λ,T)为物体的光谱发射率;C1=3.7415×108W·μm4·m-2;C2=1.43879×104μm4·K。where L(λ,T) is the radiance of the object (W·m -2 ·μm -1 sr -1 ), λ is the wavelength (μm), T is the absolute temperature (K); ε(λ,T) is Spectral emissivity of the object; C 1 =3.7415×10 8 W·μm 4 ·m −2 ; C 2 =1.43879×10 4 μm 4 ·K.
在短波段时C2/λT>>1,此时可以用维恩公式来近似代替普朗克公式When C 2 /λT>>1 in the short waveband, the Wien formula can be used to approximate the Planck formula.
本实施实例中,探测器阵列2-2选用CMOS探测器阵列有效工作波长选择工作在400-900nm,整个光谱成像单元2光谱分辨率优于5nm,光谱通道数N设置为N≥100。其中第i个通道的光谱辐射强度输出信号为:In this embodiment, the detector array 2-2 selects a CMOS detector array with an effective working wavelength of 400-900 nm, the spectral resolution of the entire spectral imaging unit 2 is better than 5 nm, and the number of spectral channels N is set to N≥100. The spectral radiation intensity output signal of the i-th channel is:
Vi=τ(λi)S(λi)L(λi,T)V i =τ(λ i )S(λ i )L(λ i ,T)
式中τ(λi)为光谱透过率,S(λi)为探测器灵敏度where τ(λ i ) is the spectral transmittance and S(λ i ) is the detector sensitivity
则第i个通道的输出信号表示为:Then the output signal of the i-th channel is expressed as:
在短波段表示为:In the short band it is expressed as:
而目标真温的求解需要构建亮度温度Ti与目标绝对温度T之间的关系,The solution of the target true temperature needs to construct the relationship between the brightness temperature T i and the target absolute temperature T,
由公式表示为:It is represented by the formula as:
建立100波长光谱发射率模型,表示为:ε(λ,T)=exp(α0+α1λ+α2λ2+…+αmλm),式中α0,α1,α2,…,αm为常数A 100-wavelength spectral emissivity model is established, expressed as: ε(λ, T)=exp(α 0 +α 1 λ+α 2 λ 2 +…+α m λ m ), where α 0 ,α 1 ,α 2 ,…,α m is a constant
此时一共可以获得100个方程式,利用线性最小二乘法对测量数据进行线性拟合计算,即可获得待测目标光谱发射率和温度。将每个光谱像元获得的温度按成像像元位置进行二维排列,便可完成温度场空间分布测量,将每个时刻获得的温度场在时间维度进行复原便可掌握待测目标温度变化规律,获得温度场的时间分布测量。At this time, a total of 100 equations can be obtained, and the measured data can be linearly fitted and calculated by the linear least squares method to obtain the spectral emissivity and temperature of the target to be measured. The temperature obtained by each spectral pixel is arranged two-dimensionally according to the position of the imaging pixel to complete the measurement of the spatial distribution of the temperature field, and the temperature change law of the target to be measured can be grasped by restoring the temperature field obtained at each moment in the time dimension. , to obtain a time distribution measurement of the temperature field.
上面仅对本发明的较佳实施例作了详细说明,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化,各种变化均应包含在本发明的保护范围之内。Only the preferred embodiments of the present invention have been described in detail above, but the present invention is not limited to the above-mentioned embodiments, and within the scope of knowledge possessed by those of ordinary skill in the art, various aspects can also be made without departing from the purpose of the present invention. Various changes should be included within the protection scope of the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210147293.0A CN114543998A (en) | 2022-02-17 | 2022-02-17 | Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210147293.0A CN114543998A (en) | 2022-02-17 | 2022-02-17 | Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114543998A true CN114543998A (en) | 2022-05-27 |
Family
ID=81675903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210147293.0A Pending CN114543998A (en) | 2022-02-17 | 2022-02-17 | Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114543998A (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110186566A (en) * | 2019-05-30 | 2019-08-30 | 上海交通大学 | Two-dimentional true temperature field imaging method and system based on the multispectral thermometric of light-field camera |
CN111458044A (en) * | 2020-05-28 | 2020-07-28 | 中国工程物理研究院流体物理研究所 | Transient temperature measurement device and method based on snapshot spectral imaging technology |
CN111664941A (en) * | 2020-07-08 | 2020-09-15 | 中北大学 | Compressed sensing spectrum measuring device and method for electro-optically regulating transmittance coding |
CN112345074A (en) * | 2020-11-11 | 2021-02-09 | 中北大学 | A chip-level spaceborne hyperspectral imaging detector and its spectral imaging method |
-
2022
- 2022-02-17 CN CN202210147293.0A patent/CN114543998A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110186566A (en) * | 2019-05-30 | 2019-08-30 | 上海交通大学 | Two-dimentional true temperature field imaging method and system based on the multispectral thermometric of light-field camera |
CN111458044A (en) * | 2020-05-28 | 2020-07-28 | 中国工程物理研究院流体物理研究所 | Transient temperature measurement device and method based on snapshot spectral imaging technology |
CN111664941A (en) * | 2020-07-08 | 2020-09-15 | 中北大学 | Compressed sensing spectrum measuring device and method for electro-optically regulating transmittance coding |
CN112345074A (en) * | 2020-11-11 | 2021-02-09 | 中北大学 | A chip-level spaceborne hyperspectral imaging detector and its spectral imaging method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108562363B (en) | An Accurate Measurement Method of Infrared Radiation Characteristic Transient Temperature Field | |
CN101476939B (en) | Double-CCD temperature field measuring apparatus and method | |
CN101943604B (en) | Temperature-measurement imaging system and measuring method thereof | |
CN103076101B (en) | Method for calibrating thermal infrared imager pixel point | |
CN110487134B (en) | A kind of explosive flame burning speed and temperature measuring device and measuring method | |
CN102736128A (en) | Method and device for processing unmanned plane optical remote sensing image data | |
CN113418613B (en) | High-temperature transient measurement system and method based on multispectral colorimetry | |
CN102538983A (en) | CCD (Charge Coupled Device) temperature measuring device | |
CN108489614A (en) | Compact-structure type full-Stokes vector polarization imaging device based on double-adhesive digital coding birefringent crystal slice | |
Schumann et al. | Infrared hyperspectral imaging Fourier transform and dispersive spectrometers: comparison of signal-to-noise-based performance | |
CN113588115A (en) | Temperature measurement method based on multispectral colorimetry | |
Anderson et al. | The physical basis of current infrared remote-sensing techniques and the interpretation of data from aerial surveys | |
CN114543998A (en) | Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging | |
CN114449079A (en) | A kind of high temperature measurement device and method based on mobile phone camera | |
Sauer et al. | Numerical assessment of uncertainty and dynamic range expansion of multispectral image-based pyrometry | |
Ivanov et al. | Thermal detectors of uncooled multi-element infrared imaging arrays. I. Thermally insulated elements | |
CN114894737B (en) | A spectral reflectance reconstruction method based on infrared images | |
Méndez-Rial et al. | A high-speed MWIR uncooled multi-aperture snapshot spectral imager for IR surveillance and monitoring | |
RU2727349C1 (en) | Method of thermography of a remote object | |
Malone et al. | Staring MWIR, LWIR and 2-color and scanning LWIR polarimetry technology | |
Wei et al. | A modified single-channel algorithm for estimating land surface temperature from UAV TIR imagery | |
CN115112244A (en) | Satellite-borne thermal infrared camera absolute radiometric calibration processing method | |
Jitsufuchi | A Multiband Uncooled Infrared Camera for Measuring Volcanic SO 2 Gas Concentration and Temperature Distributions | |
Hinnrichs et al. | Comparison of QWIP to HgCdTe detectors for gas imaging | |
CN101089658A (en) | Method of non-refrigeration focal surface infrared cloud detection |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20220527 |