CN108509845B - 基于特征拟合的变电站仪表设备巡检的视觉检测方法 - Google Patents
基于特征拟合的变电站仪表设备巡检的视觉检测方法 Download PDFInfo
- Publication number
- CN108509845B CN108509845B CN201810129335.1A CN201810129335A CN108509845B CN 108509845 B CN108509845 B CN 108509845B CN 201810129335 A CN201810129335 A CN 201810129335A CN 108509845 B CN108509845 B CN 108509845B
- Authority
- CN
- China
- Prior art keywords
- image
- value
- pointer
- calculating
- instrument
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/20—Scenes; Scene-specific elements in augmented reality scenes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/02—Recognising information on displays, dials, clocks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于特征拟合的变电站仪表设备巡检的视觉检测方法,采用基于特征拟合的仪表指针定位技术,通过指针边缘特征的环形投影、边缘统计和拟合分析,计算指针两条夹边及其交点。同时,在刻度值计算方面,选择指针最近邻的2个刻度进行线性插值,通过该方法,本发明能准确估算当前的刻度值,避免了因仪表中心定位偏差或图像倾斜而导致的刻度值计算偏差。
Description
技术领域
本发明属于机器视觉的技术领域,涉及一种变电站的仪表设备巡检的视觉检测算法,尤其涉及一种基于特征拟合的变电站仪表设备巡检的视觉检测方法。
背景技术
随着科技的不断发展和变电站管理的逐渐升级,绝大部分的变电站已经实现无人值班,但由于自动化技术和功能范围并不完善,很多设备的运行数据不能通过现有的在线检测设备简单获取,往往需要安排工作人员定期巡检。
自动巡检视觉检测技术的引入,可以显著弥补现有电网的监控技术的不足。通过变电站可见光和红外图像采集,识别和分析,可以实时诊断系统运行的重要参数,并监控异常情况,为以后的统计分析提供大数据支持。
在变电站自动巡检过程中,常需要识别各类指针式仪表,例如:压力表,油压表等,因此,仪表定位识别是自动巡检的一个核心子系统。传统的指针式仪表识别算法主要利用图像匹配或图像细化技术去定位指针并计算刻度。但在光照复杂,且干扰众多的环境下,精确的图像匹配很难实现;而基于多帧图像的识别系统,往往要采集若干张仪表图片,且对现场光照和相机定位都有较高的要求。但在动态巡检过程中难以满足此类要求,容易导致计算偏差甚至识别失败。同时,在刻度值的分析计算过程中,仪表中心的计算偏差或图像倾斜等因素往往会影响许多传统算法的计算精确度。
发明内容
本发明主要面向变电站的仪表自动巡检的视觉检测需求,针对现有技术存在巡检过程中采集图像质量不够稳定,易受到各种光照、平台振动或其它外界因素干扰等技术问题,提出一种基于特征拟合的变电站仪表设备巡检的视觉检测方法。该方法采用基于特征拟合的仪表指针定位技术,通过指针边缘特征的环形投影、边缘统计和拟合分析,计算指针两条夹边及其交点。同时,在刻度值计算方面,选择指针最近邻的2个刻度进行线性插值,从而能准确估算当前的刻度值,避免了因仪表中心定位偏差或图像倾斜而导致的刻度值计算偏差。
为了解决现有技术存在的技术问题,本发明的技术方案如下:
步骤1.采集一张变电站巡检仪表图像,具体通过CCD相机拍摄目标图像并将图像转换成灰度图像G;
步骤2.针对图像G进行预处理,去除噪声点,生成图像G1。
步骤3.利用边缘算子提取G1的边缘信息,生成图像G2。
步骤4.采用Hough算法分析图像G2中圆形目标,计算出仪表中心C。
步骤5.围绕仪表中心C,对仪表指针所在的环形区域进行投影映射,生成图像G3。
步骤6.在图像G3中,计算符合指针边缘特征的统计信息值。
步骤7:依据统计信息值,计算指针的两条边线,并计算边线交点T。
步骤8:针对图像G1,围绕仪表中心C,对仪表刻度所在的环形区域进行投影映射,生成图像G4。
步骤9:对图像G4进行二值化处理,并通过区域标记算法定位每个刻度的中心区域,确定每个刻度坐标,并映射回原始图像G1空间。
步骤10:在图像G1坐标空间中,过圆心C点,计算指针针尖T与邻近2个刻度之间夹角,并通过线性插值计算指针所指的刻度值。
步骤1所述的采集一张变电站巡检仪表的图像,具体过程如下:通过工业相机拍摄目标图像,并转换成灰度图像G。转换具体如下:
Gray=R×0.3+G×0.59+B×0.11
所述的对目标图像中的每一个像素点处理,就能将其转换成灰度图像G;其中R、G、B为一个像素点的值的基色分量;
步骤2所述的仪表图像预处理,主要是去除噪声点,一般通过中值滤波方法除去灰度图像G的噪声点,生成图像G1。具体过程如下:
2-1、选定5×5的邻域,将邻域内的像素按灰度等级进行排序,确定其中值;
2-2、将确定的中值赋予滤波后图像中的相应像素点,所述的相应像素点即进行中值滤波的那个像素点;
步骤3所述的计算G1的边缘信息,一般采用Canny或Sobel算子提取边缘信息,生成图像G2。该边缘信息主要用于定位仪表中心和指针位置。
步骤4.采用Hough算法提取图像G2中的圆形目标,并定位圆形中心C。
步骤5.围绕仪表中心C,对仪表指针所在的环形区域进行投影映射,生成图像G3。考虑仪表的形状特征,一般的环形区域可以选择在半径为(r1,r2)区域。
其中,r2的长度一般近似为指针针尖T到仪表中心C的距离。r1可以选择为r2的1/3到2/3长度。
G3的大小为(nW,nH),其中,nW=2πr2,nH=r2-r1,此时,矩形上任一点P(x,y)可以映射到仪表上的点P'(x',y'),此时,
P'距离中心点的距离为r1+y,对应弧度为x/r2;
x'=(r1+y)×cos(x/r2) (1)
y'=(r1+y)×sin(x/r2) (2)
步骤6.在图像G3中,计算符合指针边缘特征的统计信息值。其中,指针边缘特征包括指针的左边线和右边线。该特征在投影视图G3中,可以用特征三元组(α,β,d)来描述,如图2所示。考虑到投影轮廓值与X坐标相关,因此统计信息值F可以由公式(3)表示:
F=f((α,β,d),x) (3)
该公式表示为以位置x处的垂线为基准,左侧边线倾角为α,右侧倾角为β,且两侧边线在底部离垂线间距均为d/2。此时,G3中所包含的边缘点数即为统计信息值。具体计算过程如下:
6.1确定α和β搜索范围:考虑到指针夹角一般相对固定,但考虑工艺差异和拍照角度差异性,夹角可能会有一定的变化,但变化幅度较小。假设仪表指针夹角为θ(该值由用户输入,一般在3-6°),则α和β搜索的范围一般在(θ/2-2,θ/2+2)即可。
6.2左侧边缘信息统计:以G3图像底部为基准,当倾角为(θ/2-2,θ/2+2)时,以0.2°间隔计算G3图像在水平方向的投影值SL(α,x)。
6.3右侧边缘信息统计:以G3图像底部为基准,当倾角为(θ/2-2,θ/2+2)时,以0.2°间隔计算G3图像在水平方向的投影值SR(β,x)。
步骤7:依据统计信息值,计算指针的两条边线,并计算边线交点T。该过程主要基于统计信息SL,SR及间距d等参数,计算统计信息最大值,确定指针边线。具体计算过程如下:
7.1确定d的变化范围:该值与半径设置和指针形状有关,一般由用户设定,可设定在(-4,+4)之间。
7.2计算统计信息值:任选一个间距d值,此时,该值可以由公式(4)计算,采用逐一比较方法,确定取得的局部最大值Fmax(d)时(α,β,x)应满足的条件。
F(d)=SL(α,x)+SR(β,x+d) (4)
7.3在d的取值范围内,逐一比较所有的局部最大值,选择取得全局最大值时的d值和(α,β,x)值,并构成四元组,设为(α0,β0,d0,x0)。
7.4边缘线计算:
在图像G3中,左侧边缘线通过点XL(x0-d0/2,nH)和PL,其中,PL可以由公式(5)获取:
PL=(x0-d0/2-nH×tanα0,0) (5)
同理,右侧边缘线通过点XR(x0+d0/2,nH)和PR,其中,PR可以由公式(6)获取:
PR=(x0+d0/2+nH×tanβ0,0) (6)
参照公式(1)和(2),可计算G3图像空间的直线XLPL,XRPR在G1图像空间的对应直线XL'PL'和XR'PR',并计算他们在G3空间的交点T。
步骤8:针对图像G1,围绕仪表中心C,对仪表刻度所在的环形区域进行投影映射,生成图像G4。具体过程是:对图像G1的以C为中心,半径在(r3,r4)环形区域,进行投影变换,生成图像G4,如图2所示。
G4的大小为(nw,nh),其中,nw=2πr4,nh=r4-r3,
步骤9:对图像G4进行二值化处理,并通过区域标记算法定位每个刻度的中心区域,确定每个刻度坐标,并映射回原始图像G1空间。
具体过程如下:
9.1二值化处理:先对图像G4采用二值化处理,阈值可以由用户设定或根据自适应算法确定,该阈值在一般在80-160之间。
9.2计算刻度的分布:采用区域标记算法,计算所有的连通黑色区域,确定每个刻度的中心位置Mi'。
9.3计算刻度的原始坐标:参照公式(1,2),计算G4空间点Mi'在原始图像G2中的坐标Mi。
步骤10:在图像G1坐标空间中,过圆心C点,计算指针针尖T与邻近2个刻度之间夹角,并通过线性插值计算指针所指的刻度值。
具体过程如下:
10.1对所有的刻度坐标Mi,计算与指针针尖位置T最近距离的2个点M1,M2。
10.2基于圆心C,计算指针针尖T和M1和M2的夹角关系,如图3所示,Φ1为夹角∠TCM1,Φ2为夹角∠TCM2。
10.3采用插值方法计算当前的刻度值。
VT=(VM1×Φ2+VM2×Φ1)/(Φ1+Φ2) (7)
其中,VM1和VM2分别为点M1和M2对应的刻度值,该值由用户预先设定;VT为指针对应的刻度值。
本发明技术方案主要采用图像滤波、特征拟合、环形投影和线性插值等方法,具有如下有益效果:
1.在指针计算方面,基于特征拟合的仪表指针定位技术,通过指针边缘特征的环形投影、边缘统计和拟合分析,计算指针两条夹边及其交点。
2.在刻度计算方面,采用仪表指针与邻近刻度的线性插值的方法,能准确计算当前刻度值,避免了因仪表中心定位偏差或图像倾斜而导致的刻度值计算偏差。
3.在算法鲁棒性方面,通过抓取核心特征、边缘统计和线性插值等方法,算法有效避免各类干扰因素的影响,算法稳定,计算结果准确。
附图说明
图1为指针环形投影示意图;
图2为刻度投影计算示意图;
图3为刻度值计算示意图;
图4为本发明算法流程图;
具体实施方式
下面结合附图对本发明做进一步说明。
如图1、图2、图3和图4所示,一种基于特征拟合的变电站仪表设备巡检的视觉检测方法,具体包括如下步骤,参看图4:
步骤1.采集一张变电站巡检仪表图像,具体通过CCD相机拍摄目标图像并将图像转换成灰度图像G;
步骤2.针对图像G进行预处理,去除噪声点,生成图像G1。
步骤3.利用边缘算子提取G1的边缘信息,生成图像G2。
步骤4.采用Hough算法分析图像G2中圆形目标,计算出仪表中心C。
此阶段为仪表指针计算:
步骤5.围绕仪表中心C,对仪表指针所在的环形区域进行投影映射,生成图像G3,如图1所示。
步骤6.在图像G3中,计算符合指针边缘特征的统计信息值。
步骤7:依据统计信息值,计算指针的两条边线,并计算边线交点T。
此阶段为刻度计算:
步骤8:针对图像G1,围绕仪表中心C,对仪表刻度所在的环形区域进行投影映射,生成图像G4,如图2所示。
步骤9:对图像G4进行二值化处理,并通过区域标记算法定位每个刻度的中心区域,确定每个刻度坐标,并映射回原始图像G1空间。
此阶段定位刻度值:
步骤10:在图像G1坐标空间中,过圆心C点,计算指针针尖T与邻近2个刻度之间夹角,并通过线性插值计算指针所指的刻度值,如图3所示。
步骤1所述的采集一张变电站巡检仪表的图像,具体过程如下:通过工业相机拍摄目标图像,并转换成灰度图像G。转换具体如下:
Gray=R×0.3+G×0.59+B×0.11
所述的对目标图像中的每一个像素点处理,就能将其转换成灰度图像G;其中R、G、B为一个像素点的值的基色分量;
步骤2所述的仪表图像预处理,主要是去除噪声点,一般通过中值滤波方法除去灰度图像G的噪声点,生成图像G1。具体过程如下:
2-1、选定5×5的邻域,将邻域内的像素按灰度等级进行排序,确定其中值;
2-2、将确定的中值赋予滤波后图像中的相应像素点,所述的相应像素点即进行中值滤波的那个像素点;
步骤3所述的计算G1的边缘信息,一般采用Canny或Sobel算子提取边缘信息,生成图像G2。该边缘信息主要用于定位仪表中心和指针位置。
步骤4.采用Hough算法提取图像G2中的圆形目标,并定位圆形中心C。
步骤5.围绕仪表中心C,对仪表指针所在的环形区域进行投影映射,生成图像G3。考虑仪表的形状特征,一般的环形区域可以选择在半径为(r1,r2)区域,如图1所示。
其中,r2的长度一般近似为指针针尖T到仪表中心C的距离。r1可以选择为r2的1/3到2/3长度。
G3的大小为(nW,nH),其中,nW=2πr2,nH=r2-r1,此时,矩形上任一点P(x,y)可以映射到仪表上的点P'(x',y'),此时,
P'距离中心点的距离为r1+y,对应弧度为x/r2;
x'=(r1+y)×cos(x/r2) (1)
y'=(r1+y)×sin(x/r2) (2)
步骤6.在图像G3中,计算符合指针边缘特征的统计信息值。其中,指针边缘特征包括指针的左边线和右边线。该特征在投影视图G3中,可以用特征三元组(α,β,d)来描述,如图2所示。考虑到投影轮廓值与X坐标相关,因此统计信息值F可以由公式(3)表示:
F=f((α,β,d),x) (3)
该公式表示为以位置x处的垂线为基准,左侧边线倾角为α,右侧倾角为β,且两侧边线在底部离垂线间距均为d/2,如图1底部的放大区域所示。此时,G3中所包含的边缘点数即为统计信息值。具体计算过程如下:
6.1确定α和β搜索范围:考虑到指针夹角一般相对固定,但考虑工艺差异和拍照角度差异性,夹角可能会有一定的变化,但变化幅度较小。假设仪表指针夹角为θ(该值由用户输入,一般在3-6°),则α和β搜索的范围一般在(θ/2-2,θ/2+2)即可。
6.2左侧边缘信息统计:以G3图像底部为基准,当倾角为(θ/2-2,θ/2+2)时,以0.2°间隔计算G3图像在水平方向的投影值SL(α,x)。
6.3右侧边缘信息统计:以G3图像底部为基准,当倾角为(θ/2-2,θ/2+2)时,以0.2°间隔计算G3图像在水平方向的投影值SR(β,x)。
步骤7:依据统计信息值,计算指针的两条边线,并计算边线交点T。该过程主要基于统计信息SL,SR及间距d等参数,计算统计信息最大值,确定指针边线。具体计算过程如下:
7.1确定d的变化范围:该值与半径设置和指针形状有关,一般由用户设定,可设定在(-4,+4)之间。
7.2计算统计信息值:任选一个间距d值,此时,该值可以由公式(4)计算,采用逐一比较方法,确定取得的局部最大值Fmax(d)时(α,β,x)应满足的条件。
F(d)=SL(α,x)+SR(β,x+d) (4)
7.3在d的取值范围内,逐一比较所有的局部最大值,选择取得全局最大值时的d值和(α,β,x)值,并构成四元组,设为(α0,β0,d0,x0)。
7.4边缘线计算:
在图像G3中,左侧边缘线通过点XL(x0-d0/2,nH)和PL,其中,PL可以由公式(5)获取:
PL=(x0-d0/2-nH×tanα0,0) (5)
同理,右侧边缘线通过点XR(x0+d0/2,nH)和PR,其中,PR可以由公式(6)获取:
PR=(x0+d0/2+nH×tanβ0,0) (6)
参照公式(1)和(2),可计算G3图像空间的直线XLPL,XRPR在G1图像空间的对应直线XL'PL'和XR'PR',并计算他们在G3空间的交点T。
步骤8:针对图像G1,围绕仪表中心C,对仪表刻度所在的环形区域进行投影映射,生成图像G4。具体过程是:对图像G1的以C为中心,半径在(r3,r4)环形区域,进行投影变换,生成图像G4,如图2所示。
G4的大小为(nw,nh),其中,nw=2πr4,nh=r4-r3,
步骤9:对图像G4进行二值化处理,并通过区域标记算法定位每个刻度的中心区域,确定每个刻度坐标,并映射回原始图像G1空间。
具体过程如下:
9.1二值化处理:先对图像G4采用二值化处理,阈值可以由用户设定或根据自适应算法确定,该阈值在一般在80-160之间。
9.2计算刻度的分布:采用区域标记算法,计算所有的连通黑色区域,确定每个刻度的中心位置Mi'。
9.3计算刻度的原始坐标:参照公式(1,2),计算G4空间点Mi'在原始图像G2中的坐标Mi。
步骤10:在图像G1坐标空间中,过圆心C点,计算指针针尖T与邻近2个刻度之间夹角,并通过线性插值计算指针所指的刻度值。
具体过程如下:
10.1对所有的刻度坐标Mi,计算与指针针尖位置T最近距离的2个点M1,M2。
10.2基于圆心C,计算指针针尖T和M1和M2的夹角关系,如图3所示,Φ1为夹角∠TCM1,Φ2为夹角∠TCM2。
10.3采用插值方法计算当前的刻度值。
VT=(VM1×Φ2+VM2×Φ1)/(Φ1+Φ2) (7)
其中,VM1和VM2分别为点M1和M2对应的刻度值,该值由用户预先设定;VT为指针对应的刻度值。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案基础之上,本领域技术人员无需付出创造性劳动即可做出各种调整或变换方案仍然在本发明的保护范围以内。
Claims (8)
1.一种基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,包括如下步骤:
步骤1.采集变电站巡检仪表图像,并将图像转换成灰度图像G;
步骤2.对图像G进行预处理,去除噪声点,生成图像G1;
步骤3.利用边缘算子提取G1的边缘信息,生成图像G2;
步骤4.采用Hough算法分析图像G2中圆形目标,计算出仪表中心C;
步骤5.围绕仪表中心C,对仪表指针所在的环形区域进行投影映射,生成图像G3;其中,环形区域选择在半径为(r1,r2)区域,r2的长度为指针针尖T到仪表中心C的距离;r1长度为r2的1/3到2/3之间;
G3的大小为(nW,nH)的矩形,其中,nW=2πr2,nH=r2-r1,此时,矩形上任一点P(x,y)可以映射到仪表上的点P'(x',y'),
x'=(r1+y)×cos(x/r2) (1)
y'=(r1+y)×sin(x/r2) (2)
此时,P'距离中心点的距离为r1+y,对应弧度为x/r2;
步骤6.在图像G3中,计算符合指针边缘特征的统计信息值;其中,指针边缘特征包括指针的左边线和右边线;该特征在投影视图像G3中,用特征三元组(α,β,d)来描述,统计信息值F由以下公式表示:
F=f((α,β,d),x); (3)
该公式表示为以位置x处的垂线为基准,左侧边线倾角为α,右侧边线倾角为β,且两侧边线在底部离垂线间距均为d/2;此时,G3中所包含的边缘点数即为统计信息值;具体计算过程如下:
步骤6.1:确定α和β搜索范围:仪表指针夹角为θ,则α和β搜索的范围为(θ/2-2,θ/2+2);
步骤6.2:左侧边缘信息统计SL:以G3图像底部为基准,当左侧倾角为(θ/2-2,θ/2+2)时,以0.2°间隔计算G3图像在水平方向的投影值SL(α,x);
步骤6.3:右侧边缘信息统计SR:以G3图像底部为基准,当右侧倾角为(θ/2-2,θ/2+2)时,以0.2°间隔计算G3图像在水平方向的投影值SR(β,x);
步骤7:依据步骤6得到的统计信息值,计算指针的两条边线,并计算边线交点T;
步骤8:针对图像G1,围绕仪表中心C,对仪表刻度所在的环形区域进行投影映射,生成图像G4;
步骤9:对图像G4进行二值化处理,并通过区域标记算法定位每个刻度的中心区域,确定每个刻度坐标,并映射回原始图像G1空间;
步骤10:在图像G1坐标空间中,过圆心C点,计算指针针尖T与邻近2个刻度之间夹角,并通过线性插值计算指针所指的刻度值。
2.如权利要求1所述的基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,在所述步骤1中,进一步包括以下步骤:
步骤1.1:通过工业相机拍摄目标图像:
步骤1.2:对目标图像中的每一个像素点处理并将其转换成灰度图像G,转换公式为:Gray=R×0.3+G×0.59+B×0.11;其中,R、G、B为一个像素点的值的基色分量。
3.如权利要求1所述的基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,在所述步骤2中,通过中值滤波方法除去灰度图像G的噪声点,生成图像G1;进一步包括以下步骤:
步骤2.1:选定5×5的邻域,将邻域内的像素按灰度等级进行排序,确定其中值;
步骤2.2:将确定的中值赋予滤波后图像中的相应像素点,所述的相应像素点为进行中值滤波的那个像素点。
4.如权利要求1所述的基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,在所述步骤3中,采用Canny或Sobel算子提取边缘信息,该边缘信息用于定位仪表中心和指针位置。
5.如权利要求1所述的基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,在所述步骤7中,基于左侧边缘信息统计SL、右侧边缘信息统计SR及间距d,计算统计信息最大值,确定指针边线;具体计算过程如下:
步骤7.1:确定d的变化范围:该值与半径设置和指针形状有关,由用户设定,可设定在(-4,+4)之间;
步骤7.2:计算统计信息值:任选一个间距d值,此时,该值由以下公式计算,采用逐一比较方法,确定取得的局部最大值Fmax(d)时(α,β,x)应满足的条件;
F(d)=SL(α,x)+SR(β,x+d); (4)
步骤7.3:在d的取值范围内,逐一比较所有的局部最大值,选择取得全局最大值时的d值和(α,β,x)值,并构成四元组,设为(α0,β0,d0,x0);
步骤7.4:边缘线计算:
在图像G3中,左侧边缘线通过点XL(x0-d0/2,nH)和PL,其中,PL可以由以下公式获取:
PL=(x0-d0/2-nH×tanα0,0); (5)
同理,右侧边缘线通过点XR(x0+d0/2,nH)和PR,其中,PR可以由以下公式获取:
PR=(x0+d0/2+nH×tanβ0,0); (6)
由公式(1)和(2),计算G3图像空间的直线XLPL,XRPR在G1图像空间的对应直线XL'PL'和XR'PR',并计算他们在G3空间的交点T。
6.如权利要求1所述的基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,在所述步骤8中,具体过程是:对图像G1的以C为中心,半径在(r3,r4)环形区域,进行投影变换,生成图像G4,则
G4的大小为(nw,nh),其中,nw=2πr4,nh=r4-r3。
7.如权利要求1所述的基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,在所述步骤9中,具体过程如下:
步骤9.1:二值化处理:先对图像G4采用二值化处理,阈值由用户设定或根据自适应算法确定;
步骤9.2:计算刻度的分布:采用区域标记算法,计算所有的连通黑色区域,确定每个刻度的中心位置Mi';
步骤9.3:计算刻度的原始坐标:根据公式(1)和(2),计算G4空间点Mi'在原始图像G2中的坐标Mi。
8.如权利要求1所述的基于特征拟合的变电站仪表设备巡检的视觉检测方法,其特征在于,在所述步骤10中,具体过程如下:
步骤10.1:对所有的刻度坐标Mi,计算与指针针尖位置T最近距离的2个点M1,M2;
步骤10.2:基于圆心C,计算指针针尖T和M1和M2的夹角关系,其中,Φ1为夹角∠TCM1,Φ2为夹角∠TCM2;
步骤10.3:采用插值方法计算当前的刻度值;
VT=(VM1×Φ2+VM2×Φ1)/(Φ1+Φ2) (7)
其中,VM1和VM2分别为点M1和M2对应的刻度值,该值由用户预先设定;VT为指针对应的刻度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810129335.1A CN108509845B (zh) | 2018-02-08 | 2018-02-08 | 基于特征拟合的变电站仪表设备巡检的视觉检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810129335.1A CN108509845B (zh) | 2018-02-08 | 2018-02-08 | 基于特征拟合的变电站仪表设备巡检的视觉检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108509845A CN108509845A (zh) | 2018-09-07 |
CN108509845B true CN108509845B (zh) | 2020-06-23 |
Family
ID=63374547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810129335.1A Active CN108509845B (zh) | 2018-02-08 | 2018-02-08 | 基于特征拟合的变电站仪表设备巡检的视觉检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108509845B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6803940B2 (ja) * | 2019-03-26 | 2020-12-23 | 株式会社フュージョンテク | 遠隔メータ読取コンピュータ、その方法及びプログラム |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1693852A (zh) * | 2005-01-28 | 2005-11-09 | 华南理工大学 | 一种仪表指针自动检测识别方法及自动读数方法 |
JP3799408B1 (ja) * | 2005-10-18 | 2006-07-19 | 国立大学法人山口大学 | 画像処理装置および画像処理方法 |
CN103759758B (zh) * | 2014-01-26 | 2016-02-17 | 哈尔滨工业大学 | 一种基于机械角度和刻度识别的汽车仪表指针的位置检测方法 |
-
2018
- 2018-02-08 CN CN201810129335.1A patent/CN108509845B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108509845A (zh) | 2018-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104331876B (zh) | 直线检测、图像处理的方法及相关装置 | |
CN110807355B (zh) | 一种基于移动机器人的指针仪表检测与读数识别方法 | |
JP5699788B2 (ja) | スクリーン領域検知方法及びシステム | |
CN103759758B (zh) | 一种基于机械角度和刻度识别的汽车仪表指针的位置检测方法 | |
CN102176228B (zh) | 一种用于识别多个指针式仪表表盘信息的机器视觉方法 | |
CN109558871B (zh) | 一种指针式仪表读数识别方法及装置 | |
CN110910350B (zh) | 一种用于风电塔筒的螺母松动检测方法 | |
CN108007388A (zh) | 一种基于机器视觉的转盘角度高精度在线测量方法 | |
CN105303168A (zh) | 一种多视角的指针式仪表识别方法和装置 | |
CN109447061A (zh) | 基于巡检机器人的电抗器油位计识别方法 | |
JP2006120133A (ja) | アナログメータの自動読取方法および装置およびプログラム | |
CN117274722B (zh) | 基于红外图像的配电箱智能检测方法 | |
JP2009259036A (ja) | 画像処理装置、画像処理方法、画像処理プログラム、記録媒体、及び画像処理システム | |
CN109389165A (zh) | 基于巡检机器人的变压器油位计识别方法 | |
CN111354047B (zh) | 一种基于计算机视觉的摄像模组定位方法及系统 | |
CN108960236B (zh) | 基于二维码匹配与图像识别的指针式仪表识别方法 | |
CN111476246A (zh) | 应用于复杂环境下指针仪表鲁棒高效智读方法 | |
CN114005108A (zh) | 一种基于坐标变换的指针式仪表度数识别方法 | |
CN113469178A (zh) | 一种基于深度学习的电力表计识别方法 | |
CN108509845B (zh) | 基于特征拟合的变电站仪表设备巡检的视觉检测方法 | |
CN113408519B (zh) | 一种基于模板旋转匹配的指针仪表读数的方法及系统 | |
CN107767366B (zh) | 一种输电线路拟合方法及装置 | |
CN113298725A (zh) | 一种轮船图标图像叠加误差的纠正方法 | |
CN109359637B (zh) | 基于机器视觉的指针式仪表读值方法 | |
CN116880353A (zh) | 一种基于两点间隙的机床对刀方法 |
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 |