Débruitage par ondelette

[invitec925f48a]

bonjour tout le monde

Voila mon problème J'ai une image et je veux la débruiter pour cela j'ai appliqué les ondelettes :/
Tout d'abord j'ai appliqué wavedec2 et après afin d'obtenir le THR ''la valeur de thresholding'' j'ai utilisée la formule de donoho ''delta = MAD / 0.6745'' '' thr = delta * sqrt( 2 * log(n))'' sur les details horizontaux mais malheureusement je ne suis pas arrivée a une solution parfaite :/
voila mon code et une petite explication des commentaires

clear all 
I=imread('cameraman.tif');
n = prod( size(I) );
I=double(I);
Ib=I+25*randn(size(I));% ajouter un bruit
[C,S] = wavedec2(Ib,2,'bior3.7');
DH = detcoef2('all',C,S,1);% extraire seulement  les details horizontaux de niveau 1 et puis appliquer donoho
delta = median( abs(DH) ) / 0.6745;
thr = delta * sqrt(2*log(n));   
NC = wthcoef2('bior3.7',C,S,DH,thr,s) % ici je veux appliquer le soft thresholding mais :/ 
 X = waverec2(C, S, 'bior3.7'); % comment je peux reconstruire mon image après l'utilisation de  wthcoef2
 figure(2)
 imagesc(X);axis off;colormap(gray)

-----

[invitec925f48a]

c'est vous avez la reponse svp n'hesitez pas de me repondre