CN108022223B - 一种基于对数映射函数分块处理融合的色调映射方法 - Google Patents

Abstract

本发明涉及一种基于对数映射函数分块处理融合的色调映射方法，包括以下步骤：步骤100：获得一幅高动态图像；步骤200：将所述高动态图像划分为多个大小相等的局部块；步骤300：将划分的局部块分别使用对数压缩函数进行亮度压缩处理；步骤400：将所述经过亮度压缩处理的局部块使用高斯融合函数进行图像融合；步骤500：将所述图像融合得到的图像使用双边滤波进行细节增强；步骤600：将所述经过细节增强的图像进行伽马校正后输出低动态图像。

Description

一种基于对数映射函数分块处理融合的色调映射方法

技术领域

本发明涉及高动态范围图像技术领域，更具体地，涉及一种基于对数映射函数分块处理融合的色调映射方法。

背景技术

在现实世界中，自然界的亮度范围为10^-6至10⁹cd/m²，对比度高达10¹⁵。而传统的图像传感器的动态范围一般只有10²的数量级，因此在显示高动态对比度图像时会出现严重的对比度失真现象，失真现象的主要表现为细节丢失、高亮度处过曝光或者低亮度处欠曝光等特征。而高动态图像能够展现被传统低动态图像丢失但却能被人类视觉系统感知到的极暗和极亮区域的细节部分，能够正确地表现现实世界中的亮度范围。然而，由于传统的显示设备一般是24位真彩色的，每个颜色通道为8位，而高动态图像每个颜色通道为12位或者更高。因此高动态图像无法在常规显示设备中显示，因此对高动态图像进行色调映射极具意义。

色调映射将高动态图像的亮度压缩到传统显示设备可以接受的范围，同时尽可能地保留原图像的细节部分，最终使得经过映射后的低动态图像和人眼观察到的场景尽可能相同。现有的色调映射算法主要分为两类：全局色调映射算法和局部色调映射算法。

全局色调映射算法对高动态图像的所有像素点都采用相同的变换方式进行处理，而Drago映射算法则是其中的代表之一。全局色调映射算法由于对所有像素点都使用同样的算法，因此运算速度快，实现简单。但是，也正是因为对所有像素点采用相同变换方式进行处理，图像的对比度较局部色调映射算法低，容易丢失细节信息，不能很好地反映出真实场景的对应特点，对处理复杂的场景和对比度范围较大的图像映射效果不佳。

局部色调映射算法针对图像不同的区域进行不同的变换，因此对比全局色调映射算法有较好的效果，然而却有计算复杂，效率低，存在伪影现象且图像不自然的特点。

发明内容

本发明为解决现有技术提供的全局色调映射算法存在的容易丢失细节信息的技术缺陷，及现有技术提供的局部色调映射算法存在的计算复杂、存在伪影现象且图像不自然的技术缺陷，提供了一种基于对数映射函数分块处理融合的色调映射方法。

为实现上述目的，本发明采取的技术方案为：

一种基于对数映射函数分块处理融合的色调映射方法，包括以下步骤：

步骤100：获得一幅高动态图像；

步骤200：将所述高动态图像划分为多个大小相等的局部块；

步骤300：将划分的局部块分别使用对数压缩函数进行亮度压缩处理；

步骤400：将所述经过亮度压缩处理的局部块使用高斯融合函数进行图像融合；

步骤500：将所述图像融合得到的图像使用双边滤波进行细节增强；

步骤600：将所述经过细节增强的图像进行伽马校正后输出低动态图像。

优选地，所述划分的多个局部块均为独立的局部块，相邻的局部块之间不存在重合的情况。

优选地，所述划分的局部块的大小满足以下约束条件：

m≤r/2或n≤c/2

其中，m为局部块的像素点的行数，n为局部块的像素点的列数，r为高动态图像的像素点的行数，c为高动态图像的像素点的列数。

优选地，所述步骤200对高动态图像进行划分后，对各个局部块的最大亮度值和亮度进行调整，具体如下：

1)计算局部块的对数均值：

其中L_w(x,y)表示局部块上坐标为(x,y)的像素点的亮度值，δ表示一个非常小的常数，b表示对数压缩系数，N表示局部块上像素点的数量；

2)利用对数均值计算键值：L_wa＝L_a/((1.0+b-0.85)⁵)；

3)利用键值对局部块的最大亮度值和亮度进行调整：

L_amax＝L_max/L_wa

L_w'＝L_w/L_wa

其中L_amax、L_w'分别表示调整后的局部块的最大亮度值和亮度，L_max、L_w分别表示调整前的局部块的最大亮度值和亮度。

优选地，所述步骤300使用到的对数压缩函数表示如下：

优选地，所述对数压缩系数b的范围是0.4≤b≤1。

优选地，所述步骤400使用高斯融合函数进行图像融合的具体过程如下：

1)按照以下方式计算每个局部块上的像素点的高斯曲线函数的值：

2)使用高斯函数公式计算每个局部块上的像素点的权重：

3)使用每个局部块上的像素点的权重对局部块进行融合：

上述公式中，G_jk(x,y)表示第jk块局部块上坐标为(x,y)的像素点在高度为1的高斯曲面上的函数值；j表示行号，k表示列号；d为融合函数系数；(x_jk,y_jk)表示第jk块局部块的中心位置坐标；n_r、n_c分别表示第jk块局部块包含的像素点的行数和列数；I_jk(x,y)表示第jk块局部块上坐标为(x,y)的像素点的亮度值。

优选地，所述步骤S500使用双边滤波进行细节增强的具体过程如下：

1)计算高动态图像的亮度值在模板Ω中定义域核的值：

2)计算高动态图像在模板Ω中值域核的值：

3)计算经过双边滤波的图像：

4)对融合的图像进行增强：

上述公式中，Ω表示n*n的矩阵蒙版；x、y分别表示当前计算的像素点横坐标和纵坐标；k、l分别表示当前模板Ω中使用的像素点横坐标和纵坐标；L(x,y)表示高动态图像在坐标(x,y)的亮度值，L(k,l)表示模板Ω中在位置(k,l)的亮度值；σ_d表示空域核的调节系数，σ_r表示值域核的调节系数。

优选地，所述步骤600进行的伽马校正为伽马2.2校正。

与现有技术相比，本发明的优点和有益效果是：

1、本发明提供的方法通过对高动态图像进行局部块划分处理，让每个局部块的映射效果更好，最后通过融合处理，消除局部块之间的分离现象，使得图像和谐自然。和传统的全局色调映射方法相比，该映射方法对于动态范围较大的高动态图像能保留图像丰富的细节，在高亮度的地方和低亮度的地方的细节部分不容易丢失，适用场景广；

2、本发明提供的方法通过对高动态图像进行分块处理，因此属于局部色调映射算法。与传统的局部色调映射方法相比，该映射算法能够保留图像丰富的细节的同时不引入伪影，映射后的图像自然度高。

3、本发明提供的方法通过对融合后的图像使用双边滤波进行图像的细节增强，使得图像拥有更加丰富的细节，能够向用户提供更多的信息。

附图说明

图1为基于对数映射函数分块处理融合的色调映射方法的算法流程图。

图2为各个局部块经过亮度压缩后的效果图示例。

图3为融合后的低动态图像。

图4为映射完成后的低动态图像。

具体实施方式

为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合附图和具体实施例，进一步阐述本发明是如何实施的。

图1为本发明提供的基于对数映射函数分块处理融合的色调映射方法的算法流程图，参照图1，该方法步骤包括：

步骤100：首先获得一幅高动态范围(HDR)图像，对于高动态范围图像的大小和动态范围的大小均没有限制。该图像可以从数码设备(照相机、摄像机等等)直接获取或者通过一定的HDR图像合成算法使用多张低动态范围图像进行合成得到。本实施例中，采用多张低动态范围图像进行合成得到高动态图像，该图像的大小为768*512。

步骤200：将高动态图像划分为多个大小相等局部块，划分的局部块必须为m*n像素数的矩阵，同时每个局部块都应独立，不与其他局部块有重合的部分，所有的局部块组合起来是原来的高动态图像。同时，划分的局部块的大小应小于原高动态图像的大小，即：m≤r/2或n≤c/2，m为局部块的像素点的行数，n为局部块的像素点的列数，r为高动态图像的像素点的行数，c为高动态图像的像素点的列数。本实施例中，划分的局部块大小为192*128。

步骤201：对各个局部块的最大亮度值和亮度进行调整，具体如下：

1)计算局部块的对数均值：

其中L_w(x,y)表示局部块上坐标为(x,y)的像素点的亮度值，δ表示一个非常小的常数，b表示对数压缩系数；N表示局部块上像素点的数量；本实施例中，δ取值为10^-6，b取0.85。

2)利用对数均值计算键值：L_wa＝L_a/((1.0+b-0.85)⁵)；

3)利用键值对局部块的最大亮度值和亮度进行调整：

L_amax＝L_max/L_wa

L_w'＝L_w/L_wa

其中L_amax、L_w'分别表示调整后的局部块的最大亮度值和亮度，L_max、L_w分别表示调整前的局部块的最大亮度值和亮度。

步骤300：将划分的局部块分别使用对数压缩函数进行亮度压缩处理：

对于公式中有底数的对数函数，均按其底数运算，对于没有底数的对数函数，均按以2为底数运算。其中对数压缩系数b的范围是0.4≤b≤1，当输入b的值大于1时将按b＝1计算，当输入b的值小于0.4时将按b＝0.4计算。经过亮度压缩后各个局部块的效果示例如图2所示。

步骤400：将所述经过亮度压缩处理的局部块使用高斯融合函数进行图像融合：

1)按照以下方式计算每个局部块上的像素点的高斯曲线函数的值：

2)使用高斯函数公式计算每个局部块上的像素点的权重：

3)使用每个局部块上的像素点的权重对局部块进行融合：

上述公式中，G_jk(x,y)表示第jk块局部块上坐标为(x,y)的像素点在高度为1的高斯曲面上的函数值；j表示行号，k表示列号；d为融合函数系数，本实施例中，d的取值为1000；(x_jk,y_jk)表示第jk块局部块的中心位置坐标；n_r、n_c分别表示第jk块局部块包含的像素点的行数和列数；I_jk(x,y)表示第jk块局部块上坐标为(x,y)的像素点的亮度值。融合的低动态图像的示例如图3所示。

步骤500：将所述图像融合得到的图像使用双边滤波进行细节增强：

1)计算高动态图像的亮度值在模板Ω中定义域核的值：

2)计算高动态图像在模板Ω中值域核的值：

3)计算经过双边滤波的图像：

4)对融合的图像进行增强：

上述公式中，Ω表示n*n的矩阵蒙版；x、y分别表示当前计算的像素点横坐标和纵坐标；k、l分别表示当前模板Ω中使用的像素点横坐标和纵坐标；L(x,y)表示高动态图像在坐标(x,y)的亮度值，L(k,l)表示模板Ω中在位置(k,l)的亮度值；σ_d表示空域核的调节系数，σ_r表示值域核的调节系数。

步骤600：将所述经过细节增强的图像进行伽马2.2校正后输出低动态图像，若对映射效果满意，则映射流程结束，否则就重新选择合适的d值、σd和σr，返回步骤400执行。映射完成的低动态图像示例如图4所示。

最后说明，以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (8)

1.一种基于对数映射函数分块处理融合的色调映射方法，其特征在于：包括以下步骤：

步骤100：获得一幅高动态图像；

步骤200：将所述高动态图像划分为多个大小相等的局部块；

所述步骤200对高动态图像进行划分后，对各个局部块的最大亮度值和亮度进行调整，具体如下：

1)计算局部块的对数均值：

其中L_w(x,y)表示局部块上坐标为(x,y)的像素点的亮度值，δ表示一个非常小的常数，b表示对数压缩系数，N表示局部块上像素点的数量；

2)利用对数均值计算键值：L_wa＝L_a/((1.0+b-0.85)⁵)；

3)利用键值对局部块的最大亮度值和亮度进行调整：

L_amax＝L_max/L_wa

L_w'＝L_w/L_wa

其中L_amax、L_w'分别表示调整后的局部块的最大亮度值和亮度，L_max、L_w分别表示调整前的局部块的最大亮度值和亮度；

步骤300：将划分的局部块分别使用对数压缩函数进行亮度压缩处理；

步骤400：将经过亮度压缩处理的局部块使用高斯融合函数进行图像融合；

步骤500：将所述图像融合得到的图像使用双边滤波进行细节增强；

步骤600：将经过细节增强的图像进行伽马校正后输出低动态图像。

2.根据权利要求1所述的基于对数映射函数分块处理融合的色调映射方法，其特征在于：所述划分的多个局部块均为独立的局部块，相邻的局部块之间不存在重合的情况。

3.根据权利要求2所述的基于对数映射函数分块处理融合的色调映射方法，其特征在于：所述划分的局部块的大小满足以下约束条件：

m≤r/2或n≤c/2

其中，m为局部块的像素点的行数，n为局部块的像素点的列数，r为高动态图像的像素点的行数，c为高动态图像的像素点的列数。

4.根据权利要求1所述的基于对数映射函数分块处理融合的色调映射方法，其特征在于：所述步骤300使用到的对数压缩函数表示如下：

5.根据权利要求4所述的基于对数映射函数分块处理融合的色调映射方法，其特征在于：所述对数压缩系数b的范围是0.4≤b≤1。

6.根据权利要求4所述的基于对数映射函数分块处理融合的色调映射方法，其特征在于：所述步骤400使用高斯融合函数进行图像融合的具体过程如下：

1)按照以下方式计算每个局部块上的像素点的高斯曲线函数的值：

2)使用高斯函数公式计算每个局部块上的像素点的权重：

3)使用每个局部块上的像素点的权重对局部块进行融合：

上述公式中，G_jk(x,y)表示第jk块局部块上坐标为(x,y)的像素点在高度为1的高斯曲面上的函数值；j表示行号，k表示列号；d为融合函数系数；(x_jk,y_jk)表示第jk块局部块的中心位置坐标；n_r、n_c分别表示第jk块局部块包含的像素点的行数和列数；I_jk(x,y)表示第jk块局部块上坐标为(x,y)的像素点的亮度值。

7.根据权利要求6所述的基于对数映射函数分块处理融合的色调映射方法，其特征在于：所述步骤500使用双边滤波进行细节增强的具体过程如下：

1)计算高动态图像的亮度值在模板Ω中定义域核的值：

2)计算高动态图像在模板Ω中值域核的值：

3)计算经过双边滤波的图像：

4)对融合的图像进行增强：

上述公式中，Ω表示n*n的矩阵蒙版；x、y分别表示当前计算的像素点横坐标和纵坐标；k、l分别表示当前模板Ω中使用的像素点横坐标和纵坐标；L(x,y)表示高动态图像在坐标(x,y)的亮度值，L(k,l)表示模板Ω中在位置(k,l)的亮度值；σ_d表示空域核的调节系数，σ_r表示值域核的调节系数。

8.根据权利要求7所述的基于对数映射函数分块处理融合的色调映射方法，其特征在于：所述步骤600进行的伽马校正为伽马2.2校正。

Images