CN104091314B - 基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法 - Google Patents

基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法 Download PDF

Info

Publication number
CN104091314B
CN104091314B CN201410348851.5A CN201410348851A CN104091314B CN 104091314 B CN104091314 B CN 104091314B CN 201410348851 A CN201410348851 A CN 201410348851A CN 104091314 B CN104091314 B CN 104091314B
Authority
CN
China
Prior art keywords
image
turbulence
spread function
point spread
edge
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
Application number
CN201410348851.5A
Other languages
English (en)
Other versions
CN104091314A (zh
Inventor
李晖晖
钱林弘
郭雷
杨宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201410348851.5A priority Critical patent/CN104091314B/zh
Publication of CN104091314A publication Critical patent/CN104091314A/zh
Application granted granted Critical
Publication of CN104091314B publication Critical patent/CN104091314B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)

Abstract

本发明涉及一种基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,技术特征在于:首先从当前复原图像中预测出有效的边缘,然后将边缘预测信息与自然图像边缘的稀疏先验信息相结合指导点扩散函数复原,最后通过一种非盲复原算法恢复出当前的目标图像,并将此复原图像作为下一次边缘预测的输入,如此迭代循环直到求出清晰的复原图像。本发明方法结合图像的先验信息与退化图像自身包含的有效信息,能有效抑制图像复原过程中产生的伪迹并恢复出更多的细节,故复原效果更好。

Description

