以下部分文字资料整合于网络，本文仅供自己学习用！

这是一幅很绝的一维傅里叶变换动态图

一，读懂傅里叶变换

一个信号能表示成傅里叶级数的形式是有条件的，首先它必须是周期信号，第二必须是满足狄里赫利条件的周期信号。

1，关于两大域：时域与频域

（1），频域（frequency domain）是指在对函数或信号进行分析时，

分析其和频率有关部份，而不是和时间有关的部份，和时域一词相对。

（2），时域是描述数学函数或物理信号对时间的关系。

例如一个信号的时域波形可以表达信号随着时间的变化。
若考虑离散时间，时域中的函数或信号，在各个离散时间点的数值均为已知。
若考虑连续时间，则函数或信号在任意时间的数值均为已知。在研究时域的信号时，常会用示波器将信号转换为其时域的波形。

（3），两者相互间的变换
时域（信号对时间的函数）和频域（信号对频率的函数）的变换在数学上是通过积分变换实现。
对周期信号可以直接使用傅立叶变换，对非周期信号则要进行周期扩展，使用拉普拉斯变换。

2，傅立叶变换

（1），什么是傅里叶变换？
也称作傅立叶变换，表示能将满足一定条件的某个函数表示成三角函数（正弦和/或余弦函数）或者它们的积分的线性组合。在不同的研究领域，傅立叶变换具有多种不同的变体形式，如连续傅立叶变换和离散傅立叶变换。最初傅立叶分析是作为热过程的解析分析的工具被提出的。傅里叶变换是一种分析信号的方法，它可分析信号的成分，也可用这些成分合成信号。许多波形可作为信号的成分，比如正弦波、方波、锯齿波等，傅里叶变换用正弦波作为信号的成分。
傅里叶变换的实质是将一个信号分离为无穷多多正弦/复指数信号的加成，也就是说，把信号变成正弦信号相加的形式——既然是无穷多个信号相加，那对于非周期信号来说，每个信号的加权应该都是零——但有密度上的差别，你可以对比概率论中的概率密度来思考一下——落到每一个点的概率都是无限小，但这些无限小是有差别的所以，傅里叶变换之后，横坐标即为分离出的正弦信号的频率，纵坐标对应的是加权密度
（2），傅里叶变换有什么用呢？
举例说明：傅里叶变换可以将一个时域信号转换成在不同频率下对应的振幅及相位，其频谱就是时域信号在频域下的表现，而反傅里叶变换可以将频谱再转换回时域的信号。最简单最直接的应用就是时频域转换，比如在移动通信的LTE系统中，要把接收的信号从时域变成频域，就需要使用FFT。又例如对一个采集到的声音做傅立叶变化就能分出好几个频率的信号。比如南非世界杯时，南非人吹的呜呜主拉的声音太吵了，那么对现场的音频做傅立叶变化（当然是对声音的数据做），会得到一个展开式，然后找出呜呜主拉的特征频率，去掉展开式中的那个频率的sin函数，再还原数据，就得到了没有呜呜主拉的嗡嗡声的现场声音。而对图片的数据做傅立叶，然后增大高频信号的系数就可以提高图像的对比度。同样，相机自动对焦就是通过找图像的高频分量最大的时候，就是对好了。

3，快速傅里叶变换

（1），什么是快速傅里叶变换？
计算离散傅里叶变换的一种快速算法，简称FFT。快速傅里叶变换是1965年由J.W.库利和T.W.图基提出的。采用这种算法能使计算机计算离散傅里叶变换所需要的乘法次数大为减少，特别是被变换的抽样点数N越多，FFT算法计算量的节省就越显著。

（2），快速傅里叶变换有什么用？
函数或信号可以透过一对数学的运算子在时域及频域之间转换。和傅里叶变换作用一样！

（3），为什么要提出快速傅里叶变换？
人们想让计算机能处理信号 但由于信号都是连续的、无限的，计算机不能处理，于是就有了傅里叶级数、傅里叶变换，将信号由时域变到频域，把一个信号变为有很多个不同频率不同幅度的正弦信号组成，这样计算机就能处理了，但又由于傅里叶变换中要用到卷积计算，计算量很大，计算机也算不过来，于是就有了快速傅里叶变换，大大降低了运算量，使得让计算机处理信号成为可能。快速傅里叶变换是傅里叶变换的快速算法而已，主要是能减少运算量和存储开销，对于硬件实现特别有利。

