image - 傅立叶和零填充

Question

我正在使用蒙版和离散傅里叶变换过滤图像，直到现在我有这个

A=double(imread('C:\Users\samsung\Documents\Lab Imagenes\CHE.jpg','jpg'));
B=[1 4 6 4 1; 4 16 24 16 4; 6 24 36 24 6; 4 16 24 16 4; 1 4 6 4 1];
F=(1/256).*(B);
DFT_A=fftshift(fft2(A));
imshow(DFT_A);
DFT_A_F=DFT_A.*F;
figure
imshow(DFT_A_F)

但是当我想看到部分结果时，我得到了这个错误

??? Error using ==> times
Matrix dimensions must agree.

Error in ==> fourier1 at 10 
DFT_A_F=DFT_A.*F;

我知道我需要对面具进行零填充，但我不知道该怎么做，请我需要帮助谢谢！

Answer 1

在您定义 DFT_A 之后，您想要的称为 'padarray' ：

padsize= [round(0.5*size(DFT_A,1)-0.5*size(F,1)) round(0.5*size(DFT_A,2)-0.5*size(F,2))];
F = padarray(F, padsize);
DFT_A_F=DFT_A.*F;
...

但是你为什么不这样做（因为 A 是一个二维矩阵，所以如果需要的话 rgb2gray ）：

DFT_A_F = conv2(A,B,'same');

它更快，因为您不需要将所有这些零相乘，并且应该得到相同的结果。