网站优化首页付款,cdr平面设计教程,wordpress 添加新字段,微信里的小程序游戏from#xff1a;https://www.cnblogs.com/tenderwx/p/5245859.html
1.理解二维傅里叶变换的定义
1.1二维傅里叶变换
二维Fourier变换: 逆变换#xff1a; 1.2二维离散傅里叶变换
一个图像尺寸为MN的 函数的离散傅里叶变换由以下等式给出#xff1a; 其中 和。其中变量u和…fromhttps://www.cnblogs.com/tenderwx/p/5245859.html
1.理解二维傅里叶变换的定义
1.1二维傅里叶变换
二维Fourier变换: 逆变换 1.2二维离散傅里叶变换
一个图像尺寸为M×N的 函数的离散傅里叶变换由以下等式给出 其中 和。其中变量u和v用于确定它们的频率频域系统是由所张成的坐标系其中和用做频率变量。空间域是由f(x,y)所张成的坐标系。可以得到频谱系统在频谱图四角处沿和方向的频谱分量均为0。
离散傅里叶逆变换由下式给出 令R和I分别表示F的实部和需部则傅里叶频谱相位角功率谱幅度定义如下 1.3用FFT计算二维离散傅里叶变换
二维离散傅里叶变换的定义为 二维离散傅里叶变换可通过两次一维离散傅里叶变换来实现:
1作一维N点DFT对每个m做一次共M次 2作M点的DFT对每个k做一次共N次 这两次离散傅里叶变换都可以用快速算法求得,若M和N都是2的幂,则可使用基二FFT算法,所需要乘法次数为 而直接计算二维离散傅里叶变换所需的乘法次数为MNMN当M和N比较大时用用FFT运算可节约很多运算量。
1.3图像傅里叶变换的物理意义
图像的频率是表征图像中灰度变化剧烈程度的指标是灰度在平面空间上的梯度。如大面积的沙漠在图像中是一片灰度变化缓慢的区域对应的频率值很低而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域对应的频率值较高。傅里叶变换在实际中有非常明显的物理意义设f是一个能量有限的模拟信号则其傅里叶变换就表示f的频谱。从纯粹的数学意义上看傅里叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看傅里叶变换是将图像从空间域转换到频率域其逆变换是将图像从频率域转换到空间域。换句话说傅里叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数。
傅里叶逆变换是将图像的频率分布函数变换为灰度分布函数傅里叶变换以前图像未压缩的位图是由对在连续空间现实空间上的采样得到一系列点的集合通常用一个二维矩阵表示空间上各点记为zf(x,y)。又因空间是三维的图像是二维的因此空间中物体在另一个维度上的关系就必须由梯度来表示这样我们才能通过观察图像得知物体在三维空间中的对应关系。
傅里叶频谱图上我们看到的明暗不一的亮点其意义是指图像上某一点与邻域点差异的强弱即梯度的大小也即该点的频率的大小可以这么理解图像中的低频部分指低梯度的点高频部分相反。一般来讲梯度大则该点的亮度强否则该点亮度弱。这样通过观察傅里叶变换后的频谱图也叫功率图我们就可以直观地看出图像的能量分布如果频谱图中暗的点数更多那么实际图像是比较柔和的因为各点与邻域差异都不大梯度相对较小反之如果频谱图中亮的点数多那么实际图像一定是尖锐的、边界分明且边界两边像素差异较大的。
对频谱移频到原点以后可以看出图像的频率分布是以原点为圆心对称分布的。将频谱移频到圆心除了可以清晰地看出图像频率分布以外还有一个好处它可以分离出有周期性规律的干扰信号比如正弦干扰。一幅频谱图如果带有正弦干扰移频到原点上就可以看出除了中心以外还存在以另一点为中心、对称分布的亮点集合这个集合就是干扰噪音产生的。这时可以很直观的通过在该位置放置带阻滤波器消除干扰。
2.二维傅里叶变换有哪些性质
2.1二维离散傅里叶变换的性质
1分离性
二维离散傅里叶变换具有分离性 分离性质的主要优点是可借助一系列一维傅里叶变换分两步求得。第1步沿着的每一行取变换将其结果乘以1/N,取得二维函数;第2步沿着的每一列取变换再将结果乘以1/N,就得到了。这种方法是先行后列。如果采用先列后行的顺序其结果相同。
如图 对逆变换f(x,y)也可以类似地分两步进行。 2平移性
傅里叶变换和逆变换对的位移性质是指 由乘以指数项并取其乘积的傅立叶变换使频率平面的原点位移至。同样地以指数项乘以并取其反变换将空间域平面的原点位移至当N/2时指数项为 即为 这样用(xy)乘以就可以将的傅里叶变换原点移动到N*N频率方阵的中心这样才能看到整个谱图。另外对的平移不影响其傅里叶变换的幅值。
此外与连续二维傅里叶变换一样二维离散傅里叶变换也具有周期性、共轭对称性、线性、旋转性、相关定理、卷积定理、比例性等性质。这些性质在分析及处理图像时有重要意义。
2.2二维离散傅里叶变换图像性质
1、图像经过二维傅里叶变换后其变换系数矩阵具有如下性质若变换矩阵原点设在中心其频谱能量集中分布在变换系数短阵的中心附近(图中阴影区)。若所用的二维傅里叶变换矩阵的原点设在左上角那么图像信号能量将集中在系数矩阵的四个角上。这是由二维傅里叶变换本身性质决定的。同时也表明一股图像能量集中低频区域。
2、图像灰度变化缓慢的区域对应它变换后的低频分量部分图像灰度呈阶跃变化的区域对应变换后的高频分量部分。除颗粒噪音外图像细节的边缘、轮廓处都是灰度变化突变区域它们都具有变换后的高频分量特征。 3.任给一幅图像对其进行二维傅里叶变换和逆变换 原图 二维傅里叶变换 逆变换 频谱图 原图 加入高斯躁声 加入椒盐躁声 对高斯躁声中值滤波 对椒盐躁声中值滤波 对高斯躁声算术均值滤波 对椒盐躁声算术均值滤波 4.附录
4.1matlab代码
1程序一 1 [i,lcmp]imread(F:/123.jpg);%读取图像 显示图像2 3 subplot(2,2,1),imshow(i,lcmp);4 5 title(original);6 7 iiim2double(i); %将图像矩阵类型转换为double图像计算很多是不能用整型的没有这个会报错 如果不用这个就必须转化为灰度图8 9 i1 fft2(ii); %傅里叶变换
10
11 i2 fftshift(i1); %将变换的频率图像四角移动到中心原来良的部分在四角 现在移动中心便于后面的处理
12
13 i3log(abs(i2)); %显示中心低频部分加对数是为了更好的显示
14
15 subplot(2,2,2),imshow(i3,[]);
16
17 title(Fourier);
18
19 mapcolormap(lcmp); %取色谱
20
21 imwrite(i3,map,f:/ffttank.bmp); %将上面i3输入到ffttank文件中
22
23 i5 real(ifft2(ifftshift(i2))); %频域的图反变换到空域 并取实部
24
25 i6 im2uint8(mat2gray(i5)); %取其灰度图
26
27 imwrite(i6,map,f:/tank2.bmp,bmp); %利用灰度图和原来取得颜色模板 还原图像
28
29 subplot(2,2,3),imshow(i6);
30
31 title(anti-Fourier);
32
33 i7rgb2gray(i);
34
35 i8fft2(i7);%对灰色图才能归一化。因为那是2维矩阵彩色图是3维矩阵需要转化为2维灰图
36
37 mfftshift(i8); %直流分量移到频谱中心
38
39 %RRreal(m); %取傅立叶变换的实部
40
41 %IIimag(m); %取傅立叶变换的虚部
42
43 Aabs(m);%计算频谱幅值
44
45 %Asqrt(RR.^2II.^2);
46
47 A(A-min(min(A)))/(max(max(A))-min(min(A)))*225; %归一化
48
49 subplot(2,2,4),imshow(A); %显示原图像
50
51 colorbar; %显示图像的颜色条
52
53 title(FFT spectrum); %图像命名
54
55 2程序二 1 mimread(F:/123.jpg);2 3 Mrgb2gray(m); %滤波函数都是对二维灰度图Tif可直接滤波4 5 subplot(3,3,1)6 7 imshow(M);%显示原始图像8 9 title(original)
10
11 P1imnoise(M,gaussian,0.02); %加入高斯躁声
12
13 subplot(3,3,2)
14
15 imshow(P1) %加入高斯躁声后显示图像
16
17 title(gaussian noise);
18
19 P2imnoise(M,salt pepper,0.02); %加入椒盐躁声
20
21 subplot(3,3,3)
22
23 imshow(P2) %%加入椒盐躁声后显示图像
24
25 title(salt pepper noise);
26
27 gmedfilt2(P1); %对高斯躁声中值滤波
28
29 subplot(3,3,5)
30
31 imshow(g)
32
33 title(medfilter gaussian)
34
35 hmedfilt2(P2); %对椒盐躁声中值滤波
36
37 subplot(3,3,6)
38
39 imshow(h)
40
41 title(medfilter salt pepper noise)
42
43 l[1 1 1 %对高斯躁声算术均值滤波
44
45 1 1 1
46
47 1 1 1];
48
49 ll/9;
50
51 kconv2(P1,l);
52
53 subplot(3,3,8)
54
55 imshow(k,[])
56
57 title(arithmeticfilter gaussian)
58
59 %对椒盐躁声算术均值滤波
60
61 dconv2(P2,l);
62
63 subplot(3,3,9)
64
65 imshow(d,[])
66
67 title(arithmeticfilter salt pepper noise)
68
69 4.2参考文献
[1]孟凡文, 吴禄慎.基于FTP的二维傅里叶变换的研究.激光与红外. 第38卷第9期 2008年9月
[2] 董健邓国辉李金武. 基于二维傅里叶变换实现图像变换的研究. 福建电脑. 2015年第 9期
