ImageLab FFT频域理想滤波器振铃效应
floatdata2d Re,Im;
ucdata2d Ret;
SetDim(Re, 512, 512);
SetDim(Im, 512, 512);
SetDim(Ima, 512, 512);
Load(Lena.bmp, R, G, B, A);
Re=R;
FFT2D(Re,Im);
FFTShift2D(Re,Im);
@@
for (i = 0 : 1 : 512, j = 0 : 1 : 512)
if(sqrt((i-255)^2+(j-255)^2)>40)
{
Re[i][j] = 0;
Im[i][j] = 0;
}
else
{
}
@@
FFTShift2D(Re,Im);
IFFT2D(Re,Im);
Ret = Re;
Save(LenaRing.bmp, gray, Ret);