?? l5_3.m
字號:
%程序L5_3.m:同態濾波器%
clear
[image_0,map]=bmpread('L5_3.bmp'); %讀取圖像
image_1=log(image_0+1); %取自然對數
image_2=fft2(image_1); %做二維傅立葉變換
%產生高通Butterworth濾波器%
n=3;
D0=0.05*pi;
rh=0.8;
rl=0.5;
[row,col]=size(image_2);
for k=1:1:row
for l=1:1:col
D1(k,l)=sqrt((k^2+l^2));
H(k,l)=rl+(rh/(1+(D0/D1(k,l))^(2*n)));
end
end
image_3=(image_2.*H); %輸入圖像通過濾波器
image_4=ifft2(image_3); %做二維反傅立葉變換
image_5=(exp(image_4)-1); %取指數函數
figure(1)
imshow(image_0,map) %顯示原圖像
figure(2)
imshow(real(image_5),map) %顯示處理后的圖像
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -