基于无穷范数约束与最大熵原则的色调映射方法
【技术领域】
[0001] 本发明属于数字图像处理技术领域，涉及将高动态的数字图像映射至低动态范 围，可用于完成高动态数字图像的低动态显示过程。
【背景技术】
[0002] 随着数字图像技术的发展和普及，传统的低动态范围图像已经无法满足人们日益 增长的需求，高动态的数字图像视觉内容正在增多。如在航空航天领域，人们使用更为高效 的数字传感器捕获动态范围更高的航空图像；在电脑动画游戏领域，人们使用计算机图形 生成技术，产生无限动态范围的创作视觉内容。
[0003] 然而，相较于快速发展并广泛应用的计算机图形生成技术和数字传感器工艺，目 前广泛使用的却是不支持高动态数字图像的传统显示设备。因而，高动态的数字图像往往 需要映射成低动态图像，以完成在传统显示设备上显示。常用的方法是对高动态图像进行 线性映射，即直接舍去高动态图像的低位，但是，这样会舍去很多图像的细节信息，从而造 成图像显示效果的减弱。因而，如何将高动态图像映射成低动态图像，使其能够最优地在低 动态显示器上显示，成为了一个很重要的问题。
[0004] 目前为止，已提出的色调映射算法主要有以下几种：
[0005] 汤姆森特许公司提出的专利申请"基于区的色调映射"，专利申请号为 201080029221，公开了一种用于将高动态范围图像色调映射到低动态范围显示器上的方 法。该方法首先访问高动态范围图像，并将图像分段成不同区域，每个区域由矩阵表示，其 中矩阵的每个元素是像素的权重或者概率。随后，计算每个区域的曝光值并将其应用到每 个区域。最后，将不同区域融合到一起以得出最终的映射图像。然而，该算法是根据特定的 函数曲线，将高动态图像直接映射至低动态，并没有考虑到图像的细节信息以及最优显示 问题。
[0006] 惠普开发有限公司提出的专利申请"用于色调映射图像的方法"，专利申请号为 200980160696,公开了一种用于色调映射的数字图像方法。该方法首先将多个线性强度值 从线性空间映射至非线性空间，随后基于当前高位深度强度值选择边界区间值，将所选择 的边界区间值中的每一个映射到较低位深度非线性空间中，最后将被映射的所选择的边界 区间值存储到计算机可读介质上。然而，该算法是通过使用伽马函数来完成将图像从高动 态范围映射至低动态的非线性空间中，因而，其映射后，低动态图像的显示效果并不是很 好。
[0007] 综上所述，上述两种技术均是简单的采用固定形式的映射函数将高动态图像直接 映射至低动态，虽然对过亮或过暗的单峰图像有较好的显示效果，但映射函数形式单一，可 调节范围仅仅是映射函数的相关参数，无法对所有图像进行有效显示，尤其是对包括有过 亮区域和过暗区域的多峰图像进行有效显示。其次，上述显示方法并没有考虑到图像整体 的信息论意义，没有显式的最优映射表达式。最后，上述显示方法由于没有色调误差保护机 制，在映射的过程中会不加区分地丢弃所有统计特征不明显的内容，若这些内容中包含有 语义信息显著的重要内容，如小目标及隐藏目标等，则会导致重要信息的丢失。
【发明内容】
[0008] 本发明的目的在于针对上述已有技术中的不足，提出一种基于无穷范数约束与最 大熵原则的色调映射方法，以减少高动态图像映射成低动态图像时信息的损失，得到效果 最优的低动态图像。
[0009] 实现本发明目的技术思路是：根据高动态图像的系数分布，建立图论模型，采用带 无穷范数约束的最大熵原则，设计并求解最优映射函数，在保证细节信息的前提下，对图像 中的元素进行非线性映射，并最大化低动态图像系数的信息熵。具体步骤包括如下：
[0010] ⑴依据待映射图像中像素点的值，将像素点分为2m个类别，并遍历待映射 图像所有的像素点，统计出每个类别的出现概率P m(a)，其中，m为待映射图像的位宽， 0 ^ a ^ 2m-l ；
[0011] (2)根据步骤⑴所得的统计信息，建立图论模型，即从2-个类别中选取出2n-l个 节点s k，并与初始节点Stl组成节点集合S，其中，η为映射目标图像的位宽，I < k < 2 n-l，Stl =〇 ；
[0012] (3)以最大熵为原则，建立如下最优化方程：
【主权项】
1. 一种基于无穷范数约束与最大熵原则的色调映射方法，包括如下步骤： (1) 依据待映射图像中像素点的值，将像素点分为2-个类别，并遍历待映射图像所有的 像素点，统计出每个类别的出现概率Pm(a)，其中，m为待映射图像的位宽，O < a < 2m_l ; (2) 根据步骤（1)所得的统计信息，建立图论模型，即从2-个类别中选取出2 n-l个节 点sk，并与初始节点Stl组成节点集合S，其中，η为映射目标图像的位宽，I < k < 2 n-l，Stl = 〇 ； (3) 以最大熵为原则，建立如下最优化方程： 2^-1 max 2]-P (λ;.) log P (,ν;) <1> 其中，Pn(Sj) = ΣΡ>)，Sj< b < s J+1，?彡 j 彡 2n_l，\ ; (4) 根据色调失真需求，对式〈1〉中各节点sk的最大长度进行约束，即得到〈1〉的约束 条件： Sjtl-Sj^ τ 〈2> 其中，0 < j < 2η_1，τ为色调失真控制系数，根据实际需求设定； (5) 使用图论中的K-Links最长路径求解算法，在式〈2>的约束情况下，优化步骤（2) 中各节点Sk的选取，得到使式〈1>取得最大值时的最优集合S b; (6) 根据步骤（5)中求得的最优集合Sb，完成待映射图像中每个像素点的映射操作，即 将待映射图像中像素值为[s h，sh+1)的点的像素值赋为h，得到映射目标图像。
2. 根据权利要求1所述的基于无穷范数约束与最大熵原则的色调映射方法，其中所述 步骤（5)，按如下步骤进行： (5a)从步骤⑴的2-个类别中，选取出1个节点，并与初始节点s ^组成1次节点集合 S1，遍历所有的节点选取情况，求得每种情况时的1次节点最大熵e(1;i)，其中0 < i < 2m-l， e(1;i)表示选取节点为i时，2-Pn(Sj)l〇gPn( Sj)的最大值，其中Sje S1; (5b)设置变量M = 2 ; (5c)从步骤（1)的2m个类别中，选取出满足式〈2>的M个节点，并与初始节点s ^组成 M次节点集合SM，结合已经求得的M-I次节点最大摘e(M_1; ，遍历所有的节点选取情况，求得 每种情况时的M次节点最大熵e(su)，其中O < i < 2m-l，e(su)表示选取的最后一个节点为 i 时，X-Pn(Sj)IogPn(Sj)的最大值，其中 Sm; (5d)设置变量M = M+1，并判断M彡2"是否成立，若成立，则返回步骤（5c)，若不成立， 则进行步骤（5e); (5e)根据求得的2n-l次节点最大熵结果，以取得最大值时对应的选取情 况作为最优集合Sb，其中O < i < 2m-l。
【专利摘要】本发明公开了一种基于无穷范数约束与最大熵原则的色调映射方法，主要用于解决高动态图像低动态映射时信息损失较大，易造成重要语义信息细节丢失的问题。其实现步骤为：1)统计待映射图像中像素值的分布情况；2)根据统计结果建立图论模型；3)以图论模型为基础建立最优化方程；4)根据色调失真需求，设置最优化方程的约束条件；5)使用图论中K-Links最长路径求解算法，求解最优化方程；6)根据最优化方程的解，完成映射操作。本发明可有效降低图像低动态映射时信息的损失，通过加入色调误差保护机制，可确保重要语义信息的细节不会损失，且运算复杂度低，并行性好，易于实时实现，可用于高动态图像的低动态映射及显示。
【IPC分类】G06T5-00
【公开号】CN104835121
【申请号】CN201510205188
【发明人】牛毅, 石光明, 任琼巍, 李甫, 樊春晓, 焦阳, 柳林
【申请人】西安电子科技大学
【公开日】2015年8月12日
【申请日】2015年4月27日