毕设—高动态范围图像渲染算法之基于对数变换的色调映射技术（三）

最新推荐文章于 2023-04-04 09:43:58 发布
最新推荐文章于 2023-04-04 09:43:58 发布
阅读量1.8k 收藏 17
点赞数 2
分类专栏： 毕设 色调映射 文章标签： opencv
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/A_ACM/article/details/106372987
版权

这个算法是我第一个实现的色调映射算法，也是五种算法中原理最清晰，实现最简单的。（原理见：一种自适应对数映射的高对比度图像显示技术）

这个算法来源于论文《Adaptive Logarithmic Mapping For Displaying High Contrast Scenes》是一篇2003年的古老算法。其实这个算法OpenCV 3.0已经实现了，但是我还是根据自己的理解又实现了一下。

实现该算法共分为四步：

第一步：把场景亮度映射到图像亮度

其实完成这步有很多方法，上述论文中提到了两种方法：一、对于静态图像，使用场景的对数平均值作为初始缩放因子，进行压缩。二、对于需要进行交互式的图像处理，可以使用亮度的对数信息的高斯模糊后的结果，作为压缩结果，范围用户自己调，默认是15%。

但是本项目使用了另一种压缩方法，即归一化。具体步骤如下：首先我们先把提取出来的场景信息进行色彩空间转化，得到亮度值；随后我们对亮度值进行归一化处理；

	for(int i = 0; i < rows; i++){
		for(int j = 0; j < cols; j++){

			//通过这么写，不仅可以实现色调映射，还可以实现图像增强
			if(strstr(fileName,".hdr")!=NULL)
			{
				//提取场景的RGB数据
				RGB[0]=result.cols[i*cols*3+j*3+0];
				RGB[1]=result.cols[i*cols*3+j*3+1];
				RGB[2]=result.cols[i*cols*3+j*3+2];
			}
			else
			{
				RGB[2]=(double)re.at<Vec3b>(i,j)[0];
				RGB[1]=(double)re.at<Vec3b>(i,j)[1];
				RGB[0]=(double)re.at<Vec3b>(i,j)[2];
			}
			
			//色彩空间转化
			//转化成XYZ表示，虽然缩放只缩放Y，但是为了后面进行转化，这里把x，z也求出来;
			xyz[i][j].x=Mrgb[0][0]*RGB[0]+Mrgb[1][0]*RGB[1]+Mrgb[2][0]*RGB[2];
			xyz[i][j].y=Mrgb[0][1]*RGB[0]+Mrgb[1][1]*RGB[1]+Mrgb[2][1]*RGB[2];
			xyz[i][j].z=Mrgb[0][2]*RGB[0]+Mrgb[1][2]*RGB[1]+Mrgb[2][2]*RGB[2];

			//获取XYZ的比例
			double rateX=xyz[i][j].x/(xyz[i][j].x+xyz[i][j].y+xyz[i][j].z);
			double rateY=xyz[i][j].y/(xyz[i][j].x+xyz[i][j].y+xyz[i][j].z);

			//归一化
			if(strstr(fileName,".hdr")!=NULL)
				xyz[i][j].y=atan(xyz[i][j].y)*2/PI;
			else
				xyz[i][j].y/=255.0;

			//还原x，y，z的比率
			xyz[i][j].x=xyz[i][j].y/rateY*rateX;
			xyz[i][j].z=xyz[i][j].y/rateY*(1-rateX-rateY);

			Lmax = max(xyz[i][j].y,Lmax);//取出场景最亮
			Lsum += log(xyz[i][j].y+a);//对场景亮度进行对数求和
		}
	}

接下来是第二步，就是运用传说中的自适应函数：
在这里插入图片描述
这个公式看起来很唬人，其实并不难：

//对数变换，本算法的核心
xyz[i][j].y = 1.0 * log(xyz[i][j].y/Lav + 1) / log(2.0 + 8.0*pow((xyz[i][j].y/Lav / Lmax), log(bias) / log(0.5))) / log10(Lmax + 1);

其中Lav是对数平均值，简单来说除以对数平均值，映射效果会更好。为什么会更好，原因在讲基于Retinex的色调映射算法时会讲。

第三步：伽马校正

//伽马校正
RGB[0]=C*pow(RGB[0],r);
RGB[1]=C*pow(RGB[1],r);
RGB[2]=C*pow(RGB[2],r);

第四步：色彩空间还原

//转化成RGB表示
RGB[0] = Mxyz[0][0]*xyz[i][j].x + Mxyz[1][0]*xyz[i][j].y + Mxyz[2][0]*xyz[i][j].z;
RGB[1] = Mxyz[0][1]*xyz[i][j].x + Mxyz[1][1]*xyz[i][j].y + Mxyz[2][1]*xyz[i][j].z;
RGB[2] = Mxyz[0][2]*xyz[i][j].x + Mxyz[1][2]*xyz[i][j].y + Mxyz[2][2]*xyz[i][j].z;

完整代码如下：

Mat Tone_Mapping(char* fileName, double bias, double C, double r)//基于对数变换的HDR
{
	int rows,cols;
	HDRLoaderResult result;
	Mat re;
	//判断图像类型
	if(strstr(fileName,".hdr")!=NULL)
	{
		//获取图像信息
		int ret = HDRLoader::load(fileName, result);
		rows=result.height;
		cols=result.width;
	}
	else
	{
		//加载图像
		re = imread(fileName, 1);
		rows=re.rows;
		cols=re.cols;
	}

	//定义矩阵
	Mat out_image (rows,cols,CV_8UC3);//结果图像
	Mat img(rows,cols,CV_8UC3);//过程图像

	//定义变量
	double PI=acos(-1.0);
	double Ymin=INF,Ymax=-INF,Ysum=0,Yav=0;
	double Lmin=INF,Lmax=-INF,Lsum=0,Lav=0;
	for(int i = 0; i < rows; i++){
		for(int j = 0; j < cols; j++){

			//通过这么写，不仅可以实现色调映射，还可以实现图像增强
			if(strstr(fileName,".hdr")!=NULL)
			{
				RGB[0]=result.cols[i*cols*3+j*3+0];
				RGB[1]=result.cols[i*cols*3+j*3+1];
				RGB[2]=result.cols[i*cols*3+j*3+2];
			}
			else
			{
				RGB[2]=(double)re.at<Vec3b>(i,j)[0];
				RGB[1]=(double)re.at<Vec3b>(i,j)[1];
				RGB[0]=(double)re.at<Vec3b>(i,j)[2];
			}
			

			//转化成XYZ表示，虽然缩放只缩放Y，但是为了后面进行转化，这里把x，z也求出来;
			xyz[i][j].x=Mrgb[0][0]*RGB[0]+Mrgb[1][0]*RGB[1]+Mrgb[2][0]*RGB[2];
			xyz[i][j].y=Mrgb[0][1]*RGB[0]+Mrgb[1][1]*RGB[1]+Mrgb[2][1]*RGB[2];
			xyz[i][j].z=Mrgb[0][2]*RGB[0]+Mrgb[1][2]*RGB[1]+Mrgb[2][2]*RGB[2];

			//获取XYZ的比例
			double rateX=xyz[i][j].x/(xyz[i][j].x+xyz[i][j].y+xyz[i][j].z);
			double rateY=xyz[i][j].y/(xyz[i][j].x+xyz[i][j].y+xyz[i][j].z);

			//归一化
			if(strstr(fileName,".hdr")!=NULL)
				xyz[i][j].y=atan(xyz[i][j].y)*2/PI;
			else
				xyz[i][j].y/=255.0;

			//还原x，y，z的比率
			xyz[i][j].x=xyz[i][j].y/rateY*rateX;
			xyz[i][j].z=xyz[i][j].y/rateY*(1-rateX-rateY);

			Lmax = max(xyz[i][j].y,Lmax);//取出场景最亮
			Lsum += log(xyz[i][j].y+a);//对场景亮度进行对数求和
		}
	}

	Lav = exp(Lsum/(1.0*rows*cols));//求取对数平均值

	Lmax=Lmax/Lav;

	for(int i = 0; i < rows; i++){
		for(int j = 0; j < cols; j++){

			//记录XYZ空间下x，y，z的比率
			double rateX=xyz[i][j].x/(xyz[i][j].x+xyz[i][j].y+xyz[i][j].z);
			double rateY=xyz[i][j].y/(xyz[i][j].x+xyz[i][j].y+xyz[i][j].z);

			//对数变换，本算法的核心
			xyz[i][j].y = 1.0 * log(xyz[i][j].y/Lav + 1) / log(2.0 + 8.0*pow((xyz[i][j].y/Lav / Lmax), log(bias) / log(0.5))) / log10(Lmax + 1);

			//还原x，y，z的比率
			xyz[i][j].x=xyz[i][j].y/rateY*rateX;
			xyz[i][j].z=xyz[i][j].y/rateY*(1-rateX-rateY);

			//转化成RGB表示，色彩空间还原
			RGB[0] = Mxyz[0][0]*xyz[i][j].x + Mxyz[1][0]*xyz[i][j].y + Mxyz[2][0]*xyz[i][j].z;
			RGB[1] = Mxyz[0][1]*xyz[i][j].x + Mxyz[1][1]*xyz[i][j].y + Mxyz[2][1]*xyz[i][j].z;
			RGB[2] = Mxyz[0][2]*xyz[i][j].x + Mxyz[1][2]*xyz[i][j].y + Mxyz[2][2]*xyz[i][j].z;

			if (RGB[0] < 0)RGB[0] = 0; if (RGB[0]>1)RGB[0] = 1;
			if (RGB[1] < 0)RGB[1] = 0; if (RGB[1]>1)RGB[1] = 1;
			if (RGB[2] < 0)RGB[2] = 0; if (RGB[2]>1)RGB[2] = 1;

			//伽马校正
			RGB[0]=C*pow(RGB[0],r);
			RGB[1]=C*pow(RGB[1],r);
			RGB[2]=C*pow(RGB[2],r);

			if (RGB[0] < 0)RGB[0] = 0; if (RGB[0]>1)RGB[0] = 1;
			if (RGB[1] < 0)RGB[1] = 0; if (RGB[1]>1)RGB[1] = 1;
			if (RGB[2] < 0)RGB[2] = 0; if (RGB[2]>1)RGB[2] = 1;

			//修正补偿
			RGB[0] = Transform(RGB[0]);
			RGB[1] = Transform(RGB[1]);
			RGB[2] = Transform(RGB[2]);

			out_image.at<Vec3b>(i, j)[0] = (int)(RGB[2]*255);
			out_image.at<Vec3b>(i, j)[1] = (int)(RGB[1]*255);
			out_image.at<Vec3b>(i, j)[2] = (int)(RGB[0]*255);
		}
	}
	return out_image;
}

