本文介绍的是一种比较实用的低照度图像增强效果很好的方法,Xuan Dong论文《Fast efficient algorithm for enhancement of low lighting video》

概述

低光照图像取反(255 - 低光照图像),和有雾图像存在一些共性,比如在天空或者遥远的背景这些地方,亮度值都是很高的,但在近景的 RGB 三通道中至少有一个亮度值很低。作者提出直接用去雾算法对低光照图像的反转图像处理,去雾结果再取反得到低照度增强结果。

低照度图像增强算法

大气散射模型

大多透雾算法基于大气散射模型,模型如下:

I

(

x

)

=

J

(

x

)

t

(

x

)

+

A

(

1

t

(

x

)

)

I (x) = J (x)t(x) + A(1 − t(x))

I(x)=J(x)t(x)+A(1−t(x))

其中,

I

(

x

)

I(x)

I(x)表示含雾图,

J

(

x

)

J(x)

J(x)表示无雾图,

A

A

A表示大气光模型,

t

(

x

)

=

e

β

d

(

x

)

t(x)=e^{- \beta d(x)}

t(x)=e−βd(x)表示传输率矩阵,

β

\beta

β 表示大气散射系数,

d

(

x

)

d(x)

d(x)表示图像区域的景深。雾越浓,到达相机的物体反射的光越少,即传输率越小。

​上式可以改写为:

J

(

x

)

=

I

(

x

)

A

t

(

x

)

+

A

J (x) = \frac{I (x) - A}{t(x)}+A

J(x)=t(x)I(x)−A​+A

为了得到

J

(

x

)

J(x)

J(x),重点在估计

A

A

A和

t

(

x

)

t(x)

t(x)。

暗通道先验知识:

在大部分无雾图像的无天空区域,像素中至少存在一个颜色通道存在很低非常低的亮度值。这个最低的亮度值几乎等同于0。因此,对于一个观测图像J,其暗通道趋近于0,即

J

d

a

r

k

(

X

)

=

m

i

n

y

Ω

(

x

)

(

m

i

n

c

r

,

g

,

b

J

c

(

y

)

)

0

J^{dark}(X) = min_{y \in \Omega (x)}(min_{c \in {r,g,b }} J^{c}(y)) \to 0

Jdark(X)=miny∈Ω(x)​(minc∈r,g,b​Jc(y))→0

其中,

J

c

J^{c}

Jc表示彩色图像每个通道,

Ω

(

x

)

\Omega(x)

Ω(x)表示以

x

x

x为中心的窗口

透雾算法

基于大气散射模型,透雾模型的方法步骤如下:

1、从雾图I (x) 估计传输率矩阵t(x)

2、估计大气光值A

3、通过公式估计无雾图J(x)

首先求出每个像素RGB分量中的最小值,存入一副和原始图像大小相同的灰度图中,然后再对这幅灰度图进行最小值滤波,滤波的半径由窗口大小决定,一般有

W

i

n

d

o

w

S

i

z

e

=

2

R

a

d

i

u

s

+

1

WindowSize = 2 * Radius + 1

WindowSize=2∗Radius+1

将雾模型

I

(

x

)

=

J

(

x

)

t

(

x

)

+

A

(

1

t

(

x

)

)

I (x) = J (x)t(x) + A(1 − t(x))

I(x)=J(x)t(x)+A(1−t(x))

处理为:

I

c

(

x

)

A

c

=

J

c

(

x

)

A

c

t

(

x

)

+

1

t

(

x

)

\frac{I^{c} (x)}{A^{c}} = \frac{J^{c} (x)}{A^{c}}t(x) + 1 − t(x)

AcIc(x)​=AcJc(x)​t(x)+1−t(x)

假设在窗口内透射率$ t(x)

为常数,定义为

为常数,定义为

为常数,定义为\tilde t(x)

,A$值已知。对上式求两次最小值运算,得到

m

i

n

y

Ω

(

x

)

(

m

i

n

c

I

c

(

y

)

A

c

)

=

t

~

(

x

)

m

i

n

y

Ω

(

x

)

(

m

i

n

c

J

c

(

y

)

A

c

)

+

1

t

~

(

x

)

\underset {y \in \Omega (x)}{min}(\underset{c}{min} \frac{I^{c} (y)}{A^{c}}) = \tilde t(x) \underset {y \in \Omega (x)}{min}(\underset{c}{min} \frac{J^{c} (y)}{A^{c}})+ 1 − \tilde t(x)

y∈Ω(x)min​(cmin​AcIc(y)​)=t~(x)y∈Ω(x)min​(cmin​AcJc(y)​)+1−t~(x)

根据暗原色先验理论

J

d

a

r

k

(

x

)

=

m

i

n

y

Ω

(

x

)

(

m

i

n

c

J

c

(

y

)

)

=

0

J^{dark}(x) = \underset {y \in \Omega (x)}{min}(\underset{c }{min} J^{c}(y)) = 0

Jdark(x)=y∈Ω(x)min​(cmin​Jc(y))=0

推导出

m

i

n

y

Ω

(

x

)

(

m

i

n

c

J

c

(

y

)

A

c

)

=

0

\underset {y \in \Omega (x)}{min}(\underset{c }{min} \frac{J^{c}(y)}{A^{c}}) = 0

y∈Ω(x)min​(cmin​AcJc(y)​)=0

带入透射率为常数的公式,得到透射率预估值

t

~

(

x

)

=

1

m

i

n

y

Ω

(

x

)

(

m

i

n

c

I

c

(

y

)

A

c

)

\tilde t(x) = 1-\underset {y \in \Omega (x)}{min}(\underset{c }{min} \frac{I^{c}(y)}{A^{c}})

t~(x)=1−y∈Ω(x)min​(cmin​AcIc(y)​)

上式添加一个限制系数,得到

t

~

(

x

)

\tilde t(x)

t~(x):

t

~

(

x

)

=

1

ω

m

i

n

y

Ω

(

x

)

(

m

i

n

c

I

c

(

y

)

A

c

)

\tilde t(x) = 1-\omega\underset {y \in \Omega (x)}{min}(\underset{c }{min} \frac{I^{c}(y)}{A^{c}})

t~(x)=1−ωy∈Ω(x)min​(cmin​AcIc(y)​)

ω

\omega

ω取值0.95

以上假设全球达气光A值时已知的,在实际中,我们可以借助于暗通道图来从有雾图像中获取该值。具体步骤如下:

1)从暗通道图中按照亮度的大小取前0.1%的像素。

​2)在这些位置中,在原始有雾图像I中寻找对应的具有最高亮度的点的值,作为

A

A

A值。

算法参数优化

文中自适应调节

t

(

x

)

t(x)

t(x),如下式

t

(

x

)

=

