Średni filtr do wygładzania obrazów w programie Matlab

Question

Potrzebuję przetestować podstawowe techniki przetwarzania obrazów w programie Matlab. Muszę przetestować i porównać szczególnie dwa typy filtrów: filtr średni i filtr medianowy.

Aby wygładzić obraz przy użyciu filtrowania median, dostępna jest wspaniała funkcja medfilt2 z przybornika do przetwarzania obrazu. Czy istnieje jakaś podobna funkcja dla średniego filtra? Albo jak użyć funkcji filter2, aby utworzyć filtr średniej?

Jedną z najważniejszych rzeczy dla mnie jest możliwość ustawienia promienia filtra. To znaczy. za filtrem mediany, jeśli chcę [3 x 3] promienia (maskę), wystarczy użyć

imSmoothed = medfilt2(img, [3 3]); 

chciałbym osiągnąć coś podobnego do średniej filtra.

Answer 1

h = fspecial('average', n); 
filter2(h, img); 

Zobacz doc fspecial: h = fspecial('average', n) powraca filtr uśredniający. n to wektor 1 na 2 określający liczbę wierszy i kolumn w h.

Answer 2

I = imread('peppers.png'); 
H = fspecial('average', [5 5]); 
I = imfilter(I, H); 
imshow(I) 

Należy zauważyć, że filtry mogą być zastosowane do zdjęcia intensywność (2D) za pomocą matryc filter2, a na zdjęciach wielowymiarowych (obrazy macierzy RGB i 3D) imfilter używany.

Również na procesorach Intel, imfilter może użyć biblioteki Intel Integrated Performance Primitives (IPP), aby przyspieszyć wykonanie.

Answer 3

Widzę, że dobre odpowiedzi już zostały podane, ale pomyślałem, że byłoby miło po prostu podać sposób, aby wykonać średnie filtrowanie w MATLAB-ie bez specjalnych funkcji lub skrzynek narzędziowych. Jest to również bardzo dobre dla zrozumienia, jak dokładnie działa ten proces, ponieważ jest to wymagane do jawnego ustawienia jądra splotu. Średni kernel filtr jest na szczęście bardzo proste:

I = imread(...) 
kernel = ones(3, 3)/9; % 3x3 mean kernel 
J = conv2(I, kernel, 'same'); % Convolve keeping size of I 

Należy pamiętać, że w przypadku obrazów kolorowych trzeba by zastosować to do każdego z kanałów w obrazie.

Answer 4

i splot jest określona poprzez pomnożenie w przekształcić domenę:

conv2(x,y) = fftshift(ifft2(fft2(x).*fft2(y))) 

jeśli jeden kanał jest uważany za ... więcej kanałów To musi być zrobione każdy kanał

Answer 5

f=imread(...); 

h=fspecial('average', [3 3]); 
g= imfilter(f, h); 
imshow(g);