%greyscale morphological opening with 5x5 kernel

beeld=imread('blb1.jpg','jpg');
[r k dim3]=size(beeld);
A=beeld(:,:,1);
nieuw=A;


% erosion

for i=2:r-1
   for j=2:k-1
      %3x3 kernel
      if i>1 & j>1 & A(i-1,j-1)>A(i,j) & nieuw(i-1,j-1)>A(i,j)
         nieuw(i-1,j-1)=A(i,j);   
      end
      if i<r & j<k & A(i+1,j+1)>A(i,j) & nieuw(i+1,j+1)>A(i,j)
         nieuw(i+1,j+1)=A(i,j);   
      end
      if j>1 & A(i,j-1)>A(i,j) & nieuw(i,j-1)>A(i,j)
         nieuw(i,j-1)=A(i,j);   
      end
      if j<k & A(i,j+1)>A(i,j) & nieuw(i,j+1)>A(i,j)
         nieuw(i,j+1)=A(i,j);   
		end
		if i>1 & A(i-1,j)>A(i,j) & nieuw(i-1,j)>A(i,j)
         nieuw(i-1,j)=A(i,j);   
      end
      if i<r & A(i+1,j)>A(i,j) & nieuw(i+1,j)>A(i,j)
         nieuw(i+1,j)=A(i,j);   
		end
		if i>1 & j<k & A(i-1,j+1)>A(i,j) & nieuw(i-1,j+1)>A(i,j)
         nieuw(i-1,j+1)=A(i,j);   
      end
      if i<r & j>1 & A(i+1,j-1)>A(i,j) & nieuw(i+1,j-1)>A(i,j)
         nieuw(i+1,j-1)=A(i,j);   
		end
      %5x5 kernel, additional conditions
      if i>2 & j>2 & A(i-2,j-2)>A(i,j) & nieuw(i-2,j-2)>A(i,j)   
         nieuw(i-2,j-2)=A(i,j);   
      end
      if i<r-1 & j<k-1 & A(i+2,j+2)>A(i,j) & nieuw(i+2,j+2)>A(i,j)
         nieuw(i+2,j+2)=A(i,j);   
      end
      if j>2 & A(i,j-2)>A(i,j) & nieuw(i,j-2)>A(i,j)
         nieuw(i,j-2)=A(i,j);   
      end
      if j<k-1 & A(i,j+2)>A(i,j) & nieuw(i,j+2)>A(i,j)
         nieuw(i,j+2)=A(i,j);   
		end
		if i>2 & A(i-2,j)>A(i,j) & nieuw(i-2,j)>A(i,j)
         nieuw(i-2,j)=A(i,j);   
      end
      if i<r-1 & A(i+2,j)>A(i,j) & nieuw(i+2,j)>A(i,j)
         nieuw(i+2,j)=A(i,j);   
		end
		if i>2 & j<k-1 & A(i-2,j+2)>A(i,j) & nieuw(i-2,j+2)>A(i,j)
         nieuw(i-2,j+2)=A(i,j);   
      end
      if i<r-1 & j>2 & A(i+2,j-2)>A(i,j) & nieuw(i+2,j-2)>A(i,j)
         nieuw(i+2,j-2)=A(i,j);   
		end
      
      if i>2 & j>1 & A(i-2,j-1)>A(i,j) & nieuw(i-2,j-1)>A(i,j)   
         nieuw(i-2,j-1)=A(i,j);   
      end
      if i<r-1 & j<k & A(i+2,j+1)>A(i,j) & nieuw(i+2,j+1)>A(i,j)
         nieuw(i+2,j+1)=A(i,j);   
      end
      if i>2 & j<k & A(i-2,j+1)>A(i,j) & nieuw(i-2,j+1)>A(i,j)
         nieuw(i-2,j+1)=A(i,j);   
      end
      if j>1 & i<r-1 & A(i+2,j-1)>A(i,j) & nieuw(i+2,j-1)>A(i,j)
         nieuw(i+2,j-1)=A(i,j);   
		end
		if i>1 & j>2 & A(i-1,j-2)>A(i,j) & nieuw(i-1,j-2)>A(i,j)
         nieuw(i-1,j-2)=A(i,j);   
      end
      if i<r & j<k-1 & A(i+1,j+2)>A(i,j) & nieuw(i+1,j+2)>A(i,j)
         nieuw(i+1,j+2)=A(i,j);   
		end
		if i>1 & j<k-1 & A(i-1,j+2)>A(i,j) & nieuw(i-1,j+2)>A(i,j)
         nieuw(i-1,j+2)=A(i,j);   
      end
      if i<r & j>2 & A(i+1,j-2)>A(i,j) & nieuw(i+1,j-2)>A(i,j)
         nieuw(i+1,j-2)=A(i,j);   
		end

   end
