基于改进YOLOv3的避雷器红外图像故障检测方法

胡泰山, 刘浩, 刘刚, 梅琪, 马御棠, 廖民传

胡泰山, 刘浩, 刘刚, 梅琪, 马御棠, 廖民传. 基于改进YOLOv3的避雷器红外图像故障检测方法[J]. 红外技术, 2023, 45(11): 1256-1261.
引用本文: 胡泰山, 刘浩, 刘刚, 梅琪, 马御棠, 廖民传. 基于改进YOLOv3的避雷器红外图像故障检测方法[J]. 红外技术, 2023, 45(11): 1256-1261.
HU Taishan, LIU Hao, LIU Gang, MEI Qi, MA Yutang, LIAO Minchuan. Infrared Image Fault Detection Method of Arrester Based on Improved YOLOv3[J]. Infrared Technology , 2023, 45(11): 1256-1261.
Citation: HU Taishan, LIU Hao, LIU Gang, MEI Qi, MA Yutang, LIAO Minchuan. Infrared Image Fault Detection Method of Arrester Based on Improved YOLOv3[J]. Infrared Technology , 2023, 45(11): 1256-1261.

基于改进YOLOv3的避雷器红外图像故障检测方法

基金项目: 

南方电网公司科技项目 YNKJXM20191243

详细信息
    作者简介:

    胡泰山(1990-),男,湖北监利人,硕士,工程师,研究方向:电力系统过电压及绝缘配合研究工作。E-mail: artist_hts@163.com

  • 中图分类号: TP391.41

