matlab - フーリエ変換の低周波係数を抽出

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

離散フーリエ変換のエネルギー圧縮プロパティを使用して、フーリエ変換の低周波値係数の 21x21 行列を抽出するにはどうすればよいですか?

前もって感謝します！

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)

左がオリジナル、右がコア周波数成分のみを選択したもの。