用于图像的逆色调映射的方法
技术领域
本发明涉及高动态范围成像领域,并且特别地解决了扩展低动态范围内容的动态范围,以使这些内容准备好由具有特别高峰值亮度的显示设备来显示。
背景技术
显示技术的最新发展已开始以允许对要显示的颜色、亮度和对比度的扩展范围。
允许对图像内容的亮度或明度的范围的扩展技术被称为高动态范围成像(highdynamic range imaging),通常简称为HDR。HDR技术专注于捕捉、处理和显示更宽的动态范围的内容。
虽然已经出现了一些HDR显示设备,并且研发出了能够捕获具有增加的动态范围的图像的图像摄像头,但可用的HDR内容仍然有非常有限。虽然最近的发展保证了在不久的将来HDR内容的本地捕获,但这无法处理(address)现有内容。
为了为HDR显示设备准备传统的(本文称为LDR低动态范围)内容,也可以使用反向(reverse)或逆色调映射算子(iTMO,inverse tone mapping operator)。这样的算法处理图像内容中颜色的亮度信息,用于回复或重新创建原始场景的外观的目的。通常情况下,iTMOs采用传统的(即,LDR)图像作为输入,以全局的方式扩展该图像的颜色的亮度范围,以及随后局部地处理高光或亮区域,以提升图像中颜色的HDR外观。
虽然存在一些iTMO方案,这些方案专注于感知上再现原始场景的外观,并且有赖于对内容的严格假设。另外,文献中提出的多数扩展方法正朝着动态范围的极端增加进行优化。
通常,HDR成像是通过对颜色的亮度的暗值和亮值之间的动态范围的扩展、和量化步骤的数量的组合来定义的。为了实现动态范围更极端的增大,许多方法将全局扩展与增强图像的高光或其它亮区域的外观的局部处理步骤进行组合。文献中提出的已知全局扩展步骤可以是不同于逆S形的线性或分段线性。
为了增强图像中亮的局部特征,创建亮度扩展图是已知的,其中图像的每个像素与扩展值相关联以适用于该象素的亮度。在最简单的情况下,图像中的限幅区可被检测到随后被使用陡峭扩展曲线来扩张,但这样的方案不提供对图像的外观的充分控制。
需要一种新颖的iTMO,其目的是通过选择性地且动态地重映射其颜色的亮度及其对比度范围,来提高图像的视觉美观性。
发明内容
本发明的目的是通过选择性和动态地重映射图像的颜色的亮度及其对比度范围,来增强图像的视觉吸引力。基本上,不同的亮度处理被认为对图像的不同部分是必要的。当在显示侧动态范围扩展的情境中没有可用的场景相关数据时,任何亮度扩展方案需要在尽可能好地再现原始场景的外观、以及同时确保压缩的伪像或噪声不被扩展之间折衷。因此,在如下文所述的图像逆色调映射方法中,扩展在空间上是变化的,因此完全适用于图像内容。
因此,本发明的主题是用于图像的逆色调映射的方法,其中图像的颜色在将亮度从色度分离的颜色空间中表示,方法包括:
-通过对像素的空间邻域中的像素的颜色的亮度值进行低通滤波,来获取所述图像的每个像素(p)的像素扩展指数值E’(p),
-通过对像素的空间邻域中的像素的颜色的亮度值的高频率进行提取,来获取所述图像的每个像素(p)的像素亮度增强值Yenhance(p),
-将所述图像的每个像素(p)的亮度Y(p)逆色调映射到通过下述两项的乘积获得的扩展亮度Yexp(p):以该像素的亮度Y(p)为底数、以针对像素获取的像素扩展指数值E’(p)为指数的幂,和以针对该像素获取的像素亮度增强值为底数、以指数参数c为指数的幂,其中c大于或等于1。
这表示为Yexp(p)=Y(p)E’(p)x[Yenhance(p)]c
指数参数c控制由像素亮度增强值带来细节增强量。因此,更大的c值逐渐增加图像边缘的对比度。c=1.5的值是优选的。
作为亮度值的低通滤波的像素,以获得其像素扩展指数值的空间邻域进行的,这意味着每个像素将获得它自己的扩展指数。
由于在像素的空间邻域中执行了对亮度值的低通滤波来获取其像素扩展指数值,这意味着每个像素将获得其自己的扩展指数。因此,由此方法获得的亮度范围的扩展在空间上是变化的,并且适用于图像内容。
由于在像素的空间邻域中执行了对亮度值的高频的提取来获取其像素亮度增强值,这意味着每个像素将获得其自己的亮度增强,并且这样的亮度增强将有利地适用于至少部分地补偿由执行来扩展亮度的亮度值的低通滤波的空间分量造成的、对图像的细节的平滑。
优选地,低通滤波被执行,使得图像的亮区域内的像素的颜色的增强大于该图像的暗区域内的像素的颜色的增强。
优选地,所述像素扩展指数值E’(p)也通过整形经低通滤波的亮度值来获取,使得这些值被重分布,从而所获取的图像的均值扩展亮度约等于原始图像的均值亮度。
在第一变体中,如果经逆色调映射的图像意欲在具有给定的峰值亮度Dmax的显示设备上再现,则所述像素扩展指数值E’(p)也优选地通过重缩放经低通滤波的亮度值来获取,使得以所述图像的像素中最大亮度max(Y)为底数、以针对具有该最大亮度max(Y)的像素获取的像素扩展指数值为指数的幂等于所述峰值亮度Dmax。
这表示为例如其中E(p)为通过对该像素的空间邻域中的像素(p)的颜色的亮度值进行低通滤波直接获取的值,其中α是可用于控制亮度的整体扩展的重缩放参数。因此,有了项Dmax,图像的亮度处理取决于用于显示该图像的HDR显示设备的能力,特别是其峰值亮度。
该方法随后有利地适用于内容本身和显示设备的亮度范围,并且允许内容的HDR外观的提高,即使显示设备在动态范围内适度增加时被使用。
在第二变体中,所述像素扩展指数值E’(p)优选地通过经由应用到经低通滤波的亮度值的二次函数整形这些值来获取。
优选地,为了计算每个像素的像素扩展指数值,通过计算应用在像素的空间邻域中的第一高斯函数、与应用在该像素的颜色的亮度值的邻域中的第二高斯函数的乘积,来在第一变体中执行所述低通滤波。
这意味着使用应用在像素的空间域的第一高斯函数、以及应用于亮度范围域的第二高斯函数来执行滤波。
优选地,第一高斯函数的标准差σs的值大于或等于3。
第二高斯函数的标准差σr的值应优选地高到足以平滑原始图像中的纹理和噪声,但也低到足以避免该图像的对象之间的交叉边缘。标准差σr的值优选地从0.1max(Y)和0.5max(Y)之间选出,其中max(Y)是原始图像的所有像素中的最大亮度值。
优选地,为了计算像素亮度增强值,通过计算第二低通滤波除以第三低通滤波的比率执行所述亮度值的高频率的提取,其中,第三低通滤波被配置为比第二低通滤波去除图像中更多的细节。
优选地,为了计算像素的像素亮度增强值,所述第二低通滤波通过计算应用于该像素的邻域中的第三高斯函数(其具有第一标准差)、与应用于该像素的颜色的亮度邻域中的第四高斯函数的乘积来执行,并且第三低通滤波通过第二次应用所述第二低通滤波来获取,其中第二低通滤波具有高于第三高斯函数的第一标准差的第二标准差。
通过当在第二次应用所述第二低通滤波时使用第三高斯函数的较高的标准差值,相较于第二低通滤波第三低通滤波除去图像中的更多细节。因此,所获得的比例对应于高通滤波。
优选地,所述颜色的亮度被编码成小波(wavelet)分解的递增级数(1evel),每个级数具有至少一个高频系数(LH)和低频率系数(LL),在第二变体中每个像素的像素扩展指数值被计算为等于以最高级数编码所述像素(p)的颜色亮度的低频系数A(p)=LL,随后在该像素的空间邻域中执行低通滤波。
优选地,为了计算像素的像素亮度增强值,通过计算像素的亮度Y(p)与该像素的低频系数A(p)之差,来执行对位于该像素的空间邻域中的像素的亮度值的高频率的提取。
因此,只有从亮度的小波分解的(最高)级数的最低频率子带除去的高频率得以保留用于对图像的细节的增强。
优选地,计算包括将指数系数应用于该差值。该指数系数优选地等于2。
本发明的主题还是用于图像的逆色调映射的图像处理设备,其中图像的颜色在将亮度从色度分离的颜色空间中进行表示,该图像处理设备包括:
-低通滤波模块,被配置为通过对像素的空间邻域中的像素的颜色的亮度值进行低通滤波,来获取所述图像的每个像素的像素扩展指数值,
-整形/重缩放模块,被配置为对由低通滤波模块提供的经低通滤波的亮度值进行整形/重缩放,
-高频率提取模块,被配置为通过对像素的空间邻域中的像素的颜色的亮度值的高频率进行提取,来获取所述图像的每个像素的像素亮度增强值,
-逆色调映射模块,被配置为将所述图像的每个像素的亮度逆色调映射到通过下述两项的乘积获得的扩展亮度:以该像素的亮度为底数、以由整形/重缩放模块针对像素提供的像素扩展指数值为指数的幂,和以高频率提取模块针对该像素提供的像素亮度增强值为底数、以指数参数c为指数的幂,其中c大于或等于1。
附图说明
阅读下面的描述本发明将得以更清楚地理解,描述通过非限制性示例的方式给出并且参考所附附图,其中:
-图1示出了根据本发明的第一和第二实施例的不同步骤,
-图2示出了要进行逆色调映射的图像的颜色的YUV坐标的小波分解,
-图3示出了该图像的小波分解。
具体实施方式
在图中所示的各种元件的功能可以通过使用专用硬件以及能够与适当的软件结合来执行相关软件的硬件来提供。
应当理解,本发明可以以各种形式的硬件、软件、固件、专用处理器或其组合来实现。术语“处理器”不应当被解释为专指能够执行软件的硬件,并且可以隐含地包括但不限于:数字信号处理器(“DSP”)硬件、用于存储软件的只读存储器(“ROM”)、随机存取存储器(“RAM”),以及非易失性存储设备。显然本发明可以作为硬件和软件的组合来实现。此外,软件可以作为有形地体现在程序存储单元上的应用程序来实现。这样的软件可以采取要集成到另一软件的插件的形式。该应用程序可以被上传到包括任何适当架构的图像处理设备,并且由该图像处理设备执行。优选地,图像处理设备在计算机平台上实现,计算机平台具有诸如一个或多个中央处理单元(“CPU”)、随机存取存储器(“RAM”)、以及输入/输出(“I/O”)接口之类的硬件。计算机平台还可以包括操作系统和微指令代码。本文描述的各种处理和功能可以是微指令代码的一部分或应用程序的一部分,或者是它们的任意组合,其可以由CPU执行。此外,各种其它外围单元可以连接到计算机平台,例如附加的数据存储单元、显示设备、打印单元、……等等。实现根据本发明的方法实施例的图像处理设备可以是能够接收的图像的任何电子设备(例如,电视机、机顶盒、网关、蜂窝电话、平板电脑)的一部分。
用于图像的逆色调映射的图像处理设备(其中图像的颜色在将亮度从色度分离的颜色空间中进行表示)包括:
-低通滤波模块,被配置为通过对像素的空间邻域中的像素的颜色的亮度值进行低通滤波,来获取所述图像的每个像素的像素扩展指数值,
-整形/重缩放模块,被配置为对由低通滤波模块提供的经低通滤波的亮度值进行整形/重缩放,
-高频率提取模块,被配置为通过对像素的空间邻域中的像素的颜色的亮度值的高频率进行提取,来获取所述图像的每个像素的像素亮度增强值,
-逆色调映射模块,被配置为将所述图像的每个像素的亮度逆色调映射到通过下述两项的乘积获得的扩展亮度:以该像素的亮度为底数、以由整形/重缩放模块针对像素提供的像素扩展指数值为指数的幂,以及以高频率提取模块针对该像素提供的像素亮度增强值为底数、以指数参数c为指数的幂,其中c大于或等于1。
现在将使用如上所述的图像处理设备、参考图1来描述用于图像的逆色调映射的方法的第一实施例。
使用例如图像接收器,图像与和每个像素的颜色相关的所有数据一起被接收,其中与图像的每个像素相关联的颜色被编码为三个颜色坐标,即一种颜色坐标用于每个颜色通道R、G和B。
在本实施例的第一准备步骤(图1中未示出)中,如果需要的话,所接收的这些颜色的RGB颜色坐标被以本身已知的方式归一化和可选地线性化,从而所接收到的图像的颜色被表示在显示设备的RGB颜色空间中。这个RGB颜色空间可以被标准化,并且相应的显示设备是虚拟的显示设备。随后,这些颜色被转换到将亮度从色度分离的颜色空间,例如YUV颜色空间。从RGB颜色空间到YUV颜色空间的颜色转换本身是已知的,因此没有详细描述。替代地可以使用将亮度从色度分离的任何其它颜色空间,例如XYZ、Yxy、CIE、Lab。因此,亮度值Y(p)和两个色度值U(p)、V(p)与所接收的图像中的任何像素p的颜色相关联。将要描述的逆色调映射的方法的实施例的目的是通过将指数值E′(p)提供给每个亮度值Y(p),来将这些亮度值扩展到扩展亮度值。
在本实施例的第二步骤中,使用低通滤波模块,对于具有其亮度值为Y(p)的图像的每个像素p,通过对像素p的空间邻域中的像素的亮度值(并且可选地还在亮度值Y(p)的邻域中)进行低通滤波来获取中间像素扩展指数值E(p)。该低通滤波步骤优选使用高斯函数。E(p)例如通过以下等式获得:
其中fs是应用于图像的空间域的第一高斯函数,以及fr是应用于亮度范围域的第二高斯函数,其中Ω是以像素p为中心的图像的窗口的尺寸,pi是此窗口中的像素。窗口尺寸可以是例如5或7。为了计算效率,窗口尺寸的较小的值是优选的。
在根据本发明的方法的第一实施例中,随后的低通滤波是双重的(bilateral)。词语“双重的”指的是滤波既在空间域也在亮度范围域执行的事实。
优选地,第一高斯函数fs的标准差σs的值大于或等于2。
第二高斯函数fr的标准差σr的值应优选地高到足以平滑原始图像中的纹理和噪声,但也低到足以避免该图像的对象之间的交叉边缘。标准差σr的值优选地从0.1max(Y)和0.5max(Y)之间选出,其中max(Y)是原始图像的所有像素中的最大亮度值。
在第一实施例的具体实现方式中,空间高斯函数fs的标准差被设定为σs=3,并且亮度范围高斯函数fr的标准差设定为σr=0.3×max(Y)。
经由这种低通滤波获取的所有中间像素扩展指数值随后形成中间扩展指数映射E(p)。
在本实施例的第三步骤中,使用整形/重缩放模块,中间扩展指数映射是通过使用S形函数(取决于图像上像素扩展指数值的分布)被整形为最终扩展指数映射,和/或基于将用于再现经扩展的图像的显示设备的峰值亮度被重缩放。
虽然如上述计算出的像素扩展指数值以相对比例指示了每个像素的扩展,这些扩展值需要被重缩放使得它们符合一组约束。尽管有可用于再现由逆色调映射方法提供的HDR图像的HDR显示设备的增加的能力,但由该方法获得的扩展HDR图像中的平均亮度应当优选地保持在能够与原始LDR图像兼容的水平。同时,亮度的扩展应考虑到将用于再现扩展HDR图像的显示设备的峰值亮度Dmax,从而适当地扩展亮度。因此,经由该第三步骤,整形/重缩放像素扩展指数值E’(p)通过以下等式获得:
其中参数α可用来控制亮度扩展的整体行为,其中E(p)为如上文在第二步骤中计算出中间像素扩展指数值,其中max(E)是图像的所有像素中E(p)的最大值,其中项允许当Y(p)=max(Y)时有Y(p)E’(p)=Dmax。
参数α的目的是控制亮度扩展的“平面化(flat)”程度。它是的空间上变化的扩展和常量指数之间的加权平衡。α的较高值意味着亮度扩展更为局部,因此导向更极端的结果,而α的较低值导向更接近全局的扩展。值α=0.1实际上已提供了亮点扩展和对中间色调的有效管理之间的很好的权衡。
在第一实施例的具体实现方式中,参数α设定为值α=0.1,因为此值提供了亮点扩展和对中间色调的有效管理之间的良好折衷。
随后所有经整形/重缩放的像素扩展指数值形成最终扩展指数映射E’(p)。
在本实施例的第四步骤中,使用高频率提取模块,原始图像中亮度值的高空间频率被提取,以获取针对图像中亮度值为Y(p)的每个像素p的像素亮度增强值Yenhance(p)。本步骤是例如通过如下的三个子步骤执行的。
在滤波的第一子步骤,基底层Ybase计算如下,Ybase包含图像亮度的较大变化,但不包括噪声或其它小的特征:
Ybase(p)=∑i∈Ω′Y(pi)f′s(||pi-p||)f′r(||Y(pi)-Y(p)||)
使用与上述相同的高斯函数f’s,、f’r,首先在像素p的空间邻域中使用,并且第二次在亮度值Y(p)的邻域中使用,其中Ω’是以像素p为中心的图像的窗口的尺寸,具有与上述相同的优选尺寸,并且其中pi是在窗口中的像素。
在第一实施例的具体的实现方式中,空间高斯函数f’s的标准差现在设置为σ’s=10,并且亮度范围高斯函数f’r的标准差设置为σ’r=0.1×max(Y)。
在滤波的第二子步骤中,由第一子步骤中获得的基底层Ybase还使用高斯函数(此次的高斯函数具有更大的亮度范围标准差从而移除图像中更多的亮度细节)进行第二次滤波,获取Y’base如下:
其中f’s,、f’r是与上述相同的高斯函数,首先在像素p的空间邻域中使用,并且第二次在亮度值Y(p)的邻域中使用,其中Ω”是以像素p为中心的图像的窗口的尺寸,具有与上述相同的优选尺寸,并且其中pi是在窗口中的像素。
可选地,Y’base(p)可以被设置成等于如上定义的E(p),以最小化最终HDR图像中的噪声。
在第一实施例的具体实现方式中,空间高斯函数f’s的标准差现在设置为σ”s=σ’s,并且亮度范围高斯函数f’r的标准差设置为移除图像中更多的亮度细节(即,σ”r>σ’r)。例如σ”r=0.3×max(Y)。
对于图像中的亮度值的高频率的提取的第三和最后子步骤中,对于图像的每个像素,像素亮度增强值Yenhance(p)作为第一子步骤的结果Ybase(p)除以第二子步骤的结果Y’base(p)得出的比例来计算,从而Yenhance(p)=Ybase(p)/Y’base(p)。由于σ”r>σ’r,该比例对应于图像中的像素值的高频率的提取。在不脱离本发明的范围的情况下,可以使用其他方式提取图像中的亮度值的高空间频率。
随后所有最终像素亮度增强值形成亮度-增强映射,由于对亮度值的高频率的提取,当被应用于图像的亮度值时亮度-增强映射将增强其细节,
在本实施例的第五步骤中,使用逆色调映射模块,图像的每个像素p的亮度Y(p)被逆色调映到扩展亮度Yexp(p),扩展亮度Yexp(p)是通过下述两项的乘积获得的:以该像素的亮度为底数、以从上述第三步骤获得的最终像素扩展指数值E’(p)为指数的幂,以及以从上述第四步骤获得的像素亮度增强值为底数、以指数参数c为指数的幂,其中c大于或等于1。这意味着Yexp(p)=Y(p)E(p)x[Yenhance(p)]c。
指数参数c控制由像素亮度增强值带来细节增强量。因此,更大的c值逐渐增加图像边缘的对比度。c=1.5的值是优选的。
当如上所述地扩展图像的亮度时,亮度和对比度的改变可能影响图像中的饱和度和颜色的外观。在扩展其亮度范围时,图像的颜色信息可以在第六可选步骤中进行管理,以保持图像的艺术意图。优选地,颜色饱和度使用的扩展指数值作为指导(guide)进行增强。更具体地,每个像素的颜色饱和度由等于该象素的扩展指数的因子来增强。像素p的颜色饱和度例如通过调节该象素的色度值C(p)来增强,在圆柱YUV空间中计算如下:
并且调整后的色度值C’(p)作为该像素p的扩展指数E’(p)、以及该像素p的色度值C(p)的乘积来计算,从而:
C′(p)=E′(p)xC(p)
在第一实施例的具体实现方式中,将C(p)转换为C’(p)的色度缩放优选地限制为1.5倍,以避免过度饱和的亮点,例如避光爆和亮光。
C’(p)有了新的值,使用从圆柱颜色空间(例如这里的LCH)到YUV空间的惯用转换方式来计算新的色度U’(p)、V’(p):
U’(p)=cos[H(p)]×C’(p)
V’(p)=sin[(H(p)]×C’(p)
其中H(p)为从原始U’(p)、V’(p)如下计算的原始色彩:
H(p)=arctan[V(p),U(p)]
在第五或第六步的结束时,图像的每种颜色的YUV坐标随后被映射到表示YUV颜色空间中的扩展的颜色新的Y’U’V’坐标。
在第一实施例的最终的第七步骤(图1中未示出),扩展颜色的新Y’U’V’坐标以本身已知的方式转化为表示RGB颜色空间中的相同颜色的对应新的R’G’B’坐标。如果需要的话,这些R’G’B’坐标去线性化,以例如代表这样的颜色空间中扩展的颜色:在上述第一步骤中颜色曾在该颜色空间中接收。
如在图1中所述,基于这些扩展颜色的扩展图像现在准备好被发送到峰值亮度为Dmax的显示设备,从而在高动态范围内再现。
优点
通过根据本发明的方法获得的扩展图像尽可能与相同场景的HDR图像看上去近似。所得到的扩展内容相较LDR输入图像具有较高的视觉质量,即使在只能应用于适度扩展的情况下也是如此。这种亮度扩展方法增强了亮图像特征,向观看者递送高亮和光源的外观,同时保留中档值。有赖于用于再现这些扩展图像的显示设备的亮度范围,暗值可以被保留或进一步压缩,以提高图像中的全局对比度。
根据本发明的方法使用空间域中的低通滤波处理来定义用于扩展颜色的亮度的扩展指数值。该低通滤波处理平滑图像的一些细节。有利的是,对细节的移除由对空间域中高空间频率(其用于获取应用到扩展亮度的亮度增强因子)的提取来补偿。换言之,执行亮度值的高通滤波的空间分量,从而至少部分地补偿由亮度值的低通滤波的空间分量造成的、对图像的细节的平滑。
如上所述地对使亮度使用双重的低通滤波器,可以得到有利于图像的亮区域的扩展映射,同时平滑小局部细节(例如如噪声或条纹伪像)。这允许亮区的扩展,但没有放大噪声或纹理信息。
事实上,目前大多数的HDR LCD显示设备的局限性之一是,相较这些设备前端的液晶面板,这些设备随空间变化的背光具有较低的分辨率。因此,尽管这样的HDR显示设备可以达成很高的全局对比度,但在局部对比度受到影响。这在图像的小型高对比度区域(例如亮点)、或者在暗区和亮区间的尖锐过渡处尤为明显。虽然许多方案可用于增强图像的对比度,但在逆色调映射的情境中,需要特别小心以确保噪声或其它伪像不被同时增强。上文描述了用于提取图像的高空间频率的双通双重滤波处理特别适用于增强对比度而不增强噪声和伪像。
当像素扩展指数值E(p)被如上所述地整形为E’(p)时,所获得的扩展图像最好利用用于再现图像的显示设备的可用动态范围。
第一实施例的第一变体
对于图像序列(即视频内容)的逆色调映射,上述方法的实施例顺序地进行应用,序列的一幅图像之后是另一幅图像。如果此方法独立地应用于序列的每个连续图像,可能难以保证所获得的HDR序列的扩展图像之间的时间相干性。事实上,如果场景中存在移动的亮的物体,可能会产生闪烁伪像。为了避免这个问题,优选地通过确保扩展连续镜头的均值亮度保持不变来应用简单的时间处理。对于序列中的图像i(如上述第一实施例所述,已针对该图像获取了扩展亮度值Y’exp.i),经调整的扩展亮度值Y”exp.i被计算如下,并且取代Y’exp.i用于扩展图像:
为了确保在个场景中的强特征不会传播到下一场景,均值在每一镜头的开始处复位。这可以通过采用现有的镜头检测算法来实现。
第一实施例的第二变体
此第二变体涉及从如上述实施例的第二步骤处获得的中间像素扩展指数值E(p),计算经整形/重缩放的像素扩展指数E’(p)。现在根据下述等式,提出根据中间像素扩展指数值E(p)的二次函数来整形像素扩展指数值:
E’(p)=a[E(p)]2+b E(p)+c
需要注意在该等式中,E(p)不像先前的等式中那样由max(E(p))来缩放。
在第一示例中,上述等式中的参数a、b、c可以被设置为常数,并且针对具体的显示设备进行优化。
在第二示例中,上述等式中的参数a、b、c可根据Dmax来计算。下面的等式例如用于计算参数a、b、c的值:
a=pa1exp(pa2Dmax)+pa3exp(pa4Dmax) (s-1)
c=1.3 (s-3)
在等式(s-1)和(s-2)中给出的值的示例示出于表1中。
表1
在第三示例中,参数a、b、c的值取决于像素扩展指数值E(p)本身的值,从而与E(p)的暗值形成区分,并且对于E(p)的亮值更好地控制展后得到的图像的亮点和阴影。
根据本变体的所有经整形/重新调整的像素扩展指数随后形成取代先前映射的最终扩展指数映射E,(p)。
现在将公开用于图像的逆色调映射的方法的第二实施例,使用如上所述、但具体针对实现此实施例配置的图像处理设备。
与第一个实施例相比,唯一的区别涉及:
-第二步骤,对于每个像素p,通过对像素p的空间邻域中的像素的亮度值进行低通滤波来获取像素扩展指数值E(p)。
-第四步骤,提取高空间频率来获取图像的每个像素p的像素亮度增强值Yenhance(p)。
在第二实施例中,颜色作为表示将亮度从色度分离的颜色空间(例如YUV颜色空间)中的这些颜色的三色坐标组来接收。重申,可以替代地使用从色度分离的颜色空间的任何其它颜色空间。因此,亮度值Y(p)和两个色度值U(p)、V(p)被关联至图像的任意像素p的颜色。
更具体地,这些三色坐标组以压缩的格式(其为小波压缩类型)被接收,因此,颜色的亮度Y(p)被编码成小波分解的递增级数,每个级数具有至少一个高频系数LH和低频率系数LL,例如如图2所述,其中LL、LH、HL和HH子带分区指示近似图像以及分别具有水平边缘、垂直边缘和角边缘的图像,例如如图3所示。
为了创建像素扩展指数映射和像素亮度增强映射,在此实施例中仅使用了小波分解的一个级数。
为了实现方法的第二步骤,使用所谓的Haar提升小波(Haar lifting wavelet)实现方式,并且像素p的像素扩展指数值E(p)从LL子带获得。hψ(n)={1/√2;1/√2}和是卷积滤波器。
因此,这意味着在本实施例的第二步骤中,每个像素的像素扩展指数值被计算为等于编码该像素p的颜色亮度的最高级数的低频系数A(p)=LL,随后在该像素的空间邻域中执行低通滤波。
实现方法的第四步骤中,通过计算像素p的亮度Y(p)和该像素的低频系数A(p)的之差,来执行对空间域中亮度的高频率的提取。更精确地,LL子带(其是原始亮度Y(p)分辨率除以4的低频率版本A(p))被进行上采样,并且随后从该原始亮度Y(p)中减去,以如下所述地获得像素亮度增强值Yenhance(p):Yenhance(p)=((Y(p)-cA(p)↑2)+128)/255*2。
因此,只有从LL子带移除高频被保持用作创建亮度增强映射Yenhance的基础。
所有其它步骤类似于在第一实施例中描述的步骤。当然,因为所获得的扩展颜色的新Y’U’V’坐标可能不需要被转换成相应的新R’G’B’坐标,可以除去最终的第七步骤。
虽然相对于特定示例和优选实施例描述了本发明,应当理解,本发明并不限于这些示例和实施例。因此,要求保护的本发明包括来自本文所述的特定示例和优选实施例的变体,这对本技术领域技术人员是显而易见的。尽管一些具体实施例可以被单独描述并要求保护,应该理解本文所描述和要求保护实施例的各种特征可组合使用。权利要求中出现的参考数字仅仅是说明性的,并且应当对权利要求的范围没有限制作用。