{

2

t

(

x

)

,

0

<

t

(

x

)

<

0.5

1

,

0.5

<

t

(

x

)

<

1

t'(x) = \begin{cases} 2*t(x),0<t(x)<0.5\\1,0.5<t(x)<1 \end{cases}

t′(x)={2∗t(x),0<t(x)<0.51,0.5<t(x)<1​

优化

t

(

x

)

t(x)

t(x)如下,增强暗区提亮,亮区微处理,效果更好

t

(

x

)

=

{

t

(

x

)

2

2

,

0

<

t

(

x

)

<

0.5

t

(

x

)

,

0.5

<

t

(

x

)

<

1

t'(x) = \begin{cases} t(x)^2*2,0<t(x)<0.5\\t(x),0.5<t(x)<1 \end{cases}

t′(x)={t(x)2∗2,0<t(x)<0.5t(x),0.5<t(x)<1​

效果对比




图、测试图像1、原文处理效果、算法优化处理效果




图、测试图像2、原文处理效果、算法优化处理效果




图、测试图像3、原文处理效果、算法优化处理效果 ## 算法总结 简单易实现的增强算法,普世性较好,透雾算法处理部分兼容各种透雾算法。

算法实现参考

https://github.com/AomanHao/ISP_Low_Light_Image_Enhancement

我的个人博客主页,欢迎访问

我的CSDN主页,欢迎访问

我的GitHub主页,欢迎访问

我的个人博客主页,欢迎访问

我的CSDN主页,欢迎访问

我的GitHub主页,欢迎访问

【阅读笔记】低照度图像增强-《Fast efficient algorithm for enhancement of low lighting video》的更多相关文章

  1. Mysql DOC阅读笔记

    Mysql DOC阅读笔记 转自我的Github Speed of SELECT Statements 合理利用索引 隔离调试查询中花费高的部分,例如函数调用是在结果集中的行执行还是全表中的行执行 最 ...

  2. [置顶] 人工智能(深度学习)加速芯片论文阅读笔记 (已添加ISSCC17,FPGA17...ISCA17...)

    这是一个导读,可以快速找到我记录的关于人工智能(深度学习)加速芯片论文阅读笔记. ISSCC 2017 Session14 Deep Learning Processors: ISSCC 2017关于 ...

  3. 阅读笔记 1 火球 UML大战需求分析

    伴随着七天国庆的结束,紧张的学习生活也开始了,首先声明,阅读笔记随着我不断地阅读进度会慢慢更新,而不是一次性的写完,所以会重复的编辑.对于我选的这本   <火球 UML大战需求分析>,首先 ...

  4. [阅读笔记]Software optimization resources

    http://www.agner.org/optimize/#manuals 阅读笔记Optimizing software in C++   7. The efficiency of differe ...

  5. 《uml大战需求分析》阅读笔记05

    <uml大战需求分析>阅读笔记05 这次我主要阅读了这本书的第九十章,通过看这章的知识了解了不少的知识开发某系统的重要前提是:这个系统有谁在用?这些人通过这个系统能做什么事? 一般搞清楚这 ...

  6. <<UML大战需求分析>>阅读笔记(2)

    <<UML大战需求分析>>阅读笔记(2)> 此次读了uml大战需求分析的第三四章,我发现这本书讲的特别的好,由于这学期正在学习设计模式这本书,这本书就讲究对uml图的利用 ...

  7. uml大战需求分析阅读笔记01

    <<UML大战需求分析>>阅读笔记(1) 刚读了uml大战需求分析的第一二章,读了这些内容之后,令我深有感触.以前学习uml这门课的时候,并没有好好学,那时我认为这门课并没有什 ...

  8. Hadoop阅读笔记(七)——代理模式

    关于Hadoop已经小记了六篇,<Hadoop实战>也已经翻完7章.仔细想想,这么好的一个框架,不能只是流于应用层面,跑跑数据排序.单表链接等,想得其精髓,还需深入内部. 按照<Ha ...

  9. Hadoop阅读笔记(六)——洞悉Hadoop序列化机制Writable

    酒,是个好东西,前提要适量.今天参加了公司的年会,主题就是吃.喝.吹,除了那些天生话唠外,大部分人需要加点酒来作催化剂,让一个平时沉默寡言的码农也能成为一个喷子!在大家推杯换盏之际,难免一些画面浮现脑 ...

  10. Hadoop阅读笔记(五)——重返Hadoop目录结构

    常言道:男人是视觉动物.我觉得不完全对,我的理解是范围再扩大点,不管男人女人都是视觉动物.某些场合(比如面试.初次见面等),别人没有那么多的闲暇时间听你诉说过往以塑立一个关于你的完整模型.所以,第一眼 ...

随机推荐

  1. LINUX 下LAMP之yum方法部署环境

    一. 系统组成 系统由前台网站和后台保种服务两部分组成,其中后台保种可以扩展成多台设备,提高保种的资源量以及性能.下面安装以houtaiip代表后台保重服务器ip地址,qiantaiip代表前台网站i ...

  2. 链式描述线性表(C++实现)

    在链式描述中,线性表的元素在内存中的存储位置是随机的,每个元素都有一个明确的指针或链指向下一个元素的位置 chain类 在此使用单向链表实现了线性表,其中最后一个节点的指针域为NULL,即它用单向链接 ...

  3. 1 分钟给 Siri 升个级!从智Z变身 ChatSiri!

    原文链接:https://forum.laf.run/d/79/17 众所周知,Siri 是一个智 Z!那么如果能接入大火的 chatGPT,是不是就会从智 Z 变成人工智能?! 众所周知,Laf 是 ...

  4. Linux 根据名称自动kill掉当前相关进程

    ps aux | grep app | grep -v "grep" | awk '{print $2}' | xargs -r kill

  5. 【介绍】.NET新加特性介绍

    ​ 简介 当下的.Net新版本引进了几种新特性,包括全局命名空间引用.可空引用类型和顶级语句.这些特性在一定程度上改善了 .NET 平台的开发效率, 对于短小精干的小程序,这些新的特性无疑可以把开发效 ...

  6. 使用扩展函数方式,在Winform界面中快捷的绑定树形列表TreeList控件和TreeListLookUpEdit控件

    在一些字典绑定中,往往为了方便展示详细数据,需要把一些结构树展现在树列表TreeList控件中或者下拉列表的树形控件TreeListLookUpEdit控件中,为了快速的处理数据的绑定操作,比较每次使 ...

  7. C# 打开蓝牙设置界面

    蓝牙设置相关界面,以下是通过C#方式打开的几个方式,记录一下 蓝牙设置界面 1.控制面板命令bthprops.cpl 可以用控制面板 control+bthprops.cpl,也可以直接bthprop ...

  8. 在循环内调用 size() 方法的开销大吗?

    for (int i = 0; i < buildings.size(); i++) {} 和 int n = buildings.size(); for (int i = 0; i < ...

  9. Java 新的生态型应用开发框架,Solon v2.2.14 发布

    Java 新的生态型应用开发框架,Solon :更快.更小.更简单.从零开始构建,有自己的标准规范与开放生态: 150多个生态插件,可以满足各种场景开发 大量的国产框架适配,可以为应用软件国产化提供更 ...

  10. 2021-11-15:四数相加 II。给你四个整数数组 nums1、nums2、nums3 和 nums4 ,数组长度都是 n ,请你计算有多少个元组 (i, j, k, l) 能满足:0 <= i,

    2021-11-15:四数相加 II.给你四个整数数组 nums1.nums2.nums3 和 nums4 ,数组长度都是 n ,请你计算有多少个元组 (i, j, k, l) 能满足:0 <= ...