CN107292905A - 基于改进混合高斯模型的运动目标检测方法 - Google Patents
基于改进混合高斯模型的运动目标检测方法 Download PDFInfo
- Publication number
- CN107292905A CN107292905A CN201710379773.9A CN201710379773A CN107292905A CN 107292905 A CN107292905 A CN 107292905A CN 201710379773 A CN201710379773 A CN 201710379773A CN 107292905 A CN107292905 A CN 107292905A
- Authority
- CN
- China
- Prior art keywords
- gauss model
- pixel
- gauss
- model
- variance
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于改进混合高斯模型的运动目标检测方法,主要解决现有技术计算复杂度高和检测准确度低的问题。其实现步骤是:首先,采用高斯函数对输入视频集进行预处理,去除图像中的噪声;其次,使用图像的空间信息初始化高斯混合模型,自适应确定高斯模型的个数;然后,采用lambert光照模型检测全局光照变化,并对模型的权值进行分层次更新;最后,根据模型的权值和彩色图像与高斯混合模型的距离检测出运动目标图像,再对该目标图像进行后处理。本发明收敛速度快,计算复杂度低,对复杂场景适应性强,检测准确度高,可用于智能监控、智能交通、人机交互、视觉导航、虚拟现实、医学诊断和可视化手术。
Description
技术领域
本发明属于计算机视觉领域,具体地说是一种运动目标检测方法,可用于智能监控、智能交通、人机交互、视觉导航、虚拟现实、医学诊断与可视化手术。
背景技术
随着社会的发展和科技的进步,城市视频监控网络的发展越来越迅速,覆盖到了人类活动的各种场所。视频监控在很大程度上保障着城市的人员财产安全,极大地减少违法犯罪行为。视频监控主要关注的技术是对监控视频中的运动物体进行检测、跟踪、识别和分析等操作。利用这些技术,人们可以快速获得感兴趣目标的位置、轨迹以及行为等有效信息。运动目标检测是运动目标跟踪、行为识别和场景描述等技术的基础,检测的结果直接影响后续算法的准确性。因此,如何提高目标检测的准确性和鲁棒性,成为计算机视觉领域的主要研究方向之一。
目前,运动目标检测方法主要有:帧间差分法、背景减除法和光流法。其中,背景减除法是使用最广泛的方法,它通过比较当前帧和背景模型得到运动目标。因此,算法的性能高度取决于背景建模方法。高斯混合模型GMM因其简单性,多模性,实时性成为最受欢迎的背景减除算法。该算法可以较好地抑制动态背景干扰,但该方法对突然的光照变化和运动目标的阴影敏感,且初始化困难,收敛过程比较耗时。
在过去的20年里,大量文献致力于混合高斯模型的研究,主要包括模型初始化方法、模型更新方法和高斯模型个数估计等。在模型的初始化方法中,大部分的方法是基于Stauffer和Grimson提出的方法,该方法将每个像素的所有的高斯分量初始化为相同的分布,即均值设置为第一帧的相应像素的像素值,方差初始化为一个较大的值,初始权重为1。模型初始化方法也仅仅涉及到模型初始方差的计算,模型的初始均值和高斯模型权值仍然采用上述方法得到,使得模型需要较长时间才能达到稳定。模型更新方法只讨论了单个高斯模型的学习率,未考虑各高斯模型之间的关系。同时,大部分算法只能检测局部光照变化,而不能检测全局光照变化。在高斯模型个数估计上,一般采用分裂和合并的方式自适应调整高斯模型的个数,但这些算法收敛速度慢。
发明内容
针对现有混合高斯模型的不足,本发明提出了一种基于改进混合高斯模型的运动目标检测方法,以提高复杂场景中模型的收敛速度,减少计算复杂度,提高运动目标的检测精度。
实现本发明的关键技术是:使用图像的空间信息来初始化混合高斯模型,以加快模型的收敛速度,自适应确定高斯模型的个数;对模型的权值进行分层次更新,以减缓暂时静止的物体融入背景;采用lambert光照模型检测全局光照变化,使得背景模型能快速从光照变化中恢复出来。其实现步骤包括如下:
(1)读入前两帧彩色图像I1和I2,利用高斯函数分别对两幅图像的R、G、B三个通道进行平滑,得到平滑后的彩色图像为L1和L2;
(2)初始化的混合高斯模型:
(2a)初始化方差:根据平滑后彩色图像L1和图像L2差分的绝对值矩阵得到初始模型的全局方差,将像素点(x,y)与其邻域像素绝对差分的最大值作为初始模型局部方差,对初始全局方差和初始局部方差进行加权平均,得到像素点(x,y)处各高斯模型的初始方差σ1(x,y);
(2b)初始化均值和高斯模型数目:根据像素(x,y)的邻域像素的分布情况确定像素(x,y)处高斯模型数目N1(x,y)和各高斯模型的均值μk,1(x,y);
(2c)分层次初始化像素(x,y)处高斯模型的权值ωk,1(x,y),将第一个高斯模型的权值ω1,1设置为w1,1=(1/N1(x,y))+0.1,其余高斯模型的权值设置为(1-w1,1)/(N1(x,y)-1);
(2d)根据高斯模型的初始均值、方差、权重和高斯模型数目,初始化高斯混合模型;
(3)读入第t帧彩色图像It,利用高斯函数对It进行平滑,得到平滑后的第t帧彩色图像Lt;
(4)模型更新:
(4a)根据Lambert光照模型,得到第t帧彩色图像Lt的全局光照变化因子,根据全局光照变化因子改变第t-1帧构建的高斯混合模型的均值μk,t-1(x,y)和方差σk,t-1(x,y);
(4b)根据第t帧彩色图像Lt与高斯混合模型的匹配情况,更新高斯模型的数目Nt(x,y)、均值μk,t(x,y)和方差σk,t(x,y),分层次更新权重ωk,t(x,y),得到更新后的混合高斯模型;
(5)背景估计和前景检测:在像素(x,y)处,若存在高斯模型的权值ωk,t(x,y)大于门限Thω,且|Lt(x,y)-μk,t(x,y)|≤Dσk,t(x,y),则像素(x,y)为背景,即fgt(x,y)=0;否则,像素(x,y)为运动目标,即fgt(x,y)=1,得到第t帧彩色图像的前景图像fgt;
(6)利用奇异值分解、中值滤波和形态学方法对上述得到的前景图像fgt进行后处理,输出处理后的前景图像fgt,即运动目标图像;
(7)重复(3-6),直到视频结束。
本发明具有如下优点:
本发明由于引入空间信息对混合高斯模型进行初始化,加快了模型的收敛速度,减少了高斯模型的数目,提高了运动目标检测的精度;同时由于本发明在混合高斯模型的更新阶段,采用分层次更新的方法更新高斯模型的权值,延缓了暂时静止的物体融入背景的速度;此外由于本发明对全局和局部光照变化情况进行判断,根据光照变化情况,改变高斯模型的均值和方差,有效降低了运动目标误检率。
附图说明
图1本发明的实现流程图;
图2用本发明与Stauffer和Grimson提出的方法在“Waving Trees”视频序列的前三帧彩色图像上进行运动目标检测的结果图;
图3为用本发明在“Waving Trees”视频序列上得到的混合高斯模型的初始高斯模型数目结果图;
图4为用本发明与GMM|RECTGAUSS-Tex算法对“Wallflower dataset”视频集的运动目标检测结果图。
具体实施方式
以下结合附图对本发明的技术方案和效果进行进一步说明:
参照图1,本发明的实现步骤如下:
步骤1,对输入的彩色图像进行平滑处理。
现有的图像平滑方法有:均值滤波,高斯滤波,中值滤波;本发明使用高斯滤波对输入的彩色图像进行平滑处理,具体步骤如下:
读入前两帧彩色图像I1和I2,利用高斯函数分别对两幅图像的R、G、B三个通道进行平滑,得到平滑后的彩色图像为L1和L2。
所述二维高斯函数G为:
其中,σ为标准差,取值为0.5;
用高斯函数G与图像It进行卷积,得到平滑后的彩色图像Lt:
Lt=G*It,t=1,2。
其中,符号*为卷积操作。
步骤2,初始化混合高斯模型。
2.1)初始化方差:
现有的初始化方差的方法有:Stauffer和Grimson提出的方法,Evangelio和Patzold提出的方法。其中Stauffer和Grimson提出的方法,是将每个像素的所有的高斯模型初始化为相同的分布,方差初始化为一个较大的值;Evangelio和Patzold提出的方法,是根据第一帧和第二帧的绝对差分的中值作为初始方差;本发明根据平滑后彩色图像L1和图像L2差分的绝对值矩阵得到初始模型的全局方差,将像素点(x,y)与其邻域像素绝对差分的最大值作为初始模型局部方差,对初始全局方差和初始局部方差进行加权平均,得到像素点(x,y)处各高斯模型的初始方差σ1(x,y),具体步骤如下:
2.1.1)计算高斯模型的全局方差
(2.1.1.1)计算平滑后第一帧彩色图像L1和平滑后第二帧彩色图像L2的差分的绝对值矩阵diff:
diff=|L2-L1|;
(2.1.1.2)将绝对值矩阵diff分别在R、G、B三个通道上等间隔量化为P个等级,得到量化后的彩色图像Q:
其中,P=64,表示向上取整;
(2.1.1.3)分别计算彩色图像Q在R、G、B三个通道上的统计直方图分量:
其中,c={R,G,B},Qc为彩色图像Q的R、G、B三个通道之一的图像,Hc(i)为Qc在第i个量化等级上的直方图分量,P为最大量化等级,P个直方图分量构成Qc的直方图Hc,m,n表示彩色图像L1的行数和列数,δ(·)表示单位脉冲函数;
(2.1.1.4)对直方图Hc的前a个直方图分量Hc(i)求和,直到得到的和超过阈值T,得到背景变化均值ac:
其中,T=0.5;
(2.1.1.5)通过背景变化均值ac得到高斯模型在通道c上的全局方差
则初始高斯模型的全局方差为:
2.1.2)计算高斯模型的局部方差:
从平滑后第一帧彩色图像L1中选取像素(x,y)的n1×n1邻域像素(p,q),n1=5;将L1(p,q)与L1(x,y)的绝对差分的最大值作为局部方差在通道c上的局部方差为:
其中,表示像素点(p,q)在通道c上的像素值,表示像素点(x,y)在通道c上的像素值;
则初始高斯模型的局部方差为:
2.1.3)计算初始方差:
对全局方差和局部方差进行加权平均,得到像素点(x,y)处高斯模型的初始方差σ1(x,y):
其中,β=0.2表示加权系数,为了防止模型过拟合或欠拟合,设置方差的上限为40和下限为3,即初始方差值为:
2.2)初始化均值和高斯模型数目:
现有的初始化均值和高斯模型数目的方法有:Stauffer和Grimson提出的方法,该方法在每个像素点都使用相同的高斯模型数目,且每个像素处的所有高斯模型都将第一帧彩色图像在该像素点的像素值作为高斯模型的初始均值;本发明根据每个像素点的邻域像素的强度分布情况确定每个像素的初始高斯模型数目和各高斯模型的初始均值,具体步骤如下:
(2.2.1)初始化像素(x,y)处的高斯模型的个数N1(x,y)=1,将像素(x,y)对应彩色图像L1(x,y)的像素值作为第一个高斯模型的均值μ1,1(x,y),将该高斯模型的匹配次数η1,1(x,y)初始化为1;
(2.2.2)判断像素(x,y)的n1×n1邻域像素中是否存在一个像素点(p,q),且该像素点的像素值L1(p,q)与已有高斯模型均值的绝对差是否大于3倍的初始方差σ1(x,y):如果是,则在像素(x,y)处增加一个高斯模型,将像素点(p,q)的像素值L1(p,q)作为所增加的高斯模型的均值,该高斯模型的匹配次数初始化为1;否则,高斯模型的数目不变;
(2.2.3)判断高斯模型数目N1(x,y)是否小于最大高斯模型数目K=5,且像素(x,y)的n1×n1邻域像素中是否存在与已有高斯模型均值的绝对差大于3倍的初始方差σ1(x,y)的像素点,如果是,重复步骤(2.2.2);否则,循环结束,得到像素(x,y)处高斯混合模型的初始化均值μk,1(x,y)和初始化高斯模型数目N1(x,y);
2.3)初始化权值:
现有的方法有:Stauffer和Grimson提出的方法,该方法将每一个像素点处的所有高斯模型的权值初始化为相同的值;本发明采用分层次初始化高斯模型的权值ωk,1(x,y),即将第一个高斯模型的权值设置为一个较大的值,其余高斯模型的权值设置为一个较小的权值:
其中,w=(1/N1(x,y))+0.1为第一个高斯模型的权值;
2.4)根据高斯模型的初始均值、方差、权重和高斯模型数目,初始化高斯混合模型;在像素点(x,y)处的初始化高斯混合模型P(x,y)为:
其中,ωk,1(x,y)、μk,1(x,y)、分别为第k个高斯模型在1时刻的权值,均值和协方差矩阵,是高斯模型的概率密度函数。
步骤3读入第t帧彩色图像It,t>1,利用高斯函数对It进行平滑,得到平滑后的第t帧彩色图像Lt;
步骤4模型更新
现有的模型更新方法有:Stauffer和Grimson提出的方法,Kaewtrakulpong和Bowden提出的方法,Zivkovic提出的方法和Lee提出的方法。其中:
Stauffer和Grimson提出的方法选择混合高斯模型中与图像匹配的高斯模型,对该高斯模型的均值和方差进行更新,对应权值增加,其余高斯模型均值和方差保持不变,对应的权值减小;Kaewtrakulpong和Bowden提出的方法在不同的阶段使用不同的更新方程以提高模型学习速度和精度;Zivkovic提出的方法使用先验信息动态估计高斯模型的数目;Lee提出的方法根据高斯模型的匹配次数自适应确定学习率;本发明使用Stauffer和Grimson提出的方法对高斯模型的均值和方差进行更新,高斯模型权值采用分层次更新的方法,根据彩色图像与高斯混合模型的匹配情况,更新高斯模型的数目,其步骤如下:
4.1)检测全局光照变化:
(4.1.1)根据Lambert光照模型,计算平滑后第t帧彩色图像Lt和平滑后第t-1帧彩色图像Lt-1对应像素点像素值的比值,得到比值矩阵R,选取比值矩阵R中所有元素的中值作为全局光照变化因子ratio;
(4.1.2)判断|ratio-1|的值是否大于设定门限0.2,如果是,表示出现了全局光照变化,对第t-1帧构建的高斯混合模型的均值μk,t-1(x,y)和方差σk,t-1(x,y)按照全局光照变化因子进行相应的变化;否则,保持第t-1帧的高斯混合模型不变,即:
其中,k为在像素(x,y)的第k个高斯模型,k=1,…,Nt-1(x,y),Nt-1(x,y)为第t-1帧在像素点(x,y)处构建的高斯混合模型的高斯模型数目,ηk,t-1(x,y)表示第t-1帧在像素点(x,y)处构建的高斯混合模型的第k个高斯模型的匹配次数。
4.2)将平滑后的第t帧彩色图像Lt在像素点(x,y)处的像素值Lt(x,y)与所对应的Nt-1(x,y)个高斯模型的均值进行比较,判断在这Nt-1(x,y)个高斯模型中是否存在满足条件:|Lt(x,y)-μm,t-1(x,y)|<D·σm,t-1(x,y)的高斯模型m,其中,1≤m≤Nt-1(x,y),D=3表示常系数:
如果是,则第m个高斯模型与像素值Lt(x,y)匹配,更新第m个高斯模型的匹配次数ηm,t(x,y),均值μm,t(x,y),方差σm,t和权值ωm,t:
ηm,t(x,y)=ηm,t-1(x,y)+1
μm,t(x,y)=(1-ρm,t(x,y))μm,t-1(x,y)+ρm,tLt(x,y)
ωm,t=(1-α)ωm,t-1+α
其中,α=0.01表示权值的学习率,表示均值和方差的学习率;
保持高斯模型的数目不变,其余高斯模型k的均值,保持方差和匹配次数不变,将权值ωk,t更新为:
其中,Thω=0.2表示设定阈值;
否则,新建第i个高斯模型,更新高斯模型的数目Nt(x,y):
初始化第i个高斯模型均值μi,t(x,y)、方差σi,t(x,y)、权值ωi,t(x,y)和匹配次数ηi,t(x,y):
μi,t(x,y)=Lt(x,y)
σi,t(x,y)=σ1(x,y)
ηi,t(x,y)=1
其中,初始方差σ1(x,y)由平滑后彩色图像Lt和平滑后彩色图像Lt-1按照步骤2.1)的方法计算得到。
4.3)归一化高斯模型的权值ωk,t(x,y):
步骤5背景估计和前景检测。
在像素(x,y)处,判断是否若存在高斯模型的权值ωk,t(x,y)大于设定门限Thω=0.2,且是否满足|Lt(x,y)-μk,t(x,y)|≤3σk,t(x,y),如果是,则像素点(x,y)为背景,即fgt(x,y)=0;否则,像素点(x,y)为运动目标,即fgt(x,y)=1;得到第t帧彩色图像的前景图像fgt。
步骤6对上述得到的前景图像fgt进行后处理。
现有的后处理方法有:高斯滤波,中值滤波,均值滤波,形态学方法;本发明使用利用现有奇异值分解的方法检测局部光照变化,去除由于局部光照引起的误检,利用中值滤波去除噪声,利用形态学方法的闭运算消除目标内部较小的“空洞”,具体步骤如下:
6.1)检测局部光照变化
首先,使用平滑后的第t帧彩色图像Lt的像素点(x,y)的5×5邻域像素值构成邻域矩阵At(x,y),对邻域矩阵At(x,y)进行奇异值分解:
At(x,y)=UtStVt T
其中,Ut和Vt表示矩阵At(x,y)的特征向量矩阵,St=diag(λt,1,λt,2,λt,3,λt,4,λt,5)为At(x,y)的特征值矩阵,λt,j表示At(x,y)的奇异值,1≤j≤5,且λt,1≥λt,2≥λt,3≥λt,4≥λt,5;
然后,将奇异值λt,2,λtF,3,λt,4,λt,5除以λt,1,并求和,得到光照不变性特征ft(x,y):
最后,判断第t帧光照不变性特征ft(x,y)与和第t-1帧光照不变性特征ft-1(x,y)绝对差是否小于阈值Thsvd=0.05:如果是,则像素点(x,y)处发生了局部光照变化,像素点(x,y)应为背景,否则,像素点(x,y)处没有发生了局部光照变化;像素点(x,y)处前景图像fgt(x,y)为:
6.2)使用5×5的中值滤波去除在步骤6.1)所得到的前景图像fgt中的噪声,得到去除噪声后的前景图像fgt;
6.3)使用形态学的闭运算,填补轮廓线中的裂痕,消除目标内部较小的“空洞”,即先用模板矩阵B对去除噪声后的前景图像fgt进行膨胀运算;然后用模板矩阵B对膨胀后的结果进行腐蚀运算,得到预处理后的前景图像fgt,其中模板矩阵B为5×5单位矩阵;
6.4)输出处理后的前景图像fgt,即运动目标图像。
步骤7重复步骤3-步骤6,直到视频结束。
本发明的效果可通过以下实验进一步说明:
1.实验条件
实验环境:Intel Core i7 CPU 2.00Ghz,8GB内存,Matlab2014b实验平台。
2.实验内容与结果
实验1:测试本发明的混合高斯模型初始化方法
采用本发明与Stauffer和Grimson提出的初始化方法对“Wallflower dataset”中“Waving Trees”视频序列的第1-3帧进行测试,结果如图2所示,其中:
图2第一行从左到右分别为原始图像的第1-3帧;
图2第二行为采用Stauffer和Grimson提出的方法得到第1-3帧前景图像fg;
图2第三行为采用本发明得到第1-3帧前景图像fg;
从图2中可以看出,本发明的混合高斯模型初始化方法能够在初始时刻建立起稳定的背景模型。同时,该初始化方法对背景干扰不敏感,能有效提高检测精度。
实验2:测试本发明的混合高斯模型的初始化高斯模型数目
采用本发明所提出的初始化高斯模型数目的方法对“Wallflower dataset”中“Waving Trees”视频序列进行测试,结果如图3所示:
图3中颜色越深代表该像素的高斯模型个数越少,如:蓝色代表1个高斯分布,青色代表2个高斯分布,黄色代表3个高斯分布。
从图3中可以看出,在大部分图像区域,只需要1个或2个高斯模型,只有极小的区域需要3个高斯模型;而现有的混合高斯模型会在每一个像素点设定3到5个高斯模型;表明本发明能够有效降低高斯模型的数目,减少计算复杂度。
实验3:测试本发明的检测精度
采用本发明,GMM|Stauffer&Grimson方法和GMM|RECTGAUSS-Tex方法对“Wallflower dataset”视频集进行测试,三种方法对视频集中各视频的正检率(TPR)如表1所示,三种方法对视频集中各视频的误检率(FPR)如表2所示。
表1三种算法的TPR
表2三种算法的FPR
由表1可知,在六种场景中本发明都具有相对较高的TPR值。
由表2可知,除了在Camouflage和Foreground Aperture场景中本发明的FPR略高于其他两种算法外,其余场景本发明的FPR都较低。尤其是发生全局光照变化时,本发明具有较低的误检率。
实验4:测试本发明对运动目标的检测效果
用本发明与GMM|RECTGAUSS-Tex算法对“Wallflower dataset”视频集进行运动目标检测结果图,结果如图4所示,其中:
图4第一列为六种场景的原始图像;
图4第二列为手动分割的前景图像;
图4第三列为本发明得到的前景图像;
图4第4列为GMM|RECTGAUSS-Tex算法得到的前景图像;
从图4可以看出,本发明不仅对阴影、光照变化不敏感,而且能延缓暂时静止的物体融入背景,能够得到更好的检测结果。
综上所述,本发明能够有效检测复杂场景中的运动目标。
Claims (5)
1.一种基于改进混合高斯模型的运动目标检测方法,包括:
(1)读入前两帧彩色图像I1和I2,利用高斯函数分别对两幅图像的R、G、B三个通道进行平滑,得到平滑后的彩色图像为L1和L2;
(2)初始化的混合高斯模型:
(2a)初始化方差:根据平滑后彩色图像L1和图像L2差分的绝对值矩阵得到初始模型的全局方差,将像素点(x,y)与其邻域像素绝对差分的最大值作为初始模型局部方差,对初始全局方差和初始局部方差进行加权平均,得到像素点(x,y)处各高斯模型的初始方差σ1(x,y);
(2b)初始化均值和高斯模型数目:根据像素(x,y)的邻域像素的分布情况确定像素(x,y)处高斯模型数目N1(x,y)和各高斯模型的均值μk,1(x,y);
(2c)分层次初始化像素(x,y)处高斯模型的权值ωk,1(x,y),将第一个高斯模型的权值ω1,1设置为w1,1=(1/N1(x,y))+0.1,其余高斯模型的权值设置为(1-w1,1)/(N1(x,y)-1);
(2d)根据高斯模型的初始均值、方差、权重和高斯模型数目,初始化高斯混合模型;
(3)读入第t帧彩色图像It,利用高斯函数对It进行平滑,得到平滑后的第t帧彩色图像Lt;
(4)模型更新:
(4a)根据Lambert光照模型,得到第t帧彩色图像Lt的全局光照变化因子,根据全局光照变化因子改变第t-1帧构建的高斯混合模型的均值μk,t-1(x,y)和方差σk,t-1(x,y);
(4b)根据第t帧彩色图像Lt与高斯混合模型的匹配情况,更新高斯模型的数目Nt(x,y),均值μk,t(x,y)和方差σk,t(x,y),分层次更新权重ωk,t(x,y),得到更新后的混合高斯模型;
(5)背景估计和前景检测:在像素(x,y)处,若存在高斯模型的权值ωk,t(x,y)大于门限Thω,且|Lt(x,y)-μk,t(x,y)|≤Dσk,t(x,y),则像素(x,y)为背景,即fgt(x,y)=0;否则,像素(x,y)为运动目标,即fgt(x,y)=1,得到第t帧彩色图像的前景图像fgt;
(6)利用奇异值分解、中值滤波和形态学方法对上述得到的前景图像fgt进行后处理,输出处理后的前景图像fgt,即运动目标图像;
(7)重复(3-6),直到视频结束。
2.根据权利要求1所述的方法,其中步骤(2a)中根据平滑后彩色图像L1和图像L2差分的绝对值矩阵得到初始模型的全局方差,按如下步骤进行:
(2a1)计算第一帧彩色图像L1和第二帧彩色图像L2的差分的绝对值矩阵diff:
diff=|L2-L1|;
(2a2)将绝对值矩阵diff分别在R、G、B三个通道上等间隔量化为P个等级,得到量化后的彩色图像Q:
其中,P=64,表示向上取整;
(2a3)分别计算彩色图像Q在R、G、B三个通道上的统计直方图分量:
其中,c={R,G,B},Qc分别为彩色图像Q的R、G、B三个通道的量化图像,Hc(i)为Qc在第i个量化等级上的直方图分量,P为最大量化等级,P个直方图分量构成Qc的直方图Hc,m,n表示彩色图像L1的行数和列数,δ(·)表示单位脉冲函数;
(2a4)对直方图Hc的前a个直方图分量Hc(i)求和,使其和值超过设定的阈值T,得到背景变化均值ac:
其中,T=0.5;
(2a5)根据背景变化均值ac,得到高斯模型在通道c上的全局方差
则初始高斯模型的全局方差为:
3.根据权利要求1所述的方法,其中步骤(2b)中根据像素(x,y)的邻域像素的分布情况确定高斯模型数目N1(x,y)和各高斯模型的均值μk,1(x,y),按如下步骤进行:
(2b1)初始化像素(x,y)处的高斯模型的个数N1(x,y)=1,将像素(x,y)对应彩色图像L1(x,y)的像素值作为第一个高斯模型的均值μ1,1(x,y),该高斯模型的匹配次数η1,1(x,y)初始化为1;
(2b2)判断像素(x,y)的n1×n1邻域像素中是否存在一个像素点(p,q),且该像素点的像素值L1(p,q)与已有高斯模型均值的绝对差是否大于3倍初始方差σ1(x,y),如果是,则在像素(x,y)处增加一个高斯模型,即对高斯模型的数目N1(x,y)加1,将像素点(p,q)的像素值L1(p,q)作为所增加的高斯模型的均值,该高斯模型的匹配次数初始化为1;否则,高斯模型的数目不变;
(2b3)判断高斯模型数目N1(x,y)是否小于最大高斯模型数目K=5,且像素(x,y)的n1×n1邻域像素中是否存在与已有高斯模型均值的绝对差大于3倍的初始方差σ1(x,y)的像素点,如果是,重复步骤(2.2.2);否则,循环结束,得到像素(x,y)处高斯混合模型的初始化均值μk,1(x,y)和初始化高斯模型数目N1(x,y)。
4.根据权利要求1所述的方法,其中步骤(4a)中根据全局光照变化因子改变在第t-1帧构建的高斯混合模型的均值μk,t-1(x,y)和方差σk,t-1(x,y),按如下步骤进行:
(4a1)根据Lambert光照模型,计算第t帧彩色图像Lt和第t-1帧彩色图像Lt-1对应像素点像素值的比值,得到比值矩阵R,选取比值矩阵R中所有元素的中值作为全局光照变化因子ratio;
(4a2)判断|ratio-1|的值是否大于设定门限0.2,如果是,表示出现了全局光照变化,对第t-1帧构建的高斯混合模型的均值μk,t-1(x,y)和方差σk,t-1(x,y)按照全局光照变化因子进行相应的变化;否则,第t-1帧构建的高斯混合模型保持不变;即:
其中,k为在像素(x,y)的第k个高斯模型,k=1,…,Nt-1(x,y),Nt-1(x,y)为第t-1帧在像素点(x,y)处构建的高斯混合模型的高斯模型数目,ηk,t-1(x,y)表示第t-1帧在像素点(x,y)处构建的高斯混合模型的第k个高斯模型的匹配次数。
5.根据权利要求1所述的方法,其中步骤(4b)中根据第t帧彩色图像Lt与高斯混合模型的匹配情况,更新高斯模型的数目Nt(x,y),均值μk,t(x,y)和方差σk,t(x,y),分层次更新权重ωk,t(x,y),按如下步骤进行:
(4b1)将平滑后的第t帧彩色图像Lt在像素点(x,y)处的像素值Lt(x,y)与所对应的Nt-1(x,y)个高斯模型的均值进行比较,判断在这Nt-1(x,y)个高斯模型中是否存在满足条件:|Lt(x,y)-μm,t-1(x,y)|<D·σm,t-1(x,y)的高斯模型m,其中,1≤m≤Nt-1(x,y),D=3表示常系数:
如果是,则第m个高斯模型与像素值Lt(x,y)匹配,更新第m个高斯模型的匹配次数ηm,t(x,y),均值μm,t(x,y),方差σm,t和权值ωm,t:
ηm,t(x,y)=ηm,t-1(x,y)+1
μm,t(x,y)=(1-ρm,t(x,y))μm,t-1(x,y)+ρm,tLt(x,y)
ωm,t=(1-α)ωm,t-1+α
其中,α=0.01表示权值的学习率,表示均值和方差的学习率;
保持高斯模型的数目不变,其余高斯模型k的均值,保持方差和匹配次数不变,将权值ωk,t更新为:
其中,Thω=0.2表示设定阈值;
否则,新建第i个高斯模型,更新高斯模型的数目Nt(x,y):
初始化第i个高斯模型均值μi,t(x,y)、方差σi,t(x,y)、权值ωi,t(x,y)和匹配次数ηi,t(x,y):
μi,t(x,y)=Lt(x,y)
σi,t(x,y)=σ1(x,y)
ηi,t(x,y)=1
其中,初始方差σ1(x,y)由平滑后彩色图像Lt和平滑后彩色图像Lt-1按照步骤2.1)的方法计算得到。
(4b2)归一化高斯模型的权值ωk,t(x,y):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710379773.9A CN107292905B (zh) | 2017-05-25 | 2017-05-25 | 基于改进混合高斯模型的运动目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710379773.9A CN107292905B (zh) | 2017-05-25 | 2017-05-25 | 基于改进混合高斯模型的运动目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292905A true CN107292905A (zh) | 2017-10-24 |
CN107292905B CN107292905B (zh) | 2020-08-04 |
Family
ID=60094038
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710379773.9A Active CN107292905B (zh) | 2017-05-25 | 2017-05-25 | 基于改进混合高斯模型的运动目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292905B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110619645A (zh) * | 2019-09-25 | 2019-12-27 | 上海海瞩智能科技有限公司 | 一种桥吊下集装箱拖架的自动识别与定位装置及方法 |
CN110765979A (zh) * | 2019-11-05 | 2020-02-07 | 中国计量大学 | 一种基于背景建模与光控的智能led庭院灯 |
CN110991398A (zh) * | 2019-12-18 | 2020-04-10 | 长沙融创智胜电子科技有限公司 | 一种基于改进步态能量图的步态识别方法及系统 |
CN111967469A (zh) * | 2020-08-13 | 2020-11-20 | 上海明略人工智能(集团)有限公司 | 一种畸形文本矫正方法、系统及文字识别方法 |
CN112101148A (zh) * | 2020-08-28 | 2020-12-18 | 普联国际有限公司 | 一种运动目标检测方法、装置、存储介质及终端设备 |
CN113570640A (zh) * | 2021-09-26 | 2021-10-29 | 南京智谱科技有限公司 | 一种视频图像处理的方法及装置 |
CN113705817A (zh) * | 2021-08-10 | 2021-11-26 | 石家庄学院 | 基于高阶高斯混合模型的远程实时监控数据处理方法 |
CN114627079A (zh) * | 2022-03-16 | 2022-06-14 | 江苏立晶工业科技有限公司 | 基于人工智能的安防监控设备异常检测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090161981A1 (en) * | 2007-12-20 | 2009-06-25 | United States Of America As Represented By The Secretary Of The Navy | Method for Enhancing Ground-Based Detection of a Moving Object |
US20100142809A1 (en) * | 2008-12-08 | 2010-06-10 | Electronics And Telecommunications Research Institute | Method for detecting multi moving objects in high resolution image sequences and system thereof |
CN104467742A (zh) * | 2014-12-16 | 2015-03-25 | 中国人民解放军海军航空工程学院 | 基于高斯混合模型的传感器网络分布式一致性粒子滤波器 |
CN106504273A (zh) * | 2016-10-28 | 2017-03-15 | 天津大学 | 一种基于gmm运动目标检测的改进算法 |
-
2017
- 2017-05-25 CN CN201710379773.9A patent/CN107292905B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090161981A1 (en) * | 2007-12-20 | 2009-06-25 | United States Of America As Represented By The Secretary Of The Navy | Method for Enhancing Ground-Based Detection of a Moving Object |
US20100142809A1 (en) * | 2008-12-08 | 2010-06-10 | Electronics And Telecommunications Research Institute | Method for detecting multi moving objects in high resolution image sequences and system thereof |
CN104467742A (zh) * | 2014-12-16 | 2015-03-25 | 中国人民解放军海军航空工程学院 | 基于高斯混合模型的传感器网络分布式一致性粒子滤波器 |
CN106504273A (zh) * | 2016-10-28 | 2017-03-15 | 天津大学 | 一种基于gmm运动目标检测的改进算法 |
Non-Patent Citations (4)
Title |
---|
MARIN S. HALPER: "Global, Local, and Stochastic Background Modeling for Target Detection in Mixed Pixels", 《PROCEEDINGS OF SPIE》 * |
YONGQUAN ZHANG等: "A novel fast partitioning algorithm for extended target tracking using a Gaussian mixture PHD filter", 《SIGNAL PROCESSING》 * |
肖武: "基于高斯混合模型的运动目标检测算法研究和应用", 《万方数据知识服务平台》 * |
胥欣等: "混合高斯模型运动检测算法优化", 《计算机应用研究》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110619645B (zh) * | 2019-09-25 | 2022-11-25 | 上海海瞩智能科技有限公司 | 一种桥吊下集装箱拖架的自动识别与定位装置及方法 |
CN110619645A (zh) * | 2019-09-25 | 2019-12-27 | 上海海瞩智能科技有限公司 | 一种桥吊下集装箱拖架的自动识别与定位装置及方法 |
CN110765979A (zh) * | 2019-11-05 | 2020-02-07 | 中国计量大学 | 一种基于背景建模与光控的智能led庭院灯 |
CN110991398A (zh) * | 2019-12-18 | 2020-04-10 | 长沙融创智胜电子科技有限公司 | 一种基于改进步态能量图的步态识别方法及系统 |
CN111967469A (zh) * | 2020-08-13 | 2020-11-20 | 上海明略人工智能(集团)有限公司 | 一种畸形文本矫正方法、系统及文字识别方法 |
CN111967469B (zh) * | 2020-08-13 | 2023-12-15 | 上海明略人工智能(集团)有限公司 | 一种畸形文本矫正方法、系统及文字识别方法 |
CN112101148A (zh) * | 2020-08-28 | 2020-12-18 | 普联国际有限公司 | 一种运动目标检测方法、装置、存储介质及终端设备 |
CN112101148B (zh) * | 2020-08-28 | 2024-05-03 | 普联国际有限公司 | 一种运动目标检测方法、装置、存储介质及终端设备 |
CN113705817A (zh) * | 2021-08-10 | 2021-11-26 | 石家庄学院 | 基于高阶高斯混合模型的远程实时监控数据处理方法 |
CN113705817B (zh) * | 2021-08-10 | 2023-07-28 | 石家庄学院 | 基于高阶高斯混合模型的远程实时监控数据处理方法 |
CN113570640A (zh) * | 2021-09-26 | 2021-10-29 | 南京智谱科技有限公司 | 一种视频图像处理的方法及装置 |
CN113570640B (zh) * | 2021-09-26 | 2021-12-10 | 南京智谱科技有限公司 | 一种视频图像处理的方法及装置 |
CN114627079A (zh) * | 2022-03-16 | 2022-06-14 | 江苏立晶工业科技有限公司 | 基于人工智能的安防监控设备异常检测方法及系统 |
CN114627079B (zh) * | 2022-03-16 | 2023-03-10 | 江苏立晶工业科技有限公司 | 基于人工智能的安防监控设备异常检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107292905B (zh) | 2020-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107292905A (zh) | 基于改进混合高斯模型的运动目标检测方法 | |
CN106778595B (zh) | 基于高斯混合模型的人群中异常行为的检测方法 | |
Migdal et al. | Background subtraction using markov thresholds | |
Zhou et al. | Real time robust human detection and tracking system | |
Porikli et al. | Bayesian background modeling for foreground detection | |
CN111062273B (zh) | 一种遗留物品追溯检测与报警方法 | |
Cheng et al. | Scene analysis for object detection in advanced surveillance systems using Laplacian distribution model | |
CN110765964A (zh) | 基于计算机视觉的电梯轿厢内异常行为的检测方法 | |
CN103530893A (zh) | 相机抖动场景下基于背景差和运动信息的前景检测方法 | |
CN106204594A (zh) | 一种基于视频图像的弥散性运动物体的方向检测方法 | |
CN109919053A (zh) | 一种基于监控视频的深度学习车辆停车检测方法 | |
CN106650617A (zh) | 一种基于概率潜在语义分析的行人异常识别方法 | |
Wu et al. | Properties of orthogonal Gaussian-Hermite moments and their applications | |
CN106056078A (zh) | 一种基于多特征回归式集成学习的人群密度估计方法 | |
CN103617637B (zh) | 一种基于字典学习的低照度运动检测方法 | |
Wang et al. | Smoke Recognition based on Dictionary and BP Neural Network. | |
Li et al. | Combination of local binary pattern operator with sample consensus model for moving objects detection | |
CN112613430B (zh) | 一种基于深度迁移学习的步态识别方法 | |
Frejlichowski et al. | Extraction of the foreground regions by means of the adaptive background modelling based on various colour components for a visual surveillance system | |
Subudhi et al. | Moving object detection using Gaussian background model and Wronskian framework | |
Chong et al. | Visual based fall detection with reduced complexity horprasert segmentation using superpixel | |
Wen et al. | People tracking and counting for applications in video surveillance system | |
Lepisk | The use of optic flow within background subtraction | |
Benezeth et al. | Comparative study on foreground detection algorithms for human detection | |
Zou et al. | Occupancy measurement by object tracking at building entrances |
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 |