avi=aviread（‘SEM1.avi‘）;

%Seting?Parameters

T=20;

[N]=size（avi2）;
?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）;
%?????pau