tigaofenbian.m

均值滤波,加权滤波,中值滤波的MATLAB实例,用于图象及信号的处理

rgb=imread('Fig6.48(d).jpg');
[r,c,d]=size(rgb);
hr=2*r;
hc=2*c;
for i=1:r-1;
    for j=1:c-1
        a1=(rgb(i,j,1)+rgb(i,j+1,1))/2;
        rgb1(2*i-1,2*j,1)=a1;
        b1=(rgb(i,j,1)+rgb(i+1,j,1))/2;
        rgb1(2*i,2*j-1,1)=b1;
        rgb1(2*i-1,2*j-1,1)=rgb(i,j,1);
        rgb1(2*i,2*j,1)=(a1+b1)/2;
        
         a2=(rgb(i,j,2)+rgb(i,j+1,2))/2;
        rgb1(2*i-1,2*j,2)=a1;
        b2=(rgb(i,j,2)+rgb(i+1,j,2))/2;
        rgb1(2*i,2*j-1,2)=b2;
        rgb1(2*i-1,2*j-1,2)=rgb(i,j,2);
         rgb1(2*i,2*j,2)=(a2+b2)/2;
        
         a3=(rgb(i,j,3)+rgb(i,j+1,3))/2;
        rgb1(2*i-1,2*j,3)=a1;
        b3=(rgb(i,j,3)+rgb(i+1,j,3))/2;
        rgb1(2*i,2*j-1,3)=b3;
        rgb1(2*i-1,2*j-1,3)=rgb(i,j,3);
         rgb1(2*i,2*j,3)=(a3+b3)/2;
       
    end
end
%for i=1:r-1;
   % for j=1:c-1
     %   a=(rgb(i,j,1)+rgb(i+1,j,1))/2;
     %   rgb1(2*i-1,2*j-1,1)=rgb(i,j,1);
     %   rgb1(2*i,2*j,1)=a;
   % end
%
%subplot(1,2,1);
%for i=1:hr;
  % for j=1:hc
       rgb3=rgb1+0.3*rgb1;
   %end
%end
figure(1);
imshow(rgb3);
%subplot(1,2,2);
%imshow(rgb);