Finden Sie geeignete Notch-Filter, um Rauschen von Bild zu entfernen

Question

Ich möchte periodisches Rauschen aus einem Bild mit Notch-Filter beseitigen. i die folgenden Schritte aus:

I = imread('...img....'); 
ft = fftshift(fft2(I)); 
[m,n] = size(ft); 
filt = ones(m,n); 
%filt(......) = 0; % my problem is here 
ft = ft .* filt; 
ifft_ = ifft2(ifftshift(ft)); 

so i don't know what exactly to set to zero to get the proper result.

Answer 1

Werfen Sie einen Blick auf dieses Bild.

Dies ist die FFT eines Bildes mit ähnlichen Lärmproblemen zu Ihnen. Diese roten Kreise weisen auf die "Spitzen" in der Frequenzdomäne hin, die dieser Art von Rauschen zugeordnet ist. Schau dir dein Bild an und schau, ob du es finden kannst.

Diese Spitzen sind, was Sie beseitigen wollen. Mit anderen Worten, Ihr Filter sollte an diesen Stellen einen Wert von Null und an allen anderen Stellen einen Wert von eins haben. Wenn Sie dann in der folgenden Zeile multiplizieren, werden Sie diese Informationen eliminieren und das damit verbundene Rauschen drastisch reduzieren.

Answer 2

Die Grundidee hinter dem Notch-Filter besteht darin, eine kleine Anzahl von Frequenzen im FFT-Bild auf Null zu setzen, was zu Geräuschen führt.

Um sie zu finden, können Sie bei der FFT-Größe aussehen:

F = abs(ft2); 
F = log(F+1); 
imshow(F,[]) 

Es ist möglich, die Spitzen im Frequenzbild zu sehen, dass, an folgenden Stellen befinden:

pnts2del ={[88,155],[101,92],[138,71],[125,134]}; 

Wenn Sie jeden dieser Punkte nullen Umgebung werden Sie eine sauberere Bild erhalten:

G = 60; 
for ii=1:length(pnts2del) 
    pnt = pnts2del{ii}; 
    filt(pnt(1),pnt(2)-G:pnt(2)+G) = 0; 
    filt(pnt(1)-G:pnt(1)+G,pnt(2)) = 0; 
    filt(pnt(1)-1:pnt(1)+1,pnt(2)-1:pnt(2)+1) = 0; 
end 

Ergebnis (vor und nach, und FFT aktualisiert):