4，图像中傅立叶变换的物理意义

（1），图像的频率是表征图像中灰度变化剧烈程度的指标，是灰度在平面空间上的梯度。如：大面积的沙漠在图像中是一片灰度变化缓慢的区域，对应的频率值很低；而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域，对应的频率值较高。傅立叶变换在实际中有非常明显的物理意义，设f是一个能量有限的模拟信号，则其傅立叶变换就表示f的谱。从纯粹的数学意义上看，傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看，傅立叶变换是将图像从空间域转换到频率域，其逆变换是将图像从频率域转换到空间域。换句话说，傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数，傅立叶逆变换是将图像的频率分布函数变换为灰度分布函数

（2），傅立叶变换以前，图像（未压缩的位图）是由对在连续空间（现实空间）上的采样得到一系列点的集合，我们习惯用一个二维矩阵表示空间上各点，则图像可由z=f(x,y)来表示。由于空间是三维的，图像是二维的，因此空间中物体在另一个维度上的关系就由梯度来表示，这样我们可以通过观察图像得知物体在三维空间中的对应关系。为什么要提梯度？因为实际上对图像进行二维傅立叶变换得到频谱图，就是图像梯度的分布图，当然频谱图上的各点与图像上各点并不存在一一对应的关系，即使在不移频的情况下也是没有。傅立叶频谱图上我们看到的明暗不一的亮点，实际上图像上某一点与邻域点差异的强弱，即梯度的大小，也即该点的频率的大小（可以这么理解，图像中的低频部分指低梯度的点，高频部分相反）。一般来讲，梯度大则该点的亮度强，否则该点亮度弱。这样通过观察傅立叶变换后的频谱图，也叫功率图，我们首先就可以看出，图像的能量分布，如果频谱图中暗的点数更多，那么实际图像是比较柔和的（因为各点与邻域差异都不大，梯度相对较小），反之，如果频谱图中亮的点数多，那么实际图像一定是尖锐的，边界分明且边界两边像素差异较大的。对频谱移频到原点以后，可以看出图像的频率分布是以原点为圆心，对称分布的。将频谱移频到圆心除了可以清晰地看出图像频率分布以外，还有一个好处，它可以分离出有周期性规律的干扰信号，比如正弦干扰，一副带有正弦干扰，移频到原点的频谱图上可以看出除了中心以外还存在以某一点为中心，对称分布的亮点集合，这个集合就是干扰噪音产生的，这时可以很直观的通过在该位置放置带阻滤波器消除干扰

二，理性认识傅里叶变换

1，二维傅里叶fft2对简单矩阵的操作

（可以把s想象成一幅图像，傅里叶变换过程如下）

[cpp] view plain copy
  1. s=magic(2);%magic矩阵有一个神奇的特点：每行或者每列加起来都相等
  2. f=fft2(s);
  3. a=abs(f)
  4. a =
  5. 10     0
  6. 2     4
  7. >> f
  8. f =
  9. 10     0
  10. -2    -4
  11. >> s
  12. s =
  13. 1     3
  14. 4     2
  15. %3*3的矩阵的操作结果
  16. s=magic(3);
  17. f=fft2(s);
  18. a=abs(f)
  19. a =
  20. 45.0000         0         0
  21. 0.0000   15.5885    5.1962
  22. 0.0000    5.1962   15.5885
  23. >> f
  24. f =
  25. 45.0000                  0                  0
  26. 0 + 0.0000i  13.5000 + 7.7942i   0.0000 - 5.1962i
  27. 0 - 0.0000i   0.0000 + 5.1962i  13.5000 - 7.7942i
  28. >> s
  29. s =
  30. 8     1     6
  31. 3     5     7
  32. 4     9     2
  33. %4*4的矩阵操作结果
  34. s=magic(4);%magic矩阵有一个神奇的特点：每行或者每列加起来都相等
  35. f=fft2(s);
  36. a=abs(f)
  37. a =
  38. 136.0000         0         0         0
  39. 0   20.0000   11.3137   12.0000
  40. 0   45.2548         0   45.2548
  41. 0   12.0000   11.3137   20.0000
  42. >> f
  43. f =
  44. 1.0e+02 *
  45. 1.3600                  0                  0                  0
  46. 0             0.2000             0.0800 + 0.0800i        0 - 0.1200i
  47. 0             0.3200 + 0.3200i        0             0.3200 - 0.3200i
  48. 0                  0 + 0.1200i   0.0800 - 0.0800i   0.2000
  49. >> s
  50. s =
  51. 16     2     3    13
  52. 5    11    10     8
  53. 9     7     6    12
  54. 4    14    15     1

