1 引 言

近年来,稀疏表示在图像去噪中的应用越来越广泛,成为前沿研究课题。Mallat等[1]提出超完备字典,基于超完备字典的稀疏表示是图像稀疏去噪的重要研究方向,而字典的设计则是该方向上重要的研究环节。全局字典以及自适应字典是字典的两类重要分支,最近几年,超完备字典通过学习和训练而获得的方法得到很大的进展。Aharon M,Elad M等[2]提出了K-SVD(K-singular value decomposition)算法,达到了良好的去噪效果。然而传统的K-SVD算法在一些方面仍然存在不足,具体表现如下:字典在更新过程中会出现每一列无法全部更新或部分列失效的问题,原因是初始化字典的选择不恰当,最终结果会导致信号利用率低;在噪声比较强的时候,字典中含有与噪声匹配的大量噪声原子,导致字典在学习过程中无法有效地将噪声去除;K-SVD算法运算时间较长,体现在图像训练过程中会占用大量的时间;传统最速下降算法锯齿现象严重。

根据服务区域内多年来降雨强度、一次降雨历时的发生概率统计,以10年一遇大暴雨事件发生和区域面积因素,雨水调蓄池容积为[10]

湿法脱硫反应速度快，脱硫效率高，一般均高于90%，技术成熟。缺点是设备腐蚀性严重，能耗高，占地面积大，投资和运行费用高。干法和半干法的设备简单，投资和运行费用较低、无污水处理系统等。缺点是反应速度慢，脱硫率低，磨损、结垢现象比较严重，废渣处理困难。

本文针对传统的K-SVD算法的几点不足,对其进行了改进:首先,通过稀疏贝叶斯学习(sparse Bayesian learning,SBL)[3]对图像进行预处理,提高图像信号的利用率;然后,对字典进行优化,并结合Bartlett检验法[4]将噪声原子裁剪,提高去噪后图像质量;第三,针对运算时间较长的问题,研究凸优化算法中的最速下降法,将改进的最速下降算法与正交匹配追踪(orthogonal matching pursuit,OMP)相结合;第四,将最速下降算法进行一维搜索,有效地抑制了锯齿现象。实验结果表明,通过本文的图像去噪方法,图像精度和峰值信噪比(peak signal to noise ratio,PSNR)值都得到了提高。

2 图像的稀疏表示模型

2.1 K-均值聚类算法

通过K-均值(K-means)聚类算法将图像块进行聚类处理,将具有相似集合形态的图像块聚集在一起,再通过K-SVD算法对每一个相似块训练出各自的优化字典,最终得到的字典更具有针对性。

综上所述，将危机管理应用于控制呼吸内科住院患者的医院感染控制中，对于降低医院感染率，促进患者的恢复，节省住院费用有显著的效果。

K-均值聚类算法的主要方法步骤如下:首先从全体n个数据中选择k个聚类中心。根据剩下的n-k个对象与k个聚类中心的相似度,将每个对象分配到k个聚类中心中与其最为相似的聚类中去。如果类中心有变化则对其进行反复聚类计算,直到类中心不再发生变化以及平方误差函数收敛为止。平方误差函数的定义可以表示为[5,6]:

 

(1)

式中:xji为第j类第i个样本;kj为第j类的均值或称聚类中心;nj为第j类样本个数;E为数据对象及其相应的聚类中心的均方差总和。

随着我国城市化建设的迅速发展，人居环境出现较大变化，人们对原生态的旅游环境表现出较大的兴趣。保护乡村原生态旅游环境，能够为我国人民提供更加舒适和优美的旅游环境，进一步促进乡村经济的发展，使乡村旅游的市场竞争力得到有效提高。同时，对乡村原生态旅游环境进行保护，能够减少生态环境污染问题，改善当地的人居环境，使人们的生活质量得到有效提升。然而，乡村原生态旅游环境保护具有较强的复杂性及长期性，要想保证乡村原生态旅游环境保护的整体效果，使乡村原生态旅游环境保护的目标得以顺利实现，必须深入了解乡村生态环境污染现状，进而采取有效策略进行环境保护，为乡村旅游事业的可持续发展提供有利支持。

2.2 稀疏贝叶斯学习

在K-SVD算法中,字典在更新过程中经常会出现每一列不能全部更新或部分列失效的现象,为了解决这个问题,可以通过贝叶斯学习[7]对信号进行预处理,通过贝叶斯概率模型进行迭代直到信号的稀疏表示。

传统K-SVD的方法在训练含噪图像之后,得到的字典中会出现与噪声匹配的噪声原子,因此对于字典进行优化的关键步骤就是删除这些噪声原子。然而关键问题是:如何对噪声原子进行检测。文献[8]中提出利用假设检验来检测噪声原子。噪声原子采用对角、反对角、水平、垂直4个扫描方向进行对比扫描,更能提高相关性检测的能力。