基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原 方法
技术领域
本发明涉及一种基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,即基于边缘预测和稀疏正则项约束的湍流退化图像盲复原。将稀疏表示理论应用于湍流图像盲复原领域,发明成果可应用于各类军事或民用的图像处理系统中。
背景技术
大气湍流是造成天文观测图像退化的主要原因,湍流造成空气折射率的变化会导致光波振幅和相位的随机起伏,形成光强闪烁、波面畸变和光束漂移等现象。湍流对光波传播影响的本质就是改变了光波的原始波阵面,波阵面由平面畸变成随机曲面,这就使目标在成像设备焦平面的图像产生了严重模糊。观测到的图像如果不进行复原处理,就很难进行下一步的分类、识别等操作。为了解决这一难题,国内外学者提出一系列复原算法,根据点扩散函数(Point spread function,PSF)是否已知,主要将复原算法归为三大类,解卷积复原方法、盲解卷积复原方法和介于两者之间的近视解卷积复原方法。其中,解卷积方法假设退化过程的PSF是完全已知的,其中包括逆滤波,维纳滤波,卡尔曼滤波等经典复原算法,此类方法需要得到确切的退化模型,但是在实际观测活动中,造成观测目标模糊的点扩散函数变化是复杂并难于确定的,这就导致解卷积方法的实际复原效果并不理想,其实用价值不高;而盲解卷积算法则假设PSF是完全未知的,该类算法包括IBD(iterative blinddeconvolution)算法,NAS-RIF(nonnegativity and support recursive inversefiltering)算法等经典复原算法,IBD方法的计算复杂性较低,复原速度较快,但对噪声非常敏感,主要缺点是缺乏可靠性,其单值性和收敛性也是不确定的。另外,其复原对图像初始估计非常敏感,表现出不稳定性,NAS-RIF算法主要适用于支持域有限且背景较单一的空间目标湍流退化图像复原,具有比IBD算法更高的可靠性、更快的收敛性、更低的计算复杂性,但对噪声敏感,实际上还放大噪声;近似解卷积算法则是假设PSF不是完全已知也不是完全未知的,该类算法主要是将波前信息作为先验知识引入到复原算法中对图像进行复原。但由于测量设备的技术问题,波前测量噪声和误差以及是否与当前退化图像同步是一个难以解决的问题。而在实际应用中,点扩散函数往往是未知的,湍流图像复原问题被认为是盲解卷积问题,由于盲复原问题是一个严重的病态问题,病态问题的解常常是不稳定且非唯一的,因此在求解过程中常常需要加入合适的约束条件,缩小解空间,最终求得唯一稳定的解,这便是正则化技术。最初的正则化以l2范数作为约束条件,由于它的简单性,得到了广泛的应用,但是l2范数是能量的度量单位,无法证明用它求得的解是最优解。相反,在很多情况下甚至对解的结果具有误导性。近年来,研究表示自然图像的边缘具有某种稀疏特性,学者们考虑将这种稀疏特性作为先验知识引入到正则项约束中并用l0范数来衡量这种稀疏性,即将l0范数作为约束条件应用到目标函数中,它要求方程的解具有最小数目的非零项,由于解l0范数被证明是NP难问题,因此在实际应用中常使用l1范数来近似l0范数,因为l1范数最小化问题是一个凸优化问题,同时能保证所求的解是趋于稀疏的最优解,因此被广泛应用到运动模糊图像复原问题中。2011年,Krishnan研究指出l1/l2范数比l1范数更接近l0范数,用l1/l2作为约束条件,能求出更趋于真实的解,但是Krishnan利用l1/l2作为约束项所求的解有严重的伪迹。虽然近年来基于稀疏正则项约束的图像复原算法被广泛使用,但到目前为止大部分研究集中在运动模糊图像复原中。并且大部分基于稀疏正则项约束的方法直接将退化图像的梯度图像作为自然图像的边缘指导点扩散函数复原,由于退化图像存在比较严重的模糊和噪声,所求的梯度图像必然包含许多伪边缘,这些伪边缘会误导点扩散函数复原结果,从而导致复原的图像包含较多的伪迹。
以往的图像盲复原算法在退化图像噪声或模糊程度比较严重的情况下,容易使复原图像出现严重的伪迹,使复原图像难以获得令人满意的效果。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法。
技术方案
一种基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,其特征在于步骤如下:
步骤1、对湍流退化图像Y和点扩散函数k进行多尺度级数分解:首先对大小为K_size的点扩散函数进行下抽样,直到点扩散函数缩小到规定尺寸K_min,下抽样次数为L;然后,对湍流退化图像进行L次下抽样,得到L级图像,将分辨率最低的记为第1级图像,分辨率最高的记为第L级图像;
步骤2:对当前复原图像进行边缘预测:对第1级的湍流退化图像进行边缘预测,并经过双边滤波器进行滤波,得到图像XL;然后将XL经过冲击滤波器进行滤波,得到图像为XS;最后将XS映射到梯度域,以梯度阈值对XS的梯度图像进行截断处理,剔除无效的边缘,输出的边缘为▽′X={Px,Py};
步骤3、对点扩散函数k进行估计:
采用代价函数对点扩散函数k进行估计,其中λ,ψ表示权重,表示卷积算子,▽={▽x,▽y}表示梯度算子,▽x为x方向上的梯度算子,▽y为y方向上的梯度算子;
将代价函数拆分为:
边缘采用iterative shrinkage-thresholding algorithm(ISTA)算法对其进行求解;
点扩散函数k:采用unconstrained iterativere-weighted least squares算法对其进行求解,并将所求的解映射到约束k≥0,∑iki=1上;
将上述两个解▽′X和k进行交替迭代,直到达到内层循环迭代次数iter_in为止,得出估计的点扩散函数k;
步骤4、对湍流退化图像Y进行复原:采用公式选择非盲解卷积的方法对湍流图像进行复原,求得的X为当前复原图像,其中τ表示权重;
将复原图像X作为下一次迭代中边缘预测的输入图像,重复步骤2~步骤4,直到达到外层循环迭代次数iter_out为止;
对每一级图像执行步骤2~步骤4,并将上一级的复原图像X、湍流退化图像Y和点扩散函数k进行上抽样,作为下一级复原算法的初始输入,直到达到最大级数L为止。
所述梯度阈值的确定为:将图像XS的梯度方向分为4个方向,将每个方向的梯度大小按照从大到小的顺序排列,寻找一个阈值使得每个方向至少有个最大的梯度值;其中Pk为点扩散函数的总像素值。
所述双边滤波器的半宽N=2,σs=2,σr=5。
所述抽样比例为
所述点扩散函数缩小到规定尺寸K_min=5×5。
所述点扩散函数大小K_size=35×35。
有益效果
本发明提出的一种基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,对当前复原图像进行边缘预测,用预测的边缘结合稀疏比值正则项约束复原算法指导点扩散函数复原,将求得的点扩散函数代入非盲复原算法中求出复原图像,然后将当前复原的图像应用于下一次迭代的边缘预测中。循环迭代此过程,每迭代一次,所求复原图像的边缘就更接近真实图像的边缘,从而对点扩散函数估计进行修正,使其向真实的值靠近,如此迭代下去最终得出理想的复原图像。
本发明利用稀疏约束指导图像复原,将自然图像稀疏特性作为先验知识加入到正则项约束中,使解空间向该先验分布靠近,能有效的复原出符合先验特性的自然图像,对于不符合先验分布的图像则会产生严重的伪迹,而边缘预测信息能提取图像本身所含的有效信息,将边缘预测与稀疏正则项约束结合共同对解空间进行约束,从而有效的减小了复原图像的伪迹,最终使得所求的解具有较好的复原效果。
附图说明
图1:本发明方法复原图像的流程图
图2:采用不同复原算法对仿真卫星图像1进行盲复原结果
(a)源图像;(b)湍流退化图像;(c)IBD算法复原结果;(d)Shan[2008]所提算法复原结果;(e)Krishnan[2011]所提算法复原结果;(f)本发明所提算法复原结果;
图3:采用不同复原算法对仿真卫星图像2进行盲复原结果
(a)源图像;(b)湍流退化图像;(c)IBD算法复原结果;(d)Shan[2008]所提算法复原结果;(e)Krishnan[2011]所提算法复原结果;(f)本发明所提算法复原结果。
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明方法的特征在于步骤如下:
步骤1对湍流退化图像Y和点扩散函数k进行多尺度级数分解:首先对大小为K_size的点扩散函数按一定比例进行下抽样,直到点扩散函数缩小到规定尺寸K_min,记录下抽样次数L;然后,对湍流退化图像进行相同的操作得到L级图像,分辨率最低的记为第1级图像(依次类推),分辨率最高的记为第L级图像。
步骤2对当前复原图像进行边缘预测:将第1级的湍流退化图像作为边缘预测的初始输入图像,将其经过双边滤波器进行滤波,输出图像为XL;然后将XL经过冲击滤波器进行滤波,输出图像为XS;最后将XS映射到梯度域,并设置梯度阈值对XS的梯度图像进行截断处理,剔除无效的边缘,输出的边缘为▽′X={Px,Py}。
(a)双边滤波
双边滤波器用于平滑当前复原图像X的噪声并保持其边缘,它采用局部加权平均的方法获取输出图像XL的像素值,数学模型如下:
X L ( x , y ) = Σ ( i , j ∈ S x , y ) w ( i , j ) X ( i , j ) Σ ( i , j ∈ S x , y ) w ( i , j ) - - - ( 1 )
式(1)中,Sx,y表示中心点(x,y)的(2N+1)×(2N+1)大小的邻域,N表示滤波器半宽,w(i,j)为加权系数,X(i,j)为当前复原图像X的Sx,y邻域内的像素值,XL(x,y)为中心点(x,y)处的像素输出值,i表示图像像素点的行坐标,j表示图像像素点的列坐标。式(1)右边是对中心像素点邻域内像素亮度值的加权平均。对邻域内的每一个像素点X(i,j),其加权系数w(i,j)由空间邻近度因子ws和亮度像似度因子wr的乘积构成:
w s ( i , j ) = e - | i - x | 2 + | j - y | 2 2 σ s 2 - - - ( 2 )
w r ( i , j ) = e - | X ′ ( i , j ) - X ′ ( x , y ) | 2 2 σ r 2 - - - ( 3 )
因此,
wi,j=ws(i,j)wr(i,j) (4)
σs表示几何扩散参数,σr表示光度扩散参数,双边滤波器受3个参数的控制:滤波器半宽N、参数σs和σr。N越大,平滑作用越强;σs和σr分别控制着空间邻近度因子ws和亮度像似度因子wr的衰减程度。
(b)冲击滤波
冲击滤波器是增强图像特征的一个有效的工具,它从模糊的图像XL中恢复尖锐的边缘,演化方程如下:
XLt+1=XLt-sign(ΔXLt)||▽XLt||dt (5)
XLt表示t时刻的图像,Δ表示拉普拉斯算子,dt表示时间步长。其中
▿ XL = ( XL x , XL y ) , ΔXL = XL x 2 XL xx + 2 XL x XL y XL xy + XL y 2 XL yy - - - ( 6 )
XLx表示图像XL在x方向的一阶导数,XLy表示图像XL在y方向的一阶导数,XLxx表示图像XL在x方向的二阶导数,XLyy表示图像XL在y方向的二阶导数,XLxy表示图像XLx在y方向的一阶导数,最终所求得的XLt+1即为冲击滤波器的输出图像XS。
(c)梯度截断
梯度阈值的选择方法如下:将图像XS的梯度方向分为4个方向,将每个方向的梯度大小按照从大到小的顺序排列,寻找一个阈值使得每个方向至少有个最大的梯度值;其中Pk为点扩散函数的总像素值。
步骤3对点扩散函数k进行估计
采用代价函数对点扩散函数k进行估计,其中λ,ψ表示权重,表示卷积算子,▽={▽x,▽y}表示梯度算子,▽x为x方向上的梯度算子,▽y为y方向上的梯度算子。
将代价函数拆分为:
①边缘采用iterative shrinkage-thresholding algorithm(ISTA)算法对其进行求解;
②点扩散函数k:采用unconstrainediterative re-weighted least squares(IRLS)算法对其进行求解,并将所求的解映射到约束k≥0,∑iki=1上;
将上述两个公式的解▽′X和k进行交替迭代,直到达到内层循环迭代次数iter_in为止,得出估计的点扩散函数k。
步骤4对湍流退化图像Y进行复原
选择非盲解卷积的方法对湍流图像进行复原,采用公式求得的X即为当前复原图像,其中τ表示权重。将复原图像X作为下一次迭代中边缘预测的输入图像,重复步骤2~步骤4,直到达到外层循环迭代次数iter_out为止。
对每一级图像执行步骤2~步骤4算法,并将上一级的复原图像X、湍流退化图像Y和点扩散函数k进行上抽样,作为下一级复原算法的初始输入,直到达到最大级数L为止。
具体实施例:
用于实施的硬件环境是:Pentium-4 3G计算机、2GB内存、128M显卡,运行的软件环境是:Mat1ab R2012b和windows XP。采用Matlab程序设计语言实现了本发明提出的新算法。图像数据采用两幅256×256的卫星仿真图像,通过模拟大气湍流的相位屏,对卫星图像进行湍流退化模糊的仿真实验,本实验设定大气相干长度r0=0.05m,望远镜口径直径D=1.0m,将得到的点扩散函数与原始卫星图像进行卷积运算,再施加高斯随机噪声(方差为0.18)即得到实验仿真的湍流退化图像,复原时设定点扩散函数支持域为35×35。
本发明具体实施如下:
步骤1对湍流退化图像Y和点扩散函数k进行多尺度级数分解:首先对大小为K_size的点扩散函数按一定比例进行下抽样,直到点扩散函数缩小到规定尺寸K_min,记录下抽样次数L;然后,对湍流退化图像进行相同的操作得到L级图像,分辨率最低的记为第1级图像(依次类推)。本实验设定点扩散函数大小K_size=35×35,点扩散函数最小尺寸K_min=5×5,最大级数L=6,上抽样与下抽样均选择双线性插值方法,抽样比例为
步骤2对当前复原图像进行边缘预测:将第1级的湍流退化图像作为边缘预测的初始输入图像,将其经过双边滤波器进行滤波,输出图像为XL;然后将XL经过冲击滤波器进行滤波,输出图像为XS;最后将XS映射到梯度域,并设置梯度阈值对XS的梯度图像进行截断处理,剔除无效的边缘,输出的边缘为▽′X={Px,Py}。
(a)双边滤波
双边滤波器用于平滑当前复原图像X的噪声并保持其边缘,它采用局部加权平均的方法获取输出图像XL的像素值,数学模型如下:
X L ( x , y ) = Σ ( i , j ∈ S x , y ) w ( i , j ) X ( i , j ) Σ ( i , j ∈ S x , y ) w ( i , j ) - - - ( 7 )
式(7)中,Sx,y表示中心点(x,y)的(2N+1)×(2N+1)大小的邻域,N表示滤波器半宽,w(i,j)为加权系数,X(i,j)为当前复原图像X的Sx,y邻域内的像素值,XL(x,y)为中心点(x,y)处的像素输出值,i表示图像像素点的行坐标,j表示图像像素点的列坐标。式(7)右边是对中心像素点邻域内像素亮度值的加权平均。对邻域内的每一个像素点X(i,j),其加权系数w(i,j)由空间邻近度因子ws和亮度像似度因子wr的乘积构成:
w s ( i , j ) = e - | i - x | 2 + | j - y | 2 2 σ s 2 - - - ( 8 )
w r ( i , j ) = e - | X ′ ( i , j ) - X ′ ( x , y ) | 2 2 σ r 2 - - - ( 9 )
因此,
wi,j=ws(i,j)wr(i,j) (10)
σs表示几何扩散参数,σr表示光度扩散参数,双边滤波器受3个参数的控制:滤波器半宽N、参数σs和σr。N越大,平滑作用越强;σs和σr分别控制着空间邻近度因子ws和亮度像似度因子wr的衰减程度。本实验所选择的双边滤波器的半宽N=2,σs=2,σr=5。
(b)冲击滤波
冲击滤波器是增强图像特征的一个有效的工具,它从模糊的图像XL中恢复尖锐的边缘,演化方程如下:
XLt+1=XLt-sign(ΔXLt)||▽XLt||dt (11)
XLt表示t时刻的图像,Δ表示拉普拉斯算子,dt表示时间步长。其中
▿ X L = ( XL x , XL y ) , Δ X L = XL x 2 XL x x + 2 XL x XL y XL x y + XL y 2 XL y y - - - ( 12 )
XLx表示图像XL在x方向的一阶导数,XLy表示图像XL在y方向的一阶导数,XLxx表示图像XL在x方向的二阶导数,XLyy表示图像XL在y方向的二阶导数,XLxy表示图像XLx在y方向的一阶导数,最终所求得的XLt+1即为冲击滤波器的输出图像XS。本实验所选的梯度算子为▽x=[0,-1,1],▽y=[0,-1,1]T,dt=1.0,随着迭代次数的增加,外层循环每迭代一次参数σr和dt分别除以1.1;
(c)梯度截断
梯度阈值的选择方法如下:将图像XS的梯度方向分为4个方向,将每个方向的梯度大小按照从大到小的顺序排列,寻找一个阈值使得每个方向至少有个最大的梯度值;其中Pk为点扩散函数的总像素值。根据每个尺度下点扩散函数的大小可以计算出点扩散函数的总像素值Pk
步骤3对点扩散函数k进行估计
采用代价函数对点扩散函数k进行估计,其中λ,ψ表示权重,表示卷积算子,▽={▽x,▽y}表示梯度算子,▽x为x方向上的梯度算子,▽y为y方向上的梯度算子。本实验中设置权重λ=70,ψ=1×10-5
将代价函数拆分为:
①边缘采用iterative shrinkage-thresholding algorithm(ISTA)算法对其进行求解;
②点扩散函数采用unconstrainediterative re-weighted least squares(IRLS)算法对其进行求解,并将所求的解映射到约束k≥0,∑iki=1上;
将上述两个公式的解▽′X和k进行交替迭代,直到达到内层循环迭代次数iter_in为止,得出估计的点扩散函数k。本实验中设置内层循环迭代次数iter_in=40。
步骤4对湍流退化图像Y进行复原
选择非盲解卷积的方法对湍流图像进行复原,采用公式求得的X即为当前复原图像,其中τ表示权重。将复原图像X作为下一次迭代中边缘预测的输入图像,重复步骤2~步骤4,直到达到外层循环迭代次数iter_out为止,本实验中设置权重τ=1000,外层循环迭代次数iter_out=15。
对每一级图像执行步骤2~步骤4算法,并将上一级的复原图像X、湍流退化图像Y和点扩散函数k进行上抽样,作为下一级复原算法的初始输入,直到达到最大级数L为止。
选用图像的峰值信噪比(Peak Signal to Noise Ratio,PSNR)为客观评价指标,定量的评价不同图像采用不同算法的复原效果。客观评价结果如表1所示。其中峰值信噪比的定义式如式(13)所示。
PSNR=10*log(2552/MSE) (13)
式(21)中,*代表乘法,MSE代表均方误差,其定义式如式(14)所示。
M S E = Σ i = 1 m Σ j = 1 n [ f ( i , j ) - f ^ ( i , j ) ] 2 m × n - - - ( 14 )
式(14)中m和n分别表示图像的长度和宽度,f(i,j)和分别表示理想图像和被评价图像在像素点(i,j)处的灰度值。
表1不同算法复原图像的峰值信噪比

