发明名称 基于视网膜自适应模型的高动态范围图像色调映射方法
摘要 本发明公开一种基于视网膜自适应模型的高动态范围图像色调映射方法,属于数字图像处理领域。其具体操作步骤为:①输入一幅高动态范围图像;②获得该图像中每个像素p的对数域亮度L(p);③使用双边滤波器对每个像素p的对数域亮度L(p)进行双边滤波,得到其自适应因子Aopl(p);④得到每个像素p的对数域亮度L(p)经过非线性映射后的亮度Lram(p);⑤得到每个像素p分别在红色通道R、绿色通道G、蓝色通道B进行非线性映射后的红色通道值I′(pR)、绿色通道值I′(pG)、蓝色通道值I′(pB);⑥使用I′(pR)、I′(pG)、I′(pB),构造一个新的图像文件,即可实现高动态范围图像的色调映射。本发明方法同时具备了全局方法的速度快、局部方法能够较好地保留图像细节的优点。
申请公布号 CN101901475B 申请公布日期 2011.10.05
申请号 CN201010199100.3 申请日期 2010.06.12
申请人 北京理工大学 发明人 沈建冰;刘衡生
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 代理人
主权项 1.一种基于视网膜自适应模型的高动态范围图像色调映射方法,其特征在于:其具体操作步骤为:第1步:输入一幅高动态范围图像;第2步:根据该图像的红色通道R、绿色通道G、蓝色通道B的通道值,获得该图像中每个像素p的对数域亮度L(p);第3步:使用双边滤波器对每个像素p的对数域亮度L(p)进行双边滤波,得到其自适应因子A<sub>opl</sub>(p);第4步:通过公式2得到每个像素p的对数域亮度L(p)经过非线性映射后的亮度L<sub>ram</sub>(p):<maths num="0001"><![CDATA[<math><mrow><msub><mi>L</mi><mi>ram</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>m</mi><mrow><mo>(</mo><mi>L</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mrow><mo>(</mo><mi>L</mi><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>A</mi><mi>opl</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,m和s是对输入图像进行尺度平移和缩放的参数,0<m≤1,s>0;第5步:分别通过公式3~5得到每个像素p分别在红色通道R、绿色通道G、蓝色通道B进行非线性映射后的红色通道值I′(p<sub>R</sub>)、绿色通道值I′(p<sub>G</sub>)、蓝色通道值I′p<sub>B</sub>):<maths num="0002"><![CDATA[<math><mrow><msup><mi>I</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>m</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><msup><mrow><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>R</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>L</mi><mi>ram</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mi>av</mi></msub><mo>)</mo></mrow></mrow><mi>C</mi></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msup><mi>I</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>G</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>m</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>G</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>G</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>L</mi><mi>ram</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mi>av</mi></msub><mo>)</mo></mrow><mi>C</mi></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msup><mi>I</mi><mo>&prime;</mo></msup><mrow><mo>(</mo><msub><mi>p</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>m</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>s</mi><mo>)</mo></mrow><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>L</mi><mi>ram</mi></msub><mrow><mo>(</mo><mi>p</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>L</mi><mi>av</mi></msub><mo>)</mo></mrow><mi>C</mi></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,I(p<sub>R</sub>)为像素p在原始图像中的红色通道值;I(p<sub>G</sub>)为像素p在原始图像中的绿色通道值;I(p<sub>B</sub>)为像素p在原始图像中的蓝色通道值;L<sub>av</sub>为全部像素p的对数域亮度L(p)的平均值;C为参数,参数C的计算借鉴了对比度感知模型,其通过公式6得到:C=0.413+0.586×K<sup>(1.0+K)</sup>        (6)其中,K是图像的明暗关键值,<img file="FSB00000557519500021.GIF" wi="347" he="135" />L<sub>max</sub>为L(p)中的最大值,L<sub>min</sub>为L(p)中的最小值;第6步:使用第5步得到的每个像素p分别在红色通道R、绿色通道G、蓝色通道B进行非线性映射后的红色通道值I′(p<sub>R</sub>)、绿色通道值I′(p<sub>G</sub>)、蓝色通道值I′(p<sub>B</sub>),构造一个新的图像文件,即可实现高动态范围图像的色调映射。
地址 100081 北京市海淀区中关村南大街5号