基于码本和平滑转换自回归模型的黑烟车智能监控方法技术

本发明专利技术公开了一种基于码本和平滑转换自回归模型的黑烟车智能监控方法，包括如下步骤：(1)利用码本模型检测视频中的运动目标，确定目标后方关键区域；(2)提取关键区域的Tamura特征，中心对称局部二值模式直方图特征和灰度直方图特征，组合形成关键区域的三类静态特征；(3)利用平滑转换自回归模型对三类静态特征的时间序列分别进行建模，并将模型的解作为最终特征向量，用于刻画关键区域的动态特征；(4)三类特征训练得到三个SVM分类器，对三个分类结果加权融合，得到当前关键区域的识别结果，并通过连续多个关键区域识别结果的分析，对当前视频段是否有黑烟车做出判断。本发明专利技术能够提高道路监控视频中黑烟车的检测效率，降低误报率。

【技术实现步骤摘要】
基于码本和平滑转换自回归模型的黑烟车智能监控方法
本专利技术涉及图像特征提取和时序分析
，涉及一种基于码本和平滑转换自回归模型的黑烟车智能监控方法。
技术介绍
柴油货车的污染一直是机动车污染的重中之重，本专利技术讨论的黑烟车就常见于的柴油货车，如何及时发现道路上行驶的黑烟车，并根据污染程度实施淘汰或者强制维护，将非常有助于减少空气污染，打好柴油货车污染治理攻坚战，大幅降低污染物排放，赢得蓝天保卫战的胜利。目前部分地区采用车辆尾气分析装置检测黑烟车，但是还存在很多不足，一方面，车载尾气分析装置往往价格昂贵，且后续的维护和保养需要消耗大量的资金，由于车辆数量的增加，给每辆车都配置的可实施性比较困难；另一方面，对于路边安装的车辆尾气分析，检测结果往往受多种因素的影响，诸如，多车并行，跟车太近，恶劣天气(风速过大、雨雪天气)，道路环境背景，旁边车辆排放，排气管高低，设备状况(光强、噪音等)、操作人员的专业程度(安装调试)等等。
技术实现思路
本专利技术所要解决的技术问题在于，提出一种基于码本和平滑转换自回归模型的黑烟车智能监控方法，该方法能够进一步提高黑烟车检测的识别率，降低误报率。为解决上述技术问题，本专利技术提出一种基于码本模型和平滑转换自回归模型的黑烟车智能监控方法，包括如下步骤：(1)利用码本(Codebook)模型检测视频中的运动目标，确定目标后方关键区域；(2)提取关键区域的Tamura特征，中心对称局部二值模式(CS-LBP)直方图特征和灰度直方图特征，组合形成关键区域的三类静态特征；(3)利用平滑转换自回归(STAR)模型对三类静态特征的时间序列分别进行建模，并将模型的解作为最终特征向量，用于刻画关键区域的动态特征；(4)三类特征训练得到三个SVM分类器，对三个分类结果加权融合，得到当前关键区域的识别结果，并通过连续多个关键区域识别结果的分析，对当前视频段是否有黑烟车做出判断。该方法提出了一种基于码本和平滑转换自回归模型的黑烟车智能监控方法，该方法利用码本模型检测运动目标，融合Tamura特征，中心对称局部二值模式(CS-LBP)直方图特征和灰度直方图特征，并结合平滑转换自回归模型，对车辆关键区域的时序特征做出刻画，该方法能大幅度提高监控黑烟车的效率。作为本专利技术的一种改进，步骤(1)中的利用码本(Codebook)模型检测视频中的运动目标包括如下几点：(11)码本模型初始化及模型特性记C＝(c1，c2，...cL)为一个像素的码本，在初始化时，将每个像素的码本置空，即L＝0，这里的L表示每个码本中码字的个数，每个码字ci(i＝1，2，...，L.)由灰度值vi和六元组auxi两部分组成，其中，和分别表示码字中最小和最大灰度值，fi表示码字出现的频率，λi表示训练中该码字没有出现的最大时间间隔，pi和qi分别表示该码字第一次出现和最后一次出现的时间；(12)码本建模和模型更新i)设X＝{x1，x2，...，xN}是一像素的序列采样值，xt是当前像素t时刻的灰度值，对于每个像素的采样值中的每个xt，利用下述条件找出与其匹配的码字cm，其中，[Ilow，Ihigh]表示每个码字的灰度值变化范围，ε1是一个阈值，vm表示码字cm对应的灰度值；ii)如果码本为空或者找不到匹配的码字，则创建一个新的码字，如果有码字cm满足上述两个条件，把码字更新为iii)训练结束后，计算该像素每个码字没有再次出现的最大时间间隔，对于每个码字ci(i＝1，2，...，L.)，λi＝max{λi，(N-qi+pi-1)}，利用λ消除冗余码字，得到最能代表真实背景的初始码本M＝{ck|ck∈C，λk≤Tm}，Tm通常取训练帧数的一半；(13)前景目标提取对于新输入像素xt，如果在码本M中能找出满足下列条件的匹配码字cm，则xt为前景运动目标像素，否则为背景像素，其中，[Ilow，Ihigh]表示每个码字的灰度值变化范围，ε2是一个阈值，vm表示码字cm中对应的灰度值；作为本专利技术的一种改进，步骤(2)中的关键区域的Tamura特征包括以下几种：(21)特征一：粗糙度(COA)，其计算方法为，其中，m和n分别表示图像的宽和高，{E1，E2，...，Ek，...，EL}表示在k取遍集合{1，2，3，4，5}中的数值时，位置(i，j)处的水平平均强度差Ek，h(i，j)和垂直平均强度差Ek，v(i，j)的数值的集合，Ek，v(i，j)＝|Ak(i，j+2k-1)-Ak(i，j-2k-1)|Ek，h(i，j)＝|Ak(i+2k-1，j)-Ak(i-2k-1，j)|其中，F(i，j)表示位置(i，j)处的灰度值，Ak(x1，x2)表示活动窗口为2k×2k(k＝1，2，3，4，5)的像素的平均值。(22)特征二：对比度(CON)，其计算方法为，其中，σ表示标准差，μ4均值的四阶矩，参数n0一般设置为0.25。(23)特征三：方向度(DIR)，其计算方法为，其中，φ表示量化方向码，p表示直方图峰值，对任意一个峰值p，Wp表示该峰值包含的量化范围，np表示直方图中峰值的数目，φp表示Wp中最大直方图中的量化数值，HD表示直方图，通过计算梯度幅值|ΔG|大于阈值tho且角度θ在某个范围的点的个数得到，HD(k)表示点的幅值|ΔG|≥tho且角度θ满足(2k-1)π/2n≤θ≤(2k+1)π/2n的点的个数，θ和|ΔG|的计算方法为θ＝tan-1(ΔV/ΔH)+π/2|ΔG|＝(|ΔH|+|ΔV|)/2其中，ΔV和ΔH通过一个3×3的Prewitt算子得到的垂直差异值和水平差异值；(24)特征四：线性度(LIN)，其计算方法为，其中，PDd(i，j)表示n×n的局部方向共生矩阵(DCOM)在点(i，j)的元素值，该点的距离为d且沿着边界具有方向码i和方向码j，我们利用微分和阈值化计算HD和PDd，局部方向共生矩阵被定义为两个相邻单元的相对频率，通过沿着图像便捷的距离d来分离，与原方向相同则数值+1，与原方向垂直则数值-1；(25)特征五：规整度(REG)，其计算方法为，REG＝1-r(σCON+σCOR+σDIR+σLIN)其中，σXXX对应特征XXX的标准差，r是一个规范化因子；(26)特征六：粗略度(ROU)，其计算方法为，ROU＝COA+CON其中，COA和CON分别表示粗糙度和对比度；将第一类特征Tamura特征记作FTamura，即FTamura＝{CON，COR，DIR，LIN，REG，ROU}。步骤(2)中的关键区域的中心对称局部二值模式(CS-LBP)直方图特征的计算包括以下步骤：(27)中心对称局部二值模式的计算方式为，其中，R表示半径大小，N表示邻域所取像素个数，T表示阈值，CSLBPR，N，T(x，y)表示半径为R邻域点个数为N的中心对称局部二值模式特征在位置(x，y)处的数值；(28)对得到的中心对称局部二值模式图提取直方图特征；将第二类特征记作FCSLBP，即(28)提取的直方图特征。步骤(2)中的关键区域的灰度直方图特征包括以下步骤：(29)将关键区域图像转化为灰度图像，并做一个灰度值截断处理，即凡是大于阈值Th的灰度值全部赋值为Ttru，其中，Inorm(x，y)表示规范化关键区域在位置(x，y)处的灰度值；(210)本文档来自技高网...

