Filtro medio per uniformare le immagini in Matlab

Question

Ho bisogno di testare alcune tecniche di elaborazione delle immagini di base in Matlab. Ho bisogno di testare e confrontare soprattutto due tipi di filtri: filtro medio e filtro mediano.

Per uniformare l'immagine utilizzando il filtraggio mediano, esiste una grande funzione medfilt2 dalla casella degli strumenti di elaborazione delle immagini. Esiste una funzione simile per il filtro medio? O come utilizzare la funzione filter2 per creare il filtro medio?

Una delle cose più importanti per me è avere la possibilità di impostare il raggio del filtro. Cioè per filtro mediano, se voglio la [3 x 3] raggio (maschera), mi basta usare

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

Vorrei realizzare qualcosa di simile per il filtro media.

Answer 1

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

Vedi doc fspecial: h = fspecial('average', n) restituisce un filtro di media. n è un vettore 1-by-2 che specifica il numero di righe e colonne in h.

Answer 2

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

noti che i filtri possono essere applicati alle immagini di intensità (matrici 2D) utilizzando filter2, mentre sulle immagini multidimensionali (immagini RGB o matrici 3D) imfilter viene utilizzato.

Anche sui processori Intel, imfilter è possibile utilizzare la libreria Intel Integrated Performance Primitives (IPP) per accelerare l'esecuzione.

Answer 3

Vedo che sono già state fornite buone risposte, ma ho pensato che sarebbe stato utile dare un modo per eseguire il filtraggio medio in MATLAB senza utilizzare funzioni speciali o toolbox. Questo è anche molto utile per capire esattamente come funziona il processo, in quanto è necessario impostare esplicitamente il kernel della convoluzione. Il kernel del filtro media è fortunatamente molto semplice:

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

Si noti che per le immagini a colori si dovrebbe applicare questo per ciascuno dei canali nell'immagine.

Answer 4

e la circonvoluzione è definito attraverso una moltiplicazione in trasformare dominio:

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

se un canale è considerato ... per più canali questo deve essere fatto ogni canale

Answer 5

f=imread(...); 

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