蒲亦非-定义-答案-分数阶微积分在图像处理中的研究综述

第29卷第2期计算机应用研究

V01.29No.2

2012年2月

ApplicationResearchofComputers

Feb.2012

分数阶微积分在图像处理中的研究综述

黄果h,许黎岫,蒲亦非2

(1.乐山师范学院a.智能信息处理及应用实验室;b.物电学院,四川乐山614000;2.四川大学计算机学院,成

都610064)

摘要:综述了关于分数阶微积分理论在数字图像底层处理中的应用研究,具体包括:分数阶微积分和分数阶偏微分方程的基本理论及分数阶傅里叶变换的基本性质;分数阶微分滤波器的构造及在图像增强中的应用研究;分数阶积分滤波器的构造及在图像去噪中的应用研究;分数阶偏微分方程在图像去噪中的应用研究。最后,总结了分数阶微积分理论在图像处理中已取得的研究成果,并结合已有的基于分数阶微积分理论的图像底层处理模型,展望了分数阶微积分理论在图像处理中的应用前景。

关键词:分数阶微积分;分数阶偏微分方程;变分法;泛函极值;图像增强;图像去噪中图分类号:TN391

文献标志码:A

文章编号:1001.3695(2叭2)02—0414—07

doi:10.3969/j.issn.100l一3695.2012.02.003

Summaryofresearch

on

imageprocessingusingfractionalcalculus

HUANGGu01a,XULilb,PUYi.fei2

(1.Ⅱk60m£o,y矿胁扰咖眦聊备玎№砌nProc哪堍&A一池t如n,6.s如00f0,P蛳如&觑∞打口n如,肠^口nⅣ0r,mz‰妇m毋,k^口n_si-如Mn614000,醌iM;2.co朋e酽o,compM把r&如袱,.s如^∽n‰眈巧蚵,蕊e,则M610064,醌i胍)

Abstmct:The

maincontents

ofthispaperwereaboutfractionalcalculusapplicationindigitalimagelow—levelprocessing.

SpecificaⅡyincludedthefollowing:fractionalcalculusaJldfmctional—orderpartialdifEbrentialequationsandthebasicprope卜tiesoffractionalFouriertransfo册:thestmctureofthefractionaldifferentialfilteranditsapplicationtoimageenhancement;

thestmctureofthefractionalintegmtionfilteranditsapplication

to

image

denoising;theimagedenoisingmodelbased

on

fhc—

tionalpartialdi如rentialequations.FinaJly,summarizedfractionalcalculushadbeenmaderesea托hresultsandrelatedpmb—lemsinimagepmcessing.Associatedwiththeimageprocessingmodels,forecastedtheapplicationsprospectaboutfhctionalcalculusappliedtoimageprocessing.

Keywords:fhctionalcalculus;fractional-orderpartialdi矗Ierentialequations;calculus

ofvariations;functionalextreme;

image

enhancement;imagedenoising

生物、金融、流体等领域应用越来越广泛,实用价值也得以充

引言

分的体现。近几年以来,国内外许多学者研究发现,分数阶分数阶微积分理论与整数阶微积分理论几乎同时诞生,

微积分在信号分析与处理领域具有广泛的应用前景,比如在该理论最早出现在L’H6pital写给Leibniz的信件中…。近三信号的奇异性检测和提取方面具有特殊的作用H“…。分数阶微积分理论在图像底层处理中的应用是最近几年才引起百年来,分数阶微积分仅仅是数学家们的纯数学理论分析和学者关注的,到目前为止,已经取得了初步研究成果,主要包推导。Liouville等著名数学家一直致力于完善和发展分数阶括图像增强、图像去噪、图像边缘提取、图像分割、图像奇异微积分的理论研究,从而构建了分数阶微积分理论的基本系性检测等。图像处理中的整数阶微分算子包括梯度和La—统架构。直到19世纪后期,分数阶微积分理论才开始逐渐在place算子,以及墨西哥帽算子等很多经典的图像处理算子都实际工程领域中得到初步应用。Mandelbmt旧’31发现了分形可以近似看做是Rodieck…1提出的DOG感受野模型在理论学,他指出自然界和很多技术科学中存在大量分数维事实,上的仿真实现。分数阶微分算子则可以认为是DOG感受野发现在整体与部分之间存在自相似现象,并尝试将Riemann-模型的推广和延伸。基于分数阶微分理论的感受野模型较Liouville分数阶微积分理论应用于描述和分析分形媒介中的整数阶微积分理论更符合人类视觉特性。

布朗运动,从此,工程技术人员开始关注分数阶微积分。近分数阶微积分理论在图像处理领域的应用研究主要包括几十年以来,国内外许多学者研究发现分数阶微积分算子具以下方面:蒲亦非等人¨2。2o认为,选择适当阶次的分数阶微分有记忆性和非局部性,非常适合用来描述现实世界中具有记算子在增强图像过程中可以大幅提升边缘和纹理细节的同时,忆和遗传性质的材料,因此分数阶微积分理论在物理、化学、

非线性地保留了图像平滑区域的纹理信息;此外,杨柱中等

收稿日期:201l-09・16;修回日期:2011・lO一30

基金项目:国家自然科学基金资助项目(60972131);四川省教育厅科研项目(10zC019/

11zBl32);乐山市科技局重点研究计划资助项目(2011GzD046);乐山师范学院资助项目(z1162)

作者简介:黄果(1980-),男,重庆人,讲师,博士,主要研究方向为分数阶微积分及分数阶偏微分方程在图像处理中的应用研究(huangguoxuli@163.com);许黎(1982一),女,云南人,讲师,硕士,主要研究方向为分数阶微积分及混沌理论在图像处理中的应用研究;蒲亦非(1976一),男,四川人,副教授,博士(后),主要研究方向为分数阶微积分及分数阶偏微分方程在图像处理中的应用研究.

万方数据

第2期黄果,等:分数阶微积分在图像处理中的研究综述

・415・

人。2”驯利用分数阶微分算子具有弱导数性质的特点,利用分数阶梯度算子对含弱噪声的图像进行边缘检测,该方法能够有效地解决整数阶梯度算子对噪声敏感的问题,因此可以避免噪声的影响而准确地定位噪声图像的边缘;MatIlieu等人心叫提出了分数阶微分的边缘检测算子,即分数阶鲁棒轮廓(contour

r0.

busted’ordre

non

entier,cRONE)边缘检测器,说明了该检测

器的分数阶微分阶次在1<”<2时,能够有选择地检测出边缘,而在分数阶微分阶次在一1<”<1时,该检测算子能够在边缘提取的过程中克服噪声的影响;在此基础上,李远禄等人旧川提出了基于分数阶差分的滤波器并将其应用到边缘检测,该模型将数据的平滑和差分有机结合,利用分数阶微积分基本理论,设计了可以通过调节微分阶次的差分滤波器,该滤波器可以解决传统算子边缘检测出现的边缘漂移问题,并且可以抑制部分噪声;汪凯宇和刘红毅等人。28’2引利用分数阶微积分具有的记忆特性,分别将分数阶样条小波应用到图像纹理的奇异性检查和图像融合中,较整数阶微积分取得了更好的仿真效果;Liu等人∞叫提出了分数阶奇异值分解的人脸识别方法,该方法可以有效缓解面部的变化,并且在脸部出现剧烈变化时,比传统的分类方法提高其分类的性能,为高层图像处理打下坚实基础;左凯等人口1。将卡尔曼滤波和分数阶微积分理论相结合,提出了二维分数阶卡尔曼滤波器并成功应用到图像处理;汪成亮等人¨2。通过研究分数阶微积分的基本定义和相应的分数阶微分算子的实现方法,提出了将图像的梯度特征和人类视觉特征等理论引入到已有的分数阶微分算子中,由此构建了基于分数阶微分阶次自适应变化的图像增强模型;高朝邦等人”引将四元数理论和分数阶微积分理论有机结合,提出了四元数分数阶方向微分的概念,解释其物理意义和几何意义,并根据分数阶微分的特殊性质将其应用到图像增强当中,可以得到很好的效果;Bai等人旧4o以ROF去噪模型为基础,将分数阶微积分理论和偏微分方程相结合,提出了基于分数阶偏微分方程的图像去噪模型,该方法可以解决传统低阶次整数阶偏微分方程去噪模型容易产生阶梯效应,以及高阶次整数阶偏微分方程去噪模型去噪效果不佳的缺点。此后,张军等人"”驯将负指数sobolev空间的多尺度图像建模与基于分数阶微积分的图像建模有效结合,提出了统一的基于分数阶多尺度变分图像去噪模型,并初步设计了该模型涉及参数的自适应选择方法。

相关理论及分析

1.1

分数阶微积分的基本定义

从不同的应用角度去分析问题可以得到不同的分数阶微

积分定义,所以迄今为止,分数阶微积分仍然没有统一的时域定义表达式。在众多的定义表达中,主要有三种经典的分数阶微积分定义,包括Gmmwald—Letnikov、Riemann.IJiouville和ca.potu定义。

1)分数阶微积分的Grnmwald—Letnikov定义

:c=2骂∥j叁黼t一曲)weR一

(1)

2)分数阶微积分的Riemann—Liouville定义

『警

【专赢可,:茄ay。一㈨如

…eⅣ

万方数据

3)分数阶微积分的C印uto定义

细撇卜k志乃筘…川…。@’。

【鲁币%J:带如嚷”k…

f掣

…eⅣ

针对三种经典分数阶微积分定义,国内外众多学者从不同应用角度研究分析各种各样的实际问题,得到了不同形式的分数阶微积分定义表达式。分数阶微积分Riemann一“ouville定义和分数阶微积分c印uto定义都是对分数阶微积分Gnimwald.Letnikov定义的改进。分数阶微积分的Grumwald.Letnikov定义在数值实现时可以转换为卷积运算形式,因此非常适合在信号处理中的应用;分数阶微积分的Riemann—Li.ouviue定义主要应用于计算一些较为简单函数的解析解;分数阶微积分的c印uto定义适用于分数阶微分方程初边值问题的分析,因此非常适合在工程领域中应用。三种关于分数阶微积分定义之间存在着紧密的联系,在一定条件下可以相互转换。例如,如果分数阶微积分阶次w满足n一,1<”<凡,当函数,(x)的m+1阶导数连续并且满足m=L

n一1

j,分数阶微积分

Gmmwald—Letnikov的定义与分数阶微积分Riemann—Liouville的定义是完全等价的。在满足相同条件下,分数阶微积分的c叩uto定义和分数阶微积分的Gmmwald—Letnikov定义等价。当分数阶微积分阶次”为正整数和负实数时,分数阶微积分的c印uto定义和分数阶微积分的Riemann—Liouville定义满足以下关系:

:D㈧=:乒糍孵+:D荆(4)

1.2分数阶偏微分方程

分数阶偏微分方程是由传统整数阶偏微分方程拓展演化而得到的,是将传统整数阶偏微分方程中对时间或者空间的偏导数项用分数阶偏导数来替代而得到的㈤J。分数阶偏微分方程的基本形式为

!丝!!型2:拯!:吐+塑!!也

f一

5、

矿t妒x

护y

’根据偏微分阶次d和卢的不同将分数阶偏微分方程分为三种基本形式㈤““:

a)时间分数阶偏微分方程(时间偏导数阶次d为非整数,空间偏导数阶次卢为整数)。特别地,当0<d<1,口=2时,将

热传导方程(扩散方程)推广为时间分数阶热传导方程(时间分数阶扩散方程);当1<“<2,口=2时,该方程可以认为是热传导方程和波动方程的非线性插值。

b)空间分数阶偏微分方程(空间偏导数阶次口为非整数,时间偏导数阶次仅为整数)。特别地,当d=1,l<届<2时,将热传导方程(扩散方程)推广为空间分数阶热传导方程(空间分数阶扩散方程)。

c)空间一时间分数阶偏微分方程(时间偏导数阶次d为

非整数,同时空间偏导数阶次卢为非整数)。特别地,当0<

a<1,l<卢<2时,将热传导方程(扩散方程)推广为时间一空间分数阶热传导方程(时间一空间分数阶扩散方程)。

图像增强是数字图像处理的基本内容之一,以梯度算子和L印lacian算子等为代表的整数阶微分方法在图像增强中的应

2分数阶微分与图像增强

・416・

计算机应用研究第29卷

用已经相当广泛。近年来,很多学者将分数阶微分理论应用到图像增强中并且取得了较好的效果。分数阶微分可以在提高信号高频成分的同时非线性地保留信号中的甚低频,因此将该理论应用到图像增强中,可以使图像边缘更加突出,同时保留了图像平滑区域的纹理信息。分数阶微分理论应用到数字图像中的关键是增加一个自由度,即阶次。(0<”<1),通过适当调节”的大小可以构造相应的掩模算子,从而得到较好的图像增强效果。

2.1

分数阶微分对信号的作用

已知任意平方可积的能量信号以菇)Er(R),由信号处理

的基本理论可知其Fourier变换为八甜)=f。八戈)e一”出。假设信号八搿)的乃阶导数为尸(石)(n∈z+),根据Fourier变换的性质可以得到

盯^

D’弧x)铮(口D“(∞)=(妇)4以m)=d“((cJ),(∞)(6)

假设信号八戈)的分数阶”阶导数为厂(戈)(”∈月+),同样根据分数阶Fourier变换的性质可以得到

n^

D杈x)甘(功9。(∞)=(奶)”,(∞)=d”(∞班∞)

(7)

其帆

。..

b&)_l。附(乞):罕。gn(。)

一m)=(洳)”=“”(o)e‘∥‘“’

(8)

图1为分数阶微分算子在阶次不同时的幅频特性。从图1可以清楚地得出:a)分数阶微分算子对信号都有加强的作用,并且加强的幅度随频率和微分阶次的增加而非线性地急剧增强;b)2阶次微分相比1阶次微分对信号低频成分的提升较弱,但是对信号高频成分的提升却更明显;c)特别地,当微分阶次”=0时,表示分数阶微分算子不对信号进行任何改变;d)对于分数阶

微分阶次0<”<1时,在信号的低频部分(∞<1),分数阶微分算

子对信号的幅值进行了一定的提升,并且提升的幅度稍大于1阶次和2阶次微分;在信号的高频部分(∞>1),分数阶微分算子对信号的幅值进行一定的提升,但是提升的幅度也明显小于1阶次和2阶次微分。上述性质表明,分数阶微分算子具有弱导数性质,即分数阶微分算子在加强信号中高频成分的同时.对

信号的甚低频分量进行了非线性保留。此外,从信号处理的角

度分析分数阶微积分算子的物理意义可知,分数阶微分算子可以理解为广义的调幅调相,其振幅随频率呈分数阶幂指数变化,相位是频率的广义Hilbert变化。

一I三潍;

S2

