%% 读取图像、预处理
Picture1 = imread('Cameraman.tif');tryPicture1=rgb2gray(Picture1);    % 如果是RGB图像,转成灰度图
catch% 否则,不做变换
endrow = size(Picture1,1);
col= size(Picture1,2);
Picture_Gray1 = Picture1(1:row, 1:col);%% 图像傅里叶变换 及其重构
Picture_FFT1 = fft2(Picture_Gray1); % 傅里叶变换
Picture_FFT_Shift1 = fftshift(Picture_FFT1); % 对频谱进行移动,是0频率点在中心
Picture_AM_Spectrum1 = log(abs(Picture_FFT_Shift1)); % 获得傅里叶变换的幅度谱
Picture_Phase_Specture1 = log(angle(Picture_FFT_Shift1)*180/pi); % 获得傅里叶变换的相位谱Picture_Restructure1 = ifft2(abs(Picture_FFT_Shift1)); % 幅度谱重构
Picture_Restructure2 = ifft2(exp(1i*(angle(Picture_FFT_Shift1)))); % 相位谱重构
Picture_Restructure3 = ifft2(abs(Picture_FFT1).*exp(1i*(angle(Picture_FFT1)))); % 双谱重构
%% 变换结果可视化
figure
subplot(241)
imshow(Picture_Gray1),title('原始图像')
subplot(242)
imshow(Picture_FFT_Shift1,[]),title('傅里叶频谱图')
subplot(243)
imshow(Picture_AM_Spectrum1,[]),title('图像幅度谱')
% 显示图像的幅度谱,参数'[]'是为了将其值线性拉伸
subplot(244)
imshow(Picture_Phase_Specture1,[]),title('图像相位谱')
subplot(245)
imshow(Picture_Restructure1,[]),title('幅度谱重构图')
subplot(246)
imshow(Picture_Restructure2,[]),title('相位谱重构图')
subplot(247)
imshow(Picture_Restructure3,[]),title('双谱重构图')

使用matlab读取图像并通过matlab自带的fft2,fftshift,ifftshift,ifft2等函数获取该图像的傅里叶频谱图,幅度图,相位谱图,幅度图,相位谱图,幅度谱重建图以及相位谱重建相关推荐

  1. matlab 读取照片imread,利用matlab读取图像

    怎样用matlab读取20张图片并依次展示出来 指定路径下 单个文件夹data中所有图像 P = '.\data\';% 图像文件夹路径 img_path_list = dir(strcat(P,'* ...

  2. matlab读取grib三维数据,matlab读取grib2文件

    用matlab读取nc文件是简单的.但气象局的很多数据存储文件是grib2格式的.该如何读取呢?看了网上的众多方法,nctoolbox,read_grib都没有成功.最后投机取巧,用了36计中的瞒天过 ...

  3. MATLAB读取网页内链接,Matlab使用urlread()读取网页乱码问题

    2021年2月9日发(作者:television是什么意思) Matlab 使用 urlread() 读取网页乱码问题 事先声明这是本人在网上找到的( / 2012-05/ ) ,非 原创: Matl ...

  4. matlab读取地理栅格数据,基于MATLAB的时间序列和地理栅格数据分析

    248 摘 要:对MATLAB的时间序列和地理栅格数据智能化处理的特点进行介绍,并结合MATLAB的编程语言重点研究了其在地理栅格数据智能化处理和时间分析等方面的优越性能以及实践应用,能在实践工作中运 ...

  5. matlab读取正则找科学,MATLAB常用正则表达式记录

    MATLAB正则表达式 一.正则表达式 正则表达式是一串用于定义某种模式的字符.在有些情况下(例如,在解析程序输入或处理文本块时),您通常会使用正则表达式在文本中搜索与该模式匹配的一组单词. 数据处理 ...

  6. matlab读取正则找科学,Matlab正则表达式(递归、匹配上下文等)(转载)

    元字符 描述 .点 匹配任何单个字符.例如正则表达式r.t匹配这些字符串:rat.rut.r t,但是不匹配root. $ 匹配行结束符.例如正则表达式weasel$ 能够匹配字符串"He' ...

  7. matlab读取wav文件出错,MATLAB读取wav文件

    %读一个文件夹中的多个wav文件,然后拼接在一起 clc; clear all; folder = 'D:/FFOutput/xiami/'; files = dir([folder '*.wav'] ...

  8. matlab读取txt日期数据格式,matlab中如何读取带有日期和时间的txt数据文件,并放入矩阵中...

    满意答案 g100021221 2018.01.17 采纳率:59%    等级:8 已帮助:312人 假定data_zd.txt文件内容为: 2018-01-16 00:00:00 -290.27 ...

  9. php 读取优酷视频缩略图,PHP优酷土豆酷6采集入库函数(获取视频缩略图,视频swf地址,视频标题)...

    /** * 采集入库函数 * 优酷,土豆,酷6 采集 (自动获取视频缩略图,视频swf地址,视频标题) * by hkshadow * QQ 2765237 * dete: 2011-06-25 AM ...

最新文章

  1. 刘群:华为诺亚方舟NLP预训练模型工作的研究与应用 | AI ProCon 2019
  2. 让Ubuntu自动登录到X
  3. libevent中事件的添加与删除
  4. 如何实现网站文件动静分离
  5. 数据库连接客户端使用(db2,oracle,mysql)
  6. 批处理命令 For循环命令详解!
  7. Going Further with Point Pair Features: PPF, 如何走的更远?
  8. 扩展Kubernetes API
  9. Python基础第二十天:线程
  10. PFQ: a Linux kernel module for packet capturing on multi-core architectures
  11. 计算机新功能,利用win7新功能提升工作效率
  12. 如何获取到电脑所连接Wifi的密码
  13. 当360屠榜黑客奥斯卡,我们为什么要关注国家级网络安全战?
  14. android点击复制链接地址,在Android中的EditView中可点击链接和复制/粘贴菜单
  15. DDD的创新思想:开发即设计思想
  16. 【CF869E】The Untended Antiquity(哈希+二维树状数组)
  17. JS调用拨打电话功能
  18. mysql executing_MySQL数据库出现慢查询的危害
  19. Linux查看文件行数
  20. Android 关于TextView 默认显示英文字母大写问题

热门文章

  1. HUAWEI 机试题:单词接龙
  2. GaussDB(DWS) 数据仓库服务核心技术
  3. 机器人跳跃问题(字节跳动笔试题)
  4. MAC系统 WORD 如何调整自动序号的间隔距离
  5. 标准输入,标准输出和标准错误
  6. css自适应网页(大作业版)
  7. 【推荐系统技术实现】
  8. 糖果枕式包装机总体设计及横封切断装置设计(设计说明书+8张CAD图纸)
  9. 易语言与Python socket通信
  10. c++模型推理时 HWC转CHW