La traduzione di questa pagina non è aggiornata. Fai clic qui per vedere l'ultima versione in inglese.

Ricostruzione di immagini con le wavelet biortogonali

Questo esempio mostra come l’ordine in cui vengono applicati i filtri wavelet biortogonali possa influire sulla ricostruzione dell’immagine.

Generare i filtri di analisi e sintesi per la wavelet bior3.5. Caricare e visualizzare un'immagine.

[LoD,HiD,LoR,HiR] = wfilters('bior3.5');
load woman
imagesc(X)
colormap gray

I filtri di analisi LoD e HiD hanno 5 momenti di fuga. I filtri di sintesi LoR e HiR hanno 3 momenti di fuga. Eseguire una scomposizione wavelet a cinque livelli dell'immagine utilizzando i filtri di analisi.

[c1,s1] = wavedec2(X,5,LoD,HiD);

Trovare la soglia che mantiene solo i coefficienti wavelet con ampiezza nel 10 per cento superiore. Utilizzare la soglia per impostare a 0 il 90 percento dei coefficienti inferiori.

frac = 0.1;
c1sort = sort(abs(c1),'desc');
num = numel(c1);
thr = c1sort(floor(num*frac));
c1new = c1.*(abs(c1)>=thr);

Ricostruire l'immagine utilizzando i filtri di sintesi e i coefficienti con la soglia impostata. Visualizzare la ricostruzione.

X1 = waverec2(c1new,s1,LoR,HiR);
figure
imagesc(X1)
colormap gray

Eseguire una scomposizione wavelet a cinque livelli dell'immagine utilizzando i filtri di sintesi.

[c2,s2] = wavedec2(X,5,LoR,HiR);

Trovare la soglia che mantiene solo i coefficienti wavelet con ampiezza nel 10 per cento superiore. Utilizzare la soglia per impostare a 0 il 90 percento dei coefficienti inferiori

frac = 0.1;
c2sort = sort(abs(c2),'desc');
num = numel(c2sort);
thr = c2sort(floor(num*frac));
c2new = c2.*(abs(c2)>=thr);

Ricostruire l'immagine utilizzando i filtri di sintesi e i coefficienti con la soglia impostata. Visualizzare la ricostruzione. La scomposizione con un filtro che ha 3 momenti di fuga e la ricostruzione con un filtro che ha 5 momenti di fuga restituisce una ricostruzione scadente.

X2 = waverec2(c2new,s2,LoD,HiD);
figure
imagesc(X2)
colormap gray

Vedi anche

biorwavf | biorfilt | wfilters | waverec2 | wavedec2