程序代码如下: axes(handles.axes2); x=(handles.img); f=double(x); k=fft2(f); g=fftshift(k); [M,N]=size(g); nn=2; d0=25; m=fix(M/2);n=fix(N/2); for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2); if d<=d0 ...
程序代码如下:
axes(handles.axes2);
x=(handles.img);
f=double(x);
k=fft2(f);
g=fftshift(k);
[M,N]=size(g);
nn=2;
d0=25;
m=fix(M/2);n=fix(N/2);
for i=1:M
for j=1:N
d=sqrt((i-m)^2+(j-n)^2);
if d<=d0
h=0;
else h=1;
end
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
y2=ifft2(result);
y3=uint8(real(y2));
imshow(y3);
原来的图是彩色图,应该是个三维矩阵吧
f=double(rgb2gray(x));
转为灰度图处理吧