Извлечение коэффициентов низкочастотные преобразования Фурье

Question

Я применил преобразование Фурье на изображении следующим образом:

I = imread('img.jpg'); 
img = fftshift(I); 
F = fft2(img); 
magnitude = mat2gray(100*log(1+abs(fftshift(F)))); % Magnitude spectrum 
phase = mat2gray((angle(F)));      % Phase spectrum 

Использование энергии уплотнения свойства дискретного преобразования Фурье, как можно извлечь 21х21 матрицу низкого -частотные коэффициенты преобразования Фурье?

Заранее благодарен!

Answer 1

Вы могли бы попробовать что-то вроде этого:

I = imread('peppers.jpg'); 
img = fftshift(I); 
F = fftshift(fft2(img)); 

% magnitude = mat2gray(100*log(1+abs(F))); % Magnitude spectrum 
% phase = mat2gray((angle(F)));    % Phase spectrum 

[M N K] = size(F); 

L = 10; 

fsub(M,N,K)=0; 
fsub(M/2-L:M/2+L,N/2-L:N/2+L,1:K) = F(M/2-L:M/2+L,N/2-L:N/2+L,1:K); 


I2 = uint8(real(ifftshift(ifft2(ifftshift(fsub)))); 

figure 
subplot(121) 
imshow(I) 

subplot(122) 
imshow(I2) 

Панель слева является оригиналом, один справа после выбора только частотные компоненты ядра.