交通流量监测MATLAB源代码

这是我写的2011.5数学建模中B题车流量检测的matlab代码，压缩文件里面提供了当时的视频，程序可以直接运行，在车流量不是非常大的情况下统计的还是非常准确，有兴趣的朋友可以看看

%% 用再次改进的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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

     文件      12734  2011-05-02 23:23  trafficflow\增强对比度二值化结果.bmp

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

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

     文件        692  2011-05-03 00:10  trafficflow\和差.png

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

     文件        540  2011-05-03 00:10  trafficflow\背景差.png

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

     文件        407  2011-05-03 00:35  trafficflow\检测.png

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