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 PDF

Info

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
Application number
CN202210147293.0A
Other languages
Chinese (zh)
Inventor
王爽
李克武
王志斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North University of China
Original Assignee
North University of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North University of China filed Critical North University of China
Priority to CN202210147293.0A priority Critical patent/CN114543998A/en
Publication of CN114543998A publication Critical patent/CN114543998A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J5/00Radiation pyrometry, e.g. infrared or optical thermometry
    • G01J2005/0077Imaging

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Radiation Pyrometers (AREA)

Abstract

本发明属于温度场时空分布测量装置技术领域,具体涉及一种基于凝视型快照式光谱成像的温度场时空分布测量装置,包括光成像收集单元、光谱成像单元、信号采样单元和电脑PC端,所述光成像收集单元的光路方向上设置有光谱成像单元,所述光谱成像单元电性连接有信号采样单元,所述信号采样单元电性连接有电脑PC端。本发明采用光谱成像技术实现,光谱成像能够同时获取待测目标的图像信息和光谱信息,图像信息能够实现温度场空间分布测量,光谱信息能够实现每个像元对应目标的温度准确测量。

Figure 202210147293

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.

Figure 202210147293

Description

基于凝视型快照式光谱成像的温度场时空分布测量装置Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging

技术领域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次的调控光谱信号被探测器阵列探测到,构成压缩测量信号siThe 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,且

Figure BDA0003508790330000041
满足
Figure BDA0003508790330000042
多次压缩测量信号Si满足Si∈RM×1,所述
Figure BDA0003508790330000043
所述
Figure BDA0003508790330000044
代表透过率调控压缩测量矩阵;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
Figure BDA0003508790330000041
Satisfy
Figure BDA0003508790330000042
The multiple compression measurement signal S i satisfies S i ∈ R M×1 , the
Figure BDA0003508790330000043
said
Figure BDA0003508790330000044
represents the transmittance regulation compression measurement matrix;

光谱信号稀疏表示为I=Ψα,所述α为光谱信号稀疏后的K-稀疏向量,所述K-稀疏向量包含K个非零元素,所述K<<N;所述Ψ是稀疏矩阵;根据压缩感知的理论框架,将

Figure BDA0003508790330000045
改写为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
Figure BDA0003508790330000045
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:

Figure BDA0003508790330000046
Figure BDA0003508790330000046

通过寻找最优光谱稀疏信号

Figure BDA0003508790330000047
然后将待测光谱信号从稀疏信号中正确地恢复出来
Figure BDA0003508790330000048
依据每个光谱测量像元上获得的光谱信号基于普朗克辐射定律进行温度反演,根据普朗克辐射定律,绝对温度T的光谱辐射温度的公式为:By finding the optimal spectrally sparse signal
Figure BDA0003508790330000047
Then correctly recover the spectral signal to be measured from the sparse signal
Figure BDA0003508790330000048
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:

Figure BDA0003508790330000049
Figure BDA0003508790330000049

所述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:

Figure BDA00035087903300000410
Figure BDA00035087903300000410

探测器阵列第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:

Figure BDA0003508790330000051
Figure BDA0003508790330000051

在短波段第i个通道的输出信号表示为:The output signal of the ith channel in the short band is expressed as:

Figure BDA0003508790330000052
Figure BDA0003508790330000052

而目标真温的求解需要构建辐亮度温度Ti与目标绝对温度T之间的关系,由公式表示为:

Figure BDA0003508790330000053
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:
Figure BDA0003508790330000053

建立100波长光谱发射率模型,表示为:ε(λ,T)=exp(α01λ+α2λ2+…+αmλm),所述α012,…,αm为常数,此时一共获得100个方程式,利用线性最小二乘法对测量数据进行线性拟合计算,获得待测目标光谱发射率和温度。A 100-wavelength spectral emissivity model is established, expressed as: ε(λ, T)=exp(α 01 λ+α 2 λ 2 +...+α m λ m ), the α 012 ,…,α 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 signal sampling unit 3 PC terminal 4, complete data processing in the computer.

在本实施例中,光成像收集单元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,且

Figure BDA0003508790330000081
满足
Figure BDA0003508790330000082
多次压缩测量光谱信号Si满足Si∈RM×1。运用这些矩阵,测量过程可以描述为:The compressed measurement signal realized by the spectral imaging unit 2 is sampled by the signal sampling unit 3 as a discrete digital quantity S i . 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
Figure BDA0003508790330000081
Satisfy
Figure BDA0003508790330000082
The multiple-compression measurement spectral signal Si satisfies Si R M×1 . Using these matrices, the measurement process can be described as:

