Discussion :

[Unik14]

Bonjour, après avoir manipulé l'équation de la chaleur sur une image en noir et blanc pour la débruiter, j'essaye de l'appliquer sur une image couleur mais l'algo me donne un résultat étrange, une image plus lumineuse et je ne suis pas sûr que le lissage soit vraiment bien réalisé.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
Im=imread('paysage.tif');
Im_Noisy = zeros(size(Im),'uint8');
Im_denoised = zeros(size(Im),'uint8');
for i=1:3
    B=Im(:,:,i);
    B=double(B);
    noise=randn(size(B));
    B=B+2*noise;
    Im_Noisy(:,:,i) = uint8(B);
end
figure, imshow(Im_Noisy); title('Image bruitée');
dt=0.2;
T=1;
for i=1:3
    u=double(Im_Noisy(:,:,i));
    [n m] = size(u);
    for t = 0:dt:T
        u_xx = u(:, [2:n n]) - 2*u + u(:, [1 1:n-1]);
        u_yy = u(:, [2:m m]) - 2*u + u(:, [1 1:m-1]);
        u = u + dt*(u_xx + u_yy);
    end
    Im_denoised(:,:,i) = int8(u);
end
figure, imagesc(Im_denoised); title ('Image lissée');

Merci

[Invité]

Bonsoir,

Cette ligne est fausse :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

        u_yy = u(:, [2:m m]) - 2*u + u(:, [1 1:m-1]);

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

        u_yy = u([2:m m], :) - 2*u + u([1 1:m-1], :);

Sinon tu calcules le gradient deux fois horizontalement. Et tu obtiens même une erreur dans le cas d'une image non carrée.

Une autre petite erreur :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

Im_denoised(:,:,i) = uint8(u);

[Unik14]

Merci !! ça fonctionne nickel désormais !