generate_data.m

%% Regenration des variables

%chargement image
canards=imread('canards.png');
#conversion en niveaux de gris
im=double(rgb2gray(canards));
save("canards.mat",'im')

%Initialisation de la graine pour le generateur de nombres aleatoires
seed_number=2020;
rand("seed",seed_number);
randn("seed",seed_number);

%Ajout de bruit blanc Gaussian
SignalToNoiseRatio=5e2; #Rapport Signal a Bruit
sigma=norm(double(im))/SignalToNoiseRatio;
im_b=im+sigma*randn(size(im));
save("canards_gauss.mat",'im_b');

%Ajout de bruit sel
ratio = 4e-1; %ratio de pixel bruites
bruit_sel = rand(size(im))<ratio;
max_type = double(intmax(class(canards))); %get the maximum of data type
im_sel=im+bruit_sel*max_type;
save("canards_sel.mat",'im_sel');

%Ajout de bruit sel
bruit_poivre = rand(size(im))>ratio;
im_poivre=im.*bruit_poivre;
save("canards_poivre.mat",'im_poivre')

#affichage
figure(1)
subplot(221)
colormap('gray')
imagesc(im)
title('Originale')
subplot(222)
colormap('gray')
imagesc(im_b)
title('Bruit Gaussien')
subplot(223)
colormap('gray')
imagesc(im_sel)
title('Bruit Sel')
subplot(224)
colormap('gray')
imagesc(im_poivre)
title('Bruit Poivre')