车流量检测

车流量检测

%% 用再次改进的Guass迭代法求背景
% 导入视频得到用平均值求得的背景
run（‘pingjunzhi‘）

% 初始化背景
% 这里用平均值求得的背景作为初始化的背景
% 首先，将所有的点值转换为灰度，并计算其均值、方差
for i = 1:100
    mov_final（i）.cdata = rgb2gray（mov_final（i）.cdata）;
end

back=rgb2gray（back）;

for i = 1:100
    mov_final（i）.cdata = im2double（mov_final（i）.cdata）;
end

back=im2double（back）;

for i=1:vidHeight
    for j=1:vidWidth
        sum_temp=0;
        for k=1:100
            sum_temp=sum_temp+（mov_final（k）.cdata（ij）-back（ij））.^2;
        end
        squaresigma（ij）=sum_temp./100;
    end
end

% 参数初始化
alfa=0.05;
yuzhi1=0.5;
yuzhi2=2.5;

% 开始对前两幅图片做处理
mov（1:nframes） = ...
    struct（‘cdata‘ zeros（vidHeight vidWidth 3 ‘uint8‘）...
           ‘colormap‘ []）;

%对第一幅图进行处理
mov（1）.cdata=read（mov_original  1）;
mov（1）.cdata=rgb2gray（mov（1）.cdata）;
mov（1）.cdata=im2double（mov（1）.cdata）;
P=1./sqrt（2.*pi.*squaresigma）.*（-（mov（1）.cdata-back）.^2./（2.*squaresigma））;
flag=ones（vidHeightvidWidth）;
%干扰点处理
for i=2:vidHeight-1
    for j=2:vidWidth-1
        P_n（1）=mov（1）.cdata（i-1j-1）;
        P_n（2）=mov（1）.cdata（i-1j）;
        P_n（3）=mov（1）.cdata（i-1j+1）;
        P_n（4）=mov（1）.cdata（ij-1）;
        P_n（5）=mov（1）.cdata（ij+1）;
        P_n（6）=mov（1）.cdata（i+1j-1）;
        P_n（7）=mov（1）.cdata（i+1j）;
        P_n（8）=mov（1）.cdata（i+1j+1）;
        if P_n            flag（ij）=1;
        else
            flag（ij）=0;
        end
    end
end
%计算新的背景
for i= vidHeight
    for j= vidWidth
        if P（ij）<=yuzhi1&&flag（ij）==1
            m（ij1）=1;
            back（ij）=alfa*mov（1）.cdata（ij）+（1-alfa）*back（ij）;
        elseif P（ij）>yuzhi1
            m（ij1）=0;
        end
    end
end
%计算新的方差
for i=1:vidHeight
    for j=1:vidWidth
        if P（ij）<=yuzhi1&&flag（ij）==1
            squaresigma（ij）=alfa*（mov（1）.cdata（ij）-back（ij））^2+（1-alfa）*squaresigma（ij）;
        end
    end
end
    
%对第二幅图进行处理
mov（2）.cdata=read（mov_original  2）;
mov（2）.cdata=rgb2gray（mov（2）.cdata）;
mov（2）.cdata=im2double（mov（2）.cdata）;
P=1./sqrt（2.*pi.*squaresigma）.*（-（mov（2）.cdata-back）.^2./（2.*squaresigma））;
flag=ones（vidHeightvidWidth）;
%干扰点处理
for i=2:vidHeight-1
    for j=2:vidWidth-1
        P_n（1）=mov（2）.cdata（i-1j-1）;
        P_n（2）=mov（2）.cdata（i-1j）;
        P_n（3）=mov（2）.cdata（i-1j+1）;
        P_n（4）=mov（2）.cdata（ij-1）;
        P_n（5）=mov（2）.cdata（ij+1）;
        P_n（6）=mov（2）.cdata（i+1j-1）;
        P_n（7）=mov（2）.cdata（i+1j）;
        P_n（8）=mov（2）.cdata（i+1j+1）;
        if P_n            flag（ij）=1;
        else
            flag（ij）=0;
        end
    end
end
%计算新的背景
for i= vidHeight
    for j= vidWidth
        if P（ij）<=yuzhi1&&flag（ij）==1
            m（ij2）=1;
            back（ij）=alfa*mov（2）.cdata（ij）+（1-alfa）*back（ij）;
        elseif P（ij）>yuzhi1
            m（ij2）=0;
        end
    end
end
%计算新的方差
for i=1:vidHeight
    for j=1:vidWidth
        if P（ij）<=yuzhi1&&flag（ij）==1
            squaresigma（ij）=alfa*（mov（2）.cdata（ij）-back（ij））^2+（1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      14169  2011-05-02 04:56  代码\a.jpg

     文件      33935  2011-05-02 00:11  代码\aagback.png

     文件       4874  2011-05-01 20:46  代码\adadvanced_Guass.m

     文件       5263  2011-05-01 13:41  代码\Advanced_Guass.m

     文件      33935  2011-05-02 00:11  代码\agback.png

     文件      33935  2011-05-01 23:41  代码\agback1.png

     文件      12816  2011-05-02 04:58  代码\b.jpg

     文件   22984316  2011-05-01 18:50  代码\B0.avi

     文件      33935  2011-05-02 00:01  代码\back_gray.png

     文件     304182  2011-05-02 23:19  代码\back_mhf.bmp

     文件       2846  2011-05-02 20:46  代码\goodrestest1.m

     文件       1359  2011-05-02 23:06  代码\goodrestest1.rar

     文件       1225  2011-05-01 20:47  代码\pingjunzhi.m

     文件       3132  2011-05-03 00:39  代码\test1.asv

     文件       2862  2011-05-03 01:41  代码\test1.m

     文件        259  2011-05-02 00:45  代码\test2.asv

     文件       1578  2011-05-02 21:31  代码\test2.m

     文件       2917  2011-05-03 00:44  代码\test3.asv

     文件       2933  2011-05-03 00:15  代码\test3.m

     文件       2576  2011-05-03 02:46  代码\test5.m

     文件        780  2011-05-02 22:34  代码\test_BG.asv

     文件        841  2011-05-02 22:48  代码\test_BG.m

    I.A.SH.      5120  2013-05-08 12:16  代码\Thumbs.db

     文件       1683  2011-05-02 09:10  代码\tracking.asv

     文件        360  2011-05-03 00:16  代码\公共部分.png

     文件        124  2011-05-03 02:55  代码\分割结果.png

     文件        224  2011-05-03 00:10  代码\前后帧差.png

     文件     102454  2011-05-02 23:24  代码\原始图像.bmp

     文件       2995  2011-05-03 02:54  代码\原始图像.png

     文件       1746  2011-05-02 14:55  代码\参考\bftest1.m

............此处省略18个文件信息