CN101487740A - 一种三ccd温度场测量装置及方法 - Google Patents

Abstract

本发明涉及一种三CCD温度场测量装置，包括：光学镜头，用于光学成像，将待测物体的辐射聚焦在分光棱镜组的入射面上；分光棱镜组，将投射的辐射等分为三束不同方向的波段辐射，分别从三个出射面出射；三个CCD面阵传感器，分别设置于分光棱镜组的三个出射面处，对三路经过不同光谱透过率分布的滤光片的出射辐射进行成像，获得三路光谱非相关的CCD图像；数据采集分析单元，对三路CCD图像数据进行采集，并利用多光谱测温法进行温度场计算。本发明还涉及一种对应的温度场测量方法。本发明的技术方案可以实现温度场测量，应用范围广泛；且技术方案实现简单，在高温检测等工业生产领域易于推广应用。

Description

一种三CCD温度场测量装置及方法

技术领域

涉及光学测温技术，尤其涉及一种基于三个CCD(ChargeCoupled Device，电荷耦合器件)面阵传感器光谱融合的温度场测量装置及方法。

背景技术

在石油化工、冶金、钢铁、水泥、玻璃等工业生产行业的高温检测领域，辐射测温仪器具有巨大的市场需求和广阔的应用空间。例如，冶金行业的高温炉膛内部温度测量与控制对于生产过程有着重要的作用。在这些典型的应用领域，传统的热电偶接触式测温手段，由于测量的局限性以及高成本的材料消耗，目前正在逐步被价格较低、性能稳定、低消耗使用、非接触式的光学测温设备所取代，光学测温设备的应用将成为高温测量的主流趋势。

随着图像光电传感器CCD不断地更新发展，基于CCD面阵传感器的光学测温方法与技术，在温度场测量方面展现出了极大的优势。许多科研机构应用CCD传感器开展了高温辐射温度场的测量研究工作，CCD传感器在辐射测温仪器研制领域展现了很好的应用前景，现有技术已有的应用与研究现状如下所述。

一种是采用单个彩色CCD相机(或摄像机)作为光学测温仪器，利用内嵌的RGB彩色滤色器阵列实现彩色复现，提供红、绿、蓝三个颜色通道：a、直接利用三个颜色通道的波段响应，结合特定的发射率模型，进行温度场的测量计算(中国科学G辑，34(6)：639-647，2004)。b、将三个通道的波段响应测量近似地处理为单色响应测量(中国电机工程学报，20(1)：70-72，2000；仪器仪表学报，24(6)：653-656，2003)，即将颜色三基色的中心波长作为测量的三个有效波长，进而利用比色测温法，实现温度场的计算；然而，实际上有效波长并非一个常量，其是随着测量物体的辐射光谱分布不同而不同，这种有效波长的简化处理方法，会给温度计算带来误差。上述方案a和b的应用均体现了一个主要局限性，即内嵌RGB彩色滤色器阵列提供三个颜色通道实现色彩的真实复现，其是基于标准人眼的光谱三刺激值进行设计的，因此，三个通道的波段响应特性一般是固定不变的，这往往会限制仪器测温的应用，从而无法自主选择合适的三个波段光谱响应函数以实现温度的优化测量。

另一种方案是采用单个黑白CCD作为光学测温传感器，将2个具有不同单色波长的滤色片交替放置于CCD前，测量物体在两个波长下的相对辐射强度值，再根据比色法进行温度计算(IEEE TRANS-ACTIONS ON INSTRUMENTATION AND MEASUREMENT，51(5)：990-995，2002)。这种方法应用时，虽然可以较灵活的选择所需要的波长，但需要满足在2个单色滤色片交替测量的时间段内物体温度场保持稳定的条件，对于瞬态温度场测量这种要求往往难以满足，从而具有应用上的局限性。

发明内容

本发明的目的是提供一种三CCD温度场测量装置及方法，以克服现有技术中基于单CCD温度场测量技术在应用上的局限性。

为了达到上述目的，本发明的技术方案提出一种三CCD温度场测量装置，该装置包括：光学镜头、分光棱镜组、滤色片、三个CCD面阵传感器及数据采集分析单元，

所述光学镜头，用于光学成像，将待测物体的辐射聚焦在所述分光棱镜组的入射面上；

所述分光棱镜组，将投射的辐射等分为三束不同方向的波段辐射，分别从三个出射面出射；

所述三个CCD面阵传感器，相同型号，分别设置于所述分光棱镜组的三个出射面处；所述滤色片则设于每个CCD面阵传感器的焦平面与分光棱镜组的出射面之间，并且三个滤色片具有线性非相关的光谱透过率分布函数；

所述三个CCD面阵传感器，分别对三路经过所述滤光片的出射辐射进行成像，获得三路光谱非相关的CCD图像；

所述数据采集分析单元，对三路CCD图像数据进行采集，并利用多光谱测温法进行温度场计算。

上述的三CCD温度场测量装置中，所述数据采集分析单元利用多光谱测温法进行温度场计算是根据以下方程组完成：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(1\right)$

其中，上角标(i，j)表示传感器成像焦平面上任一点的坐标；

分别表示所述三个CCD面阵传感器的相对辐射强度输出值，为已知量；Φ^i，j为所述CCD面阵传感器测量的非光谱因子，为未知量；λ_min～λ_max为分光棱镜组出射的三束波段辐射的响应波段；CCD(λ)为所述CCD面阵传感器的光谱响应特性函数，为已知量；τ₁(λ)、τ₂(λ)、τ₃(λ)分别表示各滤色片的光谱透过率分布函数与光学镜头、分光棱镜组的光谱透过率曲线合并而成的综合光谱透过率函数，为已知量；T^i，j表示待测物体在点(i，j)上的温度，为未知量；I_b，λ(T^i，j)为与待测物体相同温度T^i，j下的黑体光谱功率分布函数，仅与温度T^i，j相关；ε_λ(T^i，j)为待测物体的光谱发射率函数，包含两个以下的待定参数。

上述的三CCD温度场测量装置中，所述待测物体的光谱发射率函数采用线性函数表征如下：

ε_λ(T^i，j)＝a₀+a₁λ         (2)

则方程组(1)的各方程左边为三个已知量

右边包括T^i，j、Φ^i，j·a₀、Φ^i，j·a₁三个未知量，因此方程组(1)的求解是封闭的。

上述的三CCD温度场测量装置中，根据以下步骤建立测温数据库实现温度场的实时计算：

将式(2)代入方程组(1)得到方程组(3)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{h,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CDD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(3\right)$

定义公式中仅与温度T^i，j相关的6个积分量为β_1，0、β_1，1、β_2，0、β_2，1、β_3，0、β_3，0，其中

${\mathrm{\β}}_{m,n}={\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_m\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\λ}}^n\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ},m=\mathrm{1,2,3};n=\mathrm{0,1}.---\left(4\right)$

则将方程组(3)改写为方程组(5)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{1,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{1,1}}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{2,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{2,1}}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{3,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{3,1}}\end{array}\right.---\left(5\right)$

预先对于每个温度求解所述6个积分量并建立存储温度与6个积分量对应关系的测温数据库，之后利用所述测温数据库及三个已知量

对方程组(5)进行迭代运算，反演得到温度T^i，j。

上述的三CCD温度场测量装置中，所述待测物体为温度范围为1000K～3000K的具有连续辐射特性的高温物体。

上述的三CCD温度场测量装置中，所述分光棱镜组出射的三束波段辐射的响应波段均为600nm～1000nm，则所述三个滤色片的波段响应区间也为600nm～1000nm。

本发明的技术方案还提出一种应用如上所述装置的三CCD温度场测量方法，该方法包括：

利用所述光学镜头将待测物体的辐射聚焦在所述分光棱镜组的入射面上；

利用所述分光棱镜组将投射的辐射等分为三束不同方向的波段辐射，分别从三个出射面出射；

利用所述三个CCD面阵传感器分别对三路经过所述滤光片的出射辐射进行成像，获得三路光谱非相关的CCD图像；

利用所述数据采集分析单元对三路CCD图像数据进行采集，并利用多光谱测温法进行温度场计算。

上述的三CCD温度场测量方法中，所述利用多光谱测温法进行温度场计算是根据以下方程组完成：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(1\right)$

其中，上角标(i，j)表示传感器成像焦平面上任一点的坐标；

分别表示所述三个CCD面阵传感器的相对辐射强度输出值，为已知量；Φ^i，j为所述CCD面阵传感器测量的非光谱因子，为未知量；λ_min～λ_max为分光棱镜组出射的三束波段辐射的响应波段；CCD(λ)为所述CCD面阵传感器的光谱响应特性函数，为已知量；τ₁(λ)、τ₂(λ)、τ₃(λ)分别表示各滤色片的光谱透过率分布函数与光学镜头、分光棱镜组的光谱透过率曲线合并而成的综合光谱透过率函数，为已知量；T^i，j表示待测物体在点(i，j)上的温度，为未知量；I_b，λ(T^i，j)为与待测物体相同温度T^i，j下的黑体光谱功率分布函数，仅与温度T^i，j相关；ε_λ(T^i，j)为待测物体的光谱发射率函数，包含两个以下的待定参数。

上述的三CCD温度场测量方法中，所述待测物体的光谱发射率函数采用线性函数表征如下：

ε_λ(T^i，j)＝a₀+a₁·λ        (2)

则方程组(1)的各方程左边为三个已知量

右边包括T^i，j、Φ^i，j·a₀、Φ^i，j·a₁三个未知量，因此方程组(1)的求解是封闭的。

上述的三CCD温度场测量方法中，根据以下步骤建立测温数据库实现温度场的实时计算：

将式(2)代入方程组(1)得到方程组(3)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{h,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CDD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(3\right)$

定义公式中仅与温度T^i，j相关的6个积分量为β_1，0、β_1，1、β_2，0、β_2，1、β_3，0、β_3，1’其中

${\mathrm{\β}}_{m,n}={\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_m\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\λ}}^n\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ},m=\mathrm{1,2,3};n=\mathrm{0,1}.---\left(4\right)$

则将方程组(3)改写为方程组(5)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{1,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{1,1}}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{2,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{2,1}}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{3,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{3,1}}\end{array}\right.---\left(5\right)$

预先对于每个温度求解所述6个积分量并建立存储温度与6个积分量对应关系的测温数据库，之后利用所述测温数据库及三个已知量对方程组(5)进行迭代运算，反演得到温度T^i，j。

本发明的技术方案通过三CCD光谱融合采集待测物辐射强度信息，并利用适用性更强的多光谱测温法，可以实现温度场测量，应用范围更为广泛；通过采用测温数据库，提高了温度求解速度，可应用于实时在线温度计算；且技术方案实现简单，集成系统的成本不高、性能稳定，在高温检测等工业生产领域易于推广应用。

附图说明

图1为本发明三CCD温度场测量装置实施例结构图。

具体实施方式

以下实施例用于说明本发明，但不用来限制本发明的范围。

图1为本发明三CCD温度场测量装置实施例结构图，如图所示，本实施例的温度场测量装置包括：光学镜头11、分光棱镜组12、滤色片131～133、三个CCD面阵传感器141～143及数据采集分析单元15。其中，光学镜头11用于光学成像，将待测物体10的辐射聚焦在分光棱镜组12的入射面上，其可设计为定焦距或是变焦距的镜头。分光棱镜组12是由多个棱镜拼接而成，通过棱镜镀膜设计，使入射辐射在不同的棱镜结合面上经由透射和反射，最终分解为三束不同方向的波段辐射，分别由棱镜的三个不同侧面出射出去；三束射出辐射光响应波段均为λ_min～λ_max(本实施例为600nm1000nm)，强度、光谱分布均相等。三个CCD面阵传感器141～143为同一型号的黑白CCD传感器，分别设置于分光棱镜组12的三个出射面处；滤色片131～133则设于每个CCD面阵传感器的焦平面与分光棱镜组12的出射面之间，其在600nm～1000nm波段区间具有不同的光谱透过率分布，三个光谱透过率分布函数必需是线性非相关的，以满足多光谱测量的要求；上述三个滤色片131～133的光谱透过率分布函数可根据实际测量的需要进行设计，即考虑温度测量范围、分辨精度、测温误差、优化测量、CCD响应特性等诸多因素，三个滤色片131～133的光谱透过率分布函数可具有不同的中心波长、半宽度、峰值等。续上所述，三个CCD面阵传感器141～143的曝光时间设置相同，并且通过硬件帧采集控制保障传感器同步采集，从而分别对三路经过滤光片131～133的出射辐射进行成像，由于三个滤色片的光谱非相关特性，因此将获得三路光谱非相关的CCD图像。数据采集分析单元15则对三路CCD图像数据进行采集，并利用多光谱测温法进行温度场计算。

