Débruitage par seuillage en ondelettes Résolu/Fermé

Question

Bonjour, 
salut à tous    
 
mon objectif est d'appliquer la transformation en ondelettes sur l'image bruité, ensuite appliquer l'opération de seuillage, 
 
ci-aprés mon code en matlab: 
 
clear all;close all;clc; 
%Charger et bruiter une image 
C1=double(imread('C:\Users\user\Desktop\ACTIVITEMULTI\activité_1.jpg')); 
I=rgb2gray(C1/255); 
n = size(I,1); 
sigma = 0.06 * (max(I()-min(I()); 
In = I + sigma*randn(n); 
subplot(1,2,1); 
imagesc(I); axis image; axis off; 
subplot(1,2,2); 
imagesc(In); axis image; axis off; 
colormap gray(256); 
 
%Transformée en ondelettes orthogonale 2D et seuillage 
% Transformée en ondelettes 2D 
Iw = FWT2_PO(I,Jmin,qmf); 
IwT = Iw; IwT = Iw .* (abs(Iw)>T); 
I_ond = IWT2_PO(IwT,Jmin,qmf); % transformée inverse  
% dessines des transformée originale et seuillées 
subplot(1,2,1); 
plot_wavelet(Iw, Jmin); colormap gray(256); 
subplot(1,2,2); 
plot_wavelet(IwT, Jmin); colormap gray(256); 
 
Mais je n'arrive pas à simuler mon code j'obtiens l'erreur suivante: 
??? Error using ==> plus 
Matrix dimensions must agree. 
 
Error in ==> debruitond at 7 
In = I + sigma*randn(n); 
 
j'essaye de terminer la taille de mon image I en tapant size(I) 
 
ans = 
 
   105    71 
 
JE serais reconnaissante si quelqu'un peux me déterminer l'erreur.


Configuration: Windows Vista / Safari 535.2

Char Snipeur · Answer

Salut.
Je pense que ton erreur vient de l'utilisation des min et des max. Ces fonction retournent  des vecteurs lorsque tu fournis des matrices.
Or, j'ai l'impression que tu veux que sigma soit un réel. De même, es tu sur que ton image est carré ?

tchichik · Answer

salut,
merci pour la réponse,

pour moi je cherche o faire bruité mes images par un bruit gaussien de variance variable par exemple sigma=10, 20 ou 3O puis appliqué les ondelettes en matière de débruitage tel que les ondelettes de symlet, de coifflet, de daubechies orthogonales et biorthogonale.

Pour mon image je vous ai fournie sa dimension (105, 71), elle doit etre 256 256 par exemple??

de toute facon je serais reconnaissante si vous me répondez encore une fois pour fixer d'ou vient l'erreur exactement ou si vous me prososez une autre mthode pour ce faire.
Merci

tchichik · Answer

SALUT :)

j'ai arrivé a trouvé mon ondelette adéquate ci joint mon code

%Load original image. 

X=imread(['base_originale\I' num2str(2) '.jpg']);
X=rgb2gray(X); 
X=double(X);

% Generate noisy image.

sigma =0.1 * ( max(max(X))-min(min(X)) );
x = X + sigma*randn(size(X));


% Find default values. In this case fixed form threshold
% is used with estimation of level noise, thresholding
% mode is soft and the approximation coefficients are 
% kept.
[thr,sorh,keepapp] = ddencmp('den','wv',x);


% De-noise image using global thresholding option.
xd = wdencmp('gbl',x,'Coif4',1,i,sorh,keepapp);

%figure,image(wcodemat(X)), title('Original Image');
%figure,image(wcodemat(xd)), title('De-Noised Image');

colormap(pink(255)), sm = size(X,1);
subplot(221), image(wcodemat(X,sm)), title('Original Image')
subplot(222), image(wcodemat(x,sm)), title('Noisy Image')
subplot(223), image(wcodemat(xd,sm)), title('De-Noised Image')



mon code est exécutable mais l'affichage des images n'est pas en niveau de gris 

autre chose il me donne les axes (des x et des y) numérotés.

Alors comment je peux avoir un affichage en Niveau de gris et les axes non numérotés?

Merci d'avance.