CN105631831A - 一种雾霾条件下视频图像增强方法 - Google Patents

一种雾霾条件下视频图像增强方法 Download PDF

Info

Publication number
CN105631831A
CN105631831A CN201610143521.1A CN201610143521A CN105631831A CN 105631831 A CN105631831 A CN 105631831A CN 201610143521 A CN201610143521 A CN 201610143521A CN 105631831 A CN105631831 A CN 105631831A
Authority
CN
China
Prior art keywords
value
image
frame
atmospheric light
current frame
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.)
Granted
Application number
CN201610143521.1A
Other languages
English (en)
Other versions
CN105631831B (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201610143521.1A priority Critical patent/CN105631831B/zh
Publication of CN105631831A publication Critical patent/CN105631831A/zh
Application granted granted Critical
Publication of CN105631831B publication Critical patent/CN105631831B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Processing Of Color Television Signals (AREA)
  • Picture Signal Circuits (AREA)

Abstract

本发明公开了一种视频图像去雾处理方法,能够在满足实时去雾功能的同时,有效消除去雾后视频图像闪烁。具体步骤为:通过暗通道以及自适应调整的方法求取各帧图像的大气光值;暗通道先验算法计算第一帧图像的初始透射率;构造帧间的透射率估计和邻域空间的能量函数,求解块区域最优透射率值;快速导向滤波算法细化对初始透射率;通过去雾模型求去雾后的视频图像。

Description

一种雾霾条件下视频图像增强方法
技术领域
本发明涉及数字图像处理领域,尤其涉及一种雾霾条件下视频图像增强方法,适用于机器视觉应用的前期图像增强预处理,可广泛应用于在雾霾天气条件下的道路交通监控、行车记录、安防监控、无人机地面侦察等视频图像应用领域。
背景技术
近年来,我国的大气污染越来越严重,由此产生的雾霾导致大气能见度降低,户外视频监视系统很难直接得到具有较高清晰度的图像信息,直接导致监视系统失效。为此,在摄像机得到的视频图像基础上利用软件算法进行视频图像去雾处理,增加图像清晰度具有重要意义。
在过去的十几年里,国内外学者针对单幅图像去雾开展了很多研究,并提出了多种有效的图像去雾方法。而在视频图像去雾方面,由于存在视频序列相邻帧之间高冗余性数据应用、各帧图像透射率信息获取以及视频编辑图像时空一致性保持等特殊问题,相关研究工作相对单幅图像去雾研究较为缺乏。目前已提出的视频图像去雾方法包括两种类型。一种方法是针对具有固定背景的有雾视频,将其背景进行一次去雾处理,后续帧背景采用前帧去雾后背景替换;对其余区域进行逐帧去雾。该方法运算量相对较少,但仅适用于存在固定背景的视频应用。另一种方法是采用光流法进行运动估计,然后利用视频图像前后帧的透射率信息及运动估计结果构造马尔科夫场,通过求解全局最优值来获得本帧透射率。此种方法利用了相邻帧之间高冗余性的图像数据,保持了时空一致性,但由于求解全局最优需要解算大型稀疏方程,运算量极大,不适合于具有实时性要求的视频应用场合。
发明内容
有鉴于此,本发明提供了一种雾霾条件下视频图像增强方法,通过利用相邻图像帧的数据信息,既保持了时间连续性且处理速度能够满足视频处理的实时性要求,满足实时去雾功能的同时,能够有效地消除去雾后视频图像闪烁。
为了解决上述技术问题,本发明是这样实现的:
一种雾霾条件下视频图像增强方法,包括如下步骤:
步骤一、求取基于时间连续自适应调整的大气光值,具体为:
A、针对输入视频的当前帧图像Ik,获得每一个像素点的三个独立的单通道图像R、G、B中的最小值,遍历图像Ik的各个像素点,得到暗通道图Idark
将暗通道图Idark中像素亮度值从大到小排列的前0.1%的像素点筛选出来,其它像素点不考虑,形成大气光区域图;针对图像Ik中与大气光区域图对应位置的各个像素点,获得各个像素点通道R、通道G以及通道B的平均值,则得到的三个平均值作为大气光值,该大气光值为一个有R,G,B分量的向量。
B、对视频中的每一帧图像进行如a中所述的处理,获得各帧图像原始大气光值;
C、如果视频中当前帧为第一帧,则将该帧图像的原始大气光值作为经过学习得到的当前帧大气光值Ak,此时k=0;
D、如果视频中当前帧非第一帧,则当前帧的前N帧大气光值以及当前帧原始大气光值进行加权平均,作为经学习得到的当前帧大气光值Ak;其中,所述N至少取2,且小于或等于k-1;
步骤二、估计初始透射率图,具体包括如下步骤:
S201、将图像Ik各元素的值除以步骤一得到的经过学习的当前帧大气光值Ak,得到图像IA,将图像IA划分为半径为r的窗口,求每个窗口中的各个像素点的R、G和B三通道中的最小通道值,再从该窗口中各像素点的最小通道值中选择一个最小值,得到每个窗口对应的暗通道图I′A;将该暗通道图I′A与加权系数ω的乘积后,再用1减去该乘积,最终估计得到初始的透射率
其中,半径r的取值范围为7-15个像素;加权系数ω取值范围在0.90至0.99之间;
S202、判断:如果视频中的当前帧为第一帧,则将透射率作为经过调整得到的当前帧的初始透射率;
如果视频中当前帧非第一帧,则通过时间相干性方法获取调整后的初始透射率图,所述时间相干性方法求解透射率图的方法进一步包括:
将每一帧的RGB彩色图像转换到YUV空间,得到Y通道的图像Iy;将步骤一获得的经过学习的大气光值Ak转换为标量:
A=Ar·0.3+Ag·0.59+Ab·0.11(4)
其中Ar、Ag和Ab分别表示大气光值Ak的三通道分量值;
S203、第k帧透射率tk(p)和k-1帧的透射率tk-1(p)在像素点P的关系为tk(p)=τk(p)tk-1(p),则时间相干性因子根据下式求出:
τ k ( p ) = I k ( p ) - A I k - 1 ( p ) - A - - - ( 6 )
其中,Ik(p)和Ik-1(p)表示同一像素点P在第k帧和k-1帧Y通道的图像所具有单通道值,σ2表示模型的方差;所述σ2的取值为10;
S204、根据时间相干因子τk,得到tk(p)=τk(p)tk-1(p)后构造代价函数,如下式所述:
E = Σ p ∈ Ω ( τ k t k - 1 ( p ) - t k ) 2 σ t 2 + Σ q ∈ N p ( t ~ ( q ) - t k ) 2 max ( | | I k ( p ) - I k - 1 ( q ) | | 2 , a 2 ) - - - ( 8 )
其中,tk为最终要求的透射率,为未知数;σt 2表示透射率估计的模型误差,所述σt 2取值为10;Np表示像素点P的邻域像素点;tk-1(p)表示第k-1帧中像素点P的透射率;Ik-1(q)表示第k-1帧图像中像素点p的邻域像素点的Y通道图像的值;表示像素点p的邻域像素点的透射率,由暗通道求得的初始透射率;所述a2取值为0.0052
S205、求式(8)表示的代价E最小时的透射率tk
步骤三、将输入图像的YUV空间的Y通道图像作为导向图像,初始透射率图tk作为输入图,对初始透射率tk导向滤波,得到各个区域块最终的透射率t*
步骤四、针对输入的当前帧图像的各个像素点,通过以下公式得到去雾处理后的图像:
J ( p ) = I ( p ) - A k t * ( p ) + A k - - - ( 9 )
其中,t*(p)表示任意像素点p所在位置的最终的透射率,根据步骤3的最终的透射率t*获得;I(p)为输入的RGB彩色图像的像素点p的三通道值,Ak为步骤一得到的三通道大气光值。
较佳的,在所述步骤中,当前帧的前N帧大气光值以及当前帧原始大气光值进行加权平均的具体方法如下:
找到大气光区域图中亮度最小的像素点,求该像素点3×3窗口内对应图像Ik中的像素值的差值Bstd和平均值Bmean,再求该差值Bstd与平均值Bmean之间的差值:
L=|Bstd-Bmean|(1);
最后,判断当前第k帧对应的差值Lk与第k-1帧的差值Lk-1之间的差值是否小于或等于预设的阈值Lthresh,若第k帧Lk与第k-1帧Lk-1差值小于或等于Lthresh,则将当前帧原始大气光值作为经过学习得到的当前帧大气光值Ak;否则将当前帧的前N帧大气光值以及当前帧原始大气光值进行加权平均,作为经学习得到的当前帧大气光值Ak;其中,所述Lthresh取值为0.02;
所述经学习得到的当前帧大气光值Ak的R、G和B任意一个通道的值为:
A k = A ^ 0 ( k = 0 ) 1 N + 1 Σ j = k - N k a j A ^ j ( k > 0 , | L k - L k - 1 | > L t h r e s h ) A ^ k ( k > 0 , | L k - L k - 1 | ≤ L t h r e s h ) - - - ( 2 )
其中,分别表示原始第k帧和第j帧大气光,Ak表示经过学习获得的第k帧大气光,aj表示当前帧更新权重;所述权重aj的取值原则为,与当前帧越近,取值越大,反之越小;
较佳的,所述r取15。
较佳的,ω取值范围在0.94至0.96之间。
较佳的,ω取值0.95。
较佳的,所述步骤二中,得到时间相干性因子后,再得到基于区域块的时间相干因子,具体为:
首先,计算图像Ik中任意像素点P在相邻帧出现在同一个位置的概率,一种简单的方法,即下式中的权重值wk(p):
w k ( p ) = exp ( - ( I k ( p ) - I k - 1 ( p ) ) 2 σ 2 ) - - - ( 5 )
其次,对于图像Ik中的任意像素点P,以像素点P为中心,r’为半径划分区域块Ω,r’取值为5至15个像素;采用下式得到基于块的时间相干因子
τ k ‾ = Σ p ∈ Ω w k ( p ) τ k ( p ) Σ p ∈ Ω w k ( p ) - - - ( 7 )
最后,将该基于块的时间相干因子代替时间相干因子τk代入到代价E的表达式中,最终求得代价E最小时的透射率tk
较佳的,所述步骤二中,求代价E最小时的透射率tk的方法具体包括:
设在每一个区域块中透射tk初始值为0.1,代入代价函数求代价值E1;然后将透射率tk逐次增加0.1,每次均求出对应的代价值;选出最小代价值,此时对应的tk即为满足要求的最小能量区域块透射率值;其中,累加次数小于10次;
每一个区域块均重复上述步骤,求出每一个区域块对应的最佳的tk,最终得到经过调整的整个图像透射率图。
本发明具有如下有益效果:
(1)基于时间连续的自适应调整获得大气光值,解决了现有方法存在的各帧大气光跳变问题;建立了相邻帧透射率图关系模型,通过求解区域最优透射率来保证帧间去雾画面的时间连续性,消除了去雾后视频图像的闪烁;在调整初始透射率和引导滤波时只考虑Y通道,增强了视频图像去雾的实时性。
(2)基于点的时间相干性因子是在假设像素点的场景反射率为常值的情况下得到的数值,在镜头移动和有运动物体情景下,在整幅图像中并不完全适用。引入基于块的时间相干性因子,能更准确的描述块与块之间的相似性,在运动情况下也能适用。
附图说明
图1为本发明的流程图;
图2(a)为去雾前的输入图像。
图2(b)为输入图像采用本发明的方法进行去雾后的效果。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
步骤一、基于时间连续自适应调整的大气光值求取。
a.针对输入视频的当前帧图像Ik,获得每一个像素点的三个独立的单通道图像R、G、B中的最小值,遍历图像Ik的各个像素点,得到暗通道图Idark
将暗通道图Idark中像素亮度值从大到小排列的前0.1%的像素点筛选出来,其它像素点不考虑,形成大气光区域图;针对图像Ik中与大气光区域图对应位置的各个像素点,获得各个像素点通道R、通道G以及通道B的平均值,则得到的三个平均值作为大气光值,该大气光值为一个有R,G,B分量的向量。
b.对视频中的每一帧图像进行如a中所述的处理,获得各帧图像原始大气光值。
c.如果视频中当前帧为第一帧,则将该帧图像的原始大气光值作为经过学习得到的当前帧大气光值Ak,此时k=0;
d.如果视频中当前帧非第一帧,则进行如下处理:找到大气光区域图中亮度最小的像素点,求该像素点3×3窗口内对应图像Ik中的像素值的差值Bstd和平均值Bmean,再求该差值Bstd与平均值Bmean之间的差值:
L=|Bstd-Bmean|(1);
最后,判断当前第k帧对应的差值Lk与第k-1帧的差值Lk-1之间的差值是否小于或等于预设的阈值Lthresh,若第k帧Lk与第k-1帧Lk-1差值小于或等于Lthresh,则将当前帧原始大气光值作为经过学习得到的当前帧大气光值Ak;否则将当前帧的前N帧大气光值以及当前帧原始大气光值进行加权平均,作为经学习得到的当前帧大气光值Ak;其中,所述Lthresh取值为0.02,所述N至少取2,且小于或等于k-1,即当前帧之前的所有帧。
所述经学习得到的当前帧大气光值Ak的R、G和B任意一个通道的值为:
A k = A ^ 0 ( k = 0 ) 1 N + 1 Σ j = k - N k a j A ^ j ( k > 0 , | L k - L k - 1 | > L t h r e s h ) A ^ k ( k > 0 , | L k - L k - 1 | ≤ L t h r e s h ) - - - ( 2 )
其中,分别表示原始第k帧和第j帧大气光,Ak表示经过学习获得的第k帧大气光,aj表示当前帧更新权重。
所述权重aj的取值原则为,与当前帧越近,取值越大,反之越小。本实施例中,aj取值为0.5,aj-1取值为0.3,aj-2取值为0.2。
步骤二、估计初始透射率图
S201、将图像Ik各元素的值除以步骤一得到的经过学习的当前帧大气光值Ak,得到图像IA,将图像IA划分为半径为r的窗口,求每个窗口中的各个像素点的R、G和B三通道中的最小通道值,再从该窗口中各像素点的最小通道值中选择一个最小值,得到每个窗口对应的暗通道图I′A;将该暗通道图I′A与加权系数ω的乘积后,再用1减去该乘积,最终估计得到初始的透射率
t ~ ( p ) = 1 - ω min p ∈ Ω ( min c ∈ { r , g , b } I c ( p ) A c ) - - - ( 3 )
其中,半径r的取值范围为7-15个像素,本实施例中r取15。
Ik与A0为三通道的矩阵和向量,ω为加权系数,取值0.95。根据测试结果,ω取值范围可在0.90至0.99之间,优选地,ω取值范围在0.94至0.96之间。
S202、判断:如果视频中的当前帧为第一帧,则将透射率作为经过调整得到的当前帧的初始透射率;
如果视频中当前帧非第一帧,则通过时间相干性方法获取调整后的初始透射率图,所述时间相干性方法求解透射率图的方法进一步包括:
将每一帧的RGB彩色图像转换到YUV空间,得到Y通道的图像Iy。将步骤一获得的经过学习的大气光值Ak转换为标量:
A=Ar·0.3+Ag·0.59+Ab·0.11(4)
其中Ar、Ag和Ab分别表示大气光值Ak的三通道分量值;
S203、计算图像Ik中任意像素点P在相邻帧出现在同一个位置的概率,一种简单的方法,即下式中的权重值wk(p):
w k ( p ) = exp ( - ( I k ( p ) - I k - 1 ( p ) ) 2 σ 2 ) - - - ( 5 )
其中,Ik(p)和Ik-1(p)表示同一像素点P在第k帧和k-1帧所具有的Y通道值,σ2表示模型的方差。所述σ2的取值为10。
S204、第k帧和k-1帧的透射率在像素点P的关系为tk(p)=τk(p)tk-1(p),则时间相干性因子根据下式求出:
τ k ( p ) = I k ( p ) - A I k - 1 ( p ) - A - - - ( 6 )
S205、上述计算均基于像素点,由于图像相邻像素点存在物理联系,下面计算基于块的时间相干性因子,即:对于图像Ik中的任意像素点P,以像素点P为中心,r’为半径划分区域块Ω,r’取值为5至15个像素;采用下式得到基于块的相干因子
τ k ‾ = Σ p ∈ Ω w k ( p ) τ k ( p ) Σ p ∈ Ω w k ( p ) - - - ( 7 )
S206、根据基于块的相干因子得到后构造代价函数,如下式所述:
E = Σ p ∈ Ω ( τ k ‾ t k - 1 ( p ) - t k ) 2 σ t 2 + Σ q ∈ N p ( t ~ ( q ) - t k ) 2 max ( | | I k ( p ) - I k - 1 ( q ) | | 2 , a 2 ) - - - ( 8 )
其中,tk为最终要求的透射率,为未知数;σt 2表示透射率估计的模型误差,所述σt 2取值为10。等号右侧第一项表示时间相干性,第二项表示空间相干项,用于描述邻域相干性,Np表示像素点P的邻域像素点。其中,tk-1(p)表示第k-1帧中像素点P的透射率;Ik-1(q)表示第k-1帧图像中像素点p的邻域像素点的强度;表示像素点p的邻域像素点的透射率,由暗通道求得的初始透射率;a2用来限制‖Ik(p)-Ik-1(q)‖2值的边界,所述a2取值为0.0052
S207、求式(8)表示的代价E最小时的透射率tk,即可使得透射率图具有更好的空间和时间平滑性,具体步骤包括:
设在每一个区域块中透射tk初始值为0.1,代入代价函数求代价值E1;然后将透射率tk逐次增加0.1,每次均求出Ei(1≤i≤10);选出最小值,此时对应的tk即为满足要求的最小能量区域块透射率值。
每一个区域块均重复上述步骤,求出每一个区域块对应的最佳的tk。最终得到经过调整的整个图像透射率图。
步骤三、将输入图像的YUV空间的Y通道图像作为导向图像,初始透射率图tk作为输入图,对初始透射率tk导向滤波,得到各个区域块最终的透射率t*
步骤四、针对输入的当前帧图像的各个像素点,通过以下公式得到去雾处理后的图像:
J ( p ) = I ( p ) - A k t * ( p ) + A k - - - ( 9 )
其中,t*(p)表示任意像素点p所在位置的最终的透射率,根据步骤3的最终的透射率t*获得;I(p)为输入的RGB彩色图像的像素点p的三通道值,Ak为步骤一得到的三通道大气光值。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种雾霾条件下视频图像增强方法,其特征在于,包括如下步骤:
步骤一、求取基于时间连续自适应调整的大气光值,具体为:
A、针对输入视频的当前帧图像Ik,获得每一个像素点的三个独立的单通道图像R、G、B中的最小值,遍历图像Ik的各个像素点,得到暗通道图Idark
将暗通道图Idark中像素亮度值从大到小排列的前0.1%的像素点筛选出来,其它像素点不考虑,形成大气光区域图;针对图像Ik中与大气光区域图对应位置的各个像素点,获得各个像素点通道R、通道G以及通道B的平均值,则得到的三个平均值作为大气光值,该大气光值为一个有R,G,B分量的向量。
B、对视频中的每一帧图像进行如a中所述的处理,获得各帧图像原始大气光值;
C、如果视频中当前帧为第一帧,则将该帧图像的原始大气光值作为经过学习得到的当前帧大气光值Ak,此时k=0;
D、如果视频中当前帧非第一帧,则当前帧的前N帧大气光值以及当前帧原始大气光值进行加权平均,作为经学习得到的当前帧大气光值Ak;其中,所述N至少取2,且小于或等于k-1;
步骤二、估计初始透射率图,具体包括如下步骤:
S201、将图像Ik各元素的值除以步骤一得到的经过学习的当前帧大气光值Ak,得到图像IA,将图像IA划分为半径为r的窗口,求每个窗口中的各个像素点的R、G和B三通道中的最小通道值,再从该窗口中各像素点的最小通道值中选择一个最小值,得到每个窗口对应的暗通道图I′A;将该暗通道图I′A与加权系数ω的乘积后,再用1减去该乘积,最终估计得到初始的透射率
其中,半径r的取值范围为7-15个像素;加权系数ω取值范围在0.90至0.99之间;
S202、判断:如果视频中的当前帧为第一帧,则将透射率作为经过调整得到的当前帧的初始透射率;
如果视频中当前帧非第一帧,则通过时间相干性方法获取调整后的初始透射率图,所述时间相干性方法求解透射率图的方法进一步包括:
将每一帧的RGB彩色图像转换到YUV空间,得到Y通道的图像Iy;将步骤一获得的经过学习的大气光值Ak转换为标量:
A=Ar·0.3+Ag·0.59+Ab·0.11(4)
其中Ar、Ag和Ab分别表示大气光值Ak的三通道分量值;
S203、第k帧透射率tk(p)和k-1帧的透射率tk-1(p)在像素点P的关系为tk(p)=τk(p)tk-1(p),则时间相干性因子根据下式求出:
其中,Ik(p)和Ik-1(p)表示同一像素点P在第k帧和k-1帧Y通道的图像所具有单通道值,σ2表示模型的方差;所述σ2的取值为10;
S204、根据时间相干因子τk,得到tk(p)=τk(p)tk-1(p)后构造代价函数,如下式所述:
其中,tk为最终要求的透射率,为未知数;σt 2表示透射率估计的模型误差,所述σt 2取值为10;Np表示像素点P的邻域像素点;tk-1(p)表示第k-1帧中像素点P的透射率;Ik-1(q)表示第k-1帧图像中像素点p的邻域像素点的Y通道图像的值;表示像素点p的邻域像素点的透射率,由暗通道求得的初始透射率;所述a2取值为0.0052
S205、求式(8)表示的代价E最小时的透射率tk
步骤三、将输入图像的YUV空间的Y通道图像作为导向图像,初始透射率图tk作为输入图,对初始透射率tk导向滤波,得到各个区域块最终的透射率t*
步骤四、针对输入的当前帧图像的各个像素点,通过以下公式得到去雾处理后的图像:
其中,t*(p)表示任意像素点p所在位置的最终的透射率,根据步骤3的最终的透射率t*获得;I(p)为输入的RGB彩色图像的像素点p的三通道值,Ak为步骤一得到的三通道大气光值。
2.如权利要求1所述的一种雾霾条件下视频图像增强方法,其特征在于,在所述步骤中,当前帧的前N帧大气光值以及当前帧原始大气光值进行加权平均的具体方法如下:
找到大气光区域图中亮度最小的像素点,求该像素点3×3窗口内对应图像Ik中的像素值的差值Bstd和平均值Bmean,再求该差值Bstd与平均值Bmean之间的差值:
L=|Bstd-Bmean|(1);
最后,判断当前第k帧对应的差值Lk与第k-1帧的差值Lk-1之间的差值是否小于或等于预设的阈值Lthresh,若第k帧Lk与第k-1帧Lk-1差值小于或等于Lthresh,则将当前帧原始大气光值作为经过学习得到的当前帧大气光值Ak;否则将当前帧的前N帧大气光值以及当前帧原始大气光值进行加权平均,作为经学习得到的当前帧大气光值Ak;其中,所述Lthresh取值为0.02;
所述经学习得到的当前帧大气光值Ak的R、G和B任意一个通道的值为:
其中,分别表示原始第k帧和第j帧大气光,Ak表示经过学习获得的第k帧大气光,aj表示当前帧更新权重;所述权重aj的取值原则为,与当前帧越近,取值越大,反之越小。
3.如权利要求1所述的一种雾霾条件下视频图像增强方法,其特征在于,所述r取15。
4.如权利要求1所述的一种雾霾条件下视频图像增强方法,其特征在于,ω取值范围在0.94至0.96之间。
5.如权利要求4所述的一种雾霾条件下视频图像增强方法,其特征在于,ω取值0.95。
6.如权利要求1所述的一种雾霾条件下视频图像增强方法,其特征在于,所述步骤二中,得到时间相干性因子后,再得到基于区域块的时间相干因子,具体为:
首先,计算图像Ik中任意像素点P在相邻帧出现在同一个位置的概率,一种简单的方法,即下式中的权重值wk(p):
其次,对于图像Ik中的任意像素点P,以像素点P为中心,r’为半径划分区域块Ω,r’取值为5至15个像素;采用下式得到基于块的时间相干因子
最后,将该基于块的时间相干因子代替时间相干因子τk代入到代价E的表达式中,最终求得代价E最小时的透射率tk
7.如权利要求1所述的一种雾霾条件下视频图像增强方法,其特征在于,所述步骤二中,求代价E最小时的透射率tk的方法具体包括:
设在每一个区域块中透射tk初始值为0.1,代入代价函数求代价值E1;然后将透射率tk逐次增加0.1,每次均求出对应的代价值;选出最小代价值,此时对应的tk即为满足要求的最小能量区域块透射率值;其中,累加次数小于10次;
每一个区域块均重复上述步骤,求出每一个区域块对应的最佳的tk,最终得到经过调整的整个图像透射率图。
CN201610143521.1A 2016-03-14 2016-03-14 一种雾霾条件下视频图像增强方法 Active CN105631831B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610143521.1A CN105631831B (zh) 2016-03-14 2016-03-14 一种雾霾条件下视频图像增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610143521.1A CN105631831B (zh) 2016-03-14 2016-03-14 一种雾霾条件下视频图像增强方法

Publications (2)

Publication Number Publication Date
CN105631831A true CN105631831A (zh) 2016-06-01
CN105631831B CN105631831B (zh) 2018-07-27

Family

ID=56046719

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610143521.1A Active CN105631831B (zh) 2016-03-14 2016-03-14 一种雾霾条件下视频图像增强方法

Country Status (1)

Country Link
CN (1) CN105631831B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106169176A (zh) * 2016-06-27 2016-11-30 上海集成电路研发中心有限公司 一种图像去雾方法
CN107133927A (zh) * 2017-04-21 2017-09-05 汪云飞 超像素框架下基于均值‑均方差暗通道的单幅图像去雾方法
CN107240075A (zh) * 2017-05-27 2017-10-10 上海斐讯数据通信技术有限公司 一种雾霾图像增强处理方法及系统
CN107680055A (zh) * 2017-09-26 2018-02-09 成都国翼电子技术有限公司 一种基于人工交互的航拍图像雾霾去除方法
CN109118441A (zh) * 2018-07-17 2019-01-01 厦门理工学院 一种低照度图像及视频增强方法、计算机装置及存储介质
CN109636735A (zh) * 2018-11-02 2019-04-16 中国航空工业集团公司洛阳电光设备研究所 一种基于时空一致性约束的快速视频去雾方法
CN109993714A (zh) * 2019-04-10 2019-07-09 南京航空航天大学 一种有色眩光的去除方法
CN110278346A (zh) * 2018-03-15 2019-09-24 财团法人交大思源基金会 视频除雾处理装置及方法
CN110544220A (zh) * 2019-09-05 2019-12-06 北京天地玛珂电液控制系统有限公司 一种煤矿井下视频图像智能去雾方法、存储介质及系统
CN110910319A (zh) * 2019-10-30 2020-03-24 中国医学科学院北京协和医院 一种基于大气散射模型的手术视频实时去雾增强方法
CN111192210A (zh) * 2019-12-23 2020-05-22 杭州当虹科技股份有限公司 一种自适应增强的视频去雾方法
CN111738928A (zh) * 2020-04-30 2020-10-02 南京图格医疗科技有限公司 一种基于概率优化和神经网络的内窥镜去雾方法及装置
CN114331920A (zh) * 2022-03-09 2022-04-12 浙江大华技术股份有限公司 图像处理方法及装置、存储介质、电子装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140017776A (ko) * 2012-08-01 2014-02-12 엠텍비젼 주식회사 영상 처리 장치 및 영상 내의 안개 제거 방법
CN105354806A (zh) * 2015-11-20 2016-02-24 上海熙菱信息技术有限公司 基于暗通道的快速去雾方法及系统
CN105374013A (zh) * 2014-08-15 2016-03-02 台北科技大学 基于双暗通道先验的图像能见度修复方法及图像处理装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140017776A (ko) * 2012-08-01 2014-02-12 엠텍비젼 주식회사 영상 처리 장치 및 영상 내의 안개 제거 방법
CN105374013A (zh) * 2014-08-15 2016-03-02 台北科技大学 基于双暗通道先验的图像能见度修复方法及图像处理装置
CN105354806A (zh) * 2015-11-20 2016-02-24 上海熙菱信息技术有限公司 基于暗通道的快速去雾方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KRISTOFOR B.GIBSON等: "An Investigation of Dehazing Effects on Image and Video Coding", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 *
ZHAO DONG等: "A Novel Fast Haze Removal Technique for Single Image Using Image Pyarmid", 《IMAGE PROCESSING,DEEP LEARNING VIEW PROJECT》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106169176A (zh) * 2016-06-27 2016-11-30 上海集成电路研发中心有限公司 一种图像去雾方法
CN107133927A (zh) * 2017-04-21 2017-09-05 汪云飞 超像素框架下基于均值‑均方差暗通道的单幅图像去雾方法
CN107133927B (zh) * 2017-04-21 2020-03-17 汪云飞 超像素框架下基于均值-均方差暗通道的单幅图像去雾方法
CN107240075A (zh) * 2017-05-27 2017-10-10 上海斐讯数据通信技术有限公司 一种雾霾图像增强处理方法及系统
CN107680055A (zh) * 2017-09-26 2018-02-09 成都国翼电子技术有限公司 一种基于人工交互的航拍图像雾霾去除方法
CN110278346A (zh) * 2018-03-15 2019-09-24 财团法人交大思源基金会 视频除雾处理装置及方法
CN109118441B (zh) * 2018-07-17 2022-04-12 厦门理工学院 一种低照度图像及视频增强方法、计算机装置及存储介质
CN109118441A (zh) * 2018-07-17 2019-01-01 厦门理工学院 一种低照度图像及视频增强方法、计算机装置及存储介质
CN109636735A (zh) * 2018-11-02 2019-04-16 中国航空工业集团公司洛阳电光设备研究所 一种基于时空一致性约束的快速视频去雾方法
CN109636735B (zh) * 2018-11-02 2023-03-10 中国航空工业集团公司洛阳电光设备研究所 一种基于时空一致性约束的快速视频去雾方法
CN109993714A (zh) * 2019-04-10 2019-07-09 南京航空航天大学 一种有色眩光的去除方法
CN110544220A (zh) * 2019-09-05 2019-12-06 北京天地玛珂电液控制系统有限公司 一种煤矿井下视频图像智能去雾方法、存储介质及系统
CN110544220B (zh) * 2019-09-05 2022-02-22 北京天玛智控科技股份有限公司 一种煤矿井下视频图像智能去雾方法、存储介质及系统
CN110910319B (zh) * 2019-10-30 2022-10-21 中国医学科学院北京协和医院 一种基于大气散射模型的手术视频实时去雾增强方法
CN110910319A (zh) * 2019-10-30 2020-03-24 中国医学科学院北京协和医院 一种基于大气散射模型的手术视频实时去雾增强方法
CN111192210A (zh) * 2019-12-23 2020-05-22 杭州当虹科技股份有限公司 一种自适应增强的视频去雾方法
CN111738928A (zh) * 2020-04-30 2020-10-02 南京图格医疗科技有限公司 一种基于概率优化和神经网络的内窥镜去雾方法及装置
CN111738928B (zh) * 2020-04-30 2021-04-13 南京图格医疗科技有限公司 一种基于概率优化和神经网络的内窥镜去雾方法
CN114331920A (zh) * 2022-03-09 2022-04-12 浙江大华技术股份有限公司 图像处理方法及装置、存储介质、电子装置

Also Published As

Publication number Publication date
CN105631831B (zh) 2018-07-27

Similar Documents

Publication Publication Date Title
CN105631831B (zh) 一种雾霾条件下视频图像增强方法
Mehra et al. ReViewNet: A fast and resource optimized network for enabling safe autonomous driving in hazy weather conditions
CN107301623B (zh) 一种基于暗通道和图像分割的交通图像去雾方法及系统
CN102750674B (zh) 基于自适应容差的视频图像去雾方法
CN102831591B (zh) 一种基于高斯滤波的单幅图像的实时去雾方法
CN106780356B (zh) 基于卷积神经网络和先验信息的图像去雾方法
CN109410129A (zh) 一种低光照图像场景理解的方法
CN104036466B (zh) 一种视频去雾方法与系统
CN102202163B (zh) 一种监控视频的自适应增强方法及其装置
CN103258332B (zh) 一种抗光照变化的运动目标的检测方法
CN102768760B (zh) 一种基于图像纹理的图像快速去雾方法
CN111667433B (zh) 一种基于简单线性迭代聚类优化的无人机图像去雾方法
CN107451966B (zh) 一种采用灰度图引导滤波实现的实时视频去雾方法
TWI808406B (zh) 圖像去霧方法和使用圖像去霧方法的圖像去霧設備
CN107301624A (zh) 基于区域划分及浓雾预处理的卷积神经网络去雾算法
CN105898111B (zh) 一种基于谱聚类的视频去雾方法
CN106657948A (zh) 低照度Bayer图像的增强方法及增强装置
CN114549297A (zh) 一种基于不确定分析的无监督单目深度估计方法
CN106023108A (zh) 基于边界约束和上下文正则化的图像去雾算法
CN111598814B (zh) 基于极端散射通道的单图像去雾方法
Liang et al. Learning to remove sandstorm for image enhancement
CN113379619B (zh) 去雾成像、能见度提取和景深估计的一体化处理方法
CN110544216A (zh) 基于深度学习的视频去雾系统
CN112215766B (zh) 一种图像复原与图像增强相融合的图像去雾方法及其卷积网络
CN106846260B (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