上述实施例中，分光棱镜组12、滤色片131～133、CCD面阵传感器141～143构成的分光成像系统，还需要进行光路校正，使得三个CCD面阵传感器141～143可对物体清晰成像，获得的图像能够点点对应，共同反应同一测量目标物体。至于所述的待测物体则是指温度范围在1000K～3000K之间的具有连续辐射特性的高温物体，其在CCD传感器的波段响应区间内的自发辐射强度要远远大于背景环境反射辐射强度的干扰，使得CCD传感器获得的测量信号能够直接定量反映高温物体自发辐射强度的大小。

上述的多光谱测温法为一种常用的辐射测温方法，通用的表述如下：具有连续辐射特性的物体发射率用一个关于波长的多项式函数予以描述，在多个波长下测量物体的辐射强度，结合光谱发射率模型，可以求得物体的温度。通常在有限的波段区间内，常用0阶(灰体)、1阶(线性)、2阶(二次函数)形式。在后续本实施例中将采用线性发射率函数的表述形式(灰体是其的特例形式)，发射率函数中共有两个待定系数，通过3路信号测量方程的反演计算，可以同时求得温度及两个待定系数。然而，任何包含两个或小于两个待定参数的发射率函数，均可以通过本发明技术方案获得的3路信号测量方程，计算温度数值，其仍然也归为多光谱测温法。

继续参考图1所示，应用上述本发明三CCD温度场测量装置实施例的测量方法过程具体如下所述。

首先，高温待测物体10的辐射通过光学镜头11，投射在分光棱镜组12上，分光棱镜组12将投射辐射等分为三束方向不同、强度与光谱分布相同的波段辐射，波段响应范围为600nm～1000nm。

接续，在三束分光光束的三个出射面处(棱镜上侧、下侧、右侧)，分别放置三个滤色片131～133，在600nm～1000nm波段区间滤色片131～133具有线性非相关的光谱透过率分布函数，以满足多光谱测量的要求；滤色片131～133的光谱透过率分布函数，可根据实际测量的需要进行设计，即考虑温度测量范围、分辨精度、测温误差、优化测量、CCD响应特性等诸多因素，三个滤色片的光谱透过率分布函数可具有不同的中心波长、半宽度、峰值等。例如，本实施例可以选择中心波长为750nm、850nm、950nm三个滤色片，半宽度均为10nm，峰值透过率均为75％。

然后，三个相同型号的、具有数字传输接口的黑白CCD面阵传感器141～143分别安置于滤色片131～133之后，三束光束经由滤色片131～133分别成像在CCD面阵传感器141～143上；CCD传感器141～143的曝光时间设置相同，并且通过帧控制保证三路CCD传感器同步采集。由于三个滤色片131～133的光谱非相关特性，由此CCD面阵传感器141～143将同时获得三路光谱非相关的CCD成像图像，图像灰度代表了目标物体相对辐射强度信息。

另外，分光棱镜组12、滤色片131～133、黑白CCD面阵传感器141～143共同构成的光学成像系统，需进行必要的光路校正，使得CCD面阵传感器141～143可对物体清晰成像，并且使获得的图像能够点点对应，实现非失真的融合匹配。

最后，三路图像信号传输至数据采集分析系统15，依据以下a)～e)所述的原理实现温度场的计算反演。

a).通过分光棱镜组、滤色片后，三路黑白CCD面阵传感器输出的辐射信号为：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(1\right)$

其中，上角标(i，j)表示传感器成像焦平面上任一点的坐标；

分别表示所述三个CCD面阵传感器的相对辐射强度输出值，为已知量；Φ^i，j为所述CCD面阵传感器测量的非光谱因子，为未知量；CCD(λ)为所述CCD面阵传感器的光谱响应特性函数，为已知量；τ₁(λ)、τ₂(λ)、τ₃(λ)分别表示各滤色片的光谱透过率分布函数与光学镜头、分光棱镜组的光谱透过率曲线合并而成的综合光谱透过率函数，为已知量；T^i，j，表示待测物体在点(i，j)上的温度，为未知量；I_b，λ(T^i，j)为与待测物体相同温度T^i，j下的黑体光谱功率分布函数，其仅与温度T^i，j相关。

b).以上方程组(1)中，ε_λ(T^i，j)为待测物体的光谱发射率函数。具有连续辐射特性的物体光谱发射率在数学上通常可以用多项式函数予以描述，一般而言，在一个有限的窄波段内(600nm～1000nm)应用时，采用线性函数表征光谱发射率具有很高的精度，如下式所示：

