本文内容转载自《红外技术》2021年第2期，版权归《红外技术》编辑部所有。本文内容不含参考文献，如有需要请参考原论文。
付万超，范春利，杨立
海军工程大学动力工程学院
摘要：利用红外热像仪测温需先设定被测表面的法向发射率，该发射率通常为定值。而当热像仪处于被测点的天顶角大于50︒的位置范围时，由于被测点定向发射率的变化，必造成这些点的测温误差。对于非平表面，这样的点大量存在。因此，必须对其测温结果进行修正。本文针对使用单目红外热像仪测量非平表面温度时由于各点定向发射率的变化引起的测量误差进行研究，并依据物体表面定量发射率的变化规律，给出了测量点的温度修正系数。同时，通过点云三维建模，利用热像仪的几何成像原理推导出红外热像图与实际被测表面中点与点的对应关系，给出了通过红外热像仪测量非平表面的温度分布的误差修正方法。实验证明了该方法的有效性。
关键词：红外无损检测；红外热像仪；测温误差；温度修正；非平表面
0 引言
在工业生产中，及时掌握设备的工作状态是十分必要的。好多情况下，物体表面温度可以间接地反映设备的工作状态，因此开展依据物体表面温度对工业设备运行状态进行检测和监测的理论与方法研究具有重要的意义。
对于设备表面温度分布的测量，最常用的设备为红外热像仪。通过红外热像仪进行测温，一般要求正对测量表面。然而在狭小空间或由其他原因造成可达性较差的位置进行测量时，如舰船的机舱，常常无法满足红外热成像测温的要求，因此就需要人们对特殊测量条件下所测得的温度分布进行修正。造成红外热像测温误差的因素很多，相关的研究包括红外热像仪的测温的影响因素分析及误差的计算、太阳、大气等环境的影响，周围高温物体的影响、测量距离及视角的影响、探测器温度的影响等。对于本文讨论的非平设备表面，尤其在热像仪处于被测点的天顶角的50︒范围之外时，发射率将发生明显的变化。因此，对于设为恒定表面发射率的热像仪，必将会产生显著的测温误差。因此，开展非平表面红外温度测量的修正研究具有重要的应用价值。
本文针对利用红外热像仪大角度测温或对非平表面进行测温因无法或部分无法处于表面各点定向发射率恒定的天顶角范围内时热像仪的测温误差进行研究，并基于物体表面定向发射率的变化规律，通过对所测表面进行3D建模，建立实际表面和热像图之间的坐标对应关系，提出热像仪所测温度分布的修正方法。方法通过实验进行验证。
1 红外测温及修正系数的计算
根据红外热像仪温度计算公式：
式中：ε物体表面发射率；τa为大气的光谱透射率；Tr为热像仪指示的辐射温度；εa为大气辐射率；Tu为环境温度。在3～5μm，n=9.2554;在8～12μm，n=3.9889。
由式(1)可知，温度测量值与被测物体的发射率ε密切相关。在测量复杂物体表面时，需在热像仪上将发射率设为一固定值，通常为该表面材料的法向发射率。天顶角是某个方向和天顶方向的夹角，通常记为θ（见图1）。本文中利用天顶角来描述热像仪测量物体表面中点的温度时，热像仪成像中心点相对于物体表面的位置关系。在一次测量时，热像仪所处的位置，对于非平表面上不同的点具有不同的天顶角θ。当角θ超出一定范围（如大于50︒）时，定向发射率的变化必然对红外测温的准确性产生影响。
图1 热像仪机位与其对应的天顶角
根据菲涅尔定律，光谱方向反射率表示为：
式中：
表示材料的复折射系数。
根据Snell定律，折射角与入射角满足：
根据Kirchhoff定律：
据表1中铝(AI)、黄铜(Brass)、聚乙烯(PS)、聚苯乙烯(PE)4种材料分别在3～5μm与8～12μm工作波段的复折射率可计算出其定向发射率随天顶角变化的理论值。
表1 材料各波段复折射率
式(4)的计算结果如图2所示，可观察到金属与非导电体定向发射率随天顶角的变化规律：各类材料定向发射率在被测点天顶角大于50︒附近开始出现明显的变化。
在工程测量中，需在测量前设定被测物体表面的发射率ε1，通常设为其法向发射率乱。对于相对平整的表面，可使用ε1代表所有点的发射率。当热像仪处于被测点的天顶角大于50︒的位置范围时，需采用该角度的定向发射率。对于不同角度的定向发射率εθ，根据(1)式可知实际温度T0与测量温度T1，的关系为：
图2 定向发射率随天顶角的变化
定义修正系数为：
则可得修正系数曲线如图3 所示。
由图3 观察得出，对于金属表面温度测量结果的修正系数在天顶角大于50︒时减小，在90︒附近上升迅速，而非导电体在天顶角大于70︒时迅速增加，最大时测温误差可达50%。因此，需要对红外热像图中不同位置的温度值依据其对应实际点的修正系数进行修正。
图3 温度修正系数随天顶角的变化
2 被测表面各点天顶角的计算
如前所述，要计算被测表面各个点的修正系数，必须首先计算各个点红外测温时热像仪所处的天顶角。
2.1 坐标与法向量的获取
利用3D建模或3D扫描技术可获取被测物体的空间结构，存储格式为.stl 文件。将所得.stl文件使用Meshlab软件进行网格划分处理后可获得该被测物体的点云数据，存储为.xyz文件。.xyz文件为一组六维向量，前三维是点云中各点世界坐标系下坐标
，后三维是该点的法向量坐标
。
以一锅炉为例，绘制出其点云坐标和法向量如图4所示。
图4 某锅炉点云坐标与法向量示例
本文将建模过程中的模型原点作为世界坐标系原点。可知热像仪机位中心坐标
。
2.2 天顶角的计算
由热像仪机位中心坐标与空间坐标可得测量方向向量：
则可计算天顶角为：
由空间几何关系可知，对于天顶角小于90︒的被测点，由于透视关系无法在热像仪上成像，因此可得可视部分如图5 所示。
其中坐标为(-0.75,0.12,0.1)的点为红外热像仪中心所在点，黄色部分为锅炉表面在该位置关系下热像仪的可视点，蓝色为不可视点。
图5 锅炉的红外热像仪可视点
2.3 温度修正系数的分布
由该位置关系可计算得出各点对应天顶角。由于实际工程应用中，通常在金属表面填涂一层非导电涂料，且非导电体在法向发射率相较于导电体更大，理论的光谱定向发射率与实际相差不大，因此以8～13μm的PS（ Poly styrene）材料作为被测物体材料。在图4 的位置关系下，计算得出该锅炉的修正系数c的频率分布直方图如图6 所示。经计算可得，该位置关系下修正系数大于1.2 的点占比44.1%。由此可见，对实际非平表面测温时的定向发射率进行修正是非常必要的。
图6 不同温度修正系数对应的被测点所占总被测点的比例
3 被测点在图像的映射
天顶角与成像点的实际位置和其法向量密切相关。因此，针对红外热像仪所测温度的修正需明确各成像点在世界坐标中的实际位置。
世界坐标系下被测点pw 经红外热像仪投影至平面像素点pp＝[u,v]T，在成像过程中经过两次坐标变换，分别是由世界坐标系到热像仪坐标系的变换和由热像仪坐标系到像素坐标系的变换，两个变换分别通过外参矩阵和内参矩阵来实现。
3.1 外参矩阵的获取
世界坐标pw 到热像仪坐标pi 的转换关系为：
为热像仪的外参矩阵，Ri 为热像仪的基坐标在世界坐标系的坐标，描述热像仪在世界坐标系中的姿态，I为热像仪中心在世界坐标系中的位置。
获取热像仪的外参矩阵需确定热像仪姿态，热像仪在成像过程中一般固定角度。假设热像仪在实际测量中热像仪x 轴始终与水平面平行，仅存在俯角γ与转角β。则可知当热像仪对准被测物体中心
时，热像仪中心指向被测点中心向量
为热像仪坐标z 轴方向，热像仪坐标系x 轴与水平面平行，y 轴方向垂直于x 轴和z 轴。由数学关系可得3 个方向的标准基为：
计算可得Ri＝[e1,e2,e3]为热像仪在世界坐标系旋转到当前姿态的旋转矩阵。
当热像仪角度变化时，可根据热像仪的俯角与转角求得热像仪坐标系z 轴与水平面的交点：
将上述计算过程中被测点中心
替换为p0即可根据热像仪实际姿态确定外参矩阵。
3.2 内参矩阵的获取
在世界坐标转换为热像仪坐标后，需将热像仪坐标转换为像素坐标，完成三维到二维的转变。
热像仪坐标pi 到像素坐标pp 的转换关系为：
式中：
为热像仪的内参矩阵，
为尺度因子，描述成像位置与热像仪中心的距离。
内参矩阵与热像仪的型号和图像的分辨率有关，kx、ky 为热像仪x 与y 方向的像素焦距，cx、cy 为图片主点偏移，内参矩阵均以像素为单位。
fx、fy 为热像仪x 与y 方向的毫米焦距，对于毫米焦距和像素焦距的转换关系有：
式中：ccd_size 是热像仪厂商给出的CCD 尺寸，以英寸为单位；image_x、image_y 是图像x、y 方向的像素数，与成像的分辨率有关。
3.3 热像仪标定
在工程应用中，外参矩阵主要通过机位的计算得出，内参矩阵经热像仪标定得出。
对于热像仪标定，一般使用张氏标定法。通过Matlab 软件自带的热像仪标定程序camera calibrator可对单目热像仪进行标定。可见光相机标定需打印棋盘格，对于红外热像仪，将棋盘格用金属板代替。黑色部分使用隔热胶带，白色部分使用金属板或锡纸。
在程序中输入标定板在不同角度下的红外热像图，可得到热像仪中心与标定板相对位置和热像仪内参、外参、图片尺寸等参数。本文主要使用标定得到的热像仪内参矩阵K。
4 温度修正计算流程
4.1 温度修正系数矩阵
总体的温度修正计算流程如图7 所示。
据式(7)计算出天顶角，查得该天顶角下的发射率，然后计算式(6)得到世界坐标系下各测温点的修正系数c(pw)。经热像仪成像后，测温点由世界坐标pw映射到成像平面坐标pp＝[u,v]T。pw 到pp 的转换关系由3.3 节标定的热像仪矩阵计算得出。经上述转换关系可得测温点在二维平面坐标对应的修正系数c(pp)。
热像仪拍摄的温度图像将温度以矩阵的形式存储，因此，应当将pp从离散的点云转换为与热像图格式相同的pi,j，其中，pi,j 代表像素块(i,j)的坐标。
对于修正系数在像素块范围内梯度较小的点，该像素块的修正系数取像素块内的平均值。对在热像图上像素位置为(i, j)的点，该像素的修正系数为：
当像素内修正系数梯度变化较大的，可采用该像素内修正系数最大值附近点取平均值的方法，在此像素块内凸显该形状非光滑点。
4.2 成像温度修正
对于红外热像仪，其所测温度矩阵为T0(pp)。修正后的温度矩阵为所测温度矩阵与温度修正系数矩阵做基本积（哈达玛积）：
所得温度T(pp)即修正后的温度。
5 实验验证
本文通过一个热水桶表面的温度测量实验来验证该算法有效性。首先使用Solidworks软件进行几何建模，获得被测物体几何信息，经Meshlab 软件求得热像仪可视表面的点云及各点处的法向量，再经Matlab 计算出各点温度的修正系数以及各点对应的热像图中的坐标，最后将红外热像仪直接测得的温度分布信息导入Matlab 进行修正得到水桶表面的修正后的温度。
5.1 红外测温与热像仪标定
本文选用Fluke Ti400 型号热像仪进行测温，测温表面为一聚乙烯塑料桶侧面，塑料桶底部外径为18 cm，顶部外径为20.5 cm，高为20 cm。塑料桶中为热水， 热像仪在世界坐标系中的坐标为C=[-0.75,0.12,0.10]T，热像仪正对塑料桶中心，位置关系如图8 所示。测得塑料桶外表面可视部分的红外热像图见图9。
图7 温度修正流程图
图8 热像仪机位
图9 红外热像图与温度截线
使用热像仪拍摄红外标定板（如图10 所示），获取标定板的红外热像图，将热像图输入标定程序，获得热像仪内参矩阵K。
图10 红外热像仪的标定
5.2 测温结果的修正
通过Solidworks 软件建模并将得到的水桶表面的云数据导入Matlab 软件中，可得到水桶表面与热像仪在空间中的位置关系，如图11 所示（黄色点为被测表面的可视点，蓝色为被测表面的非可视点）。本实验中，略去了接近桶口的外径较大的部分。
图11 热像仪与被测表面的位置关系
依据热像仪与可视点的空间位置关系，据式(7)可计算得各可视点的天顶角，据图2 可得不同天顶角时，所对应的发射率。然后据式(6)算得每个可视点的修正系数。如4.1 节所述，对空间中可视点进行热像仪矩阵的变换，映射至平面可得出热像图中各点对应的修正系数，见图12。进而依据式(15)得到热像图中各点修正后的温度分布。
图12 修正系数分布
为了更好地显示本文方法对红外热像仪所测得的温度分布，本节在所研究的可视表面高度方向的中心，取与图9 中相同位置的截线进行展示，修正系数如图13 所示。将所得修正系数与成像截线温度相乘得到修正后的截线温度，修正前后温度的对比如图14所示。
由图14 可见，修正前最大误差近40℃，真实温度为335℃，相对误差约为12.4%；修正后，集中在两侧大角度测量时的误差不超过5.0℃，相对误差小于1.5%。从此可见该方法的有效性。值得说明的是图13 中修正系数分布曲线并不完全对称，主要原因是该图的横坐标为非常微小的像素点，研究过程中像素点的选取没有完全对称所致。更多的平行截线结果均表明，该修正方法较好地消除了非平表面在红外热像测温时部分因为测量角度过大而引起的发射率变化所造成的误差。
图13 截线上各点的修正系数
图14 截线上温度分布修正前和修正后的对比
6 结论
针对使用红外热像仪测量非平表面温度时不同被测量点的定向发射率的变化对测温结果的影响，本文依据物体表面定量发射率的变化规律，给出了各个测量点的修正系数。同时，通过点云三维建模，利用热像仪通用几何成像原理推导出红外热像图像素与实际物体的被测表面之间各个测量点的对应关系，给出了红外热像测温过程中非平表面上的大天顶角被测点的温度测量结果的误差修正方法。实验证明了该方法的有效性。
