《数字图像处理》实验报告四
一、实验任务与要求
- 对 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…...
Android Wi-Fi 连接失败日志分析
1. Android wifi 关键日志总结 (1) Wi-Fi 断开 (CTRL-EVENT-DISCONNECTED reason3) 日志相关部分: 06-05 10:48:40.987 943 943 I wpa_supplicant: wlan0: CTRL-EVENT-DISCONNECTED bssid44:9b:c1:57:a8:90 reason3 locally_generated1解析: CTR…...

调用支付宝接口响应40004 SYSTEM_ERROR问题排查
在对接支付宝API的时候,遇到了一些问题,记录一下排查过程。 Body:{"datadigital_fincloud_generalsaas_face_certify_initialize_response":{"msg":"Business Failed","code":"40004","sub_msg…...
进程地址空间(比特课总结)
一、进程地址空间 1. 环境变量 1 )⽤户级环境变量与系统级环境变量 全局属性:环境变量具有全局属性,会被⼦进程继承。例如当bash启动⼦进程时,环 境变量会⾃动传递给⼦进程。 本地变量限制:本地变量只在当前进程(ba…...
在rocky linux 9.5上在线安装 docker
前面是指南,后面是日志 sudo dnf config-manager --add-repo https://download.docker.com/linux/centos/docker-ce.repo sudo dnf install docker-ce docker-ce-cli containerd.io -y docker version sudo systemctl start docker sudo systemctl status docker …...
Golang dig框架与GraphQL的完美结合
将 Go 的 Dig 依赖注入框架与 GraphQL 结合使用,可以显著提升应用程序的可维护性、可测试性以及灵活性。 Dig 是一个强大的依赖注入容器,能够帮助开发者更好地管理复杂的依赖关系,而 GraphQL 则是一种用于 API 的查询语言,能够提…...
土地利用/土地覆盖遥感解译与基于CLUE模型未来变化情景预测;从基础到高级,涵盖ArcGIS数据处理、ENVI遥感解译与CLUE模型情景模拟等
🔍 土地利用/土地覆盖数据是生态、环境和气象等诸多领域模型的关键输入参数。通过遥感影像解译技术,可以精准获取历史或当前任何一个区域的土地利用/土地覆盖情况。这些数据不仅能够用于评估区域生态环境的变化趋势,还能有效评价重大生态工程…...
3403. 从盒子中找出字典序最大的字符串 I
3403. 从盒子中找出字典序最大的字符串 I 题目链接:3403. 从盒子中找出字典序最大的字符串 I 代码如下: class Solution { public:string answerString(string word, int numFriends) {if (numFriends 1) {return word;}string res;for (int i 0;i &…...

什么是Ansible Jinja2
理解 Ansible Jinja2 模板 Ansible 是一款功能强大的开源自动化工具,可让您无缝地管理和配置系统。Ansible 的一大亮点是它使用 Jinja2 模板,允许您根据变量数据动态生成文件、配置设置和脚本。本文将向您介绍 Ansible 中的 Jinja2 模板,并通…...
在web-view 加载的本地及远程HTML中调用uniapp的API及网页和vue页面是如何通讯的?
uni-app 中 Web-view 与 Vue 页面的通讯机制详解 一、Web-view 简介 Web-view 是 uni-app 提供的一个重要组件,用于在原生应用中加载 HTML 页面: 支持加载本地 HTML 文件支持加载远程 HTML 页面实现 Web 与原生的双向通讯可用于嵌入第三方网页或 H5 应…...
React---day11
14.4 react-redux第三方库 提供connect、thunk之类的函数 以获取一个banner数据为例子 store: 我们在使用异步的时候理应是要使用中间件的,但是configureStore 已经自动集成了 redux-thunk,注意action里面要返回函数 import { configureS…...