1、为什么要进行 傅里叶变换,其物理意义是什么?

立叶变换是数字信号处理领域一种很重要的算法。要知道傅立叶变换算法的意义,首先要了解傅立叶原理的意义。傅立叶原理表明:任何连续测量的时序或信号,都

可以表示为不同频率的正弦波信号的无限叠加。而根据该原理创立的傅立叶变换算法利用直接测量到的原始信号,以累加方式来计算该信号中不同正弦波信号的频

率、振幅和相位。

和傅立叶变换算法对应的是反傅立

叶变换算法。该反变换从本质上说也是一种累加处理,这样就可以将单独改变的正弦波信号转换成一个信号。

因此,可以说,傅立叶变换将原来

难以处理的时域信号转换成了易于分析的频域信号(信号的频谱),可以利用一些工具对这些频域信号进行处理、加工。最后还可以利用傅立叶反变换将这些频域信

号转换成时域信号。

从现代数学的眼光来看,傅里叶变

换是一种特殊的积分变换。它能将满足一定条件的某个函数表示成正弦基函数的线性组合或者积分。在不同的研究领域,傅里叶变换具有多种不同的变体形式,如连

续傅里叶变换和离散傅里叶变换。

在数学领域,尽管最初傅立叶分析

是作为热过程的解析分析的工具,但是其思想方法仍然具有典型的还原论和分析主义的特征。"任意"的函数通过一定的分解,都能够表示为正弦函数的线性组合的形式,而正弦函数在物理上是被充分研究而相对简单

的函数类:1. 傅立叶变换是线性算子,若 赋予适当的范数,它还是酉算子;2. 傅立

叶变换的逆变换容易求出,而且形式与正变换非常类似;3. 正弦基函数是微分运算的本征函数,从而使得线性微分方程的求解可以转化为常系数的

代数方程的求解.在线性时不变杂的卷积运算为简单的乘积运算,从 而提供了计算卷积的一种简单手段;5.

离散形式的傅立叶的物理系统内,频率是个不变的性质,从而系统对于复杂激励的响应可以 通过组合其对不同频率正弦信号的响应来获取;4. 著

名的卷积定理指出:傅立叶变换可以化复变换可以利用数字计算机快速的算出(其算法称为快速傅立叶变换算法(FFT))。

正是由于上述的良好性质,傅里叶变换在物理学、数论、组合数学、信号处理、概率、统计、密码学、声学、光学等领域都有着广泛的应用。

2、图像傅立叶变换的物理意义

像的频率是表征图像中灰度变化剧烈程度的指标,是灰度在平面空间上的梯度。如:大面积的沙漠在图像中是一片灰度变化缓慢的区域,对应的频率值很低;而对于

地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域,对应的频率值较高。傅立叶变换在实际中有非常明显的物理意义,设f是一个能量有限的模拟信号,则其傅立叶变换就表示f的

谱。从纯粹的数学意义上看,傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看,傅立叶变换是将图像从空间域转换到频率域,其逆变换是将

图像从频率域转换到空间域。换句话说,傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数,傅立叶逆变换是将图像的频率分布函数变换为

灰度分布函数

傅立叶变换以

前,图像(未压缩的位图)是由对在连续空间(现实空间)上的采样得到一系列点的集合,我们习惯用一个二维矩阵表示空间上各点,则图像可由z=f(x,y)来

表示。由于空间是三维的,图像是二维的,因此空间中物体在另一个维度上的关系就由梯度来表示,这样我们可以通过观察图像得知物体在三维空间中的对应关系。

为什么要提梯度?因为实际上对图像进行二维傅立叶变换得到频谱图,就是图像梯度的分布图,当然频谱图上的各点与图像上各点并不存在一一对应的关系,即使在

不移频的情况下也是没有。傅立叶频谱图上我们看到的明暗不一的亮点,实际上图像上某一点与邻域点差异的强弱,即梯度的大小,也即该点的频率的大小(可以这

么理解,图像中的低频部分指低梯度的点,高频部分相反)。一般来讲,梯度大则该点的亮度强,否则该点亮度弱。这样通过观察傅立叶变换后的频谱图,也叫功率

