cm1=imread('twins.tif');
imshow(cm1);
r=cm1(:,:,1);
g=cm1(:,:,2);
b=cm1(:,:,3);
r=fftshift(fft2(r));
g=fftshift(fft2(g));
b=fftshift(fft2(b));
[x,y]=meshgrid(-128:127,-128:127);
z=sqrt(x.^2+y.^2);
c=(z<15);%高通>低通<
r2=r.*c;
g2=g.*c;
b2=b.*c;
cm2(:,:,1)=r2;
cm2(:,:,2)=g2;
cm2(:,:,3)=b2;
figure;fftshow(cm2,'log');
cfhi=ifft2(cm2);
figure;fftshow(cfhi,'abs');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
bl=lbutter(cm1,15,1);%butter低通
r3=r.*bl;
g3=g.*bl;
b3=b.*bl;
cfbl(:,:,1)=r3;
cfbl(:,:,2)=g3;
cfbl(:,:,3)=b3;
figure;fftshow(cfbl,'log');
cfbli=ifft2(cfbl);
figure;fftshow(cfbli,'abs');
bh=hbutter(cm1,15,1);%butter高通
r4=r.*bh;
g4=g.*bh;
b4=b.*bh;
cfbh(:,:,1)=r4;
cfbh(:,:,2)=g4;
cfbh(:,:,3)=b4;
%figure;fftshow(cfbh,'log');
cfbhi=ifft2(cfbh);
%figure;fftshow(cfbhi,'abs');
--
ID:
暱稱:
在
虧了 個妞
偷偷來了 次
被抓到了 次
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.130.19.97