三帧差分 matlab,三帧差分法

hhhhhhhh

三帧差分法的问题

大家帮忙看看这个三帧差分法的代码哪里出了问题？参照论坛上的相邻帧差法，改造成三帧差分法。

clc;

clear;

targetavi='SampleVideo.avi'; %%原始视频

resultavi='result.avi'; %%检测结果视频

%%%%%%%%%%%%%%%%%%%%%%读取视频

mov=aviread(targetavi);

%mov=mmreader(targetavi);

fnum=size(mov,2);

%%%%%%%%%%%%%%%%%%%%%%建立结果视频

aviobj = avifile(resultavi);

aviobj.Quality = 100;

aviobj.Fps = 15;

%http://www.wendangwang.compression='Indeo5';

http://www.wendangwang.compression='None';

%%%%%%%%%%%%%%%%%%%%%%帧间差分法

for i=2:fnum

x=mov(i-1).cdata(:,:,:);

y=mov(i).cdata(:,:,:);

z=mov(i+1).cdata(:,:,:);

m=rgb2gray(x);

m=medfilt2(m);

n=rgb2gray(y);

n=medfilt2(n);

o=rgb2gray(z);

o=medfilt2(o);

q=im2double(n);

w=im2double(m);

p=im2double(o);

c=q-w;%%前两帧差分

b=p-q;%%后两帧差分

t=40; %%阈值，此值可以调节

t=t/256;%%转化为double型数据

k=find(abs(c)>=t);%%find函数作用是找到图c中的值大于t的点坐标