CN114581440A - 一种基于图像识别快速定位泄漏点的方法 - Google Patents

一种基于图像识别快速定位泄漏点的方法 Download PDF

Info

Publication number
CN114581440A
CN114581440A CN202210477646.3A CN202210477646A CN114581440A CN 114581440 A CN114581440 A CN 114581440A CN 202210477646 A CN202210477646 A CN 202210477646A CN 114581440 A CN114581440 A CN 114581440A
Authority
CN
China
Prior art keywords
image
pipeline
representing
points
point
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.)
Granted
Application number
CN202210477646.3A
Other languages
English (en)
Other versions
CN114581440B (zh
Inventor
江大白
汪刚
赵建标
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Applied Technology Co Ltd
Original Assignee
China Applied Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Applied Technology Co Ltd filed Critical China Applied Technology Co Ltd
Priority to CN202210477646.3A priority Critical patent/CN114581440B/zh
Publication of CN114581440A publication Critical patent/CN114581440A/zh
Application granted granted Critical
Publication of CN114581440B publication Critical patent/CN114581440B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • G06T7/001Industrial image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/22Matching criteria, e.g. proximity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/253Fusion techniques of extracted features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于图像识别快速定位泄漏点的方法,该方法通过输入原始管道的图像数据,利用图像分割算法对原始管道图像进行分割,生成图像的不同初始化区域,对不同初始化区域的原始管道图像数据,依次进行均值、样本方差、标准化处理;利用图像极值点的检测方法对管道图像进行特征提取,利用特征融合算法对特征提取的结果,把相同管道图像聚合而不同管道图像进行分离,进行图像的相似度判断定位出泄漏点的位置,该方法是在现有技术的基础上判断相邻区域图像的相似性迅速找到管道的泄漏点,操作简单可迅速对整条管道进行检测,在管道安装时在电脑系统匹配该方法中可以实现全部泄漏过程的自动化过程,节约了人力成本。

Description

一种基于图像识别快速定位泄漏点的方法
技术领域
本发明涉及管道泄漏检测技术领域,尤其涉及一种基于图像识别快速定位泄漏点的方法。
背景技术
在工业生产过程中,经常会使用管道进行气体或者液体的输送,因为管道输送不仅运输量大、连续、迅速、经济、安全、可靠、平稳以及投资少、占地少、费用低,且可实现自动控制,所以应用十分广泛。
在管道输送的过程中,气体或者液体长时间与管道壁面接触,会造成磨损和腐蚀,或者管道的接缝处焊接强度下降,从而导致管道出现泄露的情况。一旦出现了气体或者液体的泄露,一方面会造成原料的浪费,成本的增加;另一方面,尤其是易燃易爆的化工原料,一旦有管道泄漏发生,直接或间接的影响人们的生活和环境污染,更严重的会造成人员伤亡等。
针对存在的管道泄露问题,传统的做法是定时安排专职人员去巡检,并记录相关的巡检情况,发现泄露的点位,巡检人员还需要及时地进行汇报。但是,这会存在一些问题:一、巡检人员在巡检时需要携带多种工具,往往需要几个人进行配合,会增加较多的人力成本;二、由于管道往往跨度较长,安装位置比较高,在巡检过程中,巡检人员需要花费较多的时间,巡检率较低,且容易造成遗漏。
如何搭建高效、稳定的管道泄漏检测模型,通过采集的信号进行实时检测管道状况,从而使管道运输网系统能够快速、可靠地检测是否发生泄漏的方法是管道泄漏检测的一个研究热点。现有的研究多是利用负压波法,通过压力信号的变化和泄漏点产生的负压波形成的时间差,从而实现泄漏检测,该方法针对于小泄漏检测识别率低,效果达不到预期效果。在面对管道运输网系统复杂性、稳定性等问题时,无法高效、准确的实时检测管道发生泄漏的情况,还需要进一步深入的研究。
发明内容
为了克服现有技术存在的缺点与不足,本发明提供一种基于图像识别快速定位泄漏点的方法。
本发明所采用的技术方案是,该方法步骤为:
步骤S1:输入原始管道的图像数据;
步骤S2:利用图像分割算法对原始管道图像进行分割,生成图像的不同初始化区域
Figure DEST_PATH_IMAGE001
,初始化后默认所有图像的相似度Y=100%;
步骤S3:对不同初始化区域的原始管道图像数据,依次进行均值、样本方差、标准化处理;
步骤S4:利用图像极值点的检测方法对管道图像进行特征提取,并用尺度不变特征转换描述子对管道图像的局部特征进行表示,分为四步:检测管道图像关键点、精确定位管道泄漏关键点、关键点分配主方向、计算尺度不变特征转换描述子;
步骤S5:利用特征融合算法对特征提取的结果,把相同管道图像聚合而不同管道图像进行分离;
步骤S6:选择管道图像的相邻区域,并将两块相邻区域进行匹配(Fi,Gj),计算两块相邻区域特征融合结果的相似度y(Fi,Gj),Y=Y∪y(Fi,Gj);
步骤S7:若Y≠0,令y(Fi,Gj)=max(Y),合并两块相邻区域Pz=Fi∪Gj,去除关于Fi和Gj的相似度Y=Y\y(Fi,G*)和Y=Y\y(F*,Gj);
步骤S8:迭代步骤S6和步骤S7,直到Y=0;
步骤S9:从所有区域中定位出泄漏点的位置。
进一步地,所述图像分割算法,利用高斯马尔可夫随机场模型对原始管道图像进行建模,模型的表达式为:
Figure DEST_PATH_IMAGE002
其中,B表示原始管道图像,d表示原始管道图像中一截管道的位置,γ表示相邻位置,泄漏点的位置B(d)=Bd的概率表示A(B(d)=Bd),简化为A(Bd),C表示自定义的平面;
原始管道图像B用高斯马尔可夫随机场模型模拟,Bd的表达式为:
Figure DEST_PATH_IMAGE003
其中,ηα表示高斯马尔可夫随机场的模型参数,d表示管道任一点的位置图像,a表示图像的尺寸大小,Bd+a表示管道图像的离散区域,bd表示零均值的高斯噪声,利用析取正态水平集对原始管道图像的纹理中像素点之间的相互作用进行表示,析取正态水平集的能量函数公式为:
Figure DEST_PATH_IMAGE004
其中,e1和e2分别表示原始管道图像前景区域、背景区域的灰度均值,β表示图像的全部区域集,B(n)表示能量函数的均值,dn表示对管道图像的积分,n表示单个的图像区域,f(n)表示观测图像;
所述前景区域的灰度均值表达式为:
Figure DEST_PATH_IMAGE005
其中,G表示为海维赛德函数,f(n)表示观测图像;
所述背景区域的灰度均值表达式为:
Figure DEST_PATH_IMAGE006
其中,f(n)表示观测图像,计算迭代e1,e2,直到收敛为止,输出分割后的图像,完成图像分割。
进一步地,所述均值,表达式为:
Figure DEST_PATH_IMAGE007
其中,
Figure DEST_PATH_IMAGE008
表示管道图像均值结果,x表示任一的管道图像数据点,ε表示管道图像最小批处理值,
Figure DEST_PATH_IMAGE009
表示管道图像数据点个数,h表示所有管道图像数据点灰度值的和;
所述样本方差,表达式为:
Figure DEST_PATH_IMAGE010
其中,
Figure DEST_PATH_IMAGE011
表示管道图像的样本方差值;
所述标准化处理,表达式为:
Figure DEST_PATH_IMAGE012
其中,
Figure DEST_PATH_IMAGE013
表示管道图像的标准化值
进一步地,所述特征提取,包括检测管道图像关键点、精确定位管道泄漏关键点、关键点分配主方向、计算尺度不变特征转换描述子;
检测管道图像关键点,表达式为:
Figure DEST_PATH_IMAGE014
其中,符号*表示卷积,H(j,k,φ)表示管道图像关键点结果,I(j,k,φ)表示尺度φ可变高斯函数,L(j,k)表示尺度空间,j,k表示管道图像的像素点横纵坐标;
Figure DEST_PATH_IMAGE015
的表达式为:
Figure DEST_PATH_IMAGE016
其中,φ表示尺度空间因子,利用H(j,k,φ)检测到管道图像的稳定关键点,定义DoG尺度空间为:
Figure DEST_PATH_IMAGE017
其中,M(j,k,φ)表示不同尺度的空间,p表示尺度空间系数。
进一步地,所述精确定位管道泄漏关键点,包括去除管道图像的边缘响应点和去除低对比度点,获取关键点处的拟和函数。
进一步地,所述关键点分配主方向,以关键点为中心的邻域窗口内进行采样,并计算采样点的梯度方向角,表达式为:
Figure DEST_PATH_IMAGE018
其中,j,k表示管道图像的像素点横纵坐标,
Figure DEST_PATH_IMAGE019
表示方向角,方向角每隔10°取得一个方向,在0~360°的范围内一共有36个方向,再统计所有采样点的梯度方向角直方图,其峰值则代表了该关键点的主方向。
进一步地,所述计算尺度不变特征转换描述子,步骤为:
①以每个关键点为中心,将管道图像坐标X轴旋转到尺度不变特征转换关键点的主方向上;
②以每个关键点为中心分割成子区域,在每个子区域上计算高斯加权梯度方向角直方图,绘制每个梯度方向的累加值,即可形成一个种子点。
进一步地,所述步骤S5,边缘特征融合算法:
Figure DEST_PATH_IMAGE020
其中,zw表示特征融合系数,z表示任一管道图像,w表示管道图像的尺寸大小,
Figure DEST_PATH_IMAGE021
Figure DEST_PATH_IMAGE022
分别表示管道数据集类间散度和整体散度,公式依次如下:
Figure DEST_PATH_IMAGE023
Figure DEST_PATH_IMAGE024
其中,
Figure DEST_PATH_IMAGE025
Figure DEST_PATH_IMAGE026
分别表示第x类样本与全部样本的均值,hx表示单个样本值,
Figure DEST_PATH_IMAGE027
表示正则项,R(z)表示先验知识,
Figure DEST_PATH_IMAGE028
Figure DEST_PATH_IMAGE029
表示高于0的常数,x表示任一的样本数据,r表示散度次数,n表示样本总数。
进一步地,所述相似度,建立相似度函数的表达式为:
Figure DEST_PATH_IMAGE030
其中,F,G表示两张相邻的管道图像,E表示所有经过边缘特征融合的管道图像合集,rx表示图像x的相似度系数,Sx表示图像所有的特征点合集。
本发明通过提出一种基于图像识别快速定位泄漏点的方法通过对管道图像处理判断出泄漏点,该方法是在现有技术的基础上,利用图像极值点的检测方法对管道图像进行特征提取,可以迅速的识别图像特征,并进行提取,同时利用特征融合算法完成图像的分离,最后判断相邻区域图像的相似性迅速找到管道的泄漏点,该方法操作简单可迅速对整条管道进行检测,判断出泄漏点,在管道安装时在电脑系统匹配该方法中可以实现全部泄漏过程的自动化过程,节约了人力成本。
附图说明
图1为本发明总体步骤流程图;
图2为本发明的精确定位管道泄漏关键点流程图;
图3为本发明的图像相似度计算流程图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互结合,下面结合附图和有具体实施例对本申请作进一步详细说明。
如图1所示,一种基于图像识别快速定位泄漏点的方法,该方法步骤为:
步骤S1:利用高清摄像头对管道进行拍摄,摄像头每隔5米对管道进行一次采样拍摄,对同距离管道的向阳与背阳部分同时进行采样,在电脑系统或者云端输入原始管道的图像数据,将同距离管道的向阳与背阳的图像设置为一组图像;
步骤S2:将同距离管道的向阳与背阳的一组图像组利用图像分割算法对原始管道图像进行分割,分割的目的是为了后续对管道图像进行更加细微的特征提取,将分割后的图像生成不同初始化区域
Figure 658627DEST_PATH_IMAGE001
,每个区域图像的大小相同,格式相同,初始化后默认所有图像的相似度Y=100%;
步骤S3:对不同初始化区域的原始管道图像数据,依次进行均值、样本方差、标准化处理,该步骤的目的是判断每张管道图像的数据格式与数据量是否达到特征提取的标准,为了图像的特征提取提供统一的数据格式;
步骤S4:利用图像极值点的检测方法对管道图像进行特征提取,并用尺度不变特征转换描述子对管道图像的局部特征进行表示,分为四步:检测管道图像关键点、精确定位管道泄漏关键点、关键点分配主方向、计算尺度不变特征转换描述子;其中检测管道图像关键点是依据历史经验对每张图像中容易出现泄漏点的部位首先进行检测;精确定位管道泄漏关键点是对关键点检查完成后其余图像点进行全部检测,精确定位出该图像中所有的泄漏点;关键点分配主方向,是以关键点为中心对周围所有的点进行再次检测,判断关键点一定半径范围内是否有泄漏出现,采用的是梯度方向角直方图的形式;计算尺度不变特征转换描述子,以每个关键点为中心分割成子区域,在每个子区域上计算高斯加权梯度方向角直方图,绘制每个梯度方向的累加值,即可形成一个种子点,利用种子点再次检测图像中的泄漏点;
步骤S5:根据特征融合算法对特征提取的结果,把相同管道图像聚合而不同管道图像进行分离,泄漏点的特征图像因为泄漏原因以及气候的变化会出现不同的表现方式,特征融合的目的是为了把每次出现泄漏点的特征图像进行更新,以便于在下次管道检测时,可以迅速的判断出泄漏点;
步骤S6:在管道周围安装带有红外成像的摄像头,或者采用带有红外成像无人机对管道进行视频数据采集,因为泄露的气体有可能是无色的,如果只通过传统的摄像头无法监测到,必须得使用带有红外成像的摄像设备,所获得的数据是红外热成像图像;选择管道红外热成像图像的相邻区域,并将两块相邻区域进行匹配(Fi,Gj),计算两块相邻区域特征融合结果的相似度y(Fi,Gj),Y=Y∪y(Fi,Gj),其中,Fi和Gi表示管道相邻的两块图像区域,y表示相邻管道图像的相似度,Y表示管道所有图像的相似度;
步骤S7:若Y≠0,令y(Fi,Gj)=max(Y),合并两块相邻区域Pz=Fi∪Gj,去除关于Fi和Gj的相似度Y=Y\y(Fi,G*)和Y=Y\y(F*,Gj);其中,Pz表示合并的管道图像,G*和F*表示无泄漏的两块相邻管道图像;
步骤S8:迭代步骤S6和步骤S7,直到Y=0,因为出现泄漏在图像上的表现形式为出现漏点或漏缝,当有漏点或漏缝出现时的图像与完好管道图像的相似度是完全不一样的,所以当Y=0时即可判断出泄漏的出现;
步骤S9:从所有区域中定位出泄漏点的位置,对所有的管道图像进行步骤S2至步骤S8,即可对管道进行全部的检测,找出所有的泄漏点。
在步骤S2中,图像分割算法,利用高斯马尔可夫随机场模型对原始管道图像进行建模,模型的表达式为:
Figure 106926DEST_PATH_IMAGE002
其中,B表示原始管道图像,d表示原始管道图像中一截管道的位置,γ表示相邻位置,泄漏点的位置B(d)=Bd的概率表示A(B(d)=Bd),简化为A(Bd),C表示自定义的平面;
原始管道图像B用高斯马尔可夫随机场模型模拟,Bd的表达式为:
Figure 417821DEST_PATH_IMAGE003
其中,ηα表示高斯马尔可夫随机场的模型参数,d表示管道任一点的位置图像,a表示图像的尺寸大小,Bd+a表示管道图像的离散区域,bd表示零均值的高斯噪声,利用析取正态水平集对原始管道图像的纹理中像素点之间的相互作用进行表示,析取正态水平集的能量函数公式为:
Figure 344189DEST_PATH_IMAGE004
其中,e1和e2分别表示原始管道图像前景区域、背景区域的灰度均值,β表示图像的全部区域集,B(n)表示能量函数的均值,dn表示对管道图像的积分,n表示单个的图像区域,f(n)表示观测图像;
前景区域的灰度均值表达式为:
Figure 689720DEST_PATH_IMAGE005
其中,G表示为海维赛德函数,f(n)表示观测图像;
背景区域的灰度均值表达式为:
Figure 43341DEST_PATH_IMAGE006
其中,f(n)表示观测图像,计算迭代e1,e2,直到收敛为止,输出分割后的图像,完成图像分割。
在步骤S3中,均值表达式为:
Figure 841532DEST_PATH_IMAGE007
其中,
Figure 571591DEST_PATH_IMAGE008
表示管道图像均值结果,x表示任一的管道图像数据点,ε表示管道图像最小批处理值,
Figure 506049DEST_PATH_IMAGE009
表示管道图像数据点个数,h表示所有管道图像数据点灰度值的和;
样本方差表达式为:
Figure 30571DEST_PATH_IMAGE010
其中,
Figure 316059DEST_PATH_IMAGE011
表示管道图像的样本方差值;
标准化处理表达式为:
Figure 584229DEST_PATH_IMAGE012
其中,
Figure 638773DEST_PATH_IMAGE013
表示管道图像的标准化值。
在步骤S4中特征提取算法包括,检测管道图像关键点、精确定位管道泄漏关键点、关键点分配主方向、计算尺度不变特征转换描述子;
检测管道图像关键点的表达式为:
Figure 334196DEST_PATH_IMAGE014
其中,符号*表示卷积,H(j,k,φ)表示管道图像关键点结果,I(j,k,φ)表示尺度φ可变高斯函数,L(j,k)表示尺度空间,j,k表示管道图像的像素点横纵坐标;
Figure 106980DEST_PATH_IMAGE015
的表达式为:
Figure 913262DEST_PATH_IMAGE016
其中,φ表示尺度空间因子,利用H(j,k,φ)检测到管道图像的稳定关键点,定义DoG尺度空间为:
Figure 87892DEST_PATH_IMAGE017
其中,M(j,k,φ)表示不同尺度的空间,p表示尺度空间系数。
如图2所示,精确定位管道泄漏关键点包括:去除管道图像的边缘响应点和去除低对比度点,获取关键点处的拟和函数。
去除边缘响应点,构造2x2的Hessian矩阵:
Figure DEST_PATH_IMAGE031
设ψ,μ分别为N最大、最小特征值,则:
Figure DEST_PATH_IMAGE032
其中,Tr(N)表示特征值N的离散矩阵,Det(N)表示特征值N的行列式;
令ψ=ξμ,则:
Figure DEST_PATH_IMAGE033
因为(ξ+1)2/ξ在两特征值相等时达最小,随ξ的增长而增长,取ξ=10时,当:
Figure DEST_PATH_IMAGE034
其中,Det(H)表示特征值H的行列式,成立时,将该关键点保留,反之则视作边缘响应点被剔除。
去除低对比度点,获取关键点处的拟和函数:
Figure DEST_PATH_IMAGE035
其中,T表示管道的关键点,Q表示管道不同的尺度空间矩阵,w表示空间矩阵的转置。
求导并让方程等于零,可以得到极值点:
Figure DEST_PATH_IMAGE036
对应极值点,方程的值为:
Figure DEST_PATH_IMAGE037
通常,将‖
Figure DEST_PATH_IMAGE038
‖<0.05的极值点视为低对比度的不稳定特征点,进行剔除。
关键点分配主方向,以关键点为中心的邻域窗口内进行采样,并计算采样点的梯度方向角,表达式为:
Figure 672326DEST_PATH_IMAGE018
其中,j,k表示管道图像的像素点横纵坐标,
Figure 666827DEST_PATH_IMAGE019
表示方向角,方向角每隔10°取得一个方向,在0~360°的范围内一共有36个方向,再统计所有采样点的梯度方向角直方图,其峰值则代表了该关键点的主方向。
计算尺度不变特征转换描述子步骤为:
①以每个关键点为中心,将图像坐标X轴旋转到尺度不变特征转换关键点的主方向上,以确保旋转不变性;②以每个关键点为中心取大小为16×16个像素的窗口,且分割成4×4个子区域,在每个子区域上计算8个方向的高斯加权梯度方向角直方图,越靠近关键点的像素梯度方向信息贡献越大,绘制每个梯度方向的累加值,即可形成一个种子点,因此,一共可以生成16个种子点,这样对于每个关键点就可以产生一个长度为128的数据,即最终得到一个128维的尺度不变特征转换描述子。
在步骤S5中,边缘特征融合算法:
Figure 276799DEST_PATH_IMAGE020
其中,zw表示特征融合系数,z表示任一管道图像,w表示管道图像的尺寸大小,
Figure 305935DEST_PATH_IMAGE021
Figure 608741DEST_PATH_IMAGE022
分别表示管道数据集类间散度和整体散度,公式依次如下:
Figure 90538DEST_PATH_IMAGE023
Figure 504201DEST_PATH_IMAGE024
其中,
Figure 122265DEST_PATH_IMAGE025
Figure 595971DEST_PATH_IMAGE026
分别表示第x类样本与全部样本的均值,hx表示单个样本值,
Figure 565064DEST_PATH_IMAGE027
表示正则项,R(z)表示先验知识,
Figure 516840DEST_PATH_IMAGE028
Figure 520568DEST_PATH_IMAGE029
表示高于0的常数,x表示任一的样本数据,r表示散度次数,n表示样本总数。
如图3所示,在步骤S6中的相似度计算,第一步:离线构造金字塔特征空间
①在管道图像库,随机选择N幅并分别提取它们的尺度不变特征转换特征子,记作:
Figure DEST_PATH_IMAGE039
其中,Ax表示单个管道图像的的特征子,B表示所有图像尺度不变特征转换描述子的总数。
②设置金字塔的层数E与分支数K。首先,将U的数据用K-Mean方法聚成K类,称所有的聚类中心为金字塔的第1层;然后,分别对这K类数据再用K-Mean方法又聚成K类,称这K*K个聚类中心为金字塔的第2层;依次类推,直到第E层,则停止聚类。最后,从第1层到第E层所有的聚类中心,称为金字塔特征空间,记为:
Figure DEST_PATH_IMAGE040
其中Dx表示第i层的所有聚类中心。
第二步:相似度计算
1.多分辨率直方图特征提取,设X表示某图像的尺度不变特征转换描述子集合,计算多分辨率直方图特征,记为:
Figure DEST_PATH_IMAGE041
其中Mi为第i层的直方图特征。
2.金字塔匹配函数相似度计算,计算F与G的多分辨率直方图,则F与G之间的相似度定义为:
Figure 493072DEST_PATH_IMAGE030
其中,F,G表示两张相邻的管道图像,E表示所有经过边缘特征融合的管道图像合集,rx表示图像x的相似度系数,Sx表示图像所有的特征点合集。
本发明通过提出一种基于图像识别快速定位泄漏点的方法通过对管道图像处理判断出泄漏点,该方法是在现有技术的基础上,利用图像极值点的检测方法对管道图像进行特征提取,可以迅速的识别图像特征,并进行提取,同时利用特征融合算法完成图像的分离,最后判断相邻区域图像的相似性迅速找到管道的泄漏点,该方法操作简单可迅速对整条管道进行检测,判断出泄漏点,在管道安装时在电脑系统匹配该方法中可以实现全部泄漏过程的自动化过程,节约了人力成本,可进行大规模推广进行。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解的是,在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种等效的变化、修改、替换和变型,本发明的范围由所附权利要求及其等同范围限定。

Claims (9)

1.一种基于图像识别快速定位泄漏点的方法,其特征在于,该方法步骤为:
步骤S1:输入原始管道的图像数据;
步骤S2:利用图像分割算法对原始管道图像进行分割,生成图像的不同初始化区域
Figure 117355DEST_PATH_IMAGE001
,初始化后默认所有图像的相似度Y=100%;
步骤S3:对不同初始化区域的原始管道图像数据,依次进行均值、样本方差、标准化处理;
步骤S4:利用图像极值点的检测方法对管道图像进行特征提取,并用尺度不变特征转换描述子对管道图像的局部特征进行表示,分为四步:检测管道图像关键点、精确定位管道泄漏关键点、关键点分配主方向、计算尺度不变特征转换描述子;
步骤S5:根据特征融合算法对特征提取的结果,把相同管道图像聚合而不同管道图像进行分离;
步骤S6:选择管道图像的相邻区域,并将两块相邻区域进行匹配(Fi,Gj),计算两块相邻区域特征融合结果的相似度y(Fi,Gj),Y=Y∪y(Fi,Gj);
步骤S7:若Y≠0,令y(Fi,Gj)=max(Y),合并两块相邻区域Pz=Fi∪Gj,去除关于Fi和Gj的相似度Y=Y\y(Fi,G*)和Y=Y\y(F*,Gj);
步骤S8:迭代步骤S6和步骤S7,直到Y=0;
步骤S9:从所有区域中定位出泄漏点的位置。
2.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述图像分割算法,利用高斯马尔可夫随机场模型对原始管道图像进行建模,模型的表达式为:
Figure 747051DEST_PATH_IMAGE002
其中,B表示原始管道图像,d表示原始管道图像中一截管道的位置,γ表示相邻位置,泄漏点的位置B(d)=Bd的概率表示A(B(d)=Bd),简化为A(Bd),C表示自定义的平面;
原始管道图像B用高斯马尔可夫随机场模型模拟,Bd的表达式为:
Figure 280801DEST_PATH_IMAGE003
其中,ηα表示高斯马尔可夫随机场的模型参数,d表示管道任一点的位置图像,a表示图像的尺寸大小,Bd+a表示管道图像的离散区域,bd表示零均值的高斯噪声,利用析取正态水平集对原始管道图像的纹理中像素点之间的相互作用进行表示,析取正态水平集的能量函数公式为:
Figure 430284DEST_PATH_IMAGE004
其中,e1和e2分别表示原始管道图像前景区域、背景区域的灰度均值,β表示图像的全部区域集,B(n)表示能量函数的均值,dn表示对管道图像的积分,n表示单个的图像区域,f(n)表示观测图像;
所述前景区域的灰度均值表达式为:
Figure 47079DEST_PATH_IMAGE005
其中,G表示为海维赛德函数,f(n)表示观测图像;
所述背景区域的灰度均值表达式为:
Figure 429650DEST_PATH_IMAGE006
其中,f(n)表示观测图像,计算迭代e1,e2,直到收敛为止,输出分割后的图像,完成图像分割。
3.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述均值,表达式为:
Figure 658769DEST_PATH_IMAGE007
其中,
Figure 239923DEST_PATH_IMAGE008
表示管道图像均值结果,x表示任一的管道图像数据点,ε表示管道图像最小批处理值,
Figure 106248DEST_PATH_IMAGE009
表示管道图像数据点个数,h表示所有管道图像数据点灰度值的和;
所述样本方差,表达式为:
Figure 22120DEST_PATH_IMAGE010
其中,
Figure 507459DEST_PATH_IMAGE011
表示管道图像的样本方差值;
所述标准化处理,表达式为:
Figure 959431DEST_PATH_IMAGE012
其中,
Figure 668761DEST_PATH_IMAGE013
表示管道图像的标准化值。
4.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述特征提取,包括检测管道图像关键点、精确定位管道泄漏关键点、关键点分配主方向、计算尺度不变特征转换描述子;
检测管道图像关键点,表达式为:
Figure 540771DEST_PATH_IMAGE014
其中,符号*表示卷积,H(j,k,φ)表示管道图像关键点结果,I(j,k,φ)表示尺度φ可变高斯函数,L(j,k)表示尺度空间,j,k表示管道图像的像素点横纵坐标;
Figure 626539DEST_PATH_IMAGE015
的表达式为:
Figure 510181DEST_PATH_IMAGE016
其中,φ表示尺度空间因子,利用H(j,k,φ)检测到管道图像的稳定关键点,尺度空间表达式为:
Figure 78828DEST_PATH_IMAGE017
其中,M(j,k,φ)表示不同尺度的空间,p表示尺度空间系数。
5.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述精确定位管道泄漏关键点,包括去除管道图像的边缘响应点和去除低对比度点,获取关键点处的拟和函数。
6.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述关键点分配主方向,以关键点为中心的邻域窗口内进行采样,并计算采样点的梯度方向角,表达式为:
Figure 720025DEST_PATH_IMAGE018
其中,j,k表示管道图像的像素点横纵坐标,
Figure 937380DEST_PATH_IMAGE019
表示方向角,方向角每隔10°取得一个方向,在0~360°的范围内一共有36个方向,再统计所有采样点的梯度方向角直方图,其峰值则代表了该关键点的主方向。
7.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述计算尺度不变特征转换描述子,步骤为:
①以每个关键点为中心,将管道图像坐标X轴旋转到尺度不变特征转换关键点的主方向上;
②以每个关键点为中心分割成子区域,在每个子区域上计算高斯加权梯度方向角直方图,绘制每个梯度方向的累加值,即可形成一个种子点。
8.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述特征融合算法,表达式为:
Figure 534583DEST_PATH_IMAGE020
其中,zw表示特征融合系数,z表示任一管道图像,w表示管道图像的尺寸大小,
Figure 179191DEST_PATH_IMAGE021
Figure 307684DEST_PATH_IMAGE022
分别表示管道数据集类间散度和整体散度,公式依次如下:
Figure 689249DEST_PATH_IMAGE023
Figure 281904DEST_PATH_IMAGE024
其中,
Figure 769518DEST_PATH_IMAGE025
Figure 572257DEST_PATH_IMAGE026
分别表示第x类样本与全部样本的均值,hx表示单个样本值,
Figure 865835DEST_PATH_IMAGE027
表示正则项,R(z)表示先验知识,
Figure 250680DEST_PATH_IMAGE028
Figure 863190DEST_PATH_IMAGE029
表示高于0的常数,x表示任一的样本数据,r表示散度次数,n表示样本总数。
9.如权利要求1所述的一种基于图像识别快速定位泄漏点的方法,其特征在于,所述相似度,建立相似度函数的表达式为:
Figure 28592DEST_PATH_IMAGE030
其中,F,G表示两张相邻的管道图像,E表示所有经过边缘特征融合的管道图像合集,rx表示图像x的相似度系数,Sx表示图像所有的特征点合集。
CN202210477646.3A 2022-05-05 2022-05-05 一种基于图像识别快速定位泄漏点的方法 Active CN114581440B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210477646.3A CN114581440B (zh) 2022-05-05 2022-05-05 一种基于图像识别快速定位泄漏点的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210477646.3A CN114581440B (zh) 2022-05-05 2022-05-05 一种基于图像识别快速定位泄漏点的方法

Publications (2)

Publication Number Publication Date
CN114581440A true CN114581440A (zh) 2022-06-03
CN114581440B CN114581440B (zh) 2024-01-09

Family

ID=81779279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210477646.3A Active CN114581440B (zh) 2022-05-05 2022-05-05 一种基于图像识别快速定位泄漏点的方法

Country Status (1)

Country Link
CN (1) CN114581440B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116228770A (zh) * 2023-05-09 2023-06-06 深蓝(天津)智能制造有限责任公司 一种管道泄漏的识别监测方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001337002A (ja) * 2000-05-29 2001-12-07 Tokyo Gas Co Ltd 漏洩地点推定支援システムおよび漏洩地点推定方法
CN106764465A (zh) * 2015-11-25 2017-05-31 李福霞 一种基于多传感器的管道泄漏检测装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107435817B (zh) * 2017-08-15 2019-01-25 常州大学 一种压力管道两点泄漏检测精确定位方法
CN113743637A (zh) * 2020-05-29 2021-12-03 青岛海尔电冰箱有限公司 基于协同过滤的煤气浓度预测方法、设备及冰箱

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001337002A (ja) * 2000-05-29 2001-12-07 Tokyo Gas Co Ltd 漏洩地点推定支援システムおよび漏洩地点推定方法
CN106764465A (zh) * 2015-11-25 2017-05-31 李福霞 一种基于多传感器的管道泄漏检测装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Y. ZHOU: "Study on Oil and Gas Pipeline Leakage Real-time Inspection System Based on Distributed Optical Fiber", 《KEY ENGINEERING MATERIALS》, pages 447 - 450 *
綦磊: "基于声发射信号EMD-WPD特征融合的航天器在轨泄漏辨识方法", 《震动与冲击》, pages 110 - 116 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116228770A (zh) * 2023-05-09 2023-06-06 深蓝(天津)智能制造有限责任公司 一种管道泄漏的识别监测方法及系统
CN116228770B (zh) * 2023-05-09 2023-09-01 深蓝(天津)智能制造有限责任公司 一种管道泄漏的识别监测方法及系统

Also Published As

Publication number Publication date
CN114581440B (zh) 2024-01-09

Similar Documents

Publication Publication Date Title
CN109255317B (zh) 一种基于双网络的航拍图像差异检测方法
Chen et al. Vehicle detection in high-resolution aerial images via sparse representation and superpixels
CN107133569B (zh) 基于泛化多标记学习的监控视频多粒度标注方法
CN111784746A (zh) 一种鱼眼镜头下行人多目标跟踪方法、装置及计算机系统
CN111310662B (zh) 一种基于集成深度网络的火焰检测识别方法及系统
CN112200121B (zh) 基于evm和深度学习的高光谱未知目标检测方法
CN112446370B (zh) 一种电力设备铭牌文本信息识别的方法
CN116092013A (zh) 一种用于智慧监控的危险路况识别方法
CN112258490A (zh) 基于光学和红外图像融合的低发射率涂层智能探损方法
CN113313107B (zh) 一种斜拉桥缆索表面多类型病害智能检测和识别方法
CN115861210B (zh) 一种基于孪生网络的变电站设备异常检测方法和系统
CN111860393A (zh) 一种安防系统上的人脸检测识别方法
Hendryli et al. Automatic license plate recognition for parking system using convolutional neural networks
CN114581440A (zh) 一种基于图像识别快速定位泄漏点的方法
CN115639248A (zh) 一种用于建筑外墙质量检测的系统及方法
Tang et al. Similarity and clustering of footwear prints
CN111105436B (zh) 目标跟踪方法、计算机设备及存储介质
CN112085117B (zh) 一种基于MTLBP-Li-KAZE-R-RANSAC的机器人运动监测视觉信息融合方法
CN112418262A (zh) 车辆再识别的方法、客户端及系统
CN114998792B (zh) 一种具有ai网络摄像头的安全监控方法
CN111539931A (zh) 一种基于卷积神经网络和边界限定优化的外观异常检测方法
CN116110095A (zh) 一种人脸过滤模型的训练方法、人脸识别方法及装置
CN114758135A (zh) 一种基于注意力机制的无监督图像语义分割方法
CN113688819A (zh) 一种基于标记点的目标物期望点跟踪匹配方法
CN115761211A (zh) 一种基于rgb图像与热成像耦合的石化企业泵房设备温度监测方法及装置

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