CN105606228B - 基于编码变换的双波长温度场成像设备、系统及方法 - Google Patents

基于编码变换的双波长温度场成像设备、系统及方法 Download PDF

Info

Publication number
CN105606228B
CN105606228B CN201610080112.1A CN201610080112A CN105606228B CN 105606228 B CN105606228 B CN 105606228B CN 201610080112 A CN201610080112 A CN 201610080112A CN 105606228 B CN105606228 B CN 105606228B
Authority
CN
China
Prior art keywords
optical radiation
wavelength
light
temperature
path
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
Application number
CN201610080112.1A
Other languages
English (en)
Other versions
CN105606228A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201610080112.1A priority Critical patent/CN105606228B/zh
Publication of CN105606228A publication Critical patent/CN105606228A/zh
Application granted granted Critical
Publication of CN105606228B publication Critical patent/CN105606228B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01J5/10Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors

Landscapes

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

Abstract

本申请提供了一种基于编码变换的双波长温度场成像设备,包括光辐射调制装置,配置为接收待测对象的光辐射,且加载根据预设矩阵Φ变换生成的多个掩膜,将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿不同于第一路径的第二路径射出;布置在所述第一路径上的第一滤光元件和第一探测装置;布置在所述第二路径上的第二滤光元件和第二探测装置;温度确定装置以及图像生成装置。本申请还提供了基于该设备的系统和方法。本发明将双波长测温技术、调制技术和编码技术相结合,重建出待测对象的二维红外热图像,其能够广泛地应用在深空探测、遥感、材料检测、夜视观测等相关科技领域。

Description

