法律状态公告日
法律状态信息
法律状态
2019-10-22
未缴年费专利权终止 IPC(主分类):G06T5/00 授权公告日:20171219 终止日期:20181031 申请日:20141031
专利权的终止
2017-12-19
授权
授权
2015-07-29
实质审查的生效 IPC(主分类):G06T5/00 申请日:20141031
实质审查的生效
2015-02-18
公开
公开
技术领域
本发明属于图像复原技术技术领域,涉及一种变分PDE图像保边保结构复原方法。
背景技术
图像在形成、记录、处理和传输过程中,由于成像系统、记录设备和传输介质等不完善,经常在图像中引入噪声、模糊或者部分信息丢失等质量下降问题。例如,大气湍流的扰动效应,环境条件的变化和传感元器件自身的质量原因产生的噪声干扰,被摄物与成像设备间的相对运动造成的运动模糊等。图像复原是利用图像的某些先验知识来重建图像从而改善图像质量的技术。然而,随着对图像质量和图像分辨率要求的提高,图像复原算法的代价和硬件实现的复杂度显著地增加,且实际应用中对图像复原算法通常有实时性要求,这就对现有数字图像复原技术提出了更大的挑战。
传统的图像处理技术主要分为基于概率统计的方法、基于小波理论的方法和基于变分偏微分方程的方法。最近二三十年来,基于变分偏微分方程的图像复原模型由于具有良好的数学理论和物理意义,同时可以开发高效的并行算法和软件实现,从而在图像复原领域内得到了广泛的关注,同时也预示了其广阔的应用前景。一个经典的变分偏微分方程方法是由Rudin、Osher和Fatemi提出的全变分模型(ROF模型):
>
其中K为模糊算子,f为退化图像,u为待复原图像。该模型由于含有全变分项,可以保持图像的边缘;同时是建立在能量泛函的基础上,尤其是最近几年可以利用分裂算法快速求解;因此在图像复原领域得到了广泛的应用和关注。然而ROF模型经常在图像的光滑渐变区域引入阶梯现象,可以利用高阶模型保持图像光滑渐变区域的特性来缓解阶梯现象:
>
其中
发明内容
本发明的目的在于克服上述技术存在的保持图像棱角和边缘及其光滑渐变区域等结构特征比较差、计算效率比较低和不利于处理大数据的的弱点,提供一种变分PDE图像保边保结构复原方法,该方法结合全变分型泛函空间和自适应的参数选取准则,采用将模型的目标函数转化为可高效求解多个子问题的策略,建立新的复原模型和开发出稳定高效的应用软件来复原高质量的图像,从而为遥感图像处理、医学图像处理、视频目标跟踪以及图像结构特征提取中图像分割、融合、补全、配准和识别等工作提供强有力的理论支持和技术支撑。
其具体技术方案为:
一种变分PDE图像保边保结构复原方法,包括以下步骤:
(1)输入一副大小为N×N的退化图像f(x,y),令f(i,j)表示图像在第i行第j列的像素值,利用Matlab中的diff函数定义函数g(x,y)的梯度算子▽g=(▽xg,▽yg)T和海森算子
▽2g=(▽xxg,▽xyg;▽yxg,▽yyg)T;
▽xg=[diff(g,1,2),g(:,1)-g(:,N)];
▽yg=[diff(g,1,1);g(1,:)-g(m,:)];
▽xxg=[▽xg(:,1)-▽xg(:,N);diff(▽xg,1,2)];
▽yyg=[▽yg(1,:)-▽yg(N,:);diff(▽yg,1,1)];
▽xyg=[diff(▽xg,1,1);▽xg(1,:)-▽xg(N,:)].
选取初始值u0=f,w0=0,v0=0,z0=0,α0=0,β0=0,λ0=0.
(2)求解发明的变分PDE模型模型
>
其中
>
>
由于上述问题是含有非光滑的l1范数的l1-l2优化问题,从而在数值实现中给直接求解带来了巨大的挑战。最近几年,利用算子分裂的方法在求解l1-l2问题领域得到了广泛的关注。其中一类有效的算子分裂方法是用增广拉格朗日方法,即引入约束条件
w=u,v=▽u,z=▽2u.
将其转化为下面极大极小的鞍点问题
>
其中
>
这里α,β,λ分别为对应的拉格朗日乘子,c1,c2,c3分别为对应正的罚参数,δC(w)为示性函数,即w∈C,则δC(w)=0;否则,δC(w)=+∞。
(3)极大极小的鞍点问题
●求解子问题uk+1。即考虑
>
由于wk,vk,zk,αk,βk,λk可看作常变量,因此舍弃与u无关的变量便有
>
上述优化问题为光滑优化问题,因此可求出其对应的显式解。即对应的Euler-Lagrangian方程为
(AIKTK-c1I-c2Δ+c3div2▽2)uk+1=AIKTf+αk+c1wk-divβk
-c2divvk+div2λk-c3div2zk.
其中div2为▽2的伴随矩阵。由于上述等式中选取的为循环边界差分格式,同时利用散度定理中的<p,▽q>=-<divp,q>和<s,▽2t>=<div2s,t>,因此可以利用快速FFT变化求解uk+1。即为:
>
其中F表示快速傅里叶变化,F-1表示快速傅里叶变化的逆变换。上述算子具体实现在程序中定义为:
●求解子问题wk+1。即考虑
>
舍去与w无关项,则对应的子问题为
>
该问题可以用投影梯度法快速求解,因此有
>
●求解子问题vk+1。即考虑
>
舍去与v无关项,则对应的子问题为
>
该问题对应经典的压缩阈值问题,因此其显式解为
>
●求解子问题zk+1。由于子问题z与v具有类似的形式,因此同样可写成其对应的显式解为
>
由于拉格朗日乘子对应对应的问题是平凡的,因此其更新变量即为
αk+1=αk+c1(wk+1-uk+1);
βk+1=βk+c2(vk+1-▽uk+1);
λk+1=λk+c3(zk+1-▽2uk+1).
●终止标准。由交替方向乘子法的收敛性定理可知,当>收敛到零时,算法的解收敛到初始问题的解。因此在算法中选取终止标准定为
与现有技术相比,本发明的有益效果为:
本发明中采有保持图像低阶全变分泛函空间耦合保持图像光滑渐变区域的高阶全变分泛函空间的方法建立新的变分PDE图像复原模型,从而建立的模型从理论上可以保持图像的棱角和光滑区域等结构特征;
本发明通过引入判断图像尺度自适应性的函数惩罚不同泛函空间的在模型中的权重,进而使得建立的模型能有效地刻画图像的局部信息,从而可以复原高质量的图像;
本发明中的软件包侧重于将模型转化为可快速求解的子问题,由于子问题具有可分性,且采用循环差分格式可利用快速傅里叶变换,因此可以高效求解且具有可并行性。
附图说明
图1是本发明基于耦合全变分泛函空间的图像保边保结构复原方法的流程示意图;
图2是图像去噪实验效果图,其中,图2(a)初始图像,图2(b)初始图像第255行像素值,图2(c)初始图像部分放大,图2(a1)噪声图像,图2(b1)噪声图像第255行像素值,图2(c1)噪声图像部分放大,图2(a2)ROF模型复原图像,图2(b2)复原图像第255行像素值,图2(c2)复原图像部分放大,图2(a3)高阶模型复原图像,图2(b3)复原图像第255行像素值,图2(c3)复原图像部分放大,图2(a4)建议模型复原图像,图2(b4)复原图像第255行像素值,图2(c4)复原图像部分放大;
图3是图像去模糊实验效果图,其中,图3(a)待修补图像,图3(b)ROF模型复原图像,图3(c)高阶模型复原图像,图3(d)建议模型复原图像;
图4是图像修补实验效果图,其中,图4(a)待修补图像,图4(b)ROF模型复原图像,图4(c)高阶模型复原图像,图4(d)建议模型复原图像。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合附图和具体实施例,进一步阐述本发明。
数值实验:
下述数值实验是应用Matlab7.6软件进行仿真实验,实验环境为Acer3820GT的个人笔记本电脑,Intel Core i5processing520M,2.40GHz,4GB内存。为了评估复原图像的质量,我们采用信噪比SNR来刻画:
>
其中
实施例1图像去噪实验:
如图2所示,通过对比不难发现,发明的模型能明显的保持图像的结构特征,如图像的边缘区域和光滑渐变区域。其中ROF模型复原图像的信噪比为17.7677,高阶模型的信噪比为17.9998,发明的模型信噪比为18.3694。
实施例2图像去模糊实验:
如图3所示,发明的模型可以较好地复原退化图像,尤其在保持图像的细节特征和结构特征更为有效。ROF模型复原图像的信噪比为14.3778,高阶模型的信噪比为14.4778,发明的模型信噪比为14.5118。
实施例3图像去模糊实验:
如图4所示,通过局部放大图像可以看出,ROF模型中在辣椒茎部含有明显的字母信息,高阶模型导致复原图像的边界模糊,然而建议的模型能很好的克服这两类模型的缺陷。对应地,ROF模型复原图像的信噪比为15.4055,高阶模型的信噪比为16.4336,发明的模型信噪比为16.6034。
以上所述,仅为本发明最佳实施方式,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。
机译: 用于制备适合消费的饮料的保形垫,包括制垫机和这样的保形垫的饮料制备机的组装,一种用于制造保形垫的预成型垫的方法一种方法,一种制备饮料的方法,该饮料适合使用这样的保形垫或组合件以及包括多个保形垫的零件套件
机译: 用于制备适合消费的饮料的保形垫,包括制垫机和这样的保形垫的饮料制备机的组装,一种用于制造保形垫的预成型垫的方法一种方法,一种制备饮料的方法,该饮料适合使用这种保形垫或这样的组件以及包括多个这种保形垫的零件套件。
机译: 用于制备适合消费的饮料的保形垫,包括制垫机和这样的保形垫的饮料制备机的组装,一种用于制造保形垫的预成型垫的方法一种方法,一种制备饮料的方法,该饮料适合使用这种保形垫或这样的组件以及包括多个这种保形垫的零件套件。