【技术保护点】
1.一种基于码本和平滑转换自回归模型的黑烟车智能监控方法，其特征在于，包括如下步骤：(1)利用码本(Codebook)模型检测视频中的运动目标，确定目标后方关键区域；(2)提取关键区域的Tamura特征，中心对称局部二值模式(CS‑LBP)直方图特征和灰度直方图特征，组合形成关键区域的三类静态特征；(3)利用平滑转换自回归(STAR)模型对三类静态特征的时间序列分别进行建模，并将模型的解作为最终特征向量，用于刻画关键区域的动态特征；(4)三类特征训练得到三个SVM分类器，对三个分类结果加权融合，得到当前关键区域的识别结果，并通过连续多个关键区域识别结果的分析，对当前视频段是否有黑烟车做出判断。

【技术特征摘要】
1.一种基于码本和平滑转换自回归模型的黑烟车智能监控方法，其特征在于，包括如下步骤：(1)利用码本(Codebook)模型检测视频中的运动目标，确定目标后方关键区域；(2)提取关键区域的Tamura特征，中心对称局部二值模式(CS-LBP)直方图特征和灰度直方图特征，组合形成关键区域的三类静态特征；(3)利用平滑转换自回归(STAR)模型对三类静态特征的时间序列分别进行建模，并将模型的解作为最终特征向量，用于刻画关键区域的动态特征；(4)三类特征训练得到三个SVM分类器，对三个分类结果加权融合，得到当前关键区域的识别结果，并通过连续多个关键区域识别结果的分析，对当前视频段是否有黑烟车做出判断。2.如权利要求1所述的基于码本模型和平滑转换自回归模型的黑烟车智能监控方法，其特征在于，步骤(1)中的利用码本(Codebook)模型检测视频中的运动目标包括如下几点：(11)码本模型初始化及模型特性,记C＝(c1,c2,...cL)为一个像素的码本，在初始化时，将每个像素的码本置空，即L＝0，这里的L表示每个码本中码字的个数，每个码字ci(i＝1,2,...,L.)由灰度值vi和六元组auxi两部分组成，其中，和分别表示码字中最小和最大灰度值，fi表示码字出现的频率，λi表示训练中该码字没有出现的最大时间间隔，pi和qi分别表示该码字第一次出现和最后一次出现的时间；(12)码本建模和模型更新，i)设X＝{x1,x2,...,xN}是一像素的序列采样值，xt是当前像素t时刻的灰度值，对于每个像素的采样值中的每个xt，利用下述条件找出与其匹配的码字cm，其中，[Ilow,Ihigh]表示每个码字的灰度值变化范围，ε1是一个阈值，vm表示码字cm对应的灰度值；ii)如果码本为空或者找不到匹配的码字，则创建一个新的码字，如果有码字cm满足上述两个条件，把码字更新为iii)训练结束后，计算该像素每个码字没有再次出现的最大时间间隔，对于每个码字ci(i＝1,2,...,L.),λi＝max{λi,(N-qi+pi-1)}，利用λ消除冗余码字，得到最能代表真实背景的初始码本M＝{ck|ck∈C,λk≤Tm}，Tm通常取训练帧数的一半；(13)前景目标提取，对于新输入像素xt，如果在码本M中能找出满足下列条件的匹配码字cm，则xt为前景运动目标像素，否则为背景像素，其中，[Ilow,Ihigh]表示每个码字的灰度值变化范围，ε2是一个阈值，vm表示码字cm中对应的灰度值。3.如权利要求1所述的基于码本模型和平滑转换自回归模型的黑烟车智能监控方法，其特征在于，步骤(2)中的关键区域的Tamura特征包括以下几种：(21)特征一：粗糙度(COA)，其计算方法为，其中，m和n分别表示图像的宽和高，{E1,E2,...,Ek,...,EL}表示在k取遍集合{1,2,3,4,5}中的数值时，位置(i,j)处的水平平均强度差Ek,h(i,j)和垂直平均强度差Ek,v(i,j)的数值的集合，Ek,v(i,j)＝|Ak(i,j+2k-1)-Ak(i,j-2k-1)|；Ek,h(i,j)＝|Ak(i+2k-1,j)-Ak(i-2k-1,j)|；其中，F(i,j)表示位置(i,j)处的灰度值，Ak(x1,x2)表示活动窗口为2k×2k(k＝1,2,3,4,5)的像素的平均值；(22)特征二：对比度(CON)，其计算方法为，其中，σ表示标准差，μ4均值的四阶矩，参数n0设置为0.25；(23)特征三：方向度(DIR)，其计算方法为，其中，φ表示量化方向码，p表示直方图峰值，对任意一个峰值p，Wp表示该峰值包含的量化范围，np表示直方图中峰值的数目，φp表示Wp中最大直方图中的量化数值，HD表示直方图，通过计算梯度幅值|ΔG|大于阈值tho且角度θ在某个范围的点的个数得到，HD(k)表示点的幅值|ΔG|≥tho且角度θ满足(2k-1)π/2n≤θ≤...

【专利技术属性】
技术研发人员：路小波，陶焕杰，
申请(专利权)人：东南大学，
类型：发明
国别省市：江苏,32