图,我们首先就可以看出,图像的能量分布,如果频谱图中暗的点数更多,那么实际图像是比较柔和的(因为各点与邻域差异都不大,梯度相对较小),反之,如果

频谱图中亮的点数多,那么实际图像一定是尖锐的,边界分明且边界两边像素差异较大的。对频谱移频到原点以后,可以看出图像的频率分布是以原点为圆心,对称

分布的。将频谱移频到圆心除了可以清晰地看出图像频率分布以外,还有一个好处,它可以分离出有周期性规律的干扰信号,比如正弦干扰,一副带有正弦干扰,移

频到原点的频谱图上可以看出除了中心以外还存在以某一点为中心,对称分布的亮点集合,这个集合就是干扰噪音产生的,这时可以很直观的通过在该位置放置带阻

滤波器消除干扰

另外我还想 说明以下几点:

1、图像经过二维傅立叶变换后,其变换系数矩阵表明:

若变换矩阵Fn原点设在中心,其频谱能量集中分布在变换系数短阵的中心附近(图中阴影区)。若所用的二维傅立叶变换矩阵Fn的原点设在左上角,那么图像信号能量将集中在系数矩阵的四个角上。这是由

二维傅立叶变换本身性质决定的。同时也表明一股图像能量集中低频区域。

2 、变换之后的图像在原点平移之前四角是低频,最亮,平移之后中间部分是

低频,最亮,亮度大说明低频的能量大(幅角比较大)

图像傅里叶变换及相位谱的反变换

这几天在做数字图像处理的作业,一个挺简单的题目:将一幅图像变换成傅里叶变换,显示幅度谱和相位谱,再利用相位谱进行反变换。

正变换倒是挺快的,无非就是fft2,fftshift,显示的时候也好办,运用了im2uint8(mat2gray(log(1+double(f_magnitude1)))),这是为显示而采用的对数变换,mat2gray将值限定在范围[0,1]内,im2uint8将值限定在范围[0,255]内。基本完工。

但是反变换一时把我给难到了,忘记研究公式了,于是直接把

angle()后的值拿来ifft2,结果可想而知了,与变换前一相的黑白相间的块。呜呜……把书拿来死命研究了下,才发现自己好笨哦。只要这样:

f_phase=angle(image_fft_shift);%相位谱

f_phase2=i*f_phase;

f_phase4=1.*(exp(f_phase2));

image_ifft=real(ifft2(ifftshift(f_phase4)));

OK了,完全按公式,看来太久没写程序了。利用相位谱得到的图像就是原图像的结构纹理图呢。

