《数字图像处理》实验报告四
一、实验任务与要求
- 对 Fig0403.tif 进行傅里叶变换并显示其频谱图像;fft2(x)
- 对 Fig0405.tif 图像进行填充和非填充的高斯滤波,并观察其不同;paddedsize,fft2(x,m,n)
- 由 sobel 空间滤波算子生成相应的频率域滤波算子并比较各自对 Fig0409.tif 的滤波效果;fspecial( 'sobel' )
- 任意设计一个频域的低通或高通滤波器,并对 Fig0413.tif 进行处理
二、实验报告
(一)对 Fig0403.tif 进行傅里叶变换并显示其频谱图像
1、转换图像类型为 double
f=imread('E:\Ex_4\Fig0403.tif');
f_d=im2double(f);
2、使用 fft2 进行傅里叶变换
f2=fft2(f_d);
f2=fftshift(f2); %将零频点移到频谱中间
f2=abs(f2); %取绝对值
F2=log(abs(f2)+1); %取模并进行缩放
3、效果
(二)对 Fig0405.tif 图像进行填充和非填充的高斯滤波,并观察其不同
1、转换图像类型为 double
f=imread('E:\Ex_4\Fig0405.tif');
f=im2double(f);
subplot(1,3,1)
imshow(f);title('原图');
2、填充的高斯滤波
(1)使用 paddedsize 获得填充参数
PQ=paddedsize(size(f));
(2)得到有填充图像的傅里叶变换
F=fft2(f,PQ(1),PQ(2));
(3)使用 freqz2 生成一个大小为 PQ(1)*PQ(2) 的高斯频率域滤波器
h=fspecial('gaussian',[PQ(1),PQ(2)],5);
H=freqz2(h,PQ(1),PQ(2));
(4)使用 ifftshift 进行反变换
H=ifftshift(H);
(5)滤波器乘以傅里叶变换,并使用 ifft2 进行傅里叶反变换
G=H.*F;
g=ifft2(G);
(6)将左上部的矩形修剪为原始大小,并转换回 uint8 类型
g=g(1:size(f,1),1:size(f,2));
g=im2uint8(g);
subplot(1,3,2)
imshow(g);title('填充');
3、不填充的高斯滤波
[r,c]=size(f);
F2=fft2(f); %无填充的傅里叶变换
h=fspecial('gaussian',[r,c],5);
H2=freqz2(h,r,c);
H2=ifftshift(H2);
G2=H2.*F2;
g2=ifft2(G2);
g2=im2uint8(g2);
subplot(1,3,3)
imshow(g2);title('无填充');
4、效果
可以看出,填充的高斯滤波处理图像四周边缘有都有黑边,而不填充的高斯滤波处理只有上边缘有黑边。
(三)由 sobel 空间滤波算子生成相应的频率域滤波算子并比较各自对 Fig0409.tif 的滤波效果
1、转换图像类型为 double 进行傅里叶变化
f=imread('E:\Ex_4\Fig0409.tif');
subplot(1,3,1)
imshow(f);title('原图');
f_d=im2double(f);
f2=fft2(f_d);
2、 sobel 空间域滤波器
h=fspecial('sobel'); %sobel空间滤波器
f1=imfilter(f,h);
subplot(1,3,2);
imshow(f1);title('空间域滤波算子处理');
3、sobel 频率域滤波器
(1)使用 paddedsize 获得填充参数
PQ=paddedsize(size(f));
(2)得到有填充图像的傅里叶变换
f3=fft2(f_d,PQ(1),PQ(2));
(3)使用 freqz2 转换为频率域滤波器,并进行反变换
H=freqz2(h,PQ(1),PQ(2));
H=ifftshift(H);
(4)滤波器乘以傅里叶变换,并使用 ifft2 进行傅里叶反变换
G=H.*f3;
g=ifft2(G);
g=g(1:size(f,1),1:size(f,2));
g=im2uint8(g);
subplot(1,3,3);
imshow(g);title('sobel频率域滤波算子');
4、效果
可以看出,空间域滤波和频域滤波对图像的滤波效果并无明显区别。
(四)任意设计一个频域的低通或高通滤波器,并对 Fig0413.tif 进行处理
1、对图像进行傅里叶变化并将分量移到频谱中心
f=imread('E:\Ex_4\Fig0413.tif');
subplot(1,2,1);imshow(f);title('原图');
f=im2double(f); %转换为double类型
f2=fft2(f); %傅里叶变换
f2=fftshift(f2); %将变换的原点移到频率矩形的中心
2、设计滤波器的传递函数
[R,C]=size(f2);
nn=2;
d0=50;
r=fix(R/2);
c=fix(C/2);
result=ones(R,C); %创建一个R*C的全一矩阵
for i=1:Rfor j=1:Cd=sqrt((i-r)^2+(j-c)^2);h=1/(1+0.414*(d/d0)^(2*nn));result(i,j)=h*f2(i,j);end
end
result=ifftshift(result);
f3=ifft2(result);
f4=im2uint8(real(f3));
subplot(1,2,2)
imshow(f4);title('频域滤波器处理');
3、处理效果
可以看出,产生了模糊效果
相关文章:

《数字图像处理》实验报告四
一、实验任务与要求 对 Fig0403.tif 进行傅里叶变换并显示其频谱图像;fft2(x) 对 Fig0405.tif 图像进行填充和非填充的高斯滤波,并观察其不同;paddedsize,fft2(x,m,n) 由 sobel 空间滤波算子生成相应的频率…...

算法04 模拟算法之一维数组相关内容详解【C++实现】
大家好,我是bigbigli,模拟算法我们将分为几个章节来讲,今天我们只看一维数组相关的题目 目录 模拟的概念 训练:开关灯 解析 参考代码 训练:数组变化 解析 参考代码 训练:折叠游戏 解析 参考代码 …...

【技术解码】百数SRM:如何助力企业快速优化供应链管理?
SRM应用是企业优化供应链管理的重要工具,它帮助企业全面管理供应商关系,从评估、选择到协同合作和绩效监控,确保供应链的稳定性和效率。 对于企业来说,通过全面管理供应商关系,可以降低采购风险,提升产品质…...

想要用tween实现相机的移动,three.js渲染的canvas画布上相机位置一点没动,如何解决??
🏆本文收录于「Bug调优」专栏,主要记录项目实战过程中的Bug之前因后果及提供真实有效的解决方案,希望能够助你一臂之力,帮你早日登顶实现财富自由🚀;同时,欢迎大家关注&&收藏&&…...

SQL连接与筛选:解析left join on和where的区别及典型案例分析
文章目录 前言数据库在运行时的执行顺序一、left join on和where条件的定义和作用left join on条件where条件 二、left join on和where条件的区别原理不同left join原理:where原理: 应用场景不同执行顺序不同(作用阶段不同)结果集…...
oliva-bruteforce-luks
olivaeasyLUKS v2破解、bruteforce-luks工具使用、cryptsetup使用、cap_dac_read_searcheip、mysql使用 主机发现 ┌──(kali㉿kali)-[~/桌面/OSCP] └─$ sudo netdiscover -i eth0 -r 192.168.44.148/24服务扫描 ┌──(kali㉿kali)-[~/桌面/OSCP] └─$ sudo nmap -sV -…...

图像超分辨率重建
一、什么是图像超分辨 图像超分辨是一种技术,旨在通过硬件或软件的方法提高原有图像的分辨率。这一过程涉及从一系列低分辨率的图像中获取一幅高分辨率的图像,实现了时间分辨率向空间分辨率的转换。超分辨率重建的核心思想是利用多帧图像序列的时间带宽来…...

小米上架遇到的隐私协议问题
1. 找到【APP权限设置】,点击详情,一一对照,删除没用的,新增小米商家必须要有的内容 2. APP 存在未经用户同意读取“OAID”的行为 uniapp官方文档对应内容处...
【区分vue2和vue3下的element UI Message 消息提示组件,分别详细介绍属性,事件,方法如何使用,并举例】
在 Vue 2 中,我们通常使用 Element UI 的 this.$message 方法来显示消息提示,而不是作为一个组件直接在模板中使用。然而,在 Vue 3 的 Element Plus 中,虽然 this.$message 的使用方式仍然保留,但官方文档可能更倾向于…...
架构设计 - Nginx Lua 缓存配置
摘要: web 应用业务缓存通常3级: 一级缓存:JVM 本地缓存 二级缓存:Redis集中式缓存 三级缓存:Nginx Proxy Cache 缓存 或 Nginx Lua 缓存 四级缓存:静态资源CDN缓存 页面静态化 本文主要分享 Nginx Lua 缓存配置开发 鉴于 Nginx Proxy Cache 缓存的劣势,在生产项目…...
lua的GC
关于lua的gc云风大佬在 Lua GC 的源码剖析 系列文章中讲得很清楚,这里做一下简单的记录。 分步gc lua使用的是一种三色标记清除算法(tri-color incremental mark & sweep),大体步骤如下: 初始阶段,所…...
基于python爬虫对豆瓣影评分析系统的设计与实现
基于python爬虫对豆瓣影评分析系统的设计与实现 Design and Implementation of a Python-based Web Crawler for Analyzing Douban Movie Reviews 完整下载链接:基于python爬虫对豆瓣影评分析系统的设计与实现 文章目录 基于python爬虫对豆瓣影评分析系统的设计与实现摘要第一…...

想让梦想照进现实?六西格玛绿带培训为你架起桥梁
六西格玛,这个源自摩托罗拉的质量管理方法论,如今已成为全球众多企业追求卓越的秘诀。它强调以数据为基础,通过减少变异和浪费,提高流程效率和质量,进而提升企业整体绩效。而六西格玛绿带培训,则是这个强大…...
大数据面试题之HDFS
目录 HDFS文件写入和读取流程 HDFS组成架构 介绍下HDFS,说下HDFS优缺点,以及使用场景 HDFS作用 HDFS的容错机制 HDFS的存储机制 HDFS的副本机制 HDFS的常见数据格式,列式存储格式和行存储格式异同点,列式存储优点有哪些? …...

(9)农作物喷雾器
文章目录 前言 1 必要的硬件 2 启用喷雾器 3 配置水泵 4 参数说明 前言 Copter 包括对农作物喷雾器的支持。该功能允许自动驾驶仪连接到一个 PWM 操作的泵和(可选)旋转器,根据飞行器速度控制液体肥料的流动速度。 稍微过时的视频显示了…...

智慧互联:Vatee万腾平台展现科技魅力
随着科技的迅猛发展,我们的生活正逐渐变得智能化、互联化。在这个信息爆炸的时代,一个名为Vatee万腾的平台正以其独特的魅力,引领我们走向一个更加智能的未来。 Vatee万腾,这个名字本身就充满了对科技未来的憧憬与期待。作为一家专…...

Charles抓包工具系列文章(四)-- Rewrite 重写工具
一、背景 这是一款比Map Local/Remote 还强大的工具,更加灵活,体现在以下几点: 重写request报文重写response报文header 字段的增删改query param 字段的增删改重写 body 字段改写http 响应状态status重写host/url/path 从这也可以看出其强…...

【PB案例学习笔记】-24创建一个窗口图形菜单
写在前面 这是PB案例学习笔记系列文章的第24篇,该系列文章适合具有一定PB基础的读者。 通过一个个由浅入深的编程实战案例学习,提高编程技巧,以保证小伙伴们能应付公司的各种开发需求。 文章中设计到的源码,小凡都上传到了gite…...
环境配置的相关问题
一、shap安装踩坑 遇到错误: A module that was compiled using NumPy 1.x cannot be run in NumPy 2.0.0 as it may crash. To support both 1.x and 2.x versions of NumPy, modules must be compiled with NumPy 2.0. Some module may need to rebuild instead…...

github配置可拉取项目到本地
首先配置用户名和邮箱: git config --global user.name 自己的名字git config --global user.email 自己的邮箱配置完之后检查一下: git config --global user.namegit config --global user.email如果提示的是自己配置好的名字和邮箱就Ok 然后拉取githu…...
谷歌浏览器插件
项目中有时候会用到插件 sync-cookie-extension1.0.0:开发环境同步测试 cookie 至 localhost,便于本地请求服务携带 cookie 参考地址:https://juejin.cn/post/7139354571712757767 里面有源码下载下来,加在到扩展即可使用FeHelp…...

【入坑系列】TiDB 强制索引在不同库下不生效问题
文章目录 背景SQL 优化情况线上SQL运行情况分析怀疑1:执行计划绑定问题?尝试:SHOW WARNINGS 查看警告探索 TiDB 的 USE_INDEX 写法Hint 不生效问题排查解决参考背景 项目中使用 TiDB 数据库,并对 SQL 进行优化了,添加了强制索引。 UAT 环境已经生效,但 PROD 环境强制索…...
C++.OpenGL (10/64)基础光照(Basic Lighting)
基础光照(Basic Lighting) 冯氏光照模型(Phong Lighting Model) #mermaid-svg-GLdskXwWINxNGHso {font-family:"trebuchet ms",verdana,arial,sans-serif;font-size:16px;fill:#333;}#mermaid-svg-GLdskXwWINxNGHso .error-icon{fill:#552222;}#mermaid-svg-GLd…...
JDK 17 新特性
#JDK 17 新特性 /**************** 文本块 *****************/ python/scala中早就支持,不稀奇 String json “”" { “name”: “Java”, “version”: 17 } “”"; /**************** Switch 语句 -> 表达式 *****************/ 挺好的ÿ…...
Device Mapper 机制
Device Mapper 机制详解 Device Mapper(简称 DM)是 Linux 内核中的一套通用块设备映射框架,为 LVM、加密磁盘、RAID 等提供底层支持。本文将详细介绍 Device Mapper 的原理、实现、内核配置、常用工具、操作测试流程,并配以详细的…...
iOS性能调优实战:借助克魔(KeyMob)与常用工具深度洞察App瓶颈
在日常iOS开发过程中,性能问题往往是最令人头疼的一类Bug。尤其是在App上线前的压测阶段或是处理用户反馈的高发期,开发者往往需要面对卡顿、崩溃、能耗异常、日志混乱等一系列问题。这些问题表面上看似偶发,但背后往往隐藏着系统资源调度不当…...

JVM 内存结构 详解
内存结构 运行时数据区: Java虚拟机在运行Java程序过程中管理的内存区域。 程序计数器: 线程私有,程序控制流的指示器,分支、循环、跳转、异常处理、线程恢复等基础功能都依赖这个计数器完成。 每个线程都有一个程序计数…...
C++.OpenGL (20/64)混合(Blending)
混合(Blending) 透明效果核心原理 #mermaid-svg-SWG0UzVfJms7Sm3e {font-family:"trebuchet ms",verdana,arial,sans-serif;font-size:16px;fill:#333;}#mermaid-svg-SWG0UzVfJms7Sm3e .error-icon{fill:#552222;}#mermaid-svg-SWG0UzVfJms7Sm3e .error-text{fill…...

MFC 抛体运动模拟:常见问题解决与界面美化
在 MFC 中开发抛体运动模拟程序时,我们常遇到 轨迹残留、无效刷新、视觉单调、物理逻辑瑕疵 等问题。本文将针对这些痛点,详细解析原因并提供解决方案,同时兼顾界面美化,让模拟效果更专业、更高效。 问题一:历史轨迹与小球残影残留 现象 小球运动后,历史位置的 “残影”…...

三分算法与DeepSeek辅助证明是单峰函数
前置 单峰函数有唯一的最大值,最大值左侧的数值严格单调递增,最大值右侧的数值严格单调递减。 单谷函数有唯一的最小值,最小值左侧的数值严格单调递减,最小值右侧的数值严格单调递增。 三分的本质 三分和二分一样都是通过不断缩…...