untitled.m

高斯建模实现图像背景生成

clear all;  
path='E:\毕业设计（09-02-23）\new0001\';
file_I=imfinfo([path 'new0001_1'],'bmp');
I=zeros(240,320,numel(1:5:282));
for i=1:5:282              
    q=imread([path,'new0001_', num2str(i),'.BMP']);
    I(:,:,floor(i/5)+1)=q;
end   
[h w]=size(I(:,:,1));
a=mean(I,3); 
b=std(I,1,3);
B=zeros(size(I));
C0=zeros([h w]);
C1=ones([h w]);
B(:,:,1)=a;
for i=1:h
    for j=1:w
        for k=2:56
            if  abs(I(i,j,k)-a(i,j))<5*b
                c=C0;
            else
                c=C1;
            end
            B(i,j,k)=(1-0.055)B(i,j,k-1)+0.055*c
        end
    end
end
subplot(2,2,1);imshow(B(:,:,1),[]);title('No1');
subplot(2,2,2);imshow(B(:,:,2),[]);title('No30');
subplot(2,2,3);imshow(B(:,:,30),[]);title('No56');
subplot(2,2,4);imshow(B(:,:,56),[]);title('背景图像');