Retinex 算法
目录
1、Retinex 理论及数学模型
2、Retinex 算法发展历程
2.1 基于迭代的 Retinex 算法
2.1.1 Frankle-McCann Retinex 算法
2.1.2 McCann99 Retinex 算法
2.2 基于中心环绕的 Retinex 算法
2.2.1 SSR算法(单尺度)
2.2.2 MSR算法(多尺度)
2.2.3 MSRCR算法
2.3 Retinex算法的频域处理
1、Retinex 理论及数学模型
Retinex 理论指出,入射光决定了一副图像中的所有像素点的动态范围的大小,而物体自身所固有不变的反射系数决定了图像的内在固有属性。也就是说,我们所看到的图像是照射光根据物体的反射系数所反射的光线形成的,如图所示。很显然,如果把一幅图像看作是由照射光和反射光组成的话,Retinex 图像增强的基本思想就是去除照射光的影响,保留物体自身的反射属性。
核心:从S图像中估测L分量,并去除L分量,得到原始反射分量R
retina(视网膜) + cortex(皮层) —— Retinex
提出者:埃德温·赫伯特·兰德 —— Polaroid(宝丽来)
Retinex百度百科RETINEX_百度百科
2、Retinex 算法发展历程
2.1 基于迭代的 Retinex 算法
2.1.1 Frankle-McCann Retinex 算法
基本原理:
Frankle-McCann 算法采用了一种新的基于螺旋结构的迭代分段线性比较路径,螺旋结构路径像素点间的比较是一个由远到近的比较过程,在进行完一次比较之后,下一次的做比较的两个像素点间的间距缩短为上一次比较间距的一半,并且比较路径的方向同时也按顺时针方向发生转变,就这样逐次比较直至像素点间距为 1 为止。
Frankle-McCann Retinex 算法步骤 :
(1)数据的前期转换
把原图像的像素值由整数域转到对数域,减少后续算法的运算量
对于彩色图像要先分解为R.G.B三幅灰度图再转换
因为RGB图像的像素值是(0,255),在做对数运算时为避免负值的出现,可以将原图像像素值整体加 1,即s (x, y) = log(1+S (x, y))
(2)初始化一个与原图像 S(x,y)同样大小的常数为 t 的矩阵
若m*n的图像,则矩阵包含m*n个t,保证每个像素点都进行一次迭代。
t是原图像亮度的均值,
(3)求解 S = 2P的值,P=fix[log2 min(m, n)-1]
S是目标点与两个比较点之间的最大距离
(4)计算路径上的像素点
假设 rn (x,y)是上一次迭代的结果,将此次迭代差值累加保存到相应的 rn (x,y)位置中,最终得到此次的迭代结果,然后再对两者做一个平均,最后得到输出结果。△l是目标点在此路径上的亮度差。
(5)令S=-S/2
每下一步的两个比较点与目标点的间距缩短为上一步的一半,同时方向按顺时针改变,即S=-S/2
(6)重复3.4.5三个步骤直到|S|<1
(7)迭代n次,也就是重复3.4.5.6四个步骤n次
每次迭代选取不同的初始比较点
(8)线性拉伸,彩色图像还需要三个通道的合成,然后输出显示
经过像素间的比较校正 n 次迭代之后,输出结果是以初始化值 t为中心,集中分布在 t 附近的一系列的浮点数。也就是说原图像的数据经过迭代后起到了压缩的效果,因此需要对迭代结果做线性拉伸处理,提高图像对比度。通常采用的 8 bit 图像的动态范围值是 0 到 255,拉伸公式:
其中max和min分别是迭代结果的最大和最小值,R(x,y)是最终的增强结果。
2.1.2 McCann99 Retinex 算法
基本原理:
提出了一种金字塔比较模型,由金字塔顶端到底端分辨率由低到高,依次逐层迭代。这种采用分辨率由低到高的迭代方法可以有效减少算法的运算量,节约算法的执行时间。
算法步骤:
1) 将原图像变换到对数域S;
2) 初始化(计算图像金字塔层数;初始化常数图像矩阵R作为进行迭代运算的初始值);
3) 从顶层开始,到最后一层进行8邻域比较运算,运算规则与MccCann Retinex算法相同;
4) 第n层运算结束后对第n层的运算结果进行插值,变成原来的两倍,与n+1层大小相同;
5) 当最底层计算完毕得到的即最终增强后的图像。
局限性:
读取图像的尺寸必须符合2的整数次幂,即图像长宽可表示为(col*2^n)*(row*2^n)
Col必须大于row ,并且属于集合{1,2,3,4,5}
改进算法:
适用于各种尺寸的图像,保持扩展后的图像与原图的一致性,尽量避免扩展像素点对原图的干扰,采取只对图像边界扩展的方法。
基于迭代的 Retinex 算法总结:
迭代次数对算法的影响:
算法的执行时间随着迭代次数 n 取值变大而增加
n 取 5 到 8 之间能同时兼顾图像质量和计算速度
优点:
颜色恒常性、动态范围压缩大(像素点丰富)、色彩逼真度高(图像高保真)
缺点:
增强后的图像存在“光晕”现象,即在图像色彩交界处渐变
光晕缺陷同样也存在于其他 Retinex算法之中
2.2 基于中心环绕的 Retinex 算法
基于路径比较的 Retinex 算法和 McCann 算法,像素点之间亮暗关系的比较路径实质上都是一维的,一维路径的比较方法不论在效率上还是在准确性上都有较大的缺陷。为了更加准确地估计出图像的照射分量,1986 年Land 在分析一维路径缺陷的基础上对算法做了改进,用二维区域内像素比较来取代一维路径比较,提出了基于二维路径比较的 Retinex 算法,经过后续发展之后就是我们现在所熟知的基于中心环绕的 Retinex 算法。
2.2.1 SSR算法(单尺度)
基本理论:
在计算图像中目标像素点的灰度值时,是通过以目标点为中心的区域内的像素值的加权得到的,权重的比例大小则是由环绕函数(高斯函数)来确定。
σ 越小,高斯模型会变得越尖锐,目标像素点受到周围像素影响的范围就越小;
σ 越大,则高斯模型会变得越平缓,目标像素点受到周围像素影响的范围就越大。
SSR算法σ 一般取80-100
SSR算法步骤:
2.2.2 MSR算法(多尺度)
基本理论:
采用几个不同大小的尺度参数对R.G.B三个分量分别单独做滤波后再线性加权归一化就得到了MSR算法。
MSR算法步骤:
同样在这里,对于灰度图像来说,直接对灰度值做上述步骤处理即可;对于彩色图像来说,可以将图像分解为 R、G、B 三副灰度图分别进行以上步骤处理,然后再合成彩色图像。
2.2.3 MSRCR算法
MSRCR 算法在 MSR 算法的基础上增加了一个色彩恢复的步骤
利用彩色图像的各个色彩通道间的比例关系,较好地解决了颜色恒常性的问题
2.3 Retinex算法的频域处理
SSR、MSR、MSRCR 算法都需要用高斯核对图像进行卷积计算,所有算法都要进行下式的计算
在时域中卷积,特别是对于大尺度 σ 时,需要进行的计算量过大,耗时过长。
卷积定理:时域中卷积相当于频域中乘积
假设一个序列的长度为 m,那么此序列做卷积运算是的运算复杂度是 O(m*m);如果将序列做快速傅里叶变换转换到频域再做乘积运算则运算复杂度降低为O(m*logm)。也就是说将计算搬移到频域,卷积计算变成乘法计算,计算量将会大大下降。式 3-13 可以简化为
在做图像增强之前,首先需要对原图像和高斯核做二维离散傅里叶变换
Retinex 算法相关推荐
- Retinex、log对数变换、直方图均衡化区别,边缘增强Retinex算法与拉普拉斯算法联系、均衡化与亮度调节算法、大津阈值计算
1.其中Retinex算法具有的功能:动态范围压缩(即滤掉了低频部分,提取了高频).色调再现(即还有图像色彩):具有锐化.颜色恒常性.动态范围压缩大.色彩保真度高等特点. 从算法公式上的个人理 ...
- 图像增强去雾之直方图均衡化/同态滤波/Retinex算法
图像增强去雾之直方图均衡化/同态滤波/Retinex算法 最近撸了一发图像去雾的算法,主要举四个例子,分别用了全局直方图均衡化,局部直方图均衡化,同态滤波,Retinex增强算法.感兴趣的可以一起讨论 ...
- Retinex算法解读
Retinex是一种常用的建立在科学实验和科学分析基础上的图像增强方法,它是Edwin.H.Land于1963年提出的.就跟Matlab是由Matrix和Laboratory合成的一样,Retinex ...
- 基于不均匀光照下的颜色校正——retinex算法,通态滤波算法
retinex算法原理及算法实现 Retinex是一种常用的建立在科学实验和科学分析基础上的图像增强方法,它是Edwin.H.Land于1963年提出的.就跟Matlab是由Matrix和Labora ...
- 【图像增强】基于matlab双边滤波retinex算法暗光图像增强【含Matlab源码 2305期】
⛄一.简介 1 Retinex 1.1 理论 Retinex理论始于Land和McCann于20世纪60年代作出的一系列贡献,其基本思想是人感知到某点的颜色和亮度并不仅仅取决于该点进入人眼的绝对光线, ...
- Retinex算法在暗光增强应用以及Python实现
图像暗光增强(一) Retinex简介 Single Scale Retinex(SSR) 多尺度MSR Multi-Scale Retinex 带色彩恢复的多尺度MSR,即MSRCR(Multi-S ...
- 【图像增强】双边滤波retinex算法暗光图像增强【含Matlab源码 2305期】
⛄一.简介 1 Retinex 1.1 理论 Retinex理论始于Land和McCann于20世纪60年代作出的一系列贡献,其基本思想是人感知到某点的颜色和亮度并不仅仅取决于该点进入人眼的绝对光线, ...
- 图像增强算法四种,图示与源码,包括retinex(ssr、msr、msrcr)和一种混合算法
from:http://blog.csdn.net/onezeros/article/details/6342661 两组图像:左边较暗,右边较亮 第一行是原图像,他们下面是用四种算法处理的结果 依次 ...
- 关于Retinex图像增强算法的一些新学习
关于Retinex图像增强算法的一些新学习. - Imageshop 时间 2014-06-26 16:50:00 博客园精华区 原文 http://www.cnblogs.com/Imagesh ...
- matlab Retinex图像增强算法
Retinex理论在彩色图像增强.图像去雾.彩色图像恢复方面拥有很好的效果,下面介绍一下我对该算法的理解. Retinex理论 Retinex理论始于Land和McCann于20世纪60年代作出的一系 ...
最新文章
- 用python写一个记账小程序_python实现日常记账本小程序
- 如何让word档在同一页上显示两列或更多列
- python2的xrange比range的优点_python相对于range应该更倾向于实用xrange吗
- 域名授权系统源码 网站源码授权系统_单域名授权系统
- 互联网之道,看电商的数据化管理方案
- 解决linux中xorg占用gpu问题
- Unicode16 与 UTF-8编码之间的转换
- PHP 开源搜索引擎Yioop! 0.80 发布
- 玩转大数据系列之二:数据分析与处理
- centos中多台主机免密登录_mac ssh 免用户名密码远程登录 linux 方法
- 长尾关键词生成工具-长尾关键词生成软件免费-关键词组合工具
- 2022年java知识点汇总,面试大全!超级全面,逐步完善!
- 【强档推荐】动漫初音未来Ⅱ主题
- Java文件操作——简单文件搜索
- webpack压缩图片
- 10 16 进制 转换 c语言,求一段 16进制转10进制 C语言代码。 被转换的16进制数是 0x**型,转换后为10进制数。...
- 修改mdf ldf文件权限修改方法
- 算法分析——算法的渐进效率分析 和 渐进符号大O、大Ω、大θ、小o、小ω
- 样本量太小怎么做结构方程模型?
- 让瓶子里的小人跟你互动,它是怎么做到的?