Figure BDA0003508790330000083
Figure BDA0003508790330000083

Figure BDA0003508790330000084
代表透过率调控压缩测量矩阵。本实施实例中,探测器阵列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),Ψ是稀疏矩阵。根据压缩感知的理论框架,上式可以改写为
Figure BDA0003508790330000084
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.

Figure BDA0003508790330000085
Figure BDA0003508790330000085

通过寻找最优光谱稀疏信号

Figure BDA0003508790330000086
然后将待测光谱信号从稀疏信号中正确地恢复出来
Figure BDA0003508790330000087
By finding the optimal spectrally sparse signal
Figure BDA0003508790330000086
Then correctly recover the spectral signal to be measured from the sparse signal
Figure BDA0003508790330000087

依据每个光谱测量像元上获得的光谱信号基于普朗克辐射定律进行温度反演。普朗克定律描述了物体在不同温度下的辐射亮度随波长变化的分布规律。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:

Figure BDA0003508790330000091
Figure BDA0003508790330000091

式中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.

Figure BDA0003508790330000092
Figure BDA0003508790330000092

本实施实例中,探测器阵列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:

Figure BDA0003508790330000093
Figure BDA0003508790330000093

在短波段表示为:In the short band it is expressed as:

Figure BDA0003508790330000094
Figure BDA0003508790330000094

而目标真温的求解需要构建亮度温度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:

Figure BDA0003508790330000095
Figure BDA0003508790330000095

建立100波长光谱发射率模型,表示为:ε(λ,T)=exp(α01λ+α2λ2+…+αmλm),式中α012,…,αm为常数A 100-wavelength spectral emissivity model is established, expressed as: ε(λ, T)=exp(α 01 λ+α 2 λ 2 +…+α m λ m ), where α 012 ,…,α 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)

