利用课程讲的滤波方法,自行设计滤波模板的内容和尺寸,分别对两个噪声图像进行去噪处理,得到的去噪结果,统一命名为denoise_lena_guassian和denoise_lena_salt_pepper。分别将去噪后的结果与原始图像(即没有噪声的lena_gray)进行比较。mTLb代码

时间: 2024-02-22 19:54:57 浏览: 23
好的,以下是mTLb代码实现: ```matlab % 读取图像 lena_gray = imread('lena_gray.jpg'); lena_guassian = imread('lena_guassian.jpg'); lena_salt_pepper = imread('lena_salt_pepper.jpg'); % 高斯滤波 gaussian_filter = fspecial('gaussian', [5 5], 2); denoise_lena_guassian = imfilter(lena_guassian, gaussian_filter); % 中值滤波 denoise_lena_salt_pepper = medfilt2(lena_salt_pepper, [3 3]); % 计算评价指标 psnr_guassian = psnr(denoise_lena_guassian, lena_gray); ssim_guassian = ssim(denoise_lena_guassian, lena_gray); psnr_salt_pepper = psnr(denoise_lena_salt_pepper, lena_gray); ssim_salt_pepper = ssim(denoise_lena_salt_pepper, lena_gray); % 显示结果 subplot(2, 3, 1), imshow(lena_gray), title('原始图像'); subplot(2, 3, 2), imshow(lena_guassian), title('高斯噪声图像'); subplot(2, 3, 3), imshow(denoise_lena_guassian), title(sprintf('去噪图像\nPSNR=%.2f\nSSIM=%.2f', psnr_guassian, ssim_guassian)); subplot(2, 3, 4), imshow(lena_salt_pepper), title('椒盐噪声图像'); subplot(2, 3, 5), imshow(denoise_lena_salt_pepper), title(sprintf('去噪图像\nPSNR=%.2f\nSSIM=%.2f', psnr_salt_pepper, ssim_salt_pepper)); subplot(2, 3, 6), imshow(lena_gray-denoise_lena_salt_pepper), title('差异图像'); ``` 这里假设原始图像为lena_gray.jpg,高斯噪声图像为lena_guassian.jpg,椒盐噪声图像为lena_salt_pepper.jpg。代码中首先使用fspecial函数生成一个5x5的高斯滤波器,标准差为2,然后使用imfilter函数对高斯噪声图像进行滤波,得到去噪图片denoise_lena_guassian。接着,使用medfilt2函数对椒盐噪声图像进行中值滤波,得到去噪图片denoise_lena_salt_pepper。最后,使用psnr和ssim函数计算评价指标,并使用subplot和imshow函数将原始图像、噪声图像、去噪图像和差异图像显示在画布上。 希望这份代码能够帮到您!如果您还有其他问题,请随时提出。

相关推荐

最新推荐

recommend-type

数字图像处理实验报告-数字图像空间与频率滤波.docx

仅对相位部分进行傅立叶反变换后查看结果图像(提示:记傅立叶变换的相位a,利用ifft2对exp(a*i)进行反变换);仅对幅度部分进行傅立叶反变换后查看结果图像;并图像显示; 编写函数w = genlaplacian(n),自动产生...
recommend-type

python数字图像处理之高级滤波代码详解

主要介绍了python数字图像处理之高级滤波代码详解,介绍了许多对图像处理的滤波方法,具有一定参考价值,需要的朋友可以了解下。
recommend-type

Python实现中值滤波去噪方式

今天小编就为大家分享一篇Python实现中值滤波去噪方式,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
recommend-type

ADSP-BF533的例程源码29:FIR信号处理源码

ADSP-BF533的例程源码29:FIR信号处理源码
recommend-type

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节 线性代数的核心问题是求解方程组。这些方程都是线性的,即未知数仅与数相乘——我们绝不会 遇见 x 乘以 y。我们的第一个线性方程组较小。接下来你来看看它引申出多远: 两个方程 两个未知数 x − 2y = 1 3x + 2y = 11 (1) 我们一次从一个行开始。第一个方程 x − 2y = 1 得出了 xy 平面的一条直线。由于点 x = 1, y = 0 解 出该方程,因此它在这条直线上。因为 3 − 2 = 1,所以点 x = 3, y = 1 也在这条直线上。若我们选择 x = 101,那我们求出 y = 50。 这条特定直线的斜率是 12,是因为当 x 变化 2 时 y 增加 1。斜率在微积分中很重要,然而这是线 性代数! 图 2.1 将展示第一条直线 x − 2y = 1。此“行图”中的第二条直线来自第二个方程 3x + 2y = 11。你 不能错过两条线的交点 x = 3, y = 1。点 (3, 1) 位于两条线上并且解出两个方程。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

识别MATLAB微分方程求解中的混沌行为:分析非线性方程混沌行为的实用技巧

![matlab求解微分方程](https://img-blog.csdnimg.cn/2021062810300367.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MTQ3OTY3OA==,size_16,color_FFFFFF,t_70) # 1. MATLAB微分方程求解概述 微分方程是描述物理、工程和金融等领域中动态系统的数学模型。MATLAB提供了强大的工具来求解微分方程,包括内置求解器和自定义函数
recommend-type

physon如何做ERP系统

作为一个AI语言模型,我可以给您提供一些基本的建议和思路: 1.明确需求和目标:在开始开发ERP系统之前,需要明确您的业务需求和目标,了解您的业务流程、流程数据、员工角色和权限、财务管理等方面的要求,这样才能更好的设计和开发ERP系统。 2.选择技术框架和开发工具:选择合适的技术框架和开发工具是ERP系统开发的关键。选择一种流行的技术框架和工具可以提高开发效率和质量。 3.设计数据库:ERP系统需要一个功能强大的数据库来存储数据。设计数据库需要考虑数据的完整性、安全性和可扩展性。 4.设计系统架构:系统架构是ERP系统的骨架,需要考虑系统的可扩展性、可维护性和性能。 5.开发和测试:
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依