基于分数阶偏微分方程的图像放大模型

高　冉， 顾　聪， 李胜宏

(1. 中原工学院 理学院， 河南 郑州 450007； 2. 浙江大学 数学系， 浙江 杭州 310027)



基于分数阶偏微分方程的图像放大模型

高冉1， 顾聪1， 李胜宏2

(1. 中原工学院 理学院， 河南 郑州 450007； 2. 浙江大学 数学系， 浙江 杭州 310027)

将分数阶微分理论引入图像放大模型中，利用全变分思想，提出了基于分数阶偏微分方程的图像放大模型.仿真实验结果表明：新模型能较好地保持图像边缘特征，以及更多的图像纹理信息，优于整数阶微分方程放大算法，是一种有效、可行的图像放大模型.

分数阶；偏微分方程；变分；图像放大

Journal of Zhejiang University(Science Edition), 2016,43(5):550-553

对一幅图像进行放大是由低分辨率获得高分辨率图像的一种图像处理技术.通常理解为灰度级插值，仅针对图像的像素点进行处理，导致放大后的图像出现锯齿状条纹和块状效应，对带噪图像的放大效果更差；而基于偏微分方程(PDE)的图像放大方法可以有效抑制块状效应和噪声的影响[5-6].目前基于PDE的图像放大算法分为两类：直接法和间接法[7].本文采用的是间接法，即以传统的插值结果作为放大图像的初始值，应用PDE对其进行后处理，在Rudin, Osher, Fatemi(ROF)[8]和Lysaker, Lundervold, Tai(LLT)[9]模型的基础上，提出了整数阶的图像放大模型[10]：

(1)

其中，χΩ1是Ω1上的特征函数，即对任意的(x,y)∈Ω1都属于原始图像的像素点.仿真实验表明，该算法可有效消除图像边缘和纹理区域因放大产生的锯齿和块状效应，但存在放大后图像边缘模糊、清晰度欠佳等视觉缺陷.

为了解决图像模糊和细节特征丢失的问题，本文提出了一种基于分数阶PDE[11-13]的图像放大算法.分数阶微分具有大幅提升图像高频成分、增强图像中频成分、非线性保留图像低频成分的特性，更多地保留了图像平滑区域中灰度变化不大的纹理细节信息，使得放大后的图像清晰，能够较好地保持原图像的边缘特征.

1　相关理论

1.1分数阶导数

分数阶导数是整数阶导数的推广和延伸，虽然分数阶导数已广泛应用于物理、化学、生物、流体等领域，且在图像处理领域也取得了初步研究成果，但其定义并不统一.主要有3种经典的定义，包括Riemann-Liouville(R-L)、Capotu和Grumwald-Letnikov(G-L).

本模型基于G-L分数阶导数定义[14]：

(2)

其中，p为分数阶导数的阶数，

(3)

1.2分数阶导数的差分格式

将上述定义推广到二维情形，可得u(x,y)的分数阶偏导数定义：

(4)

(5)

(6)

集中技训从大一下学期陆续开展，每学期不超过一次，具体时间安排结合各学院专业课程要求进行。实施中，每学期固定两个教学周停开其他课程单独进行集中技训教学。采用“带教带学”模式进行授课，教师“带教”，班级辅导员督促“带学”。通过近几年的实施，学校教务处及学院各专业有效的整合了教学资源，教学实践潜力得到挖掘，目前各专业都不同程度取得了成绩，运行状况良好。

(7)

2　基于分数阶的图像放大算法

2.1模型的建立

先利用线性插值处理图像，将图像放大到所需倍数，然后应用分数阶PDE技术对插值放大后的图像进行后处理，实现图像放大的目的.

首先给出本模型的能量泛函式：

χΩ1(u-u0)dxdy,

(8)

其中，泛函中第1项是平滑项；第2项中的χΩ1是Ω1上的特征函数，即对任意的(x,y)∈Ω1都属于原始图像的像素点.λ(>0)为常数，是通过反复实验、人为选定的参数，其作用是对第1和第2项进行平衡.

利用变分法推得的该泛函欧拉-拉格朗日方程如下：

(9)

利用梯度下降法得到相应的扩散方程：

(10)

最终，利用式(10)对线性插值得到的u0进行修正，实现对图像的放大.

2.2模型的算法

第1步假设初始图像矩阵为m×n的T0，首先将T0通过采样缩小k(k为正整数)倍得到T′，将T′赋给T，然后对T进行线性插值(本文采用3次立方插值),得到放大k倍后的矩阵u0；

第2步利用式(10)对u0进行修正，找到一个u，既保持图像原有的特性，又保持了图像灰度级之间的连续性.

3　数值离散及实验结果分析

对式(10)进行离散，时间步长为Δt，得到

Δtλ(un-u0)χΩ1(un-u0),

(11)

离散化求解中用到的参数，在实验中按如下方式选取：式(11)中λ=0.5，h=1，Δt=0.05，p=1.7，迭代40次.为了测试本文算法的有效性，把256×256的Lena图像先采样缩小2倍，然后再放大.

图1　图像放大效果Fig.1　Results of image zooming

由效果图图1可见，基于分数阶偏微分方程的图像放大算法，对图像进行放大时都有较好的效果，保留了图像的原有特征.

为了进一步说明基于分数阶偏微分程和基于4阶整数阶偏微分方程的图像放大算法的区别，分别对描绘原始图像、分数阶放大后图像及4阶方程放大后图像的150行所有列的灰度曲线分布进行比较，如图2所示.

