µÚ¶þ°æÊéÖеÚËÄÕÂͼ4.28ÔËÐнá¹ûºÍÊéÖÐÏÔʾµÄ²»Ò»Ñù£¬Çë¸ßÊÖ°ïÎÒ¿´Ò»Ï£¬³ÌÐòÊDz»ÊÇÓÐÎÊÌ⣿лл£¡
% ƵÓòÖÐÀÆÕÀ˹ÔöÇ¿
clear all;
close all;
clc;
f=imread('Fig4.28(a).jpg');
subplot(221),imshow(f);
f1=im2double(f);
[M,N]=size(f1);
m=fix(M/2); n=fix(N/2);
for x=1:M
for y=1:N
f2(x,y)=f1(x,y)*(-1)^(x+y);
end
end
f3=fft2(f2);
for u=1:M
for v=1:N
h(u,v)=-((u-m/2)^2+(v-n/2)^2); % ÀÆÕÀ˹Â˲¨Æ÷
f4(u,v)=h(u,v)*f3(u,v);
end
end
f5=ifft2(f4);
r=real(f5);
for x=1:M
for y=1:N
f6(x,y)=r(x,y)*(-1)^(x+y);
end
end
subplot(222),imshow(uint8(f6),[]);
% ¶ÔÀÆÕÀ˹±ä»»ºóµÄͼÏñÖØÐ±궨
mi=min(min(f6));
ma=max(max(f6));
T=ma-mi;
for i=1:M
for j=1:N
if f6(i,j)<0
f7(i,j)=0;
else
f7(i,j)=(f6(i,j)-mi)/T*255;
end
end
end
% ÏÔʾ±ê¶¨ºóµÄͼÏñ
subplot(223),imshow(uint8(f7),[]);
% ÔöÇ¿ºóµÄͼÏñ
f8=f1-f6;
subplot(224),imshow(f8,[]);
ÔËÐнá¹ûÖеĵÚÈý·ùͼ²»¶Ô£¬±³¾°²»Ó¦¸ÃÊǰ×É«µÄ£¬µ«²»Öª´íÔÚÄÄÁË¡£
![]()
Fig4.28(a).jpg
![]()
ÔËÐнá¹û
[ À´×Ô¿ÆÑмÒ×å ¿ØÖÆ&ͼÏñ¼Ò×å ]
[ Last edited by sunal on 2012-5-7 at 21:28 ] |