一种确定红外热像仪镜头平面与被测平面之间夹角的方法

摘要

本发明公开了一种确定红外热像仪镜头平面与被测平面之间夹角的方法。本发明主要原理是基于几何光学。在镜头与规则图形所在平面成角度测量时，规则图形所成的像会发生形变。所成像的形变程度随着测量角度的变化而变化。通过检测像的形变程度确定热像仪镜头平面与被测平面所成的测量角度。通过本发明使规则图形图靶附着在待测区域附近，用热像仪进行温度测量并存储下热像，检测热像上规则图形所成的像，得到规则图形的像的形变程度，进而得出热像仪与被测平面间的角度关系，以便对热像仪在一定角度上测量被测表面的情况进行角度补偿。测量时无需接触被测物表面，并且与热像仪测量温度同时进行，可以简单快速的确定热像仪的测量角度。

说明书

技术领域

本发明涉及一种角度测量方法，尤其是涉及一种确定红外热像仪镜头平面(焦平面)与被测平面之间测量角度的方法。

背景技术

红外热像仪以其非接触、快速、大面积测量温度的优点已被广泛的应用于工业现场关键点的温度测量和监测。红外热像仪镜头与被测平面之间的夹角，即红外热像仪的测量角度直接影响着红外热像仪的测温精度。但是由于工业现场的具体环境比较复杂，无法保证红外热像仪正对目标平面进行测量，从而产生热像仪测量角度引起的误差。现有的红外热像仪不具备测量热像仪与被测平面夹角的功能，传统的接触式测量方法也不适用于工业现场的角度测量，因此不能得到热像仪的测量角度，进而补偿测量角度引起的误差。本发明涉及测量红外热像仪与目标平面之间夹角的测量方法。

发明内容

鉴于上述存在的问题，本发明的目的在于提供一种测量红外热像仪镜头平面与被测平面之间所成的测量角度的方法。

本发明的主要原理是基于几何光学的原理。在镜头与规则图形图靶(圆形或正多边形)所在的被测平面成α角度测量时，规则图形图靶所成的像会发生形变。所成像的形变程度随着测量角度的变化而变化。通过检测像的形变程度确定热像仪镜头平面与被测平面所成的测量角度。

本发明采用的技术方案是：该方法依次包含以下步骤：

1)将规则图形图靶贴在被测平面测量区域；

2)将红外热像仪安装在三脚架上，放置于被测平面前方，规则图形图靶材料发射率与被测平面发射率不同，对准规则图形图靶对被测平面测量区域进行测量，并储存红外图像；

3)通过图像检测方法，得到红外图像中的六个特殊点；

4)由检测到的六个特殊点求出一个椭圆的方程；

5)将椭圆做旋转运算，将旋转后的椭圆整理成椭圆的标准方程形式，求出椭圆的长轴、短轴的值；

6)由椭圆短轴与长轴之比计算热像仪镜头平面与被测平面的夹角。

所述的规则图形图靶为多于六边的正多边形或圆形。

所述的通过图像检测方法，得到红外图像中的六个特殊点，若所述的规则图形图靶为多于六边的正多边形，则通过图像检测方法，检测得到多边形的任意六个顶点坐标；若所述的规则图形图靶为圆形，则通过图像检测方法，检测得到图形边缘上的任意六个点的坐标。

所述的由检测到的六个特殊点求出一个椭圆方程，该椭圆方程由六个公式Ax²+2Bxy+Cy²+Dx+Ey+F＝0组成的线性方程组求出，其中A、B、C、D、E、F均为椭圆方程系数。

所述的将椭圆做旋转运算包含以下步骤；

确定旋转角度旋转角度通过如下公式计算，

其中，A、B、C为第4)步骤中求出的椭圆方程的系数，

若B为0，则为0，不需要进行旋转；若B不为0，将椭圆方程做如下变换：

得到旋转后的椭圆方程。

第5)步骤中椭圆的标准方程形式为

$\frac{{(x^{\prime }-u)}^2}{a^2}+\frac{{(y^{\prime }-v)}^2}{b^2}=1$

其中x′、y′为变换后的坐标；u、v为椭圆的中心；a、b分别为椭圆的长轴或短轴，由a、b的长短决定哪个是长轴，由椭圆方程系数A、B、C、D、E、F及表示。

热像仪镜头平面与被测平面的夹角α根据下式计算得到：

Cosα＝a/b，当a为短轴，b为长轴；

Cosα＝b/a，当b为短轴，a为长轴。

本发明具有的有益效果是：

1.本发明采用非接触方式进行测量，无需对现有红外热像仪装置进行改造，可确定红外热像仪与被测物表面的测量角度；

2.本发明利用被测物表面特点，以几何光学为理论基础，在红外热像仪测温的同时测量其与被测物平面的角度关系，快速简易。

附图说明

图1是热像仪测量示意图。

图2是热像仪镜头平面与被测平面相对位置示意图。

图3是正六边形图靶及其像的示意图。

图4是圆形图靶及其像的示意图。

图中：1-红外热像仪，2-被测平面，3-三脚架，4-规则图形图靶，5-正六边形图靶，6-正六边形图靶的外接圆，7-正六边形图靶在红外图像中的像，8-正六边形图靶在红外热像中的像的外接椭圆，9-圆形图靶，10-圆形图靶在红外热像中的像。

具体实施方式

下面结合附图和具体实施例对本发明进行进一步说明。

如图1所示，红外热像仪1安装在三脚架3上，对附着在被测平面2上的规则图形图靶4进行测量。图2所示为红外热像仪镜头平面与被测平面相对位置的示意图，红外热像仪1与被测平面2之间的夹角为α。图3所示为红外热像仪镜头平面与被测平面成α角度测量时，规则图形图靶4在红外图像中的像。其中，图3a所示，正六边形图靶5具有唯一的外接圆6，正六边形图靶5在红外图像中的像如图3b所示，正六边形图靶在红外图像中的像7具有唯一的正六边形图靶在红外热像中的像的外接椭圆8；图4为圆形图靶及其像的示意图，图4a所示为圆形图靶9，图4b所示为圆形图靶在红外热像中的像10为椭圆形。

实施例1：规则图型图靶为正六边形

当红外热像仪镜头平面与被测平面成α角度测量时，此测量角度α由以下步骤确定：

第一、将正六边形图靶粘贴在被测平面测量区域。

第二、将红外热像仪安装在三脚架上，使热像仪镜头对准被测平面上的正六边形图靶，对被测平面进行测量，并存储下红外图像。

第三、通过图像检测方法，检测红外图像中正六边形图靶的像7的六个顶点的坐标(x₁，y₁)，(x₂，y₂)，(x₃，y₃)，(x₄，y₄)，(x₅，y₅)，(x₆，y₆)。

第四、由检测到的红外图像中正六边形图靶的像7的六个顶点的坐标，求出正六边形图靶的像的唯一外接椭圆8的椭圆方程，该椭圆方程的一般形式Ax²+2Bxy+Cy²+Dx+Ey+F＝0根据下列六个方程求出：

$\left(\begin{array}{c}{{\mathrm{Ax}}_1}^2+{2\mathrm{Bx}}_1{y}_1+{{\mathrm{Cy}}_1}^2+{\mathrm{Dx}}_1+{\mathrm{Ey}}_1+F=0\\ {{\mathrm{Ax}}_2}^2+{2\mathrm{Bx}}_2{y}_2+{{\mathrm{Cy}}_2}^2+{\mathrm{Dx}}_2+{\mathrm{Ey}}_2+F=0\\ {{\mathrm{Ax}}_3}^2+{2\mathrm{Bx}}_3{y}_3+{{\mathrm{Cy}}_3}^2+{\mathrm{Dx}}_3+{\mathrm{Ey}}_3+F=0\\ {{\mathrm{Ax}}_4}^2+{2\mathrm{Bx}}_4{y}_4+{{\mathrm{Cy}}_4}^2+{\mathrm{Dx}}_4+{\mathrm{Ey}}_4+F=0\\ {{\mathrm{Ax}}_5}^2+{2\mathrm{Bx}}_5{y}_5+{{\mathrm{Cy}}_5}^2+{\mathrm{Dx}}_5+{\mathrm{Ey}}_5+F=0\\ {{\mathrm{Ax}}_6}^2+{2\mathrm{Bx}}_6{y}_6+{{\mathrm{Cy}}_6}^2+{\mathrm{Dx}}_6={\mathrm{Ey}}_6+F=0\end{array}\right)$

其中A、B、C、D、E、F均为椭圆方程系数，由上式求出。

第五、旋转椭圆运算，Ax²+2Bxy+Cy²+Dx+Ey+F＝0为椭圆的一般方程，若B不为0，设椭圆旋转的角度为按照式求出旋转角

