matlab巴特沃斯高通滤波器

I1=imread('原图.jpg');
figure(1);imshow(I1);
title('原图');
m=double(I1);
f=fft2(m);
f=fftshift(f);
[N1,N2]=size(f);
n1=round(N1/2);
n2=round(N2/2);
n=2;d0=10;
for i=1:N1
    for j=1:N2
        d=sqrt((i-n1)^2+(j-n2)^2);
        h=(1/(1+(d0/d)^(2*n)))+0.5;
        y(i,j)=h*f(i,j);
    end
end
y=ifftshift(y);
A=ifft2(y);
B=uint8(real(A));
figure(2);imshow(B);
title('巴特沃斯高通滤波d0=10');
imwrite(B,'巴特沃斯高通滤波d0=10.jpg')

评论 2 您还未登录,请先 登录 后发表或查看评论
©️2022 CSDN 皮肤主题:大白 设计师:CSDN官方博客 返回首页
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值