CN103337060A - 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器 - Google Patents
一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器 Download PDFInfo
- Publication number
- CN103337060A CN103337060A CN2013102983826A CN201310298382A CN103337060A CN 103337060 A CN103337060 A CN 103337060A CN 2013102983826 A CN2013102983826 A CN 2013102983826A CN 201310298382 A CN201310298382 A CN 201310298382A CN 103337060 A CN103337060 A CN 103337060A
- Authority
- CN
- China
- Prior art keywords
- finished
- calculating
- gamma
- fractional order
- multiplier
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明所提出的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器是基于一种特殊分数阶偏微分方程去噪算法来实现对纹理图像的分数阶、非线性、多尺度、高精度去噪。该滤波器是采用缓存器、差值平方器一、-v3次方幂方器、乘法器一、λn发生器、发生器、乘法器二、乘法器三和加法器一以级联方式构成的。该滤波器具有在去除纹理图像噪声的同时,既能尽量保留平滑区域中的低频轮廓,同时又能非线性保留灰度值跃变幅度相对较大的高频边缘,而且还能非线性保留灰度值跃变幅度变化相对不大的纹理细节的显著优点。该滤波器特别适用于对富含复杂纹理细节特征的图像进行去噪的应用场合。本发明属于应用数学、数字图像处理和数字电路交叉学科的技术领域。
Description
所属领域
本发明所提出的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器是基于一种特殊分数阶偏微分方程去噪算法来实现对纹理图像的分数阶、非线性、多尺度、高精度去噪。本发明涉及的分数阶微积分的阶次v1、v2和v3不是传统的整数阶,而是非整数阶,工程应用中一般取分数或有理小数,且v3≠1,2,3。见图1,该滤波器是采用缓存器2、差值平方器一3、-v3次方幂方器4、乘法器一5、λn发生器11、发生器10、乘法器二8、乘法器三9和加法器一6以级联方式构成的。该滤波器具有在去除纹理图像噪声的同时,既能尽量保留平滑区域中的低频轮廓,同时又能非线性保留灰度值跃变幅度相对较大的高频边缘,而且还能非线性保留灰度值跃变幅度变化相对不大的纹理细节的显著优点。该滤波器特别适用于对富含复杂纹理细节特征的图像进行去噪的应用场合。本发明属于应用数学、数字图像处理和数字电路交叉学科的技术领域。
背景技术
数字图像处理理论主要包括三大类方法:随机建模、小波理论和偏微分方程方法。其中,基于偏微分方程的图像处理属于数学分析中重要的一部分,是图像处理领域中的一个重要分支。偏微分方程方法它与物理世界紧密联系在一起。著名的波动方程和热传导方程都属于整数阶偏微分方程,还有Euler方程、Poisson方程和Laplace方程等。物理学中的整数阶偏微分方程经常被应用到其他领域,如生物、金融等,并已被应用到了数字图像处理领域。关于基于整数阶偏微分方程的数字图像处理技术,一方面,该图像处理方法属于低层图像处理的范畴,其处理结果通常被当作中间结果提供给其他图像处理方法进一步使用;另一方面,随着该图像处理方法的深入研究,人们越来越深刻地挖掘图像和图像处理的本质,并试图用严格的数学理论对现存的传统图像处理方法进行改造,这对于以实用为主的传统图像处理方法是一种挑战。
目前,虽然偏微分方程已被应用到了数字图像处理领域,但是绝大多数相关研究都还仅仅局限于整数阶偏微分方程的应用,然而对于分数阶偏微分方程在数字图像处理领域中的应用在国内外都还研究甚少。整数阶偏微分方程本身来自连续域,所以它本质上能描述的是模拟图像,一旦其解的存在性和唯一性被证明了,我们就可以利用离散的数值方法对针对数字图像的整数阶偏微分方程求取其数值解。因为基于整数阶偏微分方程的数字图像处理可以与一些物理过程相联系,因此它们通常都用连续域进行描述。一般来说,整数阶偏微分方程方法与通常的滤波方法相比计算量是比较大的:需要迭代求解或者是有限差分所构造的方程组求解,整数阶偏微分方程类方法的主要优点为:第一,整数阶偏微分方程和相应的曲线(曲面)流给出分析图像的连续模型,离散的滤波表现为连续的微分算子,因而使得网格的划分、局部非线性分析易于实现。另一方面,当图像表示为连续信号,整数阶偏微分方程可以被视为在微小子邻域局部滤波的迭代,这种特性允许将已有的滤波方法进行合成与分类,并可形成新的滤波方法。第二,利用整数阶偏微分方程处理数字图像易于直接掌握和处理诸如梯度、切线、曲率以及水平集等视觉上重要的几何特征,还能有效地模拟诸如线性和非线性扩散以及信息传递机制那样的视觉上有意义的动力学过程。第三,整数阶偏微分方程领域的独特分析理论为研究更好的数字图像处理算法与有意义的理论结果,如解的存在性、唯一性等,提供了可能。特别地,最值得注意的优点在于整数阶偏微分方程方法能够获得较好的图像质量,并具有一定的稳定性。灵活多样的数值方案为图像处理方程的数值计算提供了较大的帮助。在数字图像处理领域比较有效的整数阶偏微分方程的引入可以追溯到上世纪80年代末期,在90年代得到了很长足的发展。该研究可以追溯到Nagao、Rudin等关于图像光滑和图像增强以及Koenderink对于图像结构的探索。多个经典的整数阶偏微分方程被应用到了数字图像处理当中,例如热传导方程、薛定谔方程、对流传导方程等。目前,基于整数阶偏微分方程的图像处理技术取得了一些较好的应用,例如法国宇航局已经采用了AMSS方法作为对航拍图像进行图像增强的标准方法。整数阶偏微分方程方法本身是物理学的内容,在数字图像处理中最早应用的可能要算各向同性介质中的热传导方程。若把灰度图像看成是一个各向同性介质中的温度场,那么这个温度场的热传导过程恰好对应着图像的高斯平滑过程,高斯滤波器方差参数与传导时间有关。但是由于高斯平滑是各向同性的,所以对于边缘的破坏作用很大,其应用也受到限制。1987年,Kass等利用图像边缘所需要的内部外部约束定义了一个表征轮廓曲线优劣的能量函数,其中内部约束主要考虑轮廓的光滑性和曲率,外部能量表示图像边缘轮廓的吸引。通过优化(最小化)这个能量函数,初始给定的轮廓可以收敛到邻近的图像边缘上。这种方法的物理意义明确,但是由于其考虑的对象(轮廓)是[]2的一维目标(假设是二维图像中的轮廓,若考虑三维图像,如三维医学图像,那么这个轮廓对应的是[]3中的一个二维曲面),它的描述方式和离散化均受到了一定的限制,并且其描述方式直接限制了轮廓曲线的拓扑变化,如分裂、合并等。1989年,Mumford和Shah提出了图像分割的变分模型。1992年,Chan和Vese利用模式识别中的类内距离最小的思想构造了无边缘的活动轮廓模型,之后Yezzi等同样利用类间距离最大的思想构造了一种新的活动轮廓模型。1995年,Osher等提出用水平集去描述一个与曲率有关的波前传播过程。这类方法的本质是将图像轮廓看成是一个二维函数的零水平集,那么通过研究这个二维函数的变化行为,就可以知道轮廓的变化方式。同时由于这个被研究的对象是一个二维函数,它在[]2中很容易描述和求解(相对于一维对象而言),并且一维轮廓不是直接求解的对象,所以通过二维函数的变化,使得轮廓的分裂、合并等拓扑变化的处理变得相对容易。自上世纪90年代后期起,整数阶偏微分方程开始应用到数字图像修复,即填充数字图像中丢失的部分、或移除数字图像中的障碍物等,以使得结果图像看起来像是真实的,它是图像编辑领域的一个很难的任务。基于整数阶偏微分方程的图像修复有两大方向,UCLA的Chan和Shen等利用能量优化来处理这个问题,主要是对结构图像边缘的一些性质(如简单性、曲率小等)作假设,然后构造相应的能量函数来描述,通过整数阶变分法来转化成整数阶偏微分方程进行求解;而以Bertalmio为代表的另一流派则直接考虑图像中某些性质的扩散过程,直接给出整数阶偏微分方程进行演化求偏微分方程解。这两类的方法都取得了较大的成功。另外在图像编辑领域,Poisson方程也在图像的无缝粘贴上占有主导地位。
在基于整数阶偏微分方程的图像处理中,图像去噪是其最重要的研究内容之一。基于整数阶偏微分方程的图像去噪分为两类:基于非线性扩散的方法和基于能量范函最小化的变分法。与之对应的两种基本模型是:由Perona和Malik提出的各项异性扩散(PM)模型以及Rudin,Osher和Fatemi提出的全变分(ROF)模型。PM模型使用热能的扩散过程来模拟图像的去噪过程,图像去噪的结果就是热能扩散达到平衡时的状态。用全变分来描述上述热能,就是ROF模型。在此基础之上,有学者分别将PM模型和ROF模型推广到彩色图像处理之中。有学者研究了模型中的参数选择,以及如何计算迭代求解过程的最优停止点。Rudin等人提出一种可变时间步长方法来解Euler-Lagrange方程。C.R.Vogel和M.E.Oman用不动点迭代方法来提高ROF模型的稳定性。D.C.Dobson和C.R.Vogel修改全变分形式来保证ROF模型数值计算的收敛性。A.Chambolle提出一种基于对偶公式的快速算法。J.Darbon和M.Sigelle利用水平集方法将原始问题分解为相互独立的马尔科夫随机场的优化问题,通过重建得到全局最优解。有学者提出一种迭代加权范数来求解全变分以提高计算效率。F.Catte等将原图像先经过一次高斯平滑,使PM模型具有适定性。PM模型和ROF模型都具有容易产生对比信息丢失,纹理信息丢失和阶梯效应等显著缺点。针对这些缺点,人们提出了许多改进模型。为了保持对比信息和纹理信息,有学者使用L1范数取代L2范数。S.Osher等提出一种迭代正则化方法。G.Gilboa,Y.Y.Zeevi和N.Sochen提出一种随空间变化的自适应数值保真项的方法。S.Esedoglu和S.Osher提出一种保持特定边缘的方向信息;为了消除阶梯效应,P.Blomgren提出一种全变分项随梯度变化的模型。有学者还将高阶导数引入能量范函中,或将高阶导数和原始ROF模型进行结合,或提出两阶段去噪等改进方法。上述基于整数阶偏微分方程的图像去噪改进方法,对于保持图像的对比信息和边缘信息,以及消除阶梯效应取得了一定的效果。
然而不幸的是,当我们直接将传统的基于整数阶偏微分方程的图像去噪方法应用于纹理图像去噪时,一般很难取得较好的处理效果。因为,一方面,传统的基于整数阶偏微分方程的图像去噪方法在本质上是基于整数阶微积分运算。它很难较好地处理一些非线性、非因果、非最小相位系统、非高斯、非平稳、非整数维(分形)信号和非白色的加性噪声等。如果我们直接将基于整数阶偏微分方程的图像去噪方法应用于纹理图像去噪时,会存在如下局限性:第一,整数阶偏微分方程仅根据局部信息处理图像,因此不能保持周期性的纹理特征,也不能恢复图像全局特征;第二,基于整数阶变分原理的整数阶偏微分方程模型,通过优化能量泛函实现图像处理。能量泛函实现的是局部邻域内的优化,因此处理后图像中存在块状效应;第三,仅包含前向或后向扩散的整数阶偏微分方程处理能力有限,而双向扩散的方程在扩散过程中会出现两个方向扩散信息抵消的现象,影响最终处理结果。另一方面,对于富含复杂纹理细节信息的纹理图像而言,图像的纹理细节信息对其判读的准确性显得极具价值。纹理图像去噪方法具有对比不变、纹理特征不变等特殊要求。由于常数或直流分量的整数阶微分值为零,细微波动的交流分量的整数阶微分值经过特定阈值限流后,其值亦为零,所以整数阶微分运算会对图像复杂纹理细节信息造成极大损失。当传统的基于整数阶偏微分方程的图像去噪方法直接应用于纹理图像去噪时,在低分辨率条件下,由于原始图像中的纹理细节信息本来就不够丰富和清楚,其处理结果还能够被勉强接受;然而,当分辨率逐渐增大,其处理结果图像中的复杂纹理细节信息将会被极大损失,致使我们很难对处理结果图像的纹理细节进行准确判读。因此,对于富含复杂纹理细节信息的纹理图像去噪而言,为了在去噪的同时更为有效地保持和利用复杂纹理细节信息,这就迫切要求我们提出一类能够分数阶、非线性、多尺度地处理图像复杂纹理细节特征的基于分数阶偏微分方程的纹理图像去噪方法。
发明内容
本发明所提出的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器是基于一种特殊分数阶偏微分方程去噪算法来实现对纹理图像的分数阶、非线性、多尺度、高精度去噪。本发明涉及的分数阶微积分的阶次v1、v2和v3不是传统的整数阶,而是非整数阶,工程应用中一般取分数或有理小数,且v3≠1,2,3。见图1,该滤波器是采用缓存器2、差值平方器一3、-v3次方幂方器4、乘法器一5、λn发生器11、发生器10、乘法器二8、乘法器三9和加法器一6以级联方式构成的。该滤波器具有在去除纹理图像噪声的同时,既能尽量保留平滑区域中的低频轮廓,同时又能非线性保留灰度值跃变幅度相对较大的高频边缘,而且还能非线性保留灰度值跃变幅度变化相对不大的纹理细节的显著优点。该滤波器特别适用于对富含复杂纹理细节特征的图像进行去噪的应用场合。
见图1,为了清楚说明本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的电路构成,有必要先对该基于分数阶偏微分方程的适于纹理图像去噪的高精度滤波器的数学公式推导和数值运算规则进行如下简要说明:
众所周知,分形数学理论产生了测度观的转变,分形几何否定了牛顿一莱布尼兹导数的存在性。以Hausdorff测度为基础的分形理论,虽然历经了90余年的研究至今仍然还是一种很不完善的数学理论。Hausdorff测度下的微积分数学理论的构造至今尚未能完成。目前发展比较成熟的是在欧氏测度下定义的分数阶微积分,它在数学上要求必须使用欧氏测度。在欧氏测度下,分数阶微积分最常用的是Grümwald-Letnikov定义和Riemann-Liouville定义两种。
Grümwald-Letnikov定义信号s(x)的v阶微积分为 其中,信号s(x)的持续期为[a,x],v为任意实数(包括分数),表示基于Grümwald-Letnikov定义的分数阶微分算子,Γ为Gamma函数。由分数阶微积分的Grümwald-Letnikov定义式可知,Grümwald-Letnikov定义在欧氏测度下将整数阶微积分的整数步长推广到分数步长,从而将微积分的整数阶推广到分数阶。分数阶微积分的Grümwald-Letnikov定义的计算简便易行,它仅需要与信号s(x)自身相关的的离散采样值,而不需要信号s(x)的导数与积分值。Riemann-Liouville定义信号s(x)的v阶积分(vπ0)为 vπ0。其中,表示基于Riemann-Liouville定义的分数阶微分算子。对于信号s(x)的v阶微分(v≥0),n满足n-1<v≤n。于是由Riemann-Liouville积分定义式,本发明可推导出信号s(x)的v阶微分的Riemann-Liouville定义为 0<vπn。由分数阶微分的Riemann-Liouville定义式,本发明可以推导信号s(x)的Fourier变换为 其中,i是虚数单位,ω是数字频率。当信号s(x)是因果信号时,上式可简化为FT[Dvs(x)]=(iω)vFT[s(x)]。
本发明要构造适于纹理图像去噪的分数阶偏微分方程,就必须构造分数阶Euler-Lagrange方程,而构造分数阶Green公式是构造分数阶Euler-Lagrange方程的必要前提。因此,本发明必须首先将传统的整数阶Green公式推广到分数阶,推导并构造分数阶Green公式。
见图2,本发明令Ω是以分段光滑曲线C为边界的平面单连通区域,可微积函数P(x,y)和Q(x,y)在Ω和C上连续,并存在对x和y的分数阶连续偏导数;令D1表示1阶微分算子,Dv表示v阶分数阶微分算子,I1=D-1表示1阶积分算子,Iv=D-v表示v>0阶分数阶积分算子,表示在平面Ω上的v阶分数阶曲面积分算子(将分数阶积分的Riemann-Liouville定义式从一维推广到二维),表示在曲线C的AC1B段上沿方向的v阶分数阶曲线积分算子,表示在闭曲线C上沿逆时针方向的v阶分数阶闭曲线积分算子;令区域Ω的边界C是由两曲线y=φ1(x),y=φ2(x),a≤x≤b或x=ψ1(y),x=ψ2(y),c≤y≤d所围成。
同理可得 由上述两式,本发明可以推导得分数阶Green公式为 特别地,当和互逆时(该条件要求较高,一般很难满足),即则由上式可得分数阶Green公式在特定条件下的简化表出为 由上式可见,当v1=v2=1时,可推得 传统的整数阶Green公式只是分数阶Creen公式的特例;当v1=v2=v时,可推得这是适用范围广泛的分数阶Green公式。
基于上述推导所得的分数阶Green公式,本发明可以进一步推导针对二维图像处理的分数阶泛Euler-Lagrange方程。
见图3,本发明令二维空间中的可微积数量函数为u(x,y)和可微积矢量函数为v阶分数阶微分算子 Dv是一线性算子(当v=0时,D0表示既不微分也不积分,是一个恒等算子),其中i和j分别表示在x和y方向上的单位矢量。一般而言,二维图像区域Ω是一个长方形的单连通区域,因此Ω的分段光滑边界C是一个闭合的长方形曲线。
由分数阶积分的Riemann-Liouville定义式可得 且 于是由上述分数阶Green公式可推得 恒成立。由于 且 于是本发明可以推导得由分数阶微积分的齐次性质进而可得于是,本发明可以推导得其中,符号.表示内积。与分数阶散度算子的定义类似,本发明令v阶分数阶类微分算子v阶分数阶类散度算子divPv和Pv都是线性算子,于是根据希尔伯特伴随算子理论,本发明可以推导得其中,表示v2阶分数阶内积的积分形式,(Dv)*表示Dv的v阶分数阶希尔伯特伴随算子。由上式本发明可得(Dv)*=-divPv。由上式可见,分数阶希尔伯特伴随算子(Dv)*是一线性算子。当v1=v2=1时,本发明可以推导得其中,<,>1表示一阶内积的积分形式,表示一阶散度算子,(D1)*表示D1的一阶希尔伯特伴随算子。针对数字图像而言,本发明可推导得(D1)*=-div1。于是可知,一阶希尔伯特伴随算子只是分数阶希尔伯特伴随算子的特例。当时,于是可得在上式中,由于x方向和y方向的正交性,对于任意二维数量函数u(测试函数)而言,相应的和具有随机性,由变分基本引理,要使上式成立,只需由于n为1→∞的正整数,欲使仅需使上式的等价式为于是要使上式成立,当且仅当成立。即是与相对应的分数阶Euler-Lagrange方程。
由于分数阶微积分对于所有的v均存在的原因,所以分数阶Euler-Lagrange方程与分数阶面积分的积分阶次v2无关,因此本发明在下面构造适于纹理图像去噪的分数阶偏微分方程模型的能量泛函时均不采用分数阶面积分而只采用一阶面积分的形式。
如下所述,首先,本发明利用图像分数阶微积分的特殊性质,不仅将能量范函中的梯度算子简单地从一阶推广到分数阶,而且采用分数次幂和分数阶极值来构造基于分数阶全变差的能量泛函;然后,对该能量范函求分数阶极值,而非传统的一阶极值;然后,直接运用分数阶微积分的性质来推导相应的分数阶Euler-Lagrange方程,极大地简化了分数阶偏微分方程数值实现的复杂性;最后,本发明构造出一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器。
本发明直接采用分数阶极值来构造能量泛函,令含有小参数β的曲面族s+(β-1)ξ,当β=1时为v3阶极值曲面s。本发明令Ψ1(β)=EFHTV[s+(β-1)ξ],其中Ψ2(β)是噪声[s+(β-1)ξ-s0]和无噪信号s0的交叉能量,是对[s+(β-1)ξ-s0]和s0相似性的度量。求解Ψ2(β)的极小值的过程即是求解噪声和无噪信号的相似性最小的过程。Ψ2(β)在去噪过程中起非线性保真的作用,λ为正则化参数。于是定义在曲面族s+(β-1)ξ上基于分数阶全变差的分数阶能量泛函为
对于Ψ1(β)而言,如果对β的v3阶分数阶导数存在,则当β=1时,Ψ1(β)存在v3阶分数阶极小值(分数阶驻点)。由分数阶微分算子的线性性质可得其中,对于向量而言,符号·表示内积。与基于一阶极值的传统一阶变分法不同,上式是Ψ1(β)关于β的v3阶分数阶极值,其目的在于利用图像分数阶微积分的特殊性质(既能尽量保留图像平滑区域中的低频轮廓特征,同时又能分数阶且多尺度增强图像中灰度值跃变幅度相对较大的高频边缘特征,而且还能分数阶且多尺度增强图像中灰度值跃变幅度和频率变化相对不大的高频纹理细节特征),使纹理图像在去噪的同时尽量分数阶非线性地保留其复杂纹理细节特征。
只要v为分数,可使n>v时组合数 于是由Faàde Bruno公式,本发明可推得复合函数的分数阶微分法则为 其中,a为任意常数。上式中的n=0这一项已从求和符号中分离出来了。由上式可知,复合函数的分数阶导数是一个无穷项求和。其中,Pk满足 上式中的第三个求和符号∑表示对满足 的Pk|m=1→n所有组合所对应的进行求和。于是由Gamma函数的性质可以推导得 不失一般性,为了简化计算,本发明令f(η)=η,故 于是上式可被简化为于是,当m=1时,上式中Pk必须满足Pn=1且P1=P2=Λ=Pn-1=0,于是上式进一步简化为当n取奇数(n=2k+1,k=0,1,2,3Λ)与当n取偶数(n=2k,k=1,2,3Λ)时,具有不同的表达式。于是可以分别推导得和由上述两式,本发明可以推导得 其中,本发明令对于任意二维数量函数ξ(测试函数)而言,相应的具有随机性,于是 亦具有随机性。于是本发明可以推导相应的分数阶Euler-Lagrange方程为 其中, 上式即是相对应的分数阶Euler-Lagrange方程。对于Ψ2(β)而言, 在上式中,由于试验函数ξ具有随机性,所以г(2-v3)(s-s0-ξ)+г(1-v3)ξ亦具有随机性,由变分基本引理,要使上式成立,只需成立。由于上式分别是Ψ1(β)和Ψ2(β)对β求v3阶极小值,本发明取v=v3≠1,2and3,求解分数阶极小值可得 其中, 按照分数阶差分方法实现。本发明需计算λ(t)。若图像噪声n(x,y)为白噪声,于是 当 时,上式收敛于稳定状态。于是在上式两边同时乘以(s-s0)2并在图像区域Ω上积分,上式的左边便消失,可得 式 和 所表示的分数阶偏微分方程去噪模型即为本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的数学模型。
为了使上述分数阶偏微分方程去噪模型能够完全滤除在信号的甚低频和直流部分残留的微弱噪声,本发明在数值迭代实现时,还需同时对信号的甚低频和直流部分进行低通滤波。另外,上述分数阶偏微分方程去噪模型不仅通过连乘函数以及的幂次方v2-2k-2的形式较大增强了分数阶微分阶次v2对去噪的非线性调节作用;而且通过在分母中增加г(-v3)的形式增加了分数阶微分阶次v3对去噪的非线性调节作用。另外,当v3=0时,上述分数阶偏微分方程去噪模型即为传统的位势方程(椭圆型方程);当v3=1时,上述分数阶偏微分方程去噪模型即为传统的热传导方程(抛物型方程);当v3=2时,上述分数阶偏微分方程去噪模型即为传统的波动方程(双曲型方程);当0<v3<1时,上述分数阶偏微分方程去噪模型即为传统的位势方程和传统的热传导方程之间的连续内插;当1<v3<2时,上述分数阶偏微分方程去噪模型即为传统的热传导方程和波动方程之间的连续内插。可见,在数学和物理意义上,上述分数阶偏微分方程去噪模型将传统的基于偏微分方程的图像处理从基于传统的热传导方程的各向异性扩散的普遍基本处理方法推广到了更广阔的领域。
进一步地,本发明需要数值实现上述一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的数学模型。第一,本发明需要数值实现二维数字图像在x轴和y轴方向上的分数阶微分。对于分数阶微积分的Grümwald-Letnikov定义式,当N足够的大时,可以去掉极限符号。为了提高收敛速度和收敛精度,本发明在Grümwald-Letnikov定义式中引入信号s(x)在非节点处的信号值,即 于是当v≠1时,应用拉格朗日三点插值公式对信号s(x)进行分数插值,可分别构造出数字图像在x轴和y轴的方向上的分数阶微分算子。对数字灰度图像的而言,分数阶微分算子的数值运算规则采用算子卷积的空域滤波方案。本发明选取在x轴和y轴方向上的模值最大的分数阶偏微分值作为该像素点的分数阶微分值。第二,本发明需要数值实现二维数字图像在x轴和y轴方向上的1阶微分。为了保持数值计算的稳定性,本发明采用 和 来近似一阶微分。第三,本发明需要数值实现二维数字图像对于时间t的分数阶微分。若时间等分间隔为Δt,即单位迭代时间间隔,n时刻为tn=nΔt,n=0,1,Λ(t0=0表示初始时刻)。本发明取单位迭代时间间隔Δt在(0,0.1]内取任意较小的正实数。于是n时刻的数字图像为待去噪的原始图像为s0为理想的无噪声图像,它是一个恒定值,故s0(x,y,t0)=s0(x,y,tn)。于是,本发明可将二维数字图像对于时间t的分数阶微分近似为 其中,是最佳去噪图像。另外,由于理想的无噪声图像s0(x,y,t0)事先不知道,但是每次数值迭代的去噪中间结果都是对理想的无噪声图像s0(x,y,t0)一次逼近,即故为了在数值迭代时尽量逼近s-s0,本发明令由于最佳去噪图像事先不知道,但是每次数值迭代的去噪中间结果都是对的一次逼近,即故为了在数值迭代时尽量逼近本发明令 同时为了简化计算,本发明取μ=0.005和η=1,只取k=0,1进行近似计算。于是,可以推导得式 和 的数值实现方程分别为 和 其中, 其中,在数值迭代计算时,一方面,本发明不需要预先获知或估计噪声的方差,而只需要令第一次数值迭代时的为一个较小的正数。本发明取 将带入 以启动数值迭代计算的过程,于是每一次迭代所得的都不一样,但每一都是对噪声真正方差的一次逼近;另一方面,在数值迭代计算的过程中,可能出现的情况,为了使有意义,当时,本发明取 为了使有意义,当 时,本发明取
基于上述对本发明所提出的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的数学公式推导和数值运算规则的简要说明,下面具体说明该滤波器的电路构成:
见图1,本发明所提出的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器是基于一种特殊分数阶偏微分方程去噪算法来实现对纹理图像的分数阶、非线性、多尺度、高精度去噪。本发明涉及的分数阶微积分的阶次v1、v2和v3不是传统的整数阶,而是非整数阶,工程应用中一般取分数或有理小数,且v3≠1,2,3。见图1,该滤波器是采用缓存器2、差值平方器一3、-v3次方幂方器4、乘法器一5、λn发生器11、发生器10、乘法器二8、乘法器三9和加法器一6以级联方式构成的。该滤波器具有在去除纹理图像噪声的同时,既能尽量保留平滑区域中的低频轮廓,同时又能非线性保留灰度值跃变幅度相对较大的高频边缘,而且还能非线性保留灰度值跃变幅度变化相对不大的纹理细节的显著优点。该滤波器特别适用于对富含复杂纹理细节特征的图像进行去噪的应用场合。
见图1,1是本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的输入点,即第n次迭代的数字图像的输入点。缓存器2所起的作用是将本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的输出图像进行两个单位迭代时间的缓存,从而得到第n-1次迭代的数字图像差值平方器一3的E输入点是E点是权值的输入点。差值平方器一3完成的计算是-v3次方幂方器4完成的计算是乘法器一5的F输入点是权值的输入点。乘法器一5完成的计算是 发生器10完成的计算是 λn发生器11完成的计算是 乘法器三9的H输入点是权值的输入点。乘法器三9完成的计算是乘法器二8的G输入点是权值的输入点。乘法器二8完成的计算是加法器一6完成的计算是 其中,本发明要求阶次v3满足v3≠1,2,3。本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的输出点7完成的功能是输出第n+1次迭代的数字图像
见图4,发生器10完成的计算是 差值平方器二12的I输入点是差值平方器二的权值的输入点,即第0次迭代的原始数字图像的输入点。差值平方器二12完成的计算是 发生器13完成的计算是乘法器四14完成的计算是乘法器五15完成的计算是 加法器一16完成的计算是除法器17的J输入点是权值г(1-v3)г(2-v3)的输入点。除法器17完成的计算是乘法器六18完成的计算是 λn发生器的输出点19完成的功能是输出λn值。
见图5,微分器20完成的计算是在x轴方向上的v1阶分数阶微分。微分器21完成的计算是在x轴方向上和y轴方向上同时进行v1阶分数阶微分。微分器22完成的计算是在y轴方向上的v1阶分数阶微分。v2-2k-2幂方器23完成的计算是乘法器七24完成的计算是乘法器八25完成的计算是 微分器26完成的计算是在x轴方向上的1阶微分。微分器27完成的计算是在y轴方向上的1阶微分。加法器二28完成的计算是 乘法器九32的K输入点是权值的输入点。乘法器九32完成的计算是 加法器31完成的计算是 乘法器十30的L输入点是权值的输入点。乘法器十30完成的计算是 发生器的输出点29完成的功能是输出值。
下面结合附图和实例详细说明本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的新方案:
附图说明
图1是本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的示意图。
图2是单连通区域Ω及其光滑边界曲线C的示意图。
图3是二维单连通图像区域Ω及其分段光滑边界曲线C的示意图。
图4是λn发生器的示意图。
其中,1是本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的输入点,即第n次迭代的数字图像的输入点;2是缓存器;3是差值平方器一;4是-v3次方幂方器;5是乘法器一;6是加法器一;7是本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的输出点,即第n+1次迭代的数字图像的输出点;8是乘法器二;9是乘法器三;10是发生器;11是λn发生器;12是差值平方器二:13是发生器;14是乘法器四;15是乘法器五;16是加法器一;17是除法器;18是乘法器六;19是λn发生器的输出点,即λn的输出点;20是微分器;21是微分器;22是微分器;23是v2-2k-2幂方器;24是乘法器七;25是乘法器八;26是微分器;27是微分器;28是加法器二;29是发生器的输出点,即的输出点;30是乘法器十;31是加法器;32是乘法器九;33是差值平方器三;34是加法器二;35是发生器的输出点,即的输出点。另外,3、12和33是功能和参数相同的差值平方器;5、8、9、14、15、18、24、25、30和32是功能和参数相同的乘法器;6和28是功能和参数相同的加法器;16和34是功能和参数相同的加法器。
具体实施方式
现举例介绍如下:
见图1、图4、图5和图6,在工程实际应用中,本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器涉及的分数阶微积分的阶次v1、v2和v3不是传统的整数阶,而是非整数阶,一般取分数或有理小数。另外,由上述说明可知如果本发明取阶次v1=2.25,v2=2.5,v3=0.25,取单位迭代时间间隔Δt=0..01,方差初始值取参数μ=0.005,于是可以推导得本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的具体电路参数。见图1,乘法器一5的F输入点的输入权值为乘法器二8的G输入点是输入权值为 乘法器三9的H输入点的输入权值为。见图4,除法器17的J输入点的输入权值为Γ(1-v3)Γ(2-v3)=Γ(0.75)Γ(1.75)=1.1262。见图5,乘法器九32的K输入点的输入权值为乘法器十30的L输入点的输入权值为 于是,如图1、图4、图5和图6所示,按照本说明书的发明内容中所详细说明的本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的级联电路结构及其具体电路参数,就可以方便地构造出该基于分数阶偏微分方程的适于纹理图像去噪的高精度滤波器的具体电路。在不影响准确表述的前提下,为了更加清晰明了地描述本发明的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器的具体电路,图1、图4、图5和图6未画出其中的时序控制电路及其被触发产生的时序控制信号。
Claims (6)
2.根据权利要求1所述的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器,其特征在于:该滤波器的输入点(1)输入第n次迭代的数字图像缓存器(2)将本该基于分数阶偏微分方程的适于纹理图像去噪的高精度滤波器的输出图像进行两个单位迭代时间的缓存,从而得到第n-1次迭代的数字图像差值平方器一(3)的输入权值为差值平方器一(3)完成的计算是-v3次方幂方器(4)完成的计算是乘法器一(5)的输入权值为乘法器一(5)完成的计算是 发生器(10)完成的计算是 λn发生器(11)完成的计算是 乘法器三(9)的输入权值为乘法器三(9)完成的计算是乘法器二(8)的输入权值为乘法器二(8)完成的计算是加法器一(6)完成的计算是 该基于分数阶偏微分方程的适于纹理图像去噪的高精度滤波器的输出点(7)输出第n+1次迭代的数字图像
4.根据权利要求1所述的一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器,其特征在于:其中微分器(20)完成的计算是在x轴方向上的v1阶分数阶微分。微分器(21)完成的计算是在x轴方向上和y轴方向上同时进行v1阶分数阶微分。微分器(22)完成的计算是在y轴方向上的v1阶分数阶微分。v2-2k-2幂方器(23)完成的计算是乘法器七(24)完成的计算是乘法器八(25)完成的计算是 微分器(26)完成的计算是在x轴方向上的1阶微分。微分器(27)完成的计算是在y轴方向上的1阶微分。加法器二(28)完成的计算是 乘法器九(32)的输入权值为乘法器九(32)完成的计算是 加法器(31)完成的计算是 乘法器十(30)的输入权值为乘法器十(30)完成的计算是 发生器的输出点(29)输出值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310298382.6A CN103337060B (zh) | 2013-07-17 | 2013-07-17 | 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310298382.6A CN103337060B (zh) | 2013-07-17 | 2013-07-17 | 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103337060A true CN103337060A (zh) | 2013-10-02 |
CN103337060B CN103337060B (zh) | 2015-09-23 |
Family
ID=49245210
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310298382.6A Expired - Fee Related CN103337060B (zh) | 2013-07-17 | 2013-07-17 | 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103337060B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761715A (zh) * | 2014-01-23 | 2014-04-30 | 沈阳大学 | 一种用于图像去噪的分数阶原始对偶方法 |
CN104517266A (zh) * | 2014-12-22 | 2015-04-15 | 南京信息工程大学 | 基于边缘检测算子的混合自适应图像去噪方法 |
CN104616249A (zh) * | 2014-11-24 | 2015-05-13 | 南京信息工程大学 | 一种基于曲率变分的小波变换图像去噪算法 |
CN106910215A (zh) * | 2017-03-15 | 2017-06-30 | 沈阳理工大学 | 一种基于分数阶梯度插值的超分辨率方法 |
CN106981056A (zh) * | 2017-03-31 | 2017-07-25 | 蒲亦非 | 一种基于分数阶偏微分方程的图像对比度增强滤波器 |
CN107085840A (zh) * | 2017-06-16 | 2017-08-22 | 南京信息工程大学 | 基于分数阶偏微分方程的图像去噪方法 |
CN108198139A (zh) * | 2017-12-01 | 2018-06-22 | 西安电子科技大学 | 基于传导全变分正则化的图像去噪方法 |
CN108665431A (zh) * | 2018-05-16 | 2018-10-16 | 南京信息工程大学 | 基于k-均值聚类的分数阶图像纹理增强方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917576A (zh) * | 2006-08-30 | 2007-02-21 | 蒲亦非 | 数字图像的分数阶微分滤波器 |
CN101800847A (zh) * | 2010-04-14 | 2010-08-11 | 蒲亦非 | 基于Riemann-Liouville定义的数字图像的1~2阶分数阶微分滤波器 |
CN101848319A (zh) * | 2010-04-02 | 2010-09-29 | 蒲亦非 | 高精度计算的数字图像的分数阶微积分滤波器 |
CN102262778A (zh) * | 2011-08-24 | 2011-11-30 | 重庆大学 | 基于改进的分数阶微分掩模的图像增强方法 |
CN102800051A (zh) * | 2012-03-28 | 2012-11-28 | 中南大学 | 基于分数阶微分的自适应浮选泡沫图像增强方法 |
-
2013
- 2013-07-17 CN CN201310298382.6A patent/CN103337060B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917576A (zh) * | 2006-08-30 | 2007-02-21 | 蒲亦非 | 数字图像的分数阶微分滤波器 |
CN101848319A (zh) * | 2010-04-02 | 2010-09-29 | 蒲亦非 | 高精度计算的数字图像的分数阶微积分滤波器 |
CN101800847A (zh) * | 2010-04-14 | 2010-08-11 | 蒲亦非 | 基于Riemann-Liouville定义的数字图像的1~2阶分数阶微分滤波器 |
CN102262778A (zh) * | 2011-08-24 | 2011-11-30 | 重庆大学 | 基于改进的分数阶微分掩模的图像增强方法 |
CN102800051A (zh) * | 2012-03-28 | 2012-11-28 | 中南大学 | 基于分数阶微分的自适应浮选泡沫图像增强方法 |
Non-Patent Citations (4)
Title |
---|
蒲亦非 等: "数字图像纹理细节的分数阶微分检测及其分数阶微分滤波器实现", 《中国科学E辑:信息科学》 * |
黄果 等: "分数阶微积分在图像处理中的研究综述", 《计算机应用研究》 * |
黄果 等: "基于时间-空间分数阶偏微分方程的图像去噪模型", 《系统工程与电子技术》 * |
黄果 等: "基于空间分数阶偏微分方程的图像去噪模型研究", 《四川大学学报(工程科学版)》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103761715A (zh) * | 2014-01-23 | 2014-04-30 | 沈阳大学 | 一种用于图像去噪的分数阶原始对偶方法 |
CN104616249A (zh) * | 2014-11-24 | 2015-05-13 | 南京信息工程大学 | 一种基于曲率变分的小波变换图像去噪算法 |
CN104616249B (zh) * | 2014-11-24 | 2018-01-30 | 南京信息工程大学 | 一种基于曲率变分的小波变换图像去噪方法 |
CN104517266A (zh) * | 2014-12-22 | 2015-04-15 | 南京信息工程大学 | 基于边缘检测算子的混合自适应图像去噪方法 |
CN104517266B (zh) * | 2014-12-22 | 2017-06-06 | 南京信息工程大学 | 基于边缘检测算子的混合自适应图像去噪方法 |
CN106910215A (zh) * | 2017-03-15 | 2017-06-30 | 沈阳理工大学 | 一种基于分数阶梯度插值的超分辨率方法 |
CN106910215B (zh) * | 2017-03-15 | 2022-01-07 | 沈阳理工大学 | 一种基于分数阶梯度插值的超分辨率方法 |
CN106981056B (zh) * | 2017-03-31 | 2020-01-24 | 蒲亦非 | 一种基于分数阶偏微分方程的图像对比度增强滤波器 |
CN106981056A (zh) * | 2017-03-31 | 2017-07-25 | 蒲亦非 | 一种基于分数阶偏微分方程的图像对比度增强滤波器 |
CN107085840A (zh) * | 2017-06-16 | 2017-08-22 | 南京信息工程大学 | 基于分数阶偏微分方程的图像去噪方法 |
CN108198139B (zh) * | 2017-12-01 | 2021-09-10 | 西安电子科技大学 | 基于传导全变分正则化的图像去噪方法 |
CN108198139A (zh) * | 2017-12-01 | 2018-06-22 | 西安电子科技大学 | 基于传导全变分正则化的图像去噪方法 |
CN108665431A (zh) * | 2018-05-16 | 2018-10-16 | 南京信息工程大学 | 基于k-均值聚类的分数阶图像纹理增强方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103337060B (zh) | 2015-09-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103337060B (zh) | 一种基于分数阶偏微分方程的纹理图像高精度去噪滤波器 | |
CN103427789B (zh) | 一种基于分数阶计算方程的图书馆图文信息去噪滤波器 | |
Wang et al. | A review of flux reconstruction or correction procedure via reconstruction method for the Navier-Stokes equations | |
Ryland et al. | On multivariate Chebyshev polynomials and spectral approximations on triangles | |
CN103390265B (zh) | 一种基于分数阶发展方程的纹理图像去噪滤波器 | |
CN105389839A (zh) | 基于流体分析的流体参数估计方法 | |
Sussman | Rademacher Series for $\eta $-Quotients | |
Wang | An efficient iterative method for reconstructing surface from point clouds | |
CN103325100B (zh) | 一种基于分数阶热传导方程的图像去噪滤波器 | |
Zelik et al. | Green's function asymptotics and sharp interpolation inequalities | |
Fayolle et al. | p-Laplace diffusion for distance function estimation, optimal transport approximation, and image enhancement | |
Badriev et al. | A Study of Variable Step Iterative Methods for Variational Inequalities of the Second Kind. | |
Torsello et al. | Correcting curvature-density effects in the Hamilton–Jacobi skeleton | |
Cotter et al. | Bayesian data assimilation in shape registration | |
Branding | Magnetic Dirac-harmonic maps | |
Wang et al. | Fractional variational integrators for fractional Euler–Lagrange equations with holonomic constraints | |
Shindin et al. | Analysis of a Chebyshev-type pseudo-spectral scheme for the nonlinear Schrödinger equation | |
Carlini et al. | Numerical techniques for level set models: an image segmentation perspective | |
Fayolle et al. | Rounding, filleting and smoothing of implicit surfaces | |
Bala et al. | Surface reconstruction by evaluating height from gradient values | |
Fu et al. | A wavelet multiscale iterative regularization method for the parameter estimation problems of partial differential equations | |
Ren et al. | A Walsh-function-based Finite Volume Method to Capture Discontinuity Inside Grid Cell | |
Roberts et al. | Discrete thin plate spline smoothing in 3D | |
Zhang et al. | An Improved Medical Image Denoising Algorithm Based on One-Dimensional Heat Transfer Equation | |
Pisani et al. | Lee-yang zeros and stokes phenomenon in a model with a wetting transition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150923 Termination date: 20190717 |