均值建模法

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