ε_λ(T^i，j)＝a₀+a₁·λ             (2)

光谱发射率函数中共有两个待定系数(a₀，a₁)，合并变量，方程组(1)中的三个方程中共有(T^i，j，Φ^i，j·a₀，Φ^i，j·a₁)三个未知量，因此温度的数学求解是封闭 $(V_1^{i,j},V_2^{i,j},V_3^{i,j})\⇔(T^{i,j},{\mathrm{\Φ}}^{i,j}\·a_0,{\mathrm{\Φ}}^{i,j}\·a_1),$ 从而能够继续实现基于3通道的多光谱辐射测温。

基于上述原理，在求解过程中，可继续如c)～e)所述建立测温数据库，以满足实时温度计算的需要。

c).将光谱发射率函数(2)代入方程组(1)，得到新的方程组(3)：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{h,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CDD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(3\right)$

d).定义公式中仅与温度T^i，j相关的6个积分量为β_1，0、β_1，1、β_2，0、β_2，1、β_3，0、β_3，1，其中

${\mathrm{\β}}_{m,n}={\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_m\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\λ}}^n\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ},m=\mathrm{1,2,3};n=\mathrm{0,1}.---\left(4\right)$

则将方程组(3)改写为方程组(5)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{1,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{1,1}}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{2,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{2,1}}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\mathrm{\β}}_{\mathrm{3,0}}+{\mathrm{\Φ}}^{i,j}\·a_1\·{\mathrm{\β}}_{\mathrm{3,1}}\end{array}\right.---\left(5\right)$

e).从d)中可以看出，积分运算仅与温度相关。因此，可以预先对于每个温度求解上述6个积分量并建立存储温度与6个积分量对应关系的测温数据库，之后利用所述测温数据库及三个已知量

对方程组(5)进行迭代运算，反演得到温度T^i，j。

参考上述实施例关于b)的描述，本领域的技术人员同时可以知悉，黑体、灰体测量情形(如ε(T^i，j)＝a₀)，均是光谱发射率函数式(2)的特例形式，这时式(1)的测量方程组写为：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(6\right)$

在式(6)中，仅有两个未知量(T^i，j，Φ^i，j·a₀)，利用三个通道测量量

数学求解时，未知量的数目小于方程组的数目，运用最小二乘法，可以使温度的求解误差更小。

由上述实施例可知，本发明中的光谱发射率函数用线性函数予以表述，但也不局限于此，任何包含两个或小于两个待定参数的发射率函数，均可以通过本发明获得的3路信号测量方程，计算求解逐点温度场T^i，j。例如，对于发光火焰温度场测量，具有一个待定参数的Hottel模型也已被广泛使用，用以表征火焰光谱发射率：

ε_λ＝1-exp(-KL/λ^α)              (7)

式(7)是一个半经验公式，其中，α是常数，系数K正比于粒子浓度，L是光程，其乘积KL可作为综合表征浓度的变量。这时式(1)的测量方程组可写为：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·(1-\exp (-\mathrm{KL}/{\mathrm{\λ}}^{\mathrm{\α}}))\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·(1-\exp (-\mathrm{KL}/{\mathrm{\λ}}^{\mathrm{\α}}))\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{600\mathrm{nm}}^{1000\mathrm{nm}}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·(1-\exp (-\mathrm{KL}/{\mathrm{\λ}}^{\mathrm{\α}}))\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(8\right)$

在式(8)中，同样也仅有两个未知量(T，KL)，利用三个通道测量量

可以迭代求解出温度值。具体过程与上述本发明的实施例类似，此处不再加以赘述。

综上所述，本发明基于三CCD融合的温度场测量技术，具有以下优点：

(1)与现有技术中非成像的单点辐射测温方法相比，本发明将点测量扩展到了二维场测量，获得了更为丰富的高温物体温度信息。

(2)与现有技术中CCD成像式温度场测量方法相比，本发明在测量时利用了滤色片及三通道分光成像处理方式，可以同时获取三路光谱非相关的辐射强度成像数据，三个测量通道的光谱响应函数可以根据需要进行设计，以更好地满足温度优化测量等多方面的要求。

(3)采用了适用性更强的多光谱测温法，测量物体不仅仅局限于传统的比色测温法所适用的范围，因而本发明的应用范围将更为广泛，具有更大的通用性。

(4)技术实现方案较为简单，三通道分光棱镜易于光学加工，CCD面阵传感器是成熟的商业产品，因此系统集成成本不高、性能稳定等，在高温检测等工业生产领域易于推广应用。

以上为本发明的最佳实施方式，依据本发明公开的内容，本领域的普通技术人员能够显而易见地想到一些雷同、替代方案，均应落入本发明保护的范围。

Claims (10)

1、一种三CCD温度场测量装置，其特征在于，该装置包括：光学镜头、分光棱镜组、滤色片、三个CCD面阵传感器及数据采集分析单元，

所述光学镜头，用于光学成像，将待测物体的辐射聚焦在所述分光棱镜组的入射面上；

所述分光棱镜组，将投射的辐射等分为三束不同方向的波段辐射，分别从三个出射面出射；

所述三个CCD面阵传感器，分别设置于所述分光棱镜组的三个出射面处；所述滤色片则设于每个CCD面阵传感器的焦平面与分光棱镜组的出射面之间，并且三个滤色片具有线性非相关的光谱透过率分布函数；

所述三个CCD面阵传感器，分别对三路经过所述滤光片的出射辐射进行成像，获得三路光谱非相关的CCD图像；

所述数据采集分析单元，对三路CCD图像数据进行采集，并利用多光谱测温法进行温度场计算。

2、如权利要求1所述的三CCD温度场测量装置，其特征在于，所述数据采集分析单元利用多光谱测温法进行温度场计算是根据以下方程组完成：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(1\right)$

其中，上角标(i，j)表示传感器成像焦平面上任一点的坐标；

分别表示所述三个CCD面阵传感器的相对辐射强度输出值，为已知量；Φ^i，j为所述CCD面阵传感器测量的非光谱因子，为未知量；λ_min～λ_max为分光棱镜组出射的三束波段辐射的响应波段；CCD(λ)为所述CCD面阵传感器的光谱响应特性函数，为已知量；τ₁(λ)、τ₂(λ)、τ₃(λ)分别表示各滤色片的光谱透过率分布函数与光学镜头、分光棱镜组的光谱透过率曲线合并而成的综合光谱透过率函数，为已知量；T^i，j表示待测物体在点(i，j)上的温度，为未知量；I_b，λ(T^i，j)为与待测物体相同温度T^i，j下的黑体光谱功率分布函数，仅与温度T^i，j相关；ε_λ(T^i，j)为待测物体的光谱发射率函数，包含两个以下的待定参数。

3、如权利要求2所述的三CCD温度场测量装置，其特征在于，所述待测物体的光谱发射率函数采用线性函数表征如下：

ε_λ(T^i，j)＝a₀+a₁·λ       (2)

则方程组(1)的各方程左边为三个已知量

、、

，右边包括T^i，j、Φ^i，j·a₀、Φ^i，j·a₁三个未知量，因此方程组(1)的求解是封闭的。

4、如权利要求3所述的三CCD温度场测量装置，其特征在于，根据以下步骤建立测温数据库实现温度场的实时计算：

将式(2)代入方程组(1)得到方程组(3)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\λ}\·I}_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(3\right)$

定义公式中仅与温度T^i，j相关的6个积分量为β_1，0、β_1，1、β_2，0、β_2，1、β_3，0、β_3，1，其中

${\mathrm{\β}}_{m,n}={\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_m\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\λ}}^n\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ},$ m＝1，2，3；n＝0，1.(4)

