paper 108:系统学习数字图像处理之图像复原与重建

首先,必须注意这里所限制的处理条件。

关于图像退化/复原模型

退化的图像是由成像系统的退化加上额外的噪声形成的。

1.只考虑噪声引起的退化

噪声模型,包含于空间不相关和相关两种,除了空间周期噪声,这里所讨论的都是空间不相关的噪声,如高斯,爱尔兰,瑞利,指数分布,均匀分布,脉冲(椒盐)噪声等。

针对只存在噪声引起的退化,首先要估计噪声参数,然后估计噪声模型,接着做噪声滤除。此处可以选择空间滤波方法,图像增强与复原没有区别。几类滤波器的效果特性如下:

1.均值滤波

当图像仅存在加性噪声时,采用均值滤波

算数均值滤波,模糊图像,复原效果最差

几何均值滤波,与算数均值平滑度差不多,但细节丢失更少,

以上两种适合处理高斯或均匀噪声分布。

谐波均值滤波,适用于处理高斯和均匀分布噪声处理,在白点(盐噪声)效果好,不适合胡椒噪声。

逆谐波均值滤波,根据阶数处理椒/盐噪声,不能同时消除椒盐噪声。阶数为0,变成算数均值,为-1,变成谐波均值滤波器。如果Q选择不当,会产生严重后果。

2.顺序统计滤波

中值滤波

特点:1)在相同尺寸下,比起均值滤波器引起的模糊少;

2)对单极或双极脉冲(椒盐)噪声非常有效。只要脉冲噪声的空间密度不大,经验为(小于0.2)。

3)适于处理椒盐噪声,通过多次使用小模板,可以获得很好的去噪效果,但多次应用中值滤波器,会使图像模糊。

最大值滤波

特点:1)对图像中的亮点敏感;

2)最大值滤波器对于“椒”噪声具有良好消除效果。

3)最大值滤波器可以去除”胡椒”噪声,但会从黑色物体边缘移走一些黑色像素。

最小值滤波

1) 对图像中的暗点敏感;

2) 最小值滤波器对于“盐”噪声具有良好消除效果。

3)最小值滤波器可以去除”盐”噪声,但会从亮色物体边缘移走一些白色像素。

中点滤波,

这种滤波器结合了顺序统计和求平均 (均匀),对于高斯和均匀随机分布噪声有最好的效果。

修正后的α均值滤波

当滤波器非退变为算数均值和中值滤波器的情况下,非常适合于高斯和椒盐噪声混合的情况。

3.自适应滤波

自适应滤波优于以上讨论的所有滤波器。

自适应局部噪声消除滤波,需要估计噪声方差。但要注意噪声方差和图像方差比率超过1时的处理问题,一种是超过1限制在1,这会导致滤波的非线性,但可以防止负值,另外一种是允许出现负值,但最后要重新标定灰度值,但结果是损失动态范围。

自适应中值滤波,相比传统中值滤波有更大概率的处理空间,另外可以平滑非脉冲噪声时试图保留细节,减少图像边界的细化和粗话失真。

4,频域滤波

带阻,带通,陷波已经在之前讲过,不赘述。这里只增加最佳陷波滤波器。

通常不会清楚定义干扰模式,有几种干扰同时存在的时,且干扰成分不在是单频脉冲,前面的方法就不在适用。这里采用最佳的,最小化复原估计值的局部方差。

首先要提取噪声主频成分,得到噪声函数,再次选取加权函数(或者叫调制函数),然后以某种有意义的方法使得估计值在每一点指定领域上的方差最小。

2.噪声和成像系统退化同时存在

1.先说线性,位置不变的退化

线性,输入之和的响应等于输入响应的和,位置不变,表示只与输入值有关,与位置无关。

所以,在有加性噪声,且噪声随机,与位置无关时,有结论:对于具有加性噪声的线性空间不变退化系统,可表示为退化函数与图像的卷积,再加上噪声。

许多类型的退化可近似为线性,位置不变的过程,这样可以用许多线性工具来解决图像复原问题,与位置有关的非线性技术虽然更普遍,但会带来没有已知解的困难,或者解决计算问题时非常困难,这里主要讨论线性,位置不变系统的复原技术。所以,图像去卷积通常用于表示线性图像复原,用于复原处理的滤波器,通常称为去卷积滤波器。

2.估计退化函数。

观察估计,选择一个强信号区域,这里受噪声的干扰可以忽略。接着处理该区域得到尽可能不模糊的结果。两者的DFT之比即可近似为退化函数。

试验估计,模拟出与获得退化图像的设备相似的装置,用冲激成像。

建模估计,数学建模,如大气湍流,匀速运动物体的退化函数模型等。

3.图像复原

忽略噪声后,退化模型的DFT,其中H已知,则IDFT(F)得到复原图像,这就是逆滤波。实际中,H为0或者接近0时,F变成无穷大或者非常大的数。

不忽略噪声时,不能 准确复原,因为噪声的DFT未知,如果退化是0或者接近0,噪声数值即使很小,但还是很容易决定F的估计值,因为复原公式中噪声复原引起的数值很大。

解决退化是0或者接近0的方法是,限制滤波频率,使其接近原点。即图像复原在H的原点周围的有限区域进行。

这种一般的直接逆滤波,效果较差,没有说明如何处理噪声。下面改进的方法是:

维纳滤波

它会综合考虑退化函数和噪声的统计特征。即寻找一个滤波器使其复原后的图像与原图像均方误差最小。因此,维纳滤波又叫最小均方误差复原,也叫最小二乘误差滤波器。注意,维纳滤波假设噪声和图像不相关,未污染图像和估计图像中有一个有零均值,且估计的灰度级是退化图像灰度级的线性函数。如果噪声均值为0,则维纳滤波器变成逆滤波器。信噪比,是以噪声和未退化图像的功率谱为基础的定义的。然而,未退化图像的功率谱很少已知,则常用常数K来估计。

因此维纳滤波的缺点是,

  1. 未退化图像和噪声的功率谱必须是已知的;
  2. 功率比(信噪比)常数K的估计一般还是没有合适的解。

因此,另外一种去卷积滤波器是约束最小二乘方滤波器。

它是通过将图像复原模型,用矩阵的形式表示,其核心是减少噪声敏感性问题,因此,以平滑度量的最佳复原为基础,如拉普拉斯变换的锐化。然后,利用约束条件,不断调整参数γ值。

经验是,对高噪声和中等噪声,最小二乘方要比维纳效果好,低噪声时类似,手工选择参数时,最小二乘方可能比维纳要好。参数值,可以迭代计算,但自动计算出的参数效果要比手动的差。

几何均值滤波

维纳滤波的推广,根据参数值,可以退变为维纳滤波,逆滤波,谱均衡滤波器。

时间: 2024-11-05 22:49:22

paper 108:系统学习数字图像处理之图像复原与重建的相关文章

学习数字图像处理_1

结束了数据结构的内容之后,在图书馆接着看数字图像处理的内容.教材以冈萨雷斯编写的<数字图像处理(第三版)>和<数字图像处理(matlab版)>为主.前期看的是后者,后者较薄,且以介绍matlab的各种图像处理相关函数和编程为主,原理讲的较少,前4章还能大概看懂程序,到图像复原这章由于基础知识不到位看的比较吃力,所以网上买了纯理论的第一本书.两本结合着看,前期还是以基础知识为主吧,matlab编程本身难度不大,主要是理解算法原理,虽然大部分算法matlab已经编好放在map工具箱中的

数字图像处理(MATLAB版)学习笔记——第1章 绪言

0.下定决心 当当入手数字图像处理一本,从此开此正式跨入数字图像处理大门.以前虽然多多少少接触过这些东西,也做过一些相关的事情,但感觉都不够系统,也不够专业,从今天开始,一步一步地学习下去,相信会有成果的. 1.本章知识点结构 2.书中例子实践 第一个例子,试一试手. 1 function f = twodsin(A, u0, v0, M, N) 2 f = zeros(M, N); 3 for c = 1:N 4 v0y = v0 * (c - 1); 5 for r = 1:M 6 u0x

数字图像处理学习笔记之一 DIP绪论与MATLAB基础

写在前面的话 数字图像处理系列的学习笔记是作者结合上海大学计算机学院<数字图像处理>课程的学习所做的笔记,使用参考书籍为<冈萨雷斯数字图像处理(第二版)(MATLAB版)>,同时学习过程中会参考网络学习资源.对于数字图像处理的学习不可能仅仅依靠作者所写的这一系列笔记,而是需要花时间和精力学习,本文只可作参考和交流之用.由于涉及此学科不久,在学习过程中难免存在错误,请读者不吝赐教. 数字图像处理绪论 数字图像处理(DIP)的研究目标和处理对象: DIP的研究目标是获取信息,处理对象是