1.基于凝视型快照式光谱成像的温度场时空分布测量装置,其特征在于:包括光成像收集单元(1)、光谱成像单元(2)、信号采样单元(3)和电脑PC端(4),所述光成像收集单元(1)的光路方向上设置有光谱成像单元(2),所述光谱成像单元(2)电性连接有信号采样单元(3),所述信号采样单元(3)电性连接有电脑PC端(4)。1. A temperature field spatiotemporal distribution measurement device based on staring-type snapshot spectral imaging, characterized in that: comprising a light imaging collection unit (1), a spectral imaging unit (2), a signal sampling unit (3) and a computer PC terminal (4) , a spectral imaging unit (2) is provided in the optical path direction of the optical imaging collection unit (1), the spectral imaging unit (2) is electrically connected with a signal sampling unit (3), and the signal sampling unit (3) A computer PC terminal (4) is electrically connected. 2.根据权利要求1所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置,其特征在于:所述光谱成像单元(2)包括光谱成像编码模块(2-1)和探测器阵列(2-2),所述探测器阵列(2-2)设置在光谱成像编码模块(2-1)的光路方向上。2. A temperature field spatiotemporal distribution measuring device based on staring snapshot type spectral imaging according to claim 1, wherein the spectral imaging unit (2) comprises a spectral imaging coding module (2-1) and a detection A detector array (2-2), the detector array (2-2) is arranged in the optical path direction of the spectral imaging encoding module (2-1). 3.根据权利要求1所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置,其特征在于:所述光成像收集单元(1)采用反射式望远镜结构、卡塞格林望远镜或者牛顿望远镜。3. a kind of temperature field spatiotemporal distribution measuring device based on staring snapshot type spectral imaging according to claim 1, is characterized in that: described light imaging collection unit (1) adopts reflective telescope structure, Cassegrain telescope or Newtonian telescope. 4.根据权利要求2所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置,其特征在于:所述光谱成像编码模块(2-1)采用硅基平板光子晶体,所述硅基平板光子晶体通过十字孔、圆孔或三角孔平板光子晶体空气孔隙来进行目标辐射光谱伪随机透过率编码调制。4. A temperature field spatiotemporal distribution measuring device based on staring snapshot type spectral imaging according to claim 2, wherein the spectral imaging coding module (2-1) adopts a silicon-based flat plate photonic crystal, and the The silicon-based flat photonic crystal performs pseudo-random transmittance coding modulation of the target radiation spectrum through the cross hole, the circular hole or the triangular hole flat photonic crystal air hole. 5.根据权利要求2所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置,其特征在于:所述探测器阵列(2-2)采用CMOS探测器阵列,所述探测器阵列(2-2)的波长范围为400-900nm,所述探测器阵列(2-2)的像元尺寸为5.6μm,所述探测器阵列(2-2)的像素不小于2048×1024,所述探测器阵列(2-2)的帧频不低于100fps。5 . The temperature field spatiotemporal distribution measurement device based on staring-type snapshot spectral imaging according to claim 2 , wherein the detector array ( 2 - 2 ) adopts a CMOS detector array, and the detector The wavelength range of the array (2-2) is 400-900 nm, the pixel size of the detector array (2-2) is 5.6 μm, and the pixels of the detector array (2-2) are not less than 2048×1024, The frame rate of the detector array (2-2) is not lower than 100fps. 6.根据权利要求1所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置,其特征在于:所述光谱成像单元(2)光谱分辨率大于5nm,所述光谱成像单元(2)的光谱通道数为N≥100,所述光谱成像单元(2)实际压缩测量次数M≤25,所述光谱成像单元(2)的压缩比为N:M≥4:1,所述光谱成像单元(2)的每个光谱测量像元采用5×5的硅基平板光子晶体制备。6. A temperature field spatiotemporal distribution measuring device based on staring-type snapshot spectral imaging according to claim 1, characterized in that: the spectral imaging unit (2) has a spectral resolution greater than 5 nm, and the spectral imaging unit ( 2) The number of spectral channels is N≥100, the actual compression measurement times M≤25 of the spectral imaging unit (2), the compression ratio of the spectral imaging unit (2) is N:M≥4:1, the spectral imaging unit (2) Each spectral measurement pixel of the imaging unit (2) is prepared by using a 5×5 silicon-based flat plate photonic crystal. 7.根据权利要求1-6任一项所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置的测量方法,其特征在于:包括下列步骤:7. The measurement method of a temperature field spatiotemporal distribution measuring device based on staring type snapshot spectral imaging according to any one of claims 1-6, characterized in that: comprising the following steps: S1、待测目标辐射光经光成像收集单元(1)汇聚成像在光谱成像单元(2)上;S1, the radiation light of the target to be measured is collected and imaged on the spectral imaging unit (2) through the optical imaging collection unit (1); S2、光谱成像单元(2)将光信号转变为电信号,探测器阵列(2-2)对经光谱成像编码模块(2-1)伪随机透过率编码的成像光信号进行伪随机成像探测;S2. The spectral imaging unit (2) converts the optical signal into an electrical signal, and the detector array (2-2) performs pseudo-random imaging detection on the imaging optical signal encoded by the pseudo-random transmittance of the spectral imaging encoding module (2-1) ; S3、电信号经信号采样单元(3)完成采集传输至电脑PC端(4);S3, the electrical signal is collected and transmitted to the computer PC terminal (4) through the signal sampling unit (3); S4、在电脑PC端(4)中完成光谱数据和图像数据处理,然后进一步反演获得温度场的空间分布和时间分布信息。S4. Complete the spectral data and image data processing in the computer PC terminal (4), and then further invert to obtain the spatial distribution and time distribution information of the temperature field. 8.根据权利要求7所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置的测量方法,其特征在于:所述光谱成像单元(2)获得的压缩测量信号在电脑PC端(4)复原获得待测目标光谱信号,依据每个光谱测量像元上获得的光谱信号,结合普朗克辐射定律,依次获得不少于100个与光谱发射率和温度关联方程,利用线性最小二乘法对测量数据进行线性拟合计算获得待测目标光谱发射率和温度。8. The measurement method of a temperature field spatiotemporal distribution measuring device based on staring-type snapshot spectral imaging according to claim 7, characterized in that: the compressed measurement signal obtained by the spectral imaging unit (2) is at a computer PC end (4) Restoring and obtaining 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, sequentially obtain no less than 100 equations related to spectral emissivity and temperature, using the linear minimum The square method performs linear fitting calculation on the measured data to obtain the spectral emissivity and temperature of the target to be measured. 9.根据权利要求8所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置的测量方法,其特征在于:每个所述光谱测量像元获得的温度,按光谱测量像元位置进行二维排列获得温度场空间分布测量,将每个时刻获得的温度场在时间维度进行复原便掌握待测目标温度变化规律,获得温度场的时间分布测量。9 . The method for measuring a temperature field spatiotemporal distribution measuring device based on staring-type snapshot spectral imaging according to claim 8 , wherein the temperature obtained by each of the spectral measurement pixels is measured by the spectral measurement pixel. 10 . The position is arranged in two dimensions 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 law of the target to be measured, and obtain the time distribution measurement of the temperature field. 10.根据权利要求9所述的一种基于凝视型快照式光谱成像的温度场时空分布测量装置的测量方法,其特征在于:所述依次获得不少于100个与光谱发射率和温度关联方程的方法为:待测目标光强V(λ)在被探测器阵列(2-2)探测到时,由光谱成像编码模块(2-1)的透过率调控矩阵Qi(λ)进行了调控,第i次的调控光谱信号被探测器阵列(2-2)探测到,构成压缩测量信号si10. The method for measuring a temperature field spatiotemporal distribution measuring device based on staring-type snapshot spectral imaging according to claim 9, characterized in that: the sequential acquisition of no less than 100 equations related to spectral emissivity and temperature The method is as follows: when the light intensity V (λ) of the target to be measured is detected by the detector array (2-2), it is adjusted by the transmittance control matrix Q i(λ) of the spectral imaging coding module (2-1). regulation, the i-th regulation spectral signal is detected by the detector array (2-2), forming a compressed measurement signal si as si=∫V(λ)Qi(λ)η(λ)dλ i=1,2,…,Ms i =∫V (λ) Q i(λ) η (λ) dλ i=1,2,…,M 所述η(λ)表示探测器的光谱响应转换系数,所述η(λ)由探测器阵列(2-2)量子效率和探测器阵列(2-2)外围电路积分系数决定,通过对探测器阵列(2-2)进行波长辐射定标得到;The η (λ) represents the spectral response conversion coefficient of the detector, and the η (λ) is determined by the quantum efficiency of the detector array (2-2) and the integral coefficient of the peripheral circuit of the detector array (2-2). obtained by performing wavelength radiation calibration on the device array (2-2); 所述光谱成像单元(2)实现的压缩测量信号si被信号采样单元(3)采样为离散数字量Si,将待测目标光强信号离散为N维向量表示,对应光谱通道数设为N,且
Figure FDA0003508790320000031
满足
Figure FDA0003508790320000032
多次压缩测量信号Si满足Si∈RM×1,所述
Figure FDA0003508790320000033
所述
Figure FDA0003508790320000034
代表透过率调控压缩测量矩阵;
The compressed measurement signal si realized by the spectral imaging unit (2) is sampled as a discrete digital quantity Si by the signal sampling unit (3 ) , 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 as: N, and
Figure FDA0003508790320000031
Satisfy
Figure FDA0003508790320000032
The multiple compression measurement signal S i satisfies S i ∈ R M×1 , the
Figure FDA0003508790320000033
said
Figure FDA0003508790320000034
represents the transmittance regulation compression measurement matrix;
光谱信号稀疏表示为I=Ψα,所述α为光谱信号稀疏后的K-稀疏向量,所述K-稀疏向量包含K个非零元素,所述K<<N;所述Ψ是稀疏矩阵;根据压缩感知的理论框架,将
Figure FDA0003508790320000035
改写为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
Figure FDA0003508790320000035
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:
Figure FDA0003508790320000036
Figure FDA0003508790320000036
通过寻找最优光谱稀疏信号
Figure FDA0003508790320000037
然后将待测光谱信号从稀疏信号中正确地恢复出来
Figure FDA0003508790320000038
依据每个光谱测量像元上获得的光谱信号基于普朗克辐射定律进行温度反演,根据普朗克辐射定律,绝对温度T的光谱辐射温度的公式为:
By finding the optimal spectrally sparse signal
Figure FDA0003508790320000037
Then correctly recover the spectral signal to be measured from the sparse signal
Figure FDA0003508790320000038
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:
Figure FDA0003508790320000041
Figure FDA0003508790320000041
所述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:
Figure FDA0003508790320000042
Figure FDA0003508790320000042
探测器阵列(2-2)第i个通道的光谱辐射强度输出信号为:The spectral radiation intensity output signal of the i-th channel of the detector array (2-2) 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:
Figure FDA0003508790320000043
Figure FDA0003508790320000043
在短波段第i个通道的输出信号表示为:The output signal of the ith channel in the short band is expressed as:
Figure FDA0003508790320000044
Figure FDA0003508790320000044
而目标真温的求解需要构建辐亮度温度Ti与目标绝对温度T之间的关系,由公式表示为:
Figure FDA0003508790320000045
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:
Figure FDA0003508790320000045
建立100波长光谱发射率模型,表示为:ε(λ,T)=exp(α01λ+α2λ2+…+αmλm),所述α012,…,αm为常数,此时一共获得100个方程式,利用线性最小二乘法对测量数据进行线性拟合计算,获得待测目标光谱发射率和温度。A 100-wavelength spectral emissivity model is established, expressed as: ε(λ, T)=exp(α 01 λ+α 2 λ 2 +...+α m λ m ), the α 012 ,…,α 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.
CN202210147293.0A 2022-02-17 2022-02-17 Temperature field spatiotemporal distribution measurement device based on staring snapshot spectral imaging Pending CN114543998A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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