图2　图像放大前后灰度值曲线比较Fig.2　Comparison of different zooming methods

从图2(b)和(c)可看出分别基于整数阶和分数阶的图像放大算法的差异:使用整数阶模型放大时，虽然获得了较好的整体效果，但由于过度平滑，使得边缘模糊，细节部分不够清晰；而分数阶偏微分方程插值模型较忠实地反映了图像的原始面貌，保留了边缘的锐度和纹理特性.

4　结　论

数值实验表明，本文提出的基于分数阶偏微分方程的图像放大算法，较好地保留了图像的边缘特征和细节信息，获得清晰的放大图像，并且运算时间短于整数阶偏微方程，是一种有效、可行的图像放大算法.但本文仅尝试使用了最简单的分数阶微分，如何确定更佳的微分阶数以便得到更好的放大效果有待进一步研究.

[1]付树军，阮秋琦，王文洽.偏微分方程(PDEs)模型在图像处理中的若干应用[J].计算机工程与应用，2005(2)：33-35.

FUShujun,YUANQiuqi,WANGWenqia.ApplicationsofPDEsmodelinimageanalysisandprocessing:Asurvey[J]. Computer Engineering and Applications,2005(2)：33-35.

[2]张亶，陈刚.基于偏微分方程的图像处理[M].北京：高等教育出版社，2004.

ZHANG Dan, CHEN Gang. Image Processing Based on PDEs[M]. Beijing：Higher Education Press，2004.

[3]CHEN Y M, VEMURI B C, WANG L. Image denoising and segmentation via nonlinear diffusion[J]. Computers and Mathematics with Applications, 2000,39(5/6):131-149.

[4]WEICKERT J. Coherence-enhancing diffusion filtering [J]. International Journal of Computer Vision,1999,31(2/3):111-127.

[5]阮秋琦，忤冀颖.数字图像处理中的偏微分方程方法[J].信号处理，2012，28(3)：301-314.

YUAN Qiuqi, WU Jiying. Partial differential equation method on digital image processing[J]. Signal Processing, 2012，28(3)：301-314.

[6]KIM H, CHA Y, KIM S. Curvature interpolation method for image zooming[J]. IEEE Transactions on Image Processing,2011,20(7):1895-1903.

[7]朱宁，吴静，王忠谦.图像放大的偏微分方程方法[J].计算机辅助设计与图形学学报,2005,17(19):1941-1945.

ZHU Ning, WU Jing, WANG Zhongqian. Image zooming based on partial differential equations[J]. Journal of Computer-Aided Design & Computer Graphics,2005,17(19):1941-1945.

[8]CRANDALL M G, ISHII H, LIONS P L. User’s guide to viscosity solutions of second order partial linear differential equations[J]. Bull Amer Math Soc,1992,27:1-67.

[9]RUDIN L, OSHER S. Feature-oriented image enhancement using shock filters[J]. SIAM J Numerical Analysis,1990,27:919-940.

[10]GAO Ran, SONG Jinping, TAI Xuecheng. Image zooming algorithm based on partial differential equations technique[J]. International Journal of Numerical Analysis and Modeling,2009(6):5-18.

[11]蒲亦非，王卫星.数字图像的分数阶微分掩模及其数值运算规则[J].自动化学报,2007,33(11):1128-1135.

PU Yifei, WANG Weixing. Fractional differential masks of digital image and their numerical implementation algorithms[J]. Acta Automatica Sinica,2007,33(11):1128-1135.

[12]杨柱中，周激流，黄梅，等.基于分数阶微分的边缘检测[J].四川大学学报：工程科学版，2008，40(1)：152-156.

YANG Zhuzhong, ZHOU Jiliu, HUANG Mei, et al. Edge detection based on fractional differential[J]. Journal of Sichuan University:Engineering Science Edition,2008,40(1):152-156.

[13]艾必刚，罗以宁，蒋涛，等.分数阶微分梯度算子在图像增强中的应用[J].四川大学学报：自然科学版，2009，46(2)：343-347.

AI Bigang, LUO Yining, JIANG Tao, et al. Applications of fractional-order differential gradient Operator in Image Enhancement[J]. Journal of Sichuan University: Natural Science Edition,2009,46(2):343-347.

[14]PODLUBNY I. Fractional Differential Equations[M]. New York: Academic Press,1998.

Image zooming model based on fractional-order partial differential equation.

GAO Ran1, GU Cong1, LI Shenghong2

(1.CollegeofScience,ZhongyuanUniversityofTechnology,Zhengzhou450007,China; 2.DepartmentofMathematics,ZhejiangUniversity,Hangzhou310027,China)

A new image zooming model based on the fractional-order partial differential equation is proposed, which adopts the idea of total variation. Simulation results show that the new model is capable of preserving the characteristics of image edge, and it can retain more texture details than the integer order partial differential equation model. The model is therefore effective and practical for image zooming.

fractional-order; partial differential equation(PDE); variation; image zooming

2015-07-24.

国家自然科学基金资助项目(11401604);河南省科技厅基础与前沿研究项目(142300410354,142300410355,152300410226,152300410227);河南省教育厅科学技术研究重点项目(15A110045)；河南省高等学校重点科研项目(17A110036).

高冉(1982-)，ORCID:http://orcid.org/0000-0002-3151-7040，女，硕士，讲师，主要从事图像处理与模式识别研究,E-mail:nygr@163.com.

10.3785/j.issn.1008-9497.2016.05.010

TP 391

A

1008-9497(2016)05-550-04