CN113009416B - 一种基于激光传感器阵列的绝缘子检测定位方法 - Google Patents
一种基于激光传感器阵列的绝缘子检测定位方法 Download PDFInfo
- Publication number
- CN113009416B CN113009416B CN202110377258.3A CN202110377258A CN113009416B CN 113009416 B CN113009416 B CN 113009416B CN 202110377258 A CN202110377258 A CN 202110377258A CN 113009416 B CN113009416 B CN 113009416B
- Authority
- CN
- China
- Prior art keywords
- insulator
- laser sensor
- sensors
- sensor array
- laser
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000012212 insulator Substances 0.000 title claims abstract description 187
- 238000001514 detection method Methods 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000000737 periodic effect Effects 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims description 22
- 238000005259 measurement Methods 0.000 claims description 17
- 238000011065 in-situ storage Methods 0.000 claims description 6
- 230000002441 reversible effect Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 abstract description 20
- 238000012545 processing Methods 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000013178 mathematical model Methods 0.000 description 4
- 229910052573 porcelain Inorganic materials 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000005202 decontamination Methods 0.000 description 2
- 230000003588 decontaminative effect Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/16—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using electromagnetic waves other than radio waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Abstract
本发明公开了一种基于激光传感器阵列的绝缘子检测定位方法,激光传感器阵列由五个水平安装且朝向相同的激光传感器组成,其中一个位于“十”字形刚体轴中心处,另外四个两两呈中心对称分布在“十”字形刚体轴的竖直刚体轴和水平刚体轴上。方法包括:基于竖直柱状绝缘子伞裙半径的周期性特征对绝缘子实现初步定位;基于绝缘子伞裙的周期性和波动性特征实现绝缘子鉴别并定位其绝缘子伞裙半径极大值点;基于绝缘子伞裙横截面是圆形的特征实现绝缘子中心位置的定位。本发明基于激光传感器阵列数据,实现了绝缘子空间位置的检测定位,数据处理量小,有效排除了复杂环境对绝缘子清扫工作的影响,具有良好的泛用性,适用于各式柱状绝缘子。
Description
技术领域
本发明涉及一种基于激光传感器阵列的绝缘子检测定位方法,属于电力系统运维领域。
背景技术
绝缘子在电力系统和电气设备中应用十分广泛,按照绝缘和连接形式的不同,绝缘子可分为狭义绝缘子、瓷套和套管。绝缘子清扫是变电站运维工作的重要一环。瓷套作为绝缘子的一种,是高压电气设备的外绝缘部件,保持其清洁状态尤为重要。传统的绝缘子清扫方式有人工登高清扫和人工地面清扫刷清扫两种,人工除污方式普遍存在着清洗效率低、难度大、除污质量对人员的依赖性大,清洗质量一致性不高,工作人员安全难以得到保障等缺点。因此在未来,使用机器人进行清扫的方式将变得愈加需要,其可以很好的适应变电站现场复杂的环境,极大的提升绝缘子清扫装置的自动化水平,有利于后续智能运检的推进和发展。
绝缘子清扫机器人在清扫过程中对绝缘子进行检测定位时,由于绝缘子分布场景复杂,使得机器视觉运算量过大,现有清扫机器人不能满足其运算要求,而且在清扫中还需要严格把握距离,导致机器人不能实现对绝缘子的快速检测定位。
发明内容
发明目的:本发明的目的在于提供一种基于激光传感器阵列的绝缘子检测定位方法,以解决现有绝缘子清扫机器人由于绝缘子所处物理、电磁环境复杂,机器视觉运算量过大,导致不能实现对绝缘子快速检测定位的问题。
技术方案:本发明提供一种基于激光传感器阵列的绝缘子检测定位方法,其特征在于,所述激光传感器阵列安装于“十”字形刚体轴上且随“十”字形刚体轴的移动而变化,所述激光传感器阵列由五个水平安装且朝向相同的激光传感器组成,其中一个位于“十”字形刚体轴中心处,另外四个两两呈中心对称分布在“十”字形刚体轴的竖直刚体轴和水平刚体轴上,竖直分布的两个激光传感器与中心处激光传感器的高度差为绝缘子相邻最大半径伞裙间距的整数倍,水平分布的两个激光传感器与中心处激光传感器的距离小于绝缘子最小半径,所述方法包括:
基于竖直分布的两个激光传感器和中心处激光传感器,检测待测物体侧面是否具有竖直柱状绝缘子伞裙半径的周期性特征,根据检测结果对绝缘子进行初步定位;
基于竖直分布的两个激光传感器和中心处激光传感器,对初步定位的物体检测其侧面是否同时具有绝缘子伞裙的周期性和波动性特征,根据检测结果判断是否为绝缘子,并在判定为绝缘子后定位绝缘子伞裙半径极大值点对应的高度;
基于水平分布的两个激光传感器,并根据绝缘子伞裙水平横截面的圆形特征,调整激光传感器阵列朝向,以使其对准绝缘子竖直中心轴。
进一步地,将中心处激光传感器记为1号传感器,竖直分布的两个激光传感器分别记为2号、3号传感器,所述初步定位具体包括以下步骤:
使激光传感器阵列原地转向,并读取1、2、3号传感器的测距结果,若测距结果满足如下约束条件,则初步认为检测到了绝缘子:
std({di|i=1,2,3})<α
式中:di表示第i个激光传感器的测距结果,α为初步检测判定阈值,std(·)表示求标准差。
进一步地,对初步定位的物体检测是否具有绝缘子伞裙的波动性特征具体包括以下步骤:
将激光传感器阵列缓慢下降2T,T为绝缘子相邻最大半径伞裙间距,在此过程中读取1、2、3号传感器测距结果的最大、最小值,根据下式进行绝缘子伞裙的波动性检测,若满足下式,则判断物体侧面具有绝缘子伞裙的波动性特征:
min(diM-dim)>δ;(i=1,2,3)
式中:dim(i=1,2,3)表示1、2、3号传感器在2T的下降过程中测得的最小距离,diM(i=1,2,3)表示1、2、3号传感器在2T的下降过程中所测得的最大距离,min(·)表示取最小值,δ表示设置的波动性检测阈值。
进一步地,对初步定位的物体检测是否具有绝缘子伞裙的周期性特征具体包括:
根据下降2T过程中1、2、3号传感器测得的最大、最小值,根据下式进行绝缘子伞裙的周期性检测,若满足下式,则判断物体侧面具有绝缘子伞裙的周期性特征:
||(std(dim),std(diM))||2<ξ;(i=1,2,3)
式中:std(·)表示求标准差,||·||2表示求取向量的2范数,ξ表示设置的周期性检测阈值。
进一步地,所述绝缘子伞裙半径极大值点对应的高度为在2T的下降过程中1、2、3号传感器测距结果最小值对应的高度。
进一步地,调整激光传感器阵列朝向具体包括:
根据水平方向的两个传感器的测距结果的大小关系和绝缘子伞裙水平横截面的圆形特征,确定激光传感器阵列的调整方向;并以水平方向的两个传感器的测距结果的差值确定激光传感器阵列的调整角度;
根据所述调整方向和调整角度,调整激光传感器阵列的朝向直至水平方向的两个传感器的测距结果的差值在允许范围内。
进一步地,将水平方向的两个传感器分别记为4号、5号传感器,所述调整角度根据下式计算得到:
Δθi=k(d5-d4);i=1,2,3,4,5
式中,d4,d5分别表示4、5号传感器的测距值,k为比例常数,Δθi表示激光传感器阵列的角度修正量。
进一步地,所述的一种基于激光传感器阵列的绝缘子检测定位方法,还包括:对绝缘子进行前期数学建模,具体为:
选取复合正弦函数对绝缘子进行近似拟合,拟合公式如下:
式中:Rmax、Rsec、Rmin依次表示绝缘子的最大半径、第二半径和最小半径,T表示绝缘子相邻最大半径伞裙的高度差,r表示伞裙高度为h时绝缘子横截面的圆形半径,εt表示随机误差,服从密度函数为下式的广义t分布:
其中,参数ν>0,g>0且有νg>2;Beta函数B(ν-1,g)≡Γ(ν-1)Γ(g)/Γ(ν-1+g);σt为εt的标准差;当ν=2时,广义t分布退化为t分布,自由度为2g;当ν→∞时,广义t分布退化为广义误差分布;
绝缘子表面模型建立如下:
式中,hmax和hmin分别是绝缘子的顶端高度和末端高度,(x,y)为绝缘子水平剖面的径向和轴向的坐标位置。
进一步地,所述的一种基于激光传感器阵列的绝缘子检测定位方法,还包括:对激光传感器及其阵列进行数学建模,具体为:
激光传感器的位姿由(x,y,z,θ)表示,其中(x,y,z)表示激光传感器所在的三维空间坐标,θ表示激光传感器的朝向;
激光传感阵列具备如下位姿变化能力:
前进/后退:激光传感器阵列沿着当前朝向在水平面上进行前进/后退,所有的激光传感器的高度不变、朝向不变,但在地面投影点坐标会同步改变;
转向:激光传感器阵列在原地改变自身的整体朝向,所有传感器的高度不变,朝向同步改变,安装于竖直轴的传感器位置不变,安装于水平轴的传感器绕竖直轴旋转,旋转方向与角度大小与阵列整体朝向角度变化相同;
升降:激光传感器阵列沿竖直反向进行升降,所有传感器的朝向不变、位置不变,高度同步增减。
进一步地,中心处激光传感器离地高度为100cm,竖直分布的两个激光传感器离地高度分别为100-5T cm和100+5T cm,T为绝缘子相邻最大半径伞裙间距;水平分布的两个激光传感器与中心处激光传感器的距离为10cm。
有益效果:与现有技术相比,本发明具有如下有益效果:
1)本发明基于激光传感器阵列数据,通过检测绝缘子侧面的周期性和波动性特征,实现了绝缘子空间位置的检测定位及最大伞裙半径对应高度和绝缘子竖直中心轴的定位;
2)本发明没有使用现今常用的机器视觉算法进行绝缘子检测,而是使用激光传感器阵列的测距结果进行数据分析,数据处理量小,降低了搭载算法的硬件的计算能力的需求,完全可以使用PLC等可靠性更高的工业编程方法进行实现,降低了硬件成本,并兼顾了稳定性和快速性;
3)本发明的数据来源是深度信息,基于深度信息进行后续的特征提取和从而实现绝缘子检测定位,深度信息作为输入源本质上可以更好的实现清扫过程中的绝缘子保护,安全性高;
4)本发明针对绝缘子的泛式几何特征规律进行传感器激光阵列的架构,一方面特殊严苛的筛选条件可以有效排除绝缘子所处的复杂环境影响,一方面针对绝缘子的相关特征参数在建模加以考虑,使得本发明具有良好的泛用性,适用于各式柱状绝缘子,且预留检测阈值作为调控参数也可以解决现实中检测精度和绝缘子自身尺寸等原因带来的检测误差问题。
附图说明
图1为本发明实施例的一种基于激光传感器阵列的绝缘子检测定位方法的详细流程图;
图2为激光传感器阵列结构的正视图;
图3为竖直同轴激光传感器阵列进行绝缘子侧面周期性和波动性检测示意图;
图4为水平等高激光传感器阵列进行绝缘子中心轴朝向校对示意图;
图5为典型绝缘子数学模型截面图;
图6为激光传感器测距过程数学建模示意图;
图7为1、2、3号传感器测距最大、最小值记录;
图8为绝缘子伞裙半径极大值点对应高度定位示意图;
图9为激光传感器阵列朝向校对过程;
图10为激光传感器阵列朝向校对完成结果展示图。
具体实施方式
下面结合具体实施例对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如前所述,现有技术中存在由于绝缘子所处物理、电磁环境复杂,机器视觉运算量过大,导致现有机器人不能实现对绝缘子进行快速检测定位的问题,为此,本发明提供了一种基于激光传感器阵列的绝缘子检测定位方法。
首先需要说明的是,本申请的“绝缘子”作广义理解,按照绝缘和连接形式的不同,绝缘子可分为狭义绝缘子、瓷套和套管。下面以狭义绝缘子为例进行说明。
本发明实施例的一种基于激光传感器阵列的绝缘子检测定位方法主要包括绝缘子柱状特征检测、绝缘子边缘波动性与周期性特征检测及伞裙最大半径高度定位、绝缘子中心轴定位,其具体流程如图1所示。
在本发明实施例中,首先对激光传感器阵列和绝缘子进行前期数学建模,而后对绝缘子柱状特征检测(以下简称绝缘子初步检测)、绝缘子边缘波动性与周期性特征检测及伞裙最大半径高度定位(以下简称绝缘子特征检测)、绝缘子中心轴定位(以下简称阵列朝向校对)进行具体介绍。
0)前期数学建模
1.绝缘子数学建模
根据绝缘子特征,其横截面近似视作圆形,且圆形半径与横截面离地高度呈现明显的周期性变化规律,由于本申请的检测依据与圆形半径的最大、最小值相关,与曲线形状无关,本申请选取复合正弦函数进行近似拟合,拟合公式如下:
式中:式中:Rmax、Rsec、Rmin依次表示绝缘子的最大半径、第二半径和最小半径,T表示绝缘子相邻最大半径伞裙的高度差,r表示伞裙高度为h时绝缘子横截面的圆形半径,εt表示随机误差,服从密度函数为下式的广义t分布:
其中,参数ν>0,g>0且有νg>2;Beta函数B(ν-1,g)≡Γ(ν-1)Γ(g)/Γ(ν-1+g);σt为εt的标准差;当ν=2时,广义t分布退化为t分布,自由度为2g;当ν→∞时,广义t分布退化为广义误差分布。
故据此建立绝缘子表面模型如下:
式中hmax和hmin分别是绝缘子的顶端高度和末端高度,(x,y)为绝缘子水平剖面的径向和轴向的坐标位置。
2.激光传感器及其阵列数学建模
激光传感器水平放置,故其位姿可由(x,y,z,θ)表示,其中(x,y,z)表示传感器所在的三维空间坐标,θ表示激光传感器的朝向。
所有的激光传感器均水平放置,且朝向一致,区别在于安装位置,各个激光传感器的安装方法具体为:三个安装于同一竖直刚体轴,离地高度分别为100-5T cm、100cm、100+5T cm;另外两个位于100cm高度水平面,且关于竖直轴在此高度平面投影点中心对称,且均距离竖直轴投影点10cm,安装于此三点所确定的水平刚体轴上,故所述激光传感器安装于由竖直刚体轴和水平刚体轴连接的“十”字形结构上且随“十”字形刚体轴结构移动而变化。为便于表示,规定安装于竖直轴且离地100cm的传感器为1号传感器,安装于竖直轴且离地100-5T cm的传感器为2号传感器,安装于竖直轴且离地100+5T cm的传感器为3号传感器,安装于水平轴且在激光传感器阵列朝向左侧的传感器为4号传感器,安装于水平轴且在激光传感器阵列朝向右侧的传感器为5号传感器,后续将按此方式进行表述,按此方法表述的传感器阵列模型如图2所示。
3.激光传感器阵列位姿变换约束
激光传感器阵列安装在机器人结构上,故规定其具备如下位姿变化能力,且在进行位姿变化时相关的约束介绍如下:
1)前进/后退:激光传感器阵列可以沿着当前朝向在水平面上进行前进/后退,激光传感器阵列位姿变化如下:
式中形如x0表示变量初始值,形如Δx表示变量的增加量。
2)转向:激光传感器阵列可以在原地改变自身的整体朝向,所有传感器的高度不变,朝向同步改变,但是部分传感器的位置会有变化,具体描述如下:
3)升降:竖直刚体轴应具备伸缩功能,可以进行升降操作,在进行升降时,所有传感器的朝向不变、位置不变,高度同步增减,具体位姿描述如下:
1)绝缘子初步检测:
使激光传感器阵列开始原地转向,并读取1、2、3号传感器的测距结果,若测距结果满足如下约束,则认为检测到了可能是绝缘子的器件:
std({di|i=1,2,3})<α (6)
式中:di表示第i个传感器的测距结果,α为初步检测判定阈值,std(·)表示求标准差。
若激光传感器阵列对向绝缘子,由于其同轴高度差为绝缘子周期的若干整数倍,根据绝缘子侧面截面圆半径的周期性特征,同轴激光传感器的测距结果应该近似相等,据此实现绝缘子的初步检测。
2)绝缘子特征检测
绝缘子的初步定位不能直接确认物体是否为绝缘子,因为许多器件,如墙面、电力屏柜等均会导致同轴传感器测距结果接近,故需要针对绝缘子的伞裙特征进行进一步的特征提取,从而确定激光传感器阵列所朝向器件为绝缘子,本申请将其划分为两项检测:波动性检测和周期性检测。下面将对其进行具体介绍,只有同时满足波动性和周期性特征,才能确认其为绝缘子。
1.绝缘子侧面波动性检测
绝缘子的伞裙结构势必会使得激光传感器在高度发生变化时的测距结果呈现明显的波动性,故激光传感器阵列在完成初步检测后缓慢下降2T,在此过程中读取1、2、3号激光传感器测距结果,为降低数据存储量和运算量,只保留其各自测距结果的最大、最小值,而后根据下式进行波动性检测:
min(diM-dim)>δ;(i=1,2,3) (7)
式中:dim(i=1,2,3)表示1、2、3号传感器在2T的下降过程中测得的最小距离,diM(i=1,2,3)表示所测得的最大距离,min(·)表示取最小值,δ表示设置的波动性检测阈值,一般略小于绝缘子最大半径与最小半径差。
2.绝缘子侧面周期性检测
绝缘子的伞裙结构呈现除了明显的波动性,还呈现明显的周期性,因为各个伞裙的结构近似相同,所以不同高度的激光传感器因为其相距高度差为伞裙高度的整数倍,测距结果的波动情况应近似相等,故依旧根据下降2T过程中各个传感器测距结果的最大、最小值,根据下式进行绝缘子的周期性检测:
||(std(dim),std(diM))||2<ξ;(i=1,2,3) (8)
式中:dim(i=1,2,3)表示1、2、3号传感器在2T的下降过程中测得的最小距离,diM(i=1,2,3)表示所测得的最大距离,std(·)表示求标准差,||·||2表示求取向量的2范数,ξ表示设置的周期性检测阈值。
此外,在下降2T的过程中,在读取激光传感器的测距结果的同时,也将记录1、2、3号传感器测距结果最小值对应的高度。测距结果最小意味着伞裙半径最大,记录此高度以实现绝缘子伞裙半径极大值点对应高度的定位。
3)阵列朝向校对
若完成绝缘子侧面的波动性和周期性检测,确认激光传感器已经大致对向了绝缘子,则需根据在确认所监测物体为绝缘子且完成对于绝缘子级伞裙外径极大值点的定位后,根据传感器4和5测距结果的差值控制机器人转向,从而使得机器人调整朝向对准绝缘子,为了方便起见,本申请的转向控制选取比例控制器,即转向角度△θ可由下式计算得到:
Δθi=k(d5-d4);i=1,2,3,4,5 (9)
式中d4,d5表示4、5号传感器的测距值,k为比例常数,Δθi表示激光传感器阵列的角度修正量。
下面通过具体实例对本发明的方案作进一步阐述:
1)数学模型参数设置及激光传感器补充建模
本发明在MATLAB中搭建数学模型验证算法的可行性。首先选取一个典型的绝缘子结构进行数学建模,经过相关资料学习后确立绝缘子参数如下:
表1典型绝缘子参数
而后按照公式(1)-(2)进行建模后,可得绝缘子相关的数学模型如图5所示。
在得到绝缘子模型后,由于未使用真实的激光传感器,需对激光传感器的测距过程进行数据建模,结合之前公式(3)-(5)进行的激光传感器建模结果,对激光传感器的测距过程进行建模,其求解过程如下:
式中(x0,y0)为激光传感器所在坐标,r(z0)为z0对应的绝缘子半径,可结合公式(1)与表1求得,θ为激光传感器朝向。
令公式(10)在有交点时两个根为(x1,y1)和(x2,y2),而后对其进行如下约束即可得到激光传感器射出激光打在绝缘子表面的坐标(x,y):
而后设置激光传感器阵列初始坐标信息如表2所示,由此完成激光传感器测距过程建模,其示意图如图6所示。
表2激光传感器初始位姿参数
2)绝缘子初步检测
使激光传感器阵列开始原地转向,激光传感器阵列的位姿变换符合如下约束:
而后读取1、2、3号传感器的测距结果,若测距结果满足如下约束,则认为检测到了可能是绝缘子的器件:
std({di|i=1,2,3})<α
式中:di表示第i个传感器的测距结果,α为初步检测判定阈值,本实例取0.01。
3)绝缘子波动性特征检测
在完成初步检测后,使激光传感器阵列缓慢下降20cm,在此过程中读取1、2、3号激光传感器测距结果,记录结果如图7所示,而后根据下式进行波动性检测:
min(diM-dim)>δ;(i=1,2,3)
式中:dim(i=1,2,3)表示1、2、3号传感器在20cm的下降过程中测得的最小距离,diM(i=1,2,3)表示所测得的最大距离,min(·)表示取最小值,δ表示设置的波动性检测阈值,一般略小于绝缘子最大半径与最小半径差,本实例取6cm。
4)绝缘子周期性检测
根据下降20cm过程中各个传感器测距结果的最大、最小值,根据下式进行绝缘子的周期性检测:
||(std(dim),std(diM))||2<ξ;(i=1,2,3)
式中:dim(i=1,2,3)表示1、2、3号传感器在20cm的下降过程中测得的最小距离,diM(i=1,2,3)表示所测得的最大距离,std(·)表示求标准差,||·||2表示求取向量的2范数,ξ表示设置的周期性检测阈值,本实例设置为0.1。
5)绝缘子伞裙最大半径高度对应高度定位
此外,在下降20cm的过程中,在读取激光传感器的测距结果的同时,也将记录1、2、3号传感器测距结果最小值对应的高度,以实现绝缘子伞裙半径极大值点对应高度的定位,检测结果如图8所示。
6)绝缘子朝向校对
在确认所监测物体为绝缘子且完成对于绝缘子级伞裙外径极大值点的定位后,根据传感器4和5测距结果的差值控制机器人转向,从而使得机器人调整朝向对准绝缘子,为了方便起见,本发明的转向控制选取比例控制器,即转向角度△θ可由式(9)计算得到。
式(9)中,k值的选取会影响机器人进行朝向校对时的表现,k值越大,阵列转动越剧烈,但可能导致不收敛引发振荡且可能会破坏绝缘子,k值过小,阵列收敛过慢降低了机器人清理的效率,考虑到本实例为工业模型,此处不过多探究,选定k值为0.0025。
在朝向校对过程仿真中,实时记录机器人的横向传感器阵列(1、4、5号传感器)测距结果绘制曲线如图9所示,可见在经过13次方向微调后传感器4和传感器5测距值相差不到0.1cm,符合精度要求,表明朝向校对成功。
7)校对完成,结果展示。
至此,朝向校对完成,结果如图10所示。
由以上实施例可以看出,本发明的一种基于激光传感器阵列的绝缘子检测定位方法,具有以下有益效果:
1)本发明基于激光传感器阵列数据,通过检测绝缘子侧面的周期性和波动性特征,实现了绝缘子空间位置的检测定位及最大伞裙半径对应高度和绝缘子竖直中心轴的定位,是机器人绝缘子清扫装置的必要步骤;
2)本发明没有使用现今常用的机器视觉算法进行绝缘子检测,而是使用激光传感器阵列的测距结果进行数据分析,数据处理量小,降低了搭载算法的硬件的计算能力的需求,完全可以使用PLC等可靠性更高的工业编程方法进行实现,降低了硬件成本,并兼顾了稳定性和快速性;
3)本发明的数据来源是深度信息,基于深度信息进行后续的特征提取和从而实现绝缘子检测定位,深度信息作为输入源本质上可以更好的实现清扫过程中的绝缘子保护,安全性高;
4)本发明针对绝缘子的泛式几何特征规律进行传感器激光阵列的架构,一方面特殊严苛的筛选条件可以有效排除绝缘子所处的复杂环境影响,一方面针对绝缘子的相关特征参数在建模加以考虑,使得本发明具有良好的泛用性,适用于各式柱状绝缘子,且预留检测阈值作为调控参数也可以解决现实中检测精度和绝缘子自身尺寸等原因带来的检测误差问题。
以上已以较佳实施例公布了本发明,然其并非用以限制本发明,凡采取等同替换或等效变换的方案所获得的技术方案,均落在本发明的保护范围内。
Claims (9)
1.一种基于激光传感器阵列的绝缘子检测定位方法,其特征在于,所述激光传感器阵列安装于“十”字形刚体轴上且随“十”字形刚体轴的移动而变化,所述激光传感器阵列由五个水平安装且朝向相同的激光传感器组成,其中一个位于“十”字形刚体轴中心处,另外四个两两呈中心对称分布在“十”字形刚体轴的竖直刚体轴和水平刚体轴上,竖直分布的两个激光传感器与中心处激光传感器的高度差为绝缘子相邻最大半径伞裙间距的整数倍,水平分布的两个激光传感器与中心处激光传感器的距离小于绝缘子最小半径,所述方法包括:
基于竖直分布的两个激光传感器和中心处激光传感器,检测待测物体侧面是否具有竖直柱状绝缘子伞裙半径的周期性特征,根据检测结果对绝缘子进行初步定位;
基于竖直分布的两个激光传感器和中心处激光传感器,对初步定位的物体检测其侧面是否同时具有绝缘子伞裙的周期性和波动性特征,根据检测结果判断是否为绝缘子,并在判定为绝缘子后定位绝缘子伞裙半径极大值点对应的高度;
基于水平分布的两个激光传感器,并根据绝缘子伞裙水平横截面的圆形特征,调整激光传感器阵列朝向,以使其对准绝缘子竖直中心轴;
其中,将中心处激光传感器记为1号传感器,竖直分布的两个激光传感器分别记为2号、3号传感器,所述绝缘子伞裙半径极大值点对应的高度为将激光传感器阵列缓慢下降2T的过程中1、2、3号传感器测距结果最小值对应的高度,其中T为绝缘子相邻最大半径伞裙间距。
2.根据权利要求1所述的方法,其特征在于,所述初步定位具体包括以下步骤:
使激光传感器阵列原地转向,并读取1、2、3号传感器的测距结果,若测距结果满足如下约束条件,则初步认为检测到了绝缘子:
std({di|i=1,2,3})<α
式中:di表示第i个激光传感器的测距结果,α为初步检测判定阈值,std(·)表示求标准差。
3.根据权利要求1所述的方法,其特征在于,对初步定位的物体检测是否具有绝缘子伞裙的波动性特征具体包括以下步骤:
将激光传感器阵列缓慢下降2T,在此过程中读取1、2、3号传感器测距结果的最大、最小值,根据下式进行绝缘子伞裙的波动性检测,若满足下式,则判断物体侧面具有绝缘子伞裙的波动性特征:
min(diM-dim)>δ;i=1,2,3
式中:dim,i=1,2,3表示1、2、3号传感器在2T的下降过程中测得的最小距离,diM,i=1,2,3表示1、2、3号传感器在2T的下降过程中所测得的最大距离,min(·)表示取最小值,δ表示设置的波动性检测阈值。
4.根据权利要求3所述的方法,其特征在于,对初步定位的物体检测是否具有绝缘子伞裙的周期性特征具体包括:
根据下降2T过程中1、2、3号传感器测得的最大、最小值,根据下式进行绝缘子伞裙的周期性检测,若满足下式,则判断物体侧面具有绝缘子伞裙的周期性特征:
||(std(dim),std(diM))||2<ξ;i=1,2,3
式中:std(·)表示求标准差,||·||2表示求取向量的2范数,ξ表示设置的周期性检测阈值。
5.根据权利要求1所述的方法,其特征在于,调整激光传感器阵列朝向具体包括:
根据水平方向的两个传感器的测距结果的大小关系和绝缘子伞裙水平横截面的圆形特征,确定激光传感器阵列的调整方向;并以水平方向的两个传感器的测距结果的差值确定激光传感器阵列的调整角度;
根据所述调整方向和调整角度,调整激光传感器阵列的朝向直至水平方向的两个传感器的测距结果的差值在允许范围内。
6.根据权利要求5所述的方法,其特征在于,将水平方向的两个传感器分别记为4号、5号传感器,所述调整角度根据下式计算得到:
Δθi=k(d5-d4);i=1,2,3,4,5
式中,d4,d5分别表示4、5号传感器的测距值,k为比例常数,△θi表示激光传感器阵列的角度修正量。
7.根据权利要求1所述的方法,其特征在于,还包括:对绝缘子进行前期数学建模,具体为:
选取复合正弦函数对绝缘子进行近似拟合,拟合公式如下:
式中:Rmax、Rsec、Rmin依次表示绝缘子的最大半径、第二半径和最小半径,T表示绝缘子相邻最大半径伞裙的高度差,r表示伞裙高度为h时绝缘子横截面的圆形半径,εt表示随机误差,服从密度函数为下式的广义t分布:
其中,参数ν>0,g>0且有νg>2;Beta函数B(ν-1,g)≡Γ(ν-1)Γ(g)/Γ(ν-1+g);σt为εt的标准差;当ν=2时,广义t分布退化为t分布,自由度为2g;当ν→∞时,广义t分布退化为广义误差分布;
绝缘子表面模型建立如下:
式中,hmax和hmin分别是绝缘子的顶端高度和末端高度,(x,y)为绝缘子水平剖面的径向和轴向的坐标位置。
8.根据权利要求1所述的方法,其特征在于,还包括:对激光传感器及其阵列进行数学建模,具体为:
激光传感器的位姿由(x,y,z,θ)表示,其中(x,y,z)表示激光传感器所在的三维空间坐标,θ表示激光传感器的朝向;
激光传感阵列具备如下位姿变化能力:
前进/后退:激光传感器阵列沿着当前朝向在水平面上进行前进/后退,所有的激光传感器的高度不变、朝向不变,但在地面投影点坐标会同步改变;
转向:激光传感器阵列在原地改变自身的整体朝向,所有传感器的高度不变,朝向同步改变,安装于竖直轴的传感器位置不变,安装于水平轴的传感器绕竖直轴旋转,旋转方向与角度大小与阵列整体朝向角度变化相同;
升降:激光传感器阵列沿竖直反向进行升降,所有传感器的朝向不变、位置不变,高度同步增减。
9.根据权利要求1所述的方法,其特征在于,中心处激光传感器离地高度为100cm,竖直分布的两个激光传感器离地高度分别为100-5T cm和100+5T cm;水平分布的两个激光传感器与中心处激光传感器的距离为10cm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110377258.3A CN113009416B (zh) | 2021-04-08 | 2021-04-08 | 一种基于激光传感器阵列的绝缘子检测定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110377258.3A CN113009416B (zh) | 2021-04-08 | 2021-04-08 | 一种基于激光传感器阵列的绝缘子检测定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113009416A CN113009416A (zh) | 2021-06-22 |
CN113009416B true CN113009416B (zh) | 2024-03-12 |
Family
ID=76388120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110377258.3A Active CN113009416B (zh) | 2021-04-08 | 2021-04-08 | 一种基于激光传感器阵列的绝缘子检测定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113009416B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1719269A (zh) * | 2005-04-08 | 2006-01-11 | 华北电力大学(北京) | 绝缘子带电检测仪及其实施方法 |
CN105004504A (zh) * | 2015-07-14 | 2015-10-28 | 清华大学深圳研究生院 | 复合绝缘子的起振风速评估方法及选型方法 |
CN106952249A (zh) * | 2017-02-20 | 2017-07-14 | 广东电网有限责任公司惠州供电局 | 基于交比不变性的绝缘子串中轴线提取方法 |
CN109142999A (zh) * | 2018-09-14 | 2019-01-04 | 平顶山学院 | 变电站局部放电巡检定位设备的位置检测装置及方法 |
CN109539995A (zh) * | 2018-11-19 | 2019-03-29 | 国网四川省电力公司电力科学研究院 | 一种绝缘子爬电距离自动测量装置 |
CN109801283A (zh) * | 2019-01-25 | 2019-05-24 | 西安工程大学 | 一种基于水滴偏移距离的复合绝缘子憎水性等级判定方法 |
CN110665875A (zh) * | 2019-11-11 | 2020-01-10 | 国网江苏省电力有限公司检修分公司 | 一种绝缘子清扫装置及方法 |
CN110909751A (zh) * | 2019-11-26 | 2020-03-24 | 长沙理工大学 | 用于变电站绝缘子清扫机器人的视觉识别方法、系统及介质 |
CN111161930A (zh) * | 2020-01-03 | 2020-05-15 | 西北核技术研究院 | 一种具有复合结构的真空绝缘子及其制备方法 |
CN111754465A (zh) * | 2020-06-04 | 2020-10-09 | 四川大学 | 一种绝缘子定位与掉串检测方法 |
CN111784633A (zh) * | 2020-05-26 | 2020-10-16 | 西安理工大学 | 一种面向电力巡检视频的绝缘子缺损自动检测算法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050064111A1 (en) * | 2003-09-23 | 2005-03-24 | Hiller Nathan David | Method for forming doping superlattices using standing electromagnetic waves |
US20120051378A1 (en) * | 2010-08-31 | 2012-03-01 | Aravinda Kar | Photodetection |
JP5934914B2 (ja) * | 2011-05-13 | 2016-06-15 | パナソニックIpマネジメント株式会社 | レーザアレイ光源ユニット |
WO2018146146A1 (en) * | 2017-02-08 | 2018-08-16 | Trinamix Gmbh | Detector for an optical detection of at least one object |
-
2021
- 2021-04-08 CN CN202110377258.3A patent/CN113009416B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1719269A (zh) * | 2005-04-08 | 2006-01-11 | 华北电力大学(北京) | 绝缘子带电检测仪及其实施方法 |
CN105004504A (zh) * | 2015-07-14 | 2015-10-28 | 清华大学深圳研究生院 | 复合绝缘子的起振风速评估方法及选型方法 |
CN106952249A (zh) * | 2017-02-20 | 2017-07-14 | 广东电网有限责任公司惠州供电局 | 基于交比不变性的绝缘子串中轴线提取方法 |
CN109142999A (zh) * | 2018-09-14 | 2019-01-04 | 平顶山学院 | 变电站局部放电巡检定位设备的位置检测装置及方法 |
CN109539995A (zh) * | 2018-11-19 | 2019-03-29 | 国网四川省电力公司电力科学研究院 | 一种绝缘子爬电距离自动测量装置 |
CN109801283A (zh) * | 2019-01-25 | 2019-05-24 | 西安工程大学 | 一种基于水滴偏移距离的复合绝缘子憎水性等级判定方法 |
CN110665875A (zh) * | 2019-11-11 | 2020-01-10 | 国网江苏省电力有限公司检修分公司 | 一种绝缘子清扫装置及方法 |
CN110909751A (zh) * | 2019-11-26 | 2020-03-24 | 长沙理工大学 | 用于变电站绝缘子清扫机器人的视觉识别方法、系统及介质 |
CN111161930A (zh) * | 2020-01-03 | 2020-05-15 | 西北核技术研究院 | 一种具有复合结构的真空绝缘子及其制备方法 |
CN111784633A (zh) * | 2020-05-26 | 2020-10-16 | 西安理工大学 | 一种面向电力巡检视频的绝缘子缺损自动检测算法 |
CN111754465A (zh) * | 2020-06-04 | 2020-10-09 | 四川大学 | 一种绝缘子定位与掉串检测方法 |
Non-Patent Citations (7)
Title |
---|
"基于周期性结构特征的绝缘子掉片检测算法";王海涛;《计算机测量与控制》;20180831;第26卷(第8期);213-218 * |
"基于感受野模块的绝缘子实时识别定位方法";吉志朋;《电工电气》;20200930(第9期);19-24 * |
"基于机器视觉深度学习的绝缘子前端定位方法";苏超;《广东电力》;20190930;第32卷(第9期);33-41 * |
"激光诱导击穿光谱结合CARS-PLSR法快速定量检测绝缘子污秽";王乃啸;《中国电机工程学报》;20200229;第40卷(第4期);1378-1387 * |
"盘形悬式绝缘子手持式气动清扫装置的设计";毕睿华;《电瓷避雷器》;20151231(第6期);8-13 * |
"聚合物绝缘子表面微孔阵列构筑";霍艳坤;《太赫兹科学与电子信息学报》;20190430;第17卷(第2期);332-338 * |
"超高压输电线悬垂绝缘子串激光除冰机器人构型设计研究";邝江华;《电力科学与工程》;20210131;第37卷(第1期);25-33 * |
Also Published As
Publication number | Publication date |
---|---|
CN113009416A (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109623656B (zh) | 基于厚度在线检测的移动式双机器人协同打磨装置及方法 | |
CN109176517B (zh) | 基于末端名义点约束的串联工业机器人连杆参数标定方法 | |
CN109732600B (zh) | 一种全自动顺序多站式测量系统及测量方法 | |
CN110793543B (zh) | 基于激光扫描的电力巡检机器人的定位导航精度测量装置及方法 | |
CN103759635A (zh) | 一种精度与机器人无关的扫描测量机器人检测方法 | |
CN104816307A (zh) | 工业机器人精准制孔的四点法向调平方法 | |
CN109454281A (zh) | 一种机器人铣削加工中的螺旋桨工件坐标系标定方法 | |
CN105399020A (zh) | 三维空间平面追踪控制方法、系统及高空作业设备 | |
CN108161935B (zh) | 机器人基坐标系校准的方法及装置 | |
CN109760006B (zh) | 一种基于视觉参考件的核机器人快速定位方法 | |
CN110883774A (zh) | 机器人关节角零位标定系统、方法及存储介质 | |
CN101266479A (zh) | 电极智能检测系统 | |
CA3062081A1 (en) | Autonomous modification of waterjet cutting systems | |
CN113009416B (zh) | 一种基于激光传感器阵列的绝缘子检测定位方法 | |
Ma et al. | A fast and robust seam tracking method for spatial circular weld based on laser visual sensor | |
CN114734436B (zh) | 一种机器人的编码器校准方法、装置及机器人 | |
CN105538309B (zh) | 一种有限传感能力的机器人障碍物动态识别算法 | |
KR100994742B1 (ko) | 3차원 측정기의 이동 경로에 대한 충돌 검출 및 경과점 생성 방법 | |
CN113845064B (zh) | 一种具有圆形支脚的物料承载装置的定位方法及系统 | |
CN110561400A (zh) | 圆周均匀分布零件的高效精准定位系统及定位方法 | |
JP2022176229A (ja) | 建築作業装置および建築作業方法 | |
CN114988288A (zh) | 一种2d激光器目标位置检测系统及方法 | |
CN111553083B (zh) | 共交点多孔光轴参数调整量优化计算方法、系统及介质 | |
CN112557072A (zh) | 采掘设备悬臂空间自由度的标定方法及装置 | |
Yusen et al. | A method of welding path planning of steel mesh based on point cloud for welding robot |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |