【OpenCV图像处理】9、常见的图像灰度变换

时间 2020-01-29 标签 OpenCV图像处理 9 常见的图像灰度变换

1.灰度线性变换ios

图像的灰度线性变换是图像灰度变换的一种，图像的灰度变换经过创建灰度映射来调整源图像的灰度，从而达到图像加强的目的。灰度映射一般是用灰度变换曲线来进行表示。一般来讲，它是将图像的像素值经过指定的线性函数进行变换，以此来加强或者来减弱图像的灰度，灰度线性变换的函数就是常见的线性函数。c++

g(x, y) = k · f(x, y) + d函数

设源图像的灰度值为x，则进行灰度线性变换后的灰度值为y = kx + b (0<=y<=255)，下面分别来讨论k的取值变化时线性变换的不一样效果ui

(1).|k|>1时spa

当k>1时，能够用来增长图像的对比度，图像的像素值在进行变换后所有都线性方法，加强了总体的显示效果，且通过这种变换后，图像的总体对比度明显增大，在灰度图中的体现就是变换后的灰度图明显被拉伸了。3d

(2).|k|=1时指针

当k=1时，这种状况下经常使用来调节图像的亮度，亮度的调节就是让图像的各个像素值都增长或是减小必定量。在这种状况下能够经过改变d值来达到增长或者是减小图像亮度的目的。由于当k=1，只改变d值时，只有图像的亮度被改变了，d>0时，变换曲线总体发生上移，图像的亮度增长，对应的直方图总体向右侧移动，d<0时，变换曲线总体下移，图像的亮度下降，对应的直方图发生水平左移。code

(3).0<|k|<1时orm

此时变换的效果正好与k>1时相反，即图像的总体对比度和效果都被削减了，对应的直方图会被集中在一段区域上。k值越小，图像的灰度分布也就越窄，图像看起来也就显得越是灰暗。blog

(4).k<0时

在这种状况下，源图像的灰度会发生反转，也就是原图像中较亮的区域会变暗，而较暗的区域将会变量。特别的，此时咱们令k = -1,d = 255,能够令图像实现彻底反转的效果。对应的直方图也会发生相应的变化。

相应的程序试下以下：

//实现图像的灰度线性变化
#include <iostream>
#include <opencv2\core\core.hpp>
#include <opencv2\highgui\highgui.hpp>
#include <opencv2\imgproc\imgproc.hpp>

using namespace std;
using namespace cv;

int main()
{
	Mat srcImg = imread("1234.jpg");
	if (!srcImg.data)
	{
		cout << "读入图片失败" << endl;
		return -1;
	}
	imshow("原图像", srcImg);
	double k, b;
	cout << "请输入k和b值：";
	cin >> k >> b;
	int RowsNum = srcImg.rows;
	int ColsNum = srcImg.cols;
	Mat dstImg(srcImg.size(), srcImg.type());
	//进行遍历图像像素,对每一个像素进行相应的线性变换
	for (int i = 0; i < RowsNum; i++)
	{
		for (int j = 0; j < ColsNum; j++)
		{
			//c为遍历图像的三个通道
			for (int c = 0; c < 3; c++)
			{
				//使用at操做符，防止越界
				dstImg.at<Vec3b>(i, j)[c] = saturate_cast<uchar>
					(k* (srcImg.at<Vec3b>(i, j)[c]) + b);

			}
		}
	}
	imshow("线性变换后的图像", dstImg);
	waitKey();
	return 0;
}

当k=1.2,b=50时执行程序的效果以下：

2.灰度对数变换

对数变换的基本形式为

其中，b是一个常数，用来控制曲线的弯曲程度，其中，b越小越靠近y轴，b越大越靠近x轴。表达式中的x是原始图像中的像素值，y是变换后的像素值，能够分析出，当函数自变量较低时，曲线的斜率很大，而自变量较高时，曲线的斜率变得很小。正是由于对数变换具备这种压缩数据的性质，使得它可以实现图像灰度拓展和压缩的功能。即对数变换能够拓展低灰度值而压缩高灰度级值，让图像的灰度分布更加符合人眼的视觉特征。例如进行傅里叶变换后的图像，图像中心绝对高灰度值的存在压缩了低灰度部分的动态范围，因此没法在现实的时候便显出原油的细节。这时就须要使用一个对数变换来对结果图像进行修正，通过适当的处理后，原始图像中低灰度区域的对比度将会增长，暗部细节将被加强。

使用程序进行实现以下：

//实现图像的对数变换，做用是压缩图像较亮区域的动态范围
//使用不一样的方法实现图像的对数变换
//基本公式为 y = clog(1+r)
#include <iostream>
#include <opencv2\core\core.hpp>
#include <opencv2\highgui\highgui.hpp>
#include <opencv2\imgproc\imgproc.hpp>

using namespace std;
using namespace cv;

int main()
{
	Mat srcImage = imread("1234.jpg",0);
	if (!srcImage.data)
	{
		cout << "读入图片错误~" << endl;
		return -1;
	}
	double c;
	cout << "请输入常数c:";
	cin >> c;
	Mat srcImage1(srcImage);
	imshow("原图像", srcImage);
	Mat dstImage1(srcImage.size(), srcImage.type());
	Mat dstImage2 = dstImage1.clone();
	Mat dstImage3 = dstImage1.clone();
	//使用第一种方法进行对数变换，对图像总体进行操做
	//首先计算 1+r,注意，是对每个像素点都进行加1操做
	add(srcImage, Scalar(1.0), srcImage1);
	//转换为32位的浮点数
	srcImage1.convertTo(srcImage1, CV_32F);
	//计算log(1+r)
	log(srcImage1, dstImage1);
	dstImage1 = c * dstImage1;
	//进行归一化处理
	normalize(dstImage1, dstImage1, 0, 255, NORM_MINMAX);
	//convertScaleAbs:先缩放元素再取绝对值，最后转换格式为8bit型
	//在这里不具备缩放功能，做用仅为将格式转换为8bit型
	convertScaleAbs(dstImage1, dstImage1);
	imshow("对数变换图像1", dstImage1);
	/////////////////////////////////////////////////////////////
	//使用第二种方法进行图像的对数变换，对图像的像素进行遍历
	double temp = 0.0;
	for (int i = 0; i < srcImage.rows; i++)
	{
		for (int j = 0; j < srcImage.cols; j++)
		{
			temp = (double)srcImage.at<uchar>(i, j);
			temp = c*log((double)(1 + temp));
			dstImage2.at<uchar>(i, j) = saturate_cast<uchar>(temp);
		}
	}
	//进行归一化处理
	normalize(dstImage2, dstImage2, 0, 255, NORM_MINMAX);
	convertScaleAbs(dstImage2, dstImage2);
	imshow("对数变换图像2", dstImage2);
	//////////////////////////////////////////////////////////////
	//使用第三种方法进行图像的对数变换
	//首先进行图像类型转换
	srcImage.convertTo(dstImage3, CV_32F);
	//图像矩阵元素进行加1操做
	dstImage3 = dstImage3 + 1;
	//图像对数操做
	cv::log(dstImage3, dstImage3);
	dstImage3 = c*dstImage3;
	//图像进行归一化操做
	normalize(dstImage3, dstImage3, 0, 255, NORM_MINMAX);
	convertScaleAbs(dstImage3, dstImage3);
	imshow("对数变换图像3", dstImage3);
	
	waitKey();
	return 0;
}

当c取1时，效果以下：

3.灰度幂次变换与Gamma校订

基于幂次变换的Gamma校订是图像处理中一种很是重要的非线性变换，它与对数变换相反，它是对输入图像的灰度值进行指数变换，进而校订亮度上的误差。一般Gamma校订长应用于拓展暗调的细节。一般来说，当Gamma校订的值大于1时，图像的高光部分被压缩而暗调部分被扩展；当Gamma校订的值小于1时，相反的，图像的高光部分被扩展而暗调备份被压缩。

一般状况下，最简单的Gamma校订能够用下面的幂函数来表示：

其中A是常数，函数的输入和输出都是非负数，当r=1时，为直线变换；当r<1时，低灰度区域动态范围扩大，进而图像对比度加强，高灰度值区域动态范围减少，图像对比度下降，图像总体灰度值增大，此时与图像的对数变换相似。当r>11时，低灰度区域的动态范围减少进而对比度下降，高灰度区域动态范围扩大，图像的对比度提高，图像的总体灰度值变小，Gamma校订主要应用在图像加强。目标检测和图像分析等不一样的领域。

总之，r<1的幂函数的做用是提升图像暗区域中的对比度，而下降亮区域的对比度；r>1的幂函数的做用是提升图像中亮区域的对比度，下降图像中按区域的对比度。

因此，对于灰度级总体偏暗的图像，可使用r<1的幂函数增大动态范围。对于灰度级总体偏亮的图像，可使用r>1的幂函数增大灰度动态范围。

下面使用程序进行简单的Gamma变换：

//幂次变换与Gamma灰度校订
#include <iostream>
#include <opencv2\core\core.hpp>
#include <opencv2\highgui\highgui.hpp>
#include <opencv2\imgproc\imgproc.hpp>

using namespace std;
using namespace cv;

Mat GammaTrans(Mat &srcImag, float parameter);


int main()
{
	Mat srcImage = imread("1234.jpg", 0);
	if (!srcImage.data)
	{
		cout << "读入图片失败！" << endl;
		return -1;
	}
	imshow("原始图像", srcImage);
	//初始化几组不一样的参数
	float parameter1 = 0.3;
	float parameter2 = 3.0;
	Mat dstImage1 = GammaTrans(srcImage, parameter1);
	imshow("参数1下的Gamma变换", dstImage1);
	Mat dstImage2 = GammaTrans(srcImage, parameter2);
	imshow("参数2下的Gamma变换", dstImage2);
	waitKey();
	return 0;
}