Claims (6)

1.一种基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,其特征在于步骤如下:
步骤1、对湍流退化图像Y和点扩散函数k进行多尺度级数分解:首先对大小为K_size的点扩散函数进行下抽样,直到点扩散函数缩小到规定尺寸K_min,下抽样次数为L;然后,对湍流退化图像进行L次下抽样,得到L级图像,将分辨率最低的记为第1级图像,分辨率最高的记为第L级图像;
步骤2:对当前复原图像进行边缘预测:对第1级的湍流退化图像进行边缘预测,并经过双边滤波器进行滤波,得到图像XL;然后将XL经过冲击滤波器进行滤波,得到图像为XS;最后将XS映射到梯度域,以梯度阈值对XS的梯度图像进行截断处理,剔除无效的边缘,输出的边缘为
步骤3、对点扩散函数k进行估计:
采用代价函数对点扩散函数k进行估计,其中λ,ψ表示权重,表示卷积算子,表示梯度算子,为x方向上的梯度算子,为y方向上的梯度算子;
将代价函数拆分为:
边缘 采用iterative shrinkage-thresholdingalgorithm(ISTA)算法对其进行求解;
点扩散函数k:采用unconstrained iterative re-weighted least squares算法对其进行求解,并将所求的解映射到约束k≥0,∑iki=1上;
将上述两个解和k进行交替迭代,直到达到内层循环迭代次数iter_in为止,得出估计的点扩散函数k;
步骤4、对湍流退化图像Y进行复原:采用公式选择非盲解卷积的方法对湍流图像进行复原,求得的X为当前复原图像,其中τ表示权重;
将复原图像X作为下一次迭代中边缘预测的输入图像,重复步骤2~步骤4,直到达到外层循环迭代次数iter_out为止;
对每一级图像执行步骤2~步骤4,并将上一级的复原图像X、湍流退化图像Y和点扩散函数k进行上抽样,作为下一级复原算法的初始输入,直到达到最大级数L为止。
2.根据权利要求1所述基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,其特征在于:所述梯度阈值的确定为:将图像XS的梯度方向分为4个方向,将每个方向的梯度大小按照从大到小的顺序排列,寻找一个阈值使得每个方向至少有个最大的梯度值;其中Pk为点扩散函数的总像素值。
3.根据权利要求1所述基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,其特征在于:所述双边滤波器的半宽N=2,σs=2,σr=5。
4.根据权利要求1所述基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,其特征在于:所述抽样比例为
5.根据权利要求1所述基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,其特征在于:所述点扩散函数缩小到规定尺寸K_min=5×5。
6.根据权利要求1所述基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法,其特征在于:所述点扩散函数大小K_size=35×35。
CN201410348851.5A 2014-07-22 2014-07-22 基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法 Active CN104091314B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410348851.5A CN104091314B (zh) 2014-07-22 2014-07-22 基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410348851.5A CN104091314B (zh) 2014-07-22 2014-07-22 基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法

