CN102646280A - 一种基于块累加统计的背景实时更新方法 - Google Patents
一种基于块累加统计的背景实时更新方法 Download PDFInfo
- Publication number
- CN102646280A CN102646280A CN201210099687XA CN201210099687A CN102646280A CN 102646280 A CN102646280 A CN 102646280A CN 201210099687X A CN201210099687X A CN 201210099687XA CN 201210099687 A CN201210099687 A CN 201210099687A CN 102646280 A CN102646280 A CN 102646280A
- Authority
- CN
- China
- Prior art keywords
- piece
- pixel
- background
- buf
- field picture
- 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
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于块累加统计的背景实时更新方法。该方法首先将初始背景和视频图像进行分割,并将背景每个块的Bi[N]、Fi、Ci和BUFi[N]初始化为0;接着进行背景的块累加技术;最后根据背景的块累加计数器的数值对背景进行实时更新。本发明的方法能及时更新场景中的光线变化,且计算简单、更新速度快、噪声点小,在道路畅通和车流量很大时都能更新出较完整且正确率高的背景,所更新的背景能满足实时运动目标的提取需要,具有较强的适用性。
Description
技术领域
本发明属于视频检测技术领域,具体涉及一种基于块累加统计的背景实时更新方法。
背景技术
在视频车辆检测技术中通常采用背景差分法检测车辆等运动目标,而背景差分法的关键技术是获取实时更新的背景。目前背景更新方法主要有多帧平均法、统计直方图法、像素中值训练法和IR滤波更新法。其中,多帧平均法获取的背景粗略,不适于光线的突变,并且计算量较大;统计直方图法虽能取得较好效果,但是该方法计算过于复杂、背景更新缓慢、噪声点多;像素中值训练法需要较大的存储空间,运算速度慢,不适于实时处理;IRR滤波更新法更新的背景和真实的背景间存在误差,长时间累计可能使背景失效。
发明内容
针对现有技术的缺陷或不足,本发明的目的在于提供一种基于累加统计的背景实时更新方法,以有效地实时更新背景,使得背景差法能够准确的分离运动目标。
为了实现上述目的,本发明采用的如下技术方案:
一种基于块累加统计的背景实时更新方法,其特征在于,该方法按下述步骤进行:
步骤一,初始化背景所有像素点的像素值为0,并在块坐标系X下,将背景划分成T个块,该T个块中的第i个块内像素点(m,n)的像素值为Bi(m,n),该块内所有像素点的的像素值用二维数组Bi[N]表示,该第i个块的有效背景标志为Fi,累加计数器为Ci,该第i个块内的所有像素点的像素值缓存区域为二维数组BUFi[N],缓存区域内像素点(m,n)的像素值为BUFi(m,n),将BUFi[N]中所有像素点的像素值初始化0,并且Fi和Ci初始化值均为0;
其中:T=(W/w)*(H/h),W为背景水平方向的像素,H为背景竖直方向的像素,背景的大小为W*H,w为第i个块的宽度,h为第i个块的高度;第i个块内像素点的总个数N=w*h;i=1,2,3...T;以第i个块的左下角为原点建立直角二维坐标系Y,m代表第i个块中任一像素点在坐标系Y下的横坐标,n代表第i个块中任一像素点在坐标系Y下的纵坐标,且m=0,1,2...w-1,n=0,1,2...h-1;
步骤二,在块坐标系X下,将第j+1帧图像划分为T个块,该T个块中的第i个块内像素点(m,n)的像素值为Fi,j+1(m,n),该块内所有像素点的像素值用二维数组Fi,j+1[N]表示;
同时在块坐标系X下,将第j帧图像划分为T个块,该T个块中的第i个块内像素点(m,n)的像素值为Fi,j(m,n),该块内所有像素点的像素值用二维数组Fi,j[N]表示;
其中:j≡1,2,3,...M-1,M为图像的总帧数,且M≥71;T=(W/w)*(H/h),第j+1帧图像和第j帧图像的水平方向的像素均为W,第j+1帧图像和第j帧图像竖直方向的像素均为H,第j+1帧图像和第j帧图像的大小均为W*H,第i个块的宽度为w,第i个块的高度为h;第i个块内像素点的总个数N=w*h;以第i个块的左下角为原点建立直角二维坐标系Y,m代表第i个块中任一像素点在坐标系Y下的横坐标,n代表第i个块中任一像素点在坐标系Y下的纵坐标,且m=0,1,2...w-1,n=0,1,2...h-1;
步骤三,计算第j+1帧图像中的第i个块与第j帧图像中的第i个块之间各相同坐标的像素点处的像素值差的绝对值之和SAD,其中SAD利用公式(1)计算得到,
当SAD小于或等于阈值A时,Ci加1;
当SAD大于阈值A时,Ci=0,BUFi(m,n)=0;
重复步骤三,当i=T时,执行步骤四;
步骤四,对背景中的第i个块的累加计数器Ci进行判断:
当Ci=50时,BUFi(m,n)=Fi,j+1(m,n);
当Ci=50+20*P,且Fi=0时,首先Bi(m,n)=BUFi(m,n),Fi=1,然后BUFi(m,n)=Fi,j+1(m,n);
当Ci=50+20*P,且Fi=1时,首先Bi(m,n)=(1-0.1)*Bi(m,n)+0.1*BUFi(m,n),然后BUFi(m,n)=Fi,j+1(m,n);
其中P为正整数;
重复步骤四,当i=T时,执行步骤五;
步骤五,重复执行步骤二至步骤四,当j≥M时,即可实现背景的实时更新。
所述阈值A的取值为10×块内像素点的总个数,即A=10×(w×h)。
与现有技术相比,本发明的方法能及时更新场景中的光线变化,且计算简单、更新速度快、噪声点小,在道路畅通和车流量很大时都能更新出较完整且正确率高的背景,所更新的背景能满足实时运动目标的提取需要,具有较强的适用性。
具体实施方式
本发明的基于块累加统计的背景实时更新方法,具体按下列步骤进行:
步骤一,初始化背景所有像素点的像素值为0,并在块坐标系X下,将背景划分成T个块,该T个块中的第i个块内像素点(m,n)的像素值为Bi(m,n),该块内所有像素点的的像素值用二维数组Bi[N]表示,Bi[N]中保存有N个像素点的像素值,Bi(m,n)为Bi[N]中的一个元素,该第i个块的有效背景标志为Fi,该第i个块的累加计数器为Ci,该第i个块的所有像素点的像素值缓存区域为二维数组BUFi[N],BUFi[N]中保存N个像素点的像素值,缓存区域内像素点(m,n)的像素值为BUFi(m,n),BUFi(m,n)为BUFi[N]中的一个元素,将BUFi[N]中所有像素点的像素值初始化0,并且Fi和Ci初始化值均为0;
其中:T=(W/w)*(H/h),W为背景水平方向的像素,H为背景竖直方向的像素,背景的大小为W*H,w为第i个块的宽度,h为第i个块的高度;第i个块内像素点的总数为N=w*h;i=1,2,3...T;以第i个块的左下角为原点建立直角二维坐标系Y,m代表第i个块中任一像素点在坐标系Y下的横坐标,n代表第i个块中任一像素点在坐标系Y下的纵坐标,且m=0,1,2...w-1,n=0,1,2...h-1;
步骤二,在块坐标系X下,将第j+1帧图像(当前帧图像)划分为T个块,该T个块中的第i个块内像素点(m,n)的像素值为Fi,j+1(m,n),该块内所有像素点的像素值用二维数组Fi,j+1[N]表示,Fi,j+1[N]中保存有N个像素点的像素值,Fi,j+1(m,n)为Fi,j+1[N]中的一个元素;
同时在块坐标系X下,将第j帧图像(前一帧图像)划分为T个块,该T个块中的第i个块内像素点(m,n)的像素值为Fi,j(m,n),该块内所有像素点的像素值用二维数组Fi,j[N]表示,Fi,j[N]中保存有N个像素点的像素值,Fi,j(m,n)为Fi,j[N]中一个元素;
其中:j≡1,2,3,...M-1,M为图像的总帧数,且M≥71;T=(W/w)*(H/h),第j+1帧图像和第j帧图像的水平方向的像素均为W,第j+1帧图像和第j帧图像竖直方向的像素均为H,第j+1帧图像和第j帧图像的大小均为W*H,w为第i个块的宽度,h为第i个块的高度,第i个块内像素点的总数为N=w*h;第j帧图像中的第i个块与背景中的第i个块在块坐标系X中的位置相同;以第i个块的左下角为原点建立直角二维坐标系Y,m代表第i个块中任一像素点在坐标系Y下的横坐标,n代表第i个块中任一像素点在坐标系Y下的纵坐标,且m=0,1,2...w-1,n=0,1,2...h-1;
步骤三,在第j+1帧图像中找到第i个块,并在第j帧图像中找到位置对应的第i个块,计算这两个块中对应位置像素点的像素值差的绝对值和SAD,其中SAD利用公式(1)计算得到,
当SAD小于或等于阈值A时,说明第j+1帧图像中的第i个块与第j帧图像中的第i个块匹配,背景图像中对应的第i块的累加计数器Ci加1,其中阈值A的取值为10×块内像素点的总个数,即A=10×(w×h);
当SAD大于阈值A时,说明第j+1帧中的第i个块与第j帧相应的第i个块发生了较大的变化,这两个块不匹配,背景中的第i块的累加计数器清零,即Ci=0,并将BUFi[N]中所有像素点的像素值按式BUFi(m,n)=0逐个清零;
重复步骤三,当i=T时,执行步骤四;
步骤四,对背景中的第i个块的累加计数器Ci进行判断:
当Ci=50时,将第j+1帧图像中第i个块内所有的像素点的像素值按公式(2)逐个存储在BUFi[N]中;
BUFi(m,n)=Fi,j+1(m,n) (2)
当Ci=50+20*P,且Fi=0时,将背景中的第i个块的Bi[N]中所有像素点的像素值用该块的缓存区域BUFi[N]中保存的相应像素点(即坐标相同的像素点)的像素值按公式(3)逐个更新,背景中的第i个块的标志赋值为1即Fi=1,然后将该块的缓存区域BUFi[N]中各像素点的像素点值用第j+1帧图像中第i个块内的相应像素点的像素值按公式(2)逐个更新;
Bi(m,n)=BUFi(m,n) (3)
当第Ci=50+20*P,且Fi=1时,将背景中的第i个块的Bi[N]中各像素点的像素点值按公式(4)逐个更新,然后将该块的缓存区域BUFi[N]中各像素点的像素值用第j+1帧图像中第i个块内相应像素点的像素值按公式(2)逐个更新,
Bi(m,n)=(1-0.1)*Bi(m,n)+0.1*BUFi(m,n) (4)
其中P为正整数;
重复步骤四,当i=T时,执行步骤五;
步骤五,重复执行步骤二至步骤四,当j≥M时,即可实现背景的实时更新。
以下是发明人给出的具体实施例,需要说明的是本发明并不限于以下实施例,凡是在以下实施例基础上的等同变换或者替换均属于本发明的保护范围。
实施例:
该实施例中的视频的采样频率为25帧每秒,图像大小为720*288,将每个图像划成90*48个块区域,每块区域的大小为8*6,对相邻帧对应块内所有像素求差的绝对值和,阈值A为480,遵循上述技术方案,依次对连续输入的视频图像进行处理;
当累加计数器Ci等于50时,保存该块当前帧的所有像素值;
当累加计数器等于70时,此时P=1,且该块的背景有效标志Fi为0,用保存的该块完全更新背景,得到部分有效背景,同时,用当前帧更新该块保存的像素值,且该块对应背景有效标志位Fi设置为1。
当累加计数器Ci满足Ci=50+20*P,P≥2时,此时该块对应背景有效标志位Fi=1,则按比例更新背景,即可得到实时更新的背景。
Claims (2)
1.一种基于块累加统计的背景实时更新方法,其特征在于,该方法按下述步骤进行:
步骤一,初始化背景所有像素点的像素值为0,并在块坐标系X下,将背景划分成T个块,该T个块中的第i个块内像素点(m,n)的像素值为Bi(m,n),该块内所有像素点的的像素值用二维数组Bi[N]表示,该第i个块的有效背景标志为Fi,累加计数器为Ci,该第i个块内的所有像素点的像素值缓存区域为二维数组BUFi[N],缓存区域内像素点(m,n)的像素值为BUFi(m,n),将BUFi[N]中所有像素点的像素值初始化0,并且Fi和Ci初始化值均为0;
其中:T=(W/w)*(H/h),W为背景水平方向的像素,H为背景竖直方向的像素,背景的大小为W*H,w为第i个块的宽度,h为第i个块的高度;第i个块内像素点的总个数N=w*h;i=1,2,3...T;以第i个块的左下角为原点建立直角二维坐标系Y,m代表第i个块中任一像素点在坐标系Y下的横坐标,n代表第i个块中任一像素点在坐标系Y下的纵坐标,且m=0,1,2...w-1,n=0,1,2...h-1;
步骤二,在块坐标系X下,将第j+1帧图像划分为T个块,该T个块中的第i个块内像素点(m,n)的像素值为Fi,j+1(m,n),该块内所有像素点的像素值用二维数组Fi,j+1[N]表示;
同时在块坐标系X下,将第j帧图像划分为T个块,该T个块中的第i个块内像素点(m,n)的像素值为Fi,j(m,n),该块内所有像素点的像素值用二维数组Fi,j[N]表示;
其中:j≡1,2,3,...M-1,M为图像的总帧数,且M≥71;T=(W/w)*(H/h),第j+1帧图像和第j帧图像的水平方向的像素均为W,第j+1帧图像和第j帧图像竖直方向的像素均为H,第j+1帧图像和第j帧图像的大小均为W*H,第i个块的宽度为w,第i个块的高度为h;第i个块内像素点的总个数N=w*h;以第i个块的左下角为原点建立直角二维坐标系Y,m代表第i个块中任一像素点在坐标系Y下的横坐标,n代表第i个块中任一像素点在坐标系Y下的纵坐标,且m=0,1,2...w-1,n=0,1,2...h-1;
步骤三,计算第j+1帧图像中的第i个块与第j帧图像中的第i个块之间各相同坐标的像素点处的像素值差的绝对值之和SAD,其中SAD利用公式(1)计算得到,
当SAD小于或等于阈值A时,Ci加1;
当SAD大于阈值A时,Ci=0,BUFi(m,n)=0;
重复步骤三,当i=T时,执行步骤四;
步骤四,对背景中的第i个块的累加计数器Ci进行判断:
当Ci=50时,BUFi(m,n)=Fi,j+1(m,n);
当Ci=50+20*P,且Fi=0时,首先Bi(m,n)=BUFi(m,n),Fi=1,然后BUFi(m,n)=Fi,j+1(m,n);
当Ci=50+20*P,且Fi=1时,首先Bi(m,n)=(1-0.1)*Bi(m,n)+0.1*BUFi(m,n),然后BUFi(m,n)=Fi,j+1(m,n);
其中P为正整数;
重复步骤四,当i=T时,执行步骤五;
步骤五,重复执行步骤二至步骤四,当j≥M时,即可实现背景的实时更新。
2.如权利要求1所述的基于块累加统计的背景实时更新方法,其特征在于,所述阈值A的取值为10×块内像素点的总个数,即A=10×(w×h)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210099687.XA CN102646280B (zh) | 2012-04-06 | 2012-04-06 | 一种基于块累加统计的背景实时更新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210099687.XA CN102646280B (zh) | 2012-04-06 | 2012-04-06 | 一种基于块累加统计的背景实时更新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102646280A true CN102646280A (zh) | 2012-08-22 |
CN102646280B CN102646280B (zh) | 2014-12-03 |
Family
ID=46659092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210099687.XA Expired - Fee Related CN102646280B (zh) | 2012-04-06 | 2012-04-06 | 一种基于块累加统计的背景实时更新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102646280B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106952474A (zh) * | 2017-04-12 | 2017-07-14 | 湖南源信光电科技股份有限公司 | 基于运动车辆检测的车流量统计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101465955A (zh) * | 2009-01-05 | 2009-06-24 | 北京中星微电子有限公司 | 背景更新方法和装置 |
CN101533511A (zh) * | 2009-04-14 | 2009-09-16 | 清华大学 | 更新背景图像的方法和装置 |
CN102184547A (zh) * | 2011-03-28 | 2011-09-14 | 长安大学 | 一种基于视频的车辆逆行事件检测方法 |
-
2012
- 2012-04-06 CN CN201210099687.XA patent/CN102646280B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101465955A (zh) * | 2009-01-05 | 2009-06-24 | 北京中星微电子有限公司 | 背景更新方法和装置 |
CN101533511A (zh) * | 2009-04-14 | 2009-09-16 | 清华大学 | 更新背景图像的方法和装置 |
CN102184547A (zh) * | 2011-03-28 | 2011-09-14 | 长安大学 | 一种基于视频的车辆逆行事件检测方法 |
Non-Patent Citations (3)
Title |
---|
HE ZHIWEI等: "New Method of Background Update for Video-based Vehicle Detection", 《2004 IEEE INTELLIGENT TRANSPORTATION SYSTEMS CONFERENC》, 3 October 2004 (2004-10-03), pages 580 - 584 * |
卞建勇等: "基于帧差与统计法的运动目标背景图像提取算法", 《微电子学与计算机》, vol. 28, no. 4, 30 April 2011 (2011-04-30), pages 91 - 94 * |
樊晓亮等: "基于帧间差分的背景提取与更新算法", 《计算机工程》, vol. 37, no. 22, 30 November 2011 (2011-11-30), pages 159 - 160 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106952474A (zh) * | 2017-04-12 | 2017-07-14 | 湖南源信光电科技股份有限公司 | 基于运动车辆检测的车流量统计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102646280B (zh) | 2014-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2858008B1 (en) | Target detecting method and system | |
CN104463903B (zh) | 一种基于目标行为分析的行人图像实时检测方法 | |
US8077981B2 (en) | Providing virtual inserts using image tracking with camera and position sensors | |
CN102202164B (zh) | 基于运动估计的道路视频稳定方法 | |
CN102708571B (zh) | 视频中剧烈运动的检测方法及其装置 | |
CN103226834B (zh) | 一种图像运动目标特征点快速搜索方法 | |
CN101882217B (zh) | 视频图像的目标分类方法及装置 | |
CN104820996A (zh) | 一种基于视频的自适应分块的目标跟踪方法 | |
CN106530328B (zh) | 一种基于视频图像对运动目标检测和平滑跟随的方法 | |
CN101527044A (zh) | 多视频运动目标的自动分割与跟踪方法 | |
CN102883175A (zh) | 深度图提取、判断视频场景切换及深度图边缘优化方法 | |
CN102201122B (zh) | 一种运动捕捉的数据降噪方法、系统及运动捕捉系统 | |
CN105574894A (zh) | 一种运动物体特征点跟踪结果的筛选方法和系统 | |
CN104200492A (zh) | 基于轨迹约束的航拍视频目标自动检测跟踪方法 | |
CN102930559B (zh) | 图像处理的方法及装置 | |
CN102111530A (zh) | 移动物体检测装置与方法 | |
CN102855466A (zh) | 一种基于视频图像处理的人数统计方法 | |
CN102214301B (zh) | 自适应运动相关协作多目标跟踪方法 | |
CN105374049A (zh) | 一种基于稀疏光流法的多角点跟踪方法及装置 | |
CN105139419A (zh) | 一种足球比赛视频的球员及球体跟踪方法 | |
CN110458017B (zh) | 目标跟踪的尺度估计方法以及相关装置 | |
CN102006462A (zh) | 利用运动信息的快速监控视频增强方法及其实现装置 | |
CN102646280B (zh) | 一种基于块累加统计的背景实时更新方法 | |
CN108876807B (zh) | 一种实时星载卫星图像运动对象检测跟踪方法 | |
CN105631900A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141203 Termination date: 20160406 |
|
CF01 | Termination of patent right due to non-payment of annual fee |