《数字图像处理》- 3. MATLAB实现图像复原

news/2024/11/29 18:38:38/
声明
本文仅在CSDN发布,其他均为盗版。请支持正版!

正版链接:https://blog.csdn.net/meenr/article/details/117605883

MATLAB实现图像复原

实验3、MATLAB实现图像复原

目录

  • 实验3、MATLAB实现图像复原
    • 1. 实验目的
    • 2. 实验内容
    • 3. MATLAB程序代码
    • 4. 结果及分析
    • 5. 典型的图像复原方法及其适用范围
    • 6. MATLAB程序代码
      • 途径一
      • 途径二

1. 实验目的

(1)理解退化模型。

(2)掌握常用的图像复原方法。

2. 实验内容

(1)选择一幅清晰的灰度图像,对该图像进行模糊化处理,然后分别采用逆滤波、维纳滤波和约束最小二乘方滤波对模糊图像进行复原,比较各种图像复原方法的复原效果。

(2)选择一幅清晰的灰度图像,对该图像进行模糊化处理并加入不同强度的高斯噪声,然后分别采用逆滤波、维纳滤波和约束最小二乘方滤波对退化图像进行复原,比较各种图像复原方法的复原效果。

3. MATLAB程序代码

(1)选择一幅清晰的灰度图像,对该图像进行模糊化处理,然后分别采用逆滤波、维纳滤波和约束最小二乘方滤波对模糊图像进行复原。

(2)选择一幅清晰的灰度图像,对该图像进行模糊化处理并加入高斯噪声,然后分别采用逆滤波、维纳滤波和约束最小二乘方滤波对退化图像进行复原。

4. 结果及分析

(1)本实验选择一幅清晰的灰度图像,对该图像进行模糊化处理,然后分别采用逆滤波、维纳滤波和约束最小二乘方滤波对模糊图像进行复原。原图像及其模糊化处理后的图像如图1所示。

模糊化应用点扩散函数fspecial('motion', LEN, THETA),以某一个角度逆时针旋转一定像素。函数中参数"LEN"是旋转的像素值,"THETA"为旋转角度。
在这里插入图片描述

图1 原图像与模糊图像

对模糊化的图像用噪信比滤波和自相关函数滤波处理后的图像如图2所示。
在这里插入图片描述

图2 噪信比滤波与自相关函数滤波后图像

对模糊化的图像用逆滤波和约束最小二乘方滤波处理后的图像如图3所示。
在这里插入图片描述

图3逆滤波与约束最小二乘方滤波后图像

由上述几幅图像可以看出,在没有噪声的情况下,三种复原方法的效果几乎是一样的。噪声为零时,维纳滤波,约束最小二乘方滤波可以认为就是逆滤波。

(2)选择一幅清晰的灰度图像,对该图像进行模糊化处理并加入高斯噪声,然后分别采用逆滤波、维纳滤波和约束最小二乘方滤波对退化图像进行复原。

原图像与加方差为0.001的噪声模糊图像如图4所示。

在这里插入图片描述

图4 原图像与噪声方差为0.001模糊图像

为模糊图像加方差为0.001的噪声时的四种滤波效果如下所示。其中维纳滤波的两种子方法,噪信功率比和自相关函数法滤波效果如图5所示。

在这里插入图片描述

图5 噪信比滤波与自相关函数滤波后图像

加方差为0.001的噪声时逆滤波和约束最小二乘方法的滤波效果如图6所示。

在这里插入图片描述

图6 逆滤波与约束最小二乘方滤波后图像

原图像与加方差为0.008的噪声模糊图像如图7所示。

在这里插入图片描述

图7 原图像与噪声方差为0.008模糊图像

为模糊图像加方差为0.008的噪声时的四种滤波效果如下所示。其中维纳滤波的两种子方法,噪信功率比和自相关函数法滤波效果如图8所示。

在这里插入图片描述

图8 噪信比滤波与自相关函数滤波后图像

加方差为0.008的噪声时逆滤波和约束最小二乘方法的滤波效果如图6所示。

在这里插入图片描述

图9 逆滤波与约束最小二乘方滤波后图像

由上文添加两种不同强度的噪声的滤波效果,比较复原结果可以看出,在有噪声时采用逆滤波进行复原效果很不好,在已知信噪比的情况的维纳滤波复原效果较次之,已知噪声和原图像的自相关函数的情况下的维纳滤波复原效果最佳。而对于高噪声,约束最小二乘方滤波要更好一些。

