%greyscale morphological erosion with 1x7 kernel

beeld=imread('blb1.jpg','jpg');
[r k dim3]=size(beeld);
A=beeld(:,:,1);
nieuw=A;
for i=1:r
   for j=1:k
      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 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 j>3 & A(i,j-3)>A(i,j) & nieuw(i,j-3)>A(i,j)
         nieuw(i,j-3)=A(i,j);   
      end
      if j<k-2 & A(i,j+3)>A(i,j) & nieuw(i,j+3)>A(i,j)
         nieuw(i,j+3)=A(i,j);   
		end

   end
end

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