{==茹5a

.∥/

//

墨・

一匕≤芝≥0

,一一一j曩互国龟罗舒;t帮j=

.盖矗笙孝三哆r

.q

频率w/Hz

图1分数阶微分算子的幅频特性曲线

2.2分数阶微分滤波器的构造

假设二维图像信号,(省,y)对于x和y轴方向的分数阶微分在一定的条件下是可分离的,这样可以利用分数阶微分算子分别对数字图像的八个方向进行滤波处理,然后将得到的中间结果通过相关规则结合。将图像信号,(立,y>)的持续期[o,£]

按单位间隔^=1等分,即n=I半l=[卜口】。这样可以得

r+

。1^=l

,‘

到{G—L定义下的分数阶微分算子沿x和y轴方向的数值计算

万方数据

表达式,如式(9)和(1U)所不,兵甲代表不分颧彤r微分效值计

:彤,(。,y),垒,(*,y)+(一口),(。一1,y)+L立毛型,(*一2,y)¨

算的余项,如式(11)和(12)所示。

i—二!尘jj二!二;;土j趔,(x一3,y)+R:(。,y)(9)

…。~7一"…

、‘7

:彤,(。,y),垒,(x,y)+(一口),(。,y一1)+L型妥竺业,(。,),一2)+.

R:(z)=ii—二!堕j。二!二;;_2j趟,(x,y一3)+R:(。,y)(10)

—二兰立j—=6

兰2生三土jji旦‘、’。¨7‘‘*、。’wjI兰幽,(x一4,,,)+、一。7

月;(z)=i二二旦2』乏‰m吨y)...+

,(m+1),(一口)…’…”7

—二!_=lL上{i二竺二二兰羔幽(z,,,..4)+…+

(11)、…

蒯‰职一-m)

,(m+1)厂(一")…。^~

(12)、…

图2表示分数阶微分图像增强掩模算子形,该算子将式(9)和(10)推广到图像的其余六个方向,从而得到一个基于八个方向的图像增强掩模算子(图中表格空白部分用0填充)。

该增强算子具有各向旋转不变性,增强算子形中的滤波系数

如式(13)所示。

%o=1%1=一口

%:巫掣

%:=型子坐

‘13’

%。=巫型必乏芋上业业

00%。O0

%。

00

00O阡0,0%,

O%,

000

阡名.

阡名.

阡名

阡名:

阡名.降名。阡0.%

%。

00

阡名.%.%,

阡名:

00

0OOO

%。

%。

图2分数阶微分增强掩模算子

3分数阶积分与图像去噪

众所周知,分数阶微分和分数阶积分是互为逆运算,同样图像增强和图像去噪也是相反的过程。既然分数阶微分

理论在图像增强中取得了很好的效果,因此,本章尝试将分

数阶积分理论用到图像去噪,提出了基于分数阶积分的图像去噪算法。该方法打破了传统图像去噪算法基于整数积分阶次的思想,通过设定一个较小的分数积分阶次来构建相应的图像去噪掩模,并利用迭代的方法来控制图像去噪的效果,从而在图像去噪的同时,尽量保留了图像的边缘和纹理等细节信息。

3.1

分数阶积分对信号的作用

对已知任意平方可积的能量信号以£)∈r(R),由信号处

理的基本理论可知其Fourier变换蒯∞)=k,(£)e一…出。假

第2期

黄果,等:分数阶微积分在图像处理中的研究综述

・417・

设信号八f)的n阶导数为厂(t)(nEz),根据Fourier变换的性质可以得到

Fr^

D杈f)铮(口D“(∞)=(池)”,(ccJ)=d“(∞)八∞)

(14)

这样将信号八£)的分数阶”导数为,(f)(”∈R+),同样根

据Fourier变换的性质可以得到

盯^

D弧£)甘(口D”(m)=(妇)。八∞)=d”(∞小∞)

(15)

rd”(m)=(汹)”=n”(∞)e妒(“)

其中:

‰o川计删o):孚。gIl(。)

(16)

于是,当”∈R一时,根据分数阶算子理论,即D~=,,并且”’=一”,可以得到分数阶积分算子的Fourier变换形式为

盯^

,。八£)甘(矿)”‘(∞)=(洒)”“c【,)=∥(∞)・,(∞)

(17)

r,’(甜)=d”’(埘)e‘伊I“)

其吼

k(:川计。∥(2):半。gn(。)

(18)

根据式(18)的描述,可以得出分数阶积分算子在阶次。’0,0.2,0.5,0.8,1.0,1.5,2.0时的幅频特性,如图3所示。

j3.

l\

I—v=Q.Q

;l{}S2.

¨’;

+v=1.Ov=0.5

—*一v=0.8~

胃趔\N酞

+v=1.5一罂1.

\C再N

一v=2.0

\ ̄ ̄l ̄

~~一鼍杂心

O.

1嫌遵基釜斟≮娶堇j

——寻==j!

频率w/Hz

图3分数阶积分的幅频特性曲线

从图3可以清楚地得出:a)分数阶积分对信号都有衰减作用,并且随着频率和积分次数的增加而非线性地急速衰减;b)2阶次积分相比1阶次积分对信号低频成分的提升较大,同时对信号的高频成分的削弱也明显;c)特别地,当”=0时,表示分数阶积分算子不对信号进行任何改变;d)对于积分次数0<w<1时,在信号的频率∞<l的区域,分数阶积分算子对信号的幅值进行了一定的提升,但是提升的幅度远小于1次和2次积分;在信号的频率∞>1的区域,分数阶积分算子对信号的幅值进行了一定的削弱,但是削弱的幅度也明显小于1次和2次积分。上述性质表明,分数阶积分算子在削弱信号高频部分的同时,对其中的最高频部分进行了非线性的保留;分数阶积分算子在加强信号低频部分的同时,对其中的最低频部分也进行了一定的保留。这样,分数阶积分算子在去除噪声的同时,一定程度上保留了图像的边缘和纹理等细节信息,从而使去噪后的图像不会产生严重的模糊现象。3.2分数阶积分滤波器的构造

仍然假设二维图像信号,(髫,y)对于夏轴和y轴方向的分数阶积分在一定的条件下是可分离的,这样可以利用分数阶积分算子分别对数字图像的八个方向进行滤波处理,然后将得到的中间结果通过相关规则结合。将图像信号,(石,y)的持续期

[a,£]按单位间隔危=l等分,即n=I半l=【£一o】。这样可

r.

。1^=l

,‘

以得到G—L定义下的分数阶积分算子沿x轴和y轴方向的数值计算表达式,如式(19)和(20)所示,其中尺’表示分数阶积分数值计算的余项,如式(21)和(22)所示。

:墨’,(x,y),垒,(¥,y)+F’,(x—l,),)+型掣,(x一2,),)+…+

…1、

万方数据

!:£!:!二攀,(*一3,y)+尺:(z,y)

(19)

:c

7,(*,y),垒,(。,y)+Ⅳ,,(Ⅸ,y一1)+丛妥山,(。,y一2)”.+

矗墨(z)=!:£!!±』jj』;:二必,(*一4,,,)+…+

!:£!:±』攀,(z,),一3)+R;(z,y)

(20)

R’;’(*)=竺』』£j,(m+1)厂(口’)…’…”7

高‰毗吨y)L!上工j:÷兰必,(z,,,一4)+一(21)、…

羔揣如扩m)

・+,(m+1)厂(F’)…。’…7

(22)、…

图4表示分数阶积分去噪掩模算子形,该算子将式(19)和(20)推广到图像的其余六个方向,从而得到一个基于八个方向的去噪掩模算子(图中表格空白部分用0填充)。该去噪算子具有各向旋转不变性,去噪算子形中的滤波系数如式

(23)所示。

‰=l

觋,=”

%:华

%:止岩

心3’

%。=亟必鼍≯掣业

00嘭

0O

O0孵’

形O

O0O

形形形0

彬彬形彬

OO

矿%

0O

0矿

0形

0O

0O

形O

0O%

图4分数阶积分去噪掩模算子

4分数阶偏微分方程与图像去噪

基于偏微分方程(panial

differential

equations,PDE)的图像

处理是一门崭新的交叉学科分支,融合了数学、物理学、电子信息学和计算机科学与技术等多个领域的相关理论知识。迄今为止,国内外许多学者已经建立了大量经典的基于PDE的图像处理模型,并显示出了巨大的生命力Ⅲ“J。基于整数阶PDE的图像去噪模型具有强大的数学理论基础,因此该方法具有其他去噪方法不具备的诸多优势,但是国内外许多学者在深入研究的过程中仍然发现PDE图像去噪模型存在很多不足之处。比如P—M”u和ROF模型¨2,5纠虽然能够有效地去除平滑区域的噪声,但是去噪后的图像会产生阶梯效应而出现分片光滑现象;LLT模型"4o利用了高阶PDE建模,虽然能够消除阶梯效应,但是去噪后的图像视觉效果不好;McM去噪模型p”“1虽然保护了图像的边缘信息,但是去噪后的图像会产生许多人为痕迹。为了解决各种整数阶PDE图像去噪模型存在的问题,有学者提出将分数阶微积分理论引入到PDE图像

去噪模型中,构建基于分数阶PDE的图像去噪模型。该模型

・418・

计算机应用研究第29卷

将分数阶梯度算子替代整数阶梯度算子来构建去噪模型的能量泛函。

4.1

基于分数阶Fourfer变换的图像去噪模型

因为分数阶有界变差空间曰矿(n)是传统的有界变差空

间曰y(n)的推广和延伸,By(n)空间是建立在传统梯度模值的基础上,相应地,曰矿(n)空间就应该建立在分数阶梯度模值基础上。这样可以建立基于分数阶有界变差空间日矿(门)的能量泛函来描述噪声图像。

已知噪声图像,(省,y),(戈,y)∈力,将ROF模型的正则项梯度算子V,用分数阶梯度算子V”,来替代,可以得到基于分数阶偏微分方程去噪模型的能量泛函:

E(,)=』n1D”,ld%dy+÷I|,0一,||乞

,∈BP(力),(Ⅳ,y)∈力

(24)

其中:fD”,f-√(彤,)2+(彤,)2、D=,和彤,分别表示图像,

(戈,y)关于x和y轴的分数阶导数。

为了得到上述能量泛函的EuIe卜Lagmnge方程,需要极小化上述能量泛函,因此,利用变分法的性质构造相应能量泛函。

g(8)=E(,+占妒)=

』nID。(,+印)Ickdy+{一||,0一,一即||乙

(25)

其中:占为任意实数,妒为任意函数。

利用泛函极值的基本性质,当P∈瞄(n),有g’。(0)=0,

可以得到相应的Eul矗一Lagrange方程:

』n(尚%+尚叫岫训舶棚删删(26)

为了得到最终的去噪模型,需要消除任意函数妒。因此,

采用ParseVal恒等公式』n,‘g出妙=去』n,g幽,幽z对式

(26)进行处理,可得

知(龟9尚+匆妒尚)圳旷

—]广一—广

^万万

去fn妒(,o一,)如1曲2=o

(27)

因为分数阶Fourier变换的性质有域妒:(口俐。)”9,联妒=

(丑叮『∞:)。妒,并消除系数l/(2订)可以得出

——广—r一

』n((细∥;尚郴酬;尚)岫慨一

^百百

A』n妒(毛一,)dmld∞2=O

(28)

又因为函数妒的任意性,利用变分法的基本引理可得出、

—广—广

。【(丑俐I)”面+(尼彻2)。面一A(,0—7)卸(29)。。(椭)”尚+(胁)。尚州名一?):o(29)

一一

对式(29)两边先去共轭,然后再利用Fourier逆变换可以得出相应的图像去噪的分数阶偏微分方程:

1,‘Re{[恧(尚)+谚(尚)】卜¨瑚㈣,

最后利用梯度下降法,引入时间变量f,可以得到基于

F0谢er变换的分数阶偏微分方程的去噪模型:

詈一e{[瓦(高)吲尚)肛c¨川,・,

该模型最大的难点就是如何求取分数阶微分算子的共轭

万方数据

算子D=和联。基于Fou矗er变换的分数阶偏微分方程的去噪模型是利用分数阶Fourier变换的性质来得到彤和联的数值解。

已知有噪声图像,(z,y)(省,y∈力),假设其在x和y轴方向相互独立,其石和y方向的分数阶Fourier变换为

知(m・^’训(1一≯M”沙㈨”聊)/Ⅳ

二,,(。,y)=,(。,,,)一,(。,y一1)=

‘32’

知m-‰(-一扩咧V州…蚓“

因此,由Fourier变换的性质可知,一阶差分的空域与频域的表达式对应关系为

『D;,(并,y)H(1一e一口删17Ⅳ),(∞l‘,∞2)

(33)

【D,,(并,y)+{(1一e—J2”“27Ⅳ),(∞l,m2)

同理可知,n阶差分的空域与频域的表达式对应关系为

『D:,(z,y)H(1一e—B”。17”)“,(∞l,∞2)

(34)

【D:,(z,),)H(1一e一声删27Ⅳ)“,(∞1,∞2)

由分数阶Fourier变换的基本性质可知,将整数阶次n推广到分数阶次”,即可得到基于分数阶Fourier变换的分数阶导数数值计算公式:

fD:,(z,),)H(1一e一口”’17”)”,(∞1,∞2)

(35)

【D:,(菇,y)十+(1一e一口伽’27Ⅳ)。,(∞1。,∞2)

这样,由式(35)可以得到分数阶微分算子的共轭算子:

{————.

『D:,(算,,,)++(1一e一正删17Ⅳ)。,(∞1‘,∞2)

(36)

【D:,(省,y)++(1一e一届俐27Ⅳ)”,(∞1,∞2)

已知有噪声图像,(聋,y),戈,y∈仃,并且假设噪声图像,(石,号,(戈,・)的分数阶微积分数值计算。由分数阶微积分的G—L

定义可知,图像信号,(互,・)的分数阶微分为

刚v)=(÷)“善(-1)’而器拦而如咄…s,)

当步长^=1时可以得到

:彤如)-i邑(。)“},(一,):剃加善㈢)U卜D

(38)

假设函数c”(茄)满足

c。(。)=f(一・)5(:)

算≥。

(39)

【o

膏<o

因此,可以将信号的分数阶微积分改写为信号,(石)与函

:叫,(g)=c”(z)+,(z)

(40)

于是,将式(40)代入基于分数阶曰旷(以)空间的能量泛函,可得到基于卷积积分的分数阶变分模型:

E(,)=』n

c“+,ld聋dy+—:一I|厂一,1l乙

(41)

c”},J=/(c”。,):+(c。},);2(c”+D;=』nc。(:),(*一z,y)出

(42)

为了得到式(41)的Euler_咖ge,与上节假设一样,仍然

【(c。¥,),=/nc。(z),(x,,一z)出

4.2基于分数阶基本定义的图像去噪模型

y)在x和y轴方向相互独立,因此只考虑x轴方向的图像信数c”(髫)的卷积操作形式:

其中:

第2期

黄果,等:分数阶微积分在图像处理中的研究综述

・419・

假设任蒽函数妒,从而构造函数,便兵满足式(43):

性具有特殊的生物视觉感受野模型,因此分数阶微积分理论g(s)=E(,+P妒)=』n

c”(,+印)Ickdy+争o,o一,一8妒|I乞(43)

为数字图像底层处理提供了一种崭新的研究方法。表1表示分数阶微积分理论应用于图像底层处理的类比,该表说明了仍然利用泛函极值的基本性质令g’(0)=0,可以得到

分数阶微积分理论应用于图像处理所涉及的定义、适用情。n

f。笠塑盘I磐坐兰丛尘韭出¨c”},I

““7

况、特点、优缺点情况。

Afn(,0一,)妒dxdy=O

(44)

表1

分数阶微积分应用于图像底层处理的类比

因为当(省,y)薯力时,,(菇,),)=0。所以将伫一R2时,可以得到Parseval恒等式:

^百

』R塘出dy2寿』R2,gd∞1d∞2

(45)

————F一———百一

可以根据Parseval恒等式将式(44)进一步改写为

去,n(箐等)c一:珐如,她+去,力(箐等)㈡:坫蛔幽:一

A寿』n(,0一,)妒d∞1d∞220

(46’

因为存在c”},=c。・,,并且因为函数妒的任意性,可以得到

(鬻)州‰(锗)以oAc名一2Ⅲ㈤,—r

本文介绍的各种基于分数阶微积分和分数阶偏微分方程因为c”(∞)=c”(一∞),可以得到

理论的图像处理模型还存在一些不足之处,在今后的研究工作(锗)以二训+(锗)以二㈦州J2一?)=o㈤,

中需要重点解决以下几个方面的问题:a)本文主要介绍分数阶微积分理论应用到图像去噪和增强两种基本的图像底层处利用F0urier逆变换的性质可以得到

((锗)甜㈠,)。+((锗)"㈠,),一

理,需要进一步将基于分数阶微积分的图像处理模型推广到图像分割、边缘提取、数字水印、指纹提取等图像处理的其他应用领域;b)本文介绍的基于分数阶微积分理论的图像处理模型A(如一,)=O

(49)

都只是应用于二维图像处理,因此,以后要逐步尝试将上述模慷沁)锗㈦埘胁协沁)锗㈠"胁

去掉卷积符号后可以得到

型推广到三维图像处理;c)分数阶微积分理论作为一种基本的数学分析工具,已经成功应用到了图像处理中,并取得了很A(,0一,)=0

(50)

好的应用效果,如何将分数阶微分理论和当今成熟的经典算法最后,利用梯度下降法即可以得到基于G-L定义分数阶结合,如利用神经网络、遗传算法、模糊数学等经典理论来估计偏微分方程的图像去噪模型:

分数阶微积分的阶次,从而构建性能更加完备的图像处理算詈=一(协沁)锗㈦刎胁

子;d)图像去噪和图像增强是互为逆过程,因此,可以尝试将基于分数阶偏微分方程的图像去噪模型利用反扩散原理推广』nc”(z)!i;÷等(x,,,+z)出)+A(,0一,)

得到基于分数阶偏微分方程的图像增强模型,并将该模型和其

(51)

他的一些经典图像增强方法结合得到性能更加优越的图像增强模型;e)本文介绍关于分数阶偏微分方程的图像处理模型5结束语

都是基于较弱函数空间条件下建立的,在数学完备性方面存在自从Mandelbmt发现分形现象,并将Riemann-LiouviⅡe分不足,因此需要寻找更为适合的基于分数阶微积分理论的图像数阶微积分理论应用于分析和研究分形媒介中的布朗运动描述空间对图像建模。

以来,分数阶微积分逐渐受到工程技术人员的高度关注。近分数阶微积分理论作为一种崭新的工程应用工具,已经引些年以来,分数阶微积分理论已经在化学、电磁学、控制学、起了国内外众多科技工作者的广泛关注,并在包括图像处理的材料科学和力学等众多学科领域中得到广泛应用并取得了实际应用领域中取得了较好的成果。因此可以大胆预测,随着重大的研究成果。然而,将分数阶微积分理论应用到现代信分数阶微积分理论的进一步发展和完善,数值计算方法和滤波号分析与处理中,特别是数字图像的处理,在国内外都还是器构造方法的进一步简化和改进,学术前沿的新方法和新技术一个新兴的学科分支。分数阶微积分理论适合研究具有非的进一步推广和应用,计算机硬件设备运算能力的进一步发展线性、非因果、非平稳等特征不确定的信号,因此分数阶微积

和提高,分数阶微积分理论未来必将具有更加广泛的应用分理论中许多经典算法都非常适合在现代信号分析与处理

空间。中应用。众所周知,数字图像中邻域内像素点之间的灰度值参考文献:

具有高度的自相似性,这种高度的自相似性可以近似为分形[1]

LEIBNIzG

w.corre叩ondaticedeLeibniz

avec

Hugens,vaJlzuli.

结构,并通常是以复杂的边缘和纹理等细节信息来表示。此ch帆et

leMarquisdeL’H0spital[M].P耐s:Libr.de

A.F咖ck,

外,数字图像的分数阶微分具有特殊的马赫现象,其拮抗特

1853:297—302.

万方数据

・420・计算机应用研究

第29卷

『2]

MANDELBROTBB,vanNESsJw.ThefhctionalbrDunianmo.

tions,fraction8lnoisesand

applications[J].SIAM

ReView,1968,

10(4):422—437.

[3]MANDELBROTBB.7IhefmctaJ学eometlyofnature[M].New

York:

W.H.Freem粕andCompany,1983.[4]袁晓,张红雨,虞厥邦.分数导数与数字微分器设计[J].电子

学报,2004,32(10):1658—1665.

[5]阿然,齐林,王越.分数阶Fourier变换的原理与应用[M].北京:

清华大学出版社,2004:23—49.

[6]陶然,张峰,王越.分数阶Fourier变换离散化的研究进展[J].中

国科学E辑:信息科学,2008,38(4):481.503.

[7]薛定宇,赵春娜.分数阶系统的分数阶PID控制器设计[J].控

制理论与应用,2007,24(5):771—776.[8]

PODLUBNYI.Fractional.ordersvstemsaJld

controllers[J].1EEE

T陋ns

onAutomatic

ControI,1999,44(1):208—214.

[9]zAYEDAI.OntherelationshipbetweentheFburiertransfo唧and

fhctionalFourier

t啪sfom[J].1EEE

SignaIProcessLett,1996,3

(12):310—311.

[10]廖科.分数阶微积分运算数字滤波器设计与电路实现及其应用

[D].成都:四川大学,2006.[11]RODlEcK

w.Quantitativeanalysisof

catretinal

98n对ion

cell

re—

spon8e

to

vi8ual

stimuli[J].VisionResearch,1965,5(11):583—

601.

[12]蒲亦非,袁晓,廖科,等.一种实现任意分数阶神经型脉冲振荡器

的格形模拟分抗电路[J].四川大学学报:工程科学版,2006,38

(1):128一132.

[13]陈庆利,蒲亦非,黄果,等.分数阶神经型脉冲振荡器[J].四川

大学学报:工程科学版,2011,43(1):123—128.

[14]蒲亦非,王卫星.数字图像的分数阶微分掩模及其数值运算规则

[J].自动化学报,2007,33(11):1128-1135.

[15]PUYi—fei,WANGwei.】【ing,zHOuJi-liu,封口z.Fmctional

differen—

tial印proach

to

detectingtextur■featuresofdigitalimageanditsfhc—

tionaldiⅡerentialfilterimplementation[J].

Scjence-nChinaSe—

riesF:Inforfllation

Sciences,2008,51(9):1319一1339.

[16]蒲亦非.将分数阶微分演算引入数字图像处理[J].四川大学学

报:工程科学版,2007,39(3):124—132.

[17]蒲亦非.分数阶微积分在现代信号分析与处理中的应用的研究

ID].成都:四川大学,2006.

[18]、蒲亦非,王卫星,周激流,等.数字图像纹理细节的分数阶微分

检测及其分数阶微分滤波器实现[J].中国科学E辑:信息科学,

、2008,38(12).:2252—2272.

[19],蒲亦非.将分数阶微分演算引入数字图像处理[J].四川大学学

报:工程科学版,2007,39(3):124—132.[20]Pu

Yi—fei,zH0u

Ji—liu,yuAN

xiao.FractionaldiffeIentialmask:a

fractionaI

difIbrential-based如proach

fbrmultiscale

textu南enhance-ment[J].IEEE

Trans

on

Image

Proeessing,2010,19(2):491一

