巴特沃斯高通滤波器
function im5 = freqfilter_btw_Hp(im,D0,n)% 巴特沃斯高通滤波器
% input:
% M,N:频域滤波器的尺寸
% D0:带阻滤波器的截止频率
% n :阶次
% output:
% H:M x N的矩阵,表示频域滤波器矩阵,数据类型为double,
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);
if d==0
h=0;
else
h = 1/(1 + 0.414 * (D0 / d)^(2*n));% 巴特沃斯高通滤波器
end
result(i,j) = h * im3(i,j);
end
end
result = ifftshift(result); % 反中心化
im4 = ifft2(result); % 反变换
im5 = im2uint8(real(im4)); % 滤波图像
end
页:
[1]