基于三维重建的输电线路覆冰在线监测方法_杨浩

DOI:10.3969/j
.issn.1000-1026.2012.23.019基于三维重建的输电线路覆冰在线监测方法
杨 浩1,吴 畏1,
(1.输配电装备及系统安全与新技术国家重点实验室,重庆大学,重庆市400044;
2.四川省电力公司内江电业局,四川省内江市641100
)摘要:输电线路覆冰常危害电力系统安全运行,有效监测输电线路覆冰状况能够预防由覆冰引起
的重大电力事故。文中提出一种基于图像三维重建的输电线路覆冰在线监测方法,该方法利用计算机双目视觉原理,通过2台不同位置的摄像机采集同一段线路的2幅覆冰图像,检测出图像的匹配点,利用图像间的视差计算线路覆冰在三维空间的深度,重建出线路覆冰的三维模型,最后由三维模型计算得到线路覆冰的厚度及其分布状况。应用所提出方法对2组实测图像进行三维重建分析,结果表明该方法不仅能够准确计算线路覆冰的厚度,还可反映覆冰的厚度分布特征。关键词:输电线路;覆冰;图像处理;三维重建;在线监测
收稿日期:2012-01-15;修回日期:2012-04-
05。0 引言
输电线路覆冰会引发严重危害电力系统安全运
行的事故,如倒杆(塔)、绝缘子串闪络等,造成巨大
的经济损失[
1-
3]。目前针对输电线路覆冰监测主要采用绝缘子倾角[4]、导线拉力[5]
、基于空气与冰的电阻和介电常数差异检测[6-7]及图像视频监测[8-
10]等方式。其中,前两者的覆冰厚度计算模型复杂,且覆
冰导线的重力变化、绝缘子串倾斜角和风偏角等参
数非同步采集,
计算误差较大;由于冰和空气电传导性的复杂性及测量干扰,基于电阻和介电常数差异
的检测方法不易判别冰层厚度变点[11]
,对不规则覆冰形状的测量精度还有待研究。采用图像视频监测
的方法能直观地监测导线和绝缘子的覆冰状况。文献[9
]通过对摄像机监测到的输电线路覆冰图像进行边缘检测,提取覆冰边界轮廓来测量覆冰厚度;文献[10]利用图像滤波、自适应阈值变换和基于LOG(Laplacian-Gauss)算子的边缘检测技术对远程采集的覆冰图片进行处理和覆冰厚度识别。但是,现有图像监测方法都是对采集的二维图像进行处理和识别,易受摄像机位置和角度的影响,且基于二维平面的厚度估算在精度上很难达到实际要求。
本文提出一种基于双目立体视觉技术的输电线路覆冰图像监测方法,该方法能够直观地给出覆冰厚度分布状况且测量精度高,同时能够实现远程智能监测,弥补覆冰过程监测手段的不足,对于输电线路覆冰在线监测具有重要意义。
1 输电线路覆冰图像三维重建原理
基于双目立体视觉的方法是利用视差[12-
13]原
理,由2个摄像机从不同角度获取同一段线路的覆冰图像,使覆冰线路上的同一点在2幅图像中的像点呈现不同的视觉位置,再根据摄像机成像原理得到两像点间的几何关系,进而计算出该点的三维坐标。通过获得大量特征点的三维数据重建线路覆冰的三维模型,近似于人类视觉系统的立体感知过程。线路覆冰图像三维重建原理图如图1所示,图像a和b是同一段线路的覆冰图像,由2个参数相同、位置固定的摄像机从不同角度拍摄得到
图1 覆冰图像三维重建原理图
Fig.1 Schematic diagram of icing
imag
e 3Dreconstruction图1中:Bc为摄像机投影中心的连线距离;f
为摄像机焦距,即摄像机到图像平面的距离;O1和
301—第36卷 第23期2012年12月10
日Vol.36 No.23
Dec.10,2012
O2为平面坐标原点,
即摄像机光轴与图像平面的交点。以左摄像机中心为原点建立三维坐标系
Xc-Yc-Zc(
也可用右摄像机中心为原点建立),点P为空间中覆冰线路上的一点,在左摄像机坐标系的
三维坐标为(xc,yc,zc)
,在2幅图像上的成像位置分别为Pa(x1,y1)和Pb(x2,y2)
,通过相似三角形的几何关系可得:
xc=
Bcx1
|Pa-Pb|yc=Bcy|Pa-Pb|
zc=
Bcf|Pa-Pb烅烄烆
|(1
)式中:|Pa-Pb|为视差,即点P在2幅图像中的位置偏移。
因摄像机在垂直方向上处于同一平面,故y=
y1=y2,则视差为|x1-x2|。由式(1)可见,覆冰线路上点P的空间坐标可由视差结合摄像机内外参
数计算得出。线路覆冰的三维重建,就是在覆冰图像上出大量表征线路覆冰信息的特征点,
计算特征点的空间坐标。当获得线路覆冰图像特征点的实际空间坐标后,可通过点云或网格方式重建三维模型,由模型计算覆冰厚度分布。线路覆冰图像的三维重建过程分为4个部分:图像获取、摄像机标定、特征点提取和匹配、三维重建。
2 输电线路覆冰图像三维重建方法
2.1 图像获取
采用电荷耦合器件(CCD)摄像机采集现场图像,考虑到现场测量线路覆冰的环境条件,摄像机安
装位置大多是在杆塔或固定支架上,很难垂直线路拍摄图像。试验时将摄像机安装于固定支架上,选择倾斜一定角度拍摄,如图2所示。
图2 覆冰图像获取示意图
Fig.2 Diagram of icing 
image obtain method2台摄像机采用光轴上下平行、基线共线的位置布置,
既容易得到垂直方向上的视差,又可简化后续处理。由于摄像机外部参数受2台摄像机间相对位置影响,试验中保持其相对位置不变,以确保摄像机标定求取参数的准确性。
2.2 摄像机标定
为定量描述摄像机成像过程,定义世界坐标系
Xw-Yw-Zw和摄像机坐标系Xc-Yc-Zc这2个三维
直角坐标系,以及图像物理坐标系X-Y和像素坐标
系U-V这2个平面直角坐标系[8]
如图3所示。由式(1
)可知,覆冰线路上的点的空间坐标与摄像机的内外参数有关,需通过摄像机标定确定摄像机的位置、属性等参数,从而求得覆冰线路上各点的空间坐标。其中,摄像机的内部参数包括其图像中心、焦距、镜头畸变等,外部参数是指摄像机坐标相对于世界坐标系的三维位置和方向。
图3 坐标系示意图
Fig.3 Diagram of coordinate sy
stem2.2.1 坐标系定义及转换
研究中世界坐标系选择支架与线路接点作为原点,线路及其水平、竖直方向垂线作为轴线。摄像机坐标系选择摄像机光心为原点,光轴及其水平、竖直方向垂线作为轴线。摄像机坐标系经原点平移和坐
标轴旋转能与世界坐标系重合,
故只需求取覆冰线路在摄像机坐标系的坐标及其平移旋转参数,就能得到覆冰线路世界坐标,进而计算线路覆冰尺寸。
图3中摄像机模型采用针孔模型,摄像机坐标系沿光轴方向为Zc轴并且垂直于像平面,交点O1为图像物理坐标系的原点。摄像机采集的图像以二维数组的形式储存,数组中的每一个元素即为像素,用像素坐标(u,v)
表示该像素在图像上的列数和行数。设线路上一点P的世界坐标为P(xw,yw,zw)
,摄像机坐标为P(xc,yc,zc)
,图像平面物理坐标为P(x,y)
,2个三维坐标系有如下转换关系:xc
ycz熿燀燄燅c=Rxwywz熿
燄燅w+T(2
)—
401—2012,36(23
) 
xc
yc
zc
熿燀燄
燅1=R T0T
[]
1xwywzw
燀燄燅
led显指数1(3
)式中:R为3×3阶正交旋转矩阵;T=[tXc,tYc
,tZc]T,为三维平移向量,(tXc,tYc,tZc)为世界坐标原点在摄像机坐标系的坐标。
由图3中三角形几何关系,经变换得P点像素
坐标(u,v)
满足(推导见附录A):uv熿燀燄燅
1=1
zcK1K2
xwywzw
熿燀燄燅
1(4
)K1=
fdx0u000
dy
v00 0 燀1 0熿燄燅
0(5
)K2=R T
0T[]
(6
)式中:K1为摄像机内参数矩阵,由摄像机焦距f,单位像素尺寸dx和dy及图像基准点(
物理坐标原点的像素坐标)(u0,v0)
确定;K2为摄像机外参数矩阵,包括旋转矩阵R和平移向量T,取决于摄像机坐标系和世界坐标系的相对位置。
摄像机标定的过程即是求解K1和K2的过程。2.2.2 棋盘格标定软件实现
摄像机标定方法选择张氏标定方法(
见附录B)
[14
],该方法精度较高且具有较好的鲁棒性,其原理是借助一个具有精确定位点阵的标定模板,已知其世界坐标,与单个以上不同方位的图像进行匹配,利用K1和K2矩阵计算出内外参数。成像过程中产生的畸变主要为径向畸变,如图3中P(x,y)为理想成像点,P(xq,yq)
为实际成像点,因此必须对平面坐标进行校正。在张氏标定法中的畸变模型为:
x=xq+xq[k1(x2q+y2q)+k2(x2q+y2q)2
y=yq+yq[k1(x2q+y2q)+k2(x2q+y2q)2
(7
)  u=uq+(uq-u0)[k1(x2q+y2q)+k2(x2q+y2q)2
v=vq+(vq-v0)[k1(x2q+y2
q)
+k2(x2q+y2q)2
烅烄烆
](8
)式中:xq,yq,uq,vq为校正前的坐标;x,y,u,v为校正后的坐标;k1和k2为摄像机的径向畸变系数。
采用最小二乘法,通过标定板图像上点的坐标
求解方程组可得k1和k2。
标定时使用2部相同型号摄像机,均由支架固
定且与线路的相对位置不变,
以保证内部参数和两部相机相对位置不发生变化。选取线路与支架接点
为世界坐标系原点,线路方向为Xw轴,支架和相机所在平面为Yw-Zw平面,
水平方向为Yw轴,垂直方向为Zw轴。试验采用2台维视工业相机MV-130UC,分辨率1 
280×720,像素尺寸5.2μm×5.2μm,其艾菲特定焦镜头焦距25mm,
镜头中心距离62.5mm,2号相机与线路起点在同一水平线,相距300mm。对2部相机分别标定,标定板大小为每格20mm,格数为14×9,将标定板贴于玻璃板上,垂直于线路,中央角点坐标位置为(200,0,0),变换标定板角度拍摄10张标定图像。编程环境采用
VC++,其Op
enCV库提供了张氏标定算法的函数[
15]
,流程如图4所示。图4 标定流程图
Fig
.4 Flow chart of standardization通过以上程序对2部相机进行标定,由内参数矩阵可知有效焦距和基准点像素坐标。旋转向量即
为旋转矩阵的对角线元素,可用cvRodrig
ues2(·)函数进行转换,2台相机内外参数的标定结果分别如表1和表2所示。标定误差分析时,选取标定板外围4个角点计算最大尺寸,
针织牛仔布与实际尺寸进行比较,得到在3个方向上的偏移均小于0.1mm。至此完成摄像机标定,下面将对线路覆冰图像进行处理。
表1 摄像机内参数
Table 1 Inter p
arameters of the camera摄像机f
/mm u0v0k1k2
1号25.008 4 652.451 3 347.924 0-0.152 6 0.107 12号
25.015 9 664.205 3 351.144 2-0.138 7 0.130 
1表2 摄像机外参数
Table 2 Extrinsic p
arameters of the camera摄像机旋转向量
平移向量
1号(0.711 7,0.021 3,0.720 3)(301.116 2,64.517 4,1.778 
1)2号新型环保包装材料
(0.714 3,0.021 3,0.716 9)(301.109 4,2.176 3,1.794 
3)2.3 特征点提取和匹配
在摄像机标定之后,必须对目标图像进行特征点提取和匹配。边缘检测和角点检测都属于特征点提取的方法,由于线路覆冰图像的边缘特征明显,因此边缘检测的方法是本文图像处理的首选。但线路覆冰的二维图像中不仅包含随机噪声和畸变,而且有些覆冰的边缘强度较低,需要对原始图像滤波去噪,并强化低强度的边缘特征。Canny算子作为一
501—·研制与开发· 杨 浩,等 基于三维重建的输电线路覆冰在线监测方法
种优化边缘检测算子,具有比较好的信噪比和检测
精度,但传统Canny算子在抑制噪声时容易错过低强度的边缘而将高强度噪声误检测,因此采用LOG
算子和Canny算子相结合的边缘检测方法[16
],具体实现见附录C。
提取2幅覆冰图像的特征边缘后需到匹配的像素点来计算其三维坐标,匹配的像素点越多则重建精度越高,但计算速度会降低。为保证一定的匹配精度并减少计算量,本文采用改进的基于灰度相
关系数的图像匹配算法[17]以提高匹配速度。基于
灰度相关系数的匹配算法思想是对于1号相机图片
中每个特征点,在2号相机图片中对应位置以指定大小窗口搜索,计算窗口内特征点的相关系数,系数最大值就是最佳匹配点,对于2号相机图片中特征点亦反过来计算1号相机图片中特征点的相关系
数。由于窗口内特征点的相关系数计算量太大,必须利用简化相关系数和设定阈值的方法逐步缩小搜
siv-011
索范围,在减小计算量的同时达到实时高精度匹配要求。
由于试验拍摄的覆冰线路距离摄像机较近,背景简单且灰度较低,匹配程度较高。对属于背景或畸形的匹配点,由于深度大的背景视差大,难以形成
匹配点,
深度小的背景形成的匹配点的像素坐标灰度同覆冰线路的灰度有所差异。因此,采用区域生
长方法进行剔除,取线路位置像素点作为生长点,按照灰度相似性生长,对应的覆冰匹配区域逐步扩大,囊括所有覆冰线路匹配点。
附录D图D1中含1 291对覆冰匹配点,丢掉背景和畸形匹配点306对,足够表征线路覆冰状况。但仍有少许匹配点在三维坐标重建后成为奇异点,由于重建点属于线路覆冰边缘,具有连续、聚合的特点,
当某一点的坐标值与其他点的坐标值偏差较大,可判定为奇异点将其剔除。另外还可通过改进算法尽量提高匹配特征点的数量,使还原的覆冰特征更加精确。
2.4 线路覆冰的三维重建
根据匹配特征点像素坐标,以及摄像机标定中求出的摄像机内部参数、旋转矩阵和平移矩阵,可由式(4
)经矩阵计算得到匹配点的三维坐标,并去除一些偏差较大的点。然后在VC++中运用OpenGL编程导入三维坐标参数,实现覆冰线路的三维点云
显示[18]
,并利用OpenGL提供的函数实现旋转等功能得到不同角度的线路覆冰点云模型,如附录E图
E1所示。
采用MATLAB对点云模型的Zw轴数据进行3次插值拟合出其覆冰模型,1号和2号相机拍摄的线路覆冰图见附录C图C1和图C2,
其三维拟合模型见附录E图E2。按照同样的方法,拍摄线路另一面的覆冰图像进行重建,
可得到完整的线路覆冰点云模型,
拟合出该段线路覆冰的三维网格模型。将每个点视为以Yw=0平面为底的立方体,对该覆冰模型上
冰所包络的体积进行计算得到其总体积为V1,未覆冰时线路体积为V0,即可得到该线路覆冰体积V为:
V=V1-V0(9)则覆冰重量为:
M=ρV
(10
)式中:ρ为冰的密度。
根据线路覆冰三维网格模型,设覆冰线路上含
m个覆冰截面,
每个截面边界有n个网点,则第j个覆冰截面的平均厚度Rj和整根导线的平均覆冰厚
度H分别为:
Rj=
∑ni=1
ij
-R0
(11
)H=
∑m
j=1
(12
)式中:rij为第j个截面中网点i对应的半径;R0为导线半径。
将重建所测厚度与人工所测厚度比较,得到导线平均覆冰厚度的精度为:
p=
H0-|H0-H|
H0
×1
00%(13
)式中:H0为人工测得的导线平均覆冰厚度。
3 试验结果与分析
对图5所示2段覆冰线路进行三维重建,计算该段线路覆冰体积、平均覆冰厚度及覆冰重量,与人工测数据进行比较,所得结果如表3所示
图5 2组线路覆冰试验
Fig.5 Two line icing 
tests—
601—2012,36(23
) 
表3 2组覆冰试验三维重建结果
Table 3 3Dreconstruction results of 2icing tests
线路覆冰前体积/mm3三维重建体积/
mm3
人工测得平均
厚度/mm
第1段43 080 217 180 17.32第2段43 080 255 430 18.20
线路三维重建测得
平均厚度/mm
人工测重/g三维重建测重/g
第1段16.427 38 16.13 15.67
第2段17.370 13 20.47 19.11
通过三维重建测得第1段线路的平均覆冰厚度16.43mm,精度94.86%,覆冰重量15.67g,精度97.15%;第2段线路平均覆冰厚度17.37mm,精度95.44%,覆冰重量19.11g,精度93.36%。可见基于三维重建的覆冰监测能够达到较高的精度。同时采用MATLAB分析重建模型数据,导出第1组试验以导线轴线距离为横坐标的厚度分布情况,如图6和图7所示。
图6 最大覆冰厚度沿线路分布图
Fig.6 Distribution diagram of the biggest ice cover
thickness along transmission lines
图7 平均覆冰厚度沿线路分布图
Fig.7 Distribution diagram of average ice
cover thickness along transmission lines
由图可得,在线路213.54mm处存在最大覆冰厚度,为42.61mm,对最大覆冰厚度的监测,可以防止一些异常状况的发生;线路平均覆冰厚度在225~265mm区段较大,通过监测较大覆冰厚度的
线路区间,能够直观精确地监测覆冰状况,是预防覆冰灾害的有效手段。实现覆冰厚度超过阈值的自动报警,并指示报警线路区间,能够使工作人员迅速预防或处理事故。另外,为保证实时的高精度监测,本文对三维重建算法进行简化,若监测主机的硬件资源足够,可以采用更为复杂的改进特征点提取算法,以提高监测精度。
4 结语
针对现有覆冰在线监测手段的不足,提出基于三维重建的线路覆冰图像在线监测法。试验结果表明,三维重建模型可清晰得到覆冰的三维轮廓,还原度高,能够直观反映覆冰形态。该方法在对覆冰的在线监测中引入了计算机视觉技术,但仍不够完善,针对其不足之处可以做出如下改进。
1)采用分辨率更高的工业相机采集图像,改进边缘检测和匹配算法,在保证精度的同时扩大监测区域。
2)利用云台旋转、平移摄像机扩大监测范围,但是需要先建立云台参数和摄像机参数矩阵之间的关系函数,使其能够自动确定摄像机标定结果。
3)可对三维重建模型进行三角剖分,再通过粘贴纹理还原覆冰表面,使其更为形象直观。
附录见本刊网络版(http://aeps.sgepri.sgcc.com.cn/aeps/ch/index.aspx)。
参考文献
[1]孙才新,司马文霞,舒立春.大气环境与电气外绝缘[M].北京:中国电力出版社,2002.
[2]蒋兴良.输电线路导线覆冰机理和三峡地区覆冰规律及影响因素研究[D].重庆:重庆大学,1997.
[3]张恒旭,刘玉田,张鹏飞.极端冰雪灾害下电网安全评估需求分析与框架设计[J].中国电机工程学报,2009,29(16):8-14.
ZHANG Hengxu,LIU Yutian,ZHANG Pengfei.Requirements analysis and framework design for power systemsecurity assessment considering extreme ice disasters[J].Proceedings of the CSEE,2009,29(16):8-14.
[4]马国明,李成榕,全江涛,等.架空输电线路覆冰监测光纤光栅拉力倾角传感器的研制[J].中国电机工程学报,2010,30(34):132-138.
林妙可被毒虫咬伤
MA Guoming,LI Chengrong,QUAN Jiangtao,et al.Development of fiber Bragg grating tension &tilt sensor foricing-monitoring of overhead transmission lines[J].Proceedingsof the CSEE,2010,30(34):132-138.
[5]阳林,郝艳捧,黎卫国,等.架空输电线路在线监测覆冰力学计算模型[J].中国电机工程学报,2010,30(19):100-105.
YANG Lin,HAO Yanpeng,LI Weiguo,et al.A mecha
nicalcalculation model for on-line icing-monitoring system of overheadtransmission lines[J].Proceedings of the CSEE,2010,30(19):100-105.
[6]秦建敏,张志栋,胡波,等.基于空气与冰的电阻特性差异实现电
草莓托
·研制与开发· 杨 浩,等 基于三维重建的输电线路覆冰在线监测方法

本文发布于:2024-09-21 03:22:42,感谢您对本站的认可!

本文链接:https://www.17tex.com/tex/2/142005.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:线路   摄像机   图像   监测   覆冰
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议