[发明专利]一种基于频域的反色调映射算法

说明书

技术领域

本发明涉及图像处理领域，具体为一种基于频域的反色调映射算法。

背景技术

HDR(High-Dynamic Range)图像全称是高动态范围图像，与普通图片存储格式不同。普通JPEG，BMP图片均采用8位整型数据来保存单个通道分量，如RGB，每个分量均采用8位数据来保存，单个像素总共24位。HDR图像对比JPEG的优势在于每个通道均采用浮点数来表示，每个通道采用32位来存储，使得图片的动态范围得以最大程度的保留。

HDR图像由于其高对比度等优越特性，作为下一代图像发展的一个方向被各界寄予厚望。但由于HDR显示器研制过程复杂，技术还不够成熟，目前尚未完全普及。因此当前大量的图像仍然是LDR图像，在不久的将来，HDR显示设备大规模应用后，如何将现有的大量LDR图像转换为HDR图像，是近年来研究者们越来越关注的一个问题。

解决从LDR图像获取HDR图像的方法称为反色调映射(iTMO)，可分为全局算法和局部算法。Landis提出了一种基于功率函数的全局算法，作者定义了一个阈值R来确定线性扩展的像素点。Masia等人提出基于gamma指数的扩展方法，主要针对过曝光的LDR图像，通过全局扩展使过曝光部分能够显示出较多细节。Rempel等人对图像进行线性扩展后，再定义一个模板来增强高亮度像素点，可以更好的恢复LDR图像中被截断的亮度。

如果选择了合适的参数，这些方法都可以较好的实现LDR到HDR图像的转换，但它们的不足之处在于没有关注图像细节的部分，使得经反色调映射后的HDR图像细节不够清晰。另外，Rempel虽然用模板对高亮度像素点进行亮度增强，但这个模板对所有图像固定不变，生硬而不具有灵活性。因此有必要提供一种注重图像细节纹理，又能灵活地进行亮度补偿的反色调映射算法，以满足人们对图像质量的要求。

发明内容

本发明的目的是提供一种将LDR图像扩展为HDR图像的有效方法，使图像能够较佳地在HDR设备上显示。

因此，本发明提供了一种基于频域的反色调映射算法方法，包括以下步骤：(1)把LDR图像分解成色彩通道和亮度图像；(2)LDR亮度图像通过一个双边滤波器后分离出图像的高频部分和低频部分，分别代表图像细节和基础图层；(3)对基础图层的高亮度像素点进行亮度增强(BEF)以补偿被截断的亮度，高亮度像素点通过模板来确定；(4)基础图层用指数函数进行反色调映射，扩展动态范围；(5)最后将扩展后的基础图层，图像细节，以及色彩信息合并，生成HDR图像。

步骤(2)中，用双边滤波器对图像的高频部分和低频部分进行分离，通过滤波器后的图层为低频部分，相应的，原亮度图像除以低频部分为原图像高频部分。

步骤(3)中，对高亮度像素点进行亮度增强时，用一个阈值来确定亮度增强模板，具体方法是将图像通过一个size＝max(width，height)/50的均值滤波器，从滤波后的图像中选取最大亮度值T作为阈值，图像中亮度大于T的像素点就是要进行亮度增强的对象，在一个size＝width*height的亮度增强模板中将这些像素点设为1，其余像素点设为0，我们只对模板上的非零像素点进行亮度增强，增强幅度为[0，10]效果较佳。

步骤(3)中，为避免图像在亮度增强后出现亮度变化不连贯，边界不明显以及块效应，我们对亮度增强模板进行腐蚀和膨胀处理，用亮度梯度变化阈值Tr＝0.01来限制模板膨胀边界，同时将模板通过一个双边滤波器，双边滤波器既可以起到平滑的作用，又可以保护边界。

步骤(4)中，采用的线性扩展公式为其中L_w(x)为HDR图像亮度，L_d(x)为LDR图像亮度，L_wa为HDR图像亮度几何平均值，由于HDR图像亮度未知，无法直接求出L_wa，通过实验得L_wa求解公式近似为L_wa＝|L_wa′-0.20577|，其中L_wa′为LDR图像亮度几何平均值。

与现有技术相比，本发明的基于频域的反色调映射算法将图像高频部分和低频部分分开处理，很好地保留了图像的细节，更进一步，我们可以将此细节信息储存起来，需要时可以将它移植到其他轮廓不清晰的图像中以提高图像质量；另外，本发明在对亮度图像进行动态范围扩展前增加亮度增强运算，根据图像特性自动生成亮度增强模板，并对模板进行形态学优化，以补偿由于设备限制或者图像压缩造成的亮度截断。

附图说明