高斯低通滤波器
function im5 = freqfilter_gaussian_lp(im,D0)if ~isa(im,'double')
im1 = double(im)/255;
end
im2 = fft2(im1); % 傅里叶变换
im3 = fftshift(im2);% 中心化
= size(im3);
n1 = fix(N1 / 2);
n2 = fix(N2 / 2);
for i = 1:N1
for j = 2:N2
D = sqrt((i-n1)^2+(j-n2)^2);
h = exp(-D.^2/2/D0/D0);% 高斯低通滤波器
result(i,j) = h * im3(i,j);
end
end
result = ifftshift(result); % 反中心化
im4 = ifft2(result); % 反变换
im5 = im2uint8(real(im4)); % 滤波图像
end
页:
[1]