CN109931884B - 一种长条形水嘴旋转角度非接触测量方法 - Google Patents
一种长条形水嘴旋转角度非接触测量方法 Download PDFInfo
- Publication number
- CN109931884B CN109931884B CN201910101024.9A CN201910101024A CN109931884B CN 109931884 B CN109931884 B CN 109931884B CN 201910101024 A CN201910101024 A CN 201910101024A CN 109931884 B CN109931884 B CN 109931884B
- Authority
- CN
- China
- Prior art keywords
- water nozzle
- straight line
- max
- image
- angle
- 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
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种长条形水嘴旋转角度非接触测量方法,用一种基于累计概率霍夫变换(PPHT)和加权最小二乘(WLS)的算法测量长条形水嘴旋转角度。具体包括以下步骤:A、CCD读取水嘴的图像;B、中值滤波去除噪点;C、将水嘴图像从原本的RGB空间转换到HSV空间;D、对水嘴图像进行形态学运算;E、进行Canny边缘检测;F、利用PPHT—WLS算法对水嘴图像进行分析,进行水嘴角度检测并输出角度信息。
Description
技术领域
本发明涉及一种长条形水嘴旋转角度非接触测量方法。
背景技术
利用机器视觉可以对长条形水嘴进行非接触测量,现在提出一种基于累计概率霍夫变换(PPHT)和加权最小二乘(WLS)的检测法测量水嘴旋转角度。作为传统线检测方法的霍夫变换(HT)和最小二乘(LS)方法广泛用于机器视觉,标准Hough变换(SHT)作为一种全局最优算法,具有良好的稳定性和良好的鲁棒性。然而,当要求高精度时,由于其详尽的搜索,SHT将被认为具有较慢的计算速度,存储空间的浪费和较差的实时能力。鉴于这些缺点,相关学者提出了一些改进的算法。众所周知的改进是随机Hough变换(RHT),已由Xu L等人提出。RHT利用多对一映射和随机样本选择以减轻计算和存储负荷。Kiryati等提供了霍夫变换的另一个有效版本,即概率霍夫变换(PHT),它维持SHT投票方案,但是使用相对较小的抽样。Matas J.等提出了渐进概率霍夫变换(PPHT),与PHT不同,PHT在预先选定的输入点部分上执行,PPHT通过移除相应行中的支撑点并减少投票计数来最小化检测线所需的计算量。在时间消耗方面,PPHT最快运行速度,非常适合实时线路检测。对于LS线拟合方法,考虑随机误差对线性精度的影响,获得高检测精度,可以弥补HT的精度不足。然而,当检测到的边缘受到噪声的影响,以及标识的边缘线有毛刺不在一条直线上时,LS的检测结果误差将显着增加,在这些情况下,可以使WLS进行拟合。
作为传统线检测方法的霍夫变换(HT)和最小二乘(LS)方法广泛用于机器视觉,标准Hough变换(SHT)作为一种全局最优算法,具有良好的稳定性和良好的鲁棒性。本发明提出了一种基于累计概率霍夫变换(PPHT)和加权最小二乘(WLS)的标识线检测法测量水嘴旋转角度。
发明内容
本发明提供一种长条形水嘴旋转角度非接触测量方法,用一种基于累计概率霍夫变换(PPHT)和加权最小二乘(WLS)的算法测量长条形水嘴旋转角度。
为了实现上述目的,本发明的技术方案如下:一种长条形水嘴旋转角度非接触测量方法,其特征在于所述方法包括以下步骤:
A、摄像头读取水嘴图像;
B、中值滤波去除噪点;
C、将水嘴图像从原本的RGB空间转换到HSV空间;
D、对水嘴图像进行形态学运算;
E、进行Canny边缘检测;
F、利用PPHT—WLS算法对水嘴图像进行分析,进行水嘴角度检测并输出角度信息。
进一步地,步骤A中,保证摄像头和水嘴阀芯垂直安装,且水嘴手柄完全在画面内,这样才能保证采得的水嘴画面符合处理要求。
进一步地,步骤F又包括:
(1)PPHT线检测,计算最长边;
(2)判断边缘点到最长边的距离,若距离小于阈值,则保留该边缘点,通过WLS拟合直线,得到符合误差阈值初检测的最长直线;
(3)计算新拟合的水嘴边缘线的角度,即为当前水嘴的角度。
其中步骤F具体为:设采集到的边缘检测像素点点集为M=(xi,yi)T,M中的数据点分布在最长直线附近,设定误差阈值为d,Hough变换公式为:
ρ=xicosθ+yisinθ(i=1,2,...,n)
ρ为从原点到直线的垂直距离,θ为从x轴算起的角度;
根据Hough变换公式,对M作累计概率Hough变换,并且求出最长直线拟合直线参数(θmax,ρmax)。然后找出拟合直线(θmax,ρmax)附近的数据点集M*。将最长直线(定义)方程改写为斜截式公式:
ymax=amaxx+bmax
计算M中的点到斜截式直线的距离di
若有
di<d
则
(xi,yi)∈M*(xj,yj),(j=1,2,...,n*;n*≤n)
M*为符合误差阈值初检测的最长直线斜截式方程表示方法附近的点集,n*为数据集中的点数;
该新拟合的直线即为拟合后的水嘴边缘线;
进而可求出新拟合的直线的角度变化为当前水嘴的角度θ:
本发明利用机器视觉技术可以对长条形水嘴旋转角度进行非接触测量,运用的PPHT-WLS算法与PPHT和SHT算法相比,具有更高的检测精度,同时也减少了计算时间。
附图说明
图1为本发明的流程图。
图2为Hough变换示意图。
下面结合附图和实施例对本发明作详细说明。
具体实施方式
如图所示,本发明的一种长条形水嘴旋转角度非接触测量方法,其特征在于所述方法包括以下步骤:
A、CCD摄像头读取水嘴的图像;
B、中值滤波去除噪点;
C、将水嘴图像从原本的RGB空间转换到HSV空间;
D、对水嘴图像进行形态学运算;
E、进行Canny边缘检测;
F、利用PPHT—WLS算法对水嘴图像进行分析,进行水嘴角度检测并输出角度信息。
下面对本发明的算法实现进行详细说明。
1、中值滤波
由于CCD元件本身存在缺陷以及拍摄现场噪声干扰,摄像头采集的图像中存在有噪声。噪声的存在会使Canny检测处理时产生无用的边缘信息,增加处理时间且影响检测精度,因而应首先对图像进行去噪声处理。
中值滤波是一种统计排序滤波器,其表达式为:
g(i,j)=median{f(i-k,j-l)},(k,l∈W)
其中f(i,j),g(i,j)分别为原始图像和处理后的图像,W为二维模板,通常为3*3,5*5,7*7,9*9模板。中值滤波的步骤为:
(1)将滤波模板在图像中滑动,并将模板中心与图像中的某个像素重合;
(2)读取模板中各个对应像素的像素值;
(3)将这些像素值从小到大排列;
(4)在已经进行的排序中,取大小位于最中间的像素值,并将其赋给对应模板中心位置的像素。
2、RGB空间转换到HSV空间
RGB彩色空间中三个分量有着很高的相关性,常用的数字图像处理算法中,若直接对RGB模型中的R、G、B分量进行处理,其处理过程中很可能会引起三个分量不同程度的变化,引起处理图像产生色差甚至失真,不利于后续的图像处理。本文将已经滤波的图像转到HSV颜色空间进行处理,HSV颜色空间具有两个特征,一方面亮度分量与彩色信息无关,另一方面色调和饱和度与人类所感知的颜色一一对应,这些特征使得在HSV颜色空间下处理图像所受光照的影响较小,同时精确反映图片的颜色信息。RGB颜色模型转换到HSV颜色模型的公式如下:
V=max
其中max为图像中R(红),G(绿),B(蓝)三个通道存储值中的最大值,min为最小值。RGB彩色空间转换到HSV彩色空间实则是对图像三个通道中存储的值进行一一变换,再将变换后的数值存入三个通道,此时三个通道的存储量依次为H(亮度)、S(色调)、V(饱和度)的存储量。
3、形态学处理
通常由于噪声的影响,水嘴形状在阈值得到的边界往往是很不平滑的,内部区域可能会出现一些噪声孔,背景区域上也可能随机散布着一些小的噪声孔。本文对图像采用连续的开运算和闭运算处理,改善这种情况,并且降低后续图像检测识别算法的复杂度。
本文对水嘴图像采用连续的开运算和闭运算处理,改善这种情况。
膨胀或者腐蚀操作就是将图像(或图像的一部分区域,我们称之为A)与核(我们称之为B)进行卷积,核就是进行卷积的模板。膨胀就是求局部最大值的操作,核B与图形卷积,即计算核B覆盖的区域的像素点的像素最大值,并把这个像素值最大值赋值给参考点指定的像素。腐蚀就是局部最小值的操作,核B与图形卷积,即计算核B覆盖的区域的像素点的最大值,并把这个像素值的最大值赋值给参考点指定的像素。相反闭运算是对图像先进行膨胀操作再进行腐蚀操作。
腐蚀操作用公式可表示为:
膨胀操作用公式可表示为:
开运算是先进行腐蚀再进行膨胀操作。用公式表达为:
闭运算是对图像先进行膨胀操作再进行腐蚀操作。用公式表达为:
4、Canny边缘检测
Canny算法检测性能良好,既保留重要边缘信息,又能有效抑制伪边缘,定位精度高。具体的算法如下:
(1)高斯滤波。为减小数字图像中噪声点对图像检测质量的影响,在检测环节首先对图像进行滤波处理;
(2)计算梯度幅值和方向角。利用导数算子找到图像灰度的沿着两个方向的偏导数,并求出梯度的幅值和梯度方向;
(3)非极大值抑制。遍历整个图像,使用3×3大小,包含8个方向的邻域梯度幅值阵列的所有像素沿梯度方向进行梯度幅值的插值。在每一个像素点上,比较邻域的中心像素与沿梯度方向的两个梯度幅值的插值结果的大小,找出局部极大值点,并保留幅值局部极大值点,删除局部非极大值点。
(4)双阈值法检测并连接边缘。经过非极大值处理的图像,还需使用高低阈值进行分割,得到2个阈值边缘图像。
5、霍夫变换(SHT)原理
霍夫变换是数字图象技术中一种有效的发现直线的算法。它是先把直角坐标系的目标点映射到极坐标系上进行累积,即先使直角坐标系平面上任一直线上的所有点均累积到极坐标系的同一点集中去,然后通过寻找极坐标系中点集的峰值,来发现长的直线特征。由于这种点集是通过累积统计得到的,因而能够容忍直线的间断。
霍夫变换的基本思想是点—线的对偶性,在图像空间XY里,所有过点(x,y)的直线都满足方程
y=kx+b
其中k为该直线的斜率,b为该直线在y轴上的截距,如果用参数表示则为:
ρ=xcosθ+ysinθ
其中ρ为从原点到直线的垂直距离,θ为从x轴算起的角度,这条直线在ρ-θ平面中为一点,而在x-y平面中通过一点的一簇直线变换到ρ-θ平面时,将形成一条类似正弦曲线的轨迹,也即在x-y平面上一个点对应ρ-θ平面上一条曲线。若在x-y平面上有三个共线点,它们变换到ρ-θ平面上为有一公共交点的三条曲线,公共交点的ρ-θ参数就是共线直线的参数。根据这个原理,如图2所示,可以用Hough变换抽取直线,通常将x-y平面称为图像平面,将ρ-θ平面称为参数平面。
Hough变换算法描述如下:
(1)在ρ、θ合适的最大值与最小值之间建立一个离散的参数空间;
(2)建立一个累加器A(ρ,θ),并将每一个元素置为0。
(3)对于x-y坐标系上的每一个像素点(x,y),对每一个θ的取值θ',计算ρ'=xcosθ'+ysinθ',并在相应累加器加1:A(ρ',θ')=A(ρ',θ')+1(即把矩阵A的元素A(ρ',θ')累加)
(4)对数组A(ρ,θ)进行局部峰值检测,得出被检测直线的参数ρ和θ。
由霍夫变换的原理及算法步骤可知,霍夫变换计算时间与参数离散化的区间分辨率有关,区间划分越细致,精确度越高,但遍历次数越多,时间越慢,同时计数器数量也越多,计算开销越多。并且数字图像中,像素点坐标值都以二维整数形式存在,过度细分的参数空间有可能会因为计数器无法获得足够的计数而造成检测失败。
6、累计概率霍夫变换(PPHT)原理
累计概率霍夫变换是标准霍夫变换的一种改进算法。本文使用了累计概率霍夫变换(PPHT)来对图像进行线检测。PPHT算法中的映射和直线检测是交替进行的,由于采用了随机取点映射的方法,最明显的直线特征有望被最早检测到,算法在运行一定时间后可停止运行并仍然输出有效的结果;其次,PPHT算法的运算过程也可在所有待处理点完成向参数空间的映射或被归类到某一直线中后停止。这种情况下,只有一小部分待处理点完成了映射,其余的点作为被检测到的直线上的点从待处理点集中去除而不必进行向参数空间的映射,减少了算法的运算开销。累计概率霍夫变换相比标准霍夫变换效率上有了极大提升,更适合于实时场景中的直线检测。
7、加权最小二乘法
一般最小二乘法将时间序列中的各项数据的重要性同等看待,而事实上在存在噪声或者贴条边缘有毛刺和凸起等引起的异常值时,比较合理的方法就是使用加权的方法,对比较近的点赋以较大的权数,对比较远的直线则赋以较小的权数。
假设在图像空间的特征区域中存在n个像素点,由W={(xi,yi),i=1,2,...,n}表示。配有像素点的线可以以任何方向出现。
设直线模型为
则对于给定数据点(xi,yi),其残差平方和为
其中wi表示权重
对Q求导
可解得
令
代入上式可解出
其中权重wi必须基于距离δ,本文选用Fair权函数作为权数,即
WLS方法应用于长条形水嘴边缘的线检测的精确拟合,仅在3个像素距离内设置的特征点。从权重公式中可以看出,权重随着距离的增加而减小,减小了异常值对直线拟合的影响。
8、累计概率霍夫变换(PPHT)和加权最小二乘(WLS)结合
PPHT的作用是找到图像中的所有直线,PPHT找到的图像中所有直线后,计算每一条直线的长度,长度最长的那条就是最长直线,对于长条形水嘴来说最长直线也即水嘴的边缘线。在通过PPHT确认水嘴边缘的近似位置之后,需要提取接近边缘的特征点以用于WLS的精确线性拟合。
设采集到的边缘检测像素点点集为M=(xi,yi)T,(i=1,2,3...n,n为数据集中的数据点数),其中T表示转置的意思。M中的数据点分布在最长直线附近,设定误差阈值为d。Hough变换公式为:
ρ=xi cosθ+yi sinθ(i=1,2,...,n)
根据上式,对M作累计概率Hough变换,并且求出最长直线拟合直线参数(θmax,ρmax)。然后找出拟合直线(θmax,ρmax)附近的数据点集M*。将最长直线(定义)方程改写为斜截式公式:
ymax=amaxx+bmax
计算M中的点到斜截式直线的距离
若有
di<d
则
(xi,yi)∈M*(xj,yj),(j=1,2,...,n*,n*≤n)
M*为符合误差阈值初检测的最长直线斜截式方程表示方法附近的点集,n*为数据集中的点数。最后以M*为拟合数据,可得到新的直线方程的参数以((xj)min,yj))∈M和((xj)max,yj))∈M,可确定新拟合的线段区间:
进而可求出新拟合的直线的角度变化为当前的角度θ:
Claims (2)
1.一种长条形水嘴旋转角度非接触测量方法,其特征在于所述方法包括以下步骤:
A、摄像头读取水嘴的图像;
B、中值滤波去除噪点;
C、将水嘴图像从原本的RGB空间转换到HSV空间;
D、对水嘴图像进行形态学运算;
E、进行Canny边缘检测;
F、利用PPHT—WLS算法对水嘴图像进行分析,进行水嘴角度检测并输出角度信息,其中PPHT为渐进概率霍夫变换,WLS为加权最小二乘;
步骤F又包括:
(1)PPHT线检测,计算最长边;
(2)判断边缘点到最长边的距离,若距离小于阈值,则保留该边缘点,通过WLS拟合直线,得到符合误差阈值初检测的最长直线;
(3)计算新拟合的水嘴边缘线的角度,即为当前水嘴的角度;
步骤F具体为:设采集到的边缘检测像素点点集为M=(xi,yi)T,xi,yi分别为像素点的横坐标和纵坐标,M中的数据点分布在最长直线附近,设定误差阈值为d,Hough变换公式为:
ρ=xicosθ+yisinθ,i=1,2,...,n n为正整数
ρ为从原点到直线的垂直距离,θ为从x轴算起的角度;
根据Hough变换公式,对M作累计概率Hough变换,并且求出最长直线拟合直线参数(θmax,ρmax),然后找出拟合直线(θmax,ρmax)附近的数据点集M*;将最长直线方程改写为斜截式公式:
ymax=amaxx+bmax
计算M中的点到斜截式直线的距离
若有
di<d
则
(xi,yi)∈M*(xj,yj),j=1,2,...,n*;n*≤n n*为正整数,n为正整数
M*为符合误差阈值初检测的最长直线斜截式方程表示方法附近的点集,n*为数据集中的点数;
最后以M*为拟合数据,可得到新的直线方程的参数以((xj)min,yj))∈M和((xj)max,yj))∈M,可确定新拟合的线段区间,其中(xj)min,(xj)max分别为阈值比对后得到的像素点的最大横坐标和最小横坐标:
该新拟合的直线即为拟合后的水嘴边缘线;
进而可求出新拟合的直线的角度变化为当前水嘴的角度θ:
2.如权利要求1所述的长条形水嘴旋转角度非接触测量方法,其特征在于,步骤A中,保证摄像头和水嘴阀芯垂直安装,且水嘴手柄完全在画面内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910101024.9A CN109931884B (zh) | 2019-01-31 | 2019-01-31 | 一种长条形水嘴旋转角度非接触测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910101024.9A CN109931884B (zh) | 2019-01-31 | 2019-01-31 | 一种长条形水嘴旋转角度非接触测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109931884A CN109931884A (zh) | 2019-06-25 |
CN109931884B true CN109931884B (zh) | 2021-06-04 |
Family
ID=66985426
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910101024.9A Active CN109931884B (zh) | 2019-01-31 | 2019-01-31 | 一种长条形水嘴旋转角度非接触测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109931884B (zh) |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3882891B2 (ja) * | 2001-11-05 | 2007-02-21 | 日本電気株式会社 | 電子部品検査装置、その方法およびプログラム |
CN103996197A (zh) * | 2014-05-30 | 2014-08-20 | 北京农业信息技术研究中心 | 性状采集方法及系统 |
CN104103254B (zh) * | 2014-08-12 | 2016-04-13 | 深圳市华星光电技术有限公司 | 提高wrgb色彩饱和度的方法 |
CN104279987A (zh) * | 2014-10-30 | 2015-01-14 | 盐城工学院 | 一种基于图像处理的零件角度检测方法 |
CN104616010A (zh) * | 2014-12-04 | 2015-05-13 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 一种基于Hough变换与改进型最小二乘法的直线检测方法 |
CN106157289B (zh) * | 2015-04-08 | 2019-08-09 | 株式会社理光 | 线检测方法和设备 |
CN105160309B (zh) * | 2015-08-24 | 2018-12-07 | 北京工业大学 | 基于图像形态学分割及区域生长的三车道检测方法 |
CN106934326B (zh) * | 2015-12-29 | 2020-07-07 | 同方威视技术股份有限公司 | 用于安全检查的方法、系统和设备 |
CN105700532B (zh) * | 2016-04-19 | 2018-08-10 | 长沙理工大学 | 基于视觉的变电站巡检机器人导航定位控制方法 |
CN106247969B (zh) * | 2016-09-21 | 2018-11-02 | 哈尔滨工业大学 | 一种基于机器视觉的工业磁芯元件的形变检测方法 |
CN107145890B (zh) * | 2017-05-02 | 2019-05-10 | 南通大学 | 一种远距离多视角环境下的指针式仪表盘自动读数方法 |
-
2019
- 2019-01-31 CN CN201910101024.9A patent/CN109931884B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109931884A (zh) | 2019-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112950508B (zh) | 一种基于计算机视觉的排水管道视频数据修复方法 | |
CN109859226B (zh) | 一种图形分割的棋盘格角点亚像素的检测方法 | |
CN110569857B (zh) | 一种基于质心距离计算的图像轮廓角点检测方法 | |
CN106709500B (zh) | 一种图像特征匹配的方法 | |
CN109559324B (zh) | 一种线阵图像中的目标轮廓检测方法 | |
CN107895375B (zh) | 基于视觉多特征的复杂道路线提取方法 | |
CN110414308B (zh) | 一种针对输电线路上动态异物的目标识别方法 | |
CN112200019B (zh) | 一种快速建筑夜景照明灯光故障检测方法 | |
CN111062293A (zh) | 一种基于深度学习的无人机森林火焰识别方法 | |
CN110619328A (zh) | 基于图像处理和深度学习的船舶水尺读数智能识别方法 | |
CN106815587A (zh) | 图像处理方法及装置 | |
CN114972575A (zh) | 一种基于轮廓边缘的直线拟合算法 | |
CN117496401A (zh) | 一种用于视频测量影像序列椭圆形目标点全自动识别与跟踪方法 | |
CN115880683B (zh) | 一种基于深度学习的城市内涝积水智能水位检测方法 | |
CN111241911B (zh) | 一种自适应的车道线检测方法 | |
CN109931884B (zh) | 一种长条形水嘴旋转角度非接触测量方法 | |
CN104573698A (zh) | 面向成像测量的图像纹理主方向识别方法 | |
US10115195B2 (en) | Method and apparatus for processing block to be processed of urine sediment image | |
Wu et al. | Research on crack detection algorithm of asphalt pavement | |
CN111695374B (zh) | 监控视角中斑马线区域的分割方法、系统、介质及设备 | |
CN108010076B (zh) | 一种面向密集工业棒材图像检测的端面外观建模方法 | |
US20210088456A1 (en) | Deposit detection device and deposit detection method | |
TWI384418B (zh) | 採用區域架構之影像處理方法及系統 | |
CN114511803A (zh) | 一种面向视觉跟踪任务的目标遮挡检测方法 | |
CN111784634A (zh) | 一种基于Harris-CPDA的角点检测方法 |
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 |