www.pudn.com > aasaaad.rar > first.m

imread bz_1217n3.bmp; 
a=ans; 
b=im2double(a); 
c=b; 
for i=21:3225; 
  for j=51:2100; 
     %b(i,j,1)=0; 
    %d(1)=b(i,j-2,1);      
    %d(2)=b(i,j-1,1);  
    %d(3)=b(i,j,1); 
    %d(4)=b(i,j+1,1); 
    %d(5)=b(i,j+2,1); 
   for t=1:101; 
     d(t)=b(i,j-51+t,1); 
   end 
     
    % for x=0:3; 
     %for y=x+1:4;   for y=0:4-x;       
     % if d(y)>=d(y-1); 
     %   e=d(y);  
     %   d(y)=d(y-1); 
      %  d(y-1)=e; 
     
    for x=1:100; 
     for y=x+1:101; 
      if d(y)>=d(x); 
        e=d(y);  
        d(y)=d(x); 
        d(x)=e; 
      end 
     end 
   end 
    c(i,j,1)=d(51); 
   
    %if b(i,j,1)                                      
    % c(i,j,1)=0; 
    %end 
  end 
end 
 
image(c);