Mat GammaTrans(Mat &srcImag, float parameter)
{
	//创建查表文件LUT
	unsigned char LUT[256];
	for (int i = 0; i < 256; i++)
	{
		//Gamma变换定义
		LUT[i] = saturate_cast<uchar>(pow((float)(i / 255.0), parameter)*255.0f);
	}
	Mat dstImage = srcImag.clone();
	//输入图像为单通道时，直接进行Gamma变换
	if (srcImag.channels() == 1)
	{
		MatIterator_<uchar>iterator = dstImage.begin<uchar>();
		MatIterator_<uchar>iteratorEnd = dstImage.end<uchar>();
		for (; iterator != iteratorEnd; iterator++)
			*iterator = LUT[(*iterator)];
	}
	else
	{
		//输入通道为3通道时，须要对每一个通道分别进行变换
		MatIterator_<Vec3b>iterator = dstImage.begin<Vec3b>();
		MatIterator_<Vec3b>iteratorEnd = dstImage.end<Vec3b>();
		//经过查表进行转换
		for (; iterator!=iteratorEnd; iterator++)
		{
			(*iterator)[0] = LUT[((*iterator)[0])];
			(*iterator)[1] = LUT[((*iterator)[1])];
			(*iterator)[2] = LUT[((*iterator)[2])];
		}
	}
	return dstImage;
}

4.分段线性变换

分段线性变换也是一种重要的灰度级变换。对于曝光不足，曝光过分和传感器动态范围都会形成图像表现出低对比度的特征。分段线性变换的做用是提升图像灰度级的动态范围。一般来讲，经过阶段必定比例的最亮像素和最暗像素，并使得中间亮度像素占有整个灰度级，于是可以提升图像的全局对比度。在这种状况下，一般称之为对比度拉伸，直方图裁剪，目前普遍的应用于图像后期处理中。一般使用分段函数来实现。下面先简单介绍一下对比度拉伸技术。

图像的对比度拉伸是经过扩展图像灰度级动态范围来实现的，它能够扩展对应的所有灰度范围。图像的低对比度通常是因为图像图像成像亮度不够、成像元器件参数限制或设置不当形成的。提升图像的对比度能够加强图像各个区域的对比效果，对图像中感兴趣的区域进行加强，而对图像中不感兴趣的区域进行相应的抑制做用。对比度拉伸是图像加强中的重要的技术之一。这里设点(x1,y1)与(x2,y2)是分段线性函数中折点位置坐标。常见的三段式分段线性变换函数的公式以下：

其中k1=y1/x1,k2=(y2-y1)/(x2-x1),k3=(255-y2)/(255-y1)

须要注意的是，分段线性通常要求函数是单调递增的，目的是防止图像中的灰度级不知足一一映射。

分段的灰度拉伸技术能够结合直方图处理技术，从而更加灵活地控制输出图像的直方图分布，对特定感兴趣的区域进行对比度调整，加强图像画质。对于图像灰度集中在较暗的区域，能够采用斜率k<0来进行灰度拉伸扩展；对于图像中较亮的区域，能够采用修了k<0来进行灰度拉伸压缩。

实现代码以下：

//实现对比度拉伸
#include <iostream>
#include <opencv2\core\core.hpp>
#include <opencv2\highgui\highgui.hpp>
#include <opencv2\imgproc\imgproc.hpp>

using namespace cv;
using namespace std;

int main()
{
	Mat srcImage = imread("1234.jpg",0);
	if (!srcImage.data)
	{
		cout << "读入图片错误！" << endl;
		return -1;
	}
	imshow("原始图片", srcImage);
	Mat dstImage(srcImage);
	int rowsNum = dstImage.rows;
	int colsNum = dstImage.cols;
	//图像连续性判断
	if (dstImage.isContinuous())
	{
		colsNum = colsNum*rowsNum;
		rowsNum = 1;
	}
	//图像指针操做
	uchar *pDataMat;
	int pixMax = 0, pixMin = 255;
	//计算图像像素的最大值和最小值
	for (int j = 0; j < rowsNum; j++)
	{
		pDataMat = dstImage.ptr<uchar>(j);
		for (int i = 0; i < colsNum; i++)
		{
			if (pDataMat[i]>pixMax)
				pixMax = pDataMat[i];
			if (pDataMat[i] < pixMin)
				pixMin = pDataMat[i];
		}
	}

	//进行对比度拉伸
	for (int j = 0; j < rowsNum; j++)
	{
		pDataMat = dstImage.ptr<uchar>(j);
		for (int i = 0; i < colsNum; i++)
		{
			pDataMat[i] = (pDataMat[i] - pixMin) * 255 / (pixMax - pixMin);
		}
	}
	imshow("对比度拉伸后的图像", dstImage);
	waitKey();
	return 0;
}

执行后显示效果以下

5.灰度级分层

灰度级分层，也叫作灰度级切片，做用是在整个灰度级范围内将设定窗口内的灰度和其余部分分开。从而突出图像中具备必定灰度范围的区域。大致上来讲，灰度级分层有两种类型，即：清除背景和保持背景。清除背景是将灰度窗口内的像素赋值为较亮的值，而其余部分赋值为较暗的值。通过这样的处理后产生的是二值图像，原图像的细节将所有丢失。而保持背景指的是将灰度窗口内的像素赋值为较亮的值，而其余部分的灰度保持不变。