则将方程组(3)改写为方程组(5)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_0\·{\mathrm{\β}}_{\mathrm{1,0}}+{\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_1\·{\mathrm{\β}}_{\mathrm{1,1}}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_0\·{\mathrm{\β}}_{2,0}+{\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_1\·{\mathrm{\β}}_{2,1}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_0\·{\mathrm{\β}}_{\mathrm{3,0}}+{\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_1\·{\mathrm{\β}}_{\mathrm{3,1}}\end{array}\right.---\left(5\right)$

预先对于每个温度求解所述6个积分量并建立存储温度与6个积分量对应关系的测温数据库，之后利用所述测温数据库及三个已知量对方程组(5)进行迭代运算，反演得到温度T^i，j。

5、如权利要求1所述的三CCD温度场测量装置，其特征在于，所述待测物体为温度范围为1000K～3000K的具有连续辐射特性的高温物体。

6、如权利要求1～5任一项所述的三CCD温度场测量装置，其特征在于，所述分光棱镜组出射的三束波段辐射的响应波段均为600nm～1000nm，则所述三个滤色片的波段响应区间也为600nm～1000nm。

7、一种应用如权利要求1所述装置的三CCD温度场测量方法，其特征在于，该方法包括：

利用所述光学镜头将待测物体的辐射聚焦在所述分光棱镜组的入射面上；

利用所述分光棱镜组将投射的辐射等分为三束不同方向的波段辐射，分别从三个出射面出射；

利用所述三个CCD面阵传感器分别对三路经过所述滤光片的出射辐射进行成像，获得三路光谱非相关的CCD图像；

利用所述数据采集分析单元对三路CCD图像数据进行采集，并利用多光谱测温法进行温度场计算。

8、如权利要求7所述的三CCD温度场测量方法，其特征在于，所述利用多光谱测温法进行温度场计算是根据以下方程组完成：

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\ϵ}}_{\mathrm{\λ}}\left(T^{i,j}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(1\right)$

其中，上角标(i，j)表示传感器成像焦平面上任一点的坐标；

分别表示所述三个CCD面阵传感器的相对辐射强度输出值，为已知量；Φ^i，j为所述CCD面阵传感器测量的非光谱因子，为未知量；λ_min～λ_max为分光棱镜组出射的三束波段辐射的响应波段；CCD(λ)为所述CCD面阵传感器的光谱响应特性函数，为已知量；τ₁(λ)、τ₂(λ)、τ₃(λ)分别表示各滤色片的光谱透过率分布函数与光学镜头、分光棱镜组的光谱透过率曲线合并而成的综合光谱透过率函数，为已知量；T^i，j表示待测物体在点(i，j)上的温度，为未知量；I_b，λ(T^i，j)为与待测物体相同温度T^i，j下的黑体光谱功率分布函数，仅与温度T^i，j相关；ε_λ(T^i，j)为待测物体的光谱发射率函数，包含两个以下的待定参数。

9、如权利要求8所述的三CCD温度场测量方法，其特征在于，所述待测物体的光谱发射率函数采用线性函数表征如下：

ε_λ(T^i，j)＝a₀+a₁·λ       (2)

则方程组(1)的各方程左边为三个已知量

右边包括T^i，j、Φ^i，j·a₀、Φ^i，j·a₁三个未知量。

10、如权利要求9所述的三CCD温度场测量方法，其特征在于，根据以下步骤建立测温数据库实现温度场的实时计算：

将式(2)代入方程组(1)得到方程组(3)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_1\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_2\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·\mathrm{\λ}\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·a_0\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\\ +{\mathrm{\Φ}}^{i,j}\·a_1\·{\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_3\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\λ}\·I}_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ}\end{array}\right.---\left(3\right)$

定义公式中仅与温度T^i，j相关的6个积分量为β_1，0、β_1，1、β_2，0、β_2，1、β_3，0、β_3，1，其中

${\mathrm{\β}}_{m,n}={\∫}_{{\mathrm{\λ}}_{\min }}^{{\mathrm{\λ}}_{\max }}{\mathrm{\τ}}_m\left(\mathrm{\λ}\right)\·\mathrm{CCD}\left(\mathrm{\λ}\right)\·{\mathrm{\λ}}^n\·I_{b,\mathrm{\λ}}\left(T^{i,j}\right)\mathrm{d\λ},m=\mathrm{1,2,3};n=\mathrm{0,1}.---\left(4\right)$

则将方程组(3)改写为方程组(5)，

$\left\{\begin{array}{c}{V}_1^{i,j}={\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_0\·{\mathrm{\β}}_{\mathrm{1,0}}+{\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_1\·{\mathrm{\β}}_{\mathrm{1,1}}\\ V_2^{i,j}={\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_0\·{\mathrm{\β}}_{2,0}+{\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_1\·{\mathrm{\β}}_{2,1}\\ V_3^{i,j}={\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_0\·{\mathrm{\β}}_{\mathrm{3,0}}+{\mathrm{\Φ}}^{i,j}\·{\mathrm{\α}}_1\·{\mathrm{\β}}_{\mathrm{3,1}}\end{array}\right.---\left(5\right)$

预先对于每个温度求解所述6个积分量并建立存储温度与6个积分量对应关系的测温数据库，之后利用所述测温数据库及三个已知量

对方程组(5)进行迭代运算，反演得到温度T^i，j。

Images