MATLAB代码步骤

Question

我想在过滤器索贝尔和图像之间的频域的产物，我不能让 过滤的图像什么是我的问题在此代码

[A,B]=size(f); 
%[C,D]=size(fspecial('sobel')); 
h=fspecial('gaussian'); 
[C,D]=size(h); 
hp=freqz2(h,size(f,2),size(f,1)); 
'%padarray(h,[ceil(((A-C)/2)) ceil(((B-D)/2))]); 
'%hp=hp(1:254,1:256); 
HP=fft2(hp); 
HP=ifftshift(HP); 
F=fft2(f); 
G=HP.*F; % product of transfer function and image in frequency domain 
g=real(ifft2(G)); 
imshow(g) % do not work 

Answer 1

你不需要HP=fft2(hp); ，hp已经在频域。请删除该行，并使用

HP=ifftshift(hp); 
F=fft2(f); 
G=HP.*F; 
g=real(ifft2(G)); 
imshow(g)