基于编码变换的双波长温度场成像设备、系统及方法
技术领域
本发明涉及双波长温度场成像领域,特别涉及一种基于编码变换的双波长温度场成像设备、系统及方法。
背景技术
在航空航天、冶金及汽车制造业等领域,常常需要对待测对象以及各种在线工件进行快速、实时监控,以最大限度地减少事故隐患、提高产品的安全性能及其质量。沿用传统的接触式测温仪进行测量,虽然精度高,但必须使探测器接触待测对象。可是在某些特殊的场合下(如对发动机燃烧室及高温炉中的火焰进行测温时)无法使用接触式测温仪,由此产生了非接触式测温方法。红外测温方法就属于一种非接触式测温方法,该方法通过检测物体表面发射的能量来测量温度,具有测温范围广、响应速度快和不明显破坏测温度场等特点,被广泛应用于工业各个方面。
红外测温方法主要基于黑体辐射理论,黑体是个理想化的物理模型,而自然界中实际存在的物体(测温对象),其吸收能力及辐射能力都比黑体小,称为灰体。根据普朗克辐射定律,一个绝对温度为T的黑体,单位表面积在波长λ1、λ21的附近单位波长)的间隔内向整个半球空间发射的辐射功率(简称光谱辐射度)为E0(λ,T),而灰体的光谱辐射能量的计算公式为:E(λ,T)=ε(λ,T)E0(λ,T),其中ε(λ,T)为该灰体的辐射率。
现有技术中红外测温主要经历了三个阶段的发展。
第一阶段:传统的红外测温设备,一律按黑体的热辐射定律来设计。该设计方式假设红外测温仪实际接收到的热辐射是与被测物的光谱辐射能量E(λ,T)是成比例,故在使用红外测温仪时,必须求出待测对象的辐射率ε(λ,T)数值,即进行辐射率修正。遗憾的是,该辐射率ε(λ,T)与待测对象的材料、表面状态、波长,温度以及辐射条件、环境因素等均有复杂的关系,因而很难准确测定ε(λ,T),同时由于在某些情况下待测对象的辐射率随温度变化太大,因而传统红外测温仪还存在着较大的误差。
第二阶段:科学家们为解决传统的红外测温设备存在的问题,研究出了基于单波长窄带滤波红外测温技术。可是由于存在诸如水蒸气等的周围环境对红外线的吸收,在很大幅度上影响了单波长红外测温的准确性。
第三阶段:利用双波长滤波红外测温技术来进行红外测温。双波长滤波红外测温技术的原理为:利用黑体辐射曲线中相邻两个波长对应的能量等比吸收的原理,在保证红外测温高精度测量的基础上,克服了环境对物体发射红外线吸收造成的测量误差。根据普朗克辐射定律,一个绝对温度为T的黑体,单位表面积在波长λ1、λ21的附近单位波长)的间隔内向整个半球空间发射的辐射功率(简称光谱辐射度)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)=A0 exp(B0/T)
其中,A0=C1λ-5,B0=-C2/λ,且仍然只适用于黑体。若将A0和B0视为可变参量A和B,则可推广到灰体的情况,则灰体的光谱辐射能量为:
E(T)=Aexp(B/T)。
不同于公式的是,前者只需简单改变A和B参量的值就可实现从黑体到灰体的修正,而无需确定复杂的辐射率函数ε(λ,T)。
利用相邻两波长等比吸收的原理。取2个波长的吸收能量比作为温度的函数,即可避免由于水蒸气等环境对红外线吸收的因素而造成的测量误差。
现在分别取λ1和λ2,则有:
E1(T)=A11)exp(B11)/T),
E2(T)=A22)exp(B22)/T)。
上述两式作比值得:
其中,因此,只要利用实验数据拟合确定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的绝对等分的差距很大,所以造成后期测温精度下降。此外,双波长滤波红外测温设备中的二向色镜依赖在光学平片的正反面镀上不同的膜以实现滤光,如果需要选择其它的波长,则需要更换整块二向色镜。所以双波长滤波红外测温设备测得的温度的精度不高,且使用不太方便,其适应性不广。
尽管近年来科学家也提出了多波长测温方案,用于提高测温精度,但其结构复杂,太多的波段也会导致发射率方程组的病态程度加深,不适合实际应用。
此外,红外热成像技术通常是运用光电技术检测物体辐射的特定波段的红外线光辐射,并在检测到的红外线光辐射能量与物体的表面温度之间建立相应关系,进而获得物体的红外热像图。这种红外热像图与物体表面的热分布场相对应,热图像上的不同颜色代表被测物体的不同区域具有不同的温度。
传统的红外热成像技术通常需要面阵探测器,其成像灵敏度受限于面阵探测器对于单位像素的探测灵敏度,且面阵探测带来了测量维度和测量数上的冗余,限制了其在暗场环境下的应用。
发明内容
本发明目的在于克服现有技术中的测温应用场合窄、单波长测温辐射率难以修正、误差大、多波长测温结构复杂、传统双波长测温灵敏度低、测量维度高、扩展性差等缺陷。
根据本发明的一个方面,提供了一种基于编码变换的双波长温度场成像设备,包括:
光辐射调制装置,配置为接收待测对象的光辐射,且加载预设的多个掩膜,将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿不同于第一路径的第二路径射出,所述多个掩膜由矩阵Φ变换生成;
布置在所述第一路径上的第一滤光元件,配置为接收多束所述第一光辐射,并将接收到的所述第一光辐射过滤为波长为第一波长λ1的多束光;
布置在所述第二路径上的第二滤光元件,配置为接收多束所述第二光辐射,并将接收到的所述第二光辐射过滤为波长为第二波长λ2的多束光;
布置在第一路径上的第一探测装置,配置为接收所述波长为第一波长λ1的多束光并将其转换为相应的多个第一光电信号参量;
布置在第二路径上的第二探测装置,配置为接收所述波长为第二波长λ2的多束光并将其转换为相应的多个第二光电信号参量;
温度确定装置,配置为接收来自所述第一探测装置和所述第二探测装置的多个所述第一和第二光电信号参量,并根据多个所述第一和第二光电信号参量与温度的预定关系确定出所述待测对象每一个像素点的温度值;
图像生成装置,配置为根据所述待测对象每一个像素点的温度值以及所述待测对象的二维图像反演出所述待测对象的二维红外热图像。
本实施方式采用的光辐射调制装置能够加载根据预设矩阵Φ变换生成的多个掩膜,并且能够以加载每一个掩膜时测得的第一光电信号参量和第二光电信号参量为基础,获得待测物体的二维红外热图像。
在本发明的一些实施方式中,所述光辐射调制装置加载预设的多个掩膜,所述第一探测装置接收所述波长为第一波长λ1的多束光并将其转换为相应的多个第一光电信号参量,所述第二探测装置接收所述波长为第二波长λ2的多束光并将其转换为相应的多个第二光电信号参量包括:
当所述预设矩阵Φ矩阵服从±1二值分布时:
将预设矩阵Φ拆分为两个互补的0-1矩阵H+和H-
所述光辐射调制装置加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射调制为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
所述光辐射调制装置加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从±1、0三值分布时:
将所述预设矩阵Φ拆分为两个相互独立的0-1矩阵H+和H-
所述光辐射调制装置加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
所述光辐射调制装置加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从0-1分布时:
所述光辐射调制装置顺序加载由预设矩阵Φ的每一行(或列)直接拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)i,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)i
其中,i=1,2,3,4......N,N为所述待测对象的总像素数,所述预设矩阵Φ的阶2k≥N。
在本发明的一些实施方式中,所述温度确定装置根据如下的所述预定关系确定所述待测对象每一个像素点的温度值:
当所述预设矩阵Φ服从±1二值分布时:
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
所述温度确定装置根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
当所述预设矩阵Φ服从±1、0三值分布时:
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
所述温度确定装置根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
当所述预设矩阵Φ服从0-1分布时:
所述温度确定装置根据Yi=E1(T)i/E2(T)i、结合数学模型Yi=ΦS1,并利用矩阵求逆方法计算出待测物体的二维图像S;
所述温度确定装置根据待测对象的二维图像S以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
其中,A’、B’为预设系数。
在本发明的一些实施方式中,所述光辐射调制装置包括:空间光调制器,配置为根据预设矩阵Φ变换生成的多个掩膜,以将接收到的待测对象的光辐射调制为所述第一光辐射和第二光辐射,并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出;控制元件,配置为控制所述空间光调制器依次加载由所述预设矩阵Φ变换生成的多个掩膜。
在本发明的一些实施方式中,所述空间光调制器为选自数字微镜器件、光强数字调制器或液晶光阀。
在本发明的一些实施方式中,所述第一探测装置为第一点探测器,所述第二探测装置为第二点探测器,并且
所述基于编码变换的双波长温度场成像设备还包括布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一会聚元件,和布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二会聚元件,
所述第一点探测器位于所述第一会聚元件的光焦点处;
所述第二点探测器位于所述第二会聚元件的光焦点处。
本实施方式采用两个点探测器便完成了原本需要两个面阵测温器件才能完成的工作,极大地降低了测量维度,大幅增加光通量,避免了红外光通量在维度上的分配,而将噪声压制在单像素水平,信噪比大幅提高。
在本发明的一些实施方式中,所述的基于编码变换的双波长温度场成像设备还包括:
布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一光强衰减元件,和
布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二光强衰减元件。
在本发明的一些实施方式中,所述第一滤光元件和所述第二滤光元件为中心波长相差至少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次定标时测得的波长为第二波长的光的光电信号参量。
在本发明的一些实施方式中,所述定标装置还包括将基准光源的光辐射转化为平行光辐射的扩束准直透镜。
在本发明的一些实施方式中,所述定标装置还包括将所述扩束准直透镜转化的平行光辐射射向所述光辐射调制装置的分束器。
另外,本发明还提供了一种基于编码变换的双波长温度场成像方法,包括:
利用光辐射调制装置接收待测对象的光辐射,并且加载预设的多个掩膜,而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束第一光辐射沿第一路径射出、多束第二光辐射沿第二路径射出,所述多个掩膜由矩阵Φ变换生成;
在上述的左臂方向的路径上接收多束第一光辐射,并将其过滤为波长为第一波长λ1的多束光;
在上述的右臂方向的路径上接收多束第二光辐射,并将其过滤为波长为第二波长λ2的多束光;
在上述的左臂方向的路径上接收多束波长为λ1的单波光并将其转换为多个相应的第一光电信号参量;
在上述的右臂方向的路径上接收多束波长为λ2的单波光并将其转换为多个相应的第二光电信号参量;
根据多个所述第一和第二光电信号参量与温度的预定关系确定出所述待测对象的每一个像素点的温度Ti;以及
根据所述待测对象的二维图像以及每一个像素点的温度值Ti反演出所述待测对象的二维红外热图像。
在本发明的一些实施方式中,所述光辐射调制装置加载所述多个掩膜,所述第一探测装置和所述第一探测装置根据所述掩膜探测多个所述第一光电信号参量和所述多个所述第二光电信号参量包括:
当所述预设矩阵Φ矩阵服从±1二值分布时:
将预设矩阵Φ拆分为两个互补的0-1矩阵H+和H-
加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且将接收到的光辐射调制为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从±1、0三值分布时:
将所述预设矩阵Φ拆分为两个相互独立的0-1矩阵H+和H-
加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从0-1分布时:
顺序加载由预设矩阵Φ的每一行(或列)直接拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)i,将所述第二光辐射转换为相应的第二光电信号参量E2(T)i;其中,i=1,2,3,4......N,N为所述待测对象的总像素数,所述预设矩阵Φ的阶2k≥N。
在本发明的一些实施方式中,根据多个所述第一和第二光电信号参量与温度的预定关系确定出所述待测对象的每一个像素点的温度Ti包括:
当所述预设矩阵Φ服从±1二值分布时:
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
当所述预设矩阵Φ服从±1、0三值分布时:
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
当所述预设矩阵Φ服从0-1分布时:
根据Yi=E1(T)i/E2(T)i、结合数学模型Yi=ΦS1,并利用矩阵求逆方法计算出待测物体的二维图像S;
根据待测对象的二维图像S以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
其中,A’、B’为预设系数。
在本发明的一些实施方式中,所述光辐射调制装置包括:空间光调制器,配置为根据预设矩阵Φ变换生成的多个掩膜,以将接收到的待测对象的光辐射调制为所述第一光辐射和第二光辐射,并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出;控制元件,配置为控制所述空间光调制器依次加载由所述预设矩阵Φ变换生成的多个掩膜。
在本发明的一些实施方式中,所述空间光调制器选自数字微镜器件、光强数字调制器或液晶光阀。
在本发明的一些实施方式中,所述的基于编码变换的双波长温度场成像方法进一步包括:
将所述波长为第一波长的光会聚于第一焦点,并在所述第一焦点处设置第一点探测装置用以接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量;
将所述波长为第二波长的光会聚于第二焦点,并在所述第二焦点处设置第二点探测装置用以接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量。
在本发明的一些实施方式中,所述的基于编码变换的双波长温度场成像方法进一步包括:
对所述波长为第一波长的光的强度进行衰减;及
对所述波长为第二波长的光的强度进行衰减。
在本发明的一些实施方式中,所述光辐射为红外波段的光辐射。
在本发明的一些实施方式中,所述第一波长与所述第二波长相差至少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次定标的波长为第二波长的光的光电信号参量。
本发明提供的基于编码变换的双波长温度场成像设备、系统及其方法将双波长测温技术、调制技术和编码技术相结合,适用于双波长测温领域。发明提供的基于编码变换的双波长温度场成像设备、系统及方法利用矩阵求逆方法或矩阵乘法反演出二维图像,再结合已知的拟合系数,确定二维图像上每一个像素点所对应的温度,最后重建出待测对象的二维红外热图像。凭借着显著的优势,该设备、系统及方法将替代原有的成像技术,广泛地应用在深空探测、遥感、材料检测、夜视观测等相关科技领域。
附图说明
图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为本发明一些实施方式的基于编码变换的双波长温度场成像方法的定标步骤流程示意图。
具体实施方式
为了使发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施例对发明作进一步详细的说明。虽然附图中显示了本公开示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻的理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。
参照图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无限接近时,ε11,T)≈ε22,T),则:
A′=A11)/A22)=(ε11,T)C1λ1 -5)/(ε22,T)C1λ2 -5)≈(λ12)-5
B′=B11)-B22)=-C21-(-C22),
但是,λ1与λ2越接近,对探测器的灵敏度和准确性的要求越高,另外考虑到环境噪声和探测器本征噪声(如暗计数等)对系统的测温准确性造成的较大影响,λ1与λ2在实际测量中不可能无限接近,因而两个辐射波长下的辐射率ε11,T)和ε22,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 Virtex5 application 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)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)=A11)exp(B11)/T)
Ed2(T)=A22)exp(B22)/T)
A’=A11)/A22),B’=B11)-B22)。
具体可以是利用公式拟合得出相应系数A’和B’,得到此种环境下的该辐射体的光辐射温度曲线公式;然后通过计算机程序设定好此系数下的计算公式,即可测量或者监测任何情况下的辐射温度。在实际测温过程中,利用指数形式的多项式展开,获得多项式的系数,或者利用其它复杂的拟合方程获得拟合系数。
由此,本实施方式有效地克服了红外测温中各种物体的“辐射率修正”难题,克服了测量条件复杂、现场测量条件波动或者水蒸气等因素的环境吸收造成的测量误差,实现了高精度的温度测量。
在本实施方式中,在实际测温过程中,可以通过计算机程序设定好此系数下的计算公式,即可测量或者监测任何情况下的辐射温度。可以利用指数形式的多项式展开,获得多项式的系数,或者利用其它复杂的拟合方程获得拟合系数,可在实际操作过程中实现更高精度的测量。这里实际是利用现场定标的实验方法找出实际条件下待测对象(某一具体物体,灰体而非黑体)的热辐射规律,并给与校准,再利用该规律在原有条件下进行温度监测或测温,即利用“替代法”校准该实际条件下的许多隐含参量。
本实施方式中,温度确定装置(如除法器、计算元件)可以由单片机、FPGA芯片、计算机、服务器等替换。数据运算可以通过程序模块执行。程序可以包括例程、程序、目标程序、组件、逻辑、数据结构等,它们执行特定的任务或者实现特定的抽象数据类型。计算机系统/服务器可以在分布式云计算环境中实施,分布式云计算环境中,任务是由通过通信网络链接的远程处理设备执行的。在分布式云计算环境中,程序模块可以位于包括存储设备的本地或远程计算系统存储介质上。
如图1~7所示的本发明所提供的基于光辐射的测温设备、系统和方法只能够测量待测对象(例如人体或者电灯)整体的平均温度,不能测量其每一个局部的温度并生成待测对象的二维红外热像图。
鉴于此,参照图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’接收待测对象的光辐射,并且加载根据预先设定的矩阵变换生成的多个掩膜。随着光辐射调制装置2’加载的掩膜发生变化,其将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、多束第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出。布置在第一路径上的第一滤光元件3-1接收多束所述第一光辐射,并将接收到的多束第一光辐射过滤为波长为第一波长的多束光(例如波长为λ1的单波光)。布置在第二路径上的第二滤光元件4-1接收多束第二光辐射,并将接收到的多束第二光辐射过滤为波长为第二波长的多束光(例如波长为λ2的单波光,其中λ1与λ2不等,当λ1与λ2为相邻波段的波长时,效果最佳,因为当λ1与λ2无限接近时,ε11,T)≈ε22,T),则:
A′=A11)/A22)=(ε11,T)C1λ1 -5)/(ε22,T)C1λ2 -5)≈(λ12)-5
B′=B11)-B22)=-C21-(-C22),
但是,λ1与λ2越接近,对探测器的灵敏度和准确性的要求越高,另外考虑到环境噪声和探测器本征噪声(如暗计数等)对系统的测温准确性造成的较大影响,λ1与λ2在实际测量中不可能无限接近,因而两个辐射波长下的辐射率ε11,T)和ε22,T)不可简单相消,而需要通过定标确定出预定关系。布置在第一路径上的第一探测装置3-4接收多束波长为第一波长的光并将其转换为相应的多个第一光电信号参量。布置在第二路径上的第二探测装置4-4接收多束波长为第二波长的光并将其转换为相应的多个第二光电信号参量。
温度确定装置5接收来自第一探测装置3-4和第二探测装置4-4的多个第一和多个第二光电信号参量,并根据多个所述第一和第二光电信号参量确定待测对象在波长λ1和波长λ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根据预先设置的矩阵变换生成的多个掩膜,以将接收到的待测对象的光辐射调制为所述第一光辐射和第二光辐射,并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出。控制元件2-2配置为控制所述空间光调制器依次加载由所述预设矩阵Φ变换生成的多个掩膜。例如:预设矩阵Φ,根据预设矩阵Φ变换生成掩膜,最后根据掩膜中0和1的数量控制空间光调制器2-1中微镜的翻转。
其中,所采用的预设矩阵Φ可以服从±1二值分布(诸如Hadamard矩阵),或服从±1、0三值分布,0、1分布(服从其它分布的0-1矩阵)。该部分内容下文还会继续介绍。
在本实施例中,空间光调制器可以选用如图3(a)和图3(b)所示的数字微镜器件(Digital MicromirrorDevice,DMD)、光强数字调制器或液晶光阀中的任意一种。
当采用不同的预设矩阵Φ时,光辐射调制装置2’对于预设矩阵Φ的处理方式、加载掩膜的方式均不同,温度确定装置5处理光电信号参量的方式也不同。
下面分别描述采用不同的预设矩阵Φ时,光辐射调制装置2’、温度确定装置5和图像生成装置8的工作原理:
在此之前,首先介绍待测对象的像素概念,待测对象的二维图像可以被分割成如下所示的2×4的矩阵,即待测对象的总像素N=p×q(p为横坐标像素数、q为纵坐标像素数)。
因此,实际上是以一个包含的元素个数与待测对象的总像素数相同的矩阵表示待测对象的二维图像。
其次,光辐射调制装置2’所加载的矩阵的阶2k≥N,即光辐射调制装置2’多加载的掩膜大于或等于待测对象的总像素数。
①当采用的预设矩阵Φ服从±1二值分布时,例如图9所示的8×8的Hadamard矩阵:
其中,Hadamard矩阵满足以下特点:
Hadamard矩阵Hk中的元素只有1和-1(本发明不考虑系数),其中1和-1比例为1∶1。
控制元件2-2将预设矩阵Φ拆分为两个互补的0-1矩阵H+和H-,即H-=1-H+,Φ=(H+)-(H-),则:
空间光调制器2-1交替加载互补掩膜直到第N次结束。
温度确定装置5根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1,其中:
根据ΦS1=E1(T)2i-1-E1(T)2i,可以得到下列方程组:
即X11)+X21)+X31)+X41)+X51)+X61)+X71)+X81)=E1(T)1-E1(T)2 (1);
即X11)-X21)+X31)-X41)+X51)-X61)+X71)-X81)=E1(T)3-E1(T)4 (2);
即X11)+X21)-X31)-X41)+X51)+X61)-X71)-X81)=E1(T)5-E1(T)6 (3);
即X11)-X21)-X31)+X41)+X51)-X61)-X71)+X81)=E1(T)7-E1(T)8 (4);
即X11)+X21)+X31)+X41)-X51)-X61)-X71)-X81)=E1(T)9-E1(T)10 (5);
即X11)-X21)+X31)-X41)-X51)+X61)-X71)+X81)=E1(T)11-E1(T)12 (6);
即X11)+X21)-X31)-X41)-X51)-X61)+X71)+X8(λ1)=E1(T)13-E1(T)14 (7);
即X11)-X21)-X31)+X41)-X51)+X61)+X71)-X81)=E1(T)15-E1(T)16 (8);
联立方程(1)~(8),能够解得X11)~X81),获得在波长λ1下待测对象的二维图像S1
同样地,温度确定装置5根据结合数学模型并利用矩阵求逆方法解得X12)~X82),进而获得待测对象在波长λ2下待测对象的二维图像S2
温度确定装置5根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
温度确定装置5对待测对象在波长λ1和波长λ2下的二维图像S1和S2进行定点除法运算,获得:
根据T1=B’/ln(X1/A’),计算该像素坐标下的温度T1
根据T2=B’/ln(X2/A’),计算该像素坐标下的温度T2
……
根据T8=B’/ln(X8/A’),计算该像素坐标下的温度T8
②当采用的预设矩阵Φ服从±1、0三值分布(矩阵Φ的阶2k≥N)时:
控制元件2-2将矩阵Φ拆分为两个相互独立的0-1矩阵H+和H-,即H-=1-H+,Φ=(H+)-(H-)。
即,如果预设矩阵Φ中存在0,则其在矩阵H+和矩阵H-中相应元素为均为0或者均为1。
在此之后的掩膜加载、待测对象每一个像素点的温度确定以及待测对象的二维红外热图像的生成步骤与①中的步骤完全相同,在此不再赘述。
③当采用的预设矩阵Φ为0-1分布(服从其它分布的0-1矩阵)时:
例如,所采用的预设矩阵Φ:
空间光调制器2-1顺序加载由矩阵Φ的每一行直接拉伸而得的掩膜,并且将接收到的光辐射调制为第一光辐射和第二光辐射,第一探测装置3-4将第一光辐射转换为相应的第一光电信号参量E1(T)i,第二探测装置4-4将第二光辐射转换为相应的第二光电信号参量E2(T)i,i=1,2,3,4...N;
具体而言,由Yi=E1(T)i/E2(T)i、Yi=ΦS1得ΦS=E1(T)i/E2(T)i
根据ΦS=E1(T)i/E2(T)i,可以得到下列方程组:
即X1+X2+X3+X4+X5+X6+X7+X8=E1(T)1/E2(T)1 (1);
即X1+X3+X4+X7=E1(T)2/E2(T)2 (2);
即X1+X2+X5+X6=E1(T)3/E2(T)3 (3);
即X1+X4+X5+X8=E1(T)4/E2(T)4 (4);
即X1+X2+X3+X4=E1(T)5/E2(T)5 (5);
即X1+X3+X6+X8=E1(T)6/E2(T)6 (6);
即X1+X2+X7+X8=E1(T)7/E2(T)7 (7);
即X1+X4+X6+X7=E1(T)8/E2(T)8 (8)。
联立方程(1)~(8),能够解得X1~X8,获得待测对象的二维图像:
应当能够理解,顺序加载由预设矩阵Φ的每一列直接拉伸而得的掩膜也是可行的。
温度确定装置5根据Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti
在本发明的实施方式中,温度确定装置5可以采用定义法求逆、伴随矩阵求逆法、初等变换法求逆、广义的行列初等变换法求逆、Gauss-Jordan法求逆、和化积法求逆、多项式法求逆、矩阵函数的级数展开法求逆、分解矩阵求逆、设元法求逆、递推法求逆、分块矩阵求逆、解方程组求逆、克莱姆法则求逆、行列式求逆、恒等变形法求逆、Hamilton-Caley定理求逆、三角矩阵求逆、拼接新矩阵求逆等中的任意一种以及Y=ΦS的数学模型,确定出待测对象的二维图像。
图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,下面介绍一优选实施方式的基于编码变换的双波长温度场成像设备的工作方式,其具体可以是:
待测对象(例如人体或者电灯,图中未进行标示)的光辐射可以经过透镜1射向光辐射调制装置2’。光辐射调制装置2’加载根据预先设定的矩阵变换生成的多个掩膜,且将接收的光辐射调制为双臂方向上的光辐射。在双臂中的左臂方向上分别布置有第一滤光元件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
温度确定装置5根据待测对象在波长λ1和波长λ2下的二维图像S1和S2以及所述待测对象的温度的预定关系确定出待测对象每一个像素点的温度。
图像生成装置8用于根据温度确定装置5生成的待测对象的每一个像素点的温度值以及待测对象的二维图像反演出待测对象的二维红外热图像。
在一些实施例中,为了使过滤后的单波的波长不同,第一滤光元件和所述第二滤光元件为中心波长相差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实施方式中的光辐射调制装置)接收待测对象(例如人体或者电灯)的光辐射(例如红外辐射、紫外辐射或者可见光),并且加载预设的多个掩膜,而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束第一光辐射沿第一路径(如测温设备的左臂方向的路径)射出、多束第二光辐射沿第二路径(如测温设备的右臂方向的路径)射出,所述多个掩膜由矩阵Φ变换生成;
S1202:在上述的左臂方向的路径上接收多束第一光辐射,并将其过滤为波长为第一波长的多束光(例如波长为λ1的单波光);
S1203:在上述的右臂方向的路径上接收多束第二光辐射,并将其过滤为波长为第二波长的多束光(例如波长为λ2的单波光);
S1204:在上述的左臂方向的路径上接收多束波长为λ1的单波光并将其转换为多个相应的第一光电信号参量E1(T)i
S1205:在上述的右臂方向的路径上接收多束波长为λ2的单波光并将其转换为多个相应的第二光电信号参量E2(T)i
S1206:根据多个所述第一和第二光电信号参量E1(T)i和E2(T)i(光辐射调制装置加载每一个掩膜时所获得的第一和第二光电信号参量)与温度的预定关系确定出待测对象(例如人体或者电灯)的每一个像素点的温度Ti
S1207:根据待测对象的二维图像以及每一个像素点的温度值Ti反演出待测对象的二维红外热图像。
在本实施方式中,所述光辐射调制装置根据预设矩阵Φ生成和加载掩膜的方式如下:
①当采用的预设矩阵Φ服从±1二值分布时:
a.提供2k阶的预设矩阵Φ,其中2k≥N=p×q(p为待测对象横坐标像素数、q为待测对象纵坐标像素数);
b.将预设矩阵Φ拆分为两个互补的0-1矩阵H+和H-
c.加载由H+矩阵的第i行(或第i列)H+i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,将第二光辐射转换为相应的第二光电信号参量E2(T)2i-1,i=1,2,3,4......N;
d.加载由H-矩阵的第i行(或第i列)H-i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将第一光辐射转换为相应的第一光电信号参量E1(T)2i,将第二光辐射转换为相应的第二光电信号参量E2(T)2i,i=1,2,3,4......N。
②当采用的预设矩阵Φ服从±1、0三值分布时:
a.提供2k阶的预设矩阵Φ,其中2k≥N=p×q(p为待测对象横坐标像素数、q为待测对象纵坐标像素数);
b.将预设矩阵Φ拆分为两个相互独立的0-1矩阵H+和H-
c.加载由H+矩阵的第i行(或第i列)H+i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,将第二光辐射转换为相应的第二光电信号参量E2(T)2i-1,i=1,2,3,4......N;
d.加载由H-矩阵的第i行(或第i列)H-i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将第一光辐射转换为相应的第一光电信号参量E1(T)2i,将第二光辐射转换为相应的第二光电信号参量E2(T)2i,i=1,2,3,4......N。
③当采用的预设矩阵Φ服从0-1分布(其它分布的0-1矩阵)时:
a.提供2k阶的预设矩阵Φ,其中2k≥N=p×q(p为待测对象横坐标像素数、q为待测对象纵坐标像素数);
b.顺序加载由预设矩阵Φ的每一行(或列)直接拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将第一光辐射转换为相应的第一光电信号参量E1(T)i,将第二光辐射转换为相应的第二光电信号参量E2(T)i,i=1,2,3,4......N;
特别地,当预设矩阵Φ为单位阵时,即空间光调制元件2-1每一帧的掩膜中只有一个像素,且依次以逐点扫描方式出现,则整个采样为逐点扫描采样,但将损失高通量的特性。
在本实施方式中,所述温度确定装置根据如下的预定关系确定待测对象的每个像素点的温度值:
①当采用的预设矩阵Φ服从±1二值分布时:
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti,其中,A’、B’为预设系数。
②当采用的预设矩阵Φ服从±1、0三值分布:
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
根据待测对象在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti,其中,A’、B’为预设系数。
③当采用的预设矩阵Φ为0-1分布(服从其它分布的0-1矩阵)时:
根据Yi=E1(T)i/E2(T)i、结合数学模型Yi=ΦS1,并利用矩阵求逆方法计算出待测物体的二维图像S;
根据待测对象的二维图像S以及所述Ti=B’/ln(Xi/A’)确定出待测对象每个像素点的温度值Ti,其中,A’、B’为预设系数。
二维红外热像图与物体表面的热分布相对应,也即反映被测对象的光辐射(红外热辐射)能量分布,二维红外热像图的上面的不同颜色代表被测物体不同区域的不同温度,通过二维红外热像图可以观察到被测物体的整体温度分布状况,研究被测物体的发热情况。
在一些实施方式中,该方法还包括:
将所述波长为第一波长的光会聚于第一焦点,并在所述第一焦点处设置第一点探测装置用以接收所述波长为第一波长的多束光并将其转换为多个相应的第一光电信号参量;
将所述波长为第二波长的光会聚于第二焦点,并在所述第二焦点处设置第二点探测装置用以接收所述波长为第二波长的多束光并将其转换为多个相应的第二光电信号参量。
在一些实施方式中,该方法还包括:
对所述波长为第一波长的多束光的强度进行衰减;及对所述波长为第二波长的多束光的强度进行衰减。
在本实施方式中,所述光辐射为红外线光辐射。
在本实施方式中,所述第一波长与所述第二波长相差至少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)=A11)exp(B11)/T)
Ed2(T)=A22)exp(B22)/T)
A’=A11)/A22),B’=B11)-B22)。
具体可以是利用公式拟合得出相应系数A’和B’,得到此种环境下的该辐射体的光辐射温度曲线公式;然后通过计算机程序设定好此系数下的计算公式,即可测量或者监测任何情况下的辐射温度。在实际测温过程中,利用指数形式的多项式展开,获得多项式的系数,或者利用其它复杂的拟合方程获得拟合系数。
由此,本实施方式有效地克服了红外测温中各种物体的“辐射率修正”难题,克服了测量条件复杂、现场测量条件波动或者水蒸气等因素的环境吸收造成的测量误差,实现了高精度的温度测量。
最后应说明的是,以上实施例仅用以说明本发明的技术方案而非限制性的。尽管参照实施例对本发明进行了详细说明,但本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (27)

1.一种基于编码变换的双波长温度场成像设备,其特征在于,包括:
光辐射调制装置,配置为接收待测物体的光辐射,且加载预设的多个掩膜,将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束所述第一光辐射沿第一路径射出、多束所述第二光辐射沿不同于第一路径的第二路径射出,所述多个掩膜由预设矩阵Φ变换生成;
布置在所述第一路径上的第一滤光元件,配置为接收多束所述第一光辐射,并将接收到的所述第一光辐射过滤为波长为第一波长λ1的多束光;
布置在所述第二路径上的第二滤光元件,配置为接收多束所述第二光辐射,并将接收到的所述第二光辐射过滤为波长为第二波长λ2的多束光;
布置在第一路径上的第一探测装置,配置为接收所述波长为第一波长λ1的多束光并将其转换为相应的多个第一光电信号参量;
布置在第二路径上的第二探测装置,配置为接收所述波长为第二波长λ2的多束光并将其转换为相应的多个第二光电信号参量;
温度确定装置,配置为接收来自所述第一探测装置和所述第二探测装置的多个所述第一和第二光电信号参量,并根据多个所述第一和第二光电信号参量与温度的预定关系确定出所述待测物体每一个像素点的温度值;
图像生成装置,配置为根据所述待测物体每一个像素点的温度值以及所述待测物体的二维图像反演出所述待测物体的二维红外热图像。
2.根据权利要求1所述的基于编码变换的双波长温度场成像设备,其特征在于,所述光辐射调制装置加载预设的多个掩膜,所述第一探测装置接收所述波长为第一波长λ1的多束光并将其转换为相应的多个第一光电信号参量,所述第二探测装置接收所述波长为第二波长λ2的多束光并将其转换为相应的多个第二光电信号参量包括:
当所述预设矩阵Φ矩阵服从±1二值分布时:
将预设矩阵Φ拆分为两个互补的0-1矩阵H+和H-
所述光辐射调制装置加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射调制为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
所述光辐射调制装置加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从±1、0三值分布时:
将所述预设矩阵Φ拆分为两个相互独立的0-1矩阵H+和H-
所述光辐射调制装置加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
所述光辐射调制装置加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从0-1分布时:
所述光辐射调制装置顺序加载由预设矩阵Φ的每一行或列直接拉伸变换而得的掩膜,并且所述光辐射调制装置将接收到的光辐射分为第一光辐射和第二光辐射,所述第一探测装置将所述第一光辐射转换为相应的第一光电信号参量E1(T)i,所述第二探测装置将所述第二光辐射转换为相应的第二光电信号参量E2(T)i
其中,i=1,2,3,4……N,N为所述待测物体的总像素数,所述预设矩阵Φ的阶2k≥N。
3.根据权利要求2所述的基于编码变换的双波长温度场成像设备,其特征在于,所述温度确定装置根据如下的所述预定关系确定所述待测物体每一个像素点的温度值:
当所述预设矩阵Φ服从±1二值分布时:
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
所述温度确定装置根据待测物体在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及Ti=B’/ln(Xi/A’)确定出待测物体每个像素点的温度值Ti
当所述预设矩阵Φ服从±1、0三值分布时:
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
所述温度确定装置根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
所述温度确定装置根据待测物体在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及Ti=B’/ln(Xi/A’)确定出待测物体每个像素点的温度值Ti
当所述预设矩阵Φ服从0-1分布时:
所述温度确定装置根据Yi=E1(T)i/E2(T)i、结合数学模型Yi=ΦS1,并利用矩阵求逆方法计算出待测物体的二维图像S;
所述温度确定装置根据待测物体的二维图像S以及Ti=B’/ln(Xi/A’)确定出待测物体每个像素点的温度值Ti
其中,A’、B’为预设系数。
4.根据权利要求1~3中任一项所述的基于编码变换的双波长温度场成像设备,其特征在于,所述光辐射调制装置包括:
空间光调制器,配置为根据预设矩阵Φ变换生成的多个掩膜,以将接收到的待测物体的光辐射调制为所述第一光辐射和第二光辐射,并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出;
控制元件,配置为控制所述空间光调制器依次加载由所述预设矩阵Φ变换生成的多个掩膜。
5.根据权利要求4所述的基于编码变换的双波长温度场成像设备,其特征在于,所述空间光调制器为选自数字微镜器件、光强数字调制器或液晶光阀。
6.根据权利要求4所述的基于编码变换的双波长温度场成像设备,其特征在于,所述第一探测装置为第一点探测器,所述第二探测装置为第二点探测器,并且
所述基于编码变换的双波长温度场成像设备还包括布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一会聚元件,和布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二会聚元件,
所述第一点探测器位于所述第一会聚元件的光焦点处;
所述第二点探测器位于所述第二会聚元件的光焦点处。
7.根据权利要求6所述的基于编码变换的双波长温度场成像设备,其特征在于,还包括:
布置在所述第一路径上、位于所述第一点探测器与所述空间光调制器之间的第一光强衰减元件,和
布置在所述第二路径上、位于所述第二点探测器与空间光调制器之间的第二光强衰减元件。
8.根据权利要求1~3中任一项所述的基于编码变换的双波长温度场成像设备,其特征在于,所述第一滤光元件和所述第二滤光元件为中心波长相差至少10nm的第一窄带滤光片和第二窄带滤光片,所述第一窄带滤光片和第二窄带滤光片的半高宽参数至少为10nm。
9.根据权利要求1~3中任一项所述的基于编码变换的双波长温度场成像设备,其特征在于,所述光电信号参量包括光子数、电流值、电压值、电阻值中任意一种。
10.根据权利要求1~3中任一项所述的基于编码变换的双波长温度场成像设备,其特征在于,所述光辐射为红外波段的光辐射。
11.根据权利要求10所述的基于编码变换的双波长温度场成像设备,其特征在于,
所述第一和第二点探测器选自近红外、中远红外、远红外波段的外光电效应探测器组、内光电效应探测器组、强光探测器组和弱光探测器组中的任意一种,其中,
所述外光电效应探测器组包括:雪崩二极管、真空光电管、充气光电管、光电倍增管、变像管、像增强器、摄像管;
所述内光电效应探测器组包括:本征型光电导探测器、掺杂型光电导探测器、光磁电效应探测器、光生伏特探测器;
所述强光探测器组包括:内置或者外置有模数转换器的强光探测器;
所述弱光探测器组包括:内置或者外置有计数器的弱光探测器。
12.一种基于编码变换的双波长温度场成像系统,其特征在于,包括:
如权利要求1~11中任一项所述基于编码变换的双波长温度场成像设备以及定标装置,
所述定标装置包括温度可调节的基准光源,所述基准光源配置为在一个定标阶段中,将调节成的不同温度的光辐射射向所述光辐射调制装置,由所述光辐射调制装置将接收到的光辐射等分为第一光辐射和第二光辐射,使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出,以获取多个第一和第二光电信号参量,并根据所述不同温度与所述第一和第二光电信号参量的计量值,确定出所述第一和第二光电信号参量与温度的预定关系。
13.根据权利要求12所述的基于编码变换的双波长温度场成像系统,其特征在于,所述第一和第二光电信号参量与温度的预定关系为:
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次定标时测得的波长为第二波长的光的光电信号参量。
14.根据权利要求12或13所述的基于编码变换的双波长温度场成像系统,其特征在于,所述定标装置还包括将基准光源的光辐射转化为平行光辐射的扩束准直透镜。
15.根据权利要求14所述的基于编码变换的双波长温度场成像系统,其特征在于,所述定标装置还包括将所述扩束准直透镜转化的平行光辐射射向所述光辐射调制装置的分束器。
16.一种基于编码变换的双波长温度场成像方法,其特征在于,包括:
利用光辐射调制装置接收待测物体的光辐射,并且加载预设的多个掩膜,而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束第一光辐射沿第一路径射出、多束第二光辐射沿第二路径射出,所述多个掩膜由预设矩阵Φ变换生成;
在上述的第一路径上接收多束第一光辐射,并将其过滤为波长为第一波长λ1的多束光;
在上述的第二路径上接收多束第二光辐射,并将其过滤为波长为第二波长λ2的多束光;
在上述的第一路径上接收多束波长为λ1的单波光并将其转换为多个相应的第一光电信号参量;
在上述的第二路径上接收多束波长为λ2的单波光并将其转换为多个相应的第二光电信号参量;
根据多个所述第一和第二光电信号参量与温度的预定关系确定出所述待测物体的每一个像素点的温度Ti;以及
根据所述待测物体的二维图像以及每一个像素点的温度值Ti反演出所述待测物体的二维红外热图像。
17.根据权利要求16所述的基于编码变换的双波长温度场成像方法,其特征在于,所述光辐射调制装置加载所述多个掩膜,第一探测装置和第二探测装置根据所述掩膜探测多个所述第一光电信号参量和多个所述第二光电信号参量包括:
当所述预设矩阵Φ服从±1二值分布时:
将预设矩阵Φ拆分为两个互补的0-1矩阵H+和H-
加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且将接收到的光辐射调制为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从±1、0三值分布时:
将所述预设矩阵Φ拆分为两个相互独立的0-1矩阵H+和H-
加载由H+矩阵的第i行或第i列H+i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i-1,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i-1
加载由H-矩阵的第i行或第i列H-i拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)2i,将所述第二光辐射转换为相应的第二光电信号参量E2(T)2i
当所述预设矩阵Φ服从0-1分布时:
顺序加载由预设矩阵Φ的每一行或列直接拉伸变换而得的掩膜,并且将接收到的光辐射分为第一光辐射和第二光辐射,将所述第一光辐射转换为相应的第一光电信号参量E1(T)i,将所述第二光辐射转换为相应的第二光电信号参量E2(T)i
其中,i=1,2,3,4……N,N为所述待测物体的总像素数,所述预设矩阵Φ的阶2k≥N。
18.根据权利要求17所述的基于编码变换的双波长温度场成像方法,其特征在于,根据多个所述第一和第二光电信号参量与温度的预定关系确定出所述待测物体的每一个像素点的温度Ti包括:
当所述预设矩阵Φ服从±1二值分布时:
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
根据待测物体在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及Ti=B’/ln(Xi/A’)确定出待测物体每个像素点的温度值Ti
当所述预设矩阵Φ服从±1、0三值分布时:
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ1下的二维图像S1
根据结合数学模型并利用矩阵求逆方法计算出待测物体在波长λ2下的二维图像S2
根据待测物体在波长λ1和波长λ2下的二维图像S1和S2、Xi=S1./S2以及Ti=B’/ln(Xi/A’)确定出待测物体每个像素点的温度值Ti
当所述预设矩阵Φ服从0-1分布时:
根据Yi=E1(T)i/E2(T)i、结合数学模型Yi=ΦS1,并利用矩阵求逆方法计算出待测物体的二维图像S;
根据待测物体的二维图像S以及Ti=B’/ln(Xi/A’)确定出待测物体每个像素点的温度值Ti
其中,A’、B’为预设系数。
19.根据权利要求16~18中任一项所述的基于编码变换的双波长温度场成像方法,其特征在于,所述光辐射调制装置包括:
空间光调制器,配置为根据预设矩阵Φ变换生成的多个掩膜,以将接收到的待测物体的光辐射调制为所述第一光辐射和第二光辐射,并使第一光辐射沿第一路径射出、第二光辐射沿不同于第一路径的第二路径射出;
控制元件,配置为控制所述空间光调制器依次加载由所述预设矩阵Φ变换生成的多个掩膜。
20.根据权利要求19所述的基于编码变换的双波长温度场成像方法,其特征在于,所述空间光调制器选自数字微镜器件、光强数字调制器或液晶光阀。
21.根据权利要求16~18中任一项所述的基于编码变换的双波长温度场成像方法,其特征在于,进一步包括:
将所述波长为第一波长的光会聚于第一焦点,并在所述第一焦点处设置第一点探测装置用以接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量;
将所述波长为第二波长的光会聚于第二焦点,并在所述第二焦点处设置第二点探测装置用以接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量。
22.根据权利要求16~18中任一项所述的基于编码变换的双波长温度场成像方法,其特征在于进一步包括:
对所述波长为第一波长的光的强度进行衰减;及
对所述波长为第二波长的光的强度进行衰减。
23.根据权利要求16~18中任一项所述的基于编码变换的双波长温度场成像方法,其特征在于,所述光辐射为红外波段的光辐射。
24.根据权利要求16~18中任一项所述的基于编码变换的双波长温度场成像方法,其特征在于,所述第一波长与所述第二波长相差至少10nm。
25.根据权利要求16~18中任一项所述的基于编码变换的双波长温度场成像方法,其特征在于,所述光电信号参量包括光子数、电流值、电压值、电阻值中任意一种。
26.根据权利要求16~18中任一项所述的基于编码变换的双波长温度场成像方法,其特征在于,利用光辐射调制装置接收待测物体的光辐射,并且加载预设的多个掩膜,而后将接收到的光辐射调制为多束第一光辐射和多束第二光辐射,并使多束第一光辐射沿第一路径射出、多束第二光辐射沿第二路径射出,所述多个掩膜由预设矩阵Φ变换生成步骤之前还包括定标步骤,
所述定标步骤包括:
利用光辐射调制装置接收基准光源的光辐射,将接收到的光辐射等分为第一光辐射和第二光辐射,并使所述第一光辐射沿第一路径射出、所述第二光辐射沿不同于第一路径的第二路径射出;
在所述第一路径上接收所述第一光辐射,并将接收到的第一光辐射过滤为波长为第一波长的光;
在所述第二路径上接收所述第二光辐射,并将接收到的第二光辐射过滤为波长为第二波长的光;
在所述第一路径上接收所述波长为第一波长的光并将其转换为相应的第一光电信号参量;
在所述第二路径上接收所述波长为第二波长的光并将其转换为相应的第二光电信号参量;
调节基准光源发出多个不同温度的光辐射,并获取对应的多个第一和第二光电信号参量;以及
根据所述不同温度与所述第一和第二光电信号参量的计量值,确定出所述第一和第二光电信号参量与温度的预定关系。
27.根据权利要求26所述的基于编码变换的双波长温度场成像方法,其特征在于,所述第一和第二光电信号参量与温度的预定关系为:
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次定标的波长为第二波长的光的光电信号参量。
CN201610080112.1A 2016-02-04 2016-02-04 基于编码变换的双波长温度场成像设备、系统及方法 Active CN105606228B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610080112.1A CN105606228B (zh) 2016-02-04 2016-02-04 基于编码变换的双波长温度场成像设备、系统及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610080112.1A CN105606228B (zh) 2016-02-04 2016-02-04 基于编码变换的双波长温度场成像设备、系统及方法

Publications (2)

Publication Number Publication Date
CN105606228A CN105606228A (zh) 2016-05-25
CN105606228B true CN105606228B (zh) 2018-07-17

Family

ID=55986348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610080112.1A Active CN105606228B (zh) 2016-02-04 2016-02-04 基于编码变换的双波长温度场成像设备、系统及方法

Country Status (1)

Country Link
CN (1) CN105606228B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113049589B (zh) * 2021-03-10 2021-12-03 清华大学 高温环境下可调控的图像获取方法

Citations (6)

* Cited by examiner, † Cited by third party
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 中国科学院空间科学与应用研究中心 基于二阶压缩感知的超灵敏光谱成像天文望远镜及方法
CN205642623U (zh) * 2016-02-04 2016-10-12 北京理工大学 基于编码变换的双波长温度场成像设备及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9140795B2 (en) * 2011-09-20 2015-09-22 Mesa Imaging Ag Time of flight sensor with subframe compression and method

Patent Citations (6)

* Cited by examiner, † Cited by third party
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 中国科学院空间科学与应用研究中心 基于二阶压缩感知的超灵敏光谱成像天文望远镜及方法
CN205642623U (zh) * 2016-02-04 2016-10-12 北京理工大学 基于编码变换的双波长温度场成像设备及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
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
CN105606228A (zh) 2016-05-25