Publications (2)

Publication Number Publication Date
CN104091314A CN104091314A (zh) 2014-10-08
CN104091314B true CN104091314B (zh) 2017-02-01

Family

ID=51639029

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410348851.5A Active CN104091314B (zh) 2014-07-22 2014-07-22 基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法

Country Status (1)

Country Link
CN (1) CN104091314B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104392415A (zh) * 2014-11-20 2015-03-04 杭州电子科技大学 一种结合权重因子与梯度制约的图像恢复方法
CN104574315B (zh) * 2015-01-08 2017-03-22 哈尔滨工业大学 基于光强传输矩阵的光学系统成像复原方法
CN104835126B (zh) * 2015-05-15 2016-03-02 中国人民解放军国防科学技术大学 一种基于观测信号拆分的快速稀疏图像复原方法
CN106920220B (zh) * 2017-02-24 2019-07-19 西北工业大学 基于暗原色和交替方向乘子法优化的湍流图像盲复原方法
CN108230253B (zh) * 2017-05-08 2020-11-27 北京市商汤科技开发有限公司 图像恢复方法、装置、电子设备和计算机存储介质
CN107742278B (zh) * 2017-10-25 2021-07-13 重庆邮电大学 结合l0范数和空间尺度信息的运动模糊图像盲复原方法
CN108305230A (zh) * 2018-01-31 2018-07-20 上海康斐信息技术有限公司 一种模糊图像综合处理方法和系统
CN108765288B (zh) * 2018-05-25 2022-05-31 杭州电子科技大学 一种基于边缘保持的pocs图像超分辨重建方法
CN108921132A (zh) * 2018-07-27 2018-11-30 广东电网有限责任公司 无人机机载电缆探测系统
CN110223242B (zh) * 2019-05-07 2021-08-10 北京航空航天大学 一种基于时空域残差生成对抗网络的视频湍流移除方法
CN110400280B (zh) * 2019-08-02 2023-02-03 电子科技大学 一种基于人造信标和相位屏的大气湍流退化图像复原方法
CN110796616B (zh) * 2019-10-23 2022-05-10 武汉工程大学 基于范数约束和自适应加权梯度的湍流退化图像恢复方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101487895A (zh) * 2008-01-15 2009-07-22 华晶科技股份有限公司 显示鸟瞰车辆图像的倒车雷达系统
CN103201765A (zh) * 2010-09-28 2013-07-10 马普科技促进协会 用于从所观察到的数字图像序列恢复数字图像的方法与设备

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7120308B2 (en) * 2001-11-26 2006-10-10 Seiko Epson Corporation Iterated de-noising for image recovery

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101487895A (zh) * 2008-01-15 2009-07-22 华晶科技股份有限公司 显示鸟瞰车辆图像的倒车雷达系统
CN103201765A (zh) * 2010-09-28 2013-07-10 马普科技促进协会 用于从所观察到的数字图像序列恢复数字图像的方法与设备

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
湍流退化图像复原技术研究现状及展望;王亮亮 等;《飞行器测控学报》;20091231;第28卷(第6期);第63-67页 *

