%這個程序是模擬4f系統(tǒng)濾波處理，在濾波面加菲涅爾濾波片
clc;
clear?all;

%模擬4f系統(tǒng)對圖像進行濾波處理
V=imread（‘D:\3.bmp‘）;?????%讀入字母V
V=rgb2gray（V）;
V=im2double（V）;
subplot（121）;imshow（V）;title（‘字母V‘）;
V1=fftshift（fft2（V））;?????%?將圖像頻譜低頻部分移到中央，高頻部分移到四周
subplot（122）;imshow（abs（V1）*16/max（max（abs（V1））））;title（‘字母V經(jīng)物面到達頻譜面后的頻譜圖‘）;?%字母V經(jīng)物面到達頻譜面后的頻譜圖


r=16;
a=4*sqrt（r）;
b=0;
%生成菲涅爾濾波片
F=ones（256256）;
for?k=1:30
????p=6*sqrt（k*r）;
for?x=1:256
for?y=1:256
????r1=sqrt（（x-128）^2+（y-128）^2）;
if?r1>=a&r1<=p
????F（xy）=b;
end
end
end
a=p;
b=abs（b-1）;
end
subplot（121）;imshow（F）;title（‘菲涅爾濾波片‘）;
figuresubplot（131）;imshow（F）;title（‘菲涅爾濾波片‘）;
B2=V1.*F;?????????????????%?濾波
subplot（132）;imshow（abs（B2）*64/max（max（abs（B2））））;title（‘經(jīng)頻譜面濾波后頻譜圖‘）;
C2=ifft2（B2）;???
subplot（133）;imshow（abs（C2）/max（max（abs（C2））））;title（‘濾波后的字母‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????236278??2011-01-04?18:04??光學(xué)4f系統(tǒng)菲涅爾波帶片\1-1.bmp

?????文件?????236278??2011-01-04?18:04??光學(xué)4f系統(tǒng)菲涅爾波帶片\1-2.bmp

?????文件?????196664??2010-12-31?10:49??光學(xué)4f系統(tǒng)菲涅爾波帶片\1.bmp

?????文件?????236278??2011-01-04?18:05??光學(xué)4f系統(tǒng)菲涅爾波帶片\2-1.bmp

?????文件?????236278??2011-01-04?18:05??光學(xué)4f系統(tǒng)菲涅爾波帶片\2-2.bmp

?????文件?????196664??2010-12-31?11:07??光學(xué)4f系統(tǒng)菲涅爾波帶片\2.bmp

?????文件?????236278??2011-01-04?18:06??光學(xué)4f系統(tǒng)菲涅爾波帶片\3-1.bmp

?????文件?????236278??2011-01-04?18:06??光學(xué)4f系統(tǒng)菲涅爾波帶片\3-2.bmp

?????文件?????196664??2010-12-31?11:07??光學(xué)4f系統(tǒng)菲涅爾波帶片\3.bmp

?????文件????????964??2011-01-04?18:05??光學(xué)4f系統(tǒng)菲涅爾波帶片\4f.m

?????目錄??????????0??2011-02-26?14:13??光學(xué)4f系統(tǒng)菲涅爾波帶片

-----------?---------??----------?-----??----

??????????????2008624????????????????????11