将椭圆做旋转运算，即椭圆方程坐标做如下变换：

得到旋转后的椭圆方程。

若B为0，则不需要进行旋转。

第六、将旋转后的椭圆整理成标准形式，即：

$\frac{{(x^{\prime }-u)}^2}{a^2}+\frac{{(y^{\prime }-v)}^2}{b^2}=1$的形式，

其中x′、y′为变换后的坐标；u、v为椭圆的中心；a、b分别为椭圆的长轴或短轴，由a、b的长短决定哪个是长轴，可由椭圆方程系数A、B、C、D、E、F及表示。

第七、计算夹角值，热像仪镜头平面与被测平面的夹角α可根据下式计算得到：

Cosα＝a/b，当a为短轴，b为长轴；

Cosα＝b/a，当b为短轴，a为长轴。

实施例2：规则图型图靶为圆形

当红外热像仪镜头平面与被测平面成α角度测量时，此测量角度α由以下步骤确定：

第一、将圆形图靶粘贴在被测平面测量区域。

第二、将红外热像仪安装在三脚架上，使热像仪镜头对准被测平面上的圆形图靶，对被测平面进行测量，并存储下红外图像。

第三、通过图像检测方法，检测红外图像中圆形图靶的像10的边缘曲线，此边缘曲线为椭圆曲线。取曲线上的任意六个点，坐标分别为(x₁，y₁)，(x₂，y₂)，(x₃，y₃)，(x₄，y₄)，(x₅，y₅)，(x₆，y₆)。

第四、由上述第三步中的得到的六个点的坐标，求出圆形图靶的像10的椭圆方程，该椭圆方程的一般形式Ax²+2Bxy+Cy²+Dx+Ey+F＝0根据下列六个方程求出：

$\left(\begin{array}{c}{{\mathrm{Ax}}_1}^2+{2\mathrm{Bx}}_1{y}_1+{{\mathrm{Cy}}_1}^2+{\mathrm{Dx}}_1+{\mathrm{Ey}}_1+F=0\\ {{\mathrm{Ax}}_2}^2+{2\mathrm{Bx}}_2{y}_2+{{\mathrm{Cy}}_2}^2+{\mathrm{Dx}}_2+{\mathrm{Ey}}_2+F=0\\ {{\mathrm{Ax}}_3}^2+{2\mathrm{Bx}}_3{y}_3+{{\mathrm{Cy}}_3}^2+{\mathrm{Dx}}_3+{\mathrm{Ey}}_3+F=0\\ {{\mathrm{Ax}}_4}^2+{2\mathrm{Bx}}_4{y}_4+{{\mathrm{Cy}}_4}^2+{\mathrm{Dx}}_4+{\mathrm{Ey}}_4+F=0\\ {{\mathrm{Ax}}_5}^2+{2\mathrm{Bx}}_5{y}_5+{{\mathrm{Cy}}_5}^2+{\mathrm{Dx}}_5+{\mathrm{Ey}}_5+F=0\\ {{\mathrm{Ax}}_6}^2+{2\mathrm{Bx}}_6{y}_6+{{\mathrm{Cy}}_6}^2+{\mathrm{Dx}}_6={\mathrm{Ey}}_6+F=0\end{array}\right)$

其中A、B、C、D、E、F均为椭圆方程系数，由上式求出。

第五、旋转椭圆运算，Ax²+2Bxy+Cy²+Dx+Ey+F＝0为椭圆的一般方程，若B不为0，设椭圆旋转的角度为按照式求出旋转角

将椭圆做旋转运算，即椭圆方程坐标做如下变换：

得到旋转后的椭圆方程。

若B为0，则不需要进行旋转。

第六、将旋转后的椭圆整理成标准形式，即：

$\frac{{(x^{\prime }-u)}^2}{a^2}+\frac{{(y^{\prime }-v)}^2}{b^2}=1$的形式，

其中x′、y′为变换后的坐标；u、v为椭圆的中心；a、b分别为椭圆的长轴或短轴，由a、b的长短决定哪个是长轴，可由椭圆方程系数A、B、C、D、E、F及表示。

第七、计算夹角值，热像仪镜头平面与被测平面的夹角α可根据下式计算得到：

Cosα＝a/b，当a为短轴，b为长轴；

Cosα＝b/a，当b为短轴，a为长轴。

上述描述仅为本发明的具体实施例，本发明不限于上述实施例。