function E=colour_extraction(I, k)

% supponiamo che I sia un'immagine a toni di grigio
[nx, ny]=size(I);
E=zeros(nx, ny);
soglia=5;
for i=1:nx-1
    for j=1:ny-1
   val=I(i,j);
   d1=abs(val-k);
    
   if d1<soglia
     E(i, j)=k;
        end
    end
end


