CN105675146B - 基于压缩感知的双波长三维温度场成像设备、系统及方法 - Google Patents
基于压缩感知的双波长三维温度场成像设备、系统及方法 Download PDFInfo
- Publication number
- CN105675146B CN105675146B CN201610080287.2A CN201610080287A CN105675146B CN 105675146 B CN105675146 B CN 105675146B CN 201610080287 A CN201610080287 A CN 201610080287A CN 105675146 B CN105675146 B CN 105675146B
- Authority
- CN
- China
- Prior art keywords
- wavelength
- light
- optical radiation
- path
- measured
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 74
- 230000009977 dual effect Effects 0.000 title claims abstract description 43
- 230000005855 radiation Effects 0.000 claims abstract description 596
- 239000011159 matrix material Substances 0.000 claims abstract description 239
- 238000001514 detection method Methods 0.000 claims abstract description 105
- 230000009466 transformation Effects 0.000 claims abstract description 24
- 230000000295 complement effect Effects 0.000 claims abstract description 18
- 230000003287 optical effect Effects 0.000 claims description 491
- 238000005259 measurement Methods 0.000 claims description 134
- 238000004422 calculation algorithm Methods 0.000 claims description 53
- 238000001914 filtration Methods 0.000 claims description 41
- 238000012360 testing method Methods 0.000 claims description 29
- 230000005693 optoelectronics Effects 0.000 claims description 27
- 238000009826 distribution Methods 0.000 claims description 20
- 238000013178 mathematical model Methods 0.000 claims description 20
- 230000036961 partial effect Effects 0.000 claims description 20
- 230000001678 irradiating effect Effects 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 9
- 239000004973 liquid crystal related substance Substances 0.000 claims description 8
- 230000001131 transforming effect Effects 0.000 claims description 8
- ADXSZLCTQCWMTE-UHFFFAOYSA-N 4-n-(3-bromophenyl)quinazoline-4,6,7-triamine Chemical compound C=12C=C(N)C(N)=CC2=NC=NC=1NC1=CC=CC(Br)=C1 ADXSZLCTQCWMTE-UHFFFAOYSA-N 0.000 claims description 7
- 102100031262 Deleted in malignant brain tumors 1 protein Human genes 0.000 claims description 7
- 101000844721 Homo sapiens Deleted in malignant brain tumors 1 protein Proteins 0.000 claims description 7
- ZINJLDJMHCUBIP-UHFFFAOYSA-N ethametsulfuron-methyl Chemical compound CCOC1=NC(NC)=NC(NC(=O)NS(=O)(=O)C=2C(=CC=CC=2)C(=O)OC)=N1 ZINJLDJMHCUBIP-UHFFFAOYSA-N 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 5
- 238000007906 compression Methods 0.000 claims description 5
- 230000008447 perception Effects 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 21
- 238000004861 thermometry Methods 0.000 abstract description 16
- 238000004364 calculation method Methods 0.000 abstract description 13
- 238000009529 body temperature measurement Methods 0.000 description 84
- 230000000694 effects Effects 0.000 description 40
- 230000002829 reductive effect Effects 0.000 description 25
- 238000010586 diagram Methods 0.000 description 22
- 230000008569 process Effects 0.000 description 17
- 238000005070 sampling Methods 0.000 description 16
- 238000003860 storage Methods 0.000 description 16
- 230000035945 sensitivity Effects 0.000 description 14
- 238000010521 absorption reaction Methods 0.000 description 10
- 230000007613 environmental effect Effects 0.000 description 9
- 230000003595 spectral effect Effects 0.000 description 8
- 230000008901 benefit Effects 0.000 description 6
- 238000001931 thermography Methods 0.000 description 6
- ABIFUJNCKIMWRZ-JGVFFNPUSA-N (2r,4s)-4-(3-phosphonopropyl)piperidine-2-carboxylic acid Chemical compound OC(=O)[C@H]1C[C@@H](CCCP(O)(O)=O)CCN1 ABIFUJNCKIMWRZ-JGVFFNPUSA-N 0.000 description 5
- 238000004590 computer program Methods 0.000 description 5
- 238000012937 correction Methods 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 230000005288 electromagnetic effect Effects 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 4
- 230000007935 neutral effect Effects 0.000 description 4
- 238000012634 optical imaging Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011166 aliquoting Methods 0.000 description 3
- 230000002146 bilateral effect Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000013519 translation Methods 0.000 description 3
- 230000005457 Black-body radiation Effects 0.000 description 2
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 2
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000017105 transposition Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 230000004297 night vision Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004071 soot Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/10—Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors
-
- 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
- G01J5/48—Thermography; Techniques using wholly visual means
- G01J5/485—Temperature profile
-
- 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
本申请提供了一种基于压缩感知的双波长三维温度场成像设备，包括脉冲光源；配置为接收待测对象的光辐射的光辐射调制装置，其加载预设的多个掩膜，将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，所述多个掩膜根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成；布置在所述第一路径上的第一滤光元件和第一探测装置；布置在所述第二路径上的第二滤光元件和第二探测装置；布置在第一路径或者第二路径上的时幅变换器；温度确定装置以及图像生成装置。本申请还提供了基于该设备的系统和方法。其将压缩感知理论、双波长测温技术、互补调制技术、单光子探测技术和时间相关技术相结合，重建出待测对象的三维红外热图像。
Description
技术领域
本发明涉及双波长温度场成像领域，特别涉及一种基于压缩感知的双波长三维温度场成像设备、系统及方法。
背景技术
在航空航天、冶金及汽车制造业等领域，常常需要对待测对象以及各种在线工件进行快速、实时监控，以最大限度地减少事故隐患、提高产品的安全性能及其质量。沿用传统的接触式测温仪进行测量，虽然精度高，但必须使探测器接触待测对象。可是在某些特殊的场合下(如对发动机燃烧室及高温炉中的火焰进行测温时)无法使用接触式测温仪，由此产生了非接触式测温方法。红外测温方法就属于一种非接触式测温方法，该方法通过检测物体表面发射的能量来测量温度，具有测温范围广、响应速度快和不明显破坏测温度场等特点，被广泛应用于工业各个方面。
红外测温方法主要基于黑体辐射理论，黑体是个理想化的物理模型，而自然界中实际存在的物体(测温对象)，其吸收能力及辐射能力都比黑体小，称为灰体。根据普朗克辐射定律，一个绝对温度为T的黑体，单位表面积在波长λ1、λ2(λ1的附近单位波长)的间隔内向整个半球空间发射的辐射功率(简称光谱辐射度)为E0(λ，T)，而灰体的光谱辐射能量的计算公式为：E(λ，T)＝ε(λ，T)E0(λ，T)，其中ε(λ，T)为该灰体的辐射率。
现有技术中红外测温主要经历了三个阶段的发展。
第一阶段：传统的红外测温设备，一律按黑体的热辐射定律来设计。该设计方式假设红外测温仪实际接收到的热辐射是与被测物的光谱辐射能量E(λ，T)是成比例，故在使用红外测温仪时，必须求出待测对象的辐射率ε(λ，T)数值，即进行辐射率修正。遗憾的是，该辐射率ε(λ，T)与待测对象的材料、表面状态、波长，温度以及辐射条件、环境因素等均有复杂的关系，因而很难准确测定ε(λ，T)，同时由于在某些情况下待测对象的辐射率随温度变化太大，因而传统红外测温仪还存在着较大的误差。
第二阶段：科学家们为解决传统的红外测温设备存在的问题，研究出了基于单波长窄带滤波红外测温技术。可是由于存在诸如水蒸气等的周围环境对红外线的吸收，在很大幅度上影响了单波长红外测温的准确性。
第三阶段：利用双波长滤波红外测温技术来进行红外测温。双波长滤波红外测温技术的原理为：利用黑体辐射曲线中相邻两个波长对应的能量等比吸收的原理，在保证红外测温高精度测量的基础上，克服了环境对物体发射红外线吸收造成的测量误差。根据普朗克辐射定律，一个绝对温度为T的黑体，单位表面积在波长λ1、λ2(λ1的附近单位波长)的间隔内向整个半球空间发射的辐射功率(简称光谱辐射度)E0(λ，T)，满足下式的变化关系：
其中，c为真空光速c＝2.99792458×108m/s；
h为普朗克常数，h＝6.62607004×10-34J·s；
k为玻尔兹曼常数，k＝1.3806488×10-23J/K；
C1为第一辐射常数，C1＝2πhc2＝3.741771790075259×10-16W·m2；
C2为第二辐射常数，C2＝hc/k＝1.4387770620391×10-2m·K。
而灰体的光谱辐射能量公式：
其中，E0(λ，T)为黑体发射的光谱辐射通量密度，λ为光谱辐射时的波长，T为黑体的绝对温度，单位为K，ε(λ，T)为待测对象温度为T，辐射波长为λ时的辐射率，0＜ε(λ，T)≤1。
在经典近似情况下，在红外测温仪的测温范围内满足此条件，则E0(λ，T)可近似简化为Wien公式：
若波长为定值，则上式只与温度有关，可改写为：
E0(T)＝A0exp(B0/T)
其中，A0＝C1λ-5，B0＝-C2/λ，且仍然只适用于黑体。若将A0和B0视为可变参量A和B，则可推广到灰体的情况，则灰体的光谱辐射能量为：
E(T)＝A exp(B/T)。
不同于公式的是，前者只需简单改变A和B参量的值就可实现从黑体到灰体的修正，而无需确定复杂的辐射率函数ε(λ，T)。
利用相邻两波长等比吸收的原理。取2个波长的吸收能量比作为温度的函数，即可避免由于水蒸气等环境对红外线吸收的因素而造成的测量误差。
现在分别取λ1和λ2，则有：
E1(T)＝A1(λ1)exp(B1(λ1)/T)，
E2(T)＝A2(λ2)exp(B2(λ2)/T)。
上述两式作比值得：
其中，B′＝B1(λ1)-B2(λ2)。因此，只要利用实验数据拟合确定A′和B′两个参数，即可获得被测物体的温度T与该比值X之间的关系。即，根据上述拟合系数A′和B′，可以得到此种环境下的该辐射体的温度。
图1(a)示出了现有技术中一种双波长滤波红外测温设备的示意图。图1(b)示出了图1(a)中调制盘的结构示意图。参考图1(a)和图1(b)，该双波长滤波红外测温设备的工作方法是：
待测对象发出的一束光辐射沿水平方向通过透镜9射向反射镜8，由反射镜8将光辐射反射至分光镜1(或二向色镜)。分光镜1将这束光辐射进行反射和透射，形成水平方向上的反射的第一光辐射和垂直方向上的透射的第二光辐射。水平方向上的第一光辐射经过窄带滤光片7滤为波长为第一波长的光(例如波长为λ1的光)。波长为λ1的光通过反射镜6反射为垂直方向后并射向带有电机的调制盘5。垂直方向上的第二光辐射经过反射镜2反射后形成水平方向上的第二光辐射，经过窄带滤光片3滤为波长为第二波长的光(例如波长为λ2的光)。波长为λ2的光射向带有电机的调制盘5。电机带动调制盘进行转动，波长为λ2的光可以通过调制盘5上的通孔(参见图1(b))射向光敏传感器4，波长为λ1的光可以通过调制盘5上的镜面反射后射向光敏传感器4。光敏传感器4获取波长为λ1和λ2的光的能量，再经过放大电路、计算电路进行数据处理生成待测对象的温度，该温度在显示的设备上进行显示。
本申请的发明人做了大量的实验，发现利用该双波长滤波红外测温设备测得的温度的精度尽管相对于单波长红外测温设备有了明显的提高，但仍然存在一定的误差。本发明人还发现：一方面，双波长光辐射测温的精度与光辐射分成的双波的平均度成正比。即：当入射的光辐射分成的两束光束的光子数或者能量越平均，其测温的精度越高。另一方面，该双波长滤波红外测温设备测得的温度的精度还受能量损失制约，如果光辐射在测量过程中能量损失越大，测得的温度的精度越低。本发明人根据上述发现的理论，再对照现有技术中双波长滤波红外测温设备发现：
分光镜(或者二向色镜)的反射或透射效率都不高，存在较大的能量损失，所以造成后期测温精度下降；
分光镜透射和反射的第一光辐射和第二光辐射的光子数或者能量分配比例的区间大约为，该比例离理想中的1∶1的绝对等分的差距很大，所以造成后期测温精度下降。此外，双波长滤波红外测温设备中的二向色镜依赖在光学平片的正反面镀上不同的膜以实现滤光，如果需要选择其它的波长，则需要更换整块二向色镜。所以双波长滤波红外测温设备测得的温度的精度不高，且使用不太方便，其适应性不广。
尽管近年来科学家也提出了多波长测温方案，用于提高测温精度，但其结构复杂，太多的波段也会导致发射率方程组的病态程度加深，不适合实际应用。
此外，红外热成像技术通常是运用光电技术检测物体辐射的特定波段的红外线光辐射，并在检测到的红外线光辐射能量与物体的表面温度之间建立相应关系，进而获得物体的红外热像图。这种红外热像图与物体表面的热分布场相对应，热图像上的不同颜色代表被测物体的不同区域具有不同的温度。
传统的红外热成像技术通常需要面阵探测器，其成像灵敏度受限于面阵探测器对于单位像素的探测灵敏度，且面阵探测带来了测量维度和测量数上的冗余，限制了其在暗场环境下的应用。
此外，现有的红外热成像技术仅能够生成被测物体二维红外热图像，而非三维红外热图像，因此，红外热成像技术仅能够使人们可以看到被测物体表面的温度分布状况。
发明内容
鉴于现有技术中存在的问题，本发明提供一种基于压缩感知的双波长三维温度场成像设备、系统以及方法，其主要将压缩感知理论、双波长测温技术、互补调制技术、单光子探测技术、时间相关技术相结合。
根据本发明的一个方面，提供了一种基于压缩感知的双波长三维温度场成像设备。该设备包括：
脉冲光源，配置为向待测对象发射脉冲光；
光辐射调制装置，配置为接收经由待测对象的反射的光辐射，且加载预设的多个掩膜，将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿不同于第一路径的第二路径射出，所述多个掩膜根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成；
布置在所述第一路径上的第一滤光元件，配置为接收多束所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长λ1的多束光；
布置在所述第二路径上的第二滤光元件，配置为接收多束所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长λ2的多束光；
布置在第一路径上的第一探测装置，配置为接收波长为第一波长λ1的多束光并将其转换为第一光电信号参量；
布置在第二路径上的第二探测装置，配置为接收波长为第二波长λ2的多束光并将其转换为第二光电信号参量；
布置在第一路径或者第二路径上的时幅变换器，配置为记录所述第一波长λ1的光或者第二波长λ2的光中的多束光子团到达第一探测装置或第二探测装置的一系列时间戳信息；
所述第一探测装置将波长为第一波长λ1的每一束光转换为多个与所述时间戳信息对应的第一光电信号参量；
所述第二探测装置将波长为第二波长λ2的每一束光转换为多个与所述时间戳信息对应的第二光电信号参量；
温度确定装置，配置为接收来自所述第一探测装置和第二探测装置的多个与所述时间戳信息对应的所述第一和第二光电信号参量，并根据多个与所述时间戳信息对应的所述第一和第二光电信号参量与温度的预定关系确定出所述待测对象在每一个所述时间戳信息时的每一个像素点的温度值；
图像生成装置，配置为根据所述待测对象在每一个所述时间戳信息时的每一个像素点的温度值以及待测对象在每一个所述时间戳信息时的二维图像反演出待测对象的每一个横截层面的二维红外热图像，且根据待测对象的每一个横截层面的二维红外热图像生成三维红外热图像。
在本发明的一些实施方式中，所述记录所述第一波长λ1的光或者第二波长λ2的光中的多束光子团到达第一探测装置或第二探测装置的一系列时间戳信息，包括：
利用脉冲光源照射待测对象，利用所述光辐射调制装置将经所述待测对象反射的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿第二路径射出；
在所述第一路径上接收所述第一光辐射，并将其过滤为波长为第一波长的单波光；
在所述第二路径上接收所述第二光辐射，并将其过滤为波长为第二波长的单波光；
在所述第一路径上、接收所述波长为第一波长λ1的单波光或者在所述第二路径上、接收所述波长为第二波长λ2的单波光，并记录所述第一波长λ1的单波光或者所述第二波长λ2的单波光中的第一聚集光子团至第r(r≥1)聚集光子团到达第一探测装置或第二探测装置的时刻{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}，其中，r为自然数。
在本发明的一些实施方式中，所述光辐射调制装置加载预设的多个掩膜，所述第一探测装置将波长为第一波长λ1的每一束光转换为多个与所述时间戳信息对应的第一光电信号参量，所述第二探测装置将波长为第二波长λ2的每一束光转换为多个与所述时间戳信息对应的第二光电信号参量，包括：
提供2k阶的Hk矩阵，其中2k-1≥N，N为所述待测对象的总像素数；
删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ；
将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，获得0-1分布的矩阵
以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵得到矩阵H+；
提取上述矩阵H+中M行，并保留M行中的N列作为随机测量矩阵Φ；
所述光辐射调制装置加载由H+矩阵的第i行H+i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
所述第一探测装置将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i-1[(tr-1+tr)/2]，所述第二探测装置将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i-1[(tr-1+tr)/2]；
所述光辐射调制装置加载由H+矩阵的第i行的补集H-i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
所述第一探测装置将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i[(tr-1+tr)/2]，所述第二探测装置将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i[(tr-1+tr)/2]；
其中，i＝1，2，3，4……M；M/N×100％＝5％～30％。
在本发明的一些实施方式中，所述温度确定装置接收来自所述第一探测装置和第二探测装置的多个与所述时间戳信息对应的所述第一和第二光电信号参量，并根据多个与所述时间戳信息对应的所述第一和第二光电信号参量与温度的预定关系确定出所述待测对象在每一个所述时间戳信息时的每一个像素点的温度值，包括：
所述温度确定装置根据多束波长为λ1的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第一光电信号参量E1(T)i[(tr-1+tr)/2]、结合数学模型Y1 i＝ΦS1[(tr-1+tr)/2]以及Y1 i＝E1(T)2i-1[(tr-1+tr)/2]-E1(T)2i[(tr-1+tr)/2]计算出待测对象在时间戳为[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]中的部分信号值Xi(λ1)’[(tr-1+tr)/2]；
所述温度确定装置利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1下的的全部信号值Xi(λ1)[(tr-1+tr)/2]，重构待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]；
所述温度确定装置根据根据多束波长为λ2的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第二光电信号参量E2(T)i[(tr-1+tr)/2]、结合数学模型以及计算出待测对象在时间戳为[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]中的部分信号值Xi(λ2)’[(tr-1+tr)/2]；
所述温度确定装置利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测对象在时间戳信息[(tr-1+tr)/2]、波长λ2下的的全部信号值Xi(λ2)[(tr-1+tr)/2]，重构待测对象在时间戳信息[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]；
所述温度确定装置根据待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1和波长λ2下的二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]、结合Xi[(tr-1+tr)/2]＝S1[(tr-1+tr)/2]./S2[(tr-1+tr)/2]以及Ti[(tr-1+tr)/2]＝B’/ln(Xi[(tr-1+tr)/2]/A’)确定出待测对象每个像素点的温度值Ti[(tr-1+tr)/2]，其中，A’、B’为预设系数。
在本发明的一些实施方式中，所述稀疏基Ψ采用梯度基、离散余弦变换基、小波基、傅里叶变换基、Gabor变换基中的任意一种实现。
在本发明的一些实施方式中，所述待测物体在波长λ1和波长λ2下的二维图像S1和S2的重构可以采用MP、OMP、StOMP、ReOMP、SPARCO、LARS、BCS、FBMP、SL0、PPPA、KF-CS、BPDQ、k-t FOCUSS、l1-Magic、IST、TwIST、FISTA、SOCP、RecPF、FPC、GPSR、SpaRSA、ALM、FALM、TVAL3、CoSaMP、LBI、SP、l1_ls、SMP、TFOCS、NESTA、SALSA、SPGL1、YALL1、l0重建算法、l1重建算法、l2重建算法、BPDN算法、BPDN的LAG松弛算法、LASSO算法、TV算法中的任意一种实现。
在本发明的一些实施方式中，所述光辐射调制装置包括：
空间光调制器，配置为根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成的多个掩膜，以将接收到的待测对象的光辐射调制为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出；
控制元件，配置为控制所述空间光调制器依次加载由所述测量矩阵Φ变换生成的多个掩膜。
在本发明的一些实施方式中，所述空间光调制器为选自数字微镜器件、光强数字调制器或液晶光阀。
在本发明的一些实施方式中，所述第一探测装置为第一点探测器，所述第二探测装置为第二点探测器，并且
所述基于压缩感知的双波长三维温度场成像设备还包括布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一会聚元件，和布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二会聚元件，
所述第一点探测器位于所述第一会聚元件的光焦点处；
所述第二点探测器位于所述第二会聚元件的光焦点处。
在本发明的一些实施方式中，所述的基于压缩感知的双波三维长温度场成像设备还可以包括：
布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一光强衰减元件，和
布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二光强衰减元件。
在本发明的一些实施方式中，所述第一滤光元件和所述第二滤光元件为中心波长相差至少10nm的第一窄带滤光片和第二窄带滤光片，所述第一窄带滤光片和第二窄带滤光片的半高宽参数至少为10nm。
本实施方式通过设置光线会聚收集元件可以会聚、收集光线，以方便供后续的探测元件探测。即，通过高通量测量获得超过所采用的探测装置本身的灵敏度的极限。因此，可以采用两个点探测器便完成了原本需要两个面阵测温器件才能完成的工作，极大地降低了测量维度，大幅增加光通量，避免了红外光通量在维度上的分配，而将噪声压制在单像素水平，信噪比大幅提高。此外，利用不同滤光元件(例如窄带滤光片)进行滤波得到单波长光，再将该单波长光会聚收集成以供点探测器检测，降低了无关光的干扰，提高了数据采集精度。再者，设置光强衰减元件(例如中性密度滤光片)可以衰减光线的光电信号，以防止强光对后续点探测器的损伤而降低采集的光线的精度，进而提高了后续的测温精度，同时延长了产品的使用寿命。
在本发明的一些实施方式中，所述光电信号参量包括光子数。
在本发明的一些实施方式中，所述光辐射为红外波段的光辐射。
在本发明的一些实施方式中，所述第一和第二点探测器选自雪崩二极管或光电倍增管或超导单光子探测器。
本实施方式提供的点探测器可以自由选择各种类型以满足各种需求，增强了产品的通用性能，也便于后期的维护保养。本实施方式的成像设备适用于强热辐射和弱热辐射条件，且在弱热辐射条件下可达到单光子水平，其将传统光学成像推广到温度场成像的范畴，充分利用了高通量测量来获得超过所采用的探测器件本身的灵敏度极限的系统成像灵敏度(即超灵敏)，并继承了测量维度减少的优势。
本发明还提供了一种基于压缩感知的双波长三维温度场成像系统，包括：
如上所述的基于压缩感知的双波长三维温度场成像设备以及定标装置，
所述定标装置包括温度可调节的基准光源，所述基准光源配置为在一个定标阶段中，将调节成的不同温度的光辐射射向所述光辐射调制装置，由所述光辐射调制装置将接收到的光辐射等分为第一光辐射和第二光辐射，使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出，以获取多个第一和第二光电信号参量，并根据所述不同温度与所述第一和第二光电信号参量的计量值，确定出所述第一和第二光电信号参量与温度的预定关系。
本发明的实施方式中，在实际测量被测物体之前，可以先通过定标装置以及红外辐射温度曲线公式，拟合得出此种环境下的辐射体的相应系数，然后通过计算机程序设定此系数下的计算公式，即可测量或者监测任何情况下的辐射温度；有效地克服了红外测温中各种物体的“辐射率修正”难题，克服了测量条件复杂、现场测量条件波动或者水蒸气等因素的环境吸收造成的测量误差，实现了高精度温度测量。
在本发明的一些实施方式中，所述确定调节的不同温度和获取的多个第一和第二光电信号参量的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，i为1至n的自然数；A’、B’为预设系数，Tdi为所述基准光源发出的第i个温度，Ed1(T)i为第i次定标时测得的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标时测得的波长为第二波长的光的光电信号参量。
在本发明的一些实施方式中，所述定标装置还包括将基准光源的光辐射转化为平行光辐射的扩束准直透镜。
在本发明的一些实施方式中，所述定标装置还包括将所述扩束准直透镜转化的平行光辐射射向所述光辐射等分装置的分束器。
另外，本发明还提供了一种基于压缩感知的双波长温度场三维成像方法，包括：
以脉冲光照射待测对象，记录经待测对象反射的脉冲光的第一束聚集光子团至第r(r≥1)束聚集光子团到达同一地点的时刻，获取时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}；
以脉冲光源照射待测对象，利用光辐射调制装置接收经待测对象反射的光辐射，加载预设的多个掩膜，而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿第二路径射出，所述多个掩膜由符合压缩感知RIP条件的随机测量矩阵Φ变换生成；
在所述第一路径上接收多束所述第一光辐射，并将其过滤为波长为第一波长λ1的多束光；
在所述第二路径上接收多束所述第二光辐射，并将其过滤为波长为第二波长λ2的多束光；
在所述第一路径上接收多束所述波长为λ1的单波光，并将每一束波长为λ1的单波光转换为多个与时间戳信息对应的第一光电信号参量E1(T)i[(tr-1+tr)/2]；
在所述第二路径上接收多束所述波长为λ2的单波光，并将每一束波长为λ2的单波光转换为多个与时间戳信息对应的第二光电信号参量E2(T)i[(tr-1+tr)/2]；
根据多个与时间戳信息对应的所述第一光电信号参量E1(T)i[(tr-1+tr)/2]和多个与时间戳信息对应的所述第二光电信号参量E2(T)i[(tr-1+tr)/2]与温度的预定关系确定出所述待测对象在每一个时间戳[(tr-1+tr)/2]下的每一个像素点的温度Ti[(tr-1+tr)/2]；
根据所述待测对象在每一个时间戳[(tr-1+tr)/2]下的二维图像S[(tr-1+tr)/2]以及该时间戳下每一个像素点的温度值Ti[(tr-1+tr)/2]反演出所述待测对象的每一个横截层面的二维红外热图像；
根据所述待测对象的每一个横截层面的二维红外热图像反演出所述待测对象的三维红外热图像。
在本发明的一些实施方式中，所述以脉冲光照射待测对象，记录经待测对象反射的脉冲光的第一束聚集光子团至第r(r≥1)束聚集光子团到达同一地点的时刻，获取时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}包括：
利用脉冲光源照射待测对象，利用光辐射调制装置将经待测对象反射的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿第二路径射出；
在所述第一路径上接收所述第一光辐射，并将其过滤为波长为第一波长的单波光；
在所述第二路径上接收所述第二光辐射，并将其过滤为波长为第二波长的单波光；
在所述第一路径上接收所述波长为第一波长的单波光或者在所述第二路径上接收波长为第二波长的单波光，并记录探测到所述第一波长的单波光或者所述第二波长的单波光中包括的第一聚集光子团至第r聚集光子团的时刻作为时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}，其中，r为自然数。
在本发明的一些实施方式中，所述光辐射调制装置加载预设的多个掩膜，而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿第二路径射出，所述多个掩膜由符合压缩感知RIP条件的随机测量矩阵Φ变换生成；在所述第一路径上接收多束所述波长为λ1的单波光，并将每一束波长为λ1的单波光转换为多个与时间戳信息对应的第一光电信号参量E1(T)i[(tr-1+tr)/2]；在所述第二路径上接收多束所述波长为λ2的单波光，并将每一束波长为λ2的单波光转换为多个与时间戳信息对应的第二光电信号参量E2(T)i[(tr-1+tr)/2]包括：
提供2k阶的Hk矩阵，其中2k-1≥N，N为所述待测对象的总像素数；
删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ；
将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，获得0-1分布的矩阵
以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵得到矩阵H+；
提取上述矩阵H+中M行，并保留M行中的N列作为随机测量矩阵Φ；
所述光辐射调制装置加载由H+矩阵的第i行H+i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
在所述第一路径上接收多束所述波长为λ1的单波光，将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i-1[(tr-1+tr)/2]，在所述第二路径上接收多束所述波长为λ2的单波光，将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i-1[(tr-1+tr)/2]；
所述光辐射调制装置加载由H+矩阵的第i行的补集H-i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
在所述第一路径上接收多束所述波长为λ1的单波光，将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i[(tr-1+tr)/2]，在所述第二路径上接收多束所述波长为λ2的单波光，将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i[(tr-1+tr)/2]；
其中，i＝1，2，3，4……M；M/N×100％＝5％～30％。
在本发明的一些实施方式中，所述根据多个与时间戳信息对应的所述第一光电信号参量E1(T)i[(tr-1+tr)/2]和多个与时间戳信息对应的所述第二光电信号参量E2(T)i[(tr-1+tr)/2]与温度的预定关系确定出所述待测对象在每一个时间戳[(tr-1+tr)/2]下的每一个像素点的温度Ti[(tr-1+tr)/2]包括：
根据多束波长为λ1的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第一光电信号参量E1(T)i[(tr-1+tr)/2]、结合数学模型Y1 i＝ΦS1[(tr-1+tr)/2]以及Y1 i＝E1(T)2i-1[(tr-1+tr)/2]-E1(T)2i[(tr-1+tr)/2]计算出待测对象在时间戳为[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]中的部分信号值Xi(λ1)’[(tr-1+tr)/2]；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1下的的全部信号值Xi(λ1)[(tr-1+tr)/2]，重构待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]；
根据根据多束波长为λ2的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第二光电信号参量E2(T)i[(tr-1+tr)/2]、结合数学模型以及计算出待测对象在时间戳为[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]中的部分信号值Xi(λ2)’[(tr-1+tr)/2]；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测对象在时间戳信息[(tr-1+tr)/2]、波长λ2下的的全部信号值Xi(λ2)[(tr-1+tr)/2]，重构待测对象在时间戳信息[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]；
根据待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1和波长λ2下的二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]、结合Xi[(tr-1+tr)/2]＝S1[(tr-1+tr)/2]./S2[(tr-1+tr)/2]以及Ti[(tr-1+tr)/2]＝B’/ln(Xi[(tr-1+tr)/2]/A’)确定出待测对象每个像素点的温度值Ti[(tr-1+tr)/2]，其中，A’、B’为预设系数。
在本发明的一些实施方式中，所述稀疏基Ψ采用梯度基、离散余弦变换基、小波基、傅里叶变换基、Gabor变换基中的任意一种实现。
在本发明的一些实施方式中，所述待测物体在波长λ1和波长λ2下的二维图像S1和S2的重构可以采用MP、OMP、StOMP、ReOMP、SPARCO、LARS、BCS、FBMP、SL0、PPPA、KF-CS、BPDQ、k-t FOCUSS、l1-Magic、IST、TwIST、FISTA、SOCP、RecPF、FPC、GPSR、SpaRSA、ALM、FALM、TVAL3、CoSaMP、LBI、SP、l1_ls、SMP、TFOCS、NESTA、SALSA、SPGL1、YALL1、l0重建算法、l1重建算法、l2重建算法、BPDN算法、BPDN的LAG松弛算法、LASSO算法、TV算法中的任意一种实现。
在本发明的一些实施方式中，所述光辐射调制装置包括：
空间光调制器，配置为根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成的多个掩膜，以将接收到的待测对象的光辐射调制为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出；
控制元件，配置为控制所述空间光调制器依次加载由所述测量矩阵Φ变换生成的多个掩膜。
在本发明的一些实施方式中，所述空间光调制器选自数字微镜器件、光强数字调制器或液晶光阀。
在本发明的一些实施方式中，所述的基于压缩感知的双波长三维温度场成像方法，进一步包括：
将所述波长为第一波长的光会聚于第一焦点，并在所述第一焦点处设置第一点探测装置用以接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量；
将所述波长为第二波长的光会聚于第二焦点，并在所述第二焦点处设置第二点探测装置用以接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量。
在本发明的一些实施方式中，所述的基于压缩感知的双波长三维温度场成像方法，进一步包括：
对所述波长为第一波长的光的强度进行衰减；及
对所述波长为第二波长的光的强度进行衰减。
在本发明的一些实施方式中，所述光辐射为红外波段的光辐射。
在本发明的一些实施方式中，所述第一波长与所述第二波长相差至少10nm。
在本发明的一些实施方式中，所述光电信号参量包括光子数。
在本发明的一些实施方式中，所述以脉冲光照射待测对象，记录经待测对象反射的脉冲光的第一束聚集光子团至第r(r≥1)束聚集光子团到达同一地点的时刻，获取时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}之前还包括定标步骤，
所述定标步骤包括：
利用光辐射调制装置接收基准光源的光辐射，将接收到的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出；
在所述第一路径上接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光；
在所述第二路径上接收所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光；
在所述第一路径上接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量；
在所述第二路径上接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量；
调节基准光源发出多个不同温度的光辐射，并获取对应的多个第一和第二光电信号参量；以及
根据所述不同温度与所述第一和第二光电信号参量的计量值，确定出所述第一和第二光电信号参量与温度的预定关系。
在本发明的一些实施方式中，所述确定调节的不同温度与获取的多个第一和第二光电信号参量的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量
本发明将压缩感知理论、双波长测温技术、互补调制技术、单光子探测技术和时间相关技术相结合，提出基于压缩感知的超灵敏双波长三维温度场成像设备、系统以及方法，以克服现有技术中接触式测温应用场合窄、单波长测温辐射率难以修正、误差大、多波长测温结构复杂、传统双波长测温灵敏度低、测量维度高、测量数冗余、扩展性差等缺陷。利用本发明能够重建被测对象的三维红外热图像。该图像与物体表面的热分布场相对应，即反映被测对象的红外辐射能量分布，三维红外热图像上显示的不同颜色代表被测对象的不同温度，可以观察到被测物体的整体温度分布状况，研究被测物体的发热情况。
本发明利用压缩感知原理实现亚采样和低维度测量(采样率5％～30％)，具有很好的去噪能力。并且将传统压缩光学成像推广到温度场成像的范畴。此外，本发明充分利用了高通量测量来获得超过所采用的探测器件本身的灵敏度极限的灵敏度(即超灵敏)，并继承了测量维度和测量数均减少的优势。
本发明还以Hadamard矩阵、排列矩阵、DMD互补调制相结合的方案实现了正负光强调制，使得测量矩阵满足受限等距性质(Restricted Isometry Property，RIP)条件，适合于CS的重建计算，使图像质量获得1～2个数量级的提高，其将成为测温领域的一个重要发展方向。该项技术可以广泛应用在深空探测、遥感、材料检测、夜视观测等相关科技领域。
附图说明
图1(a)为现有技术中一种双波长测温设备的结构示意图；
图1(b)为图1(a)中调制盘的结构示意图；
图2为本发明一些实施方式的基于光辐射的测温设备的结构示意图；
图3(a)为本发明一些实施方式的DMD中多个微镜结构示意图；
图3(b)为图3(a)中两片微镜结构示意图；
图4为本发明另一些实施方式的基于光辐射的测温设备的结构示意图；
图5为本发明一些实施方式的基于光辐射的测温系统的结构示意图；
图6为本发明一些实施方式的基于光辐射的测温方法的流程示意图；
图7为本发明一些实施方式的基于光辐射的测温方法的定标步骤的流程示意图；
图8为本发明一些实施方式的基于压缩感知的双波长温度场成像设备的结构示意图；
图9为8×8的Hadamard矩阵；
图10为本发明另一些实施方式的基于压缩感知的双波长温度场成像设备的结构示意图；
图11为本发明一些实施方式的基于压缩感知的双波长温度场成像系统的结构示意图；
图12为本发明一些实施方式的基于压缩感知的双波长温度场成像方法的流程示意图；
图13为本发明一些实施方式的基于压缩感知的双波长温度场成像方法的定标步骤流程示意图；
图14为本发明一些实施方式的基于压缩感知的双波长三维温度场成像设备的结构示意图；
图15为本发明另一些实施方式的基于压缩感知的双波长三维温度场成像设备的结构示意图；
图16为本发明一些实施方式的基于压缩感知的双波长三维温度场成像系统的结构示意图；
图17为本发明一些实施方式的基于压缩感知的双波长三维温度场成像方法的流程示意图；
图18为本发明一些实施方式的基于压缩感知的双波长三维温度场成像方法的定标步骤流程示意图。
具体实施方式
为了使发明的目的、技术方案和优点更加清楚，下面结合附图和具体实施例对发明作进一步详细的说明。虽然附图中显示了本公开示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻的理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
参照图2，本发明的一些实施方式提供了一种基于光辐射的测温设备的结构。该测温设备包括：光辐射等分装置2、第一滤光元件3-1、第二滤光元件4-1、第一探测装置3-4、第二探测装置4-4，以及分别与第一探测装置3-4、第二探测装置4-4连接的温度确定装置5。该测温设备的各个部件的连接关系和处理光辐射的过程可以描述如下：
待测对象(例如人体或者电灯，图中未进行标示)向光辐射等分装置2发出光辐射(例如红外辐射、紫外辐射或者可见光)。光辐射等分装置2接收待测对象的光辐射，将接收到的光辐射等分为第一光辐射和第二光辐射，并使第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出。布置在第一路径上的第一滤光元件3-1接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光(例如波长为λ1的单波光)。布置在第二路径上的第二滤光元件4-1接收第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光(例如波长为λ2的单波光，其中λ1与λ2不等，当λ1与λ2为相邻波段的波长时，效果最佳，因为当λ1与λ2无限接近时，ε1(λ1，T)≈ε2(λ2，T)，则：
A′＝A1(λ1)/A2(λ2)＝(ε1(λ1，T)C1λ1 -5)/(ε2(λ2，T)C1λ2 -5)≈(λ1/λ2)-5，
B′＝B1(λ1)-B2(λ2)＝-C2/λ1+-C2/λ2，
但是，λ1与λ2越接近，对探测器的灵敏度和准确性的要求越高，另外考虑到环境噪声和探测器本征噪声(如暗计数等)对系统的测温准确性造成的较大影响，λ1与λ2在实际测量中不可能无限接近，因而两个辐射波长下的辐射率ε1(λ1，T)和ε2(λ2，T)不可简单相消，而需要通过定标确定出预定关系。布置在第一路径上的第一探测装置3-4接收波长为第一波长的光并将其转换为相应的第一光电信号参量。布置在第二路径上的第二探测装置4-4接收波长为第二波长的光并将其转换为相应的第二光电信号参量。温度确定装置5接收来自第一探测装置3-4和第二探测装置4-4的第一和第二光电信号参量，并根据所述第一和第二光电信号参量与所述待测对象的温度的预定关系确定出待测对象的温度。
在本实施例中，光辐射等分的含义是：将接收到的光辐射的光子数或者能量按的比例区间进行分配。在本实施方式中，可以达到1∶1的平均分配的比例，此时，测温的效果最佳。由此可知，本实施例的等分的比例远远比现有技术中二向色镜中光辐射分配的比例的平均度要高很多。由于测温的精度随着平均度的提高而提高，所以，本实施例测温精度比现有技术的精度要高很多。
继续参见图2，温度确定装置5包括除法器5-1和计算元件5-2。其中，除法器5-1分别与第一探测装置3-4和第二探测装置4-4连接。计算元件5-2与除法器5-1连接。除法器5-1用于计算波长为第一波长的光的光电信号参量E1(T)和波长为第二波长的光的光电信号参量E2(T)之间的比值X。计算元件5-2用于根据第一和第二光电信号参量E1(T)和E2(T)与所述待测对象的温度的预定关系确定出待测对象的温度。
由此，本设备为运算频率较高、运算重要程度高的比值运算单独提供除法器，使得比值运算与其它逻辑运算区别开来，优化了测温设备的结构，减少了运算错误、缩短了运算的时间，提高了运算精度。此外，这种硬件结构也可以以软件模块的方式来实现。
在本实施例中，温度确定装置根据如下的预定关系确定待测对象的温度：
T＝B’/ln(X/A’)
其中，X＝E1(T)/E2(T)，A’、B’为预设系数，T为待测对象的温度，E1(T)为第一光电信号参量，E2(T)为第二光电信号参量。
再次参见图2，光辐射等分装置包括：空间光调制器2-1和控制元件2-2。其中，空间光调制器2-1根据预定控制将接收到的待测对象的光辐射等分为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出。控制元件2-2对所述空间光调制器进行预定控制(例如通过加载掩膜(数字图像处理中，掩模为二维矩阵数组)为0-1矩阵中0和1的数量来控制，该部分内容下文还会继续介绍)。
在本实施例中，空间光调制器可以选用数字微镜器件(Digital MicromirrorDevice，DMD)、光强数字调制器或液晶光阀中的任意一种。因为DMD)、光强数字调制器或液晶光阀均是现有的产品，因此，下面仅对DMD等量分配光辐射进行详细说明，其余产品不再赘述。
图3(a)示出了本发明一种实施方式的DMD中多个微镜结构示意图。图3(b)示出了图3(a)中两片微镜结构示意图。
参考如图3(a)和图3(b)，DMD包括多个微镜和与所述多个微镜对应的多个转动铰链，各个转动铰链可以根据预定控制将各个微镜向预设方向(例如与垂直方向成+12度和-12度)翻转，使得所述多个微镜中的一半微镜将接收到的待测对象的一半的光辐射沿第一路径射出，所述多个微镜中的另一半微镜将接收到的待测对象的另一半的光辐射沿不同于第一路径的第二路径射出。
在本实施例中，DMD可以选用市场上能够获得的TI(美国德州仪器)公司生产的型号为0.7XGA2XLVDS DMD的器件。控制元件可以选用FPGA(可编程门阵列芯片)，FPGA型号为：Xilinx Virtex5application FPGA。
FPGA控制DMD进行光辐射按预设比例的分配的原理如下：
DMD加电后，利用静电吸附原理，DMD中多个微镜可以呈现+12度和-12度(也有的是+10度和-10度)方向的偏转。假设DMD有1000片的微镜，当500片微镜翻转为+12度，500片微镜翻转为-12度时，此时射向DMD的一束光辐射就会反射为夹角为48度(12度×2+12度×2)的等分的两束光辐射。具体FPGA控制DMD中多少个的微镜分别向+12度和-12度的偏转，可以通过加载掩膜(数字图像处理中，掩模为二维矩阵数组)为0-1矩阵中0和1的数量来控制。例如加载0时，微镜翻转为+12度，加载1时，微镜翻转为-12度。所以想要控制两束光辐射的具体比例(可以时任何比例，例如等分的50％：50％，或者20％：80％)，只需控制矩阵中0和1的比例即可。因此，通过FPGA控制DMD可以实现将光辐射的光子数(或者光强)进行1∶1等分。由于DMD是现有产品，所以其更具体结构在此不再赘述。
图4示出了本发明另一些实施方式的基于光辐射的测温设备的结构示意图。该图示出了几种变形的实施方式。图4所示的实施例是在上述图2所示的实施例的基础之进行变形得出的。在此着重描述二者不同之处，二者相同或者相似之处不再赘述。
图4所示的基于光辐射的测温设备的第二种实施例为：
该实施例是在图2所示的第一实施例的基础上增加第一会聚元件3-3和第二会聚元件4-3。该实施例的测温设备的各部件之间的连接关系可以为：
第一探测装置3-4为第一点探测器，第二探测装置4-4为第二点探测器。第一会聚元件3-3布置在所述第一路径上、位于所述第一点探测器与所述第一滤光元件之间。第二会聚元件4-3布置在所述第二路径上、位于所述第二点探测器与所述第二滤光元件之间。第一点探测器3-4位于所述第一会聚元件的光焦点处。第二点探测器4-4位于第二会聚元件的光焦点处。
本实施例通过将点探测器布置在反射光线聚焦形成的焦点处，而不是布置在传统的焦平面后方的像平面(成像平面)处，不仅可以方便利用点探测器检测光电信号参量，而且可以增加点探测器检测到的光子数，极大地增大了信号的强度，而将散粒噪声降低到单个像素的水平，大幅度提高了测量信噪比，提高了采集数据的精确程度，进而大幅度提高了测量温度的灵敏度和精确度。本设备结构简单且可以用于不同要求测温场合，适用范围广。
图4所示的基于光辐射的测温设备的第三种实施例为：
该实施例是在上述第二实施例的基础上增加了布置在第一路径上、第一滤光元件3-1和第一会聚元件3-3之间的第一光强衰减元件3-2，和布置在第二路径上、第二滤光元件4-1和第二会聚元件4-3之间的第二光强衰减元件4-2。
由此，本装置通过设置光强衰减元件(例如中性密度滤光片)可以衰减光线的光电信号参量，以防止强光对后续点探测器的损伤，提高了测温精度，延长了产品的使用寿命。
图4所示的基于光辐射的测温设备的第四种实施例为：
在上述各实施例的基础上增加透镜1，使得待测对象的光辐射进入主光路，可以射向光辐射等分装置2。
其中，第二实施例相对于第三实施例省去了第一光强衰减元件3-2和第二光强衰减元件4-2，由此，光强衰减的作用也相应消失，但相对于现有技术，第二实施例仍然可以解决技术问题，达到相应的技术效果。第一实施例相对于第二实施例省去了第一会聚收件3-3和第二会聚元件4-3。由此，会聚的作用也相应消失，但相对于现有技术，第一实施例仍然可以解决技术问题，达到相应的技术效果。本领域的技术人员可以理解，可以根据实际测量精度或者特殊需要对上述各元件进行选择配置，或者进行多种组合配置。例如，在光强特别大的情况下，可以设置多块光强衰减元件。
参照图4，下面介绍一优选实施方式的测温装置的测温的工作方式，其具体可以是：
待测对象(例如人体或者电灯，图中未进行标示)的光辐射可以经过透镜1射向光辐射等分装置2。光辐射等分装置2将接收的光辐射平均等分(将光子数和光辐射的能量平均分配)反射为双臂方向上的两束光辐射。在双臂中的左臂方向上分别布置有第一滤光元件3-1、第一光强衰减元件3-2、第一会聚元件3-3和第一探测装置3-4。第一束光辐射经过第一滤光元件3-1过滤为波长为第一波长的光，波长为第一波长的光经过第一光强衰减元件3-2进行光强衰减，再经过第一会聚元件3-3进行会聚，第一探测装置3-4在会聚的波长为第一波长的光焦点所在的焦平面处接收波长为第一波长的光来探测波长为第一波长的光的光电信号参量。
本设备左右对称设置。同样，在双臂中的右臂方向上分别布置有第二滤光元件4-1、第二光强衰减元件4-2、第二会聚元件4-3和第二探测装置4-4。第二束光辐射经过第二滤光元件4-1过滤为波长为第二波长的光，波长为第一波长的光经过第二光强衰减元件4-2进行光强衰减，再经过第二会聚元件4-3进行会聚，第二探测装置4-4在会聚的波长为第一波长的光焦点所在的焦平面处接收波长为第二波长的光来探测波长为第二波长的光的光电信号参量。
温度确定装置5分别接收第一探测装置3-4和第二探测装置4-4探测的波长为第一波长的光的光电信号参量和波长为第二波长的光的光电信号参量，并对所述待测对象进行测温。
在一些实施例中，为了使过滤后的单波的波长不同，第一滤光元件和所述第二滤光元件为中心波长相差10nm以上的第一窄带滤光片和第二窄带滤光片，所述第一窄带滤光片和第二窄带滤光片的半高宽参数为10nm以上。
因为根据窄带带宽理论，窄带滤光片越窄效果越好；中心波长越接近越好。但是，带宽越窄，与窄带滤光片配合使用的光电信号参量探测元件能探测到的热辐射越少，外加之中心波长接近，对探测器的灵敏度和准确性提出过高的要求，而且此时环境噪声和探测器本征噪声(如暗计数等)也将对系统的测温准确性造成较大影响。所以，在衡量上述优缺点的基础上，经过大量的试验，选用的窄带滤光片的半高宽FWHM一般应在10nm以上，中心波长CWL一般相差10nm以上效果最佳。另外，本设备可以利用不同滤光元件(例如窄带滤光片)进行滤波得到单波光，再将该单波光会聚收集成以供探测器检测，降低了无关光的干扰，提高了采集温度的精度。
在一些实施例中，所述光辐射为红外线光辐射。
在一些实施例中，所述第一和第二点探测器选自近红外、中远红外、远红外波段的外光电效应探测器组、内光电效应探测器组、强光探测器组和弱光探测器组中的任意一种，其中：
所述外光电效应探测器组包括：雪崩二极管、真空光电管、充气光电管、光电倍增管、变像管、像增强器、摄像管；
所述内光电效应探测器组包括：本征型光电导探测器、掺杂型光电导探测器、光磁电效应探测器、光生伏特探测器；
所述强光探测器组包括：内置或者外置有模数转换器的强光探测器；
所述弱光探测器组包括：内置或者外置有计数器的弱光探测器。
在一些实施例中，所述光电信号参量包括光子数、电流值、电压值、电阻值中任意一种。
由此，点探测器可以自由选择各种类型以满足各种需求，增强了产品的通用性能，也便于后期的维护保养。
图5示出了本发明一些实施方式的基于光辐射的测温系统的结构示意图。该图可以有几种变形的实施方式。在此着重描述这些变形方式的不同之处，它们相同或者相似之处不再赘述。
图5所示的基于光辐射的测温系统的第一种实施例为：
参照图5，该测温系统包括：上述测温设备和定标装置6。
定标装置6可以包括温度可调节的基准光源6-1(例如，可以提供不同功率的灯泡)、扩束准直透镜6-2和分束器6-3。分束器6-3设置在待测对象7和透镜1之间。基准光源6-1、扩束准直透镜6-2和分束器6-3水平共线设置，扩束准直透镜6-2设置在基准光源6-1和分束器6-3之间。其中，基准光源6-1用于在定标阶段将调节的不同温度的光辐射射向光辐射等分装置2，并由光辐射等分装置2将接收到的光辐射等分为第一光辐射和第二光辐射，使第一光辐射沿第一路径射出(如左臂方向路径)、第二光辐射沿第二路径(如右臂方向路径)射出，以获取多个第一和第二光电信号参量，并确定调节的不同温度和获取的多个第一和第二光电信号参量的预定关系。扩束准直透镜6-2用于将基准光源的光辐射转化为平行光辐射。分束器6-3用于将扩束准直透镜转化的平行光辐射射向所述光辐射等分装置。
图5所示的基于光辐射的测温系统的第二种实施例为：
在第一实施方式的基础上减少了分束器6-3，该分束器的功能随之减少。但此实施方式仍能解决技术问题，实现相应的技术效果。
图5所示的基于光辐射的测温系统的第三种实施例为：
在第二实施方式的基础上减少了扩束准直透镜6-2，该扩束准直透镜的功能随之减少。但此实施方式仍能解决技术问题，实现相应的技术效果。
在上述实施方式中，所述确定调节的不同温度和获取的多个第一和第二光电信号参量的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，所述多个为n个，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量。
由此，本系统可以采用不同功率的标准光源进行试验测温，通过测得多组数据来确定试验中的拟合系数。其中，拟合算法利用已知试验或者真实数据，然后寻找一个模型对其规律进行模拟的过程中，求取模型中未知参数的一个过程。以确保在后期的实际测温过程中根据该拟合系数获取测得的温度。本实施方式通过多次试验测进行试验测温，提高了实际测温的精度。此外，本实施方式利用扩束准直透镜6-2将标准光源的光线变换成平行光，减少了光线因为会聚或者散射引起的误差。利用分束器6-3可以最大程度的将标准光源的光线全部传递至透镜，提高了光线的传递效率。
本领域的技术人员可以理解，可以根据实际测量精度或者特殊需要对上述各元件进行选择配置，或者进行多种组合配置。
由此，本系统通过提供定标装置，可以在正式测温之前，进行试验测温，从而制定标准，使得后面的正式测温可以参考试验测温的数据进行运算，调整测得的温度数据，进一度调高了温度的精度。
本实施例的测温系统的测温的工作方式可以参照上述的测温设备的描述。需要说明的是，在定标的试验之前，将分束器6-3、扩束准直透镜6-2和标准光源2-1移入系统中，保证分束器6-3的反射方向在透镜1和DMD2-1主轴光路上，保证标准光源的光辐射进入系统的主轴光路，且标准光源与待测对象7所在位置共轭，待定标结束后再将分束器6-3、扩束准直透镜6-2和标准光源6-1移出系统。
图6为本发明一些实施方式的基于光辐射的测温方法的流程示意图。如图6所示，该方法包括以下步骤：
S601：利用光辐射等分装置(可以利用图2和图4实施方式中的光辐射等分装置)接收待测对象(例如人体或者电灯)的光辐射(例如红外辐射、紫外辐射或者可见光)并直接将接收到的光辐射等分(此处等分的概念可以参考上述测温设备中等分的概念)为第一光辐射和第二光辐射，并使第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出；
S602：在上述的左臂方向的路径上接收第一光辐射，并将其过滤为波长为第一波长的光(例如波长为λ1的单波光)；
S603：在上述的右臂方向的路径上接收第二光辐射，并将其过滤为波长为第二波长的光(例如波长为λ2的单波光)；
S604：在上述的左臂方向的路径上接收波长为λ1的单波光并将其转换为相应的第一光电信号参量E1(T)；
S605：在上述的右臂方向的路径上接收波长为λ2的单波光并将其转换为相应的第二光电信号参量E2(T)；
S606：根据所述第一和第二光电信号参量E1(T)和E2(T)与温度的预定关系确定出待测对象(例如人体或者电灯)的温度。
在本实施例中，所述温度确定装置根据如下的预定关系确定待测对象的温度：
T＝B’/ln(X/A’)
其中，X＝E1(T)/E2(T)，A’、B’为预设系数，T为待测对象的温度，E1(T)为第一光电信号参量，E2(T)为第二光电信号参量。
在本实施例中，光辐射等分装置包括：空间光调制器和控制元件。其中，空间光调制器用于根据预定控制将接收到的待测对象的光辐射等分为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出。控制元件用于对所述空间光调制器进行预定控制。
在本实施例中，所述空间光调制器选DMD、光强数字调制器或液晶光阀。
在本实施例中，DMD包括多个微镜和与所述多个微镜对应的多个转动铰链，各个转动铰链根据所述预定控制将各个微镜向预设方向翻转，使得所述多个微镜中的一半微镜将接收到的待测对象的一半的光辐射沿第一路径射出，所述多个微镜中的另一半微镜将接收到的待测对象的另一半的光辐射沿不同于第一路径的第二路径射出。
在本发明的一些实施例中，该方法还包括：
将所述波长为第一波长的光会聚于第一焦点，并在所述第一焦点处设置第一点探测装置用以接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量；
将所述波长为第二波长的光会聚于第二焦点，并在所述第二焦点处设置第二点探测装置用以接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量。
在本发明的一些实施例中，该方法还包括：
对所述波长为第一波长的光的强度进行衰减；及对所述波长为第二波长的光的强度进行衰减。
在本发明的一些实施例中，所述光辐射为红外线光辐射。
在本发明的一些实施例中，所述第一波长与所述第二波长相差至少10nm。
在本发明的一些实施例中，所述光电信号参量包括光子数、电流值、电压值、电阻值中任意一种。
在本发明的一些实施例中，所述第一和第二点探测器选自近红外、中远红外、远红外波段的外光电效应探测器组、内光电效应探测器组、强光探测器组和弱光探测器组中的任意一种，其中，
所述外光电效应探测器组包括：雪崩二极管、真空光电管、充气光电管、光电倍增管、变像管、像增强器、摄像管；
所述内光电效应探测器组包括：本征型光电导探测器、掺杂型光电导探测器、光磁电效应探测器、光生伏特探测器；
所述强光探测器组包括：内置或者外置有模数转换器的强光探测器；
所述弱光探测器组包括：内置或者外置有计数器的弱光探测器。
在测温方法上述实施例中的技术效果与测温装置实施例中的技术效果相对应，在此不再赘述。
图7为本发明一些实施方式的定标的步骤的流程示意图。在本实施方式中，在所述利用光辐射等分装置接收待测对象的光辐射，将接收到的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出的步骤之前还包括定标步骤。如图7所示，所述定标步骤包括：
S701：利用光辐射等分装置(可以利用图2和图4实施方式中的光辐射等分装置)接收基准光源(例如温度可调节的灯泡，具体可以调节电流、电压等方式来调节)的光辐射，将接收到的光辐射等分(此处等分的概念可以参考上述测温设备中等分的概念)为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、所述第二光辐射沿不同于第一路径的第二路径(如测温设备的右臂方向的路径)射出。
S702：在所述第一路径上接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光(例如波长为λ1的单波光)。
S703：在所述第二路径上接收所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光(例如波长为λ2的单波光)。
S704：在所述第一路径上接收所述波长为λ1的单波光并将其转换为相应的第一光电信号参量E1(T)。
S705：在所述第二路径上接收所述波长为λ2的单波光并将其转换为相应的第二光电信号参量E2(T)。
S706：调节灯泡的电流或电压，使得灯泡发出光辐射的温度变为Td2，并获取对应的光电信号参量Ed1(T)2和Ed2(T)2。按此种方法可以获取灯泡发出多个不同温度Tdi的光辐射，并获取对应的多个第一和第二光电信号参量Ed1(T)i和Ed2(T)i。
S707：根据上述不同温度与第一和第二光电信号参量的计量值，确定出所述第一和第二光电信号参量与温度的预定关系。
在本实施方式中，所述确定多个第一和第二光电信号参量与调节的多个不同温度的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，所述多个为n个，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量。
具体实现方式可以是：
在定标时，标准光源以第一电流、电压、电阻的光线照射扩束准直透镜，通过分束器、透镜到达DMD。保持DMD同一帧光辐射不变，其加载的掩膜为1和0数量相等的0-1矩阵。控制元件控制DMD中各微镜的闭合、翻转，使得DMD将接收到的一帧光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出。其中：第一光辐射经过第一窄带滤光片、第一会聚元件到达第一点探测器，得到第一光电信号参量Ed1(T)1。第二光辐射经过第二窄带滤光片、第二会聚元件到达第二点探测器，得到第二光电信号参量Ed2(T)1。Ed1(T)1和Ed2(T)1经过除法器得到比值X1，X1＝Ed1(T)1/Ed2(T)1；
调节标准光源的不同等效电流或电压或电阻可计算出不同电流或电压或电阻下的一系列温度T1，T2，T3……Tdi，通过本系统测得对应的Ed1(T)i与Ed2(T)i比值X1，X2，X3……Xi，保持DMD固定一帧不变，其DMD上加载的掩膜为1和0数量相等的0-1矩阵。接着利用公式拟合得出相应系数，即得到此种环境下的辐射体(即基准光源)的光辐射温度曲线公式。即：确定拟合系数(即上述的预设系数)A’和B’，其中：
Ed1(T)i/Ed2(T)1＝A’exp(B’/Td1)
Ed1(T)2/Ed2(T)2＝A’exp(B’/Td2)
Ed1(T)3/Ed2(T)3＝A’exp(B’/Td3)
……
Ed1(T)＝A1(λ1)exp(B1(λ1)/T)
Ed2(T)＝A2(λ2)exp(B2(λ2)/T)
A’＝A1(λ1)/A2(λ2)，B’＝B1(λ1)-B2(λ2)。
具体可以是利用公式拟合得出相应系数A’和B’，得到此种环境下的该辐射体的光辐射温度曲线公式；然后通过计算机程序设定好此系数下的计算公式，即可测量或者监测任何情况下的辐射温度。在实际测温过程中，利用指数形式的多项式展开，获得多项式的系数，或者利用其它复杂的拟合方程获得拟合系数。
由此，本实施方式有效地克服了红外测温中各种物体的“辐射率修正”难题，克服了测量条件复杂、现场测量条件波动或者水蒸气等因素的环境吸收造成的测量误差，实现了高精度的温度测量。
在本实施方式中，在实际测温过程中，可以通过计算机程序设定好此系数下的计算公式，即可测量或者监测任何情况下的辐射温度。可以利用指数形式的多项式展开，获得多项式的系数，或者利用其它复杂的拟合方程获得拟合系数，可在实际操作过程中实现更高精度的测量。这里实际是利用现场定标的实验方法找出实际条件下待测对象(某一具体物体，灰体而非黑体)的热辐射规律，并给与校准，再利用该规律在原有条件下进行温度监测或测温，即利用“替代法”校准该实际条件下的许多隐含参量。
本实施方式中，温度确定装置(如除法器、计算元件)可以由单片机、FPGA芯片、计算机、服务器等替换。数据运算可以通过程序模块执行。程序可以包括例程、程序、目标程序、组件、逻辑、数据结构等，它们执行特定的任务或者实现特定的抽象数据类型。计算机系统/服务器可以在分布式云计算环境中实施，分布式云计算环境中，任务是由通过通信网络链接的远程处理设备执行的。在分布式云计算环境中，程序模块可以位于包括存储设备的本地或远程计算系统存储介质上。
如图1～7所示的本发明所提供的基于光辐射的测温设备、系统和方法只能够测量待测对象(例如人体或者电灯)整体的平均温度，不能测量其每一个局部的温度并生成待测对象的二维红外热像图。
鉴于此，本发明结合压缩感知理论，提供了基于压缩感知的双波长温度场成像设备。
压缩感知(Compressed Sensing，简称CS)也被称为压缩采样(CompressiveSampling)，稀疏采样(Sparse Sampling)，压缩传感。作为一个新的采样理论，其通过开发信号的稀疏特性，在远小于100％采样率的条件下(远低于奈奎斯特/香农采样定理的极限)，用随机采样获取信号的离散样本，然后通过非线性重建算法完美的重建信号。
压缩感知的基本步骤包括：
利用先验知识，选取合适的稀疏基Ψ，使得待测物体x经Ψ变换后得到展开系数x’是最为稀疏的；
在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的条件下，建立数学模型Y＝ΦΨx’+e，通过压缩感知算法进行凸优化，获得部分信号x’后，再由反演出全部信号x。
参照图8所示，在本发明的第一种实施例中，基于压缩感知的双波长温度场成像设备，该设备与图2所示的基于光辐射的测温设备的不同之处在于：
第一、光辐射等分装置2不必定将待测对象的光辐射进行等分处理，其可以加载预先设定的矩阵变换生成的掩膜，因此，在本实施例中，称之为光辐射调制装置2’，其根据预先设定的矩阵加载掩膜的具体规则将在后文予以详细阐述。
第二、温度确定装置5确定的为待测对象每一个像素点的温度值。
第三、增加了与温度确定装置5连接的图像生成装置8。图像生成装置8用于根据温度确定装置5生成的待测对象的每一个像素点的温度值以及待测对象的二维图像反演出待测对象的二维红外热图像。
具体而言，在第一个实施例中，基于压缩感知的双波长温度场成像设备包括：光辐射调制装置2’、第一滤光元件3-1、第二滤光元件4-1、第一探测装置3-4、第二探测装置4-4，分别与第一探测装置3-4、第二探测装置4-4连接的温度确定装置5，以及与温度确定装置5连接的图像生成装置8。
待测对象(例如人体或者电灯，图中未进行标示)向光辐射调制装置2’发出光辐射(例如红外辐射、紫外辐射或者可见光)。光辐射调制装置2’接收待测对象的光辐射，并且加载预设的多个掩膜(所述多个掩膜根据符合压缩感知理论的RIP条件的测量矩阵变换生成)。随着光辐射调制装置2’加载的掩膜发生变化，其将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、多束第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出。布置在第一路径上的第一滤光元件3-1接收多束所述第一光辐射，并将接收到的多束第一光辐射过滤为波长为第一波长的多束光(例如波长为λ1的单波光)。布置在第二路径上的第二滤光元件4-1接收多束第二光辐射，并将接收到的多束第二光辐射过滤为波长为第二波长的多束光(例如波长为λ2的单波光，其中λ1与λ2不等，当λ1与λ2为相邻波段的波长时，效果最佳，因为当λ1与λ2无限接近时，ε1(λ1，T)≈ε2(λ2，T)，则：
A′＝A1(λ1)/A2(λ2)＝(ε1(λ1，T)C1λ1 -5)/(ε2(λ2，T)C1λ2 -5)≈(λ1/λ2)-5，
B′＝B1(λ1)-B2(λ2)＝-C2/λ1-(-C2/λ2)，
但是，λ1与λ2越接近，对探测器的灵敏度和准确性的要求越高，另外考虑到环境噪声和探测器本征噪声(如暗计数等)对系统的测温准确性造成的较大影响，λ1与λ2在实际测量中不可能无限接近，因而两个辐射波长下的辐射率ε1(λ1，T)和ε2(λ2，T)不可简单相消，而需要通过定标确定出预定关系。布置在第一路径上的第一探测装置3-4接收多束波长为第一波长的光并将其转换为相应的多个第一光电信号参量。布置在第二路径上的第二探测装置4-4接收多束波长为第二波长的光并将其转换为相应的多个第二光电信号参量。
温度确定装置5接收来自第一探测装置3-4和第二探测装置4-4的多个第一和多个第二光电信号参量，并根据多个所述第一和第二光电信号参量确定待测对象在波长λ1和波长λ2下局部信号值Xi(λ1)’和Xi(λ2)’，进而结合压缩感知原理获取待测对象的二维图像S1和S2。
温度确定装置5根据待测对象在波长λ1和波长λ2下的二维图像S1和S2以及所述待测对象的温度的预定关系确定出待测对象每一个像素点的温度值。
图像生成装置8用于根据温度确定装置5生成的待测对象的每一个像素点的温度值以及待测对象的二维图像反演出待测对象的二维红外热图像。
在本实施例中，温度确定装置根据如下的预定关系确定待测对象的温度：
Ti＝B’/ln(Xi/A’)
其中，Xi＝S1./S2(矩阵的定点除法运算)，A’、B’为预设系数，Ti为待测对象的每一个像素点的温度，E1(T)i为第一光电信号参量，E2(T)i为第二光电信号参量。
再次参见图8，光辐射调制装置2’包括：空间光调制器2-1和控制元件2-2。其中，空间光调制器2-1配置为根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成的多个掩膜，以将接收到的待测对象的光辐射调制为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出。控制元件2-2配置为控制所述空间光调制器依次加载由所述测量矩阵Φ变换生成的多个掩膜。例如：根据预设矩阵变换生成掩膜，最后根据掩膜中0和1的数量控制空间光调制器2-1中微镜的翻转。
其中，所采用的预设矩阵为压缩感知理论的RIP条件的测量矩阵Φ，该测量矩阵Φ由Hadamard矩阵以及扩充排列矩阵P和Q得到，该部分内容下文还会继续介绍。
在本实施例中，空间光调制器可以选用如图3(a)和图3(b)所示的数字微镜器件(Digital MicromirrorDevice，DMD)、光强数字调制器或液晶光阀中的任意一种。
下面将详细描述测量矩阵Φ的获得方式，以及光辐射调制装置2’、温度确定装置5和图像生成装置8的工作原理：
2k阶的Hadamard矩阵满足以下特点：
矩阵Hk中的元素只有1和-1(不考虑系数)，其中1和-1比例为1∶1，且满足正交性。
而且矩阵Hk的第一行和第一列均为1，要想将矩阵Hk打乱为纯随机的矩阵，所需的排列矩阵的规模为因此打乱操作具有极大的计算复杂度。
此外，若将矩阵Hk做平移、拉伸变换使之成为0-1矩阵不能满足压缩感知测量矩阵的受限等距性质(Restricted Isometry Property，RIP)，即不适合用于压缩感知(CS)重建。因此，变换后的0-1矩阵不能直接加载在空间光调制元件2-1上。
本发明一实施方式提出的将矩阵Hk变换为满足压缩感知受限等距性质(RIP)的测量矩阵Φ的方法如下：
在此之前，首先介绍待测对象的像素概念，待测对象的二维图像可以被分割成如下所示的2×3的矩阵，其总像素N＝p×q＝3×2＝6(p为横坐标像素数、q为纵坐标像素数)。
其次，所采用的Hadamard矩阵为2k阶，且2k-1≥N，因此k最小为3，即所采用的Hadamard矩阵最小为8×8的矩阵。
为了更为清晰的表述，以图9所示的8×8的Hadamard矩阵为例。
①删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ(暂时不考虑系数)。
而(T表示转置)为对角线元素为2k-1，其余元素为-1的(2k-1)×(2k-1)的方阵，则为对角阵。因此，HΩ是近似正交的。即：
②将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，即可得到0-1分布的矩阵对于2k-1阶的0-1矩阵其每一行中有2k-1个0和2k-1-1个1，即0的数量比1的数量多一个。
③以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵得到矩阵H+(仍不满足RIP条件)，其中扩充排列矩阵P对矩阵H+的列进行打乱操作，扩充排列矩阵Q对矩阵H+的行进行打乱操作。
例如：排列矩阵P和Q分别为：
④提取上述矩阵H+中M行，并保留M行中的N列作为随机测量矩阵Φ，其中，M即可代表等效采样数。随机测量矩阵Φ满足RIP条件，适用于压缩感知(CS)重建。
例如：令M＝3，提取第2-4行和1-6列作为随机测量矩阵Φ，即：
根据上法，实际上是设计出一种基于Hadamard矩阵与排列矩阵P和Q的随机测量矩阵：Φ＝QTHΩP，其采样率为M/N×100％，且M＜＜N/2，本实施方式中采样率为5％～30％。
所采用的扩充排列矩阵P和Q均为随机生成的，而非固定不变，则每次这样产生的测量矩阵均不相同，随机性更好。
⑤空间光调制器2-1交替加载Φ矩阵的第i行拉伸变换而得的掩膜和该行的补集拉伸变换而得的掩膜。
即空间光调制器2-1交替加载互补掩膜直到第M次结束。
⑥温度确定装置5根据Y1 i＝E1(T)2i-1-E1(T)2i、结合数学模型Y1 i＝ΦS1计算出待测对象在波长λ1下的二维图像S1中的部分信号值Xi(λ1)’，根据压缩感知原理，恢复出二维图像S1中的全部信号值Xi(λ1)，构建测对象在波长λ1下的二维图像S1。
例如，根据上述加载的掩膜，获得：
即X4(λ1)+X6(λ1)＝E1(T)1-E1(T)2 (1)；
即X1(λ1)+X5(λ1)+X6(λ1)＝E1(T)3-E1(T)4 (2)；
X1(λ1)+X3(λ1)+X4(λ1)＝E1(T)5-E1(T)6 (3)；
联立方程(1)～(3)可以解得得待测对象在波长λ1下的二维图像S1：
进一步地，温度确定装置5根据压缩感知原理在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的条件下，建立数学模型Y＝ΦΨx’+e，通过压缩感知算法进行凸优化，获得部分信号x’后，再由反演出全部信号值Xi(λ1)，即获得待测对象在波长λ1下的二维图像S1：
同样地，温度确定装置5根据结合数学模型计算出待测对象在波长λ2下的二维图像S2中的部分信号值Xi(λ2)’。
进而，根据压缩感知原理在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的条件下，建立数学模型Y＝ΦΨx’+e，通过压缩感知算法进行凸优化，获得部分信号x’后，再由反演出全部信号值Xi(λ2)，即获得待测对象在波长λ2下的二维图像S2：
其中，稀疏基Ψ采用梯度基、离散余弦变换基、小波基、傅里叶变换基、Gabor变换基中的任意一种实现。
二维图像S1和S2的重构可以采用MP、OMP、StOMP、ReOMP、SPARCO、LARS、BCS、FBMP、SL0、PPPA、KF-CS、BPDQ、k-t FOCUSS、l1-Magic、IST、TwIST、FISTA、SOCP、RecPF、FPC、GPSR、SpaRSA、ALM、FALM、TVAL3、CoSaMP、LBI、SP、l1_ls、SMP、TFOCS、NESTA、SALSA、SPGL1、YALL1、l0重建算法、l1重建算法、l2重建算法、BPDN算法、BPDN的LAG松弛算法、LASSO算法、TV算法中的任意一种实现。
由于压缩感知理论为现有技术，因此在此不再赘述。相关理论可以参见：
[1]E.J.Candès，J.Romberg，and T.Tao，“Robust uncertainty principles：exact signal reconstruction from highly incomplete frequency information，”IEEE Trans.Inform.Theory52(2)，489-509(2006).
[2]E.J.Candès，J.Romberg，and T.Tao，“Stable signal recovery fromincomplete and inaccurate measurements，”Commun.Pure Appl.Math.59 1207-1223(2006).
[3]D.Donoho，“Compressed sensing，”IEEE Trans.Inform.Theory52，1289-1306(2006).
[4]E.J.Candès，“Compressive sampling，”in Proc.Int.Cong.Math.，EuropeanMathematical Society，Madrid，Spain，3，1433-1452(2006).
[5]E.J.Candès，and M.B.Wakin，“An introduction to compressivesampling，”IEEE Signal Process.Mag.25(2)，21-30(2008).
[6]D.Takhar，J.N.Laska，M.B.Wakin，M.F.Duarte，and D.Baron，“A newcompressive imaging camera architecture using optical-domain compression，”inProc.SPIE：Computational Tmaging IV，San Jose，CA，USA.
[7]J.Romberg，“Imaging via compressive sampling，”IEEE SignalProcess.Mag.25，14-20(2008).
[8]M.F.Duarte，M.A.Davenport，D.Takhar，J.N.Laska，T.Sun，K.F.Kelly，andR.G.Baraniuk，“Single-pixel imaging via compressive sampling，”IEEE SignalProcess.Mag.25(2)，83-91(2008).
⑦温度确定装置5对待测对象在波长λ1和波长λ2下的二维图像S1和S2进行定点除法运算(S1./S2)，获得：
根据Ti＝B’/ln(Xi/A’)，计算该像素坐标下的温度T1；
根据Ti＝B’/ln(Xi/A’)，计算该像素坐标下的温度T2；
……
根据Ti＝B’/ln(Xi/A’)，计算该像素坐标下的温度T6；
本发明利用两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵并由打乱后的矩阵中提取M行(保留M行中的N列)重新组成随机测量矩阵Φ。之后，在空间光调制元件上采用互补调制的策略，即交替加载随机测量矩阵Φ的第i行和该行的补集拉伸成的像素为P×Q的掩膜，其实现了正负光强调制，等价于在空间光调制元件上直接加载了±1分布的Hadamard矩阵，非传统方法中简单地将Hadamard矩阵经平移拉伸变换所得的0-1矩阵直接加载在空间光调制元件上并直接参与计算。
经过此法，使得测量矩阵满足RIP条件，并适合于CS重建计算，使得温度场图像的重建质量提高1到2个数量级。
图10示出了本发明另一些实施方式的基于压缩感知的双波长温度场成像设备的结构示意图。该图示出了几种变形的实施方式。图10实施例是在上述图8实施例的基础之进行变形得出的。在此着重描述二者不同之处，二者相同或者相似之处不再赘述。
图10所示的基于压缩感知的双波长温度场成像设备的第二种实施例为：
该实施例是在图8所示的第一种实施例的基础上增加第一会聚元件3-3和第二会聚元件4-3。该实施方式的基于压缩感知的双波长温度场成像设备的各部件之间的连接关系可以为：
第一探测装置3-4为第一点探测器，第二探测装置4-4为第二点探测器。第一会聚元件3-3布置在所述第一路径上、位于所述第一点探测器与所述第一滤光元件之间。第二会聚元件4-3布置在所述第二路径上、位于所述第二点探测器与所述第二滤光元件之间。第一点探测器3-4位于所述第一会聚元件的光焦点处。第二点探测器4-4位于第二会聚元件的光焦点处。其产生的技术效果与基于光辐射的测温设备的第二种实施例的技术效果相同，在此不再赘述。
图10所示的基于压缩感知的双波长温度场成像设备的第三种实施例为：
该实施例是在上述第二种实施例的基础上增加了布置在第一路径上、第一滤光元件3-1和第一会聚元件3-3之间的第一光强衰减元件3-2，和布置在第二路径上、第二滤光元件4-1和第二会聚元件4-3之间的第二光强衰减元件4-2。
由此，本装置通过设置光强衰减元件(例如中性密度滤光片)可以衰减光线的光电信号参量，以防止强光对后续点探测器的损伤，提高了测温精度，延长了产品的使用寿命。
图10所示的基于压缩感知的双波长温度场成像设备的第四种实施例为：
该实施例是在上述第三实施方式的基础上增加了与第一点探测器3-4和温度确定装置5连接的第一存储装置3-5，以及与第二点探测器4-4和温度确定装置5连接的第二存储装置4-5，它们分别用于存储每次测量到的第一光电信号参量和第二光电信号参量，从而减轻将温度确定装置5作为存储装置时对其造成的存储负担。
图10所示的基于压缩感知的双波长温度场成像设备的第五种实施例为：
在上述各实施例的基础上增加透镜1，使得待测对象的光辐射进入主光路，可以射向光辐射调制装置2’。
图10所示的基于压缩感知的双波长温度场成像设备的第六种实施例为：
当待测对象距离本系统的距离遥远时，在上述各实施例的基础上增加望远镜单元9，以使得待测对象的光辐射能够通过望远镜单元射向光辐射调制装置2’。
在本发明的一些实施方式中，所述望远镜单元可以是反射式望远镜、折射式望远镜、折反式望远镜中的任意一种。
在本发明的一些实施方式中，所述望远镜单元可以是伽利略望远镜、开普勒望远镜、牛顿望远镜、格里望远镜、卡塞格林望远镜、施密特-卡塞格林望远镜、马克苏托夫-卡塞格林望远镜、多镜面望远镜、双筒望远镜、空间天文望远镜中的任意一种。
本实施例中望远镜单元9采用反射式望远镜，其包括布置于光学成像透镜1和被测物体之间的凹面反射镜9-2和凸面反射镜9-1。凹面反射镜9-2位于凸面反射镜9-1的后方，且凸面反射镜9-1的凸面背向被测物体，凹面反射镜9-2的凹面朝向被测物体。
参照图10，下面介绍一优选实施方式的基于压缩感知的双波长温度场成像设备的工作方式，其具体可以是：
待测对象(例如人体或者电灯，图中未进行标示)的光辐射可以经过望远镜单元9到达透镜1，再射向光辐射调制装置2’。光辐射调制装置2’加载预设的多个掩膜(所述掩根据符合压缩感知理论的RIP条件的测量矩阵变换生成)，并且将接收的光辐射调制反射为双臂方向上的两束光辐射。在双臂中的左臂方向上分别布置有第一滤光元件3-1、第一光强衰减元件3-2、第一会聚元件3-3和第一探测装置3-4。第一束光辐射经过第一滤光元件3-1过滤为波长为第一波长的光，波长为第一波长的光经过第一光强衰减元件3-2进行光强衰减，再经过第一会聚元件3-3进行会聚，第一探测装置3-4在会聚的波长为第一波长的光焦点所在的焦平面处接收波长为第一波长的光来探测波长为第一波长的光的光电信号参量。第一存储装置3-5存储第一光电信号参量。
本设备左右对称设置。同样，在双臂中的右臂方向上分别布置有第二滤光元件4-1、第二光强衰减元件4-2、第二会聚元件4-3和第二探测装置4-4。第二束光辐射经过第二滤光元件4-1过滤为波长为第二波长的光，波长为第一波长的光经过第二光强衰减元件4-2进行光强衰减，再经过第二会聚元件4-3进行会聚，第二探测装置4-4在会聚的波长为第二波长的光焦点所在的焦平面处接收波长为第二波长的光来探测波长为第二波长的光的光电信号参量。第二存储装置4-5存储第二光电信号参量。
温度确定装置5分别接收第一存储装置3-5存储的和第二存储装置4-5存储的波长为第一波长λ1的光的光电信号参量和波长为第二波长λ2的光的光电信号参量，并利用压缩感知算法重构出待测对象在两个波长λ1和λ2下的的二维图像S1和S2。
图像生成装置8反演出待测对象二维红外热图像。
在一些实施例中，为了使过滤后的单波的波长不同，第一滤光元件和所述第二滤光元件为中心波长相差10nm以上的第一窄带滤光片和第二窄带滤光片，所述第一窄带滤光片和第二窄带滤光片的半高宽参数为10nm以上。
因为根据窄带带宽理论，窄带滤光片越窄效果越好；中心波长越接近越好。但是，带宽越窄，与窄带滤光片配合使用的光电信号参量探测元件能探测到的热辐射越少，外加之中心波长接近，对探测器的灵敏度和准确性提出过高的要求，而且此时环境噪声和探测器本征噪声(如暗计数等)也将对系统的测温准确性造成较大影响。所以，在衡量上述优缺点的基础上，经过大量的试验，选用的窄带滤光片的半高宽FWHM一般应在10nm以上，中心波长CWL一般相差10nm以上效果最佳。另外，本设备可以利用不同滤光元件(例如窄带滤光片)进行滤波得到单波光，再将该单波光会聚收集成以供探测器检测，降低了无关光的干扰，提高了采集温度的精度。
在一些实施例中，所述光辐射为红外线光辐射。
在一些实施例中，所述第一和第二点探测器选自近红外、中远红外、远红外波段的外光电效应探测器组、内光电效应探测器组、强光探测器组和弱光探测器组中的任意一种，其中，
所述外光电效应探测器组包括：雪崩二极管、真空光电管、充气光电管、光电倍增管、变像管、像增强器、摄像管；
所述内光电效应探测器组包括：本征型光电导探测器、掺杂型光电导探测器、光磁电效应探测器、光生伏特探测器；
所述强光探测器组包括：内置或者外置有模数转换器的强光探测器；
所述弱光探测器组包括：内置或者外置有计数器的弱光探测器。
在一些实施例中，所述光电信号参量包括光子数、电流值、电压值、电阻值中任意一种。
由此，点探测器可以自由选择各种类型以满足各种需求，增强了产品的通用性能，也便于后期的维护保养。
图11示出了本发明一种实施方式的基于压缩感知的双波长温度场成像系统的结构示意图。该图可以有几种变形的实施方式。在此着重描述这些变形方式的不同之处，它们相同或者相似之处不再赘述。
图11所示的基于压缩感知的双波长温度场成像系统的第一种实施例为：
参照图11，该成像系统包括：上述的基于压缩感知的双波长温度场成像设备和定标装置6。
定标装置6可以包括温度可调节的基准光源6-1(例如，可以提供不同功率的灯泡)、扩束准直透镜6-2和分束器6-3。分束器6-3设置在待测对象7和透镜1之间。基准光源6-1、扩束准直透镜6-2和分束器6-3水平共线设置，扩束准直透镜6-2设置在基准光源6-1和分束器6-3之间。其中，基准光源6-1用于在定标阶段将调节的不同温度的光辐射射向光辐射调制装置2’，并由光辐射调制装置2’将接收到的光辐射等分为第一光辐射和第二光辐射，使第一光辐射沿第一路径射出(如左臂方向路径)、第二光辐射沿第二路径(如右臂方向路径)射出，以获取多个第一和第二光电信号参量，并确定调节的不同温度和获取的多个第一和第二光电信号参量的预定关系。扩束准直透镜6-2用于将基准光源的光辐射转化为平行光辐射。分束器6-3用于将扩束准直透镜转化的平行光辐射射向所述光辐射调制装置2’。
在上述定标过程中，标准光源以第一电流、电压、电阻的光线照射扩束准直透镜，通过分束器、透镜到达DMD。保持DMD同一帧光辐射不变，其加载的掩膜为1和0数量相等的0-1矩阵。加载的掩膜中1和0的数量相等。
图11所示的基于压缩感知的双波长温度场成像系统的第二种实施例为：
在图11所示的第一实施例的基础上减少了分束器6-3，该分束器的功能随之减少。但此实施方式仍能解决技术问题，实现相应的技术效果。
图11所示的基于压缩感知的双波长温度场成像系统的第三种实施例为：
在图11所示的第二实施例的基础上减少了扩束准直透镜6-2，该扩束准直透镜的功能随之减少。但此实施方式仍能解决技术问题，实现相应的技术效果。
在上述实施方式中，所述确定调节的不同温度和获取的多个第一和第二光电信号参量的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量。
由此，本系统可以采用不同功率的标准光源进行试验测温，通过测得多组数据来确定试验中的拟合系数。其中，拟合算法利用已知试验或者真实数据，然后寻找一个模型对其规律进行模拟的过程中，求取模型中未知参数的一个过程。以确保在后期的实际测温过程中根据该拟合系数获取测得的温度。本实施方式通过多次试验测进行试验测温，提高了实际测温的精度。此外，本实施方式利用扩束准直透镜6-2将标准光源的光线变换成平行光，减少了光线因为会聚或者散射引起的误差。利用分束器6-3可以最大程度的将标准光源的光线全部传递至透镜，提高了光线的传递效率。
本领域的技术人员可以理解，可以根据实际测量精度或者特殊需要对上述各元件进行选择配置，或者进行多种组合配置。
由此，本系统通过提供定标装置，可以在正式测温之前，进行试验测温，从而制定标准，使得后面的正式测温可以参考试验测温的数据进行运算，调整测得的温度数据，进一度调高了温度的精度。
本实施方式的测温系统的测温的工作方式可以参照上述的测温设备的描述。需要说明的是，在定标的试验之前，将分束器6-3、扩束准直透镜6-2和标准光源2-1移入系统中，保证分束器6-3的反射方向在透镜1和DMD2-1主轴光路上，保证标准光源的光辐射进入系统的主轴光路，且标准光源与待测对象7所在位置共轭，待定标结束后再将分束器6-3、扩束准直透镜6-2和标准光源6-1移出系统。
图12为本发明一种实施方式的基于压缩感知的双波长温度场成像方法的流程示意图。如图12所示，该方法包括以下步骤：
S1201：利用光辐射调制装置(可以利用图8和图10实施方式中的光辐射调制装置))接收待测对象(例如人体或者电灯)的光辐射(例如红外辐射、紫外辐射或者可见光)，加载预设的多个掩膜，而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、多束第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出，所述多个掩膜由符合压缩感知RIP条件的随机测量矩阵Φ变换生成；
S1202：在上述的左臂方向的路径上接收多束第一光辐射，并将其过滤为波长为第一波长的多束光(例如波长为λ1的单波光)；
S1203：在上述的右臂方向的路径上接收多束第二光辐射，并将其过滤为波长为第二波长的多束光(例如波长为λ2的单波光)；
S1204：在上述的左臂方向的路径上接收多束波长为λ1的单波光并将其转换为多个相应的第一光电信号参量E1(T)i；
S1205：在上述的右臂方向的路径上接收多束波长为λ2的单波光并将其转换为多个相应的第二光电信号参量E2(T)i；
S1206：根据多个所述第一和第二光电信号参量与温度的预定关系确定出所述待测对象的每一个像素点的温度Ti；以及
S1207：根据所述待测对象的二维图像以及每一个像素点的温度值Ti反演出所述待测对象的二维红外热图像。
在本实施方式中，所述光辐射调制装置根据预设矩阵Φ生成和加载掩膜的方式如下：
a.提供2k阶的Hadamard矩阵(Hk)，其中且2k-1≥N＝p×q(p为待测对象横坐标像素数、q为待测对象纵坐标像素数)；
b.删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ；
c.将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，获得0-1分布的矩阵
d.以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵获得打乱后的矩阵H+；
e.提取矩阵H+中的M行，并保留该M行中的N列作为随机测量矩阵Φ，其中，M＜＜N/2。
N为待测对象的总像素数，M表示采样数(随机测量矩阵的行)。随机测量矩阵Φ满足RIP条件，适用于压缩感知(CS)重建。
f.加载由随机测量矩阵Φ的第i行拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射，将第一光辐射转换为相应的第一光电信号参量E1(T)2i-1，第将第二光辐射转换为相应的第二光电信号参量E2(T)2i-1，i＝1，2，3，4…M；
g.加载由随机测量矩阵Φ第i行的补集拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射，将第一光辐射转换为相应的第一光电信号参量E1(T)2i，将第二光辐射转换为相应的第二光电信号参量E2(T)2i，i＝1，2，3，4…M。
在本实施方式中，所述根据如下的预定关系确定待测对象的每个像素点的温度值：
如上所示，待测对象的二维图像以矩阵的形式可以表示为：
根据结合数学模型计算出待测对象在波长λ1下的二维图像S1中的部分信号值Xi(λ1)’；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出全部待测对象在波长λ1下的二维图像S1中的全部信号值Xi(λ1)，重构待测对象在波长λ1下的二维图像S1；
根据结合数学模型计算出待测对象在波长λ2下的二维图像S2中的部分信号值Xi(λ2)’；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出全部待测对象在波长λ2下的二维图像S2中的全部信号值Xi(λ2)，重构待测对象在波长λ2下的二维图像S2；
根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi＝S1./S2以及所述Ti＝B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti；
其中，A’、B’为预设系数。
所述稀疏基Ψ采用梯度基、离散余弦变换基、小波基、傅里叶变换基、Gabor变换基中的任意一种实现。
所述待测物体在波长λ1和波长λ2下的二维图像S1和S2的重构可以采用MP、OMP、StOMP、ReOMP、SPARCO、LARS、BCS、FBMP、SL0、PPPA、KF-CS、BPDQ、k-t FOCUSS、l1-Magic、IST、TwIST、FISTA、SOCP、RecPF、FPC、GPSR、SpaRSA、ALM、FALM、TVAL3、CoSaMP、LBI、SP、l1_ls、SMP、TFOCS、NESTA、SALSA、SPGL1、YALL1、l0重建算法、l1重建算法、l2重建算法、BPDN算法、BPDN的LAG松弛算法、LASSO算法、TV算法中的任意一种实现。
二维红外热像图与物体表面的热分布相对应，也即反映被测对象的光辐射(红外热辐射)能量分布，二维红外热像图的上面的不同颜色代表被测物体不同区域的不同温度，通过二维红外热像图可以观察到被测物体的整体温度分布状况，研究被测物体的发热情况。
在一些实施方式中，该方法还包括：
将所述波长为第一波长的光会聚于第一焦点，并在所述第一焦点处设置第一点探测装置用以接收多束所述波长为第一波长的光并将其转换为多个相应的第一光电信号参量；
将所述波长为第二波长的光会聚于第二焦点，并在所述第二焦点处设置第二点探测装置用以接收多束所述波长为第二波长的光并将其转换为多个相应的第二光电信号参量。
在一些实施方式中，该方法还包括：
对所述波长为第一波长的多束光的强度进行衰减；及对所述波长为第二波长的多束光的强度进行衰减。
在一些实施方式中，该方法还包括：
利用望远镜单元接收待测对象的光辐射。
在本实施方式中，所述光辐射为红外线光辐射。
在本实施方式中，所述第一波长与所述第二波长相差至少10nm。
在一些实施方式中，所述光电信号参量包括光子数、电流值、电压值、电阻值中任意一种。
在一些实施方式中，所述第一和第二点探测器选自近红外、中远红外、远红外波段的外光电效应探测器组、内光电效应探测器组、强光探测器组和弱光探测器组中的任意一种，其中，
所述外光电效应探测器组包括：雪崩二极管、真空光电管、充气光电管、光电倍增管、变像管、像增强器、摄像管；
所述内光电效应探测器组包括：本征型光电导探测器、掺杂型光电导探测器、光磁电效应探测器、光生伏特探测器；
所述强光探测器组包括：内置或者外置有模数转换器的强光探测器；
所述弱光探测器组包括：内置或者外置有计数器的弱光探测器。
在测温方法上述实施例中的技术效果与测温装置实施例中的技术效果相对应，在此不再赘述。
本发明一种实施方式的基于压缩感知的双波长温度场成像方法同样可以包括如图7所示的定标的步骤，即在步骤S1201之前还包括：
S1301：利用光辐射调制装置接收基准光源(例如温度可调节的灯泡，具体可以调节电流、电压等方式来调节)的光辐射，将接收到的光辐射等分(此处等分的概念可以参考上述测温设备中等分的概念)为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、所述第二光辐射沿不同于第一路径的第二路径(如测温设备的右臂方向的路径)射出。
S1302：在所述第一路径上接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光(例如波长为λ1的单波光)。
S1303：在所述第二路径上接收所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光(例如波长为λ2的单波光)。
S1304：在所述第一路径上接收所述波长为λ1的单波光并将其转换为相应的第一光电信号参量E1(T)。
S1305：在所述第二路径上接收所述波长为λ2的单波光并将其转换为相应的第二光电信号参量E2(T)。
S1306：调节灯泡的电流或电压，使得灯泡发出光辐射的温度变为Td2，并获取对应的光电信号参量Ed1(T)2和Ed2(T)2。按此种方法可以获取灯泡发出多个不同温度Tdi的光辐射，并获取对应的多个第一和第二光电信号参量Ed1(T)i和Ed2(T)i。
S1307：根据上述不同温度与第一和第二光电信号参量的计量值，确定出所述第一和第二光电信号参量与温度的预定关系。
在本实施方式中，所述确定多个第一和第二光电信号参量与调节的多个不同温度的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，所述多个为n个，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量。
具体实现方式可以是：
在定标时，标准光源以第一电流、电压、电阻的光线照射扩束准直透镜，通过分束器、透镜到达DMD。保持DMD同一帧光辐射不变，其加载的掩膜为1和0数量相等的0-1矩阵。控制元件控制DMD中各微镜的闭合、翻转，使得DMD将接收到的一帧光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出。其中：第一光辐射经过第一窄带滤光片、第一会聚元件到达第一点探测器，得到第一光电信号参量Ed1(T)1。第二光辐射经过第二窄带滤光片、第二会聚元件到达第二点探测器，得到第二光电信号参量Ed2(T)1。Ed1(T)1和Ed2(T)1经过除法器得到比值X1，X1＝Ed1(T)1/Ed2(T)1；
调节标准光源的不同等效电流或电压或电阻可计算出不同电流或电压或电阻下的一系列温度T1，T2，T3…Tdi，通过本系统测得对应的Ed1(T)i与Ed2(T)i比值X1，X2，X3…Xi，保持DMD固定一帧不变，其DMD上加载的掩膜为1和0数量相等的0-1矩阵。接着利用公式拟合得出相应系数，即得到此种环境下的辐射体(即基准光源)的光辐射温度曲线公式。即：确定拟合系数(即上述的预设系数)A’和B’，其中：
Ed1(T)1/Ed2(T)1＝A’exp(B’/Td1)
Ed1(T)2/Ed2(T)2＝A’exp(B’/Td2)
Ed1(T)3/Ed2(T)3＝A’exp(B’/Td3)
……
Ed1(T)＝A1(λ1)exp(B1(λ1)/T)
Ed2(T)＝A2(λ2)exp(B2(λ2)/T)
A’＝A1(λ1)/A2(λ2)，B’＝B1(λ1)-B2(λ2)。
具体可以是利用公式拟合得出相应系数A’和B’，得到此种环境下的该辐射体的光辐射温度曲线公式；然后通过计算机程序设定好此系数下的计算公式，即可测量或者监测任何情况下的辐射温度。在实际测温过程中，利用指数形式的多项式展开，获得多项式的系数，或者利用其它复杂的拟合方程获得拟合系数。
由此，本实施方式有效地克服了红外测温中各种物体的“辐射率修正”难题，克服了测量条件复杂、现场测量条件波动或者水蒸气等因素的环境吸收造成的测量误差，实现了高精度的温度测量。
需要指出的是，如图8～13所示的本发明所提供的基于压缩感知的双波长温度场成像设备、系统和方法只能够获得待测对象(例如人体或者电灯)的二维红外热图像，不能获得待测对象的三维结构，也不能生成待测对象的三维红外热图像。
鉴于此，本发明提供了一种结合时间分辨技术的基于压缩感知的双波长三维温度场成像设备。
时间分辨技术：通常利用飞行时间(Time of Flight，TOF)照相机向被测物体连续发送光脉冲，然后利用传感器接收从被测物体返回的光，通过探测光脉冲的飞行(往返)时间来得到目标物距离。
参照图14所示，本发明的基于压缩感知的双波长三维温度场成像设备的第一种实施例中，该设备与如图8所示的基于压缩感知的双波长温度场成像设备存在如下不同：
①第一探测装置3-4用于探测第一光辐射的光子数，第二探测装置4-4用于探测第二光辐射的光子数；
②增加了设置于第一探测装置3-4(或第二探测装置4-4)与温度确定装置5之间的时幅变换器11，时幅变换器11能够利用时间相关技术记录第一或者第二光辐射的光子团到达时间以确定时间戳信息；
③增加了用于照射待测对象的脉冲光源10-1；
④图像生成装置8能够将时间戳信息转化为待测对象的每一个横截层面至特定位置的距离信息，然后生成待测对象每一个横截层面上的二维红外热图像，进而结合多个横截层面的二维红外热图像生成待测对象的三维红外热图像。
具体而言，在第一个实施例中，基于压缩感知的双波长三维温度场成像设备包括：用于向待测物体发射脉冲光的脉冲光源10-1，光辐射调制装置2’、第一滤光元件3-1、第二滤光元件4-1、第一探测装置3-4、第二探测装置4-4，与第一探测装置3-4、第二探测装置4-4连接的温度确定装置5，与温度确定装置5连接的图像生成装置8，以及与第一探测装置3-4和温度确定装置5连接的时幅变换器11。
光辐射调制装置2’可以加载根据符合压缩感知RIP条件的矩阵变换生成的掩膜。
光辐射调制装置2’包括：空间光调制器1-1，配置为根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成的多个掩膜，以将接收到的待测对象的光辐射调制为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出；控制元件，配置为控制所述空间光调制器依次加载由所述测量矩阵Φ变换生成的多个掩膜。例如，根据Hadamard矩阵变换得到符合压缩感知受限等距性质(RIP)的测量矩阵，并根据测量矩阵变换生成掩膜，最后根据掩膜中0和1的数量控制空间光调制器2-1中微镜的翻转。
以脉冲光源10-1照射待测对象7，脉冲光经待测对象7反射后射向光辐射调制装置2’。光辐射调制装置2’将接收到的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出。在所述第一路径上接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光；在所述第二路径上接收所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光。时幅变换器11记录第一束聚集光子团最初到达第一探测装置3-4时刻作为初始时间节点t1，之后依次记录第二束聚集光子团最初到达第一探测装置3-4的时刻t2、第三束聚集光子团最初到达第一探测装置3-4的时刻t3、直至第r(r≥1)束聚集光子团最初到达第一探测装置3-4的时刻tr。其中，r为自然数。
第二束聚集光子团最初到达第一探测装置3-4时刻与第一束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔t2-t1；第三束聚集光子团最初到达第一探测元件3-4的时刻与第二束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔t3-t2……第r束聚集光子团最初到达第一探测装置3-4的时刻与第r-1束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔tr-tr-1，因此，在上述时间间隔内，每一束聚集光子团出现峰值的时刻为：{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}，即时间戳信息。
例如：第一束光子团在t1时刻最初达到并且持续到t2时刻，其中间时刻t1+(t2-t1)/2＝(t1+t2)/2出现第一束光子团的峰值，以此时刻(t1+t2)/2作为第一束光子团的时间戳信息。
根据上述时间戳信息生成光子到达时间的直方图。所述直方图中的峰值代表着待测对象在不同的纵向距离上的信息(即至光辐射调制装置2’具有相同距离的横截层面)，其各峰值附近形成聚集光子团。根据光子的到达时间以及光速，可以推算出待测对象7的不同横截层面至同一位置的距离(例如第一探测装置3-4或第二探测装置4-4)。若能够获得待测对象每一个横截层面的二维红外热像图，则结合时间戳信息就能够获得待测对象7完整的三维红外热像图。
以脉冲光源10-1继续照射待测对象7，脉冲光经待测对象7反射后射向光辐射调制装置2’。光辐射调制装置2’加载根据符合压缩感知RIP条件的矩阵变换生成的多个掩膜，并使第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出。布置在第一路径上的第一滤光元件3-1接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光(例如波长为λ1的单波光)。布置在第二路径上的第二滤光元件4-1接收第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光(例如波长为λ2的单波光)。
在加载每一个掩膜时，布置在第一路径上第一探测装置3-4在上述的多个时间间隔内(t2-t1、t3-t2…tr-tr-1)接收波长为第一波长的光并将其转换为相应的多个第一光电信号参量。布置在第二路径上的第二探测装置4-4在上述的多个时间间隔内(t2-t1、t3-t2、…、tr-tr-1)分别接收波长为第二波长的光并将其转换为相应的多个第二光电信号参量。
以此类推，光辐射调制装置2’每加载一个掩膜(翻转到下一帧)，第一探测装置3-4和第二探测装置4-4都在多个时间间隔内(t2-t1、t3-t2、…、tr-tr-1)测量获得与时间戳信息相关的多个第一光电信号参量和多个第二光电信号参量，直至所有掩膜加载结束。因此，在同一个时间戳信息下，都会探测获得与光辐射调制装置2’加载的所有掩膜对应的一组第一光电信号参量和一组第二光电信号参量。
温度确定装置5接收来自第一探测装置3-4和第二探测装置4-4的多组第一和多组第二光电信号参量，并根据多组所述第一光电信号参量以及时间戳信息确定待测对象在波长λ1和时间戳信息[tr-tr-1]下的二维图像S1[tr-tr-1]，根据多组所述第二光电信号参量以及时间戳信息确定待测对象在波长λ2和时间戳信息[tr-tr-1]下的二维图像S2[tr-tr-1]。
温度确定装置5根据待测对象在时间戳信息[tr-tr-1]以及波长λ1和波长λ2下的二维图像S1[tr-tr-1]和S2[tr-tr-1]以及所述待测对象的温度的预定关系确定出待测对象每每一个时间戳信息时的每一个像素点的温度Ti[tr-tr-1]。
图像生成装置8根据所述待测对象在每一个所述时间戳信息[tr-tr-1]时的每一个像素点的温度值Ti[tr-tr-1]以及待测对象在每一个所述时间戳信息时的二维图像S[tr-tr-1]反演出待测对象的每一个横截层面(将时间戳信息转化为距离信息，以获得到达特定位置不同距离的层面)的二维红外热图像，且根据待测对象的每一个横截层面的二维红外热图像生成三维红外热图像。
在本发明的实施方式中，测量矩阵Φ由Hadamard矩阵变换得到，其变换的理论模型来源于法国数学家Hadamard提出的一种2k阶矩阵方程。
Hadamard矩阵满足以下特点：
Hadamard矩阵Hk中的元素只有1和-1(不考虑系数)，其中1和-1比例为1∶1，且满足正交性。
若将其做平移和拉伸变换使之成为0-1测量矩阵，然后直接加载在空间光调制元件2-1上，其被证实不满足压缩感知测量矩阵的受限等距性质(Restricted IsometryProperty，RIP)，即不适合用于压缩感知(CS)重建。
而且矩阵Hk的第一行和第一列均为1，要想将矩阵Hk打乱为纯随机的矩阵，所需的排列矩阵的规模为因此打乱操作具有极大的计算复杂度。
本发明一实施方式提出的将Hadamard矩阵变换为满足压缩感知受限等距性质(RIP)的测量矩阵方法如下：
在此之前，首先介绍待测对象的像素概念，待测对象的二维图像可以被分割成如下所示的2×3的矩阵，其总像素N＝p×q＝3×2＝6(p为横坐标像素数、q为纵坐标像素数)。
其次，所采用的Hadamard矩阵为2k阶，且2k-1≥N，因此k最小为3，即所采用的Hadamard矩阵最小为8×8的矩阵。
为了更为清晰的表述，以图9所示的8×8的Hadamard矩阵为例。
一、时间戳信息的确定：
空间光调制元件2-1(DMD)(采用掩膜为1和0数量相等的0-1矩阵固定一帧不变)；
以脉冲光源10-1照射待测对象7，脉冲光经待测对象7反射后射向光辐射调制装置2’。
光辐射调制装置2’将接收到的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出。在所述第一路径上接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光；在所述第二路径上接收所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光。
时幅变换器11记录第一束聚集光子团最初到达第一探测装置3-4(或者第二探测装置4-4)时刻作为初始时间节点t1，之后依次记录第二束聚集光子团最初到达第一探测装置3-4的时刻t2、第三束聚集光子团最初到达第一探测装置3-4的时刻t3、直至第r(r≥1)束聚集光子团最初到达第一探测装置3-4的时刻tr。
第二束聚集光子团最初到达第一探测装置3-4时刻与第一束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔t2-t1；第三束聚集光子团最初到达第一探测装置3-4的时刻与第二束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔t3-t2……第r束聚集光子团最初到达第一探测装置3-4的时刻与第r-1束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔tr-tr-1，因此，每一束聚集光子团出现峰值的时刻，即时间戳信息为：
{(t1+t2)/2，(t2+t3)/2，……，(tr-1+tr)/2}，其中，r为自然数。
例如：第一束聚集光子团在t1和t2的中间时刻(t1+t2)/2出现峰值，第二束聚集光子团在t2和t3的中间时刻(t2+t3)/2出现峰值……第r束聚集光子团在tr和tr-1的中间时刻(tr-1+tr)/2出现峰值。
二、矩阵变换和加载
①删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ(暂时不考虑系数)。
而(T表示转置)为对角线元素为2k-1，其余元素为-1的(2k-1)×(2k-1)的方阵，则为对角阵。因此，HΩ是近似正交的。
即：
②将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，即可得到0-1分布的矩阵对于2k-1阶的0-1矩阵其每一行中有2k-1个0和2k-1-1个1，即0的数量比1的数量多一个。
③以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵得到矩阵H+(仍不满足RIP条件)，其中扩充排列矩阵P对矩阵H+的列进行打乱操作，扩充排列矩阵Q对矩阵H+的行进行打乱操作。
例如：排列矩阵P和Q分别为：
④提取上述矩阵H+中M行，并保留M行中的N列作为随机测量矩阵Φ，其中，M即可代表等效采样数。随机测量矩阵Φ满足RIP条件，适用于压缩感知(CS)重建。
例如：令M＝3，提取第2-4行和1-6列作为随机测量矩阵Φ，即：
根据上法，实际上是设计出一种基于Hadamard矩阵与排列矩阵P和Q的随机测量矩阵：Φ＝QTHΩP，其采样率为M/N×100％，且M＜＜N/2，本实施方式中采样率为5％～30％。
所采用的扩充排列矩阵P和Q均为随机生成的，而非固定不变，则每次这样产生的测量矩阵均不相同，随机性更好。
⑤空间光调制器2-1交替加载Φ矩阵的第i行拉伸变换而得的掩膜和该行的补集拉伸变换而得的掩膜。
时幅变换器11开始计时，第一探测装置3-4在时间间隔t2-t1、t3-t2、…、tr-tr-1内记录第一光辐射的多个第一光子数E1(T)2i-1(tr-tr-1)，例如在t2-t1时间间隔内的光子总数E1(T)1(t2-t1)、在t3-t2时间间隔内的光子总数E1(T)1(t3-t2)、…在tr-tr-1时间间隔内的光子总数E1(T)1(tr-tr-1)，并且根据每一段时间间隔内光子团峰值出现的时间，赋予它们时间戳信息。例如，在t2-t1时间间隔内，光子团出现峰值的时间为(t1+t2)/2，则E1(T)1(t2-t1)记为E1(T)1[(t1+t2)/2]。以此类推，记为E1(T)1[(t2+t3)/2]、…E1(T)1[(tr-1+tr)/2]。
第二探测装置4-4在时间间隔t2-t1、t3-t2、…、tr-tr-1内记录第二光辐射的多个第二光子数E2(T)2i-1，例如，E2(T)1(t2-t1)、E2(T)1(t3-t2)、…E2(T)1(tr-tr-1)，并且根据每一段时间间隔光子团峰值出现的时间，赋予它们时间戳信息，在t2-t1时间间隔内，光子团出现峰值的时间为(t1+t2)/2，则E2(T)1(t2-t1)记为E2(T)1[(t1+t2)/2]。以此类推，获取E2(T)1[(t2+t3)/2]、…E2(T)1[(tr-1+tr)/2]，其中，i＝1，2，3，4……M；
第二探测装置4-4在时间间隔t2-t1、t3-t2、…、tr-tr-1内记录第二光辐射的多个第二光子数E2(T)2i，例如，E2(T)2[(t1+t2)/2]、E2(T)2[(t2+t3)/2]、…E2(T)2[(tr-1+tr)/2]，其中i＝1，2，3，4……M；
空间光调制器2-1交替加载互补掩膜直到第M次结束。
三、图像恢复
①温度确定装置5根据Y1 i＝E1(T)2i-1[(tr-1+tr)/2]-E1(T)2i[(tr-1+tr)/2]、结合数学模型Y1 i＝ΦS1[(tr-1+tr)/2]计算出待测对象在时间戳为[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]中的部分Xi(λ1)’[(tr-1+tr)/2]值。
例如，根据上述加载的掩膜，以t2-t1时间间隔内探测获得的一系列第一光子数为例(对应的时间戳信息为[(t1+t2)/2])：
即X4(λ1)[(t1+t2)/2]+X6(λ1)[(t1+t2)/2]＝E1(T)1[(t1+t2)/2]-E1(T)2[(t1+t2)/2] (1)；
即
X1(λ1)[(t1+t2)/2]+X5(λ1)[(t1+t2)/2]+X6(λ1)[(t1+t2)/2]＝E1(T)3[(t1+t2)/2]-E1(T)4[(t1+t2)/2] (2)；
即X1(λ1)[(t1+t2)/2]+X3(λ1)[(t1+t2)/2]+X4(λ1)[(t1+t2)/2]＝E1(T)5[(t1+t2)/2]-E1(T)6[(t1+t2)/2] (3)；
进一步地，温度确定装置5根据压缩感知原理在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的条件下，建立数学模型Y＝ΦΨx’+e，通过压缩感知算法进行凸优化，获得部分信号x’后，再由反演出全部信号值Xi(λ1)[(t1+t2)/2]，即获得待测对象在时间戳信息为[(t1+t2)/2]、波长λ1下的二维图像S1：
进而，根据压缩感知原理在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的条件下，建立数学模型Y＝ΦΨx’+e，通过压缩感知算法进行凸优化，获得部分信号x’后，再由反演出全部信号值Xi(λ2)[(t1+t2)/2]，即获得待测对象在是时间戳为[(t1+t2)/2]、波长λ2下的二维图像S2：
其中，稀疏基Ψ采用梯度基、离散余弦变换基、小波基、傅里叶变换基、Gabor变换基中的任意一种实现。
二维图像S1和S2的重构可以采用MP、OMP、StOMP、ReOMP、SPARCO、LARS、BCS、FBMP、SL0、PPPA、KF-CS、BPDQ、k-t FOCUSS、l1-Magic、IST、TwIST、FISTA、SOCP、RecPF、FPC、GPSR、SpaRSA、ALM、FALM、TVAL3、CoSaMP、LBI、SP、l1_ls、SMP、TFOCS、NESTA、SALSA、SPGL1、YALL1、l0重建算法、l1重建算法、l2重建算法、BPDN算法、BPDN的LAG松弛算法、LASSO算法、TV算法中的任意一种实现。
压缩感知理论已经在上文阐述，在此不再重复说明。
②温度确定装置5对待测对象在时间戳信息为[(tr-1+tr)/2]、波长λ1和波长λ2下的二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]进行定点除法运算，获得：
例如：根据Ti＝B’/ln(Xi/A’)，计算该像素坐标下的温度T1[(t1+t2)/2]，即温度T1[(t1+t2)/2]对应的像素坐标为X1[(t1+t2)/2]、对应的时间戳信息是[(t1+t2)/2]；
根据Ti＝B’/ln(Xi/A’)，计算该像素坐标下的温度T2[(t1+t2)/2]，即温度T2[(t1+t2)/2]对应的像素坐标为X2[(t1+t2)/2]、对应的时间戳信息是[(t1+t2)/2]；
……
根据Ti＝B’/ln(Xi/A’)，计算该像素坐标下的温度T6[(t1+t2)/2]，即温度T6[(t1+t2)/2]对应的像素坐标为X6[(t1+t2)/2]、对应的时间戳信息是[(t1+t2)/2]；
以此类推，能够依次地反演出待测对象在距离光调制装置2’为D2～3、D3～4、Dr-1～r的横截层面的二维红外热图像S[(t2+t3)/2]、S[(t3+t4)/2]、S[(t4+t5)/2]……。
④图像生成装置8根据待测对象在每一个横截层面的二维红外热图像，反演出完整的三维红外热图像。
二维红外热图像S[(t1+t2)/2]对应由时间戳信息[(t1+t2)/2]推算出的距离信息D1～2、二维红外热图像S[(t2+t3)/2]对应由时间戳信息[(t2+t3)/2]推算出的距离信息D2～3、二维红外热图像S[(t3+t4)/2]对应由时间戳信息[(t3+t4)/2]推算出的距离信息D3～4……二维红外热图像S[(tr-1+tr)/2]对应由时间戳信息[(tr-1+tr)/2]推算出的距离信息Dr-1～r，因此，每一个横截层面的二维红外热图像的位置关系是确定的，根据多个二维红外热图像及其位置关系，就能够反演出待测对象完整的三维红外热图像。
本发明利用两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵并由打乱后的矩阵中提取N×M的矩阵作为随机测量矩阵Φ。之后，在空间光调制元件上采用互补调制的策略，即交替加载随机测量矩阵Φ的第i行和该行的补集拉伸成的p×q像素的掩膜，其实现了正负光强调制，等价于在空间光调制元件上直接加载了±1分布的Hadamard矩阵，非传统方法中简单地将Hadamard矩阵经平移拉伸变换所得的0-1矩阵直接加载在空间光调制元件上并直接参与计算。
经过此法，使得测量矩阵满足RIP条件，并适合于CS重建计算，使得温度场图像的重建质量提高1到2个数量级。
图15示出了本发明另一些实施方式的基于压缩感知的双波长三维温度场成像设备的结构示意图。该图示出了几种变形的实施方式。图15实施例是在上述图14实施例的基础之进行变形得出的。在此着重描述二者不同之处，二者相同或者相似之处不再赘述。
图15所示的基于压缩感知的双波长三维温度场成像设备的第二种实施例为：
可以包括扩束准直透镜10-2和反射镜10-3。脉冲光源10-1、扩束准直透镜10-2和反射镜10-3水平共线设置。扩束准直透镜10-2设置在脉冲光源10-1和反射镜10-3之间。扩束准直透镜10-2用于将脉冲光源10-1的光辐射转化为平行光辐射。反射镜6-3用于将扩束准直透镜10-2转化的平行脉冲光射向所述待测对象7。
图15所示的基于压缩感知的双波长三维温度场成像设备的第三种实施例为：
该实施例是在图13所示的第二种实施例的基础上增加第一会聚元件3-3和第二会聚元件4-3。该实施方式的基于压缩感知的双波长温度场成像设备的各部件之间的连接关系可以为：
第一探测装置3-4为第一点探测器，第二探测装置4-4为第二点探测器。第一会聚元件3-3布置在所述第一路径上、位于所述第一点探测器与所述第一滤光元件之间。第二会聚元件4-3布置在所述第二路径上、位于所述第二点探测器与所述第二滤光元件之间。第一点探测器3-4位于所述第一会聚元件的光焦点处。第二点探测器4-4位于第二会聚元件的光焦点处。其产生的技术效果与基于光辐射的测温设备的第二种实施例的技术效果相同，在此不再赘述。
图15所示的基于压缩感知的双波长三维温度场成像设备的第四种实施例为：
该实施例是在上述第三种实施例的基础上增加了布置在第一路径上、第一滤光元件3-1和第一会聚元件3-3之间的第一光强衰减元件3-2，和布置在第二路径上、第二滤光元件4-1和第二会聚元件4-3之间的第二光强衰减元件4-2。
由此，本装置通过设置光强衰减元件(例如中性密度滤光片)可以衰减光线的光电信号参量，以防止强光对后续点探测器的损伤，提高了测温精度，延长了产品的使用寿命。
图15所示的基于压缩感知的双波长三维温度场成像设备的第五种实施例为：
该实施例是在上述第四实施例的基础上增加了与第一点探测器3-4和温度确定装置5连接的第一存储装置3-5，以及与第二点探测器4-4和温度确定装置5连接的第二存储装置4-5，它们分别用于存储每次测量到的第一光电信号参量和第二光电信号参量，从而减轻将温度确定装置5作为存储装置时对其造成的存储负担。当时幅变换器11连接到第一点探测器3-4或第二点探测器4-4一侧时，第一存储装置3-5或第二存储装置4-5还可以用来储存时间戳信息。
图15所示的基于压缩感知的双波长三维温度场成像设备的第六种实施例为：
在上述各实施例的基础上增加透镜1，使得待测对象的光辐射进入主光路，可以射向光辐射调制装置2’。
图15所示的基于压缩感知的双波长三维温度场成像设备的第七种实施例为：
在上述各实施例的基础上增加望远镜单元9。
在本发明的一些实施方式中，所述望远镜单元可以是反射式望远镜、折射式望远镜、折反式望远镜中的任意一种。
在本发明的一些实施方式中，所述望远镜单元可以是伽利略望远镜、开普勒望远镜、牛顿望远镜、格里望远镜、卡塞格林望远镜、施密特-卡塞格林望远镜、马克苏托夫-卡塞格林望远镜、多镜面望远镜、双筒望远镜、空间天文望远镜中的任意一种。
本实施例中望远镜单元9采用反射式望远镜，其包括布置于光学成像透镜1和被测物体之间的凹面反射镜9-2和凸面反射镜9-1。凹面反射镜9-2位于凸面反射镜9-1的后方，且凸面反射镜9-1的凸面背向被测物体，凹面反射镜9-2的凹面朝向被测物体。
参照图15，下面介绍一优选实施方式的基于压缩感知的双波长温度场成像设备的工作方式，其具体可以是：
空间光调制元件2-1(DMD)采用掩膜为1和0数量相等的0-1矩阵(固定一帧不变)；
以脉冲光源10-1照射待测对象7；
脉冲光经待测对象7表面反射至空间光调制元件2-1(DMD)。空间光调制元件2-1(DMD)将接收的脉冲光平均等分(将光子数和光辐射的能量平均分配)反射为双臂方向上的两束脉冲光。
时幅变换器11记录第一束聚集光子团最初到达第一探测装置3-4(或第二探测装置4-4)时刻作为初始时间节点t1，之后依次记录第二束聚集光子团最初到达第一探测装置3-4的时刻t2、第三束聚集光子团最初到达第一探测装置3-4的时刻t3、直至第r(r≥1)束聚集光子团最初到达第一探测装置3-4的时刻tr。
第二束聚集光子团最初到达第一探测装置3-4时刻与第一束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔t2-t1；第三束聚集光子团最初到达第一探测装置3-4的时刻与第二束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔t3-t2……第r束聚集光子团最初到达第一探测装置3-4的时刻与第r-1束聚集光子团最初到达第一探测装置3-4的时刻的差值即为两者到达时间的间隔tr-tr-1，因此，每一束聚集光子团出现峰值的时刻，即时间戳信息为：
{(t1+t2)/2，(t2+t3)/2，……，(tr-1+tr)/2}
以脉冲光源10-1继续照射待测对象7，脉冲光经待测对象7反射后可以经过凸面反射镜9-1和凹面反射镜9-2到达透镜1，再射向光辐射调制装置2’。光辐射调制装置2’加载预设的多个掩膜(所述掩根据符合压缩感知理论的RIP条件的测量矩阵变换生成)，并且将接收的光辐射调制反射为双臂方向上的两束光辐射。在双臂中的左臂方向上分别布置有第一滤光元件3-1、第一光强衰减元件3-2、第一会聚元件3-3和第一探测装置3-4。第一束光辐射经过第一滤光元件3-1过滤为波长为第一波长的光，波长为第一波长的光经过第一光强衰减元件3-2进行光强衰减，再经过第一会聚元件3-3进行会聚，第一探测装置3-4在会聚的波长为第一波长的光焦点所在的焦平面处接收波长为第一波长的光，且在时间戳信息{(t1+t2)/2，(t2+t3)/2，……，(tr-1+tr)/2}下探测波长为第一波长的光的光电信号参量。第一存储装置3-5存储第一光电信号参量。
本设备左右对称设置。同样，在双臂中的右臂方向上分别布置有第二滤光元件4-1、第二光强衰减元件4-2、第二会聚元件4-3和第二探测装置4-4。第二束光辐射经过第二滤光元件4-1过滤为波长为第二波长的光，波长为第一波长的光经过第二光强衰减元件4-2进行光强衰减，再经过第二会聚元件4-3进行会聚，第二探测装置4-4在会聚的波长为第二波长的光焦点所在的焦平面处接收波长为第二波长的光，且在时间戳信息{(t1+t2)/2，(t2+t3)/2，……，(tr-1+tr)/2}下探测波长为第二波长的光的光电信号参量。第二存储装置4-5存储第二光电信号参量。
温度确定装置5分别接收第一存储装置3-5存储的和第二存储装置4-5存储的波长为第一波长λ1的光的光电信号参量和波长为第二波长λ2的光的光电信号参量，并利用压缩感知算法重构出待测对象在两个波长λ1和λ2下、时间戳信息为{(t1+t2)/2，(t2+t3)/2，……，(tr-1+tr)/2}时的多个二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]。
温度确定装置5根据同一个时间戳信息[(tr-1+tr)/2]下的二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]反演出时间戳信息为[(tr-1+tr)/2]的待测对象的二维图像S[(tr-1+tr)/2]。
图像生成装置8反演出多个待测对象的二维红外热图像，然后根据根据多个待测对象的二维红外热图像生成待测对象三维红外热图像。
在一些实施例中，为了使过滤后的单波的波长不同，第一滤光元件和所述第二滤光元件为中心波长相差10nm以上的第一窄带滤光片和第二窄带滤光片，所述第一窄带滤光片和第二窄带滤光片的半高宽参数为10nm以上。
因为根据窄带带宽理论，窄带滤光片越窄效果越好；中心波长越接近越好。但是，带宽越窄，与窄带滤光片配合使用的光电信号参量探测元件能探测到的热辐射越少，外加之中心波长接近，对探测器的灵敏度和准确性提出过高的要求，而且此时环境噪声和探测器本征噪声(如暗计数等)也将对系统的测温准确性造成较大影响。所以，在衡量上述优缺点的基础上，经过大量的试验，选用的窄带滤光片的半高宽FWHM一般应在10nm以上，中心波长CWL一般相差10nm以上效果最佳。另外，本设备可以利用不同滤光元件(例如窄带滤光片)进行滤波得到单波光，再将该单波光会聚收集成以供探测器检测，降低了无关光的干扰，提高了采集温度的精度。
在一些实施例中，所述光辐射为红外线光辐射。
在一些实施例中，所述第一和第二点探测器选自雪崩二极管、光电倍增管、超导单光子探测器等单光子点探测器中的任意一种实现。其可以将探测器的噪声限制在暗计数水平，仅受光场的泊松不确定性的影响。
在一些实施例中，所述光电信号参量为光子数。
由此，点探测器可以自由选择各种类型以满足各种需求，增强了产品的通用性能，也便于后期的维护保养。
图16所示的基于压缩感知的双波长三维温度场成像系统的第一种实施例为：
参照图16，该成像系统包括：上述的基于压缩感知的双波长三维温度场成像设备和定标装置6。
定标装置6可以包括温度可调节的基准光源6-1(例如，可以提供不同功率的灯泡)、扩束准直透镜6-2和分束器6-3。分束器6-3设置在待测对象7和透镜1之间。基准光源6-1、扩束准直透镜6-2和分束器6-3水平共线设置，扩束准直透镜6-2设置在基准光源6-1和分束器6-3之间。其中，基准光源6-1用于在定标阶段将调节的不同温度的光辐射射向光辐射调制装置2’，并由光辐射调制装置2’将接收到的光辐射等分为第一光辐射和第二光辐射，使第一光辐射沿第一路径射出(如左臂方向路径)、第二光辐射沿第二路径(如右臂方向路径)射出，以获取多个第一和第二光电信号参量，并确定调节的不同温度和获取的多个第一和第二光电信号参量的预定关系。扩束准直透镜6-2用于将基准光源的光辐射转化为平行光辐射。分束器6-3用于将扩束准直透镜转化的平行光辐射射向所述光辐射调制装置2’。
图16所示的基于压缩感知的双波长温度场成像系统的第二种实施例为：
在图16所示的第一实施例的基础上减少了分束器6-3，该分束器的功能随之减少。但此实施方式仍能解决技术问题，实现相应的技术效果。
图16所示的基于压缩感知的双波长温度场成像系统的第三种实施例为：
在图16所示的第二实施例的基础上减少了扩束准直透镜6-2，该扩束准直透镜的功能随之减少。但此实施方式仍能解决技术问题，实现相应的技术效果。
在上述实施方式中，所述确定调节的不同温度和获取的多个第一和第二光电信号参量的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量。
由此，本系统可以采用不同功率的标准光源进行试验测温，通过测得多组数据来确定试验中的拟合系数。其中，拟合算法利用已知试验或者真实数据，然后寻找一个模型对其规律进行模拟的过程中，求取模型中未知参数的一个过程。以确保在后期的实际测温过程中根据该拟合系数获取测得的温度。本实施方式通过多次试验测进行试验测温，提高了实际测温的精度。此外，本实施方式利用扩束准直透镜6-2将标准光源的光线变换成平行光，减少了光线因为会聚或者散射引起的误差。利用分束器6-3可以最大程度的将标准光源的光线全部传递至透镜，提高了光线的传递效率。
本领域的技术人员可以理解，可以根据实际测量精度或者特殊需要对上述各元件进行选择配置，或者进行多种组合配置。
由此，本系统通过提供定标装置，可以在正式测温之前，进行试验测温，从而制定标准，使得后面的正式测温可以参考试验测温的数据进行运算，调整测得的温度数据，进一度调高了温度的精度。
本实施方式的测温系统的测温的工作方式可以参照上述的测温设备的描述。需要说明的是，在定标的试验之前，将分束器6-3、扩束准直透镜6-2和标准光源2-1移入系统中，保证分束器6-3的反射方向在透镜1和DMD2-1主轴光路上，保证标准光源的光辐射进入系统的主轴光路，且标准光源与待测对象7所在位置共轭，待定标结束后再将分束器6-3、扩束准直透镜6-2和标准光源6-1移出系统。
图17为本发明一种实施方式的基于压缩感知的双波长温度场成像方法的流程示意图。如图17所示，该方法包括以下步骤：
S1701：以脉冲光照射待测对象，记录经待测对象反射的脉冲光的第一束聚集光子团至第r(r≥1)束聚集光子团到达同一地点的时刻，获取时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}；
S1702：继续以脉冲光源照射待测对象，利用光辐射调制装置(可以利用图14和图15实施方式中的光辐射调制装置)接收经待测对象反射的光辐射，加载预设的多个掩膜，而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、多束第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出，所述多个掩膜由符合压缩感知RIP条件的随机测量矩阵Φ变换生成；
S1703：在上述的左臂方向的路径上接收多束第一光辐射，并将其过滤为波长为第一波长的多束光(例如波长为λ1的单波光)；
S1704：在上述的右臂方向的路径上接收多束第二光辐射，并将其过滤为波长为第二波长的多束光(例如波长为λ2的单波光)；
S1705：在上述的左臂方向的路径上接收多束波长为λ1的单波光，并将每一束波长为λ1的单波光转换为多个与时间戳信息对应的第一光电信号参量E1(T)i[(tr-1+tr)/2]；
S1706：在上述的右臂方向的路径上接收多束波长为λ2的单波光，并将每一束波长为λ2的单波光转换为多个与时间戳信息对应的第二光电信号参量E2(T)i[(tr-1+tr)/2]；
S1707：根据多个与时间戳信息对应的第一光电信号参量E1(T)i[(tr-1+tr)/2]和多个与时间戳信息对应的第二光电信号参量E2(T)i[(tr-1+tr)/2]与温度的预定关系确定出所述待测对象在每一个时间戳[(tr-1+tr)/2]下的每一个像素点的温度Ti[(tr-1+tr)/2]；
S1708：根据所述待测对象在时间戳[(tr-1+tr)/2]下的二维图像以及该时间戳下每一个像素点的温度值Ti[(tr-1+tr)/2]反演出所述待测对象每一个横截层面的二维红外热图像；
S1709：根据所述待测对象的每一个横截层面的二维红外热图像反演出所述待测对象的三维红外热图像。
在本实施方式中，所述光辐射调制装置根据预设矩阵Φ生成和加载掩膜的方式如下：
a.提供2k阶的Hadamard矩阵(Hk)，其中且2k-1≥N＝p×q(p为待测对象横坐标像素数、q为待测对象纵坐标像素数)；
b.删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ；
c.将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，获得0-1分布的矩阵
d.以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵获得打乱后的矩阵H+；
e.提取矩阵H+中的M行，并保留该M行中的N列作为随机测量矩阵Φ，其中，M＜＜N/2。
N为待测对象的总像素数，M表示采样数(随机测量矩阵的行)。随机测量矩阵Φ满足RIP条件，适用于压缩感知(CS)重建。
f.加载由随机测量矩阵Φ的第i行拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射，将第一光辐射转换为与时间戳信息对应的多个第一光电信号参量E1(T)2i-1[(tr-1+tr)/2]，将第二光辐射转换为与时间戳信息对应的多个第二光电信号参量E2(T)2i-1[(tr-1+tr)/2]，i＝1，2，3，4…M；
g.加载由随机测量矩阵Φ第i行的补集拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射，将第一光辐射转换为与时间戳信息对应的多个第一光电信号参量E1(T)2i[(tr-1+tr)/2]，将第二光辐射转换为与时间戳信息对应的多个第二光电信号参量E2(T)2i[(tr-1+tr)/2]，i＝1，2，3，4…M。
在本实施方式中，所述所述时间戳信息根据如下方式获得：
利用脉冲光源照射待测对象，利用光辐射调制装置将经待测对象反射的光辐射等分为第一光辐射和第二光辐射，并使第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出；
在上述的左臂方向的路径上接收第一光辐射，并将其过滤为波长为第一波长的光(例如波长为λ1的单波光)；
在上述的右臂方向的路径上接收第二光辐射，并将其过滤为波长为第二波长的光(例如波长为λ2的单波光)；
在上述的左臂方向的路径上、接收波长为λ1的单波光或者在上述右臂方向的路径上、接收波长为λ2的单波光，并记录探测到第一束聚集光子团至第r(r≥1)束聚集光子团的时刻{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}；
在本实施方式中，所述温度确定装置根据如下的预定关系确定待测对象的每个像素点的温度值：
如上所示，待测对象的二维图像表示为：
根据多束波长为λ1的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第一光电信号参量E1(T)i[(tr-1+tr)/2]、结合数学模型Y1 i＝ΦS1[(tr-1+tr)/2]、Y1 i＝E1(T)2i-1[(tr-1+tr)/2]-E1(T)2i[(tr-1+tr)/2]计算出待测对象在时间戳为[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]中的部分Xi(λ1)’[(tr-1+tr)/2]值；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出全部待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1下的二维图像S1中的全部信号值Xi(λ1)[(tr-1+tr)/2]，重构待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]；
根据多束波长为λ2的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第一光电信号参量E2(T)i[(tr-1+tr)/2]、结合数学模型 计算出待测对象在时间戳为[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]中的部分Xi(λ2)’[(tr-1+tr)/2]值；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出全部待测对象在时间戳信息[(tr-1+tr)/2]、波长λ2下的二维图像S2中的全部信号值Xi(λ2)[(tr-1+tr)/2]，重构待测对象在时间戳信息[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]；
根据待测对象在时间戳信息[(tr-1+tr)/2]、波长λ1和波长λ2下的二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]、Xi[(tr-1+tr)/2]＝S1[(tr-1+tr)/2]./S2[(tr-1+tr)/2]以及所述Ti[(tr-1+tr)/2]＝B’/ln(Xi[(tr-1+tr)/2]/A’)确定出待测对象每个像素点的温度值Ti[(tr-1+tr)/2]，以生成时间戳信息为[(tr-1+tr)/2]的待测物体的横截层面的二维红外热图像；其中，A’、B’为预设系数。
根据多个时间戳信息下的待测物体二维红外热图像，反演出待测物体的三维红外热图像。
在一些实施方式中，该方法还包括：
将所述波长为第一波长的光会聚于第一焦点，并在所述第一焦点处设置第一点探测装置用以接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量；
将所述波长为第二波长的光会聚于第二焦点，并在所述第二焦点处设置第二点探测装置用以接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量。
在一些实施方式中，该方法还包括：
对所述波长为第一波长的光的强度进行衰减；及对所述波长为第二波长的光的强度进行衰减。
在一些实施方式中，该方法还包括：
利用望远镜单元接收待测对象的光辐射。
在本实施方式中，所述光辐射为红外线光辐射。
在本实施方式中，所述第一波长与所述第二波长相差至少10nm。
在一些实施例中，所述第一和第二点探测器选自雪崩二极管、光电倍增管、超导单光子探测器等单光子点探测器中的任意一种实现。其可以将探测器的噪声限制在暗计数水平，仅受光场的泊松不确定性的影响。
在一些实施例中，所述光电信号参量为光子数。
在测温方法上述实施例中的技术效果与测温装置实施例中的技术效果相对应，在此不再赘述。
本发明一种实施方式的基于压缩感知的双波长温度场成像方法同样可以包括如图7所示的定标的步骤，即在步骤S1701之前还包括：
S1801：利用光辐射调制装置接收基准光源(例如温度可调节的灯泡，具体可以调节电流、电压等方式来调节)的光辐射，将接收到的光辐射等分(此处等分的概念可以参考上述测温设备中等分的概念)为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、所述第二光辐射沿不同于第一路径的第二路径(如测温设备的右臂方向的路径)射出。
S1802：在所述第一路径上接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光(例如波长为λ1的单波光)。
S1803：在所述第二路径上接收所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光(例如波长为λ2的单波光)。
S1804：在所述第一路径上接收所述波长为λ1的单波光并将其转换为相应的第一光电信号参量E1(T)。
S1805：在所述第二路径上接收所述波长为λ2的单波光并将其转换为相应的第二光电信号参量E2(T)。
S1806：调节灯泡的电流或电压，使得灯泡发出光辐射的温度变为Td2，并获取对应的光电信号参量Ed1(T)2和Ed2(T)2。按此种方法可以获取灯泡发出多个不同温度Tdi的光辐射，并获取对应的多个第一和第二光电信号参量Ed1(T)i和Ed2(T)i。
S1807：根据上述不同温度与第一和第二光电信号参量的计量值，确定出所述第一和第二光电信号参量与温度的预定关系。
在本实施方式中，所述确定多个第一和第二光电信号参量与调节的多个不同温度的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，所述多个为n个，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量。
具体实现方式可以是：
在定标时，标准光源以第一电流、电压、电阻的光线照射扩束准直透镜，通过分束器、透镜到达DMD。保持DMD同一帧光辐射不变，其加载的掩膜为1和0数量相等的0-1矩阵。控制元件控制DMD中各微镜的闭合、翻转，使得DMD将接收到的一帧光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出。其中：第一光辐射经过第一窄带滤光片、第一会聚元件到达第一点探测器，得到第一光电信号参量Ed1(T)1。第二光辐射经过第二窄带滤光片、第二会聚元件到达第二点探测器，得到第二光电信号参量Ed2(T)1。Ed1(T)1和Ed2(T)1经过除法器得到比值X1，X1＝Ed1(T)1/Ed2(T)1；
调节标准光源的不同等效电流或电压或电阻可计算出不同电流或电压或电阻下的一系列温度T1，T2，T3……Tdi，通过本系统测得对应的Ed1(T)i与Ed2(T)i比值X1，X2，X3……Xi，保持DMD固定一帧不变，其DMD上加载的掩膜为1和0数量相等的0-1矩阵。接着利用公式拟合得出相应系数，即得到此种环境下的辐射体(即基准光源)的光辐射温度曲线公式。即：确定拟合系数(即上述的预设系数)A’和B’，其中：
Ed1(T)1/Ed2(T)1＝A’exp(B’/Td1)
Ed1(T)2/Ed2(T)2＝A’exp(B’/Td2)
Ed1(T)3/Ed2(T)3＝A’exp(B’/Td3)
……
Ed1(T)＝A1(λ1)exp(B1(λ1)/T)
Ed2(T)＝A2(λ2)exp(B2(λ2)/T)
A’＝A1(λ1)/A2(λ2)，B’＝B1(λ1)-B2(λ2)。
具体可以是利用公式拟合得出相应系数A’和B’，得到此种环境下的该辐射体的光辐射温度曲线公式；然后通过计算机程序设定好此系数下的计算公式，即可测量或者监测任何情况下的辐射温度。在实际测温过程中，利用指数形式的多项式展开，获得多项式的系数，或者利用其它复杂的拟合方程获得拟合系数。
以上对本发明的各种实施例进行了详细说明。本领域技术人员将理解，可在不偏离本发明范围(由所附的权利要求书限定)的情况下，对实施方案进行各种修改、改变和变化。对权利要求范围的解释应从整体解释且符合与说明一致的最宽范围，并不限于示例或详细说明中的实施范例。
Claims (33)
1.一种基于压缩感知的双波长三维温度场成像设备，其特征在于，包括：
脉冲光源，配置为向待测物体发射脉冲光；
光辐射调制装置，配置为接收经由待测物体的反射的光辐射，且加载预设的多个掩膜，将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿不同于第一路径的第二路径射出，所述多个掩膜根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成；
布置在所述第一路径上的第一滤光元件，配置为接收多束所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长λ1的多束光；
布置在所述第二路径上的第二滤光元件，配置为接收多束所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长λ2的多束光；
布置在第一路径上的第一探测装置，配置为接收波长为第一波长λ1的多束光并将其转换为第一光电信号参量；
布置在第二路径上的第二探测装置，配置为接收波长为第二波长λ2的多束光并将其转换为第二光电信号参量；
布置在第一路径或者第二路径上的时幅变换器，配置为记录所述第一波长λ1的光或者第二波长λ2的光中的多束光子团到达第一探测装置或第二探测装置的一系列时间戳信息；
所述第一探测装置将波长为第一波长λ1的每一束光转换为多个与所述时间戳信息对应的第一光电信号参量；
所述第二探测装置将波长为第二波长λ2的每一束光转换为多个与所述时间戳信息对应的第二光电信号参量；
温度确定装置，配置为接收来自所述第一探测装置和第二探测装置的多个与所述时间戳信息对应的所述第一和第二光电信号参量，并根据多个与所述时间戳信息对应的所述第一和第二光电信号参量与温度的预定关系确定出所述待测物体在每一个所述时间戳信息时的每一个像素点的温度值；
图像生成装置，配置为根据所述待测物体在每一个所述时间戳信息时的每一个像素点的温度值以及待测物体在每一个所述时间戳信息时的二维图像反演出待测物体的每一个横截层面的二维红外热图像，且根据待测物体的每一个横截层面的二维红外热图像生成三维红外热图像。
2.根据权利要求1所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述记录所述第一波长λ1的光或者第二波长λ2的光中的多束光子团到达第一探测装置或第二探测装置的一系列时间戳信息，包括：
利用脉冲光源照射待测物体，利用所述光辐射调制装置将经所述待测物体反射的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿第二路径射出；
在所述第一路径上接收所述第一光辐射，并将其过滤为波长为第一波长的单波光；
在所述第二路径上接收所述第二光辐射，并将其过滤为波长为第二波长的单波光；
在所述第一路径上、接收所述波长为第一波长λ1的单波光或者在所述第二路径上、接收所述波长为第二波长λ2的单波光，并记录所述第一波长λ1的单波光或者所述第二波长λ2的单波光中的第一聚集光子团至第r(r≥1)聚集光子团到达第一探测装置或第二探测装置的时刻{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}，其中，r为自然数。
3.根据权利要求2所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述光辐射调制装置加载预设的多个掩膜，所述第一探测装置将波长为第一波长λ1的每一束光转换为多个与所述时间戳信息对应的第一光电信号参量，所述第二探测装置将波长为第二波长λ2的每一束光转换为多个与所述时间戳信息对应的第二光电信号参量，包括：
提供2k阶的Hk矩阵，其中2k-1≥N，N为所述待测物体的总像素数；
删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ；
将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，获得0-1分布的矩阵
以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵得到矩阵H+；
提取上述矩阵H+中M行，并保留M行中的N列作为测量矩阵Φ；
所述光辐射调制装置加载由所述测量矩阵Φ的第i行H+i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
所述第一探测装置将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i-1[(tr-1+tr)/2]，所述第二探测装置将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i-1[(tr-1+tr)/2]；
所述光辐射调制装置加载由所述测量矩阵Φ的第i行的补集H-i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
所述第一探测装置将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i[(tr-1+tr)/2]，所述第二探测装置将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i[(tr-1+tr)/2]；
其中，i＝1,2,3,4……M；M/N×100％＝5％～30％。
4.根据权利要求3所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述温度确定装置接收来自所述第一探测装置和第二探测装置的多个与所述时间戳信息对应的所述第一和第二光电信号参量，并根据多个与所述时间戳信息对应的所述第一和第二光电信号参量与温度的预定关系确定出所述待测物体在每一个所述时间戳信息时的每一个像素点的温度值，包括：
所述温度确定装置根据多束波长为λ1的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第一光电信号参量E1(T)i[(tr-1+tr)/2]、结合数学模型以及计算出待测物体在时间戳为[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]中的部分信号值Xi(λ1)’[(tr-1+tr)/2]；
所述温度确定装置利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测物体在时间戳信息[(tr-1+tr)/2]、波长λ1下的的全部信号值Xi(λ1)[(tr-1+tr)/2]，重构待测物体在时间戳信息[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]；
所述温度确定装置根据根据多束波长为λ2的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第二光电信号参量E2(T)i[(tr-1+tr)/2]、结合数学模型以及计算出待测物体在时间戳为[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]中的部分信号值Xi(λ2)’[(tr-1+tr)/2]；
所述温度确定装置利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测物体在时间戳信息[(tr-1+tr)/2]、波长λ2下的的全部信号值Xi(λ2)[(tr-1+tr)/2]，重构待测物体在时间戳信息[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]；
所述温度确定装置根据待测物体在时间戳信息[(tr-1+tr)/2]、波长λ1和波长λ2下的二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]、结合Xi[(tr-1+tr)/2]＝S1[(tr-1+tr)/2]./S2[(tr-1+tr)/2]以及Ti[(tr-1+tr)/2]＝B’/ln(Xi[(tr-1+tr)/2]/A’)确定出待测物体每个像素点的温度值Ti[(tr-1+tr)/2]，其中，A’、B’为预设系数。
5.根据权利要求4所述的基于压缩感知的双波长三维温度场成像设备，其特征在于：
所述稀疏基Ψ采用梯度基、离散余弦变换基、小波基、傅里叶变换基、Gabor变换基中的任意一种实现。
6.根据权利要求5所述的基于压缩感知的双波长三维温度场成像设备，其特征在于：
所述待测物体在波长λ1和波长λ2下的二维图像S1和S2的重构可以采用MP、OMP、StOMP、ReOMP、SPARCO、LARS、BCS、FBMP、SL0、KF-CS、BPDQ、k-t FOCUSS、l1-Magic、IST、TwIST、FISTA、SOCP、RecPF、FPC、GPSR、SpaRSA、ALM、TVAL3、CoSaMP、LBI、SP、l1_ls、SMP、TFOCS、NESTA、SALSA、SPGL1、YALL1、l0重建算法、l1重建算法、l2重建算法、BPDN算法、LASSO算法、TV算法中的任意一种实现。
7.根据权利要求1～6中任一项所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述光辐射调制装置包括：
空间光调制器，配置为根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成的多个掩膜，以将接收到的待测物体的光辐射调制为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出；
控制元件，配置为控制所述空间光调制器依次加载由所述测量矩阵Φ变换生成的多个掩膜。
8.根据权利要求7所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述空间光调制器为选自数字微镜器件、光强数字调制器或液晶光阀。
9.根据权利要求7所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述第一探测装置为第一点探测器，所述第二探测装置为第二点探测器，并且
所述基于压缩感知的双波长三维温度场成像设备还包括布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一会聚元件，和布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二会聚元件，
所述第一点探测器位于所述第一会聚元件的光焦点处；
所述第二点探测器位于所述第二会聚元件的光焦点处。
10.根据权利要求9所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，还包括：
布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一光强衰减元件，和
布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二光强衰减元件。
11.根据权利要求1～6中任一项所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述第一滤光元件和所述第二滤光元件为中心波长相差至少10nm的第一窄带滤光片和第二窄带滤光片，所述第一窄带滤光片和第二窄带滤光片的半高宽参数至少为10nm。
12.根据权利要求1～6中任一项所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述光电信号参量包括光子数。
13.根据权利要求1～6中任一项所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，所述光辐射为红外波段的光辐射。
14.根据权利要求13所述的基于压缩感知的双波长三维温度场成像设备，其特征在于，
所述第一和第二点探测器选自雪崩二极管或光电倍增管或超导单光子探测器。
15.一种基于压缩感知的双波长三维温度场成像系统，其特征在于，包括：
如权利要求1～14中任一项所述基于压缩感知的双波长三维温度场成像设备以及定标装置，
所述定标装置包括温度可调节的基准光源，所述基准光源配置为在一个定标阶段中，将调节成的不同温度的光辐射射向所述光辐射调制装置，由所述光辐射调制装置将接收到的光辐射等分为第一光辐射和第二光辐射，使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出，以获取多个第一和第二光电信号参量，并根据所述不同温度与所述第一和第二光电信号参量的计量值，确定出所述第一和第二光电信号参量与温度的预定关系。
16.根据权利要求15所述的基于压缩感知的双波长三维温度场成像系统，其特征在于，所述第一和第二光电信号参量与温度的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，i为1至n的自然数；A’、B’为预设系数，Tdi为所述基准光源发出的第i个温度，Ed1(T)i为第i次定标时测得的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标时测得的波长为第二波长的光的光电信号参量。
17.根据权利要求15或16所述的基于压缩感知的双波长三维温度场成像系统，其特征在于，所述定标装置还包括将基准光源的光辐射转化为平行光辐射的扩束准直透镜。
18.根据权利要求17所述的基于压缩感知的双波长三维温度场成像系统，其特征在于，所述定标装置还包括将扩束准直透镜转化的平行光辐射射向所述光辐射调制装置的分束器。
19.一种基于压缩感知的双波长三维温度场成像方法，其特征在于，包括：
以脉冲光照射待测物体，记录经待测物体反射的脉冲光的第一束聚集光子团至第r(r≥1)束聚集光子团到达同一地点的时刻，获取时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}，其中，r为自然数；
以脉冲光源照射待测物体，利用光辐射调制装置接收经待测物体反射的光辐射，加载预设的多个掩膜，而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿第二路径射出，所述多个掩膜由符合压缩感知RIP条件的测量矩阵Φ变换生成；
在所述第一路径上接收多束所述第一光辐射，并将其过滤为波长为第一波长λ1的多束光；
在所述第二路径上接收多束所述第二光辐射，并将其过滤为波长为第二波长λ2的多束光；
在所述第一路径上接收多束所述波长为λ1的单波光，并将每一束波长为λ1的单波光转换为多个与时间戳信息对应的第一光电信号参量E1(T)i[(tr-1+tr)/2]；
在所述第二路径上接收多束所述波长为λ2的单波光，并将每一束波长为λ2的单波光转换为多个与时间戳信息对应的第二光电信号参量E2(T)i[(tr-1+tr)/2]；
根据多个与时间戳信息对应的所述第一光电信号参量E1(T)i[(tr-1+tr)/2]和多个与时间戳信息对应的所述第二光电信号参量E2(T)i[(tr-1+tr)/2]与温度的预定关系确定出所述待测物体在每一个时间戳[(tr-1+tr)/2]下的每一个像素点的温度Ti[(tr-1+tr)/2]；
根据所述待测物体在每一个时间戳[(tr-1+tr)/2]下的二维图像S[(tr-1+tr)/2]以及该时间戳下每一个像素点的温度值Ti[(tr-1+tr)/2]反演出所述待测物体的每一个横截层面的二维红外热图像；
根据所述待测物体的每一个横截层面的二维红外热图像反演出所述待测物体的三维红外热图像。
20.根据权利要求19所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述以脉冲光照射待测物体，记录经待测物体反射的脉冲光的第一束聚集光子团至第r(r≥1)束聚集光子团到达同一地点的时刻，获取时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}包括：
利用脉冲光源照射待测物体，利用光辐射调制装置将经待测物体反射的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿第二路径射出；
在所述第一路径上接收所述第一光辐射，并将其过滤为波长为第一波长的单波光；
在所述第二路径上接收所述第二光辐射，并将其过滤为波长为第二波长的单波光；
在所述第一路径上接收所述波长为第一波长的单波光或者在所述第二路径上接收波长为第二波长的单波光，并记录探测到所述第一波长的单波光或者所述第二波长的单波光中包括的第一聚集光子团至第r聚集光子团的时刻作为时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}。
21.根据权利要求20所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述光辐射调制装置加载预设的多个掩膜，而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射，并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿第二路径射出，所述多个掩膜由符合压缩感知RIP条件的测量矩阵Φ变换生成；在所述第一路径上接收多束所述波长为λ1的单波光，并将每一束波长为λ1的单波光转换为多个与时间戳信息对应的第一光电信号参量E1(T)i[(tr-1+tr)/2]；在所述第二路径上接收多束所述波长为λ2的单波光，并将每一束波长为λ2的单波光转换为多个与时间戳信息对应的第二光电信号参量E2(T)i[(tr-1+tr)/2]包括：
提供2k阶的Hk矩阵，其中2k-1≥N，N为所述待测物体的总像素数；
删除矩阵Hk的第一行和第一列后得到(2k-1)×(2k-1)阶的剩余矩阵HΩ；
将矩阵HΩ中所有元素中的-1变为0，所有1保持不变，获得0-1分布的矩阵
以两个2k-1阶的扩充排列矩阵P和Q打乱2k-1阶0-1矩阵得到矩阵H+；
提取上述矩阵H+中M行，并保留M行中的N列作为测量矩阵Φ；
所述光辐射调制装置加载由所述测量矩阵Φ的第i行H+i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
在所述第一路径上接收多束所述波长为λ1的单波光，将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i-1[(tr-1+tr)/2]，在所述第二路径上接收多束所述波长为λ2的单波光，将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i-1[(tr-1+tr)/2]；
所述光辐射调制装置加载由所述测量矩阵Φ的第i行的补集H-i拉伸变换而得的掩膜，并且将接收到的光辐射调制为第一光辐射和第二光辐射；
在所述第一路径上接收多束所述波长为λ1的单波光，将波长为第一波长λ1的每一束光转换为与所述时间戳信息对应的第一光电信号参量E1(T)2i[(tr-1+tr)/2]，在所述第二路径上接收多束所述波长为λ2的单波光，将波长为第二波长λ2的每一束光转换为与所述时间戳信息对应的第二光电信号参量E2(T)2i[(tr-1+tr)/2]；
其中，i＝1,2,3,4……M；M/N×100％＝5％～30％。
22.根据权利要求21所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述根据多个与时间戳信息对应的所述第一光电信号参量E1(T)i[(tr-1+tr)/2]和多个与时间戳信息对应的所述第二光电信号参量E2(T)i[(tr-1+tr)/2]与温度的预定关系确定出所述待测物体在每一个时间戳[(tr-1+tr)/2]下的每一个像素点的温度Ti[(tr-1+tr)/2]包括：
根据多束波长为λ1的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第一光电信号参量E1(T)i[(tr-1+tr)/2]、结合数学模型以及计算出待测物体在时间戳为[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]中的部分信号值Xi(λ1)’[(tr-1+tr)/2]；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测物体在时间戳信息[(tr-1+tr)/2]、波长λ1下的的全部信号值Xi(λ1)[(tr-1+tr)/2]，重构待测物体在时间戳信息[(tr-1+tr)/2]、波长λ1下的二维图像S1[(tr-1+tr)/2]；
根据根据多束波长为λ2的单波光在相同的时间戳信息[(tr-1+tr)/2]下的第二光电信号参量E2(T)i[(tr-1+tr)/2]、结合数学模型以及计算出待测物体在时间戳为[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]中的部分信号值Xi(λ2)’[(tr-1+tr)/2]；
利用压缩感知原理Y＝ΦΨx’+e，在已知测量值向量Y、测量矩阵Φ和稀疏基Ψ的情况下，反演出待测物体在时间戳信息[(tr-1+tr)/2]、波长λ2下的的全部信号值Xi(λ2)[(tr-1+tr)/2]，重构待测物体在时间戳信息[(tr-1+tr)/2]、波长λ2下的二维图像S2[(tr-1+tr)/2]；
根据待测物体在时间戳信息[(tr-1+tr)/2]、波长λ1和波长λ2下的二维图像S1[(tr-1+tr)/2]和S2[(tr-1+tr)/2]、结合Xi[(tr-1+tr)/2]＝S1[(tr-1+tr)/2]./S2[(tr-1+tr)/2]以及Ti[(tr-1+tr)/2]＝B’/ln(Xi[(tr-1+tr)/2]/A’)确定出待测物体每个像素点的温度值Ti[(tr-1+tr)/2]，其中，A’、B’为预设系数。
23.根据权利要求22所述的基于压缩感知的双波长三维温度场成像方法，其特征在于：
所述稀疏基Ψ采用梯度基、离散余弦变换基、小波基、傅里叶变换基、Gabor变换基中的任意一种实现。
24.根据权利要求23所述的基于压缩感知的双波长三维温度场成像方法，其特征在于：
所述待测物体在波长λ1和波长λ2下的二维图像S1和S2的重构可以采用MP、OMP、StOMP、ReOMP、SPARCO、LARS、BCS、FBMP、SL0、KF-CS、BPDQ、k-t FOCUSS、l1-Magic、IST、TwIST、FISTA、SOCP、RecPF、FPC、GPSR、SpaRSA、ALM、TVAL3、CoSaMP、LBI、SP、l1_ls、SMP、TFOCS、NESTA、SALSA、SPGL1、YALL1、l0重建算法、l1重建算法、l2重建算法、BPDN算法、LASSO算法、TV算法中的任意一种实现。
25.根据权利要求19～24中任一项所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述光辐射调制装置包括：
空间光调制器，配置为根据符合压缩感知理论的RIP条件的测量矩阵Φ变换生成的多个掩膜，以将接收到的待测物体的光辐射调制为所述第一光辐射和第二光辐射，并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出；
控制元件，配置为控制所述空间光调制器依次加载由所述测量矩阵Φ变换生成的多个掩膜。
26.根据权利要求25所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述空间光调制器选自数字微镜器件、光强数字调制器或液晶光阀。
27.根据权利要求19～24中任一项所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，进一步包括：
将所述波长为第一波长的光会聚于第一焦点，并在所述第一焦点处设置第一点探测装置用以接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量；
将所述波长为第二波长的光会聚于第二焦点，并在所述第二焦点处设置第二点探测装置用以接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量。
28.根据权利要求19～24中任一项所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，进一步包括：
对所述波长为第一波长的光的强度进行衰减；及
对所述波长为第二波长的光的强度进行衰减。
29.根据权利要求19～24中任一项所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述光辐射为红外波段的光辐射。
30.根据权利要求19～24中任一项所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述第一波长与所述第二波长相差至少10nm。
31.根据权利要求19～24中任一项所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述光电信号参量包括光子数。
32.根据权利要求19～24中任一项所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述以脉冲光照射待测物体，记录经待测物体反射的脉冲光的第一束聚集光子团至第r(r≥1)束聚集光子团到达同一地点的时刻，获取时间戳信息{(t1+t2)/2，(t2+t3)/2，(t3+t4)/2，……，(tr-1+tr)/2}之前还包括定标步骤，
所述定标步骤包括：
利用光辐射调制装置接收基准光源的光辐射，将接收到的光辐射等分为第一光辐射和第二光辐射，并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出；
在所述第一路径上接收所述第一光辐射，并将接收到的第一光辐射过滤为波长为第一波长的光；
在所述第二路径上接收所述第二光辐射，并将接收到的第二光辐射过滤为波长为第二波长的光；
在所述第一路径上接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量；
在所述第二路径上接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量；
调节基准光源发出多个不同温度的光辐射，并获取对应的多个第一和第二光电信号参量；以及
根据所述不同温度与所述第一和第二光电信号参量的计量值，确定出所述第一和第二光电信号参量与温度的预定关系。
33.根据权利要求32所述的基于压缩感知的双波长三维温度场成像方法，其特征在于，所述第一和第二光电信号参量与温度的预定关系为：
Ed1(T)i/Ed2(T)i＝A’exp(B’/Tdi)，
其中，i为1至n的自然数；
A’、B’为预设系数，Tdi为第i个基准光源的温度，Ed1(T)i为第i次定标的波长为第一波长的光的光电信号参量，Ed2(T)i为第i次定标的波长为第二波长的光的光电信号参量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610080287.2A CN105675146B (zh) | 2016-02-04 | 2016-02-04 | 基于压缩感知的双波长三维温度场成像设备、系统及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610080287.2A CN105675146B (zh) | 2016-02-04 | 2016-02-04 | 基于压缩感知的双波长三维温度场成像设备、系统及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105675146A CN105675146A (zh) | 2016-06-15 |
| CN105675146B true CN105675146B (zh) | 2018-07-17 |
Family
ID=56303625
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610080287.2A Active CN105675146B (zh) | 2016-02-04 | 2016-02-04 | 基于压缩感知的双波长三维温度场成像设备、系统及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105675146B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107846592A (zh) * | 2017-10-18 | 2018-03-27 | 中国科学院计算技术研究所 | 压缩感知图像的采集重建方法和系统 |
| CN109299569B (zh) * | 2018-10-24 | 2023-04-18 | 广州市香港科大霍英东研究院 | 一种基于相干结构的不可压缩黏性流体的大涡模拟方法 |
| CN111650600B (zh) * | 2020-05-28 | 2023-03-31 | 北京空间机电研究所 | 一种面向极微弱信号的双谱段激光成像装置 |
| JP7146146B2 (ja) * | 2020-07-22 | 2022-10-03 | 三菱電機株式会社 | 距離計測装置、距離計測方法、及び工作装置 |
| CN112082672B (zh) * | 2020-08-28 | 2021-08-03 | 上海交通大学 | 基于消光矫正时间解析的三维燃烧场温度测量方法 |
| CN113033723B (zh) * | 2021-03-08 | 2023-06-16 | 山东大学 | 环形掩膜版、光场调控方法、单像素成像方法及系统 |
| CN113506285A (zh) * | 2021-07-27 | 2021-10-15 | 西北工业大学 | 一种锅炉炉膛三维温度场检测方法、装置和计算机设备 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102393248A (zh) * | 2011-10-26 | 2012-03-28 | 中国科学院空间科学与应用研究中心 | 一种时间分辨极弱光多光谱成像系统及方法 |
| CN102510282A (zh) * | 2011-10-25 | 2012-06-20 | 中国科学院空间科学与应用研究中心 | 一种时间分辨单光子计数二维成像系统及方法 |
| CN102608619A (zh) * | 2012-03-07 | 2012-07-25 | 北京航空航天大学 | 光子计数压缩采样相控阵激光三维成像方法 |
| CN103592026A (zh) * | 2013-11-19 | 2014-02-19 | 中国科学院空间科学与应用研究中心 | 基于压缩感知和编码变换的时间飞行成像光谱系统及方法 |
| CN103968945A (zh) * | 2014-05-28 | 2014-08-06 | 中国科学院空间科学与应用研究中心 | 基于二阶压缩感知的超灵敏光谱成像天文望远镜及方法 |
| CN205642635U (zh) * | 2016-02-04 | 2016-10-12 | 北京理工大学 | 基于压缩感知的双波长三维温度场成像设备及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013041949A1 (en) * | 2011-09-20 | 2013-03-28 | Mesa Imaging Ag | Time of flight sensor with subframe compression and method |
-
2016
- 2016-02-04 CN CN201610080287.2A patent/CN105675146B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102510282A (zh) * | 2011-10-25 | 2012-06-20 | 中国科学院空间科学与应用研究中心 | 一种时间分辨单光子计数二维成像系统及方法 |
| CN102393248A (zh) * | 2011-10-26 | 2012-03-28 | 中国科学院空间科学与应用研究中心 | 一种时间分辨极弱光多光谱成像系统及方法 |
| CN102608619A (zh) * | 2012-03-07 | 2012-07-25 | 北京航空航天大学 | 光子计数压缩采样相控阵激光三维成像方法 |
| CN103592026A (zh) * | 2013-11-19 | 2014-02-19 | 中国科学院空间科学与应用研究中心 | 基于压缩感知和编码变换的时间飞行成像光谱系统及方法 |
| CN103968945A (zh) * | 2014-05-28 | 2014-08-06 | 中国科学院空间科学与应用研究中心 | 基于二阶压缩感知的超灵敏光谱成像天文望远镜及方法 |
| CN205642635U (zh) * | 2016-02-04 | 2016-10-12 | 北京理工大学 | 基于压缩感知的双波长三维温度场成像设备及系统 |
Non-Patent Citations (1)
| Title |
|---|
| Measurement dimensions compressed spectral imaging with a single point detector;Xue-Feng Liu等;《Optics Communications》;20151111;第365卷;第173-179页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105675146A (zh) | 2016-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105675146B (zh) | 基于压缩感知的双波长三维温度场成像设备、系统及方法 | |
| CN103472457B (zh) | 稀疏孔径压缩计算关联飞行时间的三维成像系统及方法 | |
| CN103472455B (zh) | 稀疏孔径压缩计算关联飞行时间四维光谱成像系统及方法 | |
| CN103884486B (zh) | 纹影测量成像系统及方法 | |
| CN103453993B (zh) | 基于稀疏孔径压缩计算关联的主动高光谱成像系统及方法 | |
| CN104121990B (zh) | 基于随机光栅的压缩感知宽波段高光谱成像系统 | |
| CN103471718B (zh) | 一种基于稀疏孔径压缩计算关联的高光谱成像系统及方法 | |
| CN103472456B (zh) | 一种基于稀疏孔径压缩计算关联的主动成像系统及方法 | |
| CN105737992B (zh) | 基于压缩感知的双波长温度场成像设备、系统及方法 | |
| Howard et al. | Hyperspectral compressive wavefront sensing | |
| CN103954357B (zh) | 压缩光谱成像系统测量矩阵的获取方法 | |
| De Borniol et al. | A 320x256 HgCdTe avalanche photodiode focal plane array for passive and active 2D and 3D imaging | |
| CN205642634U (zh) | 基于压缩感知的双波长温度场成像设备及系统 | |
| CN108663118B (zh) | 红外宽波段高光谱计算成像装置及其方法 | |
| CN205642635U (zh) | 基于压缩感知的双波长三维温度场成像设备及系统 | |
| CN205642633U (zh) | 基于光辐射的测温设备及系统 | |
| McCullough et al. | Quantum efficiency and quantum yield of an HgCdTe infrared sensor array | |
| CN105606228B (zh) | 基于编码变换的双波长温度场成像设备、系统及方法 | |
| CN205642623U (zh) | 基于编码变换的双波长温度场成像设备及系统 | |
| Volin et al. | Midwave-infrared snapshot imaging spectrometer | |
| CN105527024B (zh) | 基于光辐射的测温设备、系统及方法 | |
| Edgar et al. | 3D computational ghost imaging | |
| Brown et al. | Characterization of Earth observing satellite instruments for response to spectrally and spatially variable scenes | |
| Baron et al. | Hyper-Spectral Speckle Imaging of Resolved Targets | |
| CN118244505A (zh) | 一种宽光谱的单曝光合成孔径成像增强装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |