文档章节

基于色彩恒常( color constancy)特性的Frankle-McCann Retinex图像增强。

abcijkxyz
 abcijkxyz
发布于 2016/11/22 16:41
字数 1638
阅读 12
收藏 0

      相关随笔可见:带色彩恢复的多尺度视网膜增强算法(MSRCR)的原理、实现及应用。

     从历史的顺序上讲,本篇应该放在MSRCR之前的,只是由于现在大多论文都是描述的MSRCR,因此我也是先学习MSRCR的。

     今天,无意中找寻一些Retinex资料,搜索到一篇文章《Retinex in matlab》,原以为是MSRCR之类的matlab实现,结果仔细一看,和MSRCR算法的描述完全不同。于是找了一些资料,对Retinex的历史有了进一步了解,以下复制一些论文中的文档以供说明:

     ******************************************************************************************************************************************************************************

    文章地址:http://www.doc88.com/p-993974139685.html

    近年来,在现代生理学和心理学的启发下,美国物理学家Land等人设计出了Retinex(视网膜皮层)模型,并在符合人眼的颜色恒常性理论前提下,提出了基于Retinex的图像增强算法。这种算法建立的基础是假设人眼感知目标的亮度和颜色由环境光照和目标表面的反射光照决定,保持目标的颜色恒常性主要就是估计环境光照,并去除环境光照的影像。这样就得到了准确的目标颜色和亮度信息。由于Retinex算法具有高动态范围压缩、高色彩保真度和良好的局部细节增强等特点,引起了大量的国外学者的兴趣。为准确估计亮度分量,选取计算路径常见有一维和二维的区别。早期,Land提出选取随机路径,并对路径所经过的像素进行累加计算亮度,这种方法的缺陷在于当前像素点的亮度和随机路径上的像素亮度有关和其领域周围的像素亮度几乎无关,使得增强后的图像出现亮度不连续的现象。随后出现了两种迭代分段线性路径,即McCann99 Retinex Frankle-McCann Retinex,相比随机路径的Retinex而言,他们的增强效果较好。然而,他们的共同缺陷在于计算复杂度比较高,且迭代次数的选择对增强效果的影响至关重要。2004年,Ciurea和Funt发表了自动选择迭代次数的论文后,这一难题才得以解决。 Land分析了一维路径选择的缺陷,提出了二维路径的选择方式,即中心/环绕Retinex算法。这种方式下,当前像素点周围领域内像素亮度值作为计算依据,且Jobson等人严谨的证明了高斯卷积函数满足中心/环绕函数要求...........................................(以下就是SSR\MSR\MSRCR...等等)。

     ******************************************************************************************************************************************************************************

    由以上描述可见,Frankle-McCann Retinex采用了一维路径寻找方式估计亮度,而MSRCR是二维的,后者则更为复杂。

    关于FMR算法的原理呢,我一下子也讲不清,这里稍微对算法的一些过程进行简单的介绍吧。

    首先,算法的输入需要是[0,1]范围的浮点数,并且是按对数分布的,这个通常需要将[0,255]按对数方式量化到[0,1]范围内。

    然后按照下图所示的路径对每个点的数据进行相比、想乘、复位和平均操作。

                   

      这个我实在是说不清楚了,给几个链接大家自己看看吧:

      基于心理物理学评价和偏爱映射的高动态范围图像的色调映射算子的设计

        基于色彩恒常性的图像去雾技术

      Retinex in Matlab

      最后把处理的数据从对数域转换到[0,255]范围内(这里需要使用exp函数哦)。

      我们还是把时间放在算法效果的简要分析上。

     

     一般情况下,可以认为FMR算法只有迭代次数一个参数。由上图可见,迭代次数越小,图像的对比度越小,连读越亮,动态压缩范围越小。迭代次数越大,图像的对比度越强烈。从直方图上看,迭代次数小时,直方图会聚集在一起,随着迭代次数等增加,直方图逐渐分布到全部动态范围内。 当迭代次数特别大时,图像会和原始图像很接近。

      我们在贴一些该算法处理的结果。

          

          

          

                       原图                                                迭代两次                                                迭代10次

     可见,对一些偏暗的图像的增强效果还是很明显。

     经过测试证明,这个算法对于我们在正常光照下拍摄的照片的处理有意想不到的效果,可以产生类似HDR的效果,贴出如下:

       

       

                               原图                                                                              处理后

   另外,该算法对偏色图像也有一定的纠偏能力,比如下图。

                                                                            

      

                                 原图                                                                             处理后

   贴出一段比较粗糙的实现该函数的matlab代码:

 

function Test()
    clear all
    rgb=imread('c:\222.jpg');%需要处理的图片
    m=size(rgb,1);
    n=size(rgb,2);
    rr=zeros(m,n);
    gg=zeros(m,n);
    bb=zeros(m,n);
    for i=1:m
        for j=1:n
            rr(i,j)=logm(double(rgb(i,j,1))+eps);
            gg(i,j)=logm(double(rgb(i,j,2))+eps);
            bb(i,j)=logm(double(rgb(i,j,3))+eps);
        end
    end
    rr=rr/max(max(rr(:)));
    gg=gg/max(max(gg(:)));
    bb=bb/max(max(bb(:)));
    rrr= retinex_frankle_mccann(rr, 4);
    ggg= retinex_frankle_mccann(gg, 4);
    bbb= retinex_frankle_mccann(bb, 4);
    for i=1:m
        for j=1:n
            rrr(i,j)=round(exp(rrr(i,j)*5.54));
            ggg(i,j)=round(exp(ggg(i,j)*5.54));
            bbb(i,j)=round(exp(bbb(i,j)*5.54));
        end
    end
    rgb=cat(3,uint8(rrr),uint8(ggg),uint8(bbb));
    rgb=max(min(rgb,255),0);
    imshow(rgb);
end

function [ Retinex ] = retinex_frankle_mccann( L, nIterations )
    global RR IP OP NP Maximum
    RR = L;
    Maximum = max(L(:));                                 % maximum color value in the image
    [nrows, ncols] = size(L);
    shift = 2^(fix(log2(min(nrows, ncols)))-1);          % initial shift
    OP = Maximum*ones(nrows, ncols);                     % initialize Old Product
    while (abs(shift) >= 1)
        for i = 1:nIterations
            CompareWith(0, shift);                         % horizontal step
            CompareWith(shift, 0);                         % vertical step
        end
        shift = -shift/2;                                 % update the shift
    end
    Retinex = NP;
end

function CompareWith(s_row, s_col)
    global RR IP OP NP Maximum
    IP = OP;
    if (s_row + s_col > 0)
        IP((s_row+1):end, (s_col+1):end) = OP(1:(end-s_row), 1:(end-s_col)) + ...
            RR((s_row+1):end, (s_col+1):end) - RR(1:(end-s_row), 1:(end-s_col));
    else
        IP(1:(end+s_row), 1:(end+s_col)) = OP((1-s_row):end, (1-s_col):end) + ...
            RR(1:(end+s_row),1:(end+s_col)) - RR((1-s_row):end, (1-s_col):end);
    end
    IP(IP > Maximum) = Maximum;                          % The Reset operation
    NP = (IP + OP)/2;                                    % average with the previous Old Product
    OP = NP;                                             % get ready for the next comparison
end

    同样,提供个编译好的文件给有兴趣研究该算法的朋友看看效果:

   http://files.cnblogs.com/Imageshop/Frankle_Mccann_Retinex.zip

 

 ***************************作者: laviewpbt   时间: 2013.4.18    联系QQ:  33184777  转载请保留本行信息*************************

 

本文转载自:http://www.cnblogs.com/Imageshop/archive/2013/04/18/3029352.html

共有 人打赏支持
abcijkxyz
粉丝 63
博文 6196
码字总数 1876
作品 0
深圳
项目经理
私信 提问
SSE图像算法优化系列二十:一种快速简单而又有效的低照度图像恢复算法。

 又有很久没有动笔了,主要是最近没研究什么东西,而且现在主流的趋势都是研究深度学习去了,但自己没这方面的需求,同时也就很少有动力再去看传统算法,今天一个人在加,还是抽空分享一个简...

Imageshop
08/11
0
0
图片和图形之扩展色彩内容增强图形(20)

原文 概要 除了标准RGB(sRGB)之外,Android 8.0(API级别26)还引入了对额外色彩空间的色彩管理支持,以 在具有兼容显示的设备上呈现图形。有了这种支持,您的应用程序可以通过Java或本地代...

lichong951
05/26
0
0
用Python做图像处理

用Python做图像处理 最近在做一件比较 evil 的事情——验证码识别,以此来学习一些新的技能。因为我是初学,对图像处理方面就不太了解了,欲要利吾事,必先利吾器,既然只是做一下实验,那用...

彭博
2012/03/02
15.1K
2
15种CSS混合模式让图片产生令人惊艳的效果

如果你曾经广泛地使用过图片编辑应用软件(例如,Adobe Photoshop,Pixelmator,GIMP等),那么你可能对混合模式比较熟悉。如同名称中所暗示的,混合模式是指将上层的图像融入下层图像时采用...

欲思
2014/07/11
8.3K
1
gThumb 3.3.3 发布,Linux 图像浏览器

gThumb 3.3.3 发布,此版本添加了图像工具,改进首选项对话框 UI,修复了 GTK 3.14 等。 最值得关注的改进是添加了两个图像工具:special effects (which includes effects such as Desert,...

oschina
2015/02/12
2K
0

没有更多内容

加载失败,请刷新页面

加载更多

http协议请求头的意义

GET /day31_Http_306/index.jsp HTTP/1.1: GET请求,请求服务器路径为/hello/index.jsp,协议为1.1 请求头 1.Host:localhost:请求的主机名为localhost2.User-Agent:Mozilla/5.0(Windows NT......

潇潇程序缘
21分钟前
4
0
Netty 简单服务器 (三)

经过对Netty的基础认识,设计模型的初步了解,来写个测试,试试手感 上篇也说到官方推荐我们使用主从线程池模型,那就选择这个模型进行操作 需要操作的步骤: 需要构建两个主从线程组 写一个服务器...

_大侠__
32分钟前
5
0
day02:管道符、shell及环境变量

1、管道符:"|" 用于将前一个指令的输出作为后一个指令的输入,且管道符后面跟的是命令(针对文档的操作):cat less head tail grep cut sort wc uniq tee tr split sed awk等) [root@localho...

芬野de博客
42分钟前
11
0
Kubernetes系列——Kubernetes 组件、对象(二)

一、Kubernetes 组件 介绍了Kubernetes集群所需的各种二进制组件。 Master 组件 Master组件提供集群的管理控制中心。Master组件可以在集群中任何节点上运行。但是为了简单起见,通常在一...

吴伟祥
52分钟前
15
0
Flink-数据流编程模型

1、抽象等级 Flink提供了不同级别的抽象来开发流/批处理应用程序。 1) 低层级的抽象 最低层次的抽象仅仅提供有状态流。它通过Process函数嵌入到DataStream API中。它允许用户自由地处理来自一...

liwei2000
今天
12
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部