一种边缘保持最小化的高动态范围图像色调映射方法
【专利摘要】本发明公开了一种边缘保持最小化的高动态范围图像色调映射方法，具体按照以下步骤实施：步骤1：获取高动态范围图像的亮度图像L(x,y)；步骤2：采用导向滤波器对L(x,y)进行导向滤波；步骤3：依据简化的摄影学色调映射变换公式，对高动态范围图像的亮度图像L(x,y)进行变换；步骤4：用迭代的方法计算Zn+1(x,y)；步骤5：亮度图像Z(x,y)＝Zn+1(x,y)，计算误差e，如果e＞T，T为阈值，转入步骤4，否则，转入步骤6；步骤6：完成高动态范围图像色调映射。本发明一种边缘保持最小化的高动态范围图像色调映射方法可以将高动态范围图像在低动态范围设备上进行显示输出，既兼顾了对原图像的逼近，又能较好的保持原图像的边缘信息，显示的图像具有对比度高、细节再现效果好的特点。
【专利说明】
一种边缘保持最小化的高动态范围图像色调映射方法
技术领域
[0001] 本发明属于数字图像处理技术领域，具体涉及一种边缘保持最小化的高动态范围 图像色调映射方法。
【背景技术】
[0002] 自然场景中的动态范围通常要跨越1014数量级，随着成像技术的发展，高动态范围 图像的获取越来越方便，其被广泛应用于影视、游戏、军事、天文及医学检测等领域。与此同 时，高动态范围图像的显示设备发展还比较缓慢且价格较贵。为了能在低动态范围的显示 设备上较好的再现高动态范围图像的细节信息，常需要采用一定的方法压缩高动态范围图 像的动态范围，使之动态范围能够匹配低动态范围的显示设备，这个过程称为色调映射 (tone mapping)。
[0003] 目前主流的色调映射方法可以分为两类：（1)全局(Global)算子；（2)局部(Local) 算子。全局算子对图像中具有相同像素值的像素应用同一种变换，不考虑其周围像素;局部 映射算法在进行色调映射时，参考像素的邻域信息对其进行变换，考虑到了图像的局部相 关性，因而能取得较好的映射效果，成为目前关注的重点。目前虽然有较多的局部色调映射 方法，但其对细节的保留仍然不够好。
【发明内容】
[0004] 本发明的目的是提供一种边缘保持最小化的高动态范围图像色调映射方法，解决 了现有局部色调映射方法存在的图像细节不够清晰的问题。
[0005] 本发明所采用的技术方案是，一种边缘保持最小化的高动态范围图像色调映射方 法，具体按照以下步骤实施：
[0006] 步骤1:获取高动态范围图像的亮度图像L(x，y);
[0007] 步骤2:采用导向滤波器对L(x，y)进行导向滤波，得到滤波后的结果为b(x，y);
[0008] 步骤3:依据简化的摄影学色调映射变换公式，对高动态范围图像的亮度图像L(x， y)进行变换，变换结果为a(x，y);
[0009] 步骤4:用迭代的方法计算Zn+1(x，y)，n表示迭代的次数；
[0010] 步骤5:亮度图像Z(x，y) = Zn+1(x，y)，计算误差e，如果e>T，T为阈值，转入步骤4， 否则，转入步骤6;
[0011] 步骤6:完成高动态范围图像色调映射。
[0012] 本发明的特点还在于，
[0013] 步骤1获取高动态范围图像的亮度图像L(x，y)为：
[0014] L(x,y)=0.2126R+0.7152G+0.0722B (1)
[0015] 其中，R、G、B分别为高动态范围图像在像素点(x，y)处的红、绿、蓝颜色值。
[0016] 步骤2中在对L(x，y)进行导向滤波时需要设置两个参数:窗口半径r和正则化因子 ε〇
[0017]步骤3根据遮光-增光原理，将摄影学色调映射变换公式进行简化，对高动态范围 图像的亮度图像L(x，y)进行变换，变换结果a(x，y)为：
[0019]步骤4具体为:通过下式的最小化求得压缩后的亮度图像Z(x，y):
[0021]其中，λ为正则化参数；
[0022]迭代的计算公式为：
[0024]其中η表示迭代的次数；
[0025]迭代时取初始值Z°(X，y) = L(X，y)。
[0026] 正则化参数λ取值范围为1~2。
[0027] 步骤5具体为：
[0028]亮度图像 2(、7)=211+1&，7)，误差6为 ：
[0030]如果e > Τ，Τ为阈值，转入步骤4，否则，转入步骤6。
[0031] 步骤6具体为：
[0032]结合原高动态范围图像的亮度图像L(x，y)和R、G、B三个颜色值，利用最后得到的 低动态亮度图像Z(x，y)，完成高动态范围图像色调映射：
[0034]其中1?。、6。、8。分别表示色调映射后得到的低动态范围图像的各像素点处的红、绿、 蓝颜色值，γ表示Gamma变换的参数。
[0035]本发明的有益效果是:本发明一种边缘保持最小化的高动态范围图像色调映射方 法，可以将高动态范围图像在低动态范围设备上进行显示输出，在对高动态范围图像亮度 信息进行压缩时，既兼顾了对原图像的逼近，又考虑了保持原图像的边缘信息，因而显示的 图像具有对比度高、细节再现效果好的特点。
【附图说明】
[0036]图1是本发明高动态范围图像色调映射方法的流程图；
[0037]图2是本发明高动态范围图像色调映射方法处理后显示出的高动态范围图像； [0038]图3是采用基于双边滤波的色调映射方法处理后的高动态范围图像。
【具体实施方式】
[0039]下面结合附图和【具体实施方式】对本发明进行详细说明。
[0040] 本发明一种边缘保持最小化的高动态范围图像色调映射方法，流程图如图1所示， 具体按照以下步骤实施：
[0041] 步骤1:首先读入高动态范围图像，然后获取高动态范围图像的亮度图像L(x，y):
[0042] L(x,y)=0.2126R+0.7152G+0.0722B (1)
[0043] 其中，R、G、B分别为高动态范围图像在像素点(x，y)处的红、绿、蓝颜色值；
[0044] 步骤2:采用导向滤波器对L(x，y)进行导向滤波，滤波时的引导图就选择为L(x，y) 本身，得到滤波后的结果为b(x，y)。在对L(x，y)进行导向滤波时需要设置两个参数:窗口半 径r和正则化因子ε。窗口半径r决定了滤波过程中窗口的大小，关系到参与运算的像素的多 少，对像素最终滤波后的数值有影响。正则化因子ε决定了边缘和平滑区域的划分，该数值 的选择关系到图像中边缘和平滑区域的判断。
[0045] 步骤3:根据摄影中用到的遮光-增光原理，将摄影学色调映射变换公式进行简化， 对高动态范围图像的亮度图像L( X，y)进行变换，变换结果a (X，y)为：
[0047]步骤4:通过下式的最小化求得压缩后的亮度图像Z(x，y):
[0049] 其中λ为正则化参数，取值在1到2之间，用以确定亮度图像边缘保持的程度。
[0050] 上面的式子最小化求解的结果是亮度图像2(^7)，2(^7)是一个函数，因此是一 个变分问题，可以采用迭代的方法对其求解。迭代时取初始值迭代的计 算公式为：
[0052]其中η表示迭代的次数。
[0053]步骤5:设Z(x，y)=Zn+1(x，y)，计算误差e:
[0055] 设定一个阈值T，如果e>T，则不满足迭代的终止条件，转入步骤4;否则，迭代结束 转入步骤6。
[0056] 步骤6:结合原高动态范围图像的亮度和颜色，利用最后得到的低动态亮度图像Z (x，y)，完成尚动态范围图像色调映射：
[0058] 其中1?。、6。、8。分别表示色调映射后得到的低动态范围图像的各像素点处的红、绿、 蓝颜色值，γ表示Gamma变换的参数，取值在1.8到2.5之间。
[0059] 采用本发明一种边缘保持最小化的高动态范围图像色调映射方法，取窗口半径r =8，正则化因子￡=0.01，正则化参数\=1.5，阈值1=100,6&1111]^变换的参数丫=2.2，处理 后显示的高动态范围图像如图2所示，图3为采用基于双边滤波的色调映射方法处理后的高 动态范围图像，通过对比可以看出，采用本发明方法显示的高动态范围图像具有对比度高、 细节显示效果好的特点。
【主权项】
1. 一种边缘保持最小化的高动态范围图像色调映射方法，其特征在于，具体按照以下 步骤实施： 步骤1:获取高动态范围图像的亮度图像L(x，y); 步骤2:采用导向滤波器对L(x，y)进行导向滤波，得到滤波后的结果为b(x，y); 步骤3:依据简化的摄影学色调映射变换公式，对高动态范围图像的亮度图像L(x，y)进 行变换，变换结果为a(x，y); 步骤4:用迭代的方法计算Zn+1(x，y)，n表示迭代的次数； 步骤5:亮度图像Z(x，y)=Zn+1(x，y)，计算误差e，如果e>T，T为阈值，转入步骤4,否则， 转入步骤6; 步骤6:完成尚动态范围图像色调映射。2. 根据权利要求1所述的一种边缘保持最小化的高动态范围图像色调映射方法，其特 征在于，所述步骤1获取高动态范围图像的亮度图像L (X，y)为： L(x,y)=0.2126R+0.7152G+0.0722B (1) 其中，R、G、B分别为高动态范围图像在像素点(X，y)处的红、绿、蓝颜色值。3. 根据权利要求2所述的一种边缘保持最小化的高动态范围图像色调映射方法，其特 征在于，所述步骤2中在对L(x，y)进行导向滤波时需要设置两个参数:窗口半径r和正则化 因子ε。4. 根据权利要求2所述的一种边缘保持最小化的高动态范围图像色调映射方法，其特 征在于，所述步骤3根据遮光-增光原理，将摄影学色调映射变换公式进行简化，对高动态范 围图像的亮度图像L(x，y)进行变换，变换结果a(x，y)为：(2)〇5. 根据权利要求2所述的一种边缘保持最小化的高动态范围图像色调映射方法，其特 征在于，所述步骤4具体为:通过下式的最小化求得压缩后的亮度图像Z(x，y):其中，λ为正则化参数； 迭代的计算公式为：(4) 其中η表示迭代的次数； 迭代时取初始值6. 根据权利要求5所述的一种边缘保持最小化的高动态范围图像色调映射方法，其特 征在于，所述正则化参数λ取值范围为1~2。7. 根据权利要求5所述的一种边缘保持最小化的高动态范围图像色调映射方法，其特 征在于，所述步骤5具体为： 亮度图像2(乂，7)=211+1(1，7)，误差6为：如果e > T，T为阈值，转入步骤4，否则，转入步骤6。8.根据权利要求2所述的一种边缘保持最小化的高动态范围图像色调映射方法，其特 征在于，所述步骤6具体为： 结合原高动态范围图像的亮度图像L(x，y)和R、G、B三个颜色值，利用最后得到的低动 态亮度图像Z(x，y)，完成高动态范围图像色调映射：其中心、6。、8。分别表示色调映射后得到的低动态范围图像的各像素点处的红、绿、蓝颜 色值，丫表示Gamma变换的参数。
【文档编号】G06T5/00GK105869112SQ201610247395
【公开日】2016年8月17日
【申请日】2016年4月20日
【发明人】张二虎, 李 浩, 段敬红
【申请人】西安理工大学