针对热度图的特点，我们可以把学生生活数据信息中需要和空间相结合的数据信息在地图上附加绘制出热度图，借此反映出学生的爱好等各种数据信息。

βij协方差矩阵为

 

(2)

式中:βij为稀疏向量;Γ为k×k阶对角矩阵,γ1,γ2,…,γk为对角元素;D为过完备字典矩阵;σ为噪声标准差。

 

(3)

式中:Xij为图像子块。

2) 更新噪声:

 

首先，用字母表示数，并不是方程所独有的特征。在很多公式、规律或函数等的表达中，都会借助于字母来表示更一般的情形。比如，借助a+b=b+a来替代文字叙述——交换两个加数的位置，和不变；借助于c=4a来表示正方形的周长公式。其次，从数学发展史的角度来看，方程思想的出现早于用字母表示数。远在公元前3000多年的古埃及文化中，就出现了一元一次方程形式，但直到十六七世纪，随着笛卡尔和韦达的工作，用字母来表示数才得以规范化。

韩光曙以感谢信一事为例分析道，医务人员遵章守纪是底线，是必须做到的60分，从医院层面来讲，没有任何可责备之处。但如果在服务中能够赋予更多的关心关爱，践行医院文化和价值理念，便可以为医院加分，为患者带来更多获益和良好的就医感受。

(4)

式中:Rij为提取子图的矩阵;X为含噪图像块。

 

(5)

式中:为协方差矩阵第i对角元素。

2.3 字典噪声原子裁剪

1) 计算βij的协方差矩阵与βij:

设μi为特征向量的方差,在噪声原子4个方向上用Bartlett检验法测试特征向量的方差μi是否相等,或者至少有2个原子的方差不等,Bartlett检验统计量H的定义为:

 

(6)

式中:n为原子元素数量。

对原子库中的原子进行判断,若某个原子满足H<χ2(β;3),则认为该原子为噪声原子,并且进行删除。χ2(β;3)表示自由度为3的卡方分布上的β百分点位置。与通过挖掘原子库的冗余性来减小原子数目的传统字典优化方法对比可以发现,该方法不会对信息原子产生影响,只会从字典中删除噪声原子,从而提高了逼近质量以及去噪效果。

目前,稀疏分解方法可分为基追踪(BP),匹配追踪(MP)和正交匹配追踪(OMP)。OMP是最好的匹配追踪算法并且具有最适用的基础表达功能。式(7)可以转换为无约束最优化问题,即:

2.4 改进的最速下降正交匹配追踪

传统的最速下降算法虽然拥有很多的优点和可取之处,但是该算法依然存在缺点。在迭代过程中迭代点呈现了相应的“之”字型,此现象为锯齿现象,对收敛速度会产生影响,但是锯齿现象又是无法避免的。因此可以通过抑制锯齿现象从而可以提高收敛速度。

稀疏分解图像的目标是找到一个矩阵α并且满足:

 

(7)

式中:||α||0为α的非零元个数。

改善关系不是去讨好孩子、一味地认同他，而是要把孩子需要解决的问题和父母与孩子的亲密关系分开，比如学习和亲密关系是两回事，不要说“你不好好学习就不是我的儿子”。

 

(8)

式中:λ为正则化参数。

尽管OMP算法比其他方法快,但其仍以损失一定的精度为代价。考虑到凸优化算法具有较高的精度但收敛速度较慢的特点,本文将改进的最速下降算法与OMP算法相结合,即改进的最速下降OMP。利用改进的最速下降算法去求解式(8)。

人才是国家发展的重要资源，而教育是培养人才的主要渠道，小学教育作为基础阶段，应当受到国家和教育部门的高度重视，尤其是在信息化时代的大背景下，做好小学教育的信息化发展意义重大。

通常将式(8)看成是一种多元函数f(x)的无约束优化问题,其寻优表达式为[9,10]:

αk+1=αk+tkpk

(9)

式中:tk为优化迭代的步长;pk为第k次迭代搜索方向。

通常认为在一个点上最快的下降方向是这一点的负梯度方向。于是可以定义:

pk=-f(αk)

(10)

每次优化迭代的步长tk可以通过式(11)求得,即:

 

(11)

如果tk-1f(X(k+1))Tp(k-1)<0,则转入式(12),否则转入式(14)。

令p′k=X(k)-X′,从X(k)出发,沿着p′k搜索方向进行一维搜索,求出步长t′k,使得

 

(12)

若f(X(k)+t′kp′k)<f(X(k)+tkpk),置X(k+1)=X(k)+t′kpk,X′=X(k+1),k=k+1,转到式(10),否则,令X′=X(k)转到式(13)。

X(k+1)=X(k)+tkpk

(13)

置k=k+1,转到式(10)。

该算法通过在p=X(k+1)-X(k-1)方向上进行一维搜索,得到一个比原始算法更好的迭代点。

文献[11]对于式(8)的问题可以看作是求解迭代误差||X-Dα最小化的问题,那么问题可以等价成正定二次型函数的优化问题,如所示:

 

(14)

式中:G=DTD; b=DTX。

对式(14)求导可得:

f(Xk)=Gα-b=-DTr

(15)

pk=-f(Xk)=DTr

(16)

式中：r为残差矩阵。

 

(17)

2)更新I,I=I∪{i}

OMP的分解方法的目的是解决

 

(18)

步骤1:初始化。残差r0=X,αx是X(X是含噪图像块)的稀疏系数,Dx是X的字典,Dx∈D,ε是误差控制参数,原子下标索引集合I=φ。

步骤2:迭代L次,更新原子下标索引集I=I∪{i}。

1)选择原子i和相应的稀疏系数αx,通过集合I使得目标函数最小化

 

(19)

结合改进之后的最速下降算法在凸优化算法中的优点及贪婪迭代算法中OMP的优点。改进的最速下降OMP可以概括如下:

在本研究中，我们将PCK更多地看成是个人知识，强调教师的教学建构及在这一过程中的反思.因此，本研究采用质性研究的方法，采取个案研究策略，对具有一定代表性的某教师在新课程背景下的PCK发展现状进行深描，旨在深度地描写现状和全面地发现问题，而不是定量研究所追求的推广性.

 

(20)

4)更新残差向量

rk+1=rk-Dxtkpk

(21)

式中:为去噪后的图像;C为容差系数;σ为噪声标准差。

3 改进的K-SVD图像去噪算法

本文利用稀疏贝叶斯学习对含噪图像进行预处理,提高其信号利用率。并且结合改进的最速下降正交匹配追踪以及噪声原子裁剪对字典进行更新。具体算法步骤如下:

1.4.1 试验用药 止咳橘红颗粒，6 g/袋，生产批号151101；试验药模拟剂，6 g/袋，生产批号160101；金振口服液，10 mL/支，生产批号151017；对照药模拟剂，10 mL/支，生产批号160101；头孢呋辛酯干混悬剂，0.125 g/袋，生产批号150803。以上药物均由康芝药业股份有限公司提供。

步骤1:初始化阶。首先利用K-means聚类算法对大小为N的图像Y进行聚类,将图像分为n类(n≪N),并令有k 个超参数的向量为γ=[γ1,γ2,…,γk]T,得到图像子块Xij=RijX,i=1,2,…,k;j=1,2,…,n;Rij为图像子块矩阵。最后,初始化字典D(0)∈Rn×k为过完备DCT基,并令J=1。

步骤2:稀疏贝叶斯学习。当J≤Jmax时,针对每一类图像块Xij用稀疏贝叶斯学习不断更新噪声大小,并且求得稀疏系数向量αij。

步骤3:稀疏编码。当J>Jmax时,固定字典D,对式(22)利用改进之后的最速下降OMP算法继续求解稀疏向量

∀

(22)

步骤3:更新迭代残差,当||rk+1||2<ε停止迭代,其它情况则返回步骤2。

3)更新矩阵的最优稀疏系数αx

步骤4:字典更新。固定稀疏矩阵和D,依次对中的每个原子djl (l=1,2,…,k)进行更新。

1)定义稀疏系数不为零的图像块集合为ωl={(i,j)|αij(l)≠0}。

2)对每个稀疏系数不为零的图像块(i,j)∈ωl,计算残差得到残差矩阵El;

3)根据El=UΔVT,对其使用SVD进行分解,把U第一列替代字典的djl,用Δ(1,1)×V替换反复迭代更新,直到收敛稳定时结束字典学习;

步骤5:字典优化。进行噪声监测并删除噪声原子。

步骤6:图像去噪。在去噪过程中存在着图像块的重叠现象,因此应该对图像块的重叠部分做平均,去噪之后的图像的近似解为:

 

选取PSNR较大的图像作为最终去噪输出图像。

4 实验结果分析

实验是将本文算法与经典图像去噪算法:小波阈值去噪,基于非自适应冗余DCT字典去噪,基于传统K-SVD方法去噪在MATLAB上进行仿真对比。实验结果见图1、图2和图3。

  

图1 Boat原图与加噪图像

  

图2 不同字典对比

从仿真结果可以看出,本文方法对于含噪图像去噪效果明显,能够很好地保留图像的边缘纹理信息。与传统去噪方法进行对比之后,可以发现,本文算法的去噪效果优于传统去噪,图像精度提高。

选取2016-07—2017-07来我院就诊的需口腔正畸患者156例，按随机数字表法分为2组(n=78)，其中对照组男41例，女37例；年龄16～29岁，平均(22.64±5.24)岁；双牙弓前突45例，牙列拥挤33例。观察组男39例，女39例；年龄17～28岁，平均(22.65±5.41)岁；双牙弓前突40例，牙列拥挤38例。两组一般资料比较均无统计学差异(P>0.05)。本研究经本院伦理委员会审核批准[批准文号：(2016)第14号]。

对Barbara,Flatland和Boat 3种图像,选用不同去噪方法的峰值信噪比值比较结果见表1。

各地通常将智慧城市建设等同于建网络、装探头，或者热衷于建云计算中心、指挥中心。对如何推进符合本地特殊的应用软件开发则重视不够、投入不足，存在拿来主义的思想，将多地的先进模块引入、组合、拼装就变成自己的智慧城市建设内容。

  

图3 不同方法的去噪对比

 

表1 不同方法的峰值信噪比值比较 dB

  

图像去噪方法信噪比σ=25σ=35σ=45Barbara含噪图像20．161317．252415．0778小波23．528920．795318．9809基于DCT字典28．641626．790625．3234传统K⁃SVD29．562427．716826．1223本文方法30．470828．635927．2174Flatland含噪图像20．261517．286115．0929小波23．772621．023318．9987基于DCT字典31．021329．294827．948传统K⁃SVD32．109330．216628．787本文方法32．781831．645929．5438Boat含噪图像19．256217．654814．3698小波23．365420．321419．5987基于DCT字典31．125429．859227．6214传统K⁃SVD32．265430．525528．2136本文方法32．987531．698729．4512

通过表1可以发现,不论是稀疏性较好的Barbara图亦或是稀疏性较弱的Flatland和Boat图,利用本文方法,峰值信噪比都明显比其他方法提高不少,相对于传统的K-SVD,本文方法提高了大概1 dB。因此,本文方法相对于其他去噪方法,综合去噪能力显著提高。

5 结 论

首先通过稀疏贝叶斯学习对含噪图像进行预处理, 提高信号利用率。接着将改进之后的OMP算法,以及原子裁剪方法相结合,构成了改进之后的K-SVD算法对图像进行去噪。通过实验分析可以看出,相对于传统去噪算法,本文方法能够更好地提高图像的精度,图像边缘和纹理部分保留完好,峰值信噪比显著提高,并且去噪效率也有所提高。

[参考文献]

[1] Mallat S G, ZHANG Z F. Matching pursuit with time-frequency dictionary[J]. IEEE Transactions on Signal Processing, 1993, 41(12):3397-3415.

[2] Aharon M, Elad M, Bruckstrein A M. The K-SVD: an algorithm for designing of over-complete dictionaries for sparse representation[J]. IEEE Transaction on Image Processing, 2006, 54(11): 4311-4322.

[3] Wipf D P, Rao B D. Sparse Bayesian learning for basisselection[J]. Signal Processing IEEE Trans on, 2004, 52,(8):2153-2164.

[4] Bartlett M S. Properties of Sufficiency and Statistical Tests[J].Proceedings of the Royal Society of London. Series A,Mathematical andPhysical Sciences,1937, 160(901):268-282.

[5] 蒋帅. K-均值聚类算法研究[D]. 西安:陕西师范大学, 2010.

[6] 田源,王洪涛. 基于量子核聚类算法的图像边缘特征提取研究[J]. 计量学报，2016,37(6):582-586.

[7] ZHOU F, LI L. Modified sparse Bayesian dictionary learning K-SVD algorithm for video image sparse representation[J]. Journal of Computational Information Systems, 2015, 11(12):4567-4580.

[8] Lee J, Kim Y H, Nam J H. Adaptive noise reduction algorithms based on statistical hypotheses tests[J]. IEEE Transactions on Consumer Electronics, 2008, 54(3):1406-1414.

[9] 龙建辉, 胡锡炎, 张磊. 最速下降法解二次矩阵方程[J]. 湖南大学学报(自然科学版), 2008, 54(3):85-88.

[10] 郑慧峰,喻桑桑,王月兵,等. 基于经验模态分解和奇异值分解的振动声调制信号分析方法研究[J]. 计量学报，2016,37(4):398-401.

[11] 周灿梅. 基于压缩感知的信号重建算法研究[D].北京:北京交通大学, 2010.