資源描述:
《高通濾波-低通濾波matlab程序.doc》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、高通濾波matlab程序:clearall;I1=imread('E:matlab圖片硬幣.jpg');subplot(1,2,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=50;%濾波器截止頻率,濾波半徑fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);%計算高通濾波轉(zhuǎn)換函數(shù)ifd<=d0h=0+0.5;elseh=1+0.5;en
2、dy(i,j)=h*f(i,j);endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A));subplot(1,2,2);imshow(B);title('高通濾波');低通濾波matlab程序:clearall;I1=imread('E:matlab圖片硬幣.jpg');subplot(1,2,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=
3、2;d0=50;%濾波器截止頻率,濾波半徑fori=1:N1forj=1:N2d=sqrt((i-n1)^2+(j-n2)^2);%計算低通濾波轉(zhuǎn)換函數(shù)ifd<=d0h=1;elseh=0;endy(i,j)=h*f(i,j);endendy=ifftshift(y);A=ifft2(y);B=uint8(real(A));subplot(1,2,2);imshow(B);title('低通濾波');