5. 典型的图像复原方法及其适用范围

  1. 逆滤波:适用于无噪声图像,一般对于有噪声的图像来说,直接逆滤波的性能是较差的。
  2. 最小均方差滤波(维纳)滤波:该方法建立在图像和噪声都是随机变量的基础上。未退化图像和噪声的功率谱必须是已知的。在没有噪声的情况下,维纳滤波可以看作是逆滤波。
  3. 约束最小二乘方滤波:仅仅用噪声的均值和方差的知识,就可以实现的最佳复原算法。但要假设噪声和图像灰度值不相关。
  4. 几何均值滤波:对维纳滤波的推广。
  5. 投影重建图像,卡尔曼滤波等。图像复原算法有线性和非线性两类。线性算法通过对图像进行逆滤波来实现反卷积。非线性方法通过连续的迭代过程不断提高复原质量,直到满足预先设定的终止条件,结果往往令人满意。

6. MATLAB程序代码

感兴趣的读者可获取参照下方方式获取实验代码和实验报告等资料。优先推荐途径一,若遇途径一失效,请再尝试途径二。

途径一

优先推荐该途径
第一步:扫描下方二维码,或打开微信搜索并关注“ 2贰进制 ”公众号;
第二步:回复“ MATLAB图像复原 ”即可获取本文相关资料。
二维码

途径二

优先推荐途径一,该途径管理可能不能秒回
扫描下方二维码,加入学习交流QQ群“ 480558240 ”,联系管理员获取包括但不限于本篇内容的更多学习资料。
在这里插入图片描述
2贰进制–Echo 2020年5月
我认同兴趣是最好的老师,但是除了兴趣其次是侮辱,所以如果您觉得本文还不错,请点赞+评论+收藏,要是关注那更是对我极大地羞辱了,您的羞辱便是我前进的动力!
如果本文对你有所帮助,解决了您的困扰,可以通过赞赏来给予我更大支持:
在这里插入图片描述
此致
感谢您的阅读、点赞、评论、收藏与打赏。


http://www.ppmy.cn/news/531913.html

相关文章

数字图像处理课程设计

数字图像处理课程设计 一、前言 《数字图像处理》是一门应用型课程,为了更巩固和掌握图像处理技术的基本技能,提高实际动手能力,并通过实际编程了解图像处理算法实现的基本原理,为今后应用图像处理技术和编程技术解决实际问题奠定…

C++【STL】之反向迭代器

反向迭代器 前面在vector和list的模拟实现中都有讲到正向迭代器,今天我们就来讲解一下反向迭代器的思想和模拟实现,在某些场景下还是很实用的,下面正文直接开始。 文章目录: 反向迭代器1. 反向迭代器结构2. 反向迭代器实现2.1 多…

关于Synchronized的小结

目录 一.特性 1.既是乐观锁又是悲观锁 2.是轻量级锁,也是重量级锁 3.不是读写锁,是互斥锁 4.是可重入锁 5.非公平锁 6.加锁之后,运行完毕自动解锁释放资源 二:Synchronized使用 第一种:定义一个额外的变量来控制加锁和解锁(类似于吉祥物) 第一种:直接给类/方法上锁 三…

L2-014 列车调度 (25 分)

L2-014 列车调度 (25 分) 火车站的列车调度铁轨的结构如下图所示。 两端分别是一条入口(Entrance)轨道和一条出口(Exit)轨道,它们之间有N条平行的轨道。每趟列车从入口可以选择任意一条轨道进入,最后从出口…

427-二叉树(617.合并二叉树、700.二叉搜索树中的搜索、98. 验证二叉搜索树、530.二叉搜索树的最小绝对差)

617.合并二叉树 class Solution { public:TreeNode* mergeTrees(TreeNode* root1, TreeNode* root2) {if (root1 nullptr) return root2;if (root2 nullptr) return root1;root1->val root2->val;root1->left mergeTrees(root1->left, root2->left);r…

【卡特兰数】HJ77.火车进站

题目&#xff1a; 描述 给定一个正整数N代表火车数量&#xff0c;0<N<10&#xff0c;接下来输入火车入站的序列&#xff0c;一共N辆火车&#xff0c;每辆火车以数字1-9编号&#xff0c;火车站只有一个方向进出&#xff0c;同时停靠在火车站的列车中&#xff0c;只有后进站…

210920-车站问题

mainByHuber.cpp 文件&#xff1a; #include <iostream> #include <string.h> #include<stdlib.h> struct station_info{char station_name[100]; }; struct station_node{station_info s;station_node * next; };int findStation(station_node * head,cha…

12306一直显示服务器忙,网购春运火车票首日 12306网站又现服务器忙

①青岛火车站售票大厅内&#xff0c;购票回家的旅客已明显增多。 ②来自广西的杨先生和同伴们准备拎着行李提前回家。据了解&#xff0c;他们在船厂工作&#xff0c;由于天冷活少&#xff0c;所以提前放假了。 ③"小候鸟"在妈妈的背后等候进站上车。 本版图/记者 孙传…