1.使用paddedsize获得填充参数:PQ = paddedsize(size(f));%如果输入是彩色图像,必须要灰度化rgb2gray。
2.得到使用填充的傅里叶变化:F = fft2(f, PQ(1), PQ(2));
3. 生成一个大小为PQ(1) X PQ(2) 的滤波函数H。如果该滤波函数已居中,使用前要令H = fftshift(H)。
4.将变换乘以滤波函数:G = H.*F;
5.获得G的傅里叶逆变换的实部:g = real(ifft2(G));
6.将左上部分的矩形剪切为原来尺寸大小:g = g(1:size(f,1), 1:size(f, 2));