Infrared Image Fault Detection Method of Arrester Based on Improved YOLOv3

  • 摘要: 针对现有的金属氧化物避雷器(Metal Oxide Arrester,MOA)红外图像故障检测方法存在识别精度低、检测速度较慢的问题,提出一种基于改进YOLOv3的MOA红外图像故障检测方法。首先,以Darknet19网络代替YOLOv3原始的Darknet53网络,并在特征学习时针对样本中不同MOA长宽比例,通过K-means聚类算法对MOA图像中的目标帧进行分析,重新聚类样本中心锚点框,得到合适的锚框数目和大小。最后,利用改进YOLOv3模型完成MOA红外图像故障检测。实验结果表明,改进的YOLOv3模型识别精度达到96.3%,识别速度为6.75 ms。
    Abstract: Aiming at the problems of low recognition accuracy and slow detection speed of existing metal oxide arrester (MOA) infrared image fault detection methods, a MOA infrared image fault detection method based on improved YOLOv3 is proposed. Firstly, darknet19 network is used to replace the original darknet53 network of YOLOv3. During feature learning, the target frames in MOA images are analyzed by K-means clustering algorithm according to different MOA length width ratios in samples. The anchor frames in the center of samples are re clustered to get the appropriate number and size of anchor frames. Finally, the improved YOLOv3 model is used to complete the MOA infrared image fault detection. The experimental results show that the recognition accuracy of the improved model reaches 96.3%, and the recognition speed is 6.75ms.
  • 在工业生产中,及时掌握设备的工作状态是十分必要的。好多情况下,物体表面温度可以间接地反映设备的工作状态,因此开展依据物体表面温度对工业设备运行状态进行检测和监测的理论与方法研究具有重要的意义[1-2]

    对于设备表面温度分布的测量,最常用的设备为红外热像仪。通过红外热像仪进行测温,一般要求正对测量表面。然而在狭小空间或由其他原因造成可达性较差的位置进行测量时,如舰船的机舱,常常无法满足红外热成像测温的要求,因此就需要人们对特殊测量条件下所测得的温度分布进行修正。造成红外热像测温误差的因素很多,相关的研究包括红外热像仪的测温的影响因素分析及误差的计算[1, 3]、太阳、大气等环境的影响[1, 4-5],周围高温物体的影响[6-7]、测量距离及视角的影响[8-9]、探测器温度的影响[10]等。对于本文讨论的非平设备表面,尤其在热像仪处于被测点的天顶角的50°范围之外时,发射率将发生明显的变化。因此,对于设为恒定表面发射率的热像仪,必将会产生显著的测温误差。因此,开展非平表面红外温度测量的修正研究具有重要的应用价值。

    本文针对利用红外热像仪大角度测温或对非平表面进行测温因无法或部分无法处于表面各点定向发射率恒定的天顶角范围内时热像仪的测温误差进行研究,并基于物体表面定向发射率的变化规律,通过对所测表面进行3D建模,建立实际表面和热像图之间的坐标对应关系,提出热像仪所测温度分布的修正方法。方法通过实验进行验证。

    根据红外热像仪温度计算公式[1]

    $$ {T_0} = {\left\{ {\frac{1}{\varepsilon }\left[ {\frac{1}{{{\tau _{\rm{a}}}}}T_{\rm{r}}^n - \left( {1 - a} \right)T_{\rm{u}}^n - \frac{{{\varepsilon _{\rm{a}}}}}{{{\tau _{\rm{a}}}}}T_{\rm{a}}^n} \right]} \right\}^{\frac{1}{n}}} $$ (1)

    式中:ε为物体表面发射率; τa为大气的光谱透射率; Tr为热像仪指示的辐射温度; εa为大气辐射率; Tu为环境温度。在3~5μm,n=9.2554;在8~12μm,n=3.9889。

    由式(1)可知,温度测量值与被测物体的发射率ε密切相关。在测量复杂物体表面时,需在热像仪上将发射率设为一固定值,通常为该表面材料的法向发射率。天顶角是某个方向和天顶方向的夹角,通常记为θ(见图 1)。本文中利用天顶角来描述热像仪测量物体表面中点的温度时,热像仪成像中心点相对于物体表面的位置关系。在一次测量时,热像仪所处的位置,对于非平表面上不同的点具有不同的天顶角θ。当角θ超出一定范围(如大于50°)时,定向发射率的变化必然对红外测温的准确性产生影响。

    图  1  热像仪机位与其对应的天顶角
    Figure  1.  Position of thermal imager and corresponding zenith angle

    根据菲涅尔定律[1],光谱方向反射率表示为:

    $$ \rho (\tilde n, \theta ) = \frac{1}{2}({\left| {\frac{{\tilde n\cos \theta - \cos \theta '}}{{\tilde n\cos \theta + \cos \theta '}}} \right|^2} + {\left| {\frac{{\cos \theta - \tilde n\cos \theta '}}{{\cos \theta + \tilde n\cos \theta '}}} \right|^2}) $$ (2)

    式中:$\tilde n = n\left( \lambda \right) - i \cdot k\left( \lambda \right)$表示材料的复折射系数。

    根据Snell定律,折射角与入射角满足:

    $$ \sin \theta ' = (1/\tilde n)\theta $$ (3)

    根据Kirchhoff定律:

    $$ \varepsilon (\tilde n, \theta ) = \alpha (\tilde n, \theta ) = 1 - \rho (\tilde n, \theta ) $$ (4)

    表 1[11]中铝(Al)、黄铜(Brass)、聚乙烯(PS)、聚苯乙烯(PE)4种材料分别在3~5 μm与8~12 μm工作波段的复折射率可计算出其定向发射率随天顶角变化的理论值。

    表  1  材料各波段复折射率
    Table  1.  Complex refractive indexes of the materials in each wave bands
    Al Brass PS PE
    band/μm 3-5 8-13 3-5 8-13 3-13 3-13
    n 6.62 23.19 4.77 17.21 1.57 1.52
    k 37.77 85.04 23.73 51.85 0 0
    下载: 导出CSV 
    | 显示表格

    式(4)的计算结果如图 2所示,可观察到金属与非导电体定向发射率随天顶角的变化规律:各类材料定向发射率在被测点天顶角大于50°附近开始出现明显的变化。

    图  2  定向发射率随天顶角的变化
    Figure  2.  Variation of direction emissivity with zenith angle

    在工程测量中,需在测量前设定被测物体表面的发射率,通常设为其法向发射率ε。对于相对平整的表面,可使用ε代表所有点的发射率。当热像仪处于被测点的天顶角大于50°的位置范围时,需采用该角度的定向发射率。对于不同角度的定向发射率εθ,根据(1)式可知实际温度T0与测量温度T1的关系为:

    $$ {T_{\rm{1}}} = {\left\{ {\frac{{{\varepsilon _ \bot }}}{{{\varepsilon _\theta }}}} \right\}^{\frac{1}{n}}}{T_0} $$ (5)

    定义修正系数为:

    $$ c\left( \theta \right) = {\left\{ {\frac{{{\varepsilon _ \bot }}}{{{\varepsilon _\theta }}}} \right\}^{\frac{1}{n}}} $$ (6)

    则可得修正系数曲线如图 3所示。

    图  3  温度修正系数随天顶角的变化
    Figure  3.  Temperature correction factors vs. zenith angles

    图 3观察得出,对于金属表面温度测量结果的修正系数在天顶角大于50°时减小,在90°附近上升迅速,而非导电体在天顶角大于70°时迅速增加,最大时测温误差可达50%。因此,需要对红外热像图中不同位置的温度值依据其对应实际点的修正系数进行修正。

    如前所述,要计算被测表面各个点的修正系数,必须首先计算各个点红外测温时热像仪所处的天顶角θ

    利用3D建模或3D扫描技术可获取被测物体的空间结构,存储格式为.stl文件。将所得.stl文件使用Meshlab软件进行网格划分处理后可获得该被测物体的点云数据,存储为.xyz文件。.xyz文件为一组六维向量,前三维是点云中各点世界坐标系下坐标${p_{\rm{w}}} = {[{x_{\rm{w}}}, {y_{\rm{w}}}, {z_{\rm{w}}}]^{\rm{T}}}$,后三维是该点的法向量坐标$\mathop N\limits^ \to = {\left[ {{n_x}, {n_y}, {n_z}} \right]^{\rm{T}}}$。

    以一锅炉为例,绘制出其点云坐标和法向量如图 4所示。

    图  4  某锅炉点云坐标与法向量示例
    Figure  4.  Point cloud coordinates and normal vectors

    本文将建模过程中的模型原点作为世界坐标系原点。可知热像仪机位中心坐标\$。

    由热像仪机位中心坐标与空间坐标可得测量方向向量:$\overrightarrow {{p_{\rm{w}}}I} = {\left[ {{x_{\rm{c}}} - {x_{\rm{w}}}, {y_{\rm{c}}} - {y_{\rm{w}}}, {z_{\rm{c}}} - {z_{\rm{w}}}} \right]^{\rm{T}}}$

    则可计算天顶角为:

    $$ \theta ({p_{\rm{w}}}, I) = \langle \mathop N\limits^ \to , \overrightarrow {{p_{\rm{w}}}I} \rangle = \arccos (\frac{{\mathop N\limits^ \to \cdot \overrightarrow {{p_{\rm{w}}}I} }}{{\left| {\mathop N\limits^ \to } \right| \cdot \left| {\overrightarrow {{p_{\rm{w}}}I} } \right|}}) $$ (7)

    由空间几何关系可知,对于天顶角小于90°的被测点,由于透视关系无法在热像仪上成像,因此可得可视部分如图 5所示。

    图  5  锅炉的红外热像仪可视点
    Figure  5.  Visible points of boiler for the infrared imager

    其中坐标为(-0.75, 0.12, 0.1)的点为红外热像仪中心所在点,黄色部分为锅炉表面在该位置关系下热像仪的可视点,蓝色为不可视点。

    由该位置关系可计算得出各点对应天顶角θ。由于实际工程应用中,通常在金属表面填涂一层非导电涂料,且非导电体在法向发射率相较于导电体更大,理论的光谱定向发射率与实际相差不大,因此以8~13 μm的PS(Poly styrene)材料作为被测物体材料。在图 4的位置关系下,计算得出该锅炉的修正系数c的频率分布直方图如图 6所示。经计算可得,该位置关系下修正系数大于1.2的点占比44.1%。由此可见,对实际非平表面测温时的定向发射率进行修正是非常必要的。

    图  6  不同温度修正系数对应的被测点所占总被测点的比例
    Figure  6.  Percentage of points with different temperature correction coefficients

    天顶角与成像点的实际位置和其法向量密切相关。因此,针对红外热像仪所测温度的修正需明确各成像点在世界坐标中的实际位置。

    世界坐标系下被测点pw经红外热像仪投影至平面像素点pp=[u, v]T,在成像过程中经过两次坐标变换,分别是由世界坐标系到热像仪坐标系的变换和由热像仪坐标系到像素坐标系的变换[12],两个变换分别通过外参矩阵和内参矩阵来实现。

    世界坐标pw到热像仪坐标pi的转换关系为:

    $$ \left[ {\begin{array}{*{20}{c}} {{p_{\rm{c}}}}\\ 1 \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} R&t\\ {\rm{0}}&{\rm{1}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{p_{\rm{w}}}}\\ 1 \end{array}} \right] $$ (8)

    其中:

    $$ \left[ {\begin{array}{*{20}{c}} R&t\\ {\rm{0}}&{\rm{1}} \end{array}} \right] = {\left[ {\begin{array}{*{20}{c}} {{R_{\rm{i}}}}&I\\ {\rm{0}}&{\rm{1}} \end{array}} \right]^{ - 1}} = \left[ {\begin{array}{*{20}{c}} {{R_{\rm{i}}}^{\rm{T}}}&{ - {R_{\rm{i}}}^{\rm{T}}I}\\ 0&1 \end{array}} \right] $$ (9)

    为热像仪的外参矩阵,Ri为热像仪的基坐标在世界坐标系的坐标,描述热像仪在世界坐标系中的姿态,I为热像仪中心在世界坐标系中的位置。

    获取热像仪的外参矩阵需确定热像仪姿态,热像仪在成像过程中一般固定角度。假设热像仪在实际测量中热像仪x轴始终与水平面平行,仅存在俯角γ与转角β。则可知当热像仪对准被测物体中心$\overline {{p_{\rm{w}}}} = {[\overline {{x_{\rm{w}}}} , \overline {{y_{\rm{w}}}} , \overline {{z_{\rm{w}}}} ]^{\rm{T}}}$时,热像仪中心指向被测点中心向量$\overrightarrow {C\overline {{p_{\rm{w}}}} } = {[\overline {{x_{\rm{w}}}} - {x_{\rm{c}}}, \overline {{y_{\rm{w}}}} - {y_{\rm{c}}}, \overline {{z_{\rm{w}}}} - {z_{\rm{c}}}]^{\rm{T}}}$为热像仪坐标z轴方向,热像仪坐标系x轴与水平面平行,y轴方向垂直于x轴和z轴。由数学关系可得3个方向的标准基为:

    $$ \begin{array}{l} {e_{\rm{1}}} = {\left[ {\overline {{y_{\rm{w}}}} - {y_{\rm{c}}}, {x_{\rm{c}}} - \overline {{x_{\rm{w}}}} , 0} \right]^{\rm{T}}}/{\left( {\overline {{x_{\rm{w}}}} - {x_{\rm{c}}}} \right)^2} + {\left( {\overline {{y_{\rm{w}}}} - {y_{\rm{c}}}} \right)^2}\\ {e_{\rm{2}}} = {e_3} \times {e_1}/\left| {{e_1}} \right|\left| {{e_3}} \right|\\ {e_{\rm{3}}} = {\left[ {\overline {{x_{\rm{w}}}} - {x_{\rm{c}}}, \overline {{y_{\rm{w}}}} - {y_{\rm{c}}}, \overline {{z_{\rm{w}}}} - {z_{\rm{c}}}} \right]^{\rm{T}}}/\\ {\rm{ }}{\left( {\overline {{x_{\rm{w}}}} - {x_{\rm{c}}}} \right)^2} + {\left( {\overline {{y_{\rm{w}}}} - {y_{\rm{c}}}} \right)^2} + {\left( {\overline {{z_{\rm{w}}}} - {z_{\rm{c}}}} \right)^2} \end{array} $$ (10)

    计算可得Ri=[e1, e2, e3]为热像仪在世界坐标系旋转到当前姿态的旋转矩阵。

    当热像仪角度变化时,可根据热像仪的俯角与转角求得热像仪坐标系z轴与水平面的交点:

    $$ {p_0} = {\left[ {\left( {{x_{\rm{c}}} + {z_{\rm{c}}}/\tan \gamma } \right)\sin \beta , \left( {{y_{\rm{c}}} + {z_{\rm{c}}}/\tan \gamma } \right)\cos \beta , 0} \right]^{\rm{T}}} $$ (11)

    将上述计算过程中被测点中心$\overline {{p_{\rm{w}}}} $替换为p0即可根据热像仪实际姿态确定外参矩阵。

    在世界坐标转换为热像仪坐标后,需将热像仪坐标转换为像素坐标,完成三维到二维的转变。

    热像仪坐标pi到像素坐标pp的转换关系为:

    $$ \left[ {\begin{array}{*{20}{c}} u\\ v\\ 1 \end{array}} \right] = \frac{1}{{{z_{{p_{\rm{c}}}}}}}\left[ {\begin{array}{*{20}{c}} {{k_x}}&0&{{c_x}}\\ 0&{{k_y}}&{{c_y}}\\ 0&0&1 \end{array}} \right]{p_{\rm{c}}} $$ (12)

    式中:$\mathit{\boldsymbol{K}} = \left[ {\begin{array}{*{20}{c}} {{k_x}}&0&{{c_x}}\\ 0&{{k_y}}&{{c_y}}\\ 0&0&1 \end{array}} \right]$为热像仪的内参矩阵,$\frac{1}{{{z_{{p_{\rm{c}}}}}}}$为尺度因子,描述成像位置与热像仪中心的距离。

    内参矩阵与热像仪的型号和图像的分辨率有关,kxky为热像仪xy方向的像素焦距,cxcy为图片主点偏移,内参矩阵均以像素为单位。

    fxfy为热像仪xy方向的毫米焦距,对于毫米焦距和像素焦距的转换关系有:

    $$ \begin{array}{l} {f_x}{\rm{ }} = {\rm{ }}{k_{\rm{u}}}{\rm{ }}*{\rm{ d}}x\\ {f_y}{\rm{ }} = {\rm{ }}{k_{\rm{v}}}{\rm{ }}*{\rm{ d}}y\\ \;{\rm{d}}x{\rm{ }} = {\rm{ (}}1/{\rm{ccd}}\_{\rm{size }}*{\rm{ }}24.5{\rm{ }}*\\ \quad \;{\rm{image}}\_x{\rm{ }}/{\rm{ image}}\_y{\rm{) }}/{\rm{ image}}\_x\\ \;{\rm{d}}y{\rm{ }} = {\rm{ (}}1/{\rm{ccd\_size }}*{\rm{ }}24.5{\rm{ }}*\\ \quad \;{\rm{image}}\_x{\rm{ }}/{\rm{ image\_}}y{\rm{) }}/{\rm{ image}}\_y \end{array} $$ (13)

    式中:ccd_size是热像仪厂商给出的CCD尺寸,以英寸为单位; image_x、image_y是图像xy方向的像素数,与成像的分辨率有关。

    在工程应用中,外参矩阵主要通过机位的计算得出,内参矩阵经热像仪标定得出。

    对于热像仪标定,一般使用张氏标定法[13]。通过Matlab软件自带的热像仪标定程序camera calibrator可对单目热像仪进行标定。可见光相机标定需打印棋盘格,对于红外热像仪,将棋盘格用金属板代替。黑色部分使用隔热胶带,白色部分使用金属板或锡纸。

    在程序中输入标定板在不同角度下的红外热像图,可得到热像仪中心与标定板相对位置和热像仪内参、外参、图片尺寸等参数。本文主要使用标定得到的热像仪内参矩阵K

    总体的温度修正计算流程如图 7所示。

    图  7  温度修正流程图
    Figure  7.  Flow chart of temperature compensation

    据式(7)计算出天顶角,查得该天顶角下的发射率,然后计算式(6)得到世界坐标系下各测温点的修正系数c(pw)。经热像仪成像后,测温点由世界坐标pw映射到成像平面坐标pp=[u, v]Tpwpp的转换关系由3.3节标定的热像仪矩阵计算得出。经上述转换关系可得测温点在二维平面坐标对应的修正系数c(pp)。

    热像仪拍摄的温度图像将温度以矩阵的形式存储,因此,应当将pp从离散的点云转换为与热像图格式相同的pi, j,其中,pi, j代表像素块(i, j)的坐标。

    对于修正系数在像素块范围内梯度较小的点,该像素块的修正系数取像素块内的平均值。对在热像图上像素位置为(i, j)的点,该像素的修正系数为:

    $$ {c_{i, j}}\left( {{p_{\rm{p}}}} \right) = \frac{{\sum\limits_{k = 1}^m {{c_0}\left( {u, v} \right)} }}{m}, \{ u \in [i, i + 1), v \in [j, j + 1)\} $$ (14)

    当像素内修正系数梯度变化较大的,可采用该像素内修正系数最大值附近点取平均值的方法,在此像素块内凸显该形状非光滑点。

    $$ {c_{i, j}}\left( {{p_p}} \right) = \max \left( {{c_0}\left( {u, v} \right)} \right), \{ u \in [i, i + 1), v \in [j, j + 1)\} $$

    对于红外热像仪,其所测温度矩阵为T0(pp)。修正后的温度矩阵为所测温度矩阵与温度修正系数矩阵做基本积(哈达玛积):

    $$ {T_{i, j}}\left( {{p_{\rm{p}}}} \right) = {T_{0i, j}}\left( {{p_{\rm{p}}}} \right) * {c_{i, j}}\left( {{p_{\rm{p}}}} \right) $$ (15)

    所得温度T(pp)即修正后的温度。

    本文通过一个热水桶表面的温度测量实验来验证该算法有效性。首先使用Solidworks软件进行几何建模,获得被测物体几何信息,经Meshlab软件求得热像仪可视表面的点云及各点处的法向量,再经Matlab计算出各点温度的修正系数以及各点对应的热像图中的坐标,最后将红外热像仪直接测得的温度分布信息导入Matlab进行修正得到水桶表面的修正后的温度。

    本文选用Fluke Ti400型号热像仪进行测温,测温表面为一聚乙烯塑料桶侧面,塑料桶底部外径为18 cm,顶部外径为20.5 cm,高为20 cm。塑料桶中为热水,热像仪在世界坐标系中的坐标为$C{\rm{ = [}} - {\rm{0}}{\rm{.75, 0}}{\rm{.12, 0}}{\rm{.10}}{{\rm{]}}^{\rm{T}}}$,热像仪正对塑料桶中心,位置关系如图 8所示。测得塑料桶外表面可视部分的红外热像图见图 9

    图  8  热像仪机位
    Figure  8.  Position of thermal imager
    图  9  红外热像图与温度截线
    Figure  9.  Infrared thermal image and temperature crossing line

    使用热像仪拍摄红外标定板(如图 10所示),获取标定板的红外热像图,将热像图输入标定程序,获得热像仪内参矩阵K

    图  10  红外热像仪的标定
    Figure  10.  Calibration of infrared thermal imager

    通过Solidworks软件建模并将得到的水桶表面的云数据导入Matlab软件中,可得到水桶表面与热像仪在空间中的位置关系,如图 11所示(黄色点为被测表面的可视点,蓝色为被测表面的非可视点)。本实验中,略去了接近桶口的外径较大的部分。

    图  11  热像仪与被测表面的位置关系
    Figure  11.  Positions of thermal imager and measured surface

    依据热像仪与可视点的空间位置关系,据式(7)可计算得各可视点的天顶角,据图 2可得不同天顶角时,所对应的发射率。然后据式(6)算得每个可视点的修正系数。如4.1节所述,对空间中可视点进行热像仪矩阵的变换,映射至平面可得出热像图中各点对应的修正系数,见图 12。进而依据式(15)得到热像图中各点修正后的温度分布。

    图  12  修正系数分布
    Figure  12.  Correct factors distribution

    为了更好地显示本文方法对红外热像仪所测得的温度分布,本节在所研究的可视表面高度方向的中心,取与图 9中相同位置的截线进行展示,修正系数如图 13所示。将所得修正系数与成像截线温度相乘得到修正后的截线温度,修正前后温度的对比如图 14所示。

    图  13  截线上各点的修正系数
    Figure  13.  Correction factors of points in cross-line
    图  14  截线上温度分布修正前和修正后的对比
    Figure  14.  Comparison of temperature distributions before and after correction of points in the cross-line

    图 14可见,修正前最大误差近40℃,真实温度为335℃,相对误差约为12.4%;修正后,集中在两侧大角度测量时的误差不超过5.0℃,相对误差小于1.5%。从此可见该方法的有效性。值得说明的是图 13中修正系数分布曲线并不完全对称,主要原因是该图的横坐标为非常微小的像素点,研究过程中像素点的选取没有完全对称所致。更多的平行截线结果均表明,该修正方法较好地消除了非平表面在红外热像测温时部分因为测量角度过大而引起的发射率变化所造成的误差。

    针对使用红外热像仪测量非平表面温度时不同被测量点的定向发射率的变化对测温结果的影响,本文依据物体表面定量发射率的变化规律,给出了各个测量点的修正系数。同时,通过点云三维建模,利用热像仪通用几何成像原理推导出红外热像图像素与实际物体的被测表面之间各个测量点的对应关系,给出了红外热像测温过程中非平表面上的大天顶角被测点的温度测量结果的误差修正方法。实验证明了该方法的有效性。

  • 图  1   不同状态MOA红外图像

    Figure  1.   MOA infrared images in different states

    图  2   MOA红外故障检测模型框架图

    Figure  2.   Framework diagram of MOA infrared fault detection model

    图  3   中值滤波处理后图像灰度值结果

    Figure  3.   Grayscale value result of the image after median filtering processing

    图  4   K-means聚类分析结果

    Figure  4.   K-means cluster analysis result YOLO recognition

    图  5   修改后的YOLOv3网络结构图

    Figure  5.   Modified network structure of YOLOv3

    图  6   YOLOv3-Darknet53训练损失变化

    Figure  6.   Loss trend of YOLOv3-Darknet53 training

    图  7   YOLOv3-Darknet19训练损失变化

    Figure  7.   Loss trend of YOLOv3-Darknet19 training

    图  8   不同网络训练损失对比

    Figure  8.   Comparison of different network training loss

    图  9   MOA红外图像检测结果

    Figure  9.   MOA Infrared image detection results

    表  1   Darknet19网络结构

    Table  1   Network structure of Darknet19

    Type Fitlers Size/stride Output
    Convolutional 32 3×3/1 256×256
    Maxpool 2×2/2 128×128
    Convolutional 64 3×3/1 128×128
    Maxpool 2×2/2 64×64
    Convolutional 128 3×3/1 64×64
    Convolutional 64 1×1/1 64×64
    Convolutional 128 3×3/1 64×64
    Maxpool 2×2/2 32×32
    Convolutional 256 3×3/1 32×32
    Convolutional 128 1×1/1 32×32
    Convolutional 256 3×3/1 32×32
    Maxpool 2×2/2 16×16
    Convolutional 512 3×3/1 16×16
    Convolutional 256 1×1/1 16×16
    Convolutional 512 3×3/1 16×16
    Convolutional 256 1×1/1 16×16
    Convolutional 512 3×3/1 16×16
    Maxpool 2×2/2 8×8
    Convolutional 1024 3×3/1 8×8
    Convolutional 512 1×1/1 8×8
    Convolutional 1024 3×3/1 8×8
    Convolutional 512 1×1/1 8×8
    Convolutional 1024 3×3/1 8×8
    下载: 导出CSV

    表  2   YOLOv3的先验框尺寸

    Table  2   Prior box size of YOLOv3

    Scale prediction Feature map size Anchor value
    Scale 1 13×13 (105, 153) (49, 250)(32, 132)
    Scale 2 26×26 (144, 43) (33, 78)(47, 56)
    Scale 3 52×52 (52, 25) (22, 34)(16, 15)
    下载: 导出CSV

    表  3   YOLOv3-Darknet19训练环境

    Table  3   Training environment of YOLOv3-Darknet19

    Term configuration
    CPU Core i9 10900K
    GPU INVIDIA RTX 3080
    Operating system Microsoft Windows10
    Parallel computing library Cuda10.0+Cudnn7.4.1.5
    Image processing Python3.7
    Deep learning framework Darknet19
    下载: 导出CSV
  • [1] 任大江, 叶海鹏, 李建萍, 等. 一起500 kV金属氧化锌避雷器故障原因分析[J]. 电瓷避雷器, 2020(3): 127-132. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ202302023.htm

    REN Dajiang, YE Haipeng, LI Jianping, et al. Analysis of the causes of a 500 kV metal zinc oxide arrester fault[J]. Insulators and Surge Arresters, 2020(3): 127-132. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ202302023.htm

    [2] 张明龙, 钱健, 王健, 等. 基于配电避雷器全电流监测的系统过电压告警研究[J]. 高压电器, 2020, 56(8): 256-260, 267. DOI: 10.13296/j.1001-1609.hva.2020.08.039

    ZHANG Minglong, QIAN Jian, WANG Jian, et al. Research on over-voltage alarm based on full current monitoring of distribution arrester[J]. High Voltage Apparatus, 2020, 56(8): 256-260, 267. DOI: 10.13296/j.1001-1609.hva.2020.08.039

    [3] 何贵先, 行鸿彦, 季鑫源, 等. 金属氧化物避雷器在线监测的谐波校正及研究[J]. 电子测量与仪器学报, 2017, 31(10): 1549-1554. https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201710004.htm

    HE Guixian, XING Hongyan, JI Xinyuan, et al. Harmonic correction and study for MOA on-line monitoring algorithm[J]. Journal of Electronic Measurement and Instrumentation, 2020(3): 127-132. https://www.cnki.com.cn/Article/CJFDTOTAL-DZIY201710004.htm

    [4] 陈登义, 周雪会, 李鹏. 电涌保护器老化监测指标对比研究[J]. 电瓷避雷器, 2016(5): 90-93. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201605017.htm

    CHEN Dengyi, ZHOU Xuehui, LI Peng, et al. Comparative study about index of SPD aging monitoring[J]. High Voltage Apparatus, 2016(5): 90-93. https://www.cnki.com.cn/Article/CJFDTOTAL-DCPQ201605017.htm

    [5] 陈丹, 傅中君, 柳益君, 等. 一种氧化锌避雷器阻性电流的提取方法[J]. 电测与仪表, 2019, 56(13): 117-122. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ201913021.htm

    CHEN Dan, FU Zhongjun, LIU Yijun, et al. A method for the resistive current extraction of metal oxide surge arresters[J]. Electrical Measurement & Instrumentation, 2019, 56(13): 117-122. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ201913021.htm

    [6] 李艳鹏, 晋涛, 梁基重, 等. 红外技术在避雷器绝缘缺陷检测中的应用[J]. 电测与仪表, 2019, 56(20): 87-90. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ201920015.htm

    LI Yanpeng, JIN Tao, LIANG Jichong, et al. The application of infrared diagnosis in detection of lightning arrester insulation defects[J]. Electrical Measurement & Instrumentation, 2019, 56(20): 87-90. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ201920015.htm

    [7] 李红光, 于若男, 丁文锐. 基于深度学习的小目标检测研究进展[J/OL]. 航空学报: 1-19, [2020-12-19]. http://kns.cnki.net/kcms/detail/11.1929.V.20201026.0947.004.html.

    LI Hongguang, YU Ruonan, DING Wenrui. Research development of small object tracking based on deep learning[J/OL]. Acta Aeronautica et Astronautica Sinica: 1-19, [2020-12-19]. http://kns.cnki.net/kcms/detail/11.1929.V.20201026.0947.004.html.

    [8] 黄继鹏, 史颖欢, 高阳. 面向小目标的多尺度Faster-RCNN检测算法[J]. 计算机研究与发展, 2019, 56(2): 319-327. https://www.cnki.com.cn/Article/CJFDTOTAL-JFYZ201902008.htm

    HUANG Jipeng, SHI Yinhuan, GAO Yang. Multi-scale faster-RCNN algorithm for small object detection[J]. Journal of Computer Research and Development, 2019, 56(2): 319-327. https://www.cnki.com.cn/Article/CJFDTOTAL-JFYZ201902008.htm

    [9] 赵振兵, 崔雅萍, 戚银城, 等. 基于改进的R-FCN航拍巡线图像中的绝缘子检测方法[J]. 计算机科学, 2019, 46(3): 159-163. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJA201903024.htm

    ZHAO Zhenbing, CUI Yaping, QI Yincheng, et al. Detection method of insulator in aerial inspection image based on modified R-FCN[J]. Computer Science, 2019, 46(3): 159-163. https://www.cnki.com.cn/Article/CJFDTOTAL-JSJA201903024.htm

    [10] 李文璞, 谢可, 廖逍, 等. 基于Faster RCNN变电设备红外图像缺陷识别方法[J]. 南方电网技术, 2019, 13(12): 79-84. https://www.cnki.com.cn/Article/CJFDTOTAL-NFDW202110012.htm

    LI Wenpu, XIE Ke, LIAO Yao, et al. Intelligent diagnosis method of infrared image for transformer equipment based on improved faster RCNN[J]. Southern Power System Technology, 2019, 13(12): 79-84. https://www.cnki.com.cn/Article/CJFDTOTAL-NFDW202110012.htm

    [11] 刘云鹏, 裴少通, 武建华, 等. 基于深度学习的输变电设备异常发热点红外图片目标检测方法[J]. 南方电网技术, 2019, 13(2): 27-33. https://www.cnki.com.cn/Article/CJFDTOTAL-NFDW201902006.htm

    LIU Yunpeng, PEI Shaotong, WU Jianhua, et al. Deep learning based target detection method for abnormal hot spots infrared images of transmission and transformation equipment[J]. Southern Power System Technology, 2019, 13(2): 27-33. https://www.cnki.com.cn/Article/CJFDTOTAL-NFDW201902006.htm

    [12] 赵琰, 刘荻, 赵凌君. 基于YOLOv3的复杂环境红外弱小目标检测[J]. 航空兵器, 2019, 26(6): 29-34. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906006.htm

    ZHAO Yan, LIU Di, ZHAO Lingjun. Infrared dim and small target detection based on YOLOv3 in complex environment[J]. Aero Weaponry, 2019, 26(6): 29-34. https://www.cnki.com.cn/Article/CJFDTOTAL-HKBQ201906006.htm

    [13] 罗元, 王薄宇, 陈旭. 基于深度学习的目标检测技术的研究综述[J]. 半导体光电, 2020, 41(1): 1-10. https://www.cnki.com.cn/Article/CJFDTOTAL-BDTG202001001.htm

    LUO Yuan, WANG Boyu, CHEN Xu. Research progress of target detection technology based on deep learning[J]. Semiconductor Optoelectronics, 2020, 41(1): 1-10. https://www.cnki.com.cn/Article/CJFDTOTAL-BDTG202001001.htm

    [14] 王永平, 张红民, 彭闯, 等. 基于YOLO v3的高压开关设备异常发热点目标检测方法[J]. 红外技术, 2020, 42(10): 983-987. http://hwjs.nvir.cn/article/id/hwjs202010011

    WANG Yongping, ZHANG Hongmin, PENG Chuang, et al. The target detection method for abnormal heating point of high-voltage switchgear based on YOLO v3[J]. Infrared Technology, 2020, 42(10): 983-987. http://hwjs.nvir.cn/article/id/hwjs202010011

    [15] 王芳, 李传强, 伍博, 等. 基于多尺度特征融合的红外小目标检测方法[J]. 红外技术, 2021, 43(7): 688-695. http://hwjs.nvir.cn/article/id/483e7824-ea0f-4381-96b7-22257615d475

    WANG Fang, LI Chuanqiang, WU Bo, et al. Infrared small target detection method based on multi-scale feature fusion[J]. Infrared Technology, 2021, 43(7): 688-695. http://hwjs.nvir.cn/article/id/483e7824-ea0f-4381-96b7-22257615d475

    [16] 刘云鹏, 张喆, 裴少通, 等. 基于深度学习的红外图像中劣化绝缘子片的分割方法[J]. 电测与仪表, 2022, 59(9): 63-68. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ202209009.htm

    LIU Yunpeng, ZHANG Zhe, PEI Shaotong, et al. Faulty insulator segmentation method in infrared image based on deep learning[J]. Electrical Measurement & Instrumentation, 2022, 59(9): 63-68. https://www.cnki.com.cn/Article/CJFDTOTAL-DCYQ202209009.htm

    [17] 中华人民共和国国家发展和改革委员会. 带电设备红外诊断应用规范[S]. DL/T 664-2008, [2008-11-01].

    National Development and Reform Commission. Application Rules of infrared diagnosis for live electrical equipment[S]. DL/T 664-2008, [2008-11-01].

    [18] 汪烈兵, 姜雄飞, 石春光, 等. 基于图像滤波与Hough变换的红外弱小目标检测[J]. 红外技术, 2020, 42(7): 683-687. http://hwjs.nvir.cn/article/id/hwjs202007012

    WANG Liebing, JIANG Xiongfei, SHI Chunguang, et al. Infrared small target detection based on image filtering and hough transform[J]. Infrared Technology, 2020, 42(7): 683-687. http://hwjs.nvir.cn/article/id/hwjs202007012

    [19]

    Redmon J, Farhadi A. YOLOv3: An incremental improvement [J/OL]//Computer Science, 2018 https://arxiv.org/abs/1804.02767.

    [20]

    CHENG G, LIU L. Survey of image segmentation methods based on clustering[C]//IEEE International Conference on Information Technology (ICIBA), 2020: 1111-1115.

  • 期刊类型引用(2)

    1. 宋宏伟,刘塑,李海英,于洪有,史圣兵,郭秀丽,李昊岚,张亚平,王向前. 中波红外640×512(25μm)探测器组件可靠性验证试验. 红外技术. 2024(08): 879-882 . 本站查看
    2. 何邦盛,王忠华. 基于结构张量筛选和局部对比度分析的空中红外小目标检测算法. 红外技术. 2023(12): 1169-1176 . 本站查看

    其他类型引用(1)

图(9)  /  表(3)
计量
  • 文章访问数:  150
  • HTML全文浏览量:  44
  • PDF下载量:  36
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-06-06
  • 修回日期:  2021-08-19
  • 刊出日期:  2023-11-19

目录

/

返回文章
返回
x 关闭 永久关闭

尊敬的专家、作者、读者:

端午节期间因系统维护,《红外技术》网站(hwjs.nvir.cn)将于2024年6月7日20:00-6月10日关闭。关闭期间,您将暂时无法访问《红外技术》网站和登录投审稿系统,给您带来不便敬请谅解!

预计6月11日正常恢复《红外技术》网站及投审稿系统的服务。您如有任何问题,可发送邮件至编辑部邮箱(irtek@china.com)与我们联系。

感谢您对本刊的支持!

《红外技术》编辑部

2024年6月6日