function  J=salepepe(I, q)
% aggiunge rumore casuale di tipo sale e pepe ad I
% q stabilisce la percentuale di rumore che vogliamo aggiungere.

if ~isa(I, 'double')
    I=double(I);
end
J=I;
n_p=rand(size(I)); % si stabilisce il numero di pixel da 
n_i=find(n_p<=q);  % alterare con il rimore
n=zeros(size(I));
n(n_i)=sign(randn(size(n_i))); % n matrice del rumore: nei pixel n_i ci metto -1 o 1; 
Imax=max(I(:));
Imin=min(I(:));
salt=find(n==-1); % trova in n i pixel che valcgono -1
pepper=find(n==1); % trova in n i pixel che valgono 1
J(salt)=Imax;
J(pepper)=Imin;
nsalt=length(salt)
npepper=length(pepper)