matlab巴特沃斯高通滤波器
生活随笔
收集整理的這篇文章主要介紹了
matlab巴特沃斯高通滤波器
小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
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')
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')
總結(jié)
以上是生活随笔為你收集整理的matlab巴特沃斯高通滤波器的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: chome(谷歌浏览器)上传文件崩溃/上
- 下一篇: 如何查看一个avi文件的播放帧率?和帧数