python 图像的离散傅立叶变换实例


Posted in Python onJanuary 02, 2020

图像(MxN)的二维离散傅立叶变换可以将图像由空间域变换到频域中去,空间域中用x,y来表示空间坐标,频域由u,v来表示频率,二维离散傅立叶变换的公式如下:

python 图像的离散傅立叶变换实例

在python中,numpy库的fft模块有实现好了的二维离散傅立叶变换函数,函数是fft2,输入一张灰度图,输出经过二维离散傅立叶变换后的结果,但是具体实现并不是直接用上述公式,而是用快速傅立叶变换。结果需要通过使用abs求绝对值才可以进行可视化,但是视觉效果并不理想,因为傅立叶频谱范围很大,所以要用log对数变换来改善视觉效果。

在使用log函数的时候,要写成log(1 + x) 而不是直接用log(x),这是为了避开对0做对数处理。

另外,图像变换的原点需要移动到频域矩形的中心,所以要对fft2的结果使用fftshift函数。最后也可以使用log来改善可视化效果。

代码如下:

import numpy as np
import matplotlib.pyplot as plt

img = plt.imread('photo.jpg')

#根据公式转成灰度图
img = 0.2126 * img[:,:,0] + 0.7152 * img[:,:,1] + 0.0722 * img[:,:,2]

#显示原图
plt.subplot(231),plt.imshow(img,'gray'),plt.title('original')

#进行傅立叶变换,并显示结果
fft2 = np.fft.fft2(img)
plt.subplot(232),plt.imshow(np.abs(fft2),'gray'),plt.title('fft2')

#将图像变换的原点移动到频域矩形的中心,并显示效果
shift2center = np.fft.fftshift(fft2)
plt.subplot(233),plt.imshow(np.abs(shift2center),'gray'),plt.title('shift2center')

#对傅立叶变换的结果进行对数变换,并显示效果
log_fft2 = np.log(1 + np.abs(fft2))
plt.subplot(235),plt.imshow(log_fft2,'gray'),plt.title('log_fft2')

#对中心化后的结果进行对数变换,并显示结果
log_shift2center = np.log(1 + np.abs(shift2center))
plt.subplot(236),plt.imshow(log_shift2center,'gray'),plt.title('log_shift2center')

运行结果:

python 图像的离散傅立叶变换实例

python 图像的离散傅立叶变换实例

根据公式实现的二维离散傅立叶变换如下:

import numpy as np
import matplotlib.pyplot as plt
PI = 3.141591265
img = plt.imread('temp.jpg')

#根据公式转成灰度图
img = 0.2126 * img[:,:,0] + 0.7152 * img[:,:,1] + 0.0722 * img[:,:,2]

#显示原图
plt.subplot(131),plt.imshow(img,'gray'),plt.title('original')

#进行傅立叶变换,并显示结果
fft2 = np.fft.fft2(img)
log_fft2 = np.log(1 + np.abs(fft2))
plt.subplot(132),plt.imshow(log_fft2,'gray'),plt.title('log_fft2')

h , w = img.shape
#生成一个同样大小的复数矩阵
F = np.zeros([h,w],'complex128')
for u in range(h):
 for v in range(w):
  res = 0
  for x in range(h):
   for y in range(w):
    res += img[x,y] * np.exp(-1.j * 2 * PI * (u * x / h + v * y / w))
  F[u,v] = res
log_F = np.log(1 + np.abs(F))
plt.subplot(133),plt.imshow(log_F,'gray'),plt.title('log_F')

直接根据公式实现复杂度很高,因为是四重循环,时间复杂度为python 图像的离散傅立叶变换实例 ,所以实际用的时候需要用快速傅立叶变换来实现

以上这篇python 图像的离散傅立叶变换实例就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python检测某个变量是否有定义的方法
May 20 Python
Python脚本实现自动发带图的微博
Apr 27 Python
python如何在终端里面显示一张图片
Aug 17 Python
Python正则表达式教程之三:贪婪/非贪婪特性
Mar 02 Python
matplotlib在python上绘制3D散点图实例详解
Dec 09 Python
详解python selenium 爬取网易云音乐歌单名
Mar 28 Python
python爬虫简单的添加代理进行访问的实现代码
Apr 04 Python
Python3 批量扫描端口的例子
Jul 25 Python
python多线程同步之文件读写控制
Feb 25 Python
python中@property和property函数常见使用方法示例
Oct 21 Python
Python实现把多维数组展开成DataFrame
Nov 30 Python
最简单的matplotlib安装教程(小白)
Jul 28 Python
Python加密模块的hashlib,hmac模块使用解析
Jan 02 #Python
在win64上使用bypy进行百度网盘文件上传功能
Jan 02 #Python
pytorch实现onehot编码转为普通label标签
Jan 02 #Python
pytorch标签转onehot形式实例
Jan 02 #Python
Python socket聊天脚本代码实例
Jan 02 #Python
解决Pytorch训练过程中loss不下降的问题
Jan 02 #Python
Pytorch evaluation每次运行结果不同的解决
Jan 02 #Python
You might like
PHP 已经成熟
2006/12/04 PHP
PHP 图片文件上传实现代码
2010/12/29 PHP
php实现基于pdo的事务处理方法示例
2017/07/21 PHP
laravel 如何实现引入自己的函数或类库
2019/10/15 PHP
Aster vs KG BO3 第三场2.18
2021/03/10 DOTA
javascript的document.referrer浏览器支持、失效情况总结
2014/07/18 Javascript
javascript在当前窗口关闭前检测窗口是否关闭
2014/09/29 Javascript
jQuery层级选择器用法分析
2015/02/10 Javascript
浅析javascript异步执行函数导致的变量变化问题解决思路
2016/05/13 Javascript
Active控件问题小结(附解决办法)
2016/06/09 Javascript
js实现数组去重方法及效率?Ρ? target=
2017/02/14 Javascript
常用的js方法合集
2017/03/10 Javascript
使用selenium抓取淘宝的商品信息实例
2018/02/06 Javascript
如何编写一个d.ts文件的步骤详解
2018/04/13 Javascript
JS原生瀑布流效果实现
2019/04/26 Javascript
用Angular实现一个扫雷的游戏示例
2020/05/15 Javascript
使用Vue实现一个树组件的示例
2020/11/06 Javascript
浅谈Python中的数据类型
2015/05/05 Python
python实现合并两个数组的方法
2015/05/16 Python
Python中的字典与成员运算符初步探究
2015/10/13 Python
用不到50行的Python代码构建最小的区块链
2017/11/16 Python
Django model select的多种用法详解
2019/07/16 Python
Python实现京东抢秒杀功能
2021/01/25 Python
中国领先的专业演出票务网:永乐票务
2016/08/29 全球购物
Harrods英国:世界领先的奢侈品百货商店
2020/09/23 全球购物
PHP数据运算类型都有哪些
2013/11/05 面试题
超市5.1促销活动
2014/01/15 职场文书
社区健康教育实施方案
2014/03/18 职场文书
村干部培训方案
2014/05/02 职场文书
青安岗事迹材料
2014/05/14 职场文书
员工安全责任书范本
2014/07/24 职场文书
企业领导对照检查材料
2014/08/20 职场文书
精神文明建设汇报材料
2014/12/24 职场文书
员工辞职信怎么写
2015/02/27 职场文书
幼儿园元旦主持词
2015/07/06 职场文书
Python编程根据字典列表相同键的值进行合并
2021/10/05 Python