.一{

511.

[21]黄果,蒲亦非,陈庆利,等.非整数步长的分数阶微分滤波器在图

像增强中的应用[J]..四川大学学报:工程科学版,2011,43(1):

129一136.

[龙q,黄果,蒲亦非,陈庆利,等.纂于分数阶积分的图像去噪[".+系统

工程与电子技术,2011,33(4):925—932.

[23]杨柱中,周激流,黄梅,等.基于分数阶微分的边缘检测[J].四川

大学学报:工程科学版,2008,40(1):152-157.

[24]杨柱中,周教流,晏祥玉,等.基于分数阶微分的图像增强[J].计

算机辅助设计与图形学学报。2008,20(3):343—348.

万方数据

[25]杨柱中,周激流,黄梅,等.用分数阶微分提取图像边缘[J].计算

机工程与应用,2007,43(35):15-18.[26]MATHIEU

B,MELcHl0RP,0uTSALOuPA,doZ.Fmctionaldif-

ferentiationfbredgedetec60n[J].SignaIPorcessing,2003,83

(11):2421—2432.

[27]李远禄,于盛林.分数阶差分及边缘检测[J].光电工程,2006,33

(12):70一74.

[28]汪凯宇,肖亮,韦志辉.基于图像整体变分和分数阶奇异性提取

的图像恢复模型[J].南京理工大学学报:自然科学版,2003,27

(4):400-4J04.

[29]刘红毅,韦志辉.基于分数阶样条小波与IHs变换的图像融合

[J].南京理工大学学报:自然科学版,2006,30(1):82—84.[30]LIuJun,cHEN

song—can,TANxiao—yang.Fractionalordersingular

value

deco‘叩ositionrepresentationforfacereco印ition[J].

Pattem

Recognition,2007,41(1):168-182.

[31]左凯,孙同景,李振华,等.二维分数阶卡尔曼滤波及其在图像处

理中的应用[J].电子与信息学报,2010,32(12):3027—3031.[32]汪成亮,兰利彬,周尚波.自适应分数阶微分在图像纹理增强中的

应用[J].重庆大学学报,20ll,34(2):32—37.

[33]高朝邦,周激流.基于四元数分数阶方向微分的图像增强[J].自

动化学报,2011,37(2):150一159.[34]BAIJian,FENG

xiang—chu.Fractional-orderanisotmpicdi饪hsionfbr

image

denoising[J].IEEE

Trans

on

ImageProcessing,2007,16

(10):2492-2502.

[35]张军,韦志辉.sAR图像去噪的分数阶多尺度变分PDE模型及

自适应算法[J].电子与信息学报,2010,32(7):1654—1659.[36]张军,韦志辉.一种基于卷积积分的图像去噪变分方法[J].中国

图象图形学报,2008,13(9):1673二1677.

[37]张军,韦志辉.分数阶图像去噪变分模型及投影算法[J].计算

机工程与应用,2009,45(5):1—6.

[38]张军.基于分数阶变分的图像建模与去噪算法研究[D].南京:

南京理工大学,2009.[39]zHANG

Jun,w_EIzhi—hui.AclassoffbcⅡonal-order舢lti—scale

v耐ational

modelsand

altematingprojectionalgofithmforimagede—

noising[J].App¨ed

MathematicaI

ModeⅢng,20“,35(5):2516—

2528.

[40]沈淑君.分数阶微分方程的数值解方法及误差分析[D].厦门:厦

门大学,2005.

[41]庄平辉,刘发旺.空间一时间分数阶扩散方程的显式差分近似

[J].高等学校计算数学学报,2005,27(51):223-228.[42]MAINARDI

F,LucHKOY,PAGNINIG,PARADIsIG.Thefunda-

mentalsolutionofthespace—timefractionaldiffusion[J].

FracCalc

AppI

Anal,2001,4(2):153—192.

[43]HuANGFeng—hui,uu

Fa-wang.The

fund砌ental

solutionofthe

space—tirrlefhctionaladvection1dispersionequation[J].

JournaIof

App¨edMathematics&Compu“ng,2005,18(1—2):339-350.[44]鸭cHuMPE砒正D,DERICHE

R.Vector—valuedimageregIllariza一

“.tion:whhPDEs:。acDmmon。fI’舭单Iework

fbrdi]日ferent

applications[J].

。IEEETrans

On

Pattern

AnaIysisandMachineInte¨igence,

2005,27(4):506—517.[45]GILBOA

G,socHENN,zEEVIYY,Variationaldenoisingofpar【一

ly

texturedimagesbyspatiaIlyvaryi“g

constmints[J].

JEEETrans

on

Image

Processing,2006,15(8):228l一2289.

[46]朱立新,王平安,夏德深.引入耦合梯度保真项的非线性扩散图

像去噪方法[J].计算机研究与发展,2咖,44(8):1390-1398.

(下转第426页)

・426・

tionwith

计算机应用研究

mobile

第29卷

on

phone:e%ctsoftouchkeysizeandt叫chkeyloca.

pleasure

theIntemetshoppingexpe^ence[J].JournaI

ofRetai-

tion[J].IntematjonaI

(1):68-76.

JournaIofIndustriaIErgonomics,20lO,40Iing,2002,78(1):31—40.[30]FRANKA

M,cHRISTOPHM.Mobile

Webusabilityevaluation—com.

[19]CHANG

dren

on

on

C.u8abilitycomparisonofpen-basedi“putforyoungchil—mobile

biningthemodifiedthinkaloudmethodwiththetestingofemotional,

devices[c]//Proc

ofIEEEIntemationalCor迁erence

co邸itiVeandconatiVeaspects0ftheusageof

WebapPlication

on

SensorNetworks,Ubiquitous,andTmstworthyComputing.2008:[C]//Pmc

ofthe9t}I

IEEE/AClsIntemationalCoIlference

Com.

531.536.

puterandInfoHnation

s.Mobile

Science.Washin群onDC:IEEEComputerSo—

[20]sHREsTHAwebbrowsing:usabiIitystudy[c]//Procof

on

ciety,2010:367・372.

the4thIntemationalConferenceandSystemsandthe1

st

MobileTechnology,Applications,

on

[31]BANDuRA

nitive

A.Socialfoundationsofthoughtandaction:asocialco哥

clif!fj:Prentice—HallInc,1986.

IntemationalSymposiumComputerHumantheory[M].En毋ewood

IntemctioninMobileTechnoIogy.NewYork:ACMPress,2007:187一194.

[32]s‘I’APLEsD,HuLLANDJ,HIGGINSC.A

self-emcacytheoIy

ex.

planationforthemanagement0fremoteworke玛invinual

H,KIM

o唧i豫tions

[21]KIMJ,LEEY.Anempiricalstudyof

on

u8econtextsinthe

[J].

JoumaIof

Computer-MediatedCommunlcation,1998,3

mobileIntemet,focusingtheusabilityofinf0Ⅱnation

architecture

(4):758—776.

[J].Jnformation[22]Lu

Systems

Front_erS,2005,7(2):175—186.

on

[33]HsU

M,cHIu

c.Intemetself-emcacyaJldelectronicserviceaccep.

H,suP.Factorsafbctingpurchaseintentionmobileshopping

tance[J].DecisionSupportSystems,2004,38(3):369.381.[34]LuARN

tion

to

websites[J].IntemetResearch,2009,19(4):442—458.

P,uNH.T0ward

mobile

an

understanding0fthebehavioralinten—Computersin

HumanBehavior,

[23]HALLR,HANNAP.Theimpact0f

on

webpagetext-backgroundc010r

aestheticsandbehavioraIin—

use

baIIking[J].

combinations

readability,Tetention,2005,21(6):873—891.

tention[J].BehaVior&InformationTechnoIogy,2004,23(3):

183—195.

[35]TslA0usIsAs,GIAGLIsGM.Anempiricalassessmentofenvimn—

theusabilityof

on

mentalfactorsthatinnuence

sK,HAsSANMA.AsmaItframeworkfor

mobile

website[c]//

Wash.

[24]GuIRGuIS

and

web

content

Proc0fthe9山IntemationalConferenceMobileBusiness.

resources

ad印tationinmobiledevices[c]//Procofthe12thIn・

0n

in昏onDC:IEEECo。叩uter

Society,2010:161—167.cmwding’sinnuence

on

temationalConferenceAdv蛐cedCommunication

Technology.

[36]sE0s,JuNGs.

Per℃eivedcon8umer’s

20lO:487—492.

shopping

behavior[J].KoreanJoumalofDist—buaonResearch,

[25]KAIKKONENA.usabilityproblemsintoday’smobileofthe2nd

Intemetpor-

Mobile7rech—

2005(1):23-32.

tals[c]//Proc

Intemationalcorlference

on

[37]scHMIDT

than

A,BEIGLM,GELLERsENH.7Ihereismore

tocontext

nology,ApplicationsaIldSystems.2005:1—7.location[J].Computer&Graphics,1999,23(6):893.901,

D,PARKs,以叫.AframeworkforevallIatingthe

[26]AcTON

T,GOLDEN

w,GuDEAs,越of.usab订ityaIldacceptance

ofthe9th

[38]HE0J,HAMB

insmall—scTeen

inf0啪ationsystems[c]//PmcEur叩eaIlusabilityofmobileusability263.275.

phones

based

on

multi—level,hieraI℃hicalHlodelofwith

CollaborativeElectronicrence.2004.

CommerceTechnologyand

Research

Co出-

factors[J].InteractingComputerS,2009,21(4):

[27]NAH

F.A

to

study

on

tolerablewaitingtime:howIong

are

web

users

[39]LESEMANN

luatemobile

on

E,NATALIE

w,SVENK.Combiningmethods

toeva—

willi“g

wait[J].BehaVior&|nformationTechnoIogy,2004,23

usability[c]//Pmc

ofthe9thIntemationalconference

M0bile

DevicesaJld

Services.

(3):153一163.HumanComputerIntemctionwith

[28]TARAsEwIcH

P.

MobilecommeI℃eoppor【unitiesandchallenges:

2007:444—447.

designingmobilecommercethe

applications[J].

Communicationsof

[40]DuHH,TANG,cHENV.usabilityevaluationformobiledevice:a

0f

laboratory

aJld

field

ACM,2003,46(12):1・60.

s,KAHNB.cross・categorydfects

ofinducedamusaland

comp面son

tests[J].

JoumaIof

UsabⅢty

[29]MENONStudjes,2005,1(I):4—16.

268.

(上接第420页)

[47]suNxiao—li,wANG

Jun—ping,SONGGuo—xiang.Anisotropicdifh-

[53]0sHERS,RuDINLI,FA7rEMIE.Nonlineartotalvariationbased

sionfilteringmethodwithand

fidelity

tem[J].SystemsEngineering

noiseremoval

algorithms[J].Physica

D,NonIinearPhenomena,

EIectronics,2008,30(8):1421—1423.1992,60(3):259-268.

[48]陈利霞,冯象初,王卫卫,等.加权变分的图像去噪算法[J].系

统工程与电子技术,2010,32(2):392—395.[49]u£E

SH,sE0

[54]LYsAKERM,LUNDERVOLDA,TAI

xc.

Noiseremovalusing

to

fourth—orderpanialdifferentialequationwith印plication8

medical

on

JK.NoiseremovalwithGausscurvature.drivendifhJ—

Trans

on

ma驴eticreson锄ceimagesinspaceandtime[J].ImageProcessing,2003,12(12):1579-1590.

1EEETrans

sion[J].IEEE

909.

Image

Processing,2005,14(7):904—

[55]KAssM,wITKIN

models[J].

A,TERz0POuL0sD.snakes:active

contour

[50]王新楼,乔明,邹谋炎.一种基于偏微分方程的sAR图像去噪方

法[J].电子与信息学报,2005,27(9):1365—1369.[51]PERONA

P,MALIK

InternationaIJournaIof

ComputerVision,1988,1

(4):321—331.

J.scale—spaceandedgedetectionusinganiso—

IEEE

Trans

on

[56]xuchen—yallg,JERRY

[J].1EEE[57]xu

Trans

on

LP.snal【es

sh印esandgmI|ientvect叫flow

tI、叩icdi】6fusion[J].

chine

PatternAnaIysisandMa-

lmage

Processing,1998,7(3):359—369.

vector

Inte¨igence,1990,12(7):629.639.

L,0sHERs,FA7IEMIE.。Nonlineartotalvariationbased

chen—ya“g,JERRYLP.Generalized目mdient

nowexte卜

[52]RuDING

nalforcesforactive131一】39.

contours[J].SignalProcessing,1998,71(2):

noisel-emoval

algorithms[J].PhysicaID,1992,60(1—4):259—

万方数据

分数阶微积分在图像处理中的研究综述

作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:

黄果, 许黎, 蒲亦非, HUANG Guo, XU Li, PU Yi-fei

黄果,HUANG Guo(乐山师范学院智能信息处理及应用实验室,四川乐山,614000), 许黎,XU Li(乐山师范学院物电学院,四川乐山,614000), 蒲亦非,PU Yi-fei(四川大学计算机学院,成都,610064)计算机应用研究

Application Research of Computers2012,29(2)4次

参考文献(57条)

1. LEIBNIZ G W Correspondatice de Leibniz avec Hugens,van Zulichem et le Marquis de L' Hospital 18532. MANDELBROT B B;Van NESS J W The fractional brounian motions,fractional noises and applications 1968(04)3. MANDELBROT B B The fractal geometry of nature 1983

4. 袁晓;张红雨;虞厥邦 分数导数与数字微分器设计[期刊论文]-电子学报 2004(10)5. 陶然;齐林;王越 分数阶Fourier变换的原理与应用 2004

6. 陶然;张峰;王越 分数阶Fourier变换离散化的研究进展[期刊论文]-中国科学E辑 2008(04)7. 薛定宇;赵春娜 分数阶系统的分数阶PID控制器设计[期刊论文]-控制理论与应用 2007(05)8. PODLUBNY I Fractional-order systems and controllers[外文期刊] 1999(01)

9. ZAYED A I On the relationship between the Fourier transform and fractional Fourier transform 1996(12)10. 廖科 分数阶微积分运算数字滤波器设计与电路实现及其应用[学位论文] 2006

11. RODIECK R W Quantitative analysis of cat retinal ganglion cell response to visual stimuli 1965(11)12. 蒲亦非;袁晓;廖科 一种实现任意分数阶神经型脉冲振荡器的格形模拟分抗电路[期刊论文]-四川大学学报(工程科学版)2006(01)

13. 陈庆利;蒲亦非;黄果 分数阶神经型脉冲振荡器[期刊论文]-四川大学学报(工程科学版) 2011(01)14. 蒲亦非;王卫星 数字图像的分数阶微分掩模及其数值运算规则[期刊论文]-自动化学报 2007(11)

15. PU Yi-fei;WANG Wei-xing;ZHOU Ji-liu Fractional differential approach to detecting textural features of digitalimage and its fractional differential filter implementation[期刊论文]-Science in China Series F:InformationSciences 2008(09)

16. 蒲亦非 将分数阶微分演算引入数字图像处理[期刊论文]-四川大学学报(工程科学版) 2007(03)17. 蒲亦非 分数阶微积分在现代信号分析与处理中的应用的研究[学位论文] 2006

18. 蒲亦非;王卫星;周激流 数字图像纹理细节的分数阶微分检测及其分数阶微分滤波器实现[期刊论文]-中国科学E辑 2008(12)19. 蒲亦非 将分数阶微分演算引入数字图像处理[期刊论文]-四川大学学报(工程科学版) 2007(03)

20. PU Yi-fei;ZHOU Ji-liu;YUAN Xiao Fractional differential mask:a fractional differential-based approach formultiscale texture enhancement 2010(02)

21. 黄果;蒲亦非;陈庆利 非整数步长的分数阶微分滤波器在图像增强中的应用[期刊论文]-四川大学学报(工程科学版) 2011(01)22. 黄果;蒲亦非;陈庆利 基于分数阶积分的图像去噪[期刊论文]-系统工程与电子技术 2011(04)23. 杨柱中;周激流;黄梅 基于分数阶微分的边缘检测[期刊论文]-四川大学学报(工程科学版) 2008(01)24. 杨柱中;周激流;晏祥玉 基于分数阶微分的图像增强[期刊论文]-计算机辅助设计与图形学学报 2008(03)25. 杨柱中;周激流;黄梅 用分数阶微分提取图像边缘[期刊论文]-计算机工程与应用 2007(35)

26. MATHIEU B;MELCHIOR P;OUTSALOUP A Fractional differentiation for edge detection[外文期刊] 2003(11)27. 李远禄;于盛林 分数阶差分及边缘检测[期刊论文]-光电工程 2006(12)

28. 汪凯宇;肖亮;韦志辉 基于图像整体变分和分数阶奇异性提取的图像恢复模型[期刊论文]-南京理工大学学报(自然科学版)2003(04)

29. 刘红毅;韦志辉 基于分数阶样条小波与IHS变换的图像融合[期刊论文]-南京理工大学学报(自然科学版) 2006(01)30. LIU Jun;CHEN Song-can;TAN Xiao-yang Fractional order singular value decomposition representation for facerecognition 2007(01)

31. 左凯;孙同景;李振华 二维分数阶卡尔曼滤波及其在图像处理中的应用[期刊论文]-电子与信息学报 2010(12)32. 汪成亮;兰利彬;周尚波 自适应分数阶微分在图像纹理增强中的应用[期刊论文]-重庆大学学报 2011(02)33. 高朝邦;周激流 基于四元数分数阶方向微分的图像增强[期刊论文]-自动化学报 2011(02)

34. BAI Jian;FENG Xiang-chu Fractional-order anisotropic diffusion for image denoising[外文期刊] 2007(10)35. 张军;韦志辉 SAR图像去噪的分数阶多尺度变分PDE模型及自适应算法[期刊论文]-电子与信息学报 2010(07)36. 张军;韦志辉 一种基于卷积积分的图像去噪变分方法[期刊论文]-中国图象图形学报 2008(09)37. 张军;韦志辉 分数阶图像去噪变分模型及投影算法[期刊论文]-计算机工程与应用 2009(05)38. 张军 基于分数阶变分的图像建模与去噪算法研究[学位论文] 2009

39. ZHANG Jun;WEI Zhi-hui A class of fractional-order multi-scale variational models and alternating projectionalgorithm for image denoising 2011(05)

40. 沈淑君 分数阶微分方程的数值解方法及误差分析 2005

41. 庄平辉;刘发旺 空间-时间分数阶扩散方程的显式差分近似 2005(51)

42. MAINARDI F;LUCHKO Y;PAGNINI G;PARADISI G The fundamental solution of the space-time fractional diffusion2001(02)

43. HUANG Feng-hui;LIU Fa-wang The fundamental solution of the space-time fractional advection dispersion equation 2005(1-2)

44. TSCHUMPERLE D;DERICHE R Vector-valued image regularization with PDEs:a common framework for differentapplications [外文期刊] 2005(04)

45. GILBOA G;SOCHEN N;ZEEVI Y V Variational denoising of partly textured images by spatially varying constraints[外文期刊] 2006(08)

46. 朱立新;王平安;夏德深 引入耦合梯度保真项的非线性扩散图像去噪方法[期刊论文]-计算机研究与发展 2007(08)

47. SUN Xiao-li;WANG Jun-ping;SONG Guo-xiang Anisotropic diffusion filtering method with a fidelity term[期刊论文]

-Journal of Systems Engineering and Electronics 2008(08)

48. 陈利霞;冯象初;王卫卫 加权变分的图像去噪算法[期刊论文]-系统工程与电子技术 2010(02)

49. LEE S H;SEO J K Noise removal with Gauss curvature-driven diffusion[外文期刊] 2005(07)50. 王新楼;乔明;邹谋炎 一种基于偏微分方程的SAR图像去噪方法[期刊论文]-电子与信息学报 2005(09)

51. PERONA P;MALIK J Scale-space and edge detection using anisotropic diffusion[外文期刊] 1990(07)52. RUDING L;OSHER S;FATEMI E Nonlinear total variation based noise removal algorithms 1992(1-4)53. OSHER S;RUDIN L I;FATEMI E Nonlinear total variation based noise removal algorithms 1992(03)

54. LYSAKER M;LUNDERVOLD A;TAI X C Noise removal using fourth-order partial differential equation withapplications to medical magnetic resonance images in space and time[外文期刊] 2003(12)55. KASS M;WITKIN A;TERZOPOULOS D Snakes:active contour models 1988(04)56. XU Chen-yang;JERRY L P Snakes shapes and gradient vector flow 1998(03)

57. XU Chen-yang;JERRY L P Generalized gradient vector flow external forces for active contours[外文期刊] 1998(02)

引证文献(4条)

1. 勾荣 基于Riemann-Liouville分数阶微分的图像增强[期刊论文]-制造业自动化

2013(16)2. 勾荣 基于 G-L 分数阶微分的图像边缘检测[期刊论文]-计算机与现代化

2013(11)3. 勾荣 基于分数阶微分的图像增强算法[期刊论文]-电子科技 2013(12)

4. 魏文力. 李宝树. 崔克彬. 徐雪涛 基于改进的最大熵航拍输电线路图像复原方法[期刊论文]-科学技术与工程 2013(32)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_jsjyyyj201202003.aspx

第29卷第2期计算机应用研究

V01.29No.2

2012年2月

ApplicationResearchofComputers

Feb.2012

分数阶微积分在图像处理中的研究综述

黄果h,许黎岫,蒲亦非2

(1.乐山师范学院a.智能信息处理及应用实验室;b.物电学院,四川乐山614000;2.四川大学计算机学院,成

都610064)

摘要:综述了关于分数阶微积分理论在数字图像底层处理中的应用研究,具体包括:分数阶微积分和分数阶偏微分方程的基本理论及分数阶傅里叶变换的基本性质;分数阶微分滤波器的构造及在图像增强中的应用研究;分数阶积分滤波器的构造及在图像去噪中的应用研究;分数阶偏微分方程在图像去噪中的应用研究。最后,总结了分数阶微积分理论在图像处理中已取得的研究成果,并结合已有的基于分数阶微积分理论的图像底层处理模型,展望了分数阶微积分理论在图像处理中的应用前景。

关键词:分数阶微积分;分数阶偏微分方程;变分法;泛函极值;图像增强;图像去噪中图分类号:TN391

文献标志码:A

文章编号:1001.3695(2叭2)02—0414—07

doi:10.3969/j.issn.100l一3695.2012.02.003

Summaryofresearch

on

imageprocessingusingfractionalcalculus

HUANGGu01a,XULilb,PUYi.fei2

(1.Ⅱk60m£o,y矿胁扰咖眦聊备玎№砌nProc哪堍&A一池t如n,6.s如00f0,P蛳如&觑∞打口n如,肠^口nⅣ0r,mz‰妇m毋,k^口n_si-如Mn614000,醌iM;2.co朋e酽o,compM把r&如袱,.s如^∽n‰眈巧蚵,蕊e,则M610064,醌i胍)

Abstmct:The

maincontents

ofthispaperwereaboutfractionalcalculusapplicationindigitalimagelow—levelprocessing.

SpecificaⅡyincludedthefollowing:fractionalcalculusaJldfmctional—orderpartialdifEbrentialequationsandthebasicprope卜tiesoffractionalFouriertransfo册:thestmctureofthefractionaldifferentialfilteranditsapplicationtoimageenhancement;

thestmctureofthefractionalintegmtionfilteranditsapplication

to

image

denoising;theimagedenoisingmodelbased

on

fhc—

tionalpartialdi如rentialequations.FinaJly,summarizedfractionalcalculushadbeenmaderesea托hresultsandrelatedpmb—lemsinimagepmcessing.Associatedwiththeimageprocessingmodels,forecastedtheapplicationsprospectaboutfhctionalcalculusappliedtoimageprocessing.

Keywords:fhctionalcalculus;fractional-orderpartialdi矗Ierentialequations;calculus

ofvariations;functionalextreme;

image

enhancement;imagedenoising

生物、金融、流体等领域应用越来越广泛,实用价值也得以充

引言

分的体现。近几年以来,国内外许多学者研究发现,分数阶分数阶微积分理论与整数阶微积分理论几乎同时诞生,

微积分在信号分析与处理领域具有广泛的应用前景,比如在该理论最早出现在L’H6pital写给Leibniz的信件中…。近三信号的奇异性检测和提取方面具有特殊的作用H“…。分数阶微积分理论在图像底层处理中的应用是最近几年才引起百年来,分数阶微积分仅仅是数学家们的纯数学理论分析和学者关注的,到目前为止,已经取得了初步研究成果,主要包推导。Liouville等著名数学家一直致力于完善和发展分数阶括图像增强、图像去噪、图像边缘提取、图像分割、图像奇异微积分的理论研究,从而构建了分数阶微积分理论的基本系性检测等。图像处理中的整数阶微分算子包括梯度和La—统架构。直到19世纪后期,分数阶微积分理论才开始逐渐在place算子,以及墨西哥帽算子等很多经典的图像处理算子都实际工程领域中得到初步应用。Mandelbmt旧’31发现了分形可以近似看做是Rodieck…1提出的DOG感受野模型在理论学,他指出自然界和很多技术科学中存在大量分数维事实,上的仿真实现。分数阶微分算子则可以认为是DOG感受野发现在整体与部分之间存在自相似现象,并尝试将Riemann-模型的推广和延伸。基于分数阶微分理论的感受野模型较Liouville分数阶微积分理论应用于描述和分析分形媒介中的整数阶微积分理论更符合人类视觉特性。

布朗运动,从此,工程技术人员开始关注分数阶微积分。近分数阶微积分理论在图像处理领域的应用研究主要包括几十年以来,国内外许多学者研究发现分数阶微积分算子具以下方面:蒲亦非等人¨2。2o认为,选择适当阶次的分数阶微分有记忆性和非局部性,非常适合用来描述现实世界中具有记算子在增强图像过程中可以大幅提升边缘和纹理细节的同时,忆和遗传性质的材料,因此分数阶微积分理论在物理、化学、

非线性地保留了图像平滑区域的纹理信息;此外,杨柱中等

收稿日期:201l-09・16;修回日期:2011・lO一30

基金项目:国家自然科学基金资助项目(60972131);四川省教育厅科研项目(10zC019/

11zBl32);乐山市科技局重点研究计划资助项目(2011GzD046);乐山师范学院资助项目(z1162)

作者简介:黄果(1980-),男,重庆人,讲师,博士,主要研究方向为分数阶微积分及分数阶偏微分方程在图像处理中的应用研究(huangguoxuli@163.com);许黎(1982一),女,云南人,讲师,硕士,主要研究方向为分数阶微积分及混沌理论在图像处理中的应用研究;蒲亦非(1976一),男,四川人,副教授,博士(后),主要研究方向为分数阶微积分及分数阶偏微分方程在图像处理中的应用研究.

万方数据

第2期黄果,等:分数阶微积分在图像处理中的研究综述

・415・

人。2”驯利用分数阶微分算子具有弱导数性质的特点,利用分数阶梯度算子对含弱噪声的图像进行边缘检测,该方法能够有效地解决整数阶梯度算子对噪声敏感的问题,因此可以避免噪声的影响而准确地定位噪声图像的边缘;MatIlieu等人心叫提出了分数阶微分的边缘检测算子,即分数阶鲁棒轮廓(contour

r0.

busted’ordre

non

entier,cRONE)边缘检测器,说明了该检测

器的分数阶微分阶次在1<”<2时,能够有选择地检测出边缘,而在分数阶微分阶次在一1<”<1时,该检测算子能够在边缘提取的过程中克服噪声的影响;在此基础上,李远禄等人旧川提出了基于分数阶差分的滤波器并将其应用到边缘检测,该模型将数据的平滑和差分有机结合,利用分数阶微积分基本理论,设计了可以通过调节微分阶次的差分滤波器,该滤波器可以解决传统算子边缘检测出现的边缘漂移问题,并且可以抑制部分噪声;汪凯宇和刘红毅等人。28’2引利用分数阶微积分具有的记忆特性,分别将分数阶样条小波应用到图像纹理的奇异性检查和图像融合中,较整数阶微积分取得了更好的仿真效果;Liu等人∞叫提出了分数阶奇异值分解的人脸识别方法,该方法可以有效缓解面部的变化,并且在脸部出现剧烈变化时,比传统的分类方法提高其分类的性能,为高层图像处理打下坚实基础;左凯等人口1。将卡尔曼滤波和分数阶微积分理论相结合,提出了二维分数阶卡尔曼滤波器并成功应用到图像处理;汪成亮等人¨2。通过研究分数阶微积分的基本定义和相应的分数阶微分算子的实现方法,提出了将图像的梯度特征和人类视觉特征等理论引入到已有的分数阶微分算子中,由此构建了基于分数阶微分阶次自适应变化的图像增强模型;高朝邦等人”引将四元数理论和分数阶微积分理论有机结合,提出了四元数分数阶方向微分的概念,解释其物理意义和几何意义,并根据分数阶微分的特殊性质将其应用到图像增强当中,可以得到很好的效果;Bai等人旧4o以ROF去噪模型为基础,将分数阶微积分理论和偏微分方程相结合,提出了基于分数阶偏微分方程的图像去噪模型,该方法可以解决传统低阶次整数阶偏微分方程去噪模型容易产生阶梯效应,以及高阶次整数阶偏微分方程去噪模型去噪效果不佳的缺点。此后,张军等人"”驯将负指数sobolev空间的多尺度图像建模与基于分数阶微积分的图像建模有效结合,提出了统一的基于分数阶多尺度变分图像去噪模型,并初步设计了该模型涉及参数的自适应选择方法。

相关理论及分析

1.1

分数阶微积分的基本定义

从不同的应用角度去分析问题可以得到不同的分数阶微

积分定义,所以迄今为止,分数阶微积分仍然没有统一的时域定义表达式。在众多的定义表达中,主要有三种经典的分数阶微积分定义,包括Gmmwald—Letnikov、Riemann.IJiouville和ca.potu定义。

1)分数阶微积分的Grnmwald—Letnikov定义

:c=2骂∥j叁黼t一曲)weR一

(1)

2)分数阶微积分的Riemann—Liouville定义

『警

【专赢可,:茄ay。一㈨如

…eⅣ

万方数据

3)分数阶微积分的C印uto定义

细撇卜k志乃筘…川…。@’。

【鲁币%J:带如嚷”k…

f掣

…eⅣ

针对三种经典分数阶微积分定义,国内外众多学者从不同应用角度研究分析各种各样的实际问题,得到了不同形式的分数阶微积分定义表达式。分数阶微积分Riemann一“ouville定义和分数阶微积分c印uto定义都是对分数阶微积分Gnimwald.Letnikov定义的改进。分数阶微积分的Grumwald.Letnikov定义在数值实现时可以转换为卷积运算形式,因此非常适合在信号处理中的应用;分数阶微积分的Riemann—Li.ouviue定义主要应用于计算一些较为简单函数的解析解;分数阶微积分的c印uto定义适用于分数阶微分方程初边值问题的分析,因此非常适合在工程领域中应用。三种关于分数阶微积分定义之间存在着紧密的联系,在一定条件下可以相互转换。例如,如果分数阶微积分阶次w满足n一,1<”<凡,当函数,(x)的m+1阶导数连续并且满足m=L

n一1

j,分数阶微积分

Gmmwald—Letnikov的定义与分数阶微积分Riemann—Liouville的定义是完全等价的。在满足相同条件下,分数阶微积分的c叩uto定义和分数阶微积分的Gmmwald—Letnikov定义等价。当分数阶微积分阶次”为正整数和负实数时,分数阶微积分的c印uto定义和分数阶微积分的Riemann—Liouville定义满足以下关系:

:D㈧=:乒糍孵+:D荆(4)

1.2分数阶偏微分方程

分数阶偏微分方程是由传统整数阶偏微分方程拓展演化而得到的,是将传统整数阶偏微分方程中对时间或者空间的偏导数项用分数阶偏导数来替代而得到的㈤J。分数阶偏微分方程的基本形式为

!丝!!型2:拯!:吐+塑!!也

f一

5、

矿t妒x

护y

’根据偏微分阶次d和卢的不同将分数阶偏微分方程分为三种基本形式㈤““:

a)时间分数阶偏微分方程(时间偏导数阶次d为非整数,空间偏导数阶次卢为整数)。特别地,当0<d<1,口=2时,将

热传导方程(扩散方程)推广为时间分数阶热传导方程(时间分数阶扩散方程);当1<“<2,口=2时,该方程可以认为是热传导方程和波动方程的非线性插值。

b)空间分数阶偏微分方程(空间偏导数阶次口为非整数,时间偏导数阶次仅为整数)。特别地,当d=1,l<届<2时,将热传导方程(扩散方程)推广为空间分数阶热传导方程(空间分数阶扩散方程)。

c)空间一时间分数阶偏微分方程(时间偏导数阶次d为

非整数,同时空间偏导数阶次卢为非整数)。特别地,当0<

a<1,l<卢<2时,将热传导方程(扩散方程)推广为时间一空间分数阶热传导方程(时间一空间分数阶扩散方程)。

图像增强是数字图像处理的基本内容之一,以梯度算子和L印lacian算子等为代表的整数阶微分方法在图像增强中的应

2分数阶微分与图像增强

・416・

计算机应用研究第29卷

用已经相当广泛。近年来,很多学者将分数阶微分理论应用到图像增强中并且取得了较好的效果。分数阶微分可以在提高信号高频成分的同时非线性地保留信号中的甚低频,因此将该理论应用到图像增强中,可以使图像边缘更加突出,同时保留了图像平滑区域的纹理信息。分数阶微分理论应用到数字图像中的关键是增加一个自由度,即阶次。(0<”<1),通过适当调节”的大小可以构造相应的掩模算子,从而得到较好的图像增强效果。

2.1

分数阶微分对信号的作用

已知任意平方可积的能量信号以菇)Er(R),由信号处理

的基本理论可知其Fourier变换为八甜)=f。八戈)e一”出。假设信号八搿)的乃阶导数为尸(石)(n∈z+),根据Fourier变换的性质可以得到

盯^

D’弧x)铮(口D“(∞)=(妇)4以m)=d“((cJ),(∞)(6)

假设信号八戈)的分数阶”阶导数为厂(戈)(”∈月+),同样根据分数阶Fourier变换的性质可以得到

n^

D杈x)甘(功9。(∞)=(奶)”,(∞)=d”(∞班∞)

(7)

其帆

。..

b&)_l。附(乞):罕。gn(。)

一m)=(洳)”=“”(o)e‘∥‘“’

(8)

图1为分数阶微分算子在阶次不同时的幅频特性。从图1可以清楚地得出:a)分数阶微分算子对信号都有加强的作用,并且加强的幅度随频率和微分阶次的增加而非线性地急剧增强;b)2阶次微分相比1阶次微分对信号低频成分的提升较弱,但是对信号高频成分的提升却更明显;c)特别地,当微分阶次”=0时,表示分数阶微分算子不对信号进行任何改变;d)对于分数阶

微分阶次0<”<1时,在信号的低频部分(∞<1),分数阶微分算

子对信号的幅值进行了一定的提升,并且提升的幅度稍大于1阶次和2阶次微分;在信号的高频部分(∞>1),分数阶微分算子对信号的幅值进行一定的提升,但是提升的幅度也明显小于1阶次和2阶次微分。上述性质表明,分数阶微分算子具有弱导数性质,即分数阶微分算子在加强信号中高频成分的同时.对

信号的甚低频分量进行了非线性保留。此外,从信号处理的角

度分析分数阶微积分算子的物理意义可知,分数阶微分算子可以理解为广义的调幅调相,其振幅随频率呈分数阶幂指数变化,相位是频率的广义Hilbert变化。

一I三潍;

S2

{==茹5a

.∥/

//

墨・

一匕≤芝≥0

,一一一j曩互国龟罗舒;t帮j=

.盖矗笙孝三哆r

.q

频率w/Hz

图1分数阶微分算子的幅频特性曲线

2.2分数阶微分滤波器的构造

假设二维图像信号,(省,y)对于x和y轴方向的分数阶微分在一定的条件下是可分离的,这样可以利用分数阶微分算子分别对数字图像的八个方向进行滤波处理,然后将得到的中间结果通过相关规则结合。将图像信号,(立,y>)的持续期[o,£]

按单位间隔^=1等分,即n=I半l=[卜口】。这样可以得

r+

。1^=l

,‘

到{G—L定义下的分数阶微分算子沿x和y轴方向的数值计算

万方数据

表达式,如式(9)和(1U)所不,兵甲代表不分颧彤r微分效值计

:彤,(。,y),垒,(*,y)+(一口),(。一1,y)+L立毛型,(*一2,y)¨

算的余项,如式(11)和(12)所示。

i—二!尘jj二!二;;土j趔,(x一3,y)+R:(。,y)(9)

…。~7一"…

、‘7

:彤,(。,y),垒,(x,y)+(一口),(。,y一1)+L型妥竺业,(。,),一2)+.

R:(z)=ii—二!堕j。二!二;;_2j趟,(x,y一3)+R:(。,y)(10)

—二兰立j—=6

兰2生三土jji旦‘、’。¨7‘‘*、。’wjI兰幽,(x一4,,,)+、一。7

月;(z)=i二二旦2』乏‰m吨y)...+

,(m+1),(一口)…’…”7

—二!_=lL上{i二竺二二兰羔幽(z,,,..4)+…+

(11)、…

蒯‰职一-m)

,(m+1)厂(一")…。^~

(12)、…

图2表示分数阶微分图像增强掩模算子形,该算子将式(9)和(10)推广到图像的其余六个方向,从而得到一个基于八个方向的图像增强掩模算子(图中表格空白部分用0填充)。

该增强算子具有各向旋转不变性,增强算子形中的滤波系数

如式(13)所示。

%o=1%1=一口

%:巫掣

%:=型子坐

‘13’

%。=巫型必乏芋上业业

00%。O0

%。

00

00O阡0,0%,

O%,

000

阡名.

阡名.

阡名

阡名:

阡名.降名。阡0.%

%。

00

阡名.%.%,

阡名:

00

0OOO

%。

%。

图2分数阶微分增强掩模算子

3分数阶积分与图像去噪

众所周知,分数阶微分和分数阶积分是互为逆运算,同样图像增强和图像去噪也是相反的过程。既然分数阶微分

理论在图像增强中取得了很好的效果,因此,本章尝试将分

数阶积分理论用到图像去噪,提出了基于分数阶积分的图像去噪算法。该方法打破了传统图像去噪算法基于整数积分阶次的思想,通过设定一个较小的分数积分阶次来构建相应的图像去噪掩模,并利用迭代的方法来控制图像去噪的效果,从而在图像去噪的同时,尽量保留了图像的边缘和纹理等细节信息。

3.1

分数阶积分对信号的作用

对已知任意平方可积的能量信号以£)∈r(R),由信号处

理的基本理论可知其Fourier变换蒯∞)=k,(£)e一…出。假

第2期

黄果,等:分数阶微积分在图像处理中的研究综述

・417・

设信号八f)的n阶导数为厂(t)(nEz),根据Fourier变换的性质可以得到

Fr^

D杈f)铮(口D“(∞)=(池)”,(ccJ)=d“(∞)八∞)

(14)

这样将信号八£)的分数阶”导数为,(f)(”∈R+),同样根

据Fourier变换的性质可以得到

盯^

D弧£)甘(口D”(m)=(妇)。八∞)=d”(∞小∞)

(15)

rd”(m)=(汹)”=n”(∞)e妒(“)

其中:

‰o川计删o):孚。gIl(。)

(16)

于是,当”∈R一时,根据分数阶算子理论,即D~=,,并且”’=一”,可以得到分数阶积分算子的Fourier变换形式为

盯^

,。八£)甘(矿)”‘(∞)=(洒)”“c【,)=∥(∞)・,(∞)

(17)

r,’(甜)=d”’(埘)e‘伊I“)

其吼

k(:川计。∥(2):半。gn(。)

(18)

根据式(18)的描述,可以得出分数阶积分算子在阶次。’0,0.2,0.5,0.8,1.0,1.5,2.0时的幅频特性,如图3所示。

j3.

l\

I—v=Q.Q

;l{}S2.

¨’;

+v=1.Ov=0.5

—*一v=0.8~

胃趔\N酞

+v=1.5一罂1.

\C再N

一v=2.0

\ ̄ ̄l ̄

~~一鼍杂心

O.

1嫌遵基釜斟≮娶堇j

——寻==j!

频率w/Hz

图3分数阶积分的幅频特性曲线

从图3可以清楚地得出:a)分数阶积分对信号都有衰减作用,并且随着频率和积分次数的增加而非线性地急速衰减;b)2阶次积分相比1阶次积分对信号低频成分的提升较大,同时对信号的高频成分的削弱也明显;c)特别地,当”=0时,表示分数阶积分算子不对信号进行任何改变;d)对于积分次数0<w<1时,在信号的频率∞<l的区域,分数阶积分算子对信号的幅值进行了一定的提升,但是提升的幅度远小于1次和2次积分;在信号的频率∞>1的区域,分数阶积分算子对信号的幅值进行了一定的削弱,但是削弱的幅度也明显小于1次和2次积分。上述性质表明,分数阶积分算子在削弱信号高频部分的同时,对其中的最高频部分进行了非线性的保留;分数阶积分算子在加强信号低频部分的同时,对其中的最低频部分也进行了一定的保留。这样,分数阶积分算子在去除噪声的同时,一定程度上保留了图像的边缘和纹理等细节信息,从而使去噪后的图像不会产生严重的模糊现象。3.2分数阶积分滤波器的构造

仍然假设二维图像信号,(髫,y)对于夏轴和y轴方向的分数阶积分在一定的条件下是可分离的,这样可以利用分数阶积分算子分别对数字图像的八个方向进行滤波处理,然后将得到的中间结果通过相关规则结合。将图像信号,(石,y)的持续期

[a,£]按单位间隔危=l等分,即n=I半l=【£一o】。这样可

r.

。1^=l

,‘

以得到G—L定义下的分数阶积分算子沿x轴和y轴方向的数值计算表达式,如式(19)和(20)所示,其中尺’表示分数阶积分数值计算的余项,如式(21)和(22)所示。

:墨’,(x,y),垒,(¥,y)+F’,(x—l,),)+型掣,(x一2,),)+…+

…1、

万方数据

!:£!:!二攀,(*一3,y)+尺:(z,y)

(19)

:c

7,(*,y),垒,(。,y)+Ⅳ,,(Ⅸ,y一1)+丛妥山,(。,y一2)”.+

矗墨(z)=!:£!!±』jj』;:二必,(*一4,,,)+…+

!:£!:±』攀,(z,),一3)+R;(z,y)

(20)

R’;’(*)=竺』』£j,(m+1)厂(口’)…’…”7

高‰毗吨y)L!上工j:÷兰必,(z,,,一4)+一(21)、…

羔揣如扩m)

・+,(m+1)厂(F’)…。’…7

(22)、…

图4表示分数阶积分去噪掩模算子形,该算子将式(19)和(20)推广到图像的其余六个方向,从而得到一个基于八个方向的去噪掩模算子(图中表格空白部分用0填充)。该去噪算子具有各向旋转不变性,去噪算子形中的滤波系数如式

(23)所示。

‰=l

觋,=”

%:华

%:止岩

心3’

%。=亟必鼍≯掣业

00嘭

0O

O0孵’

形O

O0O

形形形0

彬彬形彬

OO

矿%

0O

0矿

0形

0O

0O

形O

0O%

图4分数阶积分去噪掩模算子

4分数阶偏微分方程与图像去噪

基于偏微分方程(panial

differential

equations,PDE)的图像

处理是一门崭新的交叉学科分支,融合了数学、物理学、电子信息学和计算机科学与技术等多个领域的相关理论知识。迄今为止,国内外许多学者已经建立了大量经典的基于PDE的图像处理模型,并显示出了巨大的生命力Ⅲ“J。基于整数阶PDE的图像去噪模型具有强大的数学理论基础,因此该方法具有其他去噪方法不具备的诸多优势,但是国内外许多学者在深入研究的过程中仍然发现PDE图像去噪模型存在很多不足之处。比如P—M”u和ROF模型¨2,5纠虽然能够有效地去除平滑区域的噪声,但是去噪后的图像会产生阶梯效应而出现分片光滑现象;LLT模型"4o利用了高阶PDE建模,虽然能够消除阶梯效应,但是去噪后的图像视觉效果不好;McM去噪模型p”“1虽然保护了图像的边缘信息,但是去噪后的图像会产生许多人为痕迹。为了解决各种整数阶PDE图像去噪模型存在的问题,有学者提出将分数阶微积分理论引入到PDE图像

去噪模型中,构建基于分数阶PDE的图像去噪模型。该模型

・418・

计算机应用研究第29卷

将分数阶梯度算子替代整数阶梯度算子来构建去噪模型的能量泛函。

4.1

基于分数阶Fourfer变换的图像去噪模型

因为分数阶有界变差空间曰矿(n)是传统的有界变差空

间曰y(n)的推广和延伸,By(n)空间是建立在传统梯度模值的基础上,相应地,曰矿(n)空间就应该建立在分数阶梯度模值基础上。这样可以建立基于分数阶有界变差空间日矿(门)的能量泛函来描述噪声图像。

已知噪声图像,(省,y),(戈,y)∈力,将ROF模型的正则项梯度算子V,用分数阶梯度算子V”,来替代,可以得到基于分数阶偏微分方程去噪模型的能量泛函:

E(,)=』n1D”,ld%dy+÷I|,0一,||乞

,∈BP(力),(Ⅳ,y)∈力

(24)

其中:fD”,f-√(彤,)2+(彤,)2、D=,和彤,分别表示图像,

(戈,y)关于x和y轴的分数阶导数。

为了得到上述能量泛函的EuIe卜Lagmnge方程,需要极小化上述能量泛函,因此,利用变分法的性质构造相应能量泛函。

g(8)=E(,+占妒)=

』nID。(,+印)Ickdy+{一||,0一,一即||乙

(25)

其中:占为任意实数,妒为任意函数。

利用泛函极值的基本性质,当P∈瞄(n),有g’。(0)=0,

可以得到相应的Eul矗一Lagrange方程:

』n(尚%+尚叫岫训舶棚删删(26)

为了得到最终的去噪模型,需要消除任意函数妒。因此,

采用ParseVal恒等公式』n,‘g出妙=去』n,g幽,幽z对式

(26)进行处理,可得

知(龟9尚+匆妒尚)圳旷

—]广一—广

^万万

去fn妒(,o一,)如1曲2=o

(27)

因为分数阶Fourier变换的性质有域妒:(口俐。)”9,联妒=

(丑叮『∞:)。妒,并消除系数l/(2订)可以得出

——广—r一

』n((细∥;尚郴酬;尚)岫慨一

^百百

A』n妒(毛一,)dmld∞2=O

(28)

又因为函数妒的任意性,利用变分法的基本引理可得出、

—广—广

。【(丑俐I)”面+(尼彻2)。面一A(,0—7)卸(29)。。(椭)”尚+(胁)。尚州名一?):o(29)

一一

对式(29)两边先去共轭,然后再利用Fourier逆变换可以得出相应的图像去噪的分数阶偏微分方程:

1,‘Re{[恧(尚)+谚(尚)】卜¨瑚㈣,

最后利用梯度下降法,引入时间变量f,可以得到基于

F0谢er变换的分数阶偏微分方程的去噪模型:

詈一e{[瓦(高)吲尚)肛c¨川,・,

该模型最大的难点就是如何求取分数阶微分算子的共轭

万方数据

算子D=和联。基于Fou矗er变换的分数阶偏微分方程的去噪模型是利用分数阶Fourier变换的性质来得到彤和联的数值解。

已知有噪声图像,(z,y)(省,y∈力),假设其在x和y轴方向相互独立,其石和y方向的分数阶Fourier变换为

知(m・^’训(1一≯M”沙㈨”聊)/Ⅳ

二,,(。,y)=,(。,,,)一,(。,y一1)=

‘32’

知m-‰(-一扩咧V州…蚓“

因此,由Fourier变换的性质可知,一阶差分的空域与频域的表达式对应关系为

『D;,(并,y)H(1一e一口删17Ⅳ),(∞l‘,∞2)

(33)

【D,,(并,y)+{(1一e—J2”“27Ⅳ),(∞l,m2)

同理可知,n阶差分的空域与频域的表达式对应关系为

『D:,(z,y)H(1一e—B”。17”)“,(∞l,∞2)

(34)

【D:,(z,),)H(1一e一声删27Ⅳ)“,(∞1,∞2)

由分数阶Fourier变换的基本性质可知,将整数阶次n推广到分数阶次”,即可得到基于分数阶Fourier变换的分数阶导数数值计算公式:

fD:,(z,),)H(1一e一口”’17”)”,(∞1,∞2)

(35)

【D:,(菇,y)十+(1一e一口伽’27Ⅳ)。,(∞1。,∞2)

这样,由式(35)可以得到分数阶微分算子的共轭算子:

{————.

『D:,(算,,,)++(1一e一正删17Ⅳ)。,(∞1‘,∞2)

(36)

【D:,(省,y)++(1一e一届俐27Ⅳ)”,(∞1,∞2)

已知有噪声图像,(聋,y),戈,y∈仃,并且假设噪声图像,(石,号,(戈,・)的分数阶微积分数值计算。由分数阶微积分的G—L

定义可知,图像信号,(互,・)的分数阶微分为

刚v)=(÷)“善(-1)’而器拦而如咄…s,)

当步长^=1时可以得到

:彤如)-i邑(。)“},(一,):剃加善㈢)U卜D

(38)

假设函数c”(茄)满足

c。(。)=f(一・)5(:)

算≥。

(39)

【o

膏<o

因此,可以将信号的分数阶微积分改写为信号,(石)与函

:叫,(g)=c”(z)+,(z)

(40)

于是,将式(40)代入基于分数阶曰旷(以)空间的能量泛函,可得到基于卷积积分的分数阶变分模型:

E(,)=』n

c“+,ld聋dy+—:一I|厂一,1l乙

(41)

c”},J=/(c”。,):+(c。},);2(c”+D;=』nc。(:),(*一z,y)出

(42)

为了得到式(41)的Euler_咖ge,与上节假设一样,仍然

【(c。¥,),=/nc。(z),(x,,一z)出

4.2基于分数阶基本定义的图像去噪模型

y)在x和y轴方向相互独立,因此只考虑x轴方向的图像信数c”(髫)的卷积操作形式:

其中:

第2期

黄果,等:分数阶微积分在图像处理中的研究综述

・419・

假设任蒽函数妒,从而构造函数,便兵满足式(43):

性具有特殊的生物视觉感受野模型,因此分数阶微积分理论g(s)=E(,+P妒)=』n

c”(,+印)Ickdy+争o,o一,一8妒|I乞(43)

为数字图像底层处理提供了一种崭新的研究方法。表1表示分数阶微积分理论应用于图像底层处理的类比,该表说明了仍然利用泛函极值的基本性质令g’(0)=0,可以得到

分数阶微积分理论应用于图像处理所涉及的定义、适用情。n

f。笠塑盘I磐坐兰丛尘韭出¨c”},I

““7

况、特点、优缺点情况。

Afn(,0一,)妒dxdy=O

(44)

表1

分数阶微积分应用于图像底层处理的类比

因为当(省,y)薯力时,,(菇,),)=0。所以将伫一R2时,可以得到Parseval恒等式:

^百

』R塘出dy2寿』R2,gd∞1d∞2

(45)

————F一———百一

可以根据Parseval恒等式将式(44)进一步改写为

去,n(箐等)c一:珐如,她+去,力(箐等)㈡:坫蛔幽:一

A寿』n(,0一,)妒d∞1d∞220

(46’

因为存在c”},=c。・,,并且因为函数妒的任意性,可以得到

(鬻)州‰(锗)以oAc名一2Ⅲ㈤,—r

本文介绍的各种基于分数阶微积分和分数阶偏微分方程因为c”(∞)=c”(一∞),可以得到

理论的图像处理模型还存在一些不足之处,在今后的研究工作(锗)以二训+(锗)以二㈦州J2一?)=o㈤,

中需要重点解决以下几个方面的问题:a)本文主要介绍分数阶微积分理论应用到图像去噪和增强两种基本的图像底层处利用F0urier逆变换的性质可以得到

((锗)甜㈠,)。+((锗)"㈠,),一

理,需要进一步将基于分数阶微积分的图像处理模型推广到图像分割、边缘提取、数字水印、指纹提取等图像处理的其他应用领域;b)本文介绍的基于分数阶微积分理论的图像处理模型A(如一,)=O

(49)

都只是应用于二维图像处理,因此,以后要逐步尝试将上述模慷沁)锗㈦埘胁协沁)锗㈠"胁

去掉卷积符号后可以得到

型推广到三维图像处理;c)分数阶微积分理论作为一种基本的数学分析工具,已经成功应用到了图像处理中,并取得了很A(,0一,)=0

(50)

好的应用效果,如何将分数阶微分理论和当今成熟的经典算法最后,利用梯度下降法即可以得到基于G-L定义分数阶结合,如利用神经网络、遗传算法、模糊数学等经典理论来估计偏微分方程的图像去噪模型:

分数阶微积分的阶次,从而构建性能更加完备的图像处理算詈=一(协沁)锗㈦刎胁

子;d)图像去噪和图像增强是互为逆过程,因此,可以尝试将基于分数阶偏微分方程的图像去噪模型利用反扩散原理推广』nc”(z)!i;÷等(x,,,+z)出)+A(,0一,)

得到基于分数阶偏微分方程的图像增强模型,并将该模型和其

(51)

他的一些经典图像增强方法结合得到性能更加优越的图像增强模型;e)本文介绍关于分数阶偏微分方程的图像处理模型5结束语

都是基于较弱函数空间条件下建立的,在数学完备性方面存在自从Mandelbmt发现分形现象,并将Riemann-LiouviⅡe分不足,因此需要寻找更为适合的基于分数阶微积分理论的图像数阶微积分理论应用于分析和研究分形媒介中的布朗运动描述空间对图像建模。

以来,分数阶微积分逐渐受到工程技术人员的高度关注。近分数阶微积分理论作为一种崭新的工程应用工具,已经引些年以来,分数阶微积分理论已经在化学、电磁学、控制学、起了国内外众多科技工作者的广泛关注,并在包括图像处理的材料科学和力学等众多学科领域中得到广泛应用并取得了实际应用领域中取得了较好的成果。因此可以大胆预测,随着重大的研究成果。然而,将分数阶微积分理论应用到现代信分数阶微积分理论的进一步发展和完善,数值计算方法和滤波号分析与处理中,特别是数字图像的处理,在国内外都还是器构造方法的进一步简化和改进,学术前沿的新方法和新技术一个新兴的学科分支。分数阶微积分理论适合研究具有非的进一步推广和应用,计算机硬件设备运算能力的进一步发展线性、非因果、非平稳等特征不确定的信号,因此分数阶微积

和提高,分数阶微积分理论未来必将具有更加广泛的应用分理论中许多经典算法都非常适合在现代信号分析与处理

空间。中应用。众所周知,数字图像中邻域内像素点之间的灰度值参考文献:

具有高度的自相似性,这种高度的自相似性可以近似为分形[1]

LEIBNIzG

w.corre叩ondaticedeLeibniz

avec

Hugens,vaJlzuli.

结构,并通常是以复杂的边缘和纹理等细节信息来表示。此ch帆et

leMarquisdeL’H0spital[M].P耐s:Libr.de

A.F咖ck,

外,数字图像的分数阶微分具有特殊的马赫现象,其拮抗特

1853:297—302.

万方数据

・420・计算机应用研究

第29卷

『2]

MANDELBROTBB,vanNESsJw.ThefhctionalbrDunianmo.

tions,fraction8lnoisesand

applications[J].SIAM

ReView,1968,

10(4):422—437.

[3]MANDELBROTBB.7IhefmctaJ学eometlyofnature[M].New

York:

W.H.Freem粕andCompany,1983.[4]袁晓,张红雨,虞厥邦.分数导数与数字微分器设计[J].电子

学报,2004,32(10):1658—1665.

[5]阿然,齐林,王越.分数阶Fourier变换的原理与应用[M].北京:

清华大学出版社,2004:23—49.

[6]陶然,张峰,王越.分数阶Fourier变换离散化的研究进展[J].中

国科学E辑:信息科学,2008,38(4):481.503.

[7]薛定宇,赵春娜.分数阶系统的分数阶PID控制器设计[J].控

制理论与应用,2007,24(5):771—776.[8]

PODLUBNYI.Fractional.ordersvstemsaJld

controllers[J].1EEE

T陋ns

onAutomatic

ControI,1999,44(1):208—214.

[9]zAYEDAI.OntherelationshipbetweentheFburiertransfo唧and

fhctionalFourier

t啪sfom[J].1EEE

SignaIProcessLett,1996,3

(12):310—311.

[10]廖科.分数阶微积分运算数字滤波器设计与电路实现及其应用

[D].成都:四川大学,2006.[11]RODlEcK

w.Quantitativeanalysisof

catretinal

98n对ion

cell

re—

spon8e

to

vi8ual

stimuli[J].VisionResearch,1965,5(11):583—

601.

[12]蒲亦非,袁晓,廖科,等.一种实现任意分数阶神经型脉冲振荡器

的格形模拟分抗电路[J].四川大学学报:工程科学版,2006,38

(1):128一132.

[13]陈庆利,蒲亦非,黄果,等.分数阶神经型脉冲振荡器[J].四川

大学学报:工程科学版,2011,43(1):123—128.

[14]蒲亦非,王卫星.数字图像的分数阶微分掩模及其数值运算规则

[J].自动化学报,2007,33(11):1128-1135.

[15]PUYi—fei,WANGwei.】【ing,zHOuJi-liu,封口z.Fmctional

differen—

tial印proach

to

detectingtextur■featuresofdigitalimageanditsfhc—

tionaldiⅡerentialfilterimplementation[J].

Scjence-nChinaSe—

riesF:Inforfllation

Sciences,2008,51(9):1319一1339.

[16]蒲亦非.将分数阶微分演算引入数字图像处理[J].四川大学学

报:工程科学版,2007,39(3):124—132.

[17]蒲亦非.分数阶微积分在现代信号分析与处理中的应用的研究

ID].成都:四川大学,2006.

[18]、蒲亦非,王卫星,周激流,等.数字图像纹理细节的分数阶微分

检测及其分数阶微分滤波器实现[J].中国科学E辑:信息科学,

、2008,38(12).:2252—2272.

[19],蒲亦非.将分数阶微分演算引入数字图像处理[J].四川大学学

报:工程科学版,2007,39(3):124—132.[20]Pu

Yi—fei,zH0u

Ji—liu,yuAN

xiao.FractionaldiffeIentialmask:a

fractionaI

difIbrential-based如proach

fbrmultiscale

textu南enhance-ment[J].IEEE

Trans

on

Image

Proeessing,2010,19(2):491一

.一{

511.

[21]黄果,蒲亦非,陈庆利,等.非整数步长的分数阶微分滤波器在图

像增强中的应用[J]..四川大学学报:工程科学版,2011,43(1):

129一136.

[龙q,黄果,蒲亦非,陈庆利,等.纂于分数阶积分的图像去噪[".+系统

工程与电子技术,2011,33(4):925—932.

[23]杨柱中,周激流,黄梅,等.基于分数阶微分的边缘检测[J].四川

大学学报:工程科学版,2008,40(1):152-157.

[24]杨柱中,周教流,晏祥玉,等.基于分数阶微分的图像增强[J].计

算机辅助设计与图形学学报。2008,20(3):343—348.

万方数据

[25]杨柱中,周激流,黄梅,等.用分数阶微分提取图像边缘[J].计算

机工程与应用,2007,43(35):15-18.[26]MATHIEU

B,MELcHl0RP,0uTSALOuPA,doZ.Fmctionaldif-

ferentiationfbredgedetec60n[J].SignaIPorcessing,2003,83

(11):2421—2432.

[27]李远禄,于盛林.分数阶差分及边缘检测[J].光电工程,2006,33

(12):70一74.

[28]汪凯宇,肖亮,韦志辉.基于图像整体变分和分数阶奇异性提取

的图像恢复模型[J].南京理工大学学报:自然科学版,2003,27

(4):400-4J04.

[29]刘红毅,韦志辉.基于分数阶样条小波与IHs变换的图像融合

[J].南京理工大学学报:自然科学版,2006,30(1):82—84.[30]LIuJun,cHEN

song—can,TANxiao—yang.Fractionalordersingular

value

deco‘叩ositionrepresentationforfacereco印ition[J].

Pattem

Recognition,2007,41(1):168-182.

[31]左凯,孙同景,李振华,等.二维分数阶卡尔曼滤波及其在图像处

理中的应用[J].电子与信息学报,2010,32(12):3027—3031.[32]汪成亮,兰利彬,周尚波.自适应分数阶微分在图像纹理增强中的

应用[J].重庆大学学报,20ll,34(2):32—37.

[33]高朝邦,周激流.基于四元数分数阶方向微分的图像增强[J].自

动化学报,2011,37(2):150一159.[34]BAIJian,FENG

xiang—chu.Fractional-orderanisotmpicdi饪hsionfbr

image

denoising[J].IEEE

Trans

on

ImageProcessing,2007,16

(10):2492-2502.

[35]张军,韦志辉.sAR图像去噪的分数阶多尺度变分PDE模型及

自适应算法[J].电子与信息学报,2010,32(7):1654—1659.[36]张军,韦志辉.一种基于卷积积分的图像去噪变分方法[J].中国

图象图形学报,2008,13(9):1673二1677.

[37]张军,韦志辉.分数阶图像去噪变分模型及投影算法[J].计算

机工程与应用,2009,45(5):1—6.

[38]张军.基于分数阶变分的图像建模与去噪算法研究[D].南京:

南京理工大学,2009.[39]zHANG

Jun,w_EIzhi—hui.AclassoffbcⅡonal-order舢lti—scale

v耐ational

modelsand

altematingprojectionalgofithmforimagede—

noising[J].App¨ed

MathematicaI

ModeⅢng,20“,35(5):2516—

2528.

[40]沈淑君.分数阶微分方程的数值解方法及误差分析[D].厦门:厦

门大学,2005.

[41]庄平辉,刘发旺.空间一时间分数阶扩散方程的显式差分近似

[J].高等学校计算数学学报,2005,27(51):223-228.[42]MAINARDI

F,LucHKOY,PAGNINIG,PARADIsIG.Thefunda-

mentalsolutionofthespace—timefractionaldiffusion[J].

FracCalc

AppI

Anal,2001,4(2):153—192.

[43]HuANGFeng—hui,uu

Fa-wang.The

fund砌ental

solutionofthe

space—tirrlefhctionaladvection1dispersionequation[J].

JournaIof

App¨edMathematics&Compu“ng,2005,18(1—2):339-350.[44]鸭cHuMPE砒正D,DERICHE

R.Vector—valuedimageregIllariza一

“.tion:whhPDEs:。acDmmon。fI’舭单Iework

fbrdi]日ferent

applications[J].

。IEEETrans

On

Pattern

AnaIysisandMachineInte¨igence,

2005,27(4):506—517.[45]GILBOA

G,socHENN,zEEVIYY,Variationaldenoisingofpar【一

ly

texturedimagesbyspatiaIlyvaryi“g

constmints[J].

JEEETrans

on

Image

Processing,2006,15(8):228l一2289.

[46]朱立新,王平安,夏德深.引入耦合梯度保真项的非线性扩散图

像去噪方法[J].计算机研究与发展,2咖,44(8):1390-1398.

(下转第426页)

・426・

tionwith

计算机应用研究

mobile

第29卷

on

phone:e%ctsoftouchkeysizeandt叫chkeyloca.

pleasure

theIntemetshoppingexpe^ence[J].JournaI

ofRetai-

tion[J].IntematjonaI

(1):68-76.

JournaIofIndustriaIErgonomics,20lO,40Iing,2002,78(1):31—40.[30]FRANKA

M,cHRISTOPHM.Mobile

Webusabilityevaluation—com.

[19]CHANG

dren

on

on

C.u8abilitycomparisonofpen-basedi“putforyoungchil—mobile

biningthemodifiedthinkaloudmethodwiththetestingofemotional,

devices[c]//Proc

ofIEEEIntemationalCor迁erence

co邸itiVeandconatiVeaspects0ftheusageof

WebapPlication

on

SensorNetworks,Ubiquitous,andTmstworthyComputing.2008:[C]//Pmc

ofthe9t}I

IEEE/AClsIntemationalCoIlference

Com.

531.536.

puterandInfoHnation

s.Mobile

Science.Washin群onDC:IEEEComputerSo—

[20]sHREsTHAwebbrowsing:usabiIitystudy[c]//Procof

on

ciety,2010:367・372.

the4thIntemationalConferenceandSystemsandthe1

st

MobileTechnology,Applications,

on

[31]BANDuRA

nitive

A.Socialfoundationsofthoughtandaction:asocialco哥

clif!fj:Prentice—HallInc,1986.

IntemationalSymposiumComputerHumantheory[M].En毋ewood

IntemctioninMobileTechnoIogy.NewYork:ACMPress,2007:187一194.

[32]s‘I’APLEsD,HuLLANDJ,HIGGINSC.A

self-emcacytheoIy

ex.

planationforthemanagement0fremoteworke玛invinual

H,KIM

o唧i豫tions

[21]KIMJ,LEEY.Anempiricalstudyof

on

u8econtextsinthe

[J].

JoumaIof

Computer-MediatedCommunlcation,1998,3

mobileIntemet,focusingtheusabilityofinf0Ⅱnation

architecture

(4):758—776.

[J].Jnformation[22]Lu

Systems

Front_erS,2005,7(2):175—186.

on

[33]HsU

M,cHIu

c.Intemetself-emcacyaJldelectronicserviceaccep.

H,suP.Factorsafbctingpurchaseintentionmobileshopping

tance[J].DecisionSupportSystems,2004,38(3):369.381.[34]LuARN

tion

to

websites[J].IntemetResearch,2009,19(4):442—458.

P,uNH.T0ward

mobile

an

understanding0fthebehavioralinten—Computersin

HumanBehavior,

[23]HALLR,HANNAP.Theimpact0f

on

webpagetext-backgroundc010r

aestheticsandbehavioraIin—

use

baIIking[J].

combinations

readability,Tetention,2005,21(6):873—891.

tention[J].BehaVior&InformationTechnoIogy,2004,23(3):

183—195.

[35]TslA0usIsAs,GIAGLIsGM.Anempiricalassessmentofenvimn—

theusabilityof

on

mentalfactorsthatinnuence

sK,HAsSANMA.AsmaItframeworkfor

mobile

website[c]//

Wash.

[24]GuIRGuIS

and

web

content

Proc0fthe9山IntemationalConferenceMobileBusiness.

resources

ad印tationinmobiledevices[c]//Procofthe12thIn・

0n

in昏onDC:IEEECo。叩uter

Society,2010:161—167.cmwding’sinnuence

on

temationalConferenceAdv蛐cedCommunication

Technology.

[36]sE0s,JuNGs.

Per℃eivedcon8umer’s

20lO:487—492.

shopping

behavior[J].KoreanJoumalofDist—buaonResearch,

[25]KAIKKONENA.usabilityproblemsintoday’smobileofthe2nd

Intemetpor-

Mobile7rech—

2005(1):23-32.

tals[c]//Proc

Intemationalcorlference

on

[37]scHMIDT

than

A,BEIGLM,GELLERsENH.7Ihereismore

tocontext

nology,ApplicationsaIldSystems.2005:1—7.location[J].Computer&Graphics,1999,23(6):893.901,

D,PARKs,以叫.AframeworkforevallIatingthe

[26]AcTON

T,GOLDEN

w,GuDEAs,越of.usab订ityaIldacceptance

ofthe9th

[38]HE0J,HAMB

insmall—scTeen

inf0啪ationsystems[c]//PmcEur叩eaIlusabilityofmobileusability263.275.

phones

based

on

multi—level,hieraI℃hicalHlodelofwith

CollaborativeElectronicrence.2004.

CommerceTechnologyand

Research

Co出-

factors[J].InteractingComputerS,2009,21(4):

[27]NAH

F.A

to

study

on

tolerablewaitingtime:howIong

are

web

users

[39]LESEMANN

luatemobile

on

E,NATALIE

w,SVENK.Combiningmethods

toeva—

willi“g

wait[J].BehaVior&|nformationTechnoIogy,2004,23

usability[c]//Pmc

ofthe9thIntemationalconference

M0bile

DevicesaJld

Services.

(3):153一163.HumanComputerIntemctionwith

[28]TARAsEwIcH

P.

MobilecommeI℃eoppor【unitiesandchallenges:

2007:444—447.

designingmobilecommercethe

applications[J].

Communicationsof

[40]DuHH,TANG,cHENV.usabilityevaluationformobiledevice:a

0f

laboratory

aJld

field

ACM,2003,46(12):1・60.

s,KAHNB.cross・categorydfects

ofinducedamusaland

comp面son

tests[J].

JoumaIof

UsabⅢty

[29]MENONStudjes,2005,1(I):4—16.

268.

(上接第420页)

[47]suNxiao—li,wANG

Jun—ping,SONGGuo—xiang.Anisotropicdifh-

[53]0sHERS,RuDINLI,FA7rEMIE.Nonlineartotalvariationbased

sionfilteringmethodwithand

fidelity

tem[J].SystemsEngineering

noiseremoval

algorithms[J].Physica

D,NonIinearPhenomena,

EIectronics,2008,30(8):1421—1423.1992,60(3):259-268.

[48]陈利霞,冯象初,王卫卫,等.加权变分的图像去噪算法[J].系

统工程与电子技术,2010,32(2):392—395.[49]u£E

SH,sE0

[54]LYsAKERM,LUNDERVOLDA,TAI

xc.

Noiseremovalusing

to

fourth—orderpanialdifferentialequationwith印plication8

medical

on

JK.NoiseremovalwithGausscurvature.drivendifhJ—

Trans

on

ma驴eticreson锄ceimagesinspaceandtime[J].ImageProcessing,2003,12(12):1579-1590.

1EEETrans

sion[J].IEEE

909.

Image

Processing,2005,14(7):904—

[55]KAssM,wITKIN

models[J].

A,TERz0POuL0sD.snakes:active

contour

[50]王新楼,乔明,邹谋炎.一种基于偏微分方程的sAR图像去噪方

法[J].电子与信息学报,2005,27(9):1365—1369.[51]PERONA

P,MALIK

InternationaIJournaIof

ComputerVision,1988,1

(4):321—331.

J.scale—spaceandedgedetectionusinganiso—

IEEE

Trans

on

[56]xuchen—yallg,JERRY

[J].1EEE[57]xu

Trans

on

LP.snal【es

sh印esandgmI|ientvect叫flow

tI、叩icdi】6fusion[J].

chine

PatternAnaIysisandMa-

lmage

Processing,1998,7(3):359—369.

vector

Inte¨igence,1990,12(7):629.639.

L,0sHERs,FA7IEMIE.。Nonlineartotalvariationbased

chen—ya“g,JERRYLP.Generalized目mdient

nowexte卜

[52]RuDING

nalforcesforactive131一】39.

contours[J].SignalProcessing,1998,71(2):

noisel-emoval

algorithms[J].PhysicaID,1992,60(1—4):259—

万方数据

分数阶微积分在图像处理中的研究综述

作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:

黄果, 许黎, 蒲亦非, HUANG Guo, XU Li, PU Yi-fei

黄果,HUANG Guo(乐山师范学院智能信息处理及应用实验室,四川乐山,614000), 许黎,XU Li(乐山师范学院物电学院,四川乐山,614000), 蒲亦非,PU Yi-fei(四川大学计算机学院,成都,610064)计算机应用研究

Application Research of Computers2012,29(2)4次

参考文献(57条)

1. LEIBNIZ G W Correspondatice de Leibniz avec Hugens,van Zulichem et le Marquis de L' Hospital 18532. MANDELBROT B B;Van NESS J W The fractional brounian motions,fractional noises and applications 1968(04)3. MANDELBROT B B The fractal geometry of nature 1983

4. 袁晓;张红雨;虞厥邦 分数导数与数字微分器设计[期刊论文]-电子学报 2004(10)5. 陶然;齐林;王越 分数阶Fourier变换的原理与应用 2004

6. 陶然;张峰;王越 分数阶Fourier变换离散化的研究进展[期刊论文]-中国科学E辑 2008(04)7. 薛定宇;赵春娜 分数阶系统的分数阶PID控制器设计[期刊论文]-控制理论与应用 2007(05)8. PODLUBNY I Fractional-order systems and controllers[外文期刊] 1999(01)

9. ZAYED A I On the relationship between the Fourier transform and fractional Fourier transform 1996(12)10. 廖科 分数阶微积分运算数字滤波器设计与电路实现及其应用[学位论文] 2006

11. RODIECK R W Quantitative analysis of cat retinal ganglion cell response to visual stimuli 1965(11)12. 蒲亦非;袁晓;廖科 一种实现任意分数阶神经型脉冲振荡器的格形模拟分抗电路[期刊论文]-四川大学学报(工程科学版)2006(01)

13. 陈庆利;蒲亦非;黄果 分数阶神经型脉冲振荡器[期刊论文]-四川大学学报(工程科学版) 2011(01)14. 蒲亦非;王卫星 数字图像的分数阶微分掩模及其数值运算规则[期刊论文]-自动化学报 2007(11)

15. PU Yi-fei;WANG Wei-xing;ZHOU Ji-liu Fractional differential approach to detecting textural features of digitalimage and its fractional differential filter implementation[期刊论文]-Science in China Series F:InformationSciences 2008(09)

16. 蒲亦非 将分数阶微分演算引入数字图像处理[期刊论文]-四川大学学报(工程科学版) 2007(03)17. 蒲亦非 分数阶微积分在现代信号分析与处理中的应用的研究[学位论文] 2006

18. 蒲亦非;王卫星;周激流 数字图像纹理细节的分数阶微分检测及其分数阶微分滤波器实现[期刊论文]-中国科学E辑 2008(12)19. 蒲亦非 将分数阶微分演算引入数字图像处理[期刊论文]-四川大学学报(工程科学版) 2007(03)

20. PU Yi-fei;ZHOU Ji-liu;YUAN Xiao Fractional differential mask:a fractional differential-based approach formultiscale texture enhancement 2010(02)

21. 黄果;蒲亦非;陈庆利 非整数步长的分数阶微分滤波器在图像增强中的应用[期刊论文]-四川大学学报(工程科学版) 2011(01)22. 黄果;蒲亦非;陈庆利 基于分数阶积分的图像去噪[期刊论文]-系统工程与电子技术 2011(04)23. 杨柱中;周激流;黄梅 基于分数阶微分的边缘检测[期刊论文]-四川大学学报(工程科学版) 2008(01)24. 杨柱中;周激流;晏祥玉 基于分数阶微分的图像增强[期刊论文]-计算机辅助设计与图形学学报 2008(03)25. 杨柱中;周激流;黄梅 用分数阶微分提取图像边缘[期刊论文]-计算机工程与应用 2007(35)

26. MATHIEU B;MELCHIOR P;OUTSALOUP A Fractional differentiation for edge detection[外文期刊] 2003(11)27. 李远禄;于盛林 分数阶差分及边缘检测[期刊论文]-光电工程 2006(12)

28. 汪凯宇;肖亮;韦志辉 基于图像整体变分和分数阶奇异性提取的图像恢复模型[期刊论文]-南京理工大学学报(自然科学版)2003(04)

29. 刘红毅;韦志辉 基于分数阶样条小波与IHS变换的图像融合[期刊论文]-南京理工大学学报(自然科学版) 2006(01)30. LIU Jun;CHEN Song-can;TAN Xiao-yang Fractional order singular value decomposition representation for facerecognition 2007(01)

31. 左凯;孙同景;李振华 二维分数阶卡尔曼滤波及其在图像处理中的应用[期刊论文]-电子与信息学报 2010(12)32. 汪成亮;兰利彬;周尚波 自适应分数阶微分在图像纹理增强中的应用[期刊论文]-重庆大学学报 2011(02)33. 高朝邦;周激流 基于四元数分数阶方向微分的图像增强[期刊论文]-自动化学报 2011(02)

34. BAI Jian;FENG Xiang-chu Fractional-order anisotropic diffusion for image denoising[外文期刊] 2007(10)35. 张军;韦志辉 SAR图像去噪的分数阶多尺度变分PDE模型及自适应算法[期刊论文]-电子与信息学报 2010(07)36. 张军;韦志辉 一种基于卷积积分的图像去噪变分方法[期刊论文]-中国图象图形学报 2008(09)37. 张军;韦志辉 分数阶图像去噪变分模型及投影算法[期刊论文]-计算机工程与应用 2009(05)38. 张军 基于分数阶变分的图像建模与去噪算法研究[学位论文] 2009

39. ZHANG Jun;WEI Zhi-hui A class of fractional-order multi-scale variational models and alternating projectionalgorithm for image denoising 2011(05)

40. 沈淑君 分数阶微分方程的数值解方法及误差分析 2005

41. 庄平辉;刘发旺 空间-时间分数阶扩散方程的显式差分近似 2005(51)

42. MAINARDI F;LUCHKO Y;PAGNINI G;PARADISI G The fundamental solution of the space-time fractional diffusion2001(02)

43. HUANG Feng-hui;LIU Fa-wang The fundamental solution of the space-time fractional advection dispersion equation 2005(1-2)

44. TSCHUMPERLE D;DERICHE R Vector-valued image regularization with PDEs:a common framework for differentapplications [外文期刊] 2005(04)

45. GILBOA G;SOCHEN N;ZEEVI Y V Variational denoising of partly textured images by spatially varying constraints[外文期刊] 2006(08)

46. 朱立新;王平安;夏德深 引入耦合梯度保真项的非线性扩散图像去噪方法[期刊论文]-计算机研究与发展 2007(08)

47. SUN Xiao-li;WANG Jun-ping;SONG Guo-xiang Anisotropic diffusion filtering method with a fidelity term[期刊论文]

-Journal of Systems Engineering and Electronics 2008(08)

48. 陈利霞;冯象初;王卫卫 加权变分的图像去噪算法[期刊论文]-系统工程与电子技术 2010(02)

49. LEE S H;SEO J K Noise removal with Gauss curvature-driven diffusion[外文期刊] 2005(07)50. 王新楼;乔明;邹谋炎 一种基于偏微分方程的SAR图像去噪方法[期刊论文]-电子与信息学报 2005(09)

51. PERONA P;MALIK J Scale-space and edge detection using anisotropic diffusion[外文期刊] 1990(07)52. RUDING L;OSHER S;FATEMI E Nonlinear total variation based noise removal algorithms 1992(1-4)53. OSHER S;RUDIN L I;FATEMI E Nonlinear total variation based noise removal algorithms 1992(03)

54. LYSAKER M;LUNDERVOLD A;TAI X C Noise removal using fourth-order partial differential equation withapplications to medical magnetic resonance images in space and time[外文期刊] 2003(12)55. KASS M;WITKIN A;TERZOPOULOS D Snakes:active contour models 1988(04)56. XU Chen-yang;JERRY L P Snakes shapes and gradient vector flow 1998(03)

57. XU Chen-yang;JERRY L P Generalized gradient vector flow external forces for active contours[外文期刊] 1998(02)

引证文献(4条)

1. 勾荣 基于Riemann-Liouville分数阶微分的图像增强[期刊论文]-制造业自动化

2013(16)2. 勾荣 基于 G-L 分数阶微分的图像边缘检测[期刊论文]-计算机与现代化

2013(11)3. 勾荣 基于分数阶微分的图像增强算法[期刊论文]-电子科技 2013(12)

4. 魏文力. 李宝树. 崔克彬. 徐雪涛 基于改进的最大熵航拍输电线路图像复原方法[期刊论文]-科学技术与工程 2013(32)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_jsjyyyj201202003.aspx


相关文章

  • 分数阶控制系统
  • 2008年5月第15卷增刊 文章编号:167127848(2008) S 020052203 控制工程C ontrol Engineering of China May 2008V ol . 15, S 0 分数阶控制系统 吴 娟 (重庆邮 ...查看


  • _浙江省教师招聘考试小学数学考试大纲 1
  • 浙江省小学数学教师招聘考试说明 一.考试性质 浙江省教师招聘考试是为全省教育行政部门招聘教师而进行的选拔性考试, 其目的是为教育 行政部门录用教师提供智育方面的参考.各地根据考生的考试成绩,结合面试情况,按已确定的招聘计划,从教师应有的素质 ...查看


  • 浙江省教师招聘考试小学数学考试大纲
  • 浙江省小学数学教师招聘考试说明 一.考试性质 浙江省教师招聘考试是为全省教育行政部门招聘教师而进行的选拔性考试, 其目的是 为教育行政部门录用教师提供智育方面的参考.各地根据考生的考试成绩,结合面试情况, 按已确定的招聘计划,从教师应有的素 ...查看


  • 尔雅公开课数学文化答案
  • 数学文化(一) 12002年,为中国少年数学论坛活动题词"数学好玩"的是(). A.邓东皋B.钱学森C.齐民友D.陈省身 正确答案: D 2"数学文化"一词最早进入官方文件,是出现在中华人民共和国教育 ...查看


  • 中学数学思想方法论-在线作业
  • 1.第4题 著名的英国数学家罗素是( )学派的代表人物. A. 直觉主义学派 B. 形式主义学派 C. 逻辑主义学派 D. 以上答案都不对 您的答案:C 题目分数:1 此题得分:1.0 2.第5题 韦达定理和距离公式,解析几何中有关方程等命 ...查看


  • 商务智能技术发展和应用研究综述
  • 1 商务智能技术发展和应用研究综述 商务智能技术发展和应用研究综述 x 斌 1 1. 北京航空航天大学 www.wenshan.me 摘 要:首先介绍了商务智能技术的发展历史,数据仓库.联机分析处理.数据挖掘概念.其次介绍了商务智能在当前大 ...查看


  • 例谈概率论与微积分的联系及相互间的应用
  • 第4卷第3期2008年7月 沈阳工程学院学报(自然科学版) JournalofShenyangInstituteofEngineering(NaturalScience) Vd.4No.3Jul.2008 例谈概率论与微积分的联系及相互间的 ...查看


  • 从记数法到复数域:数系理论的历史
  • 从记数法到复数域:数系理论的历史 纪志刚,上海交通大学 一. 记数法.位置制和零 人类在进化的蒙昧时期,就具有了一种"识数"的才能,心理学家称这种才能为"数觉"(perception of numbe ...查看


  • 分数阶Lü系统中的混沌及其控制
  • 计算机科学2007V01.34N0.12 分数阶Lii系统中的混沌及其控制*) 武相军1王兴元2 (河南大学计算中心 开封475004)1 (大连理工大学电信学院计算机系大连116024)2 摘要研究了分数阶Ln系统的混沌动力学行为.数值模 ...查看


热门内容