这些步骤没有很明显的前后之分，比如在我的项目中，色彩空间还原就放在了伽马校正的前面。

上面用到的色彩空间转化，是RGB到XYZ，其中Y代表亮度，在映射过程中就使用了亮度这一个变量。色彩空间还原那肯定是XYZ到RGB了。

空间转化矩阵如下：
在这里插入图片描述
在这里插入图片描述
运行效果如下：
在这里插入图片描述在这里插入图片描述在这里插入图片描述

A comparative review of tone-mapping algorithms for high dynamic range video（高动态范围视频色调映射算法的比较研究）（—）
ISP行业的忠实粉丝
05-18 1690
研读论文A comparative review of tone-mapping algorithms for high dynamic range video（高动态范围视频色调映射算法的比较研究）（—） 目录 研读论文A comparative review of tone-mapping algorithms for high dynamic range video（高动态范围视频色调映射算法的比较研究）（—） 0、摘要 1、介绍 2. 一般的背景 2.1.LDR和HDR图像之间的差异
A comparative review of tone-mapping algorithms for high dynamic range video（高动态范围视频色调映射算法的比较研究）（二）
ISP行业的忠实粉丝
05-19 2210
研读论文A comparative review of tone-mapping algorithms for high dynamic range video（高动态范围视频色调映射算法的比较研究）（二） 目录 研读论文A comparative review of tone-mapping algorithms for high dynamic range video（高动态范围视频色调映射算法的比较研究）（二） 5. Video tone-mapping 5.1.背景 5.2.时间相干性
论文研究-亮度分区自适应对数色调映射算法.pdf
07-22
色调映射可将高动态范围图像显示在低动态范围显示器上。常用的对数全局色调映射算法由于压缩范围有限容易引起细节丢失，针对高动态范围图像的显示问题给出一种基于亮度分区的自适应对数色调映射算法。将高动态范围图像由RGB颜色空间转换为XYZ颜色空间以提取图像亮度信息，然后将亮度图分为高、中、低三个照度区域。根据区域亮度属性实施对数色调映射实现动态范围局部压缩，并进行融合处理以消除区域交界处的显示效果。同时采用双边滤波技术进行细节补偿。实验结果表明，此算法能有效压缩动态范围并再现真实场景信息，同时可以保留丰富的细节。
毕设—基于OpenCV的高动态范围图像渲染算法研究与实现（一）
A_ACM的博客
05-22 1268
本论文共实现了五种色调映射算法——基于对数变换的色调映射技术、基于直方图均衡化的色调映射技术、基于Retinex的色调映射技术、基于双边滤波的色调映射技术、基于亮度分区与引导滤波的色调映射技术。 在实现这五种算法之前我们需要实现一个功能，就是对HDRI文件的读取，色调映射技术就是把无法在显示器上直接显示的HDR图像转化为可以直接显示的RGB图像，那显然HDRI文件的读取是必不可少的。 读取完成之后，我们把文件中表示图像长宽和色彩信息的数据用结构体记录下来，之后我们对这些数据进行色调映射。 最后就是界面设计。
HDR高动态范围光照渲染
weixin_30496751的博客
07-27 917
1 什么是HDR HDR是High-Dynamic Range的缩写，实时渲染中称之为高动态范围光照渲染，是一种表达超过了显示器所能表现的亮度范围的图像映射技术。 2 为什么要用HDR 通常，显示器能够显示RGB分量在（0,255）之间的像素，也就是我们最终呈现在显示器上的backbuffer format最多只是R8G8B8,每个通道8bit, 这样能够表现出的明暗...
高动态范围
qq_22393417的专栏
11-20 1606
高动态范围 高动态范围(HDR)是一项技术，这项技术可以产生图像，并且所产生图像的亮度的动态范围比使用标准动态范围(SDR)来得更广，它强调对颜色和光泽的真实绘制。 HDR的工作原理 通常在渲染过程中，像素的红绿蓝三个通道的值均是使用0到1之间的8位值进行存储，对于显示设备而言，0表示没有强度，1表示最大强度。这个限制了 颜色值的范围，并不能准确地反映现实生活中人们感知灯光的方式，以及在非常亮或者非常暗的情况下会导致不真实的图像的产生。 在HDR渲染过程中，像素值是浮点型。这意味着像素值可以有更大的选择范围
基于自适应对数映射的局部对比度增强算法
大熊
11-14 3215
基于自适应对数映射的局部对比度增强算算法 光照环境不是很好的环境下所拍摄的图片，要么不是太亮，要么就是太暗，不能很好的兼顾亮区和暗区的细节。解决方式有很多种，一种是使用宽动态的sensor，获取一帧长曝光图像和一帧短曝光图像，再将两帧图像进行融合，这种合成的方式可以获取到比较高动态范围图像。当然使用宽动态的sensor成本会比一般的非宽动态sensor高；另...
高动态范围图像的非配对学习算法
14657高动态范围图像色调映射的非配对学习算法Yael Vinker Inbar Huberman-Spiegeldah Raanan Fattal{yael.vinker，inbar.huberman1，raanan.fattal} @ mail.huji.ac.il以色列耶路撒冷希伯来大学计算机科学与工程...
基于生成对抗网络的无监督域自适应方法和图像渲染数据集的研究
1基于生成对抗网络的无监督KonstantinosBousmalis GoogleBrain英国伦敦konstantinos@google.comNathanSilbermanGoogle Research纽约，NYnsilberman@google.comDavid Dohan，Google Brain，加利福尼亚州山景城ddohan@...
高动态范围（High-Dynamic Range，简称HDR）
吴建明wujianming_110117
05-19 1万+
高动态范围（High-Dynamic Range，简称HDR） 一．HDR介绍 高动态范围（High-Dynamic Range，简称HDR），又称宽动态范围技术，是在非常强烈的对比下让摄像机看到影像的特色而运用的一种技术。 当在强光源（日光、灯具或反光等）照射下的高亮度区域及阴影、逆光等相对亮度较低的区域在图像中同时存在时，摄像机输出的图像会出现明亮区域因曝光过度成为白色，而黑暗区域因曝光不足成为黑色，严重影响图像质量。摄像机在同一场景中对最亮区域及较暗区域的表现是存在局限的，这种局限就是通常所讲的“动态
毕设—高动态范围图像渲染算法之基于Retinex的色调映射技术（七）
A_ACM的博客
06-12 1542
这次和前几次一样，也是只介绍算法思路和实现步骤，再给出本人项目中Retinex的代码，Retinex的原理就不做解释。参考的论文是《Adaptive Local Tone Mapping Based on Retinex for High Dynamic Range Images》 本篇参考论文把算法分为两个部分，第一个部分是全局色调映射；第二部分是局部色调映射。全局自适应这部分原理既简单又容易实现，局部自适应这部分原理不简单但是容易实现。
Tone Mapping---色调映射算法
热门推荐
think ofu
12-18 2万+
1. 【背景 - 高动态图像】 图像的动态范围： 图像的 Max(亮度值) / Min(亮度值) 高动态图像 ： 动态范围大的图像 高动态 VS 普通 ： 高动态图像(位数 > 8位) 普通灰度图像(位数 = 8位 一般情况) 显示器的灰度Only 8位 怎么办 -处理 ： 将高动态图像的颜色 -> 变换 才能显示出来； 另外高动态图像的灰度值分布的很...
Tone Mapping(色调映射) 认识
梁桥江的专栏
02-21 5634
转自：https://www.cnblogs.com/cjhd/p/7530440.html Tone Mapping(色调映射) 认识 目录 【背景 - 高动态图像】 【起源】 【定义】 【举个例子】 【研究分类】 详细 1. 【背景 - 高动态图像】 图像的动态范围： 图像的 Max(亮度值) / Min(亮度值) 高动态图像 ： 动态范围大的图像 高动态 VS 普通 ： 高动态图像(位数 > 8位) 普通灰度图像(位数 = 8位 一般情况)...
局部色调映射(Local Tone Mapping)
程序员深度学习
04-23 5821
重建视觉外观是色调映射的终极目标。色调映射算法在降低高动态图像(HDR)范围的同时着力保护捕捉到的原始图像的外观。色调映射算子分两种策略，一种是全局的，另一种是局部的。 1. 全局映射算子 每一个像素点将会根据它的全图特征和亮度信息进行映射，不管其空间位置几何。全局算子一个比较典型的例子就是色调曲线。全局色调映射在处理12位(12-bit)深度的图像的时候是完全OK的，当图像的动态范围特别高的时候，那就不行了。这是因为所有的像素点都采取同一种方式进行处理，根本就没有管它是在较亮区域还是较暗区域。这样的话，.
Chapter4、色调映射的一般方法
qq_55621259的博客
06-09 1844
对第四章的一些内容进行翻译提炼
ISP-全局色调映射(GTM)算法简单实现
最新发布
notgoodenouf的博客
04-04 1348
色调映射算法，色彩空间转换LAB
Tone-mapping high dynamic range images by novel histogram adjustment
u013049912的博客
08-20 1848
abstract本文提出了一种新的直方图调整方法，用于显示高动态范围图像。 我们首先提出了一种基于全局直方图调整的色调映射算子，该算子能够很好地再现高动态范围图像的全局对比度。 然后利用全局色调映射算子对图像进行分割，在局部区域进行自适应对比度调整，以再现局部对比度，保证较好的质量。 我们证明了我们的方法是快速的，易于使用和一组固定的参数值产生了很好的结果，广泛的图像。 1. Introduction场景、图像或成像设备的动态范围定义为亮度或信号电平的最高与最低之比。 因此，现实世界的动态范围约为14个数量
Chapter3、色调映射
qq_55621259的博客
08-04 2221
色调映射大体上说就是使HDR图像能够在我们的低动态范围的显示器（LDR）显示出来，并且尽量符合人眼视觉体验。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值