end


A = nieuw;


% dilation
for i=1:r
   for j=1:k
      %3x3 kernel
      if i>1 & j>1 & A(i-1,j-1)<A(i,j) & nieuw(i-1,j-1)<A(i,j)   
         nieuw(i-1,j-1)=A(i,j);   
      end
      if i<r & j<k & A(i+1,j+1)<A(i,j) & nieuw(i+1,j+1)<A(i,j)
         nieuw(i+1,j+1)=A(i,j);   
      end
      if j>1 & A(i,j-1)<A(i,j) & nieuw(i,j-1)<A(i,j)
         nieuw(i,j-1)=A(i,j);   
      end
      if j<k & A(i,j+1)<A(i,j) & nieuw(i,j+1)<A(i,j)
         nieuw(i,j+1)=A(i,j);   
		end
		if i>1 & A(i-1,j)<A(i,j) & nieuw(i-1,j)<A(i,j)
         nieuw(i-1,j)=A(i,j);   
      end
      if i<r & A(i+1,j)<A(i,j) & nieuw(i+1,j)<A(i,j)
         nieuw(i+1,j)=A(i,j);   
		end
		if i>1 & j<k & A(i-1,j+1)<A(i,j) & nieuw(i-1,j+1)<A(i,j)
         nieuw(i-1,j+1)=A(i,j);   
      end
      if i<r & j>1 & A(i+1,j-1)<A(i,j) & nieuw(i+1,j-1)<A(i,j)
         nieuw(i+1,j-1)=A(i,j);   
      end
      %5x5 kernel, additional conditions
      if i>2 & j>2 & A(i-2,j-2)<A(i,j) & nieuw(i-2,j-2)<A(i,j)   
         nieuw(i-2,j-2)=A(i,j);   
      end
      if i<r-1 & j<k-1 & A(i+2,j+2)<A(i,j) & nieuw(i+2,j+2)<A(i,j)
         nieuw(i+2,j+2)=A(i,j);   
      end
      if j>2 & A(i,j-2)<A(i,j) & nieuw(i,j-2)<A(i,j)
         nieuw(i,j-2)=A(i,j);   
      end
      if j<k-1 & A(i,j+2)<A(i,j) & nieuw(i,j+2)<A(i,j)
         nieuw(i,j+2)=A(i,j);   
		end
		if i>2 & A(i-2,j)<A(i,j) & nieuw(i-2,j)<A(i,j)
         nieuw(i-2,j)=A(i,j);   
      end
      if i<r-1 & A(i+2,j)<A(i,j) & nieuw(i+2,j)<A(i,j)
         nieuw(i+2,j)=A(i,j);   
		end
		if i>2 & j<k-1 & A(i-2,j+2)<A(i,j) & nieuw(i-2,j+2)<A(i,j)
         nieuw(i-2,j+2)=A(i,j);   
      end
      if i<r-1 & j>2 & A(i+2,j-2)<A(i,j) & nieuw(i+2,j-2)<A(i,j)
         nieuw(i+2,j-2)=A(i,j);   
		end
      
      if i>2 & j>1 & A(i-2,j-1)<A(i,j) & nieuw(i-2,j-1)<A(i,j)   
         nieuw(i-2,j-1)=A(i,j);   
      end
      if i<r-1 & j<k & A(i+2,j+1)<A(i,j) & nieuw(i+2,j+1)<A(i,j)
         nieuw(i+2,j+1)=A(i,j);   
      end
      if i>2 & j<k & A(i-2,j+1)<A(i,j) & nieuw(i-2,j+1)<A(i,j)
         nieuw(i-2,j+1)=A(i,j);   
      end
      if j>1 & i<r-1 & A(i+2,j-1)<A(i,j) & nieuw(i+2,j-1)<A(i,j)
         nieuw(i+2,j-1)=A(i,j);   
		end
		if i>1 & j>2 & A(i-1,j-2)<A(i,j) & nieuw(i-1,j-2)<A(i,j)
         nieuw(i-1,j-2)=A(i,j);   
      end
      if i<r & j<k-1 & A(i+1,j+2)<A(i,j) & nieuw(i+1,j+2)<A(i,j)
         nieuw(i+1,j+2)=A(i,j);   
		end
		if i>1 & j<k-1 & A(i-1,j+2)<A(i,j) & nieuw(i-1,j+2)<A(i,j)
         nieuw(i-1,j+2)=A(i,j);   
      end
      if i<r & j>2 & A(i+1,j-2)<A(i,j) & nieuw(i+1,j-2)<A(i,j)
         nieuw(i+1,j-2)=A(i,j);   
		end
   end
end


imwrite(nieuw,'blb1openingvierk5.jpg','jpg');