python傅里叶变换相位谱图_[转载]图像傅立叶变换及相位谱相关推荐

  1. 图像傅立叶变换的原理和物理意义

    图像傅立叶变换的原理和物理意义 图像的傅立叶变换,原始图像由N行N列构成,N必须是基2的,把这个N*N个包含图像的点称为实部,另外还需要N*N个点称为虚部,因为FFT是基于复数的,如下图所示: (// ...

  2. python 二维强度图_荐 python数据分析matplotlib库使用之二维图形绘制

    本篇内容会在后期不定时更新 什么是matplotlib matplotlib是最流行的python底层绘图库,主要做数据可视化图表. 为什么要学习matplotlib 能将数据进行可视化,更直观的呈现 ...

  3. python课程思维导图_零基础Python学习思维导图,记得收藏

    原标题:零基础Python学习思维导图,记得收藏 本文主要涵盖了 Python 编程的核心知识(暂不包括标准库及第三方库). 1.按顺序依次展示了以下内容的一系列思维导图:基础知识,数据类型(数字,字 ...

  4. python 网页樱花动态图_如何用Python实现动态图?

    GIF(Graphics Interchange Format,图形交换格式)是一种位图图像格式, GIF格式的图像文件具有如下特点: (1)GIF格式图像文件的扩展名是".gif" ...

  5. python怎么变成动图_教你用 Python 生成 GIF 动图 !

    最近啊 ,微信订阅号改变频繁 ,很多读者后台说 :小詹啊 ,我总是容易错过你公号的消息 ,现在没有置顶功能很难过啊 ! 不止你们难过 ,订阅号的作者恐怕更难过 !现在人人公众号时代 ,大家关注的公众号 ...

  6. sin傅里叶变换公式_全面解析傅立叶变换(非常详细)

    前言 第一部分. DFT 第一章.傅立叶变换的由来 第二章.实数形式离散傅立叶变换(Real DFT) 从头到尾彻底理解傅里叶变换算法.下 第三章.复数 第四章.复数形式离散傅立叶变换 前言:&quo ...

  7. 语谱图(四) Mel spectrogram 梅尔语谱图

    1. 语谱图的产生 前文介绍了语谱图的产生 这里在简单小结如下, 信号预处理,预加重, 分帧, 加窗, 进行STFT 变换, 生成频谱图: 对单帧信号的频谱进行映射, (2,3 步骤如下图) 将映射后 ...

  8. 图像傅立叶变换的物理意义

    原文:http://blog.csdn.net/dadaadao/article/details/6093882 傅立叶变换可以看做拉普拉斯变换的特殊形式.拉氏变换就是将原时域函数乘上一个与 σ相关的 ...

  9. python程序画漂亮图_用python画图代码:正弦图像、多轴图等案例

    python画图 正弦图像: #coding:utf-8 import numpy as np import matplotlib.pyplot as plt x=np.linspace(0,10,1 ...

  10. python画车辆轨迹图_如何利用 Python 绘制酷炫的 车辆轨迹 — 速度时空图?三维数据用二维图像呈现...

    说明:本文系交通攻城狮原创文章,如需转载请私信联系,侵权必究. 2020,第 30 期,编程笔记 建议直接阅读精编版:如何利用 Python 绘制酷炫的 车辆轨迹 - 速度时空图?三维数据用二维图像呈 ...

最新文章

  1. 语言全局变量跨文件_C 语言课上不会讲的一些问题
  2. 跟我一起写 Makefile(三)
  3. 视觉平衡与物理平衡_设计中的视觉平衡
  4. 互联网时代下,看租赁电商如何玩转消费金融
  5. 密歇根安娜堡大学的计算机科学教授,美国密歇根大学安娜堡分校读机械硕士在美国好就业吗?...
  6. 最新xampp下安装mod_security,mod_evasive,mod-cband笔记
  7. php设置路径别名,react设置文件路径别名的具体方法你知道么
  8. 网易回应“公开叫卖网易邮箱账号”:不涉及用户敏感信息
  9. linux是否有安装java_Linux 安装 Java
  10. OpenCV2 图像处理与计算机视觉(一)—— 去除一幅二值化图像中的椒盐噪声
  11. 2020-09-02 微机原理与接口课程复习题
  12. c语言数据结构用矩阵存储图,数据结构之---C语言实现图的数组(邻接矩阵)存储表示...
  13. 面试中的小知识(网络 以及 测试)
  14. hcna华为认证网络工程师
  15. 马斯克喊话库克:昨天你对我爱答不理,今天我让你高攀不起
  16. VS2017设置透明主题
  17. 估算活动持续时间-规划管理
  18. Vue 状态管理与与SSR详解
  19. hand crafted feature:histogram(直方图)
  20. [高数][高昆轮][高等数学上][第一章-函数与极限]03.函数的极限

热门文章

  1. java数字转大写_java实现数字转大写的方法
  2. 高斯过程回归python_高斯过程回归在pythony中的实现(n个样本,n个目标)
  3. 手机浏览器的扫一扫功能在哪里,有什么作用
  4. Ubuntu 18.04 安装驱动
  5. 树莓派远程监控+运动目标检测
  6. retinex 的水下图像增强算法_Retinex图像增强算法
  7. SVN 版本回退 命令行
  8. opensips(1)——安装opensips详细流程
  9. AcrelEMS-BP生物制药工厂能效管理系统
  10. chm打开秒退_无法打开chm文件