2，二维快速傅里叶fft2对灰度图的操作

[cpp] view plain copy
  1. <span style="font-size:12px;">clear all;
  2. close all;
  3. I=imread('peppers.png');%读入图像
  4. J=rgb2gray(I);%将图像转换为灰度图
  5. K=fft2(J);%对图像进行二维快速傅里叶变换
  6. K=fftshift(K);%将频谱转移到中心，其实就是在傅里叶变换时乘以了某个因子
  7. L=abs(K/256);%取模
  8. %显示图片
  9. figure;
  10. subplot(131);
  11. imshow(I);title('原图像')
  12. subplot(132);
  13. imshow(J);title('被转换为灰度图后')
  14. subplot(133);
  15. imshow(uint8(L));title('二维傅里叶变换后的频谱图')</span>

[cpp] view plain copy
  1. clear all; close all;
  2. I=imread('peppers.png');%读入图像
  3. J=rgb2gray(I);%转换为灰度图
  4. J=imrotate(J,45,'bilinear');    %将图像旋转45度角
  5. K=fft2(J);%对图像进行二维傅里叶比变换
  6. K=fftshift(K);%转移频谱中心
  7. L=abs(K/256);%取模
  8. figure;
  9. subplot(121);
  10. imshow(J);title('原图像被旋转45度角')
  11. subplot(122);
  12. imshow(uint8(L));title('傅里叶变换后的图像')

[cpp] view plain copy
  1. clear all; close all;
  2. I=imread('peppers.png');
  3. J=rgb2gray(I);
  4. J=imnoise(J, 'gaussian', 0, 0.01);%加入高斯噪声
  5. K=fft2(J);
  6. K=fftshift(K);
  7. L=abs(K/256);
  8. figure;
  9. subplot(131);
  10. imshow(I);title('原图像')
  11. subplot(132);
  12. imshow(J);title('高斯噪声污染')
  13. subplot(133);
  14. imshow(uint8(L));title('污染后进行福利叶变换')

[cpp] view plain copy
  1. clear all; close all;
  2. I=imread('onion.png');
  3. J=rgb2gray(I);
  4. K=fft2(J);%傅里叶变换
  5. L=fftshift(K);%转移频谱中心
  6. L=abs(L/256);%取模
  7. M=ifft2(K);%二维傅里叶反变换
  8. figure;
  9. subplot(131);
  10. imshow(J);title('原灰度图像')
  11. subplot(132);
  12. imshow(uint8(L));title('傅里叶变换后的图像')
  13. subplot(133);
  14. imshow(uint8(M));title('反变换后的图像')

参考资源：

【1】http://blog.csdn.net/struggle_for_m/article/details/51207370

【2】动态图来源于“维基百科”

【3】https://zhuanlan.zhihu.com/p/19763358?columnSlug=wille

本文转自：本文为EbowTang原创文章

数字图像处理，读懂频域处理的“傅里叶变换”相关推荐

  1. Matlab数字图像处理 实验1、数字图像的傅里叶变换

    子曰:桃李不言,下自成蹊. 数字图像的傅里叶变换数字图像的傅里叶变换数字图像的傅里叶变换 Part.0实验目的仅限于学习交流Part.0 实验目的^{仅限于学习交流}Part.0实验目的仅限于学习交流 ...

  2. 数字图像处理及图像频域

    数字图像处理方法的重要性源于两个主要的应用领域:改善图示信息以便人们解释:为存储.传输和表示而对图像数据进行处理,以便于机器自动理解. 一幅数字图像可定义为一个二维函数f(x,y),其中x和y是空间( ...

  3. 【MATLAB】数字图像处理--图像的频域处理-相位谱重构

    效果 %% clc clear alldir = 'image/56.jpg';img = imread(dir); imggray = rgb2gray(img);%灰度处理 imgf = fft2 ...

  4. 数字图像处理：python对图像做傅里叶变换，理想低通滤波器，理想高通滤波器

    对Lena图像进行以下处理(可使用编程语言自带函数,以实验报告形式提交word文件): (1)做傅里叶变换,并显示频谱图像: (2)使用理想低通滤波器,截止频率半径分别选取10.20.50.80,分别 ...

  5. 沃尔什哈达玛变换Matlab,哈达玛变换矩阵-数字图像处理.ppt

    哈达玛变换矩阵-数字图像处理 3.1 二维离散傅里叶变换(DFT) 3.1.1 二维连续傅里叶变换 二维连续函数 f (x, y)的傅里叶变换定义如下: 设 是独立变量 的函数,且在 上绝对可积,则定 ...

  6. JAVA公路车几何图_单车基械匠：读懂自行车几何角度图，是你成为老鸟的关键一步...

    大家好,欢迎大家收看新一期的单车基械匠. 学会阅读自行车的几何角度图,可以帮助你了解你大致了解你要购买的车对操控的处理方式,骑行感觉和舒适性.当然,尺寸是否合适也是一大关键性难题,而所有这些都可以在车 ...

  7. 数字图像处理报告：实验3 同态滤波、频域滤波、傅里叶变换性质、DCT变换性质

    实验3 同态滤波.频域滤波.傅里叶变换性质.DCT变换性质 一.实验主题: 学习同态滤波.频域滤波.傅里叶变换性质.DCT变换性质,掌握其原理,并学会编程,实现这些方法与性质. 二.实验目的: 1.掌 ...

  8. 基于skimage的数字图像处理——频域滤波（傅里叶变换）

    注意:文章代码基于notebook实现,也可修改后用于pycharm等编译器. 文章目录 前言 一.频域滤波与傅里叶变换 二.skimage代码实现 1.简单的傅里叶变换: 2.不同种类的正弦波及其傅 ...

  9. 10句话读懂图像频域滤波——不能不知道的信号与系统基本理论

    10句话读懂图像频域滤波--不能不知道的信号与系统基本理论 今天的图像处理依靠各类方便易用的工具箱与函数库似乎已然成为上手就能用的应用科学.但没有那种算法是普适的,知其所以然才能真正理解原理,深刻的理 ...

最新文章

  1. Poj2480欧拉函数
  2. 安装数据库windows 安装redmine 详解
  3. 【C/C++】递归算法
  4. 信恒支付源码-第四方支付源码
  5. 【今日CV 计算机视觉论文速览 第106期】Fri, 26 Apr 2019
  6. RESTful API标准设计教程
  7. OSX Finder显示隐藏文件
  8. nekohtml和xpath
  9. 三星手机刷机后显示无服务器,三星Odin刷机工具连接不上_Odin工具识别不了手机的问题...
  10. 徐培成大数据第一季java基础-徐培成-专题视频课程
  11. golang安装配置
  12. 文档计算机无法分页,word文档总是重新分页 word文档老是分页显示怎么解决
  13. Program Slicing
  14. WORD图片无法同时选中/WORD图片无法组合
  15. HZOI20190819模拟26题解
  16. Flutter 使用自定义 fluro 路由转场动画实现个性化页面切换
  17. jQuery - fadeIn淡入效果
  18. 【EasyX图形库】c语言贪吃蛇-音乐-图片-界面比较美观
  19. 信息的可再生性举例_举例说说信息的可处理性
  20. 几个有用的谷歌Chrome插件

热门文章

  1. 32MB块意味着比特币现金为大量采用做好准备
  2. TypeScript + React 学习render props
  3. ubuntu Gitolite管理git server代码库权限
  4. C中位域结合对齐(pragma)的例子及解释
  5. Dropbox 官方中文版！最优秀实用的免费跨平台文件网络同步网盘云存储服务
  6. MongoDB索引类型
  7. 今日英语：out of the box
  8. Vue Router的详细教程
  9. 小白入坑 Web 渗透测试必备指南
  10. 如何使用KingShard构建MySQL集群