CN103745238B - 基于AdaBoost和主动形状模型的受电弓识别方法 - Google Patents
基于AdaBoost和主动形状模型的受电弓识别方法 Download PDFInfo
- Publication number
- CN103745238B CN103745238B CN201310566693.6A CN201310566693A CN103745238B CN 103745238 B CN103745238 B CN 103745238B CN 201310566693 A CN201310566693 A CN 201310566693A CN 103745238 B CN103745238 B CN 103745238B
- Authority
- CN
- China
- Prior art keywords
- pantograph
- image
- adaboost
- shape model
- active shape
- 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
Abstract
本发明属于计算机数字图像处理与模式识别技术领域,主要涉及电动牵引机车使用的受电弓在线自动识别方法,具体涉及一种基于AdaBoost和主动形状模型的受电弓在线自动识别方法。基本流程包括:通过受电弓在线拍照系统采集若干受电弓图像,组成受电弓图像学习样本集;基于样本训练学习生成受电弓主动形状模型和AdaBoost级联分类器;采用AdaBoost级联分类器对新获取、待识别的图像进行受电弓检测;结合受电弓检测结果和学习构建的受电弓主动形状模型,精确匹配和识别受电弓形状;最后在精确匹配结果基础上进行定量检测和分析。本发明能有效对受电弓碳滑板厚度进行线上自动定量检测,可快速地进行车辆维护,节约受电弓检测成本。
Description
技术领域
本发明属于计算机数字图像处理与模式识别技术领域,主要涉及电动牵引机车使用的受电弓在线自动识别方法,具体涉及一种基于AdaBoost和主动形状模型的受电弓识别方法。
背景技术
电力牵引机车受电弓是电气化铁路电力机车从接触网上受取电源的装置。受电弓弓头破损、碳滑板过度磨耗不仅影响电力机车的正常供电,甚至会造成轨道交通的中断。因此作为直接与接触网接触的取流关键设备-受电弓碳滑板的状态极为关键,它关系到机车能否稳定、可靠地从接触线取电。实际管理中,相关部门需要对受电弓碳滑板运行状态进行定期检测,以保证机车的稳定、安全运行。随着我国轨道交通和高速铁路的飞速发展,对受电弓的可靠运行提出了更高的要求,实现对受电弓运行状态的在线自动识别具有重要意义。
目前国内外受电弓状态的检测方法主要包括入库静态检测、机载检测系统和在线定点式检测3种方式。
入库静态检测要求接触网必须在受电弓检测前断电,且通常只能在机车入库停车后才能检测,不能实现受电弓的在线自动检测。在实际运作中需要投入较大规模的人力、成本较高。详见文献:谢力.基于图像处理的受电弓状态检测技术研究:[硕士学位论文]:西南交通大学,2009。
机载受电弓检测方法属于一种在线自动检测方法。主要包括2种方式,其一是通过在机车上安装光、电、力学等多种传感器采集受电弓位置、碳滑板的磨耗,以及受电弓弓头左右偏移量状态等信息对受电弓进行检测,但传感器获取的受电弓运行状态信息容易受机车运行过程中产生的震动的干扰;其二是将光纤内埋式磨损传感器嵌入受电弓滑板中,当受电弓滑板收到接触网的冲击作用而产生磨损时或者产生缺失时,传感器给出相应的磨损信号。该方法的主要缺点是对材料的要求比较高,要在受电弓滑板中内嵌传感器,材料工艺要求高,受电弓滑板材料结构改变可能会降低受电弓的使用寿命,提高了成本。
在线定点式自动检测方式,国内外主要有基于超声波传感器的检测和基于图像的检测方法。超声波传感器检测方法多见于国外机车系统,该方法存在的显著不足是超声波信号容易受超声波工作距离、超声波信号强弱和受电弓运行中的姿态的影响。
近年来,国内外出现了许多基于在线获取图像的受电弓检测方法。该方法核心包括拍照和识别2部分,首先完成受电弓图像的在线采集,然后采用数字图像处理方法对采集到的图像进行仔细分析,最终获取受电弓在线运行状态。该方法的优势非常明显,可以在机车行驶过程中做到在线自动识别;缺点是精度依赖于拍照系统采集图像的质量,对受电弓在线拍照系统精度要求比较高。国内目前已有厂商和科研单位正在研发基于图像的受电弓在线运行状态的识别系统,存在的问题主要集中于:1)受电弓图像采集质量较低,不利于后续图像处理与分析;2)对机车运行速度有限制,难以实现真正的在线识别;3)采集到的受电弓图像背景复杂、不统一,图像处理方法难以进行精确的受电弓定量识别,相关方法亟需进一步改进。详见南京大学于2012年公开的专利《高速机车受电弓滑板磨损自动检测装置》,专利公开号:CN102507600A。
基于可变模型对目标物体特征点定位的研究取得了很大成就,特别是对于图像中的目标物体变化很大的情况。可变模型的基本思想是:建立一个目标物体的通用模型,对于任意给定图像,如果图像中有与该模型相同的目标物体,则可以根据物体图像的内部和外部特征与模型之间的差异,通过调整模型参数将模型进行变化,其形状和纹理能够与目标物体在一定误差内实现匹配,详见文献:Cootes T.F.Deformable Object Modelling andMatching.In:Kimmel R,Klette R,Sugimoto A,eds.Computer Vision-Accv2010,PtI.Berlin:Springer-Verlag Berlin;2011:1-10.
英国曼彻斯特大学(University of Manchester)的Tim.Cootes等人提出的主动形状模型ASM(Active Shape Models)是实现上述思想的一个典型代表。该方法已成为当前主流的目标定位方法,在目标对象的定位中,特别是在人脸检测、手势识别和医学人体器官图像检测方面得到了广泛的应用,详见文献:Cootes T.F.,C.J.Taylor,D.H.Cooper,etal.Active Shape Models-Their Training and Application.Computer Vision andImage Understanding.1995,61(1):38-59.
ASM是一种基于统计学的灰度和形状分离的可变形模型,允许待定位目标有一定程度的形变,同时采用灰度和梯度信息指导形状模型收敛,其收敛速度较快。其主要方法是通过对样本图像中的目标物体进行关键特征点手动标定,进而形成样本集目标点分布模型(Point Distribution Model),然后对点分布模型进行训练、统计分析建立ASM模型,再以此模型作为依据,在测试图像中与目标特征进行快速匹配,找到被定位目标的位置。其优点是容易选择目标轮廓特征点作为建模基础,但ASM的匹配精度与匹配起始位置密切相关,因此获取ASM匹配起始位置的目标定位方法是影响ASM最终匹配结果优劣的一个重要因素。实际应用ASM方法进行目标精确定位时,需要首先进行目标的初略定位,否则匹配过程容易陷入局部最优,导致匹配不到实际的目标上。因此,本发明在受电弓进行ASM精确识别之前,采用基于AdaBoost的级联分类器对受电弓进行检测,保证ASM识别精度。
AdaBoost算法是一种分类器算法,是由Freund和Schapire在1995年提出的,详见文献:Freund,Y.and Schapire,R.E.1995.A decision-theoretic generalization ofon-line learning and an application to boosting.In Computational LearningTheory:Eurocolt95,Springer-Verlag,pp.23-37。其基本思想是利用大量的分类能力一般的弱分类器(weaker classifier)通过一定方法组合(boost)起来,构成一个分类能力很强的强分类器。理论证明,只要每个简单分类器分类能力比随机猜测要好,当简单分类器个数趋向于无穷时,强分类器的错误率将趋于零。详见文献:龙伶敏,基于Adaboost的人脸检测方法及眼睛定位算法研究[D],电子科技大学,2008。
基于Harr-like特征的AdaBoost级联分类器最初应用于人脸检测,详见文献:PaulViola and Michael J.Jones.Robust Real-Time Face Detection.InternationalJournal of Computer Vision,Vol.57,pp.137-154,May2004。该方法提供了一种通用的目标检测框架,.使用者只需建立某目标的训练样本集,就可以方便地完成特定目标物检测模型的建立和检测。详见北京航空航天大学于2011年公开的专利《基于冗余特征消减AdaBoost分类器的物体检测方法》,专利公开号:CN102129565A。
随着我国高速铁路建设的不断发展,机车的行驶速度越来越快,对受电弓运行质量的要求越来越高,同时也对受电弓在线自动识别的精确性、稳定性提出了更高的要求。针对目前国内基于图像的受电弓在线识别方法精度差的问题,迫切需要研发一种新型的受电弓运行状态在线自动识别方法。
发明内容
(一)发明目的
本发明的目的是:针对基于图像的受电弓识别,提供一种基于AdaBoost和主动形状模型的受电弓在线自动识别方法,在机车正常行驶过程中,能够对受电弓碳滑板厚度进行定量检测。
(二)技术解决方案
本发明提供了一种基于AdaBoost和主动形状模型的受电弓自动识别方法,包括以下步骤:
步骤10、获取受电弓图像学习样本集;
步骤20、AdaBoost级联分类器训练;
步骤30、采用AdaBoost级联分类器对新获取、待识别的图像进行受电弓检测;
步骤40、构建受电弓主动形状模型;
步骤50、结合检测结果和主动形状模型的受电弓精确匹配和识别;
步骤60、受电弓定量检测和分析。
所述步骤10中获取受电弓图像学习样本集包括两套:AdaBoost学习样本集:用于训练AdaBoost级联分类器;主动形状模型学习样本集:用于训练受电弓主动形状模型。具体通过安装于电力牵引车接触网上方的受电弓在线拍照系统采集若干受电弓图像,形成受电弓图像学习样本集。其中,AdaBoost学习样本集包括正样本和负样本;主动形状模型学习样本集至少包含50张受电弓图像,且各受电弓图像的分辨率保持一致。
所述步骤20中AdaBoost级联分类器训练,具体步骤是:
1、收集受电弓正样本和负样本。正样本中必须包含待检测的受电弓目标,且正样本归一化为同一大小;负样本中不能包含受电弓,图像尺寸大于正样本;
2、Harr-like特征选择与计算。选择边缘(Edge)和线性(Line)的扩展Harr-like特征作为检测受电弓的特征;采用积分图(Intergral Image)快速计算特征;
3、强分类器训练。基于步骤10所述的AdaBoost正样本和负样本集,采用AdaBoost算法训练多个弱分类器,并线性加权组成强分类器,具体训练过程描述如为:
给定受电弓正样本和负样本集(x1,y1),(x2,y2),…,(xi,yi),…,(xn,yn),其中xi表示第i个样本,yi=0时表示其为负样本(非受电弓),yi=1时表示其为正样本(受电弓)。n为正样本和负样本总数。
初始化权重。对于正样本:wi=1/2l,l为正样本个数;对于负样本:wi=1/2m,m为负样本个数;
当t=1,...,T,进行如下迭代处理:
(1)权重归一化:
(2)对每个特征f,训练一个弱分类器h(xi,f,p,θ),并计算该弱分类器的加权错误率,找出一个具有最小错误率εt的弱分类器作为最优弱分类器:
(3)权重更新:
其中,ei=0表示样本xi检测正确,ei=1表示样本检测错误。
最终的强分类器:
其中,
4、级联分类器训练。训练多个强分类器,组成级联分类器。具体步骤是:
(1)确定每层,即每个强分类器的最大误识率(false positive rate)f,最小要达到的检测率(detection rate)d;
(2)指定最终级联分类器的误识率为Ftarget;
设P代表正样本集,N代表负样本集,初始误识率F0=1.0,初始检测率D0=1.0,i=0;
(3)当Fi>Ftarget,进行如下处理:
i=i+1;
ni=0;Fi=Fi-1;
当Fi>f×Fi-1,进行如下处理:
●ni=ni+1;
●利用AdaBoost算法在P和N上训练具有ni个特征数(弱分类器)的强分类器;
●衡量当前级联分类器误识率Fi和检测率Di;
●低第i层的强分类器阈值,直到当前级联分类器的检测率达到d×Di-1;
将负样本集N设置为空集;
利用当前的级联分类器检测受电弓负样本图像,将误识的图像放入负样本集N。
所述步骤30采用AdaBoost级联分类器对新获取、待识别的图像进行受电弓检测,具体采用多尺度检测机制对受电弓进行检测,步骤是:
1、计算待检测图像的积分图,便于快速计算各个检测窗口的特征;
2、初始化检测窗口,大小与级联分类器训练采用的正样本一致,对输入图像进行多种规模的、各个位置的窗口检测;
3、等比放大检测窗口,以检测不同大小的子图像;
4、区域合并:将同一受电弓在不同的尺度和邻近的位置上检测到的多次结果进行合并,得到合适的受电弓位置。
所述步骤40中构建受电弓主动形状模型,具体步骤是:
采用形态学滤波运算方法对步骤10采集的受电弓图像逐一进行图像降噪预处理;
在每张图像中选取受电弓轮廓的角点、边界点作为特征点,通过人工手动的方式,对受电弓特征点进行标记,且每张受电弓图像中的特征标记点必须相互对应、数量保持一致;
将标记完成的受电弓特征点的图像坐标保存于文本文件中,不同受电弓图像的特征点分开保存,且受电弓图像文件名称与对应的特征点文件名称保持一致。图像坐标指以图像左上角为坐标原点,水平向右方向为X轴,垂直向下方向为Y轴的坐标体系的像素坐标。此时,所有受电弓图像可以采用特征点的点分布模型(Pointdistribution model,PDM)对受电弓形状进行描述,即受电弓图像i的形状可以通过其所有特征点数学表示为:N为受电弓图像的特征点总数。受电弓图像学习样本集可以表示为:M为受电弓图像总数目;
基于受电弓图像和相应的特征点,根据主动形状模型(Active Shape Model,ASM)算法原理,建立受电弓主动形状模型。具体步骤包括:
1、采用Generalized Procrustes analysis(GPA)方法对所有受电弓图像的点分布模型进行对齐(Align)。对齐步骤如下:
(1)对于受电弓形状xi,i=2,3,…,M,逐个进行旋转、缩放和平移,并与形状x1做对齐,得到变换后的形状集合
(2)计算变换后的所有受电弓图像形状的平均值m;
其中:
(3)将平均形状m进行旋转、缩放和平移,与样本做对齐;
(4)将进行旋转、缩放和平移,与调整后的平均形状进行对齐匹配;
(5)如果平均形状收敛,停止。否则跳至第(2)步。
最终收敛的判定是依据重新对齐的各个受电弓形状与平均形状之间的差别,即寻找到变换(Ti),使得下式取值最小。
∑|m-Ti(xi)|2
所述受电弓形状之间通过旋转、缩放和平移进行对齐描述为:以两个受电弓形状为例,每个形状有N个坐标对:
首先定义一个变换T,T是由一个平移(tx,ty),旋转θ和尺度s组成,则形状x2的变换可以表示为:
设
将形状x2经过变换T对齐到x1,最佳的变换可以通过最小化下面的表达式得到:
E=[x1-Rx2-(tx,ty)T]T[x1-Rx2-(tx,ty)T]
这个最小化是最小二乘法的常规应用:E对未知变量θ,s,tx,ty的偏微分可以计算出来,并令它们为零,联立这些方程,既可求解得变换T。
2、ASM模型建立:假设目前已经存在经过对齐处理的M个训练形状每个形状可以由N对坐标给出:平均形状设为:然后计算协方差矩阵(2N×2N维):训练形状在某些方向上的变化是描述受电弓形状的重要性质,可以从协方差矩阵S的特征向量中得到,亦即,求解线性方程:Spi=λipi。协方差矩阵S的特征向量为(对应的特征值已按降序排列):
P=(p1p2p3…p2N)
则对于任何向量X,存在向量b(形状模型参数),满足
也可表示为:
特征值较大的特征向量描述了训练形状变化最大的方向,在描述“合理”的形状与平均形状偏差有多大时,p2N,p2N-1,…方向上的贡献是微不足道的。因此可以设:
Pt=(p1p2p3…pt) t≤2N
bt=(b1,b2,…,bt)T
可以得到估计
如果X是与训练集相关的合理形状,对于足够大的t,该估计可以很好地拟合真实形状。
向量bt定义了一组可变模型参数,不同的bt可以拟合出不同变化的形状。研究表明,bi在训练集上的方差将会与特征值λi相关;相应地,对于较好的形状,bi通常要求满足下面的不等式
所述步骤50中结合受电弓检测结果和主动形状模型的受电弓精确匹配和识别,具体采用单分辨率搜索算法精确匹配受电弓形状,具体步骤包括:
1、根据步骤40建立受电弓主动形状模型过程中生成的平均形状和待检测图像中受电工的初始位置,初始化受电弓形状,表示如下:
2、在初始化受电弓形状的每一个标记点处,在边界附近沿边界的法向进行搜索,确定具有最高梯度的像素点,给该点打上最佳目标位置的标志,将标记点向这个标志点移动,如果没有明显的新目标点,标记点位置不作移动。
3、通过上述标记点移动后,形状发生了改变,发生改变的形状与初始化受电弓形状之间存在一个位移向量由上述知:发生位移后可以表示为:
可得进一步推导得
4、重复第2、3步,直到姿态参数变化可以忽略不计。
所述步骤60中受电弓定量检测和分析,指在精确匹配和识别结果基础上对受电弓碳滑板厚度进行定量检测和分析。具体步骤是:
根据受电弓精确匹配结果,获取受电弓厚度;
由于受电弓支架厚度与受电弓类型相关,且对于特定的受电弓类型,受电弓支架厚度是固定已知的,可作为先验知识。因此,碳滑板厚度可以通过受电弓厚度与受电弓支架厚度相减得到。
(三)技术效果
本发明与现有的技术方案相比具有如下的优点及有益效果:本发明对受电弓进行图像检测仅需一张图像,不必要求拍照系统对同一受电弓成像多次;由于受电弓图像具有复杂多变的背景,不能简单地通过边缘提取等常规方法对受电弓进行定量检测。本发明提出的基于AdaBoost和主动形状模型的受电弓自动识别方法,事先通过受电弓形状的学习,结合受电弓检测结果,能保证受电弓在运行状态下有较大角度、大小变化情况下的检测准确率。实验结果表明:本发明提出的基于AdaBoost和主动形状模型的受电弓自动识别方法准确率比边缘提取等常规方法明显提高,且能更好地为后续的受电弓的定量检测分析打下基础。
附图说明
图1是本发明实施例的基于AdaBoost和主动形状模型的受电弓自动识别方法流程图
图2是部分AdaBoost学习样本集(左为正样本,右为负样本)
图3是部分主动形状模型学习样本集
图4是检测受电弓的Harr-like特征
图5是受电弓检测结果
图6是受电弓特征点标记示意图
图7是ASM单分辨率搜索算法示意图
图8是受电弓厚度定量检测结果.
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
如图1所示,是本发明实施例的基于AdaBoost和主动形状模型的受电弓自动识别方法流程图,本实施例包括如下步骤:
步骤10:获取受电弓图像学习样本集;
本实施例受电弓图像学习样本集包括两套:AdaBoost学习样本集:用于训练AdaBoost级联分类器;主动形状模型学习样本集:用于训练受电弓主动形状模型。具体通过安装于电力牵引车接触网上方的受电弓在线拍照系统采集若干受电弓图像,形成受电弓图像学习样本集。其中,AdaBoost学习样本集包括正样本和负样本,正样本数目为756,负样本为3215,图2是部分AdaBoost学习样本集;主动形状模型学习样本集至少包含50张受电弓图像,且各受电弓图像的分辨率保持一致,图3是部分主动形状模型学习样本集。
步骤20、AdaBoost级联分类器训练;
本实例AdaBoost级联分类器训练步骤是:
21:收集受电弓正样本和负样本。正样本中必须包含待检测的受电弓目标,且正样本归一化为同一大小,本实例归一化为80×25;负样本中不能包含受电弓,图像尺寸大于正样本;
22:Harr-like特征选择与计算。选择边缘(Edge)和线性(Line)的扩展Harr-like特征作为检测受电弓的特征,如图4所示;采用积分图(Intergral Image)快速计算特征;
23:强分类器训练。基于步骤10所述的AdaBoost正样本和负样本集,采用AdaBoost算法训练多个弱分类器,并线性加权组成强分类器,具体训练过程描述如为:
231:给定受电弓正样本和负样本集(x1,y1),(x2,y2),…,(xi,yi),…,(xn,yn),其中xi表示第i个样本,yi=0时表示其为负样本(非受电弓),yi=1时表示其为正样本(受电弓)。n为正样本和负样本总数;
232:初始化权重。对于正样本:wi=1/2l,l为正样本个数;对于负样本:wi=1/2m,m为负样本个数;
233:当t=1,...,T,进行如下迭代处理:
2331:权重归一化:
2332:对每个特征f,训练一个弱分类器h(xi,f,p,θ),并计算该弱分类器的加权错误率,找出一个具有最小错误率εt的弱分类器作为最优弱分类器:
2333:权重更新:
其中,ei=0表示样本xi检测正确,ei=1表示样本检测错误。
234:最终的强分类器:
其中,
24:级联分类器训练。训练多个强分类器,组成级联分类器。具体步骤是:
241:确定每层,即每个强分类器的最大误识率(false positive rate)f,最小要达到的检测率(detection rate)d;
242:指定最终级联分类器的误识率为Ftarget;
设P代表正样本集,N代表负样本集,初始误识率F0=1.0,初始检测率D0=1.0,i=0;
243:当Fi>Ftarget,进行如下处理:
2431:i=i+1;
2432:ni=0;Fi=Fi-1;
2433:当Fi>f×Fi-1,进行如下处理:
24331:ni=ni+1;
24332:利用AdaBoost算法在P和N上训练具有ni个特征数(弱分类器)的强分类器;
24333:衡量当前级联分类器误识率Fi和检测率Di;
24334:低第i层的强分类器阈值,直到当前级联分类器的检测率达到d×Di-1;
2434:将负样本集N设置为空集;
2435:利用当前的级联分类器检测受电弓负样本图像,将误识的图像放入负样本集N。
步骤30、采用AdaBoost级联分类器对新获取、待识别的图像进行受电弓检测;
本实例采用多尺度检测机制对受电弓进行检测,步骤是:
31:计算待检测图像的积分图,便于快速计算各个检测窗口的特征;
32:初始化检测窗口,大小与级联分类器训练采用的正样本一致,对输入图像进行多种规模的、各个位置的窗口检测;
33:等比放大检测窗口,以检测不同大小的子图像;
34:区域合并:将同一受电弓在不同的尺度和邻近的位置上检测到的多次结果进行合并,得到合适的受电弓位置,如图5所示。
步骤40:构建受电弓主动形状模型;
本实施例构建受电弓主动形状模型步骤是:
41:采用形态学滤波运算方法对步骤10采集的受电弓图像逐一进行图像降噪预处理;
42:在每张图像中选取受电弓轮廓的角点、边界点作为特征点,通过人工手动的方式,对受电弓特征点进行标记,且每张受电弓图像中的特征标记点必须相互对应、数量保持一致;
43:将标记完成的受电弓特征点的图像坐标保存于文本文件中,不同受电弓图像的特征点分开保存,且受电弓图像文件名称与对应的特征点文件名称保持一致,标记结果如图6所示。图像坐标指以图像左上角为坐标原点,水平向右方向为X轴,垂直向下方向为Y轴的坐标体系的像素坐标。此时,所有受电弓图像可以采用特征点的点分布模型(Pointdistribution model,PDM)对受电弓形状进行描述,即受电弓图像i的形状可以通过其所有特征点数学表示为:N为受电弓图像的特征点总数。受电弓图像学习样本集可以表示为:M为受电弓图像总数目;
44:基于受电弓图像和相应的特征点,根据主动形状模型(Active Shape Model,ASM)算法原理,建立受电弓主动形状模型。具体步骤包括:
441:采用Generalized Procrustes analysis(GPA)方法对所有受电弓图像的点分布模型进行对齐(Align)。对齐步骤如下:
4411:对于受电弓形状xi,i=2,3,…,M,逐个进行旋转、缩放和平移,并与形状x1做对齐,得到变换后的形状集合
4412:计算变换后的所有受电弓图像形状的平均值m;
其中:
4413:将平均形状m进行旋转、缩放和平移,与样本做对齐;
4414:将进行旋转、缩放和平移,与调整后的平均形状进行对齐匹配;
4415:如果平均形状收敛,停止。否则跳至第4412步。
最终收敛的判定是依据重新对齐的各个受电弓形状与平均形状之间的差别,即寻找到变换(Ti),使得下式取值最小。
∑|m-Ti(xi)|2
所述受电弓形状之间通过旋转、缩放和平移进行对齐描述为:以两个受电弓形状为例,每个形状有N个坐标对:
首先定义一个变换T,T是由一个平移(tx,ty),旋转θ和尺度s组成,则形状x2的变换可以表示为:
设
将形状x2经过变换T对齐到x1,最佳的变换可以通过最小化下面的表达式得到:
E=[x1-Rx2-(tx,ty)T]T[x1-Rx2-(tx,ty)T]
这个最小化是最小二乘法的常规应用:E对未知变量θ,s,tx,ty的偏微分可以计算出来,并令它们为零,联立这些方程,既可求解得变换T。
442:ASM模型建立:假设目前已经存在经过对齐处理的M个训练形状每个形状可以由N对坐标给出:平均形状设为:然后计算协方差矩阵(2N×2N维):训练形状在某些方向上的变化是描述受电弓形状的重要性质,可以从协方差矩阵S的特征向量中得到,亦即,求解线性方程:Spi=λipi。协方差矩阵S的特征向量为(对应的特征值已按降序排列):
P=(p1p2p3…p2N)
则对于任何向量X,存在向量b(形状模型参数),满足
也可表示为:
特征值较大的特征向量描述了训练形状变化最大的方向,在描述“合理”的形状与平均形状偏差有多大时,p2N,p2N-1,…方向上的贡献是微不足道的。因此可以设:
Pt=(p1p2p3…pt) t≤2N
bt=(b1,b2,…,bt)T
可以得到估计
如果X是与训练集相关的合理形状,对于足够大的t,该估计可以很好地拟合真实形状。
向量bt定义了一组可变模型参数,不同的bt可以拟合出不同变化的形状。研究表明,bi在训练集上的方差将会与特征值λi相关;相应地,对于较好的形状,bi通常要求满足下面的不等式
步骤50:结合检测结果和主动形状模型的受电弓精确匹配和识别;
本实施例结合受电弓检测结果和主动形状模型进行受电弓精确匹配和识别,具体采用单分辨率搜索算法精确匹配受电弓形状,具体步骤包括:
51:根据步骤40建立受电弓主动形状模型过程中生成的平均形状和待检测图像中受电工的初始位置,初始化受电弓形状,表示如下:
52:在初始化受电弓形状的每一个标记点处,在边界附近沿边界的法向进行搜索,确定具有最高梯度的像素点,给该点打上最佳目标位置的标志,将标记点向这个标志点移动,如果没有明显的新目标点,标记点位置不作移动,算法示意如图7所示。
53:通过上述标记点移动后,形状发生了改变,发生改变的形状与初始化受电弓形状之间存在一个位移向量由上述知:发生位移后可以表示为:
可得进一步推导得
54:重复第52、53步,直到姿态参数变化可以忽略不计。
步骤60:受电弓定量检测和分析。
本实例在精确匹配和识别结果基础上对受电弓碳滑板厚度进行定量检测和分析。具体步骤是:
61:根据受电弓精确匹配结果,获取受电弓厚度;
62:由于受电弓支架厚度与受电弓类型相关,且对于特定的受电弓类型,受电弓支架厚度是固定已知的,可作为先验知识。因此,碳滑板厚度可以通过受电弓厚度与受电弓支架厚度相减得到,如图8所示。
最后,通过设定受电弓碳滑板厚度磨损的标准,与检测结果进行对比,若判断检测受电弓碳滑板磨损超出标准,则发出报警并提示更换受电弓。
实验结果表明,通过本技术方案,可以对复杂背景下、不同角度和不同尺度变化的受电弓作出更加精确的定位和识别。相比现有的相关算法,本发明在保证精确率的同时,增强了算法的鲁棒性,实现了受电弓的线上自动检测。
本发明跳出了传统多图像、安装检测装置的受电弓在线检测设计思路,提出一种全新的基于AdaBoost和主动形状模型的受电弓自动识别方法。该方法只需在线获取一张受电弓图像,首先通过AdaBoost级联分类器检测受电弓,然后结合检测结果,加入基于主动形状模型的学习算法,从而能够自适应地把具有形状可变性的受电弓精确检测出来,可更加快速地进行车辆维护,节约受电弓检测成本。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变型,这些改进和变型也应视为本发明的保护范围。
Claims (5)
1.一种基于AdaBoost和主动形状模型的受电弓识别方法,其特征在于,包括:
步骤10、获取受电弓图像学习样本集;
步骤20、AdaBoost级联分类器训练;
步骤30、采用AdaBoost级联分类器对新获取、待识别的图像进行受电弓检测;
步骤40、构建受电弓主动形状模型;
步骤50、结合检测结果和主动形状模型的受电弓精确匹配和识别;
步骤60、受电弓定量检测和分析,
所述步骤10中获取受电弓图像学习样本集包括两套:AdaBoost学习样本集:用于训练AdaBoost级联分类器;主动形状模型学习样本集:用于训练受电弓主动形状模型,具体通过安装于电力牵引车接触网上方的受电弓在线拍照系统采集若干受电弓图像,形成受电弓图像学习样本集,其中,AdaBoost学习样本集包括正样本和负样本;主动形状模型学习样本集至少包含50张受电弓图像,且各受电弓图像的分辨率保持一致;
所述步骤20中AdaBoost级联分类器训练,具体步骤是:
一、收集受电弓正样本和负样本,正样本中必须包含待检测的受电弓目标,且正样本归一化为同一大小;负样本中不能包含受电弓,图像尺寸大于正样本;
二、Harr-like特征选择与计算,选择边缘和线性的扩展Harr-like特征作为检测受电弓的特征;采用积分图快速计算特征;
三、强分类器训练,基于步骤10所述的AdaBoost正样本和负样本集,采用AdaBoost算法训练多个弱分类器,并线性加权组成强分类器,具体训练过程描述如为:
给定受电弓正样本和负样本集(x1,y1),(x2,y2),…,(xi,yi),…,(xn,yn),其中xi表示第i个样本,yi=0时表示其为负样本,yi=1时表示其为正样本,n为正样本和负样本总数;
初始化权重,对于正样本:wi=1/2l,l为正样本个数;对于负样本:wi=1/2m,m为负样本个数;
当t=1,...,T,进行如下迭代处理:
(1)权重归一化:
(2)对每个特征f,训练一个弱分类器h(xi,f,p,θ),并计算该弱分类器的加权错误率,找出一个具有最小错误率εt的弱分类器作为最优弱分类器:
(3)权重更新:
其中,ei=0表示样本xi检测正确,ei=1表示样本检测错误;
最终的强分类器:
其中,
四、级联分类器训练,训练多个强分类器,组成级联分类器,具体步骤是:
(1)确定每层,即每个强分类器的最大误识率f,最小要达到的检测率d;
(2)指定最终级联分类器的误识率为Ftarget;
设P代表正样本集,N代表负样本集,初始误识率F0=1.0,初始检测率D0=1.0,i=0;
(3)当Fi>Ftarget,进行如下处理:
i=i+1;
ni=0;Fi=Fi-1;
当Fi>f×Fi-1,进行如下处理:
ni=ni+1;
利用AdaBoost算法在P和N上训练具有ni个特征数的强分类器;
衡量当前级联分类器误识率Fi和检测率Di;
低第i层的强分类器阈值,直到当前级联分类器的检测率达到d×Di-1;
将负样本集N设置为空集;
利用当前的级联分类器检测受电弓负样本图像,将误识的图像放入负样本集N。
2.如权利要求1所述的基于AdaBoost和主动形状模型的受电弓识别方法,其特征在于,所述步骤30采用AdaBoost级联分类器对新获取、待识别的图像进行受电弓检测,具体采用多尺度检测机制对受电弓进行检测,步骤是:
计算待检测图像的积分图,便于快速计算各个检测窗口的特征;
初始化检测窗口,大小与级联分类器训练采用的正样本一致,对输入图像进行多种规模的、各个位置的窗口检测;
等比放大检测窗口,以检测不同大小的子图像;
区域合并:将同一受电弓在不同的尺度和邻近的位置上检测到的多次结果进行合并,得到合适的受电弓位置。
3.如权利要求1所述的基于AdaBoost和主动形状模型的受电弓识别方法,其特征在于,所述步骤40中构建受电弓主动形状模型,具体步骤是:
采用形态学滤波运算方法对步骤10采集的主动形状模型学习样本集逐一进行图像降噪预处理;
在每张图像中选取受电弓轮廓的角点、边界点作为特征点,通过人工手动的方式,对受电弓特征点进行标记,且每张受电弓图像中的特征标记点必须相互对应、数量保持一致;
将标记完成的受电弓特征点的图像坐标保存于文本文件中,不同受电弓图像的特征点分开保存,且受电弓图像文件名称与对应的特征点文件名称保持一致,图像坐标指以图像左上角为坐标原点,水平向右方向为X轴,垂直向下方向为Y轴的坐标体系的像素坐标;
基于受电弓图像和相应的特征点,根据主动形状模型算法原理,建立受电弓主动形状模型,该模型以文件的形式存在。
4.如权利要求1所述的基于AdaBoost和主动形状模型的受电弓识别方法,其特征在于,所述步骤50中结合受电弓检测结果和主动形状模型的受电弓精确匹配和识别,具体采用单分辨率搜索算法精确匹配受电弓形状。
5.如权利要求1所述的基于AdaBoost和主动形状模型的受电弓识别方法,其特征在于,所述步骤60中受电弓定量检测和分析,指在精确匹配和识别结果基础上对受电弓碳滑板厚度进行定量检测和分析,具体步骤是:
根据受电弓精确匹配结果,获取受电弓厚度;
由于受电弓支架厚度与受电弓类型相关,且对于特定的受电弓类型,受电弓支架厚度是固定已知的,可作为先验知识, 因此,碳滑板厚度可以通过受电弓厚度与受电弓支架厚度相减得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310566693.6A CN103745238B (zh) | 2013-11-15 | 2013-11-15 | 基于AdaBoost和主动形状模型的受电弓识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310566693.6A CN103745238B (zh) | 2013-11-15 | 2013-11-15 | 基于AdaBoost和主动形状模型的受电弓识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103745238A CN103745238A (zh) | 2014-04-23 |
CN103745238B true CN103745238B (zh) | 2018-04-24 |
Family
ID=50502255
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310566693.6A Active CN103745238B (zh) | 2013-11-15 | 2013-11-15 | 基于AdaBoost和主动形状模型的受电弓识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103745238B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105550710A (zh) * | 2015-12-17 | 2016-05-04 | 西安电子科技大学 | 基于非线性拟合的接触网运行异常状态智能检测方法 |
CN108694349B (zh) * | 2017-04-07 | 2020-03-06 | 成都唐源电气股份有限公司 | 一种基于线阵相机的受电弓图像提取方法及装置 |
CN107633266B (zh) * | 2017-09-07 | 2020-07-28 | 西安交通大学 | 一种电力机车接触网受电弓电弧检测方法 |
CN108009574B (zh) * | 2017-11-27 | 2022-04-29 | 成都明崛科技有限公司 | 一种轨道扣件检测方法 |
CN108229337B (zh) * | 2017-12-14 | 2021-03-30 | 创新先进技术有限公司 | 一种数据处理的方法、装置及设备 |
CN110095984A (zh) * | 2019-04-23 | 2019-08-06 | 西南交通大学 | 一种基于状态估计的高速铁路受电弓多目标鲁棒控制方法 |
CN113324864B (zh) * | 2020-02-28 | 2022-09-20 | 南京理工大学 | 一种基于深度学习目标检测的受电弓碳滑板磨耗检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009098598A2 (en) * | 2008-02-08 | 2009-08-13 | Administrador De Infraestructuras Ferroviarias (Adif) | Detection system for the state of train pantographs |
CN102129565A (zh) * | 2011-03-01 | 2011-07-20 | 北京航空航天大学 | 基于冗余特征消减AdaBoost分类器的物体检测方法 |
-
2013
- 2013-11-15 CN CN201310566693.6A patent/CN103745238B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009098598A2 (en) * | 2008-02-08 | 2009-08-13 | Administrador De Infraestructuras Ferroviarias (Adif) | Detection system for the state of train pantographs |
CN102129565A (zh) * | 2011-03-01 | 2011-07-20 | 北京航空航天大学 | 基于冗余特征消减AdaBoost分类器的物体检测方法 |
Non-Patent Citations (3)
Title |
---|
《基于AdaBoost和ASM算法的人眼定位》;郑玉铎等;《北京信息科技大学学报》;20130630;第28卷(第3期);全文 * |
基于主动形状模型的受电弓自动检侧方法;岳安志等;《第十三届中国体视学与图像分析学术会议论文集》;20130930;第313页第一栏第3段、第二栏第3段,图1,第314页第一栏第2段、第二栏第1-2段、图2-3,第315页,第316页第二栏第1段、第6-7段 * |
基于学习算法受电弓瓷瓶定位检测方法的研究;王士成;《中国优秀硕士学位论文全文数据库 信息科技辑》;20120415;第38页第2段,第39页第图4-2,第40页第2-3段、图4-3,第45页第2段,第46页第2段, 第47页第1段,第48页第1段、第3段,第49页第1段,图5-7 * |
Also Published As
Publication number | Publication date |
---|---|
CN103745238A (zh) | 2014-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103745238B (zh) | 基于AdaBoost和主动形状模型的受电弓识别方法 | |
CN103116756B (zh) | 一种人脸检测与跟踪方法及装置 | |
CN106682603B (zh) | 一种基于多源信息融合的实时驾驶员疲劳预警系统 | |
CN105046717B (zh) | 一种鲁棒性的视频目标对象跟踪方法 | |
CN103837087B (zh) | 基于主动形状模型的受电弓自动检测方法 | |
CN103279768B (zh) | 一种基于增量学习人脸分块视觉表征的视频人脸识别方法 | |
CN105844621A (zh) | 一种印刷品质量检测方法 | |
CN107564002A (zh) | 塑料管表面缺陷检测方法、系统及计算机可读存储介质 | |
Liu et al. | An improved surface roughness measurement method for micro-heterogeneous texture in deep hole based on gray-level co-occurrence matrix and support vector machine | |
CN105243386A (zh) | 人脸活体判断方法以及系统 | |
CN102254183B (zh) | 一种基于AdaBoost算法的人脸检测方法 | |
CN101726498B (zh) | 基于视觉仿生的铜带表面质量智能检测装置及方法 | |
CN106650805A (zh) | 一种视觉目标跟踪方法及装置 | |
CN108898623A (zh) | 目标跟踪方法及设备 | |
CN114037703B (zh) | 基于二维定位和三维姿态解算的地铁阀门状态检测方法 | |
CN110427943A (zh) | 一种基于r-cnn的智能电表数值识别方法 | |
CN115993096A (zh) | 一种高层建筑变形测量方法 | |
CN108647722A (zh) | 一种基于过程尺寸特征的锌矿品位软测量方法 | |
CN110751005B (zh) | 融合深度感知特征和核极限学习机的行人检测方法 | |
CN117288094B (zh) | 基于激光传感器的掘进机实时定位系统 | |
CN105825215A (zh) | 一种基于局部近邻嵌入核函数的仪表定位方法及其使用载体 | |
CN106683116B (zh) | 一种基于支持向量机的粒子滤波联合跟踪方法 | |
CN101021948A (zh) | 人体对称性运动图像中关节的自动识别装置及方法 | |
CN206019585U (zh) | 钢轨磨耗检测系统 | |
CN112329893A (zh) | 基于数据驱动的异源多目标智能检测方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |