本发明专利技术涉及一种基于拉普拉斯图特征向量的稀疏编码的图像去噪方法，利用能够表示图像整体结构信息的拉普拉斯图的特征向量，通过将其参与到对图像块稀疏编码的过程，构建了一种新的基于拉普拉斯图特征向量的稀疏编码模型。对于该种新稀疏编码模型的求解与优化，则采用多轮迭代优化的算法，并在每次局部优化时，将这种新稀疏编码模型转化为已知的双稀疏模型进行快速求解。与此同时，本发明专利技术还利用一种较优策略来选择图特征向量的数量，即通过控制失控率的取值来自适应地选择不同图像在不同噪声方差下合适的图特征向量数量，以提高去噪性能。实验表明，通过本发明专利技术可以有效的进行图像去噪。
【技术实现步骤摘要】
本专利技术涉及，属于图像 处理领域，具体涉及到基于拉普拉斯图特征向量的稀疏编码技术。
技术介绍
图像去噪技术是指含噪图像经过降噪处理得到高质量图像，图像去噪在图像处理 中起着重要作用。随着稀疏表示技术的发展，人们提出了许多基于稀疏表示的图像去噪算 法，与传统图像去噪算法相比，基于稀疏表示的算法图像去噪，能够有效表示图像的结构， 充分降低噪声。在基于拉普拉斯图特征向量的稀疏编码的图像去噪算法中，充分利用由图 像块关系构建的拉普拉斯图的特征向量表示图像结构信息，进行降噪处理。通过拉普拉斯 图特征向量既可以有效表示图像的整体结构，也可作为一组基函数来重构图像。故本专利技术 充分利用拉普拉斯图特征向量的特性，通过将其应用于稀疏表示，进一步增强图像去噪性 能。
技术实现思路
本专利技术涉及。 ，包括输入含噪图像采 样，目标函数构造及稀疏优化和去噪图像输出过程， 包含以下步骤： (1)、对输入含噪图像进行采样，获得N个大小为M1XM2像素（长为M i个像素，宽 为M2个像素）的图像块，记录各采样图像块在含噪图像中的位置，对第i个M i X 112大小的图 像块像素矩阵B1，按列进行堆砌，构成一个（M1M2)Xl的列向量 71。同时将所有列向量yi，i =1，2, . . .，N，进行顺序排列，构成输入含噪图像块矩阵Y = Iiy1 y2 · . . yN]; (2)、根据所有列向量yi，i = 1，2，...，N，在M1M2维空间的几何结构中，构建所有 列向量Y1之间关系图，计算其对应的图的拉普拉斯矩阵L，进而计算矩阵L对应的特征向量 Iii-L' (3)、构造基于稀疏编码理论的去噪模型，并对参数X进行优化求解。该 去噪模型的目标函数为：' i = 1，2，...，Ν，对应的约束条件为：其中，D = 为字典，dk为字典中的第k 个原子，k = 1，2, K.，，K为字典中的原子总个数，X = 为稀疏系数矩阵，Xi为稀疏系数，畫为X对应的优化稀疏系数矩阵。I I. I I。为〇范数，I I. I I 2为2范数，ε = YM1M2O 2为残差，γ为一加权系数，σ 2为噪声方差。71为输入含噪图像块矩阵Y的第1列 的列向量。C1是乘法矩阵C = UUt的第1列的列向量，U = 为根据步骤2矩 阵L计算的对应特征向量中前M个特征向量； (4)、利用字典D和优化稀疏系数矩阵X，获得优化的去噪图像块矩阵 歹=，免为Y1对应的去噪列向量。其f的求解公式为f = (5)、去噪图像的输出，将步骤（4)中获得的去噪图像块矩阵f的列向量龙 转化为M1XM2大小的图像块像素矩阵爲，并取代步骤（1)中含噪图像相应位置的图 像块B1，同时记录含噪图像中每个像素点上需要覆盖的去噪图像块个数^ i，j为 像素位置坐标，及对应的估计图像块的像素值為;^k= 1，2，...。利用获得的 覆盖估计图像块个数W1Y对#；进行权值相加，得到降噪后图像对应位置的像素值為;α分别为含噪图像的像素值和对应的加权系 数。 上述步骤⑶所述的基于稀疏编码理论的去噪模型：目标函数为约束条件为，进行X的优化求解，其按如下步骤进行： (a)、在对该基于稀疏编码理论的去噪模型优化求解X时，采用多轮迭代优化的算 法，每次迭代时依次对一个稀疏系数Xi，i = 1，2, ...，N，进行局部优化。对稀疏系数乂;进 行优化时，首先固定其他稀疏系数{X]} # i取值，求解仅与X 1有关的对应稀疏编码理论的去 噪模型。该去噪模型的目标函数为I约束条件为其中， Cu为步骤（3)中所述的列向量c :的第i行元素 ，e := y「Σ# Axpu; (b)、对步骤（a)中的约束公式采用如下公式进行松弛转换： 。其中，残差矩阵E= ，II. I Ip^frobenius范数，为 C的第i个行向量； (c)、对步骤（a)的目标函数 和（b)中的松弛后的约束条件进行等价，采用如下公式：目标函数为约束条件 为其中，归一化参数 W fgi/lkIL，ε r = Νε -Tr(ETE) + (g,χ) TETEg' y Tr(.)为矩阵的迹。此时，目标函数可利用标准的基于误差约束 的稀疏编码模型求得矣，并求得本次局部优化后的稀疏系数為，.即. 上述步骤（3)所述的基于稀疏编码理论的去噪模型，在每次迭代时依次对一个稀 疏系数Xl，i = 1，2，...，N，进行局部优化时，可能遇到某些稀疏系数无解的情况，即在目标 函数为，约束条件为采无解。此时，则按如下步骤进 行对吳的求解： (a)此时将目标函数，约束条件·的去噪模型替 换为如下去噪模型进行优化：目标函数 ；，约束条件为I Ix1I I。彡丁， 其中τ稀疏度，并进行对X1的优化求解。 (b)对X1的优化求解时，将目标函数约束条件 Ix1I I。彡τ转化为目标函数约束条件I Ix1I I。彡τ。可利用标 准的基于稀疏度约束的稀疏编码模型求得氧，并求得本次局部优化后的稀疏系数％，即 上述步骤（3)所述的基于稀疏编码理论的去噪模型，对于乘法矩阵C = 1]1/中U = 的特征向量个数M的选取采用较优取值算法。在特征向量个数M取值时，首 先定义失控率函数？_= NS/NX100%，其中，Ns为在最后一轮迭代中需按基于稀疏编码理 论的构造的去噪模型优化求解算法求解稀疏系数的总个数。对特征向量个数M从0至N2进 行搜索，主要通过线性搜索的算法进行搜索。当P_< T p，Tp为一预设门限值时，则定义此 时的搜索值M。为较优的M取值。 本专利技术利用能够表示图像整体结构信息的拉普拉斯图的特征向量，通过将其参与 到对图像块稀疏编码的过程，构建了一种新的基于拉普拉斯图特征向量的稀疏编码模型。 对于该种新稀疏编码模型的求解与优化，则采用多轮迭代优化的算法，并在每次局部优化 时，将这种新稀疏编码模型转化为已知的双稀疏模型进行快速求解。与此同时，该算法还利 用一种较优策略来选择图特征向量的数量，即通过控制失控率的取值来自适应地选择不同 图像在不同噪声方差下合适的图特征向量数量，以提高去噪性能。通过该方法可以有效的 进行图像去噪。【附图说明】 图1是在噪声偏差〇 = 20时的降噪性能比较图像； 图2是在噪声偏差〇 = 40时的降噪性能比较图像； 图3是在噪声偏差〇 = 60时的降噪性能比较图像； 图4是在噪声偏差〇 = 20时的PSNR和失控率的变化图； 图5是在噪声偏差〇 = 40时的PSNR和失控率的变化图； 图6是在噪声偏差〇 = 60时的PSNR和失控率的变化图。【具体实施方式】 下面对本专利技术的一种基于拉普拉斯图特征向量的稀疏编码图像去噪算法作进一 步阐述。[002当前第1页1 2 3 本文档来自技高网...
【技术保护点】
一种基于拉普拉斯图特征向量的稀疏编码的图像去噪方法，其特征在于，其包含以下步骤：(1)、对输入含噪图像进行采样，获得N个大小为M1×M2像素的图像块，M1×M2像素表示长为M1个像素，宽为M2个像素，记录各采样图像块在含噪图像中的位置，对第i个M1×M2大小的图像块像素矩阵Bi，按列进行堆砌，构成一个(M1M2)×1的列向量yi，i＝1,2,...,N，同时将所有列向量yi，i＝1,2,...,N，进行顺序排列，构成输入含噪图像块矩阵Y＝[y1 y2 ... yN]；(2)、根据所有列向量yi，i＝1,2,...,N，在M1M2维空间的几何结构中，构建所有列向量yi之间关系图，计算其对应的图的拉普拉斯矩阵L，进而计算矩阵L对应的特征向量(3)、构造基于稀疏编码理论的去噪模型，并对参数X进行优化求解；该去噪模型的目标函数为：i＝1,2,...,N，对应的约束条件为：l＝1,2,...,N；其中，D＝[d1 d2 ... dK]为字典，dk为字典中的第k个原子，k＝1,2,...,K，K为字典中的原子总个数，X＝[x1 x2 ... xN]为稀疏系数矩阵，xi为稀疏系数，为X对应的优化稀疏系数矩阵，||·||0为0范数，||·||2为2范数，ε＝γM1M2σ2为残差，γ为一加权系数，σ2为噪声方差，yl为输入含噪图像块矩阵Y的第l列的列向量，cl是乘法矩阵C＝UUT的第l列的列向量，U＝[u1 u2 ... uM]为根据步骤2矩阵L计算的对应特征向量中前M个特征向量；(4)、利用字典D和优化稀疏系数矩阵获得优化的去噪图像块矩阵Y~=y~1y~2...y~N,]]>为yi对应的去噪优化的列向量，其的求解公式为Y~=DX~C;]]>(5)、去噪图像的输出，将步骤(4)中获得的去噪图像块矩阵的列向量转化为M1×M2大小的图像块像素矩阵并取代步骤(1)中含噪图像相应位置的图像块Bi，同时记录含噪图像中每个像素点上需要覆盖的去噪图像块个数wi,j，i,j为像素位置坐标，及对应的估计图像块的像素值k＝1,2,...,wi,j，利用获得的覆盖估计图像块个数wi,j，对进行权值相加，得到降噪后图像对应位置的像素值P^i,jk=(Σwi,jP^i,j+αPi,j)/(wi,j+α),]]>和α分别为含噪图像的像素值和对应的加权系数。...
【技术特征摘要】
【专利技术属性】
技术研发人员：汤一彬，张凤，张燕，高远，韩庆邦，徐宁，
申请(专利权)人：河海大学常州校区，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。