基于PDE的图像降噪方法的研究综述
- 格式:doc
- 大小:200.50 KB
- 文档页数:6
图像去噪技术研究背景及方法1背景21世纪是信息化的时代,信息的形式不再是单纯的语音,而是发展到包括数据、文字、图像、视频等在内的多媒体形式。
据统计,人类接受外界的信息中有70%来自于图像。
图像处理技术在人类生产和生活的方方面面起到了越来越重要的作用。
例如,人们在网络上浏览、下载、共享丰富的图像、视频等多媒体信息;医生根据核磁共振扫描图像对患者进行疾病诊断;科学家根据卫星遥感图像对矿产进行定位和预测等。
图像的广泛应用对图像的表示方法、处理模型和算法提出了新的发展要求。
近年来,在数学分析、计算机视觉、模式识别和统计分析等学科中,分别独立地发展着一种彼此极其相似的理论,人们称之为:多尺度几何分析(MGA,Multiscale Geometric Analysis)。
发展多尺度几何分析的目的就是为了检测、表示、处理某些高维空间数据,这些空间的主要特点是:其中数据的某些重要特征集中体现于其低维子集中(如曲线、面等)。
在指纹图像去噪中,特别考虑减小各种元器件带来的噪声,主要包括了高斯白噪声和交沿噪声,其中交沿噪声对指纹图像影响较为突出。
傅立叶分析揭示了时域与频域之间的内在联系,反应了信号在“整个”时间范围内的“全部”频谱成分,是研究周期现象不可或缺的工具。
然而,傅立叶变换虽然具有很强的频域局部化能力,但并不具有时域局部化能力,而后一点,对于很多信号处理工作而言,特别是对于涉及非平稳信号处理的任务而言,是至关重要的。
小波分析的理论和方法是从傅立叶分析演变而来的。
小波变换以牺牲部分频域定位性能来取得时频局部性能的折衷,不仅能够提供较精确的频域定位,也能够提供较精确的时域定位。
我们所面对的真实物理信号,更多地表现出非平稳的特性,而小波变换恰恰是处理非平稳信号的有力工具。
小波理论的兴起,得益于其对信号的时频局部分析能力及其对一维有界变差函数的最优逼近性能,也得益于Mallat和Meyer等人引入的多分辨分析概念,以及Mallat提出的快速小波变换实现方法。
毕业论文基于小波变换的图像去噪方法的研究学生姓名: 学号:学系 专 指导教师:2011年 5 月基于小波变换的图像去噪方法的研究摘要图像是人类传递信息的主要媒介。
然而,图像在生成和传输的过程中会受到各种噪声的干扰,对信息的处理、传输和存储造成极大的影响。
寻求一种既能有效地减小噪声,又能很好地保留图像边缘信息的方法,是人们一直追求的目标。
小波分析是局部化时频分析,它用时域和频域联合表示信号的特征,是分析非平稳信号的有力工具。
它通过伸缩、平移等运算功能对信号进行多尺度细化分析,能有效地从信号中提取信息。
随着小波变换理论的完善,小波在图像去噪中得到了广泛的应用,与传统的去噪方法相比小波分析有着很大的优势,它能在去噪的同时保留图像细节,得到原图像的最佳恢复。
本文对基于小波变换的图像去噪方法进行了深入的研究分析,首先详细介绍了几种经典的小波变换去噪方法。
对于小波变换模极大值去噪法,详细介绍了其去噪原理和算法,分析了去噪过程中参数的选取问题,并给出了一些选取依据;详细介绍了小波系数相关性去噪方法的原理和算法;对小波变换阈值去噪方法的原理和几个关键问题进行了详细讨论。
最后对这些方法进行了分析比较,讨论了它们各自的优缺点和适用条件,并给出了仿真实验结果。
在众多基于小波变换的图像去噪方法中,运用最多的是小波阈值萎缩去噪法。
传统的硬阈值函数和软阈值函数去噪方法在实际中得到了广泛的应用,而且取得了较好的效果。
但是硬阈值函数的不连续性导致重构信号容易出现伪吉布斯现象;而软阈值函数虽然整体连续性好,但估计值与实际值之间总存在恒定的偏差,具有一定的局限性。
鉴于此,本文提出了一种基于小波多分辨率分析和最小均方误差准则的自适应阈值去噪算法。
该方法利用小波阈值去噪基本原理,在基于最小均方误差算法LMS和Stein无偏估计的前提下,引出了一个具有多阶连续导数的阈值函数,利用其对阈值进行迭代运算,得到最优阈值,从而得到更好的图像去噪效果。
基于变分PDE的地震资料空间自适应保边缘去噪方法王德华;高静怀;张丽丽【摘要】地震资料去噪是地震资料处理的一个基本环节,其信噪比的好坏会直接影响地震资料的可靠性及地质解释的精度.本文主要针对已有方法在去噪的同时不能很好地保持边缘信息这一问题,通过构造空间自适应边缘检测函数,建立了基于分数阶偏微分方程的自适应变分去噪模型;然后,推导了求解新模型的分数阶Euler-Lagrange方程,并给出其离散格式;最后,将本文提出的方法用于合成地震数据及实际地震数据去噪,通过与传统去噪方法相比较,验证了新方法的有效性与实用性.【期刊名称】《工程数学学报》【年(卷),期】2019(036)001【总页数】10页(P33-42)【关键词】地震资料去噪;边缘保持;变分正则化;偏微分方程【作者】王德华;高静怀;张丽丽【作者单位】西安工业大学理学院,西安710021;西安交通大学电信学院,西安710049;西安工业大学理学院,西安710021【正文语种】中文【中图分类】O29;P6311 引言随着我国油气勘探程度的不断提高,勘探目标逐渐向深层和地表复杂地区转移,面临的油气藏多为地层岩性油气藏,这对地震资料处理解释工作提出更高的要求,而地震资料的信噪比问题解决的好坏,将会直接影响地震资料的可靠性以及地质解释的精度.因此,在地震勘探的各个环节中,压制噪声干扰、提高资料的信噪比是一项贯彻始终的关键任务.常规地震资料处理中使用的去噪技术有小波变换去噪[1]、f-x 域预测滤波[2]、KL 变换[3]、SVD 分解[4]等,小波变换具有很强的去数据相关特性,通过将信号能量投影在小波域少数大系数上,而噪声能量分布于整个小波域内,进行阈值去噪;f-x 预测滤波会对相干信号进行增强,但由于高频段信噪比较低,求取的预测算子受噪声影响较大,从而使得滤波处理后的高频段有效信号发生畸变,降低信号的保真度;KL 变换与SVD 分解主要利用相邻道信号在相同时刻的相关性,所以对水平同相轴有较好的去噪效果,而处理倾斜或者弯曲同相轴时效果不佳.这些方法大多利用了信号的空间相干特性,以牺牲横向分辨率达到提高信噪比的目的,容易使倾斜和弯曲同相轴受到衰减,并且会模糊和压制一些细微的信号结构,如小断裂、细河道等,甚至还能引起较大断距的断层两侧同相轴的错误连接,给精细的地震结构解释带来困难.因而,能够较好地保护地震数据的细微结构与断层、裂缝等边缘结构的叠后噪声衰减技术受到地球物理工作者的重视[5-7].偏微分方程(PDEs)是一类重要的数学模型,具有各向异性扩散性能,并且整个扩散过程在局部信息的作用下进行,不同方程的处理特性由具体的扩散项和扩散方向决定,方程直接使用图像中的几何特征(比如梯度、曲率、水平集等)控制扩散项及扩散方向,因此偏微分方程处理图像可在平滑同质区域的同时保持区域边界等几何特征[8],这是传统的统计方法与小波技术所缺少的.在计算机视觉和图像处理领域的许多问题,例如图像的去噪、增强、分割、跟踪等等都可以根据变分原理转化为某个泛函的最优化问题.习惯上,称这种从变分模型推导出的偏微分方程为变分偏微分方程.基于变分PDE的图像处理模型主要有两大类:一类是Mumford-shah 模型[9]及其改进模型[10,11],另一类是基于全变差(TotalVariation,TV)的模型[12]及其改进模型[13,14].1992 年,Rudin 等(ROF)提出用BV 空间对图像进行建模,建立的基于TV的非线性变分模型式中BV(Ω)表示有界变差函数空间,L2(Ω)表示平方可积函数空间,第一项为正则化项,第二项为保真项,u0 表示含噪图像,u表示去噪后图像,λ >0 为正则化参数.该模型的优点是去噪同时可以保护边缘,缺点是去噪后的图像产生阶梯效应,去噪同时滤掉了纹理等细节信息[15].为了克服ROF 去噪模型的“阶梯效应”这一缺陷,一些学者通过引入高阶PDEs 来处理信号能量的缓变区域,从而增强信号平滑区域的连续性,其中的代表性工作之一就是Lysaker 等在2003 年提出了一种基于四阶PDE的去噪模型(简称LLT 模型)[16]式中,参数λ >0 用来调节模型中光滑项与拟合项之间的权重,模型中其它参数的意义与模型(1)中的对应参数的意义相同.LLT 模型大大缓解了ROF 模型带来的“阶梯效应”,但是该模型与ROF 模型的一个共同缺陷是:在使用二者进行去噪时,包含在残差项u0 −u中的纹理信息和噪声不加区分地一起被滤除,因此,在去噪的同时不能有效地保持图像中的纹理、细节等有用信息.2008 年,Zhang 等人在文献[17]中,根据分数阶导数算子的滤波特性,通过改变ROF 模型中的正则化项,将其中的整数阶导数改为分数阶导数,给出下列分数阶变分去噪模型式中表示图像u(x,y)关于变量x的α阶导数;表示图像u(x,y)关于变量y的α阶导数;其他变量的含义和上述ROF 模型中相应变量的含义相同.注意到,当α=1 时,上述模型(3)简化为ROF 模型.本文主要针对上述分数阶变分去噪模型的缺点,通过构造自适应边缘检测函数,提出了一种能够更好地保持数据的边缘信息、抑制“阶梯效应”的自适应分数阶变分PDE 去噪模型,推导了模型求解的Euler-Lagrange 方程,并将模型应用于合成地震数据、实际地震数据去噪,验证了方法的有效性.本文第2 部分为模型的建立与分析,第3 部分给出模型求解的数值算法,第4 部分为数值实验,最后,给出本文的结论与展望.2 模型的建立与分析为了在去噪的同时,既能较好的保持地震数据的纹理等细节特征也能保护其中的边缘信息,下面我们首先构造自适应边缘检测函数,进而建立新的空间自适应分数阶变分去噪模型.2.1 自适应边缘检测函数的构造1990 年,Perona 和Malik 在文献[18]中首先提出了基于非线性扩散的边缘检测方法.文中考虑了如下PDE在该扩散方程中,g表示扩散函数,是一个非负递减函数,而且满足0.这里,保持边缘的思路是在边缘附近减小滤波程度,从而构造了如下边缘检测函数然而,这种边缘检测函数不适用于噪声存在的情况.1992 年,Lions 等人在文献[19]中针对上述边缘检测函数的缺陷,在方程(4)中引入高斯滤波算子,考虑如下PDE从而,边缘检测函数可以构造为该边缘检测函数能够识别含噪数据中的边缘,但其中的高斯滤波算子会导致边缘的错位[20].本文针对含噪数据,首先采用TV 扩散流[21]对其进行滤波,该非线性扩散过程可以在压制噪声的同时保持边缘信息.在此基础上,构造边缘检测函数如下。
小波图像去噪文献综述Overview of Wavelet Image Denoising0.简要说明本课题主要研究方向是通过小波变换的方法对图像进行去噪处理,以期得到能反应实际特征的容易分析的图像。
本综述介绍了小波图像去噪的原理,主要方法,并对小波图像去噪的未来进行了分析与展望。
根据本课题的特点,需要检索的文献在小波变换的研究方面主要检索1992年以来的国内外的期刊和学位论文,而在小波去噪及小波变换发展方面主要检索2000年以来的国内外的期刊以期得到最新的信息。
1.前言在信号数据采集及传输时,不仅能采集或接收到与所研究的问题相关的有效信号,同时也会观测到各种类型的噪声。
在实际应用中,为降低噪声的影响,不仅应研究信号采集的方式方法及仪器的选择,更重要的是对已采集或接收的信号寻找最佳的降噪处理方法。
对于信号去噪方法的研究可谓是信号处理中一个永恒的话题。
传统的去噪方法是将被噪声污染的信号通过一个滤波器,滤除掉噪声频率成分。
但对于瞬间信号、宽带噪声信号、非平稳信号等,采用传统方法具有一定的局限性。
其次还有傅里叶(Fourier)变换也是信号处理中的重要手段。
这是因为信号处理中牵涉到的绝大部分都是语音或其它一维信号,这些信号可以近似的认为是一个高斯过程,同时由于信号的平稳性假设,傅立叶交换是一个很好的信号分析工具。
但也有其不足之处,给实际应用带来了困难。
小波变换是继Fourier 变换后的一重大突破,它是一种窗口面积恒定、窗口形状可变(时间域窗口和频率域窗口均可改变)的时频局域化分析方法,它具有这样的特性;在低频段具有较高的频率分辨率及较低的时间分辨率,在高频段具有较高的时间分辨率及较低的频率分辨率,实现了时频窗口的自适应变化,具有时频分析局域性。
小波变换的一个重要应用就是图像信号去噪。
将小波变换用于信号去噪,它能在去噪的同时而不损坏信号的突变部分。
在过去的十多年,小波方法在信号和图像去噪方面的应用引起学者广泛的关注。
基于人工智能的图像去噪与图像增强算法研究图像去噪与图像增强是计算机视觉领域中重要的研究方向之一。
近年来,随着人工智能技术的快速发展,基于人工智能的图像去噪与图像增强算法已经取得了一系列令人瞩目的成果。
本文将对基于人工智能的图像去噪与图像增强算法进行研究,探讨其原理、方法和应用。
一、图像去噪算法研究在实际应用中,图像中常常受到各种噪声的影响,如高斯噪声、椒盐噪声等。
图像去噪的目标是通过算法将图像中的噪声去除,使得图像更加清晰和可识别。
1. 基于深度学习的图像去噪算法深度学习是人工智能领域的热门技术之一,其强大的特征提取和学习能力使得其在图像去噪领域取得了突破性的进展。
深度学习图像去噪算法可以学习到图像中的潜在噪声分布,并通过神经网络实现去噪的过程。
2. 基于小波变换的图像去噪算法小波变换是一种时间-频率分析方法,在图像处理领域具有很大的应用潜力。
基于小波变换的图像去噪算法通过将图像转换到小波域中,并对小波系数进行滤波处理,实现去除图像中的噪声。
3. 基于稀疏表示的图像去噪算法稀疏表示是一种数学工具,广泛应用于图像信号处理领域。
基于稀疏表示的图像去噪算法通过对图像进行稀疏表示,利用一个稀疏的表示矩阵来恢复原始图像。
二、图像增强算法研究图像增强是指通过一系列的算法和技术手段,对图像进行处理,改善图像的质量和视觉效果,使其更加清晰、具有更多细节和更好的对比度。
1. 基于深度学习的图像增强算法深度学习在图像增强领域也发挥着重要的作用。
基于深度学习的图像增强算法往往采用卷积神经网络结构,通过学习图像的特征,并通过非线性映射函数实现图像的增强。
2. 基于直方图均衡的图像增强算法直方图均衡是一种经典的图像增强算法,通过调整图像的灰度分布,使得图像的对比度更加均衡。
该算法简单高效,适用于多种图像增强场景。
3. 基于边缘增强的图像增强算法边缘是图像中的重要特征之一,基于边缘的图像增强算法旨在提高图像的边缘信息,使得图像在各种分析任务中能够更好地展现和突出边缘特征。
图像增强算法综述①靳阳阳, 韩现伟, 周书宁, 张世超(河南大学 物理与电子学院, 开封 475001)通讯作者: 韩现伟摘 要: 图像增强算法主要是对成像设备采集的图像进行一系列的加工处理, 增强图像的整体效果或是局部细节,从而提高整体与部分的对比度, 抑制不必要的细节信息, 改善图像的质量, 使其符合人眼的视觉特性. 首先, 本文从图像增强算法的基本原理出发, 归纳了直方图均衡图像增强、小波变换图像增强、偏微分方程图像增强、分数阶微分的图像增强、基于Retinex 理论的图像增强和基于深度学习的图像增强算法, 并讨论了它们的改进算法. 然后,从视觉效果、对比度、信息熵等方面对几种算法进行了定性和定量的对比, 分析了它们的优势和劣势. 最后, 对图像增强算法的未来发展趋势作了简单的展望.关键词: 图像增强; 直方图均衡; 小波变换; 微分方程; Retinex 理论; 深度学习引用格式: 靳阳阳,韩现伟,周书宁,张世超.图像增强算法综述.计算机系统应用,2021,30(6):18–27. /1003-3254/7956.htmlReview on Image Enhancement AlgorithmsJIN Yang-Yang, HAN Xian-Wei, ZHOU Shu-Ning, ZHANG Shi-Chao(School of Physics and Electronics, Henan University, Kaifeng 475001, China)Abstract : Image enhancement algorithm mainly process the captured images to enhance the overall effect or local details,increasing the overall and partial contrast while suppressing unwanted details. As a result, the quality of the images is improved, conforming to the visual perception of the human eye. Firstly, according to the basic principles of image enhancement algorithms, this study analyzes those based on histogram equalization, wavelet transform, partial differential equations, fractional-order differential equations, the Retinex theory and deep learning, and their improved algorithms.Then, the qualitative and quantitative comparisons between image enhancement algorithms are carried out with regard to visual effect, contrast, and information entropy to indentify the advantages and disadvantages of them. Finally, the future development trend of image enhancement algorithms is briefly predicted.Key words : image enhancement; histogram equalization; wavelet transform; differential equation; Retinex theory; deep learning在全球信息化大幅发展的时代, 对于这个世界的认识越来越依靠于信息的爆炸性传递. 大部分人认识世界的主要途径还是眼睛的可视性, 人眼所看到的一切都可以化作图像的形式. 图像的获取、生成、压缩、存储、变换过程自然会受到各种状况的影响, 例如获取图像时会因为天气原因, 不同光照条件, 图像亮度也有着细微的变化, 同样由于仪器设备的质量, 参数的设置, 人员的操作都会使图像质量在一定程度上的损伤, 影响图像的质量. 图像增强算法的出现, 无疑是对受损的图像做一个“修补”的工作, 以此来满足各样的需求. 图像增强的目的是为了适应人眼的视觉特性,且易于让机器来进行识别. 近些年来, 图像增强的发展计算机系统应用 ISSN 1003-3254, CODEN CSAOBNE-mail: Computer Systems & Applications,2021,30(6):18−27 [doi: 10.15888/ki.csa.007956] ©中国科学院软件研究所版权所有.Tel: +86-10-62661041① 收稿时间: 2020-10-12; 修改时间: 2020-11-05; 采用时间: 2020-11-17; csa 在线出版时间: 2021-06-01涉及了很多领域, 其中包括了遥感卫星成像领域、医学影像领域、影视摄影等各领域[1].要想真正地实现图像增强的效果, 首先对于整个图像来讲, 要提高图像部分和整体的对比度, 细节也不能忽略; 其次应提高图像的信噪比, 抑制噪声的产生,对“降质”的图像处理; 然后是对于增强过的图像来讲,避免出现局部增强不适, 影响人眼的观看模式.下面我们将列出几类典型的且应用范围比较广的图像增强算法以及改进的算法. 直方图均衡(HE)技术原理是对原图像的灰度直方图从比较集中的某个灰度区间转换为全部灰度区域内的均匀分布[2]; 由此算法进行转化的局部直方图均衡化[3], 符合图像局部特性; Kim 等提出的保持亮度的双直方图均衡算法(BBHE)[4],最大亮度双直方图均衡(MMBEBHE)算法有效地保持图像亮度[5]; 迭代阈值的双直方图均衡算法(IBBHE)[6]用迭代的方法达到增强对比度和亮度保持的效果; 彩色图像直方图均衡算法[7], 运算复杂度很低, 合并图像的视觉效果很好. 基于偏微分方程(PDE)的增强方法是把图像作为水平集或高维空间中的曲面, 再根据曲线和曲面演化逐步来增强图像的对比度[8]; 基于全变分模型插值的图像增强方法[9], 保留原图像的细节, 提高了对比度; 基于HE的偏微分方程增强方法, 在梯度域增强对比度基础上[10]提出新梯度变换函数. 小波变换中增强本质是图像信号分解为不同频段图像分量[11]; 小波变换图像多聚集模糊增强方法[12], 增强后的图像较为清晰; 基于离散余弦变换(DCT)和离散小波变换(DWT)的图像增强方法, 提高图像的质量, 同时减少计算复杂度和内存使用量[13]; 基于小波分析和伪彩色处理的图像增强方法[14], 在降噪增强的同时进一步提高图像分辨率. 基于量子力学偏微分方程的缺陷图像增强的研究[15]. 基于PDE的红外图像增强, 很好改进了传统对比度增强方法的不足[16]; 基于PDE平滑技术是一种新兴的图像增强滤波技术, 实质性、开创性的研究在图像增强滤波中引入的尺度空间理论[17]. 基于LBPV (Local Binary Pattern Variance)的分数阶微分图像增强算法[18],在图像纹理和细节方面处理效果比现有分数阶算法效果更好; 自适应分数阶微分理论指纹图像增强算法改进了传统分数阶微分形式, 提高了计算精度[19]. 基于多尺度Retinex的HSV彩色快速图像增强算法, 在HSV 颜色模型中有与Multi-Scale Retinex (MSR)等同的结果, 处理时间短[20]; 基于多尺度Retinex的数字射线照相增强算法, 改善对比度, 抑制噪声[21]; MSR与颜色恢复(MSRCR)算法增强的图像在复杂的情况下进行识别物体[22]; 基于变分Retinex方法的图像增强, 良好结合了MSRCR和变分方法的优点, 保证图像自然度[23].近年来, 基于深度学习的图像处理算法迎来了一个新的时代[24]. Hu等利用超分辨卷积神经网络(SRCNN)方法提高了风云卫星亮温图像的峰值信噪比, 结果较传统方法更精细[25]; Li等利用深度学习来增强低光图像, 提出利用深度的卷积神经网络进行学习, 提高图像质量[26].1 图像增强算法的介绍1.1 直方图均衡算法直方图均衡化算法, 简言之就是对图像直方图的每个灰度级来进行统计[3]. 实现归一化的处理, 再对每一灰度值求累积分布的结果, 可求得它的灰度映射表,由灰度映射表, 可对原始图像中的对应像素来进行修正, 生成一个修正后的图像.1.1.1 传统标准直方图均衡算法f HE传统直方图均衡算法是通过图像灰度级的映射,在变换函数作用下, 呈现出相对均匀分布的输出图像灰度级, 增强了图像的对比度. 该算法是相对于图1中n=1, 均衡函数为的简化模型[27], 即:f HEX k= {X0,X1,···,X L−1}其中, 函数代表直方图均衡过程, 其大致过程为: 已知输入和输出图像为X和Y, 总灰度级为L, 则存在, 均衡后输出和输入图之间有如下变换关系:c(X k)其中, 展现的累积概率分布表示函数输入图像灰度级.图1 全局均衡算法的模型L=∞如果输入图像看作一个连续随机变量, 即,则输出图像自然是一个随机变量, 输出图像灰度级均衡后的概率分布将趋于均匀, 则输出图像的亮度均值为:2021 年 第 30 卷 第 6 期计算机系统应用得到均衡后图像的均值分布与原图像无关, 由此可知其不能有效保持原始图像的亮度, 由于原图像各灰度级概率密度的差异, 简并现象的产生明显变多.1.1.2 保持亮度的双直方图均衡算法BBHE 实质是利用两个独立的子图像的直方图等价性[4]. 两个子图像的直方图等价性是根据输入图像的均值对其进行分解得到, 其约束条件是得到均衡化后的子图像在输入均值附近彼此有界作为基于图像均值进行的分割, 均衡后图像均值偏离原始图像均值的现象不会出现, 达到了亮度保持的目的, 其算法流程如下:G mean 1)计算输入图像均值, 根据均值将原始直方图分为左右两个子直方图.P L (i )P R (i )2)分别计算左右两个子直方图的灰度分布概率直方图和, 即:N L N R 其中, 和分别表示左右两个子直方图的总像素数,L 表示图像总灰度级数.cd f L (i )cd f R (i )3)计算左右两子直方图的累积分布直方图和, 即:tab L (i )tab R (i )4)计算左右两个映射表和, 合并之后得到最终的映射表tab , 其中round 表示四舍五入取整, 即:对于一些低照度和高亮的图像, 均值会处于较低和较高的地方, 若此时基于均值进行分割并分别均衡的话, 很大程度上会导致一个有大量数据的子直方图在小范围内进行均衡的情况出现, 另一个只有少量数据的子直方图却在较宽的范围内均衡.1.2 小波变换图像增强算法19世纪80年代Morlet 提出小波变换的概念, 数学家Merey 在十几年后提出小波基构造思想, 随着Mallat 的加入, 两个人共同建立了小波变换算法. 通过小波逆变换将同态滤波处理的低频分量和经自应阈值噪、改进模糊增强的高频分量得到增强处理后的红外图像[28].1.2.1 标准小波变换图像增强小波理论具有低熵和多分辨率的性质, 处理小波系数对降噪有一定作用, 噪声主要在高通系数中呈现,对高低通子带均需要增强对比度和去噪处理. 标准小波变换图像增强(WT)将图像分解为1个低通子图像和3个具有方向性的高通子图像, 高通子图像包括水平细节图像、垂直细节图像和对角细节图像[29]. 小波变换最大的特点是能较好地用频率表示某些特征的局部特征, 而且小波变换的尺度可以不同[30].1.2.2 改进后的小波变换图像增强算法针对传统方法对图像多聚焦模糊特征进行增强会出现图像不清晰、细节丢失现象, 小波变换图像多聚焦模糊特征增强方法, 利用背景差分法将目标图像的前景区域提取出来, 背景区域亮度会随时间发生变化,进而完成背景区域特征更新; 根据全局像素点熵值和预设阈值校正加强模糊特征, 突出小波变换图像边界局部纹理细节信息, 完成增强变换. 基于小波变换域的医学图像增强方法[31], 是基于Shearlet 变换改进的Gamma 校正, 采用改进的伽玛校正对低频进行处理, 利用模糊对比函数增强图像细节, 增强图像的对比度.二进小波变换简单的对信号尺度参数实现了离散化, 不过仍具备和连续小波变换同样的平移不变特性.利用二进小波变换将指纹图像分解[32], 步骤如下:1)首先将获取的指纹图像进行尺度的分解, 这样得到的频率分量为一低三高;2)对低频分量进行直方图均衡;3)对3个高频分量先进行高斯拉普拉斯掩膜锐化, 得到锐化后的图像;4)直方图均衡后的低频分量和处理后的3个高频分量进行二进小波逆变换重构, 得到增强后的图像.1.3 偏微分方程图像增强算法u (x 1,x 2,···,x n )关于未知函数的偏微分方程是形如式(11)的等式:计算机系统应用2021 年 第 30 卷 第 6 期x =(x 1,x 2,···,x n )Du =u x 1,u x 2,···,u x n 其中, , , F 是关于x 和未知函数u 加上u 的有限多个偏导数的基础函数. 偏微分方程(Partial Differential Equation, PDE)是微分方程的一种, 如果一个微分方程出现多元函数的偏导数, 这种方程就是偏微分方程[33].1.3.1 标准偏微分方程图像增强V l o (p )V l (p )l o V l o (p )V l (p )l o l o l o 假设和分别为两幅图像和l 的对比度场, 若与在每一点上具有相同的梯度方向,但前者大小均大于后者, 则图像应该比l 具有更高的对比度, 可以将看作l 的增强图像. 实际上, 从图像l 到图像的过程就是标准PDE 图像增强实现的过程,可以由以下式子来描述它们的关系:V l o (p )式中, 为增强后图像的对比度场; k 为增强因子,一般情况下k >1, 过大的话会增大噪声. 对于式(12),图像l 是已知的, 其解为:φl o (p )式中, 是一个与坐标无关的常数. 可看到两幅图像之间的动态范围存在k 倍的差距. 对于可在计算机屏幕上显示的数字图像, 其动态范围为0 ~ 255. 我们要做到先要对的对比度场进行约束, 之后开始按照步骤运算, 最后才能得到比较准确的数据.1.3.2 改进的偏微分方程增强方法∇u max ∥∇u ∥min为避免增强图像梯度场同时造成噪声的危害加剧,寻找一种比较适合的增强方法. 定义原图像的数值梯度函数为, 梯度模的最大值为, 最小值为, 增强之后的图像梯度为S [10]:∥∇u ∥[min ∥∇u ∥,max ∥∇u ∥][0,max ∥∇u ∥]式中, 表示梯度场的方向信息. 经过改进的梯度函数梯度场从的区域内映射到内. 原本纹理突显出来的同时保留梯度值较大的边缘.基于量子力学偏微分方程的缺陷图像增强研究方法[15]. 航空材料缺陷的图像增强对缺陷的定性和定量性能起着至关重要的作用, 由于复合材料分布不均匀,将导致缺陷成像对比度不高, 会让识别和量化的难度加大. 算法主要分为两个步骤: 首先是根据量子力学理论, 计算图像边缘的量子概率; 在此基础上, 建立融合各向异性量子概率的偏微分方程来增强航空材料缺陷图像. 此算法可以在有效抑制噪声和减少成像不均匀性的同时, 更好保留缺陷的特征, 增强图像的对比度.1.4 分数阶微分方程增强算法近些年, 分数阶微积分在多领域都有了突破性进展[34]. 分数阶微分不仅可以提升图像中的高频分量, 还可以以一种非线性形式保留图像中低频分量所带有的性能. 常用的分数阶微分定义有G-L 、R-L 、Caputo 三种定义, 其中最常用的是采用非整型分数阶微积分的G-L 定义[35].1.4.1 图像增强的分数阶微分算子构造m ×n 让图像像素邻域中任一像素与对应系数进行乘法运算, 得到的结果再进行和运算, 得到像素点所在位置的回复, 当邻域的大小为, 要求的系数会很多. 这些系数被排列成一个矩阵, 称为滤波器、模板或者掩模[36].f (x ,y )在整数阶微分方程的增强算子中, 有一类是拉普拉斯算子, 对任一二元连续函数来讲, 其拉氏变换可表示为:f (x ,y )f (x ,y )f (x ,y )x ∈[x 1,x 2]y∈[y 1,y 2]n x =[x 2−x 1]n y =[y 2−y 1]由于在图像中, 两个相邻像素点之间灰度产生差异的距离最小, 因此图像在它的x 和y 方向上灰度值的变化只能以像素之间的最小距离为单位来进行数值度量和分析, 所以的最小等分间隔只能设为: h =1, 如果图像中x 和y 方向的持续区间分别为和, 则最大等分份数分别为和.将上式拉普拉斯变换写成离散的表示形式, 对x 方向和y 方向重新定义, 得到它的二阶微分表示:根据以上定义, 可以得到:拉氏算子还要对处理前后的图像完成进一步的叠加, 其方式如下:2021 年 第 30 卷 第 6 期计算机系统应用在雾天图像中应用算子增强图像, 边缘轮廓还有纹理部分的效果会很容易看到, 不过若是图像像素中某一范围灰度变化不明显, 细节可能受到损失. 因此,构建图像增强的分数阶微分算子, 将整数阶微分扩展到分数阶微分上并且应用于图像增强中[37].1.4.2 改进的分数阶微分算子增强图像相比传统的分数阶微分算法的不足, 提出新的改进算法, 在极端条件处理拍摄的交通图像时, 具有良好效果. 上文提到的指纹图像增强算法, 对传统形式加以改造, 在计算精度上有所提升, 进而构造了更加高精度的分数阶微分掩模. 通过对像素周围的纹理对比从而逐点选择微分阶, 明确的选择了具有二阶精度的分数阶微分形式来构造IRH 算子, 并对算子结构进行相应的改进, 之后利用图像的梯度信息和局部统计信息, 结合中心像素对相邻像素的影响, 建立自适应分数阶微分的自适应函数, 此法保留了指纹纹线和图像纹理细节, 对于降噪起到很好的作用.1.5 Retinex 图像增强算法S (x ,y )L (x ,y )R (x ,y )S (x ,y )L (x ,y )Retinex 是retina(视网膜)和cortexv(大脑皮层)组成的, Retinex 算法由美国物理学家提出[38]. Retinex 理论的基础是人类视觉系统的色彩恒常性, 人类视觉感知系统的色知觉存在“先入为主”的特性, 即光源条件发生改变, 视网膜接收到的彩色信息也会被人们的大脑驳回. Retinex 理论的依据就是是原始图像可以分解为照射图像和反射图像, 最重要的就是让摆脱的影响, 以便得到图像的反射属性.1.5.1 经典的Retinex 图像增强对数域进行操作可以把乘法运算变成简单的加法运算, 进而出现了多种Retinex 算法. 经典的有: 单尺度Retinex 算法(SSR)、多尺度Retinex 算法(MSR)和带色彩恢复的多尺度Retinex 算法(MSMCR)等[39].针对运算速度缓慢的问题, 在1986年, Jobson 等[40]将高斯低通滤波与Retinex 结合, 改进了Land 提出的中心环绕Retinex 算法(Center/Surround Retinex), 提出了单尺度Retinex(SSR)算法. 在SSR 算法中, Jobson 等创新的使用高斯函数与图像进行卷积的方式来近似实现了入射分量的表达. 它的数学表达式如式(20)表示:I i (x ,y )i ∈(R ,G ,B )G (x ,y ,c )∗L i (x ,y )其中, 表示原始图像的第i 个通道分量的像素值,颜色通道中的一个, 表示中心环绕函数, 是一种卷积操作表示, 入射分量的表达可以借用Jobson 等的成果, 则可以看做入射图像的第i 个通道分量. SSR 的实现过程如式(21)至式(23)所示:由于SSR 算法处理要对图像细节对比度和色彩的保留做到很好的发展, 而尺度c 又相对难做到极好的运用, MSR 算法的出现, 在很大程度上解决了这一问题, 起到了平衡图像色彩和细节的良好效果.1.5.2 改进的Retinex 图像增强Retinex 算法对于图像增强的效果需要经过精确且复杂的计算, 最后的结果精确度越高, 增强效果将会更好. 文献[20]中基于多尺度Retinex 的HSV 彩色快速图像增强算法. 在HSV 模型中用多尺度Retinex 进行图像增强, 由于颜色转换的非线性, 计算起来非常复杂. 使用亮度校正的MSR 算法基于HSV 颜色模型和修正的V 频道输出图像的RGB 分量的线性形式减少30–75%的平均处理时间, MSR 算法在Haar 小波变换低频区域应用亮度校正的处理速度有很明显优势, 平均加速度接近3倍. 文献[22,23]中介绍了MSRCR 算法. 由于传统均值移位算法有不少的不足, 改进后, 对要增强的图像可以在情况复杂下进行识别物体, 增强对比度的同时, 光晕现象的产生被消灭, 噪声得到抑制,保证图像自然度. 基于Retinex 提出一种自适应的图像增强方法, 其中包括如下4个步骤: (1)用引导滤波器估计其照度分量; (2)提取图像的反射分量; (3)对反射分量进行颜色恢复校正; (4)后处理. 由于雾霾和照度较低, 自然生成的图像质量比较差, 而此法不管是在定量还是定性上都突出了更好的优势. 此算法最终的结果图像具有清晰的对比度和生动自然的颜色[41].1.6 基于深度学习的图像增强算法在当今社会经济科技奋进之时, 深度学习的发展可谓是如日中天, 特别是在图像增强方面.1.6.1 卷积神经网络图像增强算法神经网络(neural networks)最基本的组成结构是计算机系统应用2021 年 第 30 卷 第 6 期神经元(neuron), 神经元概念源于生物神经网络[42]. 卷积神经网络(Convolutional Neural Network, CNN)在传统神经网络基础上, 引入了卷积(convolution)和池化(pooling), CNN 的建筑灵感来自于视觉感知[43]. CNN 是深度学习领域最重要的网络之一, CNN 在计算机视觉和自然语言处理等诸领域都有很大成就. 卷积神经网络的特性比较突出, 除了可以实现权值共享外, 可调的参数相对来说不多, 对二维图像这类的, 它的平移、倾斜、缩放包括其他形变都拥有着极高的不变性.CNN 相比于一般的神经网络, 具有很大优势[44]: (1)局部连接. 每个神经元只与少数神经元相连, 有效地减少了参数, 加快了收敛速度; (2)重量共享. 一组连接可做到同时分享相同的权值, 进一步降低了所需的参数;(3)降采样降维. 池化层利用图像部分相关的依据对图像进行降采样, 降低运算数据量, 留存有效的信息值.卷积神经网络大致包含4部分, 卷积层、池化层、全连接层以及反卷积层, 各自具有不同作用, 承担独自的工作. 深度越深, 网络性能越好; 随着深度增加, 网络性能逐渐饱和.1.6.2 基于深度学习图像增强的改进算法f o=F (g )F (g )Hu 等基于深度学习方法增强MMSI 亮温图像, 设计卷积神经网络重建风云四号卫星MMSI 的亮温图像和风云三号卫星微波成像仪亮温图像[25]. 在根据SRCNN进行实现映射函数, 式中, g 为监测的天线温度的图像, 可用于复原, 使其尽可能接近地面真实高分辨率亮温图像f . 映射函数F 的完成可以依据学习思想, 构建一种卷积神经网络, 为了让观测图像数据重新构建为理想的高分辨数据, 需要对卷积神经网络进行一系列特征变换, 此过程即达成卷积核的卷积操作.相比古老的插值方法而言, SRCNN 方法除了提高图像的峰值信噪比之外, 在提高图像细节较古老的方法也有很大的提高.2 图像增强算法的评价和对比2.1 各种算法增强效果的分析通过对论文文献研究比对, 以及对于其中的经典算法以及改进的算法, 对应用广泛的上述6大类图像增强算法进行较概括的研究分析.图2是几种不同算法得到的增强图像. 从增强图像的效果来看, HE 增强效果是对图像的动态范围进行拉大实现的, 增强效果随动态范围增加而变差. BBHE算法均衡后的图像在增强对比度的同时很好保持原图像的平均亮度. IBBHE 根据各子图像的直方图分别进行独立的均衡化处理, IBBHE 增强效果更好. WT 算法增强图像细节信息, 但是增加了噪声. 小波变换图像多聚集模糊增强方法, 对图像增强后, 图像较为清晰, 细节没有丢失, 效果较好. PDE 和TVPDE 算法放大了图像对比度场, 增强后图像都有较高对比度[45]. 自适应分数阶微分可以很好降噪. SSR 和MSR 算法去除了图像中照度分量影响, 还原景物本身的亮度信息, MSRCR 处理后的图像比原图像细节增加了, 亮度有所提高, 颜色有一定矫正, 对颜色的恢复存在失真现象. 基于深度学习的图像增强算法通过复杂的神经网络, 进行大量的训练, 得到的模型同时减少了训练时间, 取得了更好的精度.2.2 算法增强效果的对比对一幅图像的增强效果来讲, 需要对图像对比度和信息熵来进行评价和比较, 可以对图像有很好认识.图像对比度的计算公式:I i ,j 其中, 为中心像素点的灰度值, N 为图像局部块内像素点的个数. 为了计算一幅完整图像的对比度, 需要对图像中所有部分块对比度总体的平均值来表示.图像的信息熵公式如下:p (k )式中, 为灰度级k 的概率密度, M 为最大的灰度级.表1中为第一幅图通过不同算法得到的图像质量的客观结果评价, 评价指标为对比度和信息熵. 通过对文献中算法的研究以及本文中对增强算法的分析对比, 我们得到表2中对不同算法优缺点的总结.3 增强算法发展趋势及有意义的研究方向根据上文所介绍的不同图像增强算法及实验分析对比结果, 可预见未来的图像增强算法发展将有以下特点: 超分辨率、多维化、智能化和超高速.1)超分辨率, 对获得的低分辨率图像进行增强从而得到超高分辨率的图像, 重点是对采集分辨率以及显示分辨率做进一步的提升, 突破技术壁垒限制, 向时空感知超分辨率迈进.2021 年 第 30 卷 第 6 期计算机系统应用。
基于PDE方程的图像处理技术研究一、引言图像处理技术是数字图像处理领域的重要研究方向之一。
随着计算机科学和技术的快速发展,基于PDE(偏微分方程)的图像处理技术在图像降噪、边缘检测、图像增强和图像分割等方面取得了显著的成果。
本文将系统地探讨基于PDE方程的图像处理技术在不同应用领域的研究现状和进展。
二、图像降噪图像降噪是图像处理领域的重要研究课题之一。
基于PDE方程的图像降噪技术能够通过去除图像中的噪声,提高图像的质量和清晰度。
常见的PDE方程图像降噪方法有扩散滤波模型、非线性扩散滤波模型和全变分模型等。
其中,扩散滤波模型通过图像的梯度流来消除噪声,而非线性扩散滤波模型能够更好地保持图像的边缘信息。
全变分模型则通过最小化图像的总变差来实现降噪效果。
三、边缘检测图像中的边缘信息对于图像分析和理解至关重要。
基于PDE方程的图像边缘检测技术能够提取图像中的边缘信息,帮助用户更好地理解图像内容。
常用的PDE方程边缘检测方法有Canny算子、Marr-Hildreth算子和基于水平集的方法等。
Canny算子通过寻找图像中灰度值变化最大的位置来检测边缘。
Marr-Hildreth算子则利用高斯滤波和拉普拉斯算子来提取图像中的边缘。
基于水平集的方法通过优化演化方程来检测图像中的边缘。
四、图像增强图像增强是指通过一系列的处理手段,改善图像的视觉效果和质量。
基于PDE方程的图像增强技术能够增强图像的对比度和细节,使得图像更加清晰和易于理解。
常用的PDE方程图像增强方法有Retinex算法、双曲正切函数和直方图均衡化等。
Retinex算法通过模拟人眼的调节特性,提高图像的动态范围和对比度。
双曲正切函数则通过变换函数来显著增强图像的细节。
直方图均衡化则通过调整图像像素的统计分布来增强图像的对比度。
五、图像分割图像分割是将图像分成若干个不同区域的过程,能够帮助用户更好地理解图像内容和提取感兴趣的目标。
基于PDE方程的图像分割技术能够通过优化演化方程,将图像像素划分到不同的区域中。
华中科技大学博士学位论文基于小波及其统计特性的图像去噪方法研究姓名:侯建华申请学位级别:博士专业:控制科学与工程指导教师:田金文20070601华中科技大学博士学位论文摘要图像在获取或传输过程中不可避免地会受到噪声污染, 图像中的噪声严重影响了后续的图像处理工作,如图像分割、编码、特征提取和目标检测等。
为了提高图像的质量以及后续更高层次处理的需要, 对图像进行去噪就成为图像预处理中一项非常重要的工作。
图像去噪的目的就是从被噪声污染的含噪图像中恢复出原始的“干净”图像, 即在滤除噪声的同时尽可能的保留重要的图像特征与细节。
传统图像去噪方法在降噪与保细节折中方面难以令人满意;小波变换作为一种新的时频分析方法,具有多尺度、多分辨分析的特点,为信号处理提供了一种新的、强有力的手段,在图像去噪领域得到了成功的应用。
目前, 基于小波的去噪方法已经成为图像去噪和恢复的重大分支,而根据图像小波系数的统计特性,研究基于模型的去噪方法,是目前图像去噪领域中的主要研究方向,无论在理论上还是在实际应用中都具有重要意义。
本文以小波分析理论为工具, 对小波域图像去噪理论与方法进行了系统、深入的研究,主要工作包括以下四部分:1、小波图像去噪方法研究综述本文前两章作为全文的基础,对基于小波的图像去噪方法进行了全面的研究总结。
首先综述了图像去噪技术的发展现状,特别是小波图像去噪方法的研究进展。
针对目前小波图像去噪领域尚未有一个较全面的分类方法, 本文以该领域发展的三个阶段为线索, 将小波图像去噪算法进行了新的分类并划分为四类, 并对每种类型中代表性的算法做了分析讨论。
阈值去噪是小波去噪研究中一类非常重要的方法, 对此进行了系统、深入的分析,在阈值选择这一核心问题上,对最具代表性的阈值结合具体的算法在原理与方法上做了清晰的阐述, 并对这些典型算法分别在正交小波变换基和平移不变小波基下进行了全面的实验仿真和分析讨论。
从实验结果和性能分析中得到了一些有意义的结论。
基于PDE的图像降噪方法的研究综述
摘要:回顾了三种传统的PDE图像降噪方法,即j函数、基于散度的PDE 和基于拉普拉斯算子的扩散PDE这三种图像降噪方法,并指出可将扩散方程与自由薛定锷方程结合形成新的基于复扩散的图像降噪处理方法。
关键词:函数;拉普拉斯算子;PDE;复扩散
1引言
实际的图像常受一些误差的影响而发生退化,这种退化就称之为“噪声”,可将噪声视为“妨碍人们感觉器官对所接收的信源信息理解的因素”[1~3]。
噪声一般由其概率特征来描述。
根据处理的目的一般有两种不同的降噪方法:一是平滑,其目的在于抑制噪声或其他小的波动,但平滑也会使含有图像重要信息的明显边缘变模糊。
二是梯度算子,梯度算子的作用是在图像中显现这些位置。
但如果在图像中使用梯度算子,也会同时抬高噪声水平。
很明显,平滑和梯度算子的目标相互抵触,因此需要寻找解决上述问题的方法。
自上个世纪80年代末,偏微分方程(PDE)开始用于解决传统的计算机视觉问题,像图像降噪问题等。
PDE可以描述高度非线性和迭代性的连续对象的演变如曲线、曲面或者向量场等。
用PDE表示的曲面、曲线或者向量场来模型化并进行演变已成为一种自然的图像处理方式。
2 图像降噪处理及变分PDE方法
2.1 图像降噪处理
可将图像降噪处理视为一个输入—输出系统,即是
噪声图像或者模糊的图像。
I是降噪和增强后的图像。
是图像降噪处理器,可以表示为减少模糊和降噪。
u0 = Ku+n,u0是原始图像,K是模糊核,n是加性噪声。
u是降噪和增强后的图像。
由此可见,图像降噪处理解决两个问题:一是如何模型化和表示原始图像
I0;二是模型化图像降噪处理器。
原始图像模型化的方法将影响图像降噪处理的方法。
对原始图像模型化和表示的方法有很多种,应用PDE方法进行图像处理时,往往基于Mumford和Shah的对象边缘模型和Rudin、Osher和Fatemi的BV模型。
其中BV模型[4]常用于数字图像处理。
该模型假设理想的图像是有界的总体变分
,因此,图像的降噪处理器可以用变分优化的方法加以实现。
例如:通
过计算下面的公式来获得:满足
其中,认为白噪声是高斯函数的近似。
该方程可以用变分方法加以
解决。
2.2 变分PDE方法
变分计算是一种寻找偏微分方程(该方程与函数积分有关)极值的有效方法,如:
(1)
可以通过Euler-Lagrange方程来求解:
(2)
方程(2)可以使用梯度减小的迭代方法来求解。
实际上,可以认为方程(2)是方程E(I)的梯度。
从初始条件I0开始,沿着梯度的相反方向将会求得E(I)的最小值Imin:
(3)
一般情况下,考虑到方程E(I)是非凸的,I0的选择很重要,理想的情况下,选择I0接近于方程E(I)的全局最小值。
不同的I0将得到不同的结果(不同的局部最小值)。
注意到,这种PDE演变用时间变量t进行参数化。
它描述了函数
I获得E(I)最小值的连续变化过程。
当获得最小值时,。
Euler-Lagrange方程将PDE演变与梯度下降联系起来。
一般情况下,更关注于梯度下降而不是方程的最小值。
变分PDE方法与其它方法相比,无论在理论和计算上都具有明显的优势:一是可以直接处理重要的几何特征如梯度、正切、曲线和level set,也可以有效地模拟可视化意义上的动态过程如线性、非线性扩散和信息传输机制;二是在计算方面,可以从大量已存在的关于PDE数值分析和计算的文献中改进计算的效率。
3 PDE图像降噪方法
早期主要应用线性扩散PDE方法进行图像处理,但这种方法的局限性在于对局部图像特征和噪声的一致滤波。
因此Persona和Malik提出非线性扩散过程,这种扩散是通过变量的扩散来控制平滑的效果。
3.1 P-M方法
P-M方法的扩散系数是图像梯度的减函数。
这种操作有选择地对不含大梯度的区域进行低通滤波。
P-M方法使用如下形式的非线性扩散方程进行图像处理[5]:
(4)
其中是梯度的减函数,在接近边缘的地方(高梯度)接近于0,在平滑区域(低梯度)接近于1。
根据“最小-最大”原理,在一维情况下任意时
刻都不会产生新的局部极小值和极大值。
当图像随着时间发生演变时,全局最小值(最大值)总是受输入数据的约束,并且是一个非减(非增)函数。
P-M方法可以满足这些条件,保证了 PDE的稳定性。
P-M方法的扩散系数。
K是区分平滑区域和边缘区域的固
定的梯度阈值。
3.2 P-M型方程的分解
为了更好地理解P-M型方程的扩散过程,可以将方程(4)分解为如下类型[6]
其中,如果使用上述P-M的系数,则
和是I在互相垂直方向上和的二阶导数。
,
在二维情况下,
,处处正切于等值线。
由分解可知,P-M的系数当,。
这就意味着:即使,在图像的某些点上也存在着逆扩散。
而逆扩散增强噪声的不稳定过程。
3.3 函数方法
上述P-M方法实际上可以用于求解如下的偏微分方程(5)的极小值。
(5)
该方程更一般的解可以描述如下:
(6)
不同的函数将导致不同的降噪效果。
目前,已提出了一些函数(如下表
所示)。
一般来说,函数应满足如下条件:
一是,
二是,
三是。
在上述函数中只有Minimal surfaces和Green的j函数才满足上述要求。
因此,这两个函数用于降噪的效果比其它函数效果要好一些。
3.4 基于散度的PDE降噪方法
为了能依据更准确的图像特征来表达,基于散度的PDE降噪方法将用其它的表达式来代替,这样就可以更方便地设计基于PDE的图像降噪处理器。
一般的处理方法可表述如下[7]:,是扩散张量。
其中,
,是的特征向量,,是它的特征值。
,
3.5 基于拉普拉斯算子的扩散PDE方法
基于拉普拉斯算子的扩散PDE方法可以表述如下:,其中
,,分别是的二阶导数。
其基本思想是考虑到函数方
法和基于散度的PDE方法对噪声较为敏感。
为了克服这个缺点,不同的学者提出了一些较为有效的算法。
一种常见的算法是将图像信号的二阶导数用低通滤波器作卷积运算。
如[8]
(7)
在很多情况下,上式的缺陷用中等波长的高斯函数作卷积运算不能消除噪声产生的拐点,而长波高斯函数可以消除噪声,但又会使拐点的位置不准确。
因此,基于对平滑区域降噪而对边缘进行增强的思想,Alvarez和Mazorra 提出了如下方程[9]:
(8)
其中c是正常数。
另外,Kornprobst提出了如下方程[10]:
(9)
其中,当时,
,否则为0。
Coulon和Arridge提出了如下方程[13]:
(10)
其中
4 新的PDE图像降噪方法
最近,Gilboa等将实数域中的分析扩展到复数域[12],并将扩散方程与自由薛定锷方程联系起来。
认为线性复扩散同时具有前向扩散和逆扩散的特征,可以在保留希望得到的数学和视觉特征的基础上获得稳定的流。
同时,从理论和数值计算角度来看,当复扩散系数接近于实轴时,虚部可以视为边缘探测算子。
他们提出了一些线性和非线性复扩散算法。
5 小结
针对现有的PDE图像去噪处理方法,分为j函数、基于散度的PDE和基于拉普拉斯算子的扩散PDE加以总结,并指出将分析扩展到复数域将是一个很有前景的基于PDE图像去噪处理的方向。