Similar Documents

Publication Publication Date Title
CN105675146B (zh) 基于压缩感知的双波长三维温度场成像设备、系统及方法
CN105737992B (zh) 基于压缩感知的双波长温度场成像设备、系统及方法
Fu et al. Measurement performance of an optical CCD-based pyrometer system
CN110530529A (zh) 红外热成像设备的检测系统
CN205642634U (zh) 基于压缩感知的双波长温度场成像设备及系统
CN205642633U (zh) 基于光辐射的测温设备及系统
JP2023535745A (ja) 光放射線を測定するための分光計装置及び方法
Yao et al. Non-uniformity calibration method of space-borne area CCD for directional polarimetric camera
CN105606228B (zh) 基于编码变换的双波长温度场成像设备、系统及方法
CN205642635U (zh) 基于压缩感知的双波长三维温度场成像设备及系统
CN205642623U (zh) 基于编码变换的双波长温度场成像设备及系统
CN113758573A (zh) 基于发射率迭代的三波长比色红外测温系统、方法及装置
JP4324693B2 (ja) 光検出器の分光応答度測定装置、その測定方法及び光源の分光放射照度校正方法
CN105527024B (zh) 基于光辐射的测温设备、系统及方法
Yang et al. Radiometric calibration algorithm for high dynamic range infrared imaging system
JP2004045306A (ja) 放射率分布測定方法および装置
Bonefačić et al. Two-color temperature measurement method using BPW34 PIN photodiodes
Coles et al. An automated system to measure the quantum efficiency of CCDs for astronomy
CN113418613A (zh) 一种基于多光谱比色的高温瞬态测量系统和方法
KR20220034887A (ko) 방사선 모니터링을 위한 방법 및 장치
Samedov et al. Filter-radiometer-based realization of candela and establishment of photometric scale at UME
Mahavarkar et al. Designing a sky-scanning Fabry–Perot interferometer system with a large size integrating sphere
Richwine et al. A comprehensive model for bolometer element and uncooled array design and imaging sensor performance prediction
Eppeldauer et al. 4. Transfer Standard Filter Radiometers: Applications to Fundamental Scales
Sametoglu New traceability chains in the photometric and radiometric measurements at the National Metrology Institute of Turkey

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