




版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、各類濾波器的MATLA醍序清單各類濾波器的 MATLAB 程序一、 理 想低通濾波器IA=imread('lena.bmp'); f1,f2=freqspace(size(IA),'meshgrid');Hd=ones(size(IA);r=sqrt(f1.2+f2.2);Hd(r>0.2)=0;Y=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(
2、uint8(Ia);figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')二、理想高通濾波器IA=imread('lena.bmp'); f1,f2=freqspace(size(IA),'meshgrid');Hd=ones(size(IA);r=sqrt(f1.2+f2.2);Hd(r<0.2)=0;Y=fft2(double(IA);Y=fftshift(Y
3、);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')三、 Butterworth 低通濾波器IA=imread('lena.bmp'); f1,f2=f
4、reqspace(size(IA),'meshgrid'); D=0.3;r=f1.2+f2.2;n=4;for i=1:size(IA,1)for j=1:size(IA,2)t=r(i,j)/(D*D);Hd(i,j)=1/(tn+1);end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia); figuresur
5、f(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')四、 Butterworth 高通濾波器IA=imread('lena.bmp');f1,f2=freqspace(size(IA),'meshgrid');D=0.3;r=f1.2+f2.2;n=4;for i=1:size(IA,1)for j=1:size(IA,2)t=(D*D)/r(i,j);Hd(i,j)=1/(tn+1)
6、;endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figure五、高 斯低通濾波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');D=100/size(IA,1);r=f1
7、.2+f2.2;Hd=ones(size(IA);for i=1:size(IA,1)for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=exp(-t);endendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figuresurf(Hd,'Facecolor','interp',
8、39;Edgecolor','none','Facelighting','phong')六、高 斯高通濾波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D=0.3;r=f1.2+f2.2;for i=1:size(IA,1)for j=1:size(IA,2) t=r(i,j)/(D*D); Hd(i,j)=1-exp(-t);en
9、dendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(uint8(Ia);figure七、梯 形低通濾波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D0=0.1
10、;D1=0.4; r=sqrt(f1.2+f2.2); Hd=zeros(size(IA);Hd(r<D0)=1;for i=1:size(IA,1)for j=1:size(IA,2) if r(i,j)>=D0 & r(i,j)<=D1 Hd(i,j)=(D1-r(i,j)/(D1-D0); end end end Y=fft2(double(IA); Y=fftshift(Y); Ya=Y.*Hd; Ya=ifftshift(Ya); Ia=real(ifft2(Ya); figure subplot(2,2,1),imshow(uint8(IA); subpl
11、ot(2,2,2),imshow(uint8(Ia); figure surf(Hd,'Facecolor','interp','Edgecolor','none','Facelighting','phong')八、梯 形高通濾波器IA=imread('lena.bmp');IB=imread('babarra.bmp'); f1,f2=freqspace(size(IA),'meshgrid');%D=100/size(IA,1);D0=0.1;D1
12、=0.4;r=sqrt(f1.2+f2.2);Hd=ones(size(IA);Hd(r<D1)=0;for i=1:size(IA,1)for j=1:size(IA,2)if r(i,j)>=D0 & r(i,j)<=D1Hd(i,j)=(D0-r(i,j)/(D0-D1);endendendY=fft2(double(IA);Y=fftshift(Y);Ya=Y.*Hd;Ya=ifftshift(Ya);Ia=real(ifft2(Ya);figuresubplot(2,2,1),imshow(uint8(IA); subplot(2,2,2),imshow(u
13、int8(Ia);figure九、用 其他方法編寫的理想低通、理想高通、Butterworth 低通、同態(tài)濾波程序1、 理想低通i1=imread('lena.bmp');i2=imnoise(i1,'salt & pepper',0.1);f=double(i2);k=fft2(f);g=fftshift(k);N1,N2=size(g);d0=50;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1for j=1:N2 d=sqrt(i-u0)2+(j-v0)2);if d<=d0h=1;elseh=0;e
14、ndy(i,j)=g(i,j)*h;endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figuresubplot(2,2,1),imshow(uint8(i1);subplot(2,2,2),imshow(uint8(i2);subplot(2,2,3),imshow(uint8(E2);2、 理想高通i1=imread('lena.bmp');i2=imnoise(i1,'salt & pepper',0.1); f=double(i2);k=fft2(f);g=fftshift(k);N1,N2=size(g)
15、;n=2;d0=10;u0=floor(N1/2)+1;v0=floor(N2/2)+1;for i=1:N1for j=1:N2d=sqrt(i-u0)2+(j-v0)2);if d<=d0h=0;else h=1;endy(i,j)=g(i,j)*h;endendy=ifftshift(y);E1=ifft2(y);E2=real(E1);figure subplot(2,2,1),imshow(uint8(i1); subplot(2,2,2),imshow(uint8(i2); subplot(2,2,3),imshow(uint8(E2);3、 Butterworth 低通i1
16、=imread('lena.bmp'); i2=imnoise(i1,'salt & pepper',0.1); f=double(i2);k=fft2(f); g=fftshift(k);N1,N2=size(g);n=2; d0=50; u0=floor(N1/2)+1; v0=floor(N2/2)+1;for i=1:N1for j=1:N2 d=sqrt(i-u0)2+(j-v0)2); h=1/(1+(d/d0)(2*n); y(i,j)=g(i,j)*h;end end y=ifftshift(y); E1=ifft2(y); E2=real(E1); figure subplot(2,2,1),imshow(uint8(i1); subplot(2,2,2),imshow(uint8(i2); subplot(2,2,3),imshow(uint8(E2);4、 同態(tài)濾波I=rgb2gray(imread('fabric00.b
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年度二手車寄售代理服務合同
- 2025年度吊車租賃與施工現(xiàn)場管理服務協(xié)議
- 2025年車貸居間服務收費標準合同范本
- 建立適應性強的財務體系計劃
- 學生領導力培養(yǎng)課程計劃
- 教學工作創(chuàng)新方向計劃
- 年度戰(zhàn)略性談判的技巧分享計劃
- 市場調研與分析工作安排計劃
- 特色課程開發(fā)與實施方案計劃
- 高危作業(yè)的安全保障措施優(yōu)化計劃
- 元宇宙視域下非遺保護與傳播途徑探究
- 2025年買賣個人房屋合同(4篇)
- 2025代運營合同范本
- 武漢2025年湖北武漢理工大學管理人員招聘筆試歷年參考題庫附帶答案詳解
- 家庭燃氣和煤氣防火安全
- 第十一章《功和機械能》達標測試卷(含答案)2024-2025學年度人教版物理八年級下冊
- 大學學生宿舍管理員工作培訓
- 初三物理常識試卷單選題100道及答案
- 初中新課標培訓課件
- 辦公用品價格清單
- 公司銀行貸款申請書范文
評論
0/150
提交評論