ImageLab 频域处理图像达到增强效果
floatdata2d Re,Im;
ucdata2d Ret,Ret2;
floatdata2d Retf,Ret2f;
SetDim(Re, 512, 512);
SetDim(Im, 512, 512);
SetDim(Ima, 512, 512);
Load(Lena.bmp, R, G, B, A);
Re=R;
Retf=Re;
FFT2D(Re,Im);
FFTShift2D(Re,Im);
@@
for (i = 0 : 1 : 512, j = 0 : 1 : 512)
if(sqrt((i-255)^2+(j-255)^2)<50)
{
Re[i][j] = 0;
Im[i][j] = 0;
}
else
{
}
@@
FFTShift2D(Re,Im);
IFFT2D(Re,Im);
Ret = Re;
##Save(lenafft.raw, Re, flt);
##ToUc(Ret,Ret2,0.0,0.0,0.0);
Ret2f=Ret;
Operand(Ret2f,0.9,multi);
Operand(Ret2f,Retf,add);
Ret2=Ret2f;
Save(Lenafft.bmp, gray, Ret2);