简单的求取下灰度图像的幅度谱和相位谱并进行双谱重构:

直接上代码:

clear all

Picture = imread('E:\others\Picture\Library.jpg');

Picture_Gray = rgb2gray(Picture);%灰度处理

Picture_FFT = fft2(Picture_Gray);%傅里叶变换

Picture_FFT_Shift = fftshift(Picture_FFT);%对频谱进行移动,是0频率点在中心

Picture_AM_Spectrum = log(abs(Picture_FFT_Shift));%获得傅里叶变换的幅度谱

Picture_Phase_Specture = log(angle(Picture_FFT_Shift)*180/pi);%获得傅里叶变换的相位谱

Picture_Restructure = ifft2(abs(Picture_FFT).*exp(j*(angle(Picture_FFT))));%双谱重构

figure(1)

subplot(221)

imshow(Picture_Gray)

title('原图像')

subplot(222)

imshow(Picture_AM_Spectrum,[])%显示图像的幅度谱,参数'[]'是为了将其值线性拉伸

title('图像幅度谱')

subplot(223)

imshow(Picture_Phase_Specture,[]);

title('图像相位谱')

subplot(224)

imshow(Picture_Restructure,[]);

title('双谱重构图')

图像的幅度谱代表的是图像各像素点的亮度信息,即该像素应该显示什么颜色,但是做出来的幅度谱却不知道每一点在原图像中具体是哪一点,即幅度谱虽然存储了各个像素点的幅值信息,但是原像素点的位置已经被打乱,所以仅凭幅度谱是没有办法重构原图像的。幅度谱的中心是低频部分,越亮的地方代表的幅度越大。幅度谱中“十”字形亮线表示原图像中水平和垂直方向的分量较其他方向要多,因为在人们周围的自然场景中水平和垂直的线条出现的可能性较大。

来看下仅有幅度谱重构出来的原图像:

Picture_Restructure = ifft2(abs(Picture_FFT));%幅度谱重构

可以看到仅有幅度谱重构出来的图像啥也不是,,,

而相位谱记录的是所有点的相位信息,看起来相位谱是一团噪声,这也说明相位信息是以一种更为隐蔽的方式出现在人们面前的,但它非常重要,因为相位信息中携带者图像的位置信息,没有它将无法从品频谱还原出原图像。

相位谱重构原图像:

Picture_Restructure = ifft2(exp(j*(angle(Picture_FFT))));%相位谱重构

可以看到仅有相位谱重构出来的图像还是可以看到一些轮廓信息的,利用相位谱记录的位置信息和幅度谱记录的亮度信息,就可以用双谱重构的方法恢复出原图像。

————————————————

版权声明:本文为CSDN博主「逸凌123」的原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接及本声明。

原文链接:https://blog.csdn.net/qq_36554582/article/details/88701865

python傅里叶变换相位谱图_图像傅里叶变换的幅度谱、相位谱以及双谱重构原图像...相关推荐

  1. 图像傅里叶变换的幅度谱、相位谱以及双谱重构原图像

    简单的求取下灰度图像的幅度谱和相位谱并进行双谱重构: 直接上代码: clear all Picture = imread('E:\others\Picture\Library.jpg');Pictur ...

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

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

  3. python opencv屏幕找图_使用Python+OpenCV进行图像模板匹配(Match Template)实例-找到百度首页按钮并点击...

    意图:准备一张小图,在电脑屏幕上找到小图坐标,并点击. 1  安装 opencv 和 numpy: pip3 install opencv-python 上述命令将 opencv 和 numpy都安装 ...

  4. python如何做四象图_(Python)从零开始,简单快速学机器仿人视觉Opencv—第十四节:图像梯度...

    事先说明 图像梯度原理:简单来说就是求导 OpenCV提供了三种不同的梯度滤波器,或者说高通滤波器:Sobel,Scharr和Laplacian.Sobel和Scharr是求一阶或二阶导数.Schar ...

  5. python画三维折线图_使用Python的matplotlib画折线图,柱状图,三维图

    因为在各种场景下需要各种实验数据的对比图像,有的中还要求dpi,这些在Python中的matplotlib中都可以实现,下面是总结的各种画图命令. 打包文件:https://download.csdn ...

  6. python绘制三维矢量图_【知乎专栏】让你的论文图表更sexy

    原标题:[知乎专栏]让你的论文图表更sexy 有些同学在写论文的时候光注意文字的阐述,忽略掉了一个重要因素--图表.插入的图表不仅无意义,有时候甚至造成图表影响了文章内容表达.下面干货满满,十分钟教你 ...

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

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

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

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

  9. 图卷积神经网络的数学原理——谱图理论和傅里叶变换初探

    Math of GCN 一.warmup 1.Graph vs Image ①Graph是相比于Image更加广义的一种拓扑结构. ②Image是Grape在欧式空间的一种特例. 2.符号含义 ①A: ...

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

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

最新文章

  1. 时间:2014年3月27日文件和目录操作函数
  2. ArcGIS桌面基本操作 和 ArcSDE数据库基本概念
  3. 记录一下LabelImg实际绘制框颜色的位置
  4. 学习笔记11-C语言-指针
  5. java的向下转型_终于搞明白向下转型的作用了,还不懂的进来看下.
  6. console.log()不显示结果_提醒低端电子显示屏易致视疲劳,OLED屏幕表现略好
  7. OpenCV+Python实现视频文件裁剪功能
  8. 人生苦短,喝点python性能鸡汤
  9. 2017 ACM/ICPC(北京)总结
  10. linux时间同步和登陆数据库shell脚本
  11. 基于SSM的ERP管理系统和仓库管理系统
  12. ps剪贴蒙版教程(ps创建剪贴蒙版步骤)
  13. windows 10屏蔽系统自动更新
  14. rna聚类分析_新技术助力单细胞RNA测序数据聚类分析
  15. 飞机加油游戏 --- 我的易术
  16. 年轻人转行需三思,未来5-10年什么工作真吃香?
  17. Java分布式锁那点事
  18. 【小飞象出海】Facebook海外企业广告账户如何申请基本流程
  19. WindowsTool
  20. C语言二叉树叶子节点的求法

热门文章

  1. jmeter录制脚本后请求太多_使用Jmeter录制web脚本
  2. 部队计算机操作使用教案,计算机基础教案2(键盘鼠标操作).doc
  3. 软件项目管理:使用PERT评价不确定性的方法
  4. 企业工资管理系统论文
  5. 二级路由dhcp关闭连不上wifi_如何做到让家里WiFi真正全覆盖的几个布线方案
  6. 态路小课堂丨光缆知识
  7. 用计算机运行搜索Ip的方法,怎样查ip地址 几种查ip地址的方法【图文】
  8. php 连接 mysql 8.0
  9. 笔记本电脑 联想 Thinkpad E420 无法打开摄像头怎么办
  10. 2020最新Javaweb视频教程-Javaweb从入门到精通【JSP】