[数字图像处理]图像复原--逆滤波

1.逆滤波的问题点 图像的老化,可以视为以下这样的一个过程.一个是退化函数的影响(致使图片模糊,褪色等),一个可加性噪声的影响. 用算式表示为 前几篇博文,主要是介绍可加性噪声的去除.本博文,主要介绍图像的逆滤波,即退化函数的去除.然而,逆滤波在空间域内的处理是很不方便的. 简单的来考虑,加法的逆运算是减法,乘法的逆运算的除法,微分的逆运算是积分(严密一点说是不定积分).那么就可以得到一个简单的结论了,要出去卷积的话,肯定需要用到卷积的逆运算.卷积的逆运算是---------反卷积,额,好像是一

数字图像处理—学习总结

VC++数字图像处理 杨淑滢 本文对该老师的代码进行总结,代码小白顺便补充一下自己比较困惑的知识点,方便后续复习. 对于BMP格式的图像结构分为以下四个部分: 位图文件头结构BITMAPFILEHEADER 位图信息头结构BITMAPINFOHEADER 位图颜色表RGBQUAD 位图像素数据部分 对于位图的图像分为四类: 单色图像:有两个值:0,255其中0表示黑,255表示白(也可以进行归一化处理为0和1). 灰度图像:带有颜色表,相当于索引表.RGB的分量值相同. 伪彩色图像:RGB的分量

大数据分析与机器学习系统学习路线指导(1)

3.开始专业课程的学习,请根据自己的学习能力和基础,逐渐选择学习课程,一年完成学习即可!三个专业方向课程详解如下: 数据分析师专业方向 学习路线建议 课程名称 难易程度 基础课 大数据的统计学基础 ★ 基础课 大数据的矩阵计算基础 ★ 工具课 深度玩转Excel ★ 工具课 R语言数据分析.展现与实例 ★ 工具课 SPSS数据分析入门与提高 ★★ 工具课 数据分析与SAS ★★ 应用课 基于软件学习数据挖掘算法与案例 ★ 应用课 BI实战之cognos 篇 ★★ 应用课 量化投资 ★★ 应用课

数字图像处理 简单介绍

 数字图像处理发展概况 数字图像处理发展概况数字图像处理(Digital Image Processing)又称为计算机图像处理,它是指将图像信号转换成数字信号并利用计算机对其进行处理的过程.数字图像处理最早出现于20世纪50年代,当时的电子计算机已经发展到一定水平,人们開始利用计算机来处理图形和图像信息.数字图像处理作为一门学科大约形成于20世纪60年代初期.早期的图像处理的目的是改善图像的质量,它以人为对象,以改善人的视觉效果为目的.图像处理中,输入的是质量低的图像,输出的是改善质量后的图像

【数字图像处理之(一)】数字图像处理与相关领域概述

数字图像(Digital Image) 一副图像可以定义为一个二维函数f(x, y),这里的x和y是空间坐标,而在任意坐标(x, y)处的幅度f被称为这一坐标位置图像的亮度或灰度.当x.y和f的幅值都是有限的离散值时,称为数字图像.注意,数字图像由有限数量的元素组成,每个元素都有特殊的位置和数值.这些元素称为画像元素.图像元素和像素,像素是定义数字图像元素时使用最广泛的术语. --Digital Image Pricessing Using MATLAB( Rafacel C. Gonzalez

数字图像处理原理与实践(MATLAB版)勘误表

本文系<数字图像处理原理与实践(MATLAB版)>一书的勘误表. [内容简单介绍]本书全面系统地介绍了数字图像处理技术的理论与方法,内容涉及几何变换.灰度变换.图像增强.图像切割.图像去噪.小波变换.形态学处理.多尺度融合.偏微分方程应用.正交变换与图像压缩.边缘及轮廓检測.图像复原.图像去雾.多尺度空间构建与特征匹配等15大核心话题.全部算法均配有完整的MATLAB实现代码.并以此为基础具体介绍了MATLAB中与图像处理有关的近200个函数的用法,便于读者学习与实践.此外,本书还提供了丰富的