CN106447632B - 一种基于稀疏表示的raw图像去噪方法 - Google Patents
一种基于稀疏表示的raw图像去噪方法 Download PDFInfo
- Publication number
- CN106447632B CN106447632B CN201610846774.5A CN201610846774A CN106447632B CN 106447632 B CN106447632 B CN 106447632B CN 201610846774 A CN201610846774 A CN 201610846774A CN 106447632 B CN106447632 B CN 106447632B
- Authority
- CN
- China
- Prior art keywords
- dictionary
- training sample
- rarefaction representation
- matrix
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012549 training Methods 0.000 claims abstract description 61
- 239000011159 matrix material Substances 0.000 claims abstract description 47
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 16
- 230000001965 increasing effect Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims 1
- 238000012545 processing Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000010790 dilution Methods 0.000 description 3
- 239000012895 dilution Substances 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- AHLBNYSZXLDEJQ-FWEHEUNISA-N orlistat Chemical compound CCCCCCCCCCC[C@H](OC(=O)[C@H](CC(C)C)NC=O)C[C@@H]1OC(=O)[C@H]1CCCCCC AHLBNYSZXLDEJQ-FWEHEUNISA-N 0.000 description 2
- 108091092878 Microsatellite Proteins 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011210 chromatographic step Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于稀疏表示的RAW图像去噪方法,包括以下步骤:步骤一、将待去噪RAW图像分解成为固定大小相互重叠的RAW矩形图像块;步骤二、将步骤一中得到的RAW矩形图像块按照不同颜色通道,重新排列成G1RBG2颜色层;步骤三、将步骤二中得到的G1RBG2颜色层按照G1、R、B、G2的顺序依次拉伸后得到列向量形式的图像块向量,将每个图像块向量从左至右拼接成G1RBG2矩阵,G1RBG2矩阵的每一列为一个训练样本;步骤四、从步骤三中的G1RBG2矩阵中随机选择部分训练样本作为训练样本集T,采用K‑SVD方法对训练样本集T进行学习得出字典D;步骤五、采用改进后的OMP算法对步骤四中得出的所述字典D进行图像重建后,形成无噪RAW图像;本方法解决了RAW图像去噪过程中计算复杂的问题。
Description
【技术领域】
本发明属于图像处理技术领域,特别涉及一种基于稀疏表示的RAW图像去噪方法。
【背景技术】
微小卫星相对于大卫星,具有质量轻,成本低,设计与制造周期短等特点。对地观测作为一项空间应用的重要任务,有着广泛的应用,如植被监测和灾难监测。在该领域,微小卫星已经成功得到应用。图像处理***是对地观测任务的关键组成部分。
但是,由于复杂的太空环境,相机拍摄的图像受到了严重的噪声影响。例如光线强度较弱将会带来很大的噪声,相机通过增加输入信号的感光度(增加ISO)来放大感光元件的电压,会造成噪声增强。数字相机中的图像前处理技术在原本噪声源的基础上还会带来特定的方法噪声,例如白平衡、色彩校正、伽马校正和边缘增强会改变噪声源的性质或增加噪声的强度。如果直接压缩图像,传输到地面处理,那么方法噪声、量化噪声和传输过程的冲击噪声会使得噪声的模型更加复杂。
目前传统的信号稀疏表示去噪方法对已知噪声模型(如加性高斯白噪声)的灰度或者彩色图像去噪。自然信号在冗余字典中稀疏分解,是一种有效的信号表示方法。在图像去噪领域,稀疏表示很好地将信号分解成稀疏向量,通过稀疏向量恢复无噪图像。
但实际的问题是,经过白平衡、去马赛克、颜色变换三种主要的图像前处理技术,数字相机输出YUV或者RGB图像噪声模型是未知的。而且图像具有信号相关、颜色相关和空间相关的噪声特性。如果简单地利用传统稀疏表示方法去噪,那么就会出现去噪效果差的问题。另一方面,原始图像(或称为:RAW图像)的噪声模型是已知的,即加性高斯白噪声。RAW图像的格式一般采用Bayer颜色滤波阵列,即图像中像素总数的一半是绿色(G),而总数的四分之一分别为红色(R)和蓝色(B),像素首先按GRGR排列,下一行像素按BGBG排列。如果直接处理RAW图像,相应的传统稀疏表示图像处理方法有限制,因为RAW图像数据不符合一般的图像假设,如:平滑性、分段不变性等。
【发明内容】
本发明的目的是提供一种基于稀疏表示的RAW图像去噪方法,以解决在处理未去噪RAW图像过程中,不能有效地估计出噪声的参数,计算复杂的问题。
本发明采用以下技术方案:一种基于稀疏表示的RAW图像去噪方法,其特征在于,包括以下步骤:
步骤一、将待去噪RAW图像分解成为固定大小相互重叠的RAW矩形图像块;
步骤二、将步骤一中得到的RAW矩形图像块按照不同颜色通道,重新排列成G1RBG2颜色层;
步骤三、将步骤二中得到的G1RBG2颜色层按照G1、R、B、G2的顺序依次拉伸后得到列向量形式的图像块向量,将每个图像块向量从左至右拼接成G1RBG2矩阵,G1RBG2矩阵的每一列为一个训练样本;
步骤四、从步骤三中的G1RBG2矩阵中随机选择部分训练样本作为训练样本集T,采用K-SVD方法对训练样本集T进行学习得出字典D;
步骤五、采用改进后的OMP算法对步骤四中得出的字典D进行图像重建后,形成无噪RAW图像。
进一步地,步骤四中采用K-SVD方法对训练样本集T进行学习得出字典D的具体步骤如下:
阶段一:初始化阶段,即从训练样本集T中随机选取i个n维样本ti用来组成字典D,并始字典D进行J’次迭代,J’为迭代总次数;
阶段二:固定字典D阶段,即采用OMP算法找出训练样本集T中的每个训练样本的稀疏表示αi;
阶段三:更新字典D阶段,即将固定字典D阶段得到的每个训练样本的稀疏表示组成稀疏矩阵A,通过稀疏矩阵A对阶段二中得出的字典D进行处理,依次更新字典D中的每一列,得出训练完毕的字典D。
进一步地,步骤五采用如下步骤具体实现:
5.1)设置各个训练样本的初始参数,并计算残差与训练完毕的字典D中每列dl的内积,l为字典D的列数,找出最大内积所对应的下标p;
5.2)更新索引集,更新残差,并将迭代变量的数值增加1;
5.3)当更新后的残差<重构误差阈值时,结束迭代,且当i个训练样本ti均都已找到稀疏表示αi时,采用加权求和将重叠的图像块组合,得出无噪RAW图像,当i个训练样本ti未全部找到稀疏表示αi,则返回步骤5.1),继续执行;
当更新后的残差≥重构误差阈值时,返回步骤5.2),继续执行。
进一步地,步骤5.1)采用如下步骤具体实现:
找到训练样本集T中每个训练样本ti的稀释表示,i=1,2,…,i,针对其中的一个样本,设置每个训练样本ti的初始化残差r(0),设置重构误差阈值为n(Cσ)2,设置索引集Ω0=0;
计算残差r(L-1)与训练完毕的字典D中每列dl的改进的内积,找出最大的内积所对应的下标p,p∈{1,2,…,k},即:
式中,dl为字典D的每一列,l=1,2,…,k;r(L-1)为OMP算法中迭代(L-1)次的残差;Bn为n×n的全1矩阵;γ用于调节颜色相关性;a为γ的分解因子;n表示图像块向量的大小,I为单位矩阵。
进一步地,步骤5.2)采用如下步骤具体实现:
更新索引集Ω(L)(L)=p,更新列dp构成集D(L)=D(:,Ω(L)(1:L)),D(L)为迭代L次时的字典D;
通过最小二乘法计算得到稀疏逼近βi (L)=argmin||ti-D(L)αi||2,
更新残差r(L)=ti-D(L)βi (L),更新残差之后,迭代变量L的数值加1。
进一步地,步骤5.3)采用如下步骤具体实现:
5.3.1)当r(L)<n(Cσ)2,得到样本ti的稀疏表示:αi=Sβi (L),结束迭代;并判断i个训练样本ti是否都找到稀疏表示αi;
当i个训练样本ti都已找到稀疏表示αi,则通过公式
采用加权求和将重叠的图像块组合,得出无噪RAW图像,式中,取λ=30/σ;
当未全部找到i个样本ti的稀疏表示,则返回步骤5.1),继续执行;
5.3.2)当r(L)≥n(Cσ)2时,则返回步骤5.2),继续执行。
本发明的有益效果是:通过对OMP算法进行改进,使其适应于RAW图像的重建,通过该方法能够很好地估计出来噪声的参数,而且图像去噪计算复杂度减少至少1/3,本方法主要是针对在太空中低照度下,低信噪比成像***去噪方法,不仅在太空中噪声较大的恶劣环境中具有较好的去噪性能,而且在成像过程中的RAW图像阶段去噪,能有效地在源头上抑制噪声,为后续图像处理和压缩提供必要的准备。
【附图说明】
图1为本发明一种基于稀疏表示的RAW图像去噪方法的流程图;
图2为本发明一种基于稀疏表示的RAW图像去噪方法中将图像块分离成G1RBG2颜色层步骤的示意图;
图3为本发明一种基于稀疏表示的RAW图像去噪方法将G1RBG2颜色层重组成G1RBG2矩阵步骤的示意图;
图4为本发明一种基于稀疏表示的RAW图像去噪方法中通过G1RBG2矩阵训练字典并重组图像的示意图。
图中,1.RAW矩形图像块,2.G1RBG2颜色层,3.G1RBG2矩阵,4.字典D,5.无噪RAW图像。
【具体实施方式】
下面结合附图和具体实施方式对本发明进行详细说明。
本发明公开了一种基于稀疏表示的RAW图像去噪方法,包括以下步骤:
步骤一、将待去噪RAW图像以给定步长分解成为固定大小相互重叠的RAW矩形图像块1;
步骤二、将步骤一中得到的RAW矩形图像块1按照不同颜色通道,重新排列成G1RBG2颜色层2;
步骤三、将步骤二中得到的G1RBG2颜色层2按照G1、R、B、G2的顺序依次拉伸后得到列向量形式的图像块向量,将每个图像块向量从左至右拼接成G1RBG2矩阵3,G1RBG2矩阵3的每一列为一个训练样本;
步骤四、从步骤三中的G1RBG2矩阵3中随机选择部分训练样本作为训练样本集T,采用K-SVD方法对训练样本集T进行学习得出字典D4。
步骤五、采用改进后的OMP算法对步骤四中得出的字典D4进行图像重建后,形成无噪RAW图像。
步骤一中的RAW矩形图像块1的大小为b*b,即RAW矩形图像块1中相邻两边的长度为均b个像素点,一般的b取值为8~16之间的偶数,优选的b取值为8;为了使图像块重叠,步长取小于b的偶数,优选的步长取值为2;且每个RAW图像块的左上角均是同一颜色分量。
步骤二中的G1RBG2颜色层2的颜色格式为:G1、R、B、G2;
步骤四中的采用K-SVD方法得出字典D4,其具体通过公式:
式中:n表示图像块向量的大小,n=4*4*4=64,k表示字典D4的大小,即列向量的个数,k=256,C表示噪声增益,取值为1.15,αij为稀疏表示向量;RijX为图像块向量;σ为噪声的标准差,D为字典D4,i为字典D4中的列数,j为字典D4的列数。
步骤四中,采用K-SVD方法对训练样本集T进行学习得到字典D4,目的是找出更好的字典稀释来表示样本数据。
此步骤分为三个阶段进行实施:
阶段一:初始化阶段,具体步骤如下:
从训练样本集T中随机选区i个n维样本ti,i=1,2,…,i;用ti来初始化初始字典D4的每一列dl,l表示字典D4的列数,l=1,2,…,l,即初始字典D4的每一列设置为i个训练样本的前l个;
初始字典D=[d1,d2,…dl]=[t1,t2,…tl];
然后,归一化初始字典D4,即将字典D4进行迭代,
设置迭代次数变量J,初始状态下J=1,每迭代一次后,次数变量J的值增加1,总共迭代J’次,J’为迭代总次数;
当迭代次数增加,去躁效果会提升,迭代次数越多,噪声标准差越小,而且根据零均值高斯噪声模型,采用经典PCA主成分分析方法估计标准差参数;
所以本方法设置:当初始标准差小于等于5时,J’=10;当标准差大于5时,J’=15;
阶段二:稀疏编码阶段,具体步骤如下:
当初始化阶段完成后,进入稀疏编码阶段,即固定字典D4,此步骤采用OMP算法来实现,目的是为了找到训练样本集T中每个训练样本的稀释表示。
4.1)依次针对训练样本集T的每一个样本:
初始化残差r(0)等于训练样本ti,设置稀疏表示迭代次数L=1,总共迭代L’=3次,即稀疏度为3,设置索引集Ω0=0;
其中,r(0)就是没有迭代的初始残差,r(1)为第一次迭代计算生成的更新生成的残差,直到更新完r(3),停止迭代;
4.2)计算前一次迭代中的残差r(L-1)与字典D4中每列dl的内积,通过r(L-1)和dl的内积值,找出最大的内积所对应的下标p,p∈{1,2,…,k},即:p=argmaxl=1,...,k|<r(L-1),dl>|;
将索引集更新为Ω(L)(L)=p,更新所选择的列dp构成集合D(L)=D(:,Ω(L)(1:L)),D(L)为迭代L次后的字典D4;
通过最小二乘法计算得到稀疏逼近αi (L)=argmin||ti-D(L)αi||2,αi为ti的稀疏表示;
更新本次迭代中生成的残差r(L)=ti-D(L)αi (L)之后,再次进行迭代,将迭代变量的次数更新为L+1次;
4.3)将L+1与L’进行比较;
4.3.1)如果L+1>L’,则结束迭代,并判断i个训练样本ti是否全部找到稀疏表示αi;
当i个训练样本ti全部找到稀疏表示αi,则进入字典更新阶段;
当i个训练样本ti未全部找到稀疏表示αi,返回步骤4.2);
4.3.2)如果L+1≤L’,继续返回步骤4.2);
阶段三:字典更新阶段,将上述稀疏编码阶段得到的稀疏表示αi,组成k行i列的稀疏矩阵A,其元素为αi,通过稀疏矩阵A处理字典D4,依次更新字典D4中的每一个原子列dl,l=1,2,…,l,具体步骤如下:
4.4)找到使用了当前需更新原子列dl的稀疏矩阵A的列元素集合{(l,i)|αli≠0},即对应于dl原子列的是稀疏矩阵A中的l行,找到l行中不为零元素所在的列序号i构成的集合(l,i)∈wl;
对于该集合wl中所有的元素i,选择稀疏矩阵A和训练样本集T的第i列分别组合成新的稀疏矩阵A’和样本矩阵T’,其中i为小于g的非0整数,再去除新的稀疏矩阵和样本矩阵中的第l行后,得出稀疏矩阵Al和样本矩阵Tl;
然后,将字典D4中去除第l行得到新的字典D4,通过公式El=Tl-DlAl得出误差矩阵El;
然后,应用奇异值分解SVD分解误差矩阵El=UΔVT,得出U矩阵和V矩阵,更新列dl为U矩阵的首列,更新系数值为V矩阵的首列,然后乘以Δ(1,1);
4.5)判断字典D4的k列原子是否依次更新完;
4.5.1)如果字典D4的k列原子没有更新完,则返回步骤4.4),继续执行;
4.5.2)如果字典D4的k列原子已更新完,判断迭代次数变量J的大小;
若J>J’,完成字典D4的更新,得出更新后的字典D4,并结束迭代;否则,返回步骤4.2),继续执行;
步骤五中,OMP算法是基于l0惩罚算法,实际就是寻找信号在字典D4的稀疏表示,即为下列优化问题:
该式的优化目标只有αij,因为字典D4在步骤四中已经训练完成。
通过优化完αij之后,再通过公式得出无噪图像。
在步骤五中,由于经典的OMP算法导致RAW图像空间的稀疏趋于R=B=G,使得彩色图像趋于灰度图像,无法保证重建图像块保持与原图像块具有相同的平均色彩,因此,采用对OMP算法进行改进,使其适应于RAW图像的重建。
通过改进后的OMP算法,计算所有n维图像块的样本集Tn×m的稀疏表示,固定步骤四中训练好的字典D4,其具体步骤如下:
5.1)找到每个训练样本ti的稀释表示,i=1,2,…,i,针对其中的一个样本,首先,初始化残差r(0)等于训练样本ti,设置重构误差阈值为n(Cσ)2,设置索引集Ω0=0;
然后,计算残差r(L-1)与字典D4中每列dl的改进的内积,找出最大的内积所对应的下标p,p∈{1,2,…,k},即:
式中,dl为字典D4的每一列,l为字典D4的列数,l=1,2,…,l;r(L-1)为OMP算法中迭代(L-1)次的残差,Bn为n×n的全1矩阵,γ用于调节颜色相关性;a为γ的分解因子,n表示图像块向量的大小,I为单位矩阵。。
即每个图像块和字典D4中的原子列都乘以尺度矩阵(I+(a/n)K),然后在计算它们的内积。
值得注意的是,由于字典D4乘以尺度矩阵(I+(a/n)K)后,各个列向量需要重新归一化,即用(I+(a/n)K)DS处理字典D4,得到归一化之后的字典D4,式中,S代表对角矩阵,矩阵S的对角线元素对应为(I+(a/n)K)D各列2范数的转置:
式中,diag表示把向量元素对角化,sum表示对矩阵各列求和,“./”和“.*”分别表示矩阵点除和点乘;
5.2)更新索引集Ω(L)(L)=p,更新所选择的列dp构成集D(L)=D(:,Ω(L)(1:L));
通过最小二乘法计算得到稀疏逼近βi (L)=argmin||ti-D(L)αi||2;更新残差r(L)=ti-D(L)βi (L),更新残差之后,迭代变量L的数值加1;
5.3)如果r(L)<n(Cσ)2,得到样本ti的稀疏表示:αi=Sβi (L),结束迭代;并判断i个训练样本ti是否都找到稀疏表示αi;
5.3.1)如果i个训练样本ti都找到稀疏表示αi,则通过公式
采用加权求和将重叠的图像块组合,估计去噪后的RAW图像;式中,λ表示拉格朗日乘子,它与噪声水平有关,噪声越大λ越小,因为噪声小的图像对重建图像影响大,噪声大的图像对重建图像影响应该小,取λ=30/σ。
5.3.2)如果未全部找到i个样本ti的稀疏表示,则返回步骤5.1),继续执行;
5.4)若r(L)≥n(Cσ)2,则返回步骤5.2),继续执行。
Claims (6)
1.一种基于稀疏表示的RAW图像去噪方法,其特征在于,包括以下步骤:
步骤一、将待去噪RAW图像分解成为固定大小相互重叠的RAW矩形图像块(1);
步骤二、将所述步骤一中得到的RAW矩形图像块(1)按照不同颜色通道,重新排列成G1RBG2颜色层(2);
步骤三、将所述步骤二中得到的G1RBG2颜色层(2)按照G1、R、B、G2的顺序依次拉伸后得到列向量形式的图像块向量,将每个所述图像块向量从左至右拼接成G1RBG2矩阵(3),所述G1RBG2矩阵(3)的每一列为一个训练样本;
步骤四、从所述步骤三中的G1RBG2矩阵(3)中随机选择部分所述训练样本作为训练样本集T,采用K-SVD方法对所述训练样本集T进行学习得出字典D(4);
步骤五、采用改进后的OMP算法对所述步骤四中得出的所述字典D(4)进行图像重建后,形成无噪RAW图像(5)。
2.如权利要求1所述的一种基于稀疏表示的RAW图像去噪方法,其特征在于,所述步骤四中采用K-SVD方法对所述训练样本集T进行学习得出字典D(4)的具体步骤如下:
阶段一:初始化阶段,即从训练样本集T中随机选取i个n维样本ti用来组成字典D(4),并使字典D(4)进行J’次迭代,J’为迭代总次数;
阶段二:固定字典D(4)阶段,即采用OMP算法找出训练样本集T中的每个训练样本的稀疏表示αi;
阶段三:更新字典D(4)阶段,即将固定字典D(4)阶段得到的每个训练样本的稀疏表示组成稀疏矩阵A,通过稀疏矩阵A对所述阶段二中得出的字典D(4)进行处理,依次更新字典D(4)中的每一列,得出训练完毕的字典D(4)。
3.如权利要求1所述的一种基于稀疏表示的RAW图像去噪方法,其特征在于,所述步骤五采用如下步骤具体实现:
5.1)设置各个训练样本的初始参数,并计算残差与训练完毕的字典D(4)中每列dl的内积,l为字典D(4)的列数,找出最大内积所对应的下标p;
5.2)更新索引集,更新残差,并将迭代变量的数值增加1;
5.3)当更新后的残差<重构误差阈值时,结束迭代,且当i个训练样本ti均都已找到稀疏表示αi时,采用加权求和将重叠的图像块组合,得出无噪RAW图像(5),当i个训练样本ti未全部找到稀疏表示αi,则返回步骤5.1),继续执行;
当更新后的残差≥重构误差阈值时,返回步骤5.2),继续执行。
4.如权利要求3所述的一种基于稀疏表示的RAW图像去噪方法,其特征在于,所述步骤5.1)采用如下步骤具体实现:
找到所述训练样本集T中每个训练样本ti的稀疏表示,i=1,2,…,i,针对其中的一个样本,设置每个训练样本ti的初始化残差r(0),设置重构误差阈值为n(Cσ)2,C表示噪声增益,σ为噪声的标准差,设置索引集Ω0=0;
计算残差r(L-1)与训练完毕的字典D(4)中每列dl的改进的内积,找出最大的内积所对应的下标p,p∈{1,2,…,k},即:
式中,dl为字典D(4)的每一列,l=1,2,…,k;r(L-1)为OMP算法中迭代(L-1)次的残差;Bn为n×n的全1矩阵;γ用于调节颜色相关性;a为γ的分解因子;n表示图像块向量的大小,I为单位矩阵。
5.如权利要求3或4所述的一种基于稀疏表示的RAW图像去噪方法,其特征在于,所述步骤5.2)采用如下步骤具体实现:
更新索引集Ω(L)(L)=p,更新列dp构成集D(L)=D(:,Ω(L)(1:L)),D(L)为迭代L次时的字典D(4);
通过最小二乘法计算得到稀疏逼近βi (L)=argmin||ti-D(L)αi||2,
更新残差r(L)=ti-D(L)βi (L),更新残差之后,迭代变量L的数值加1。
6.如权利要求3或4所述的一种基于稀疏表示的RAW图像去噪方法,其特征在于,所述步骤5.3)采用如下步骤具体实现:
5.3.1)当r(L)<n(Cσ)2,C表示噪声增益,σ为噪声的标准差,得到样本ti的稀疏表示:αi=Sβi (L),结束迭代;并判断i个训练样本ti是否都找到稀疏表示αi;
当i个训练样本ti都已找到稀疏表示αi,则通过公式
采用加权求和将重叠的图像块组合,得出无噪RAW图像(5),式中,取λ=30/σ,S为对角矩阵,βi (L)为稀疏逼近,X为无噪图像,Y为去噪前图像,I为单位矩阵,Rij为图像块的提取矩阵,为Rij的转置,D为字典D(4);
当未全部找到i个样本ti的稀疏表示,则返回步骤5.1),继续执行;
5.3.2)当r(L)≥n(Cσ)2时,则返回步骤5.2),继续执行。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610846774.5A CN106447632B (zh) | 2016-09-23 | 2016-09-23 | 一种基于稀疏表示的raw图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610846774.5A CN106447632B (zh) | 2016-09-23 | 2016-09-23 | 一种基于稀疏表示的raw图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106447632A CN106447632A (zh) | 2017-02-22 |
CN106447632B true CN106447632B (zh) | 2019-04-02 |
Family
ID=58166916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610846774.5A Active CN106447632B (zh) | 2016-09-23 | 2016-09-23 | 一种基于稀疏表示的raw图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106447632B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107833194A (zh) * | 2017-11-21 | 2018-03-23 | 长沙全度影像科技有限公司 | 一种结合raw图像去噪的简单透镜图像复原方法 |
CN108550110B (zh) * | 2018-03-22 | 2021-03-02 | 影石创新科技股份有限公司 | 一种基于raw图像的全景拼接方法、装置及电子设备 |
CN111212198B (zh) * | 2018-11-22 | 2022-04-08 | 中移(杭州)信息技术有限公司 | 一种视频去噪的方法和设备 |
CN110189354B (zh) * | 2019-04-18 | 2021-12-28 | 北京迈格威科技有限公司 | 图像处理方法与图像处理器、图像处理设备及介质 |
CN112966131B (zh) * | 2021-03-02 | 2022-09-16 | 中华人民共和国成都海关 | 一种海关数据风控类型识别方法、海关智能化风险布控方法、装置、计算机设备及存储介质 |
CN116109519B (zh) * | 2023-04-04 | 2023-07-28 | 南京邮电大学 | 图像去噪方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1780405A (zh) * | 2005-05-26 | 2006-05-31 | 智多微电子(上海)有限公司 | 一种从原始拜尔插值到全彩色的转换方法和模块 |
CN102081795A (zh) * | 2011-01-14 | 2011-06-01 | 天津大学 | 基于稀疏表示的自动去块效应方法 |
CN103279932A (zh) * | 2013-06-05 | 2013-09-04 | 北京工业大学 | 一种二维合成稀疏模型和基于该模型的字典训练方法 |
CN104866905A (zh) * | 2015-04-27 | 2015-08-26 | 北京工业大学 | 一种基于beta过程的非参数稀疏张量字典的学习方法 |
-
2016
- 2016-09-23 CN CN201610846774.5A patent/CN106447632B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1780405A (zh) * | 2005-05-26 | 2006-05-31 | 智多微电子(上海)有限公司 | 一种从原始拜尔插值到全彩色的转换方法和模块 |
CN102081795A (zh) * | 2011-01-14 | 2011-06-01 | 天津大学 | 基于稀疏表示的自动去块效应方法 |
CN103279932A (zh) * | 2013-06-05 | 2013-09-04 | 北京工业大学 | 一种二维合成稀疏模型和基于该模型的字典训练方法 |
CN104866905A (zh) * | 2015-04-27 | 2015-08-26 | 北京工业大学 | 一种基于beta过程的非参数稀疏张量字典的学习方法 |
Non-Patent Citations (2)
Title |
---|
A Patch Aware Multiple Dictionary Framework for Demosaicing;Meiqing Zhang et.al;《 Computer Vision —ACCV 2014 》;20150416(第1期);第1-16页 |
一种基于压缩感知的图像去马赛克算法;康彬 等;《南京邮电大学学报(自然科学版)》;20140831;第34卷(第4期);第39-43页 |
Also Published As
Publication number | Publication date |
---|---|
CN106447632A (zh) | 2017-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106447632B (zh) | 一种基于稀疏表示的raw图像去噪方法 | |
CN109741256B (zh) | 基于稀疏表示和深度学习的图像超分辨率重建方法 | |
Li et al. | Underwater scene prior inspired deep underwater image and video enhancement | |
CN103093444B (zh) | 基于自相似性和结构信息约束的图像超分辨重建方法 | |
US20220301114A1 (en) | Noise Reconstruction For Image Denoising | |
CN110210282B (zh) | 一种基于非凸低秩稀疏分解的运动目标检测方法 | |
CN108269244B (zh) | 一种基于深度学习和先验约束的图像去雾*** | |
CN113658057B (zh) | 一种Swin Transformer微光图像增强方法 | |
WO2019186407A1 (en) | Systems and methods for generative ensemble networks | |
CN107085835B (zh) | 基于四元数加权核范数最小的彩色图像去噪方法 | |
CN116152120A (zh) | 一种融合高低频特征信息的低光图像增强方法及装置 | |
CN113284061B (zh) | 一种基于梯度网络的水下图像增强方法 | |
Parmar et al. | Spatio-spectral reconstruction of the multispectral datacube using sparse recovery | |
CN114820352A (zh) | 一种高光谱图像去噪方法、装置及存储介质 | |
Zhou et al. | IACC: cross-illumination awareness and color correction for underwater images under mixed natural and artificial lighting | |
CN116703752A (zh) | 融合近红外的Transformer结构的图像去雾方法及装置 | |
Li et al. | Efficient visual computing with camera raw snapshots | |
CN117422653A (zh) | 一种基于权重共享和迭代数据优化的低光照图像增强方法 | |
US11997246B2 (en) | Trained artificial intelligence model for raw to RGB image transformation | |
CN116843559A (zh) | 一种基于图像处理与深度学习的水下图像增强方法 | |
CN113191970B (zh) | 一种正交色彩传递网络及方法 | |
CN115861749A (zh) | 一种基于窗***叉注意力的遥感图像融合方法 | |
CN112862906B (zh) | 基于神经网络的颜色空间转换方法 | |
Bai et al. | Convolutional sparse coding for demosaicking with panchromatic pixels | |
CN111028159B (zh) | 一种图像条带噪声抑制方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |