频率域图像增强
用傅里叶变换表示的函数特征可以完全通过傅里叶反变换进行重建而不丢失任何信息。
吉布斯现象Gibbs phenomenon（又叫吉布斯效应）：将具有不连续点的周期函数（如矩形脉冲)进行傅立叶级数展开后，选取有限项进行合成。当选取的项数越多，在所合成的波形中出现的峰起越靠近原信号的不连续点。当选取的项数很大时，该峰起值趋于一个常数，大约等于总跳变值的9%。这种现象称为吉布斯现象。
一、傅里叶变换
l = fft2(x); %快速傅里叶变换l = fft2(x,m,n);
x为输入图像；m和n分别用于将x的第一和第二维规整到指定的长度。当m和n均为2的整数次幂时算法的执行速度要比m和n均为素数时快。
l1 = abs(l); %计算l的幅度谱l2 = angle(l); %计算l的相位谱Y = fftshift(l); %频谱平移l = ifft2(x);%快速傅里叶逆变换l =iff2(x,m,n);
%% fftshift 对数变换,所应用的图片本身很简单，就只有黑白2种颜色clcclearf = imread('.\images\dipum_images_ch04\Fig0403(a)(image).tif');imshow(f)title('原始图像')imfinfo('.\images\dipum_images_ch04\Fig0403(a)(image).tif');%此处如果用Imfinfo(f)就会报错fft%没有居中的傅里叶频谱F=fft2(f);%进行二维快速傅里叶变换，其结果和DFT的一样，只是计算机的计算速度变快了而已，因而叫fftS=abs(F);%求傅里叶变换后的幅值figure,subplot(121),imshow(S,[]),title('傅里叶频谱图像1');%title函数一定要放在坐标显示的下一句才有效。subplot(122),imshow(S),title('傅里叶频谱图像2')；%当没有第二个参数时，显示的图像为竖线加一些孤立的黑点%居中的傅里叶频谱Fc=fftshift(F);%将频谱图像原点移至图像矩形中间S1=abs(Fc);figure,subplot(121),imshow(S1,[]);%加了第二个参数后显示的图像正常%使用对数后视觉增强后的傅里叶频谱S2=log(1+S1);subplot(122),imshow(S2,[]);
原始图像： 傅里叶频谱图：
居中频谱和对数后频谱：
理想低通：
%%%=============理想低通========================function out = imidealflpf(f,freq)%imidealflpf函数 构造高斯频域低通滤波器[M,N] = size(l);out = ones(M,N);for i = 1:Mfor j = 1:Nif(sqrt(((i-M/2)^2+(j-N/2)^2))>freq)out(i,i) = 0;endendend
高斯低通：
%%%=============高斯低通========================function out = imgausslpf(l,sigma)%imgaussflpf函数 构造高斯频域低通滤波器[M,N] = size(l);out = ones(M,N);for i = 1:Mfor j = 1:Nout(i,j) = exp(-((i-M/2)^2+(j-N/2)^2)/2/sigma^2);endend