Also Published As

Publication number Publication date
CN104091314A (zh) 2014-10-08

Similar Documents

Publication Publication Date Title
CN104091314B (zh) 基于边缘预测和稀疏比值正则约束的湍流退化图像盲复原方法
Li et al. Image dehazing using residual-based deep CNN
Mao et al. Non rigid geometric distortions correction-application to atmospheric turbulence stabilization
Chan et al. An augmented Lagrangian method for total variation video restoration
Xie et al. Removing turbulence effect via hybrid total variation and deformation-guided kernel regression
JP6727642B2 (ja) 学習型アルゴリズムによるフォーカス補正処理方法
CN110796616A (zh) 基于分数阶微分算子的l0范数约束和自适应加权梯度的湍流退化图像恢复方法
Flamary Astronomical image reconstruction with convolutional neural networks
Bajpai et al. Analysis of image enhancement techniques used in remote sensing satellite imagery
Chen et al. Hyperspectral image denoising with weighted nonlocal low-rank model and adaptive total variation regularization
Fu et al. An anisotropic Gaussian filtering model for image de-hazing
Ljubenović et al. CNN-based deblurring of terahertz images
Wei et al. Fast space-varying convolution using matrix source coding with applications to camera stray light reduction
Gu et al. A novel total generalized variation model for image dehazing
De Dravo et al. Stress for dehazing
CN114298950A (zh) 一种基于改进的GoDec算法的红外与可见光图像融合方法
CN110246152A (zh) Piv图像处理方法及系统
CN114066786A (zh) 一种基于稀疏和滤波器的红外和可见光图像融合方法
Chen et al. Blind restoration for nonuniform aerial images using nonlocal Retinex model and shearlet-based higher-order regularization
Wong et al. Regularization-based modulation transfer function compensation for optical satellite image restoration using joint statistical model in curvelet domain
Rawshan Habib et al. Study on Advanced Image Processing Techniques for Remote Sensor Data Analysis
Xu et al. Geometric distortion correction of long-range imaging containing moving objects
Kaur et al. Study of Image Enhancement Techniques in Image Processing: A
Yu et al. A novel nighttime dehazing model integrating Retinex algorithm and atmospheric scattering model
Gunturk et al. Fundamentals of image restoration

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant