帧差法matlab代码

mov=aviread（‘C:\Documents and Settings\Administrator\桌面\照片\cup.avi‘）;
temp=size（mov）; 
fnum=temp（2）;
for i=1:fnum
    strtemp=strcat（int2str（i）‘.‘‘JPG‘）;
    imwrite（mov（i）.cdata（:::）strtemp）;
end          %  从视频中提出每一帧图像
o=1;
e=‘.jpg‘;
for i=1:15
    u=o-1;v=o-2;
    m=int2str（o）;
    n=int2str（u）;h=int2str（v）;
    s=strcat（me）;%%把字符串b与m连接后在连接e得到图像文件存储的位置
    m=imread（s）;%%从S处把图像读取出来
    m=rgb2gray（m）;%%将图像M灰度化
   
m=medfilt2（m[33]）;
if（o>=3）
        s=strcat（ne）;
        n=imread（s）;
        n=rgb2gray（n）;
        s=strcat（he）;
        h=imread（s）;
        h=rgb2gray（h）;
        
        n=medfilt2（n[33]）;
        h=medfilt2（h[33]）;

        q=im2double（m）;%%将图像数组转换为double型
        w=im2double（n）; g=im2double（h）;
        c=q-w;j=w-g;
        
        th=10/255;
        
        k=find（abs（c）>=th）;
        c（k）=1;
        k=find（abs（c）        c（k）=0;  
        c=bwareaopen（c15）;
          
        se90=strel （‘line‘390）;se0=strel （‘line‘30）;