avi=aviread('temp4.avi');
%Seting Parameters
T=20;
[N]=size(avi,2);
im_mean=double(avi(1).cdata);
for k=1:N
avi(k).cdata=rgb2gray(avi(k).cdata);
end
% figure;
% colormap(gray);
%
% %Defferitial image between coutinous two images
% for k=2:N
% imagesc((avi(k-1).cdata-avi(k).cdata)>T);
% pause(.2)
% end
%
%
% %Accumulating deffertial image.......
% int_img=0;
% figure;
% colormap(gray);
% for k=2:N
% int_img=int_img+double((avi(1).cdata-avi(k).cdata)>T);
% imagesc(int_img>1);
% pause(.2)
% end
%Deferitial image with a reference images
im_mean=double(avi(1).cdata);
colormap(gray);
for k=2:N
im_mean=im_mean+double(avi(k).cdata);
end
im_mean=uint8(im_mean/N);
imagesc(im_mean);
figure;
colormap(gray);
for k=1:N
imagesc((im_mean-avi(k).cdata)>T);
pause(.2)
end
%Accumulating deffertial image.......
int_img=double((im_mean-avi(k).cdata)>T);
figure;
colormap(gray);
for k=1:N
int_img=int_img+double((im_mean-avi(k).cdata)>T);
imagesc(int_img);
pause(.2)
end
