数字视频处理（二）——直方图均衡化

---

  解决方案资源管理器界面如下：

main.cpp

#include <iostream>
#include <cstdio>
#include <fstream>
#include "YHistogram.h"

using namespace std;

#define width 500
#define height 500


int main()
{
	//读入文件
	ifstream inYUV("seed.yuv", ios::binary);
	ofstream outYUV("HisSeed.yuv", ios::binary);
	if (!inYUV) { cout << "error to open file1!" << endl; }
	if (!outYUV) { cout << "error to open file2" << endl; }

	int TotalSize = width * height * 3;
	int YSize = width * height;
	unsigned char* in_YUV = new unsigned char[TotalSize];
	unsigned char* in_Y = new unsigned char[YSize];
	unsigned char* out_Y = new unsigned char[YSize];
	unsigned char* out_YUV = new unsigned char[TotalSize];

	inYUV.read((char*)in_YUV, TotalSize);
	
	//for (int i = 0; i < TotalSize; i++)
	//{
	//	int y = in_YUV[i];
	//}
	//提取Y分量
	EXTRACT(in_YUV, in_Y, width, height);
	
	//直方图均衡化
	Histogram(in_Y, out_Y, width, height);
	
	//UV分量的计算
	Y2YUV(out_Y, out_YUV, width, height);
	
	//写入YUV文件
	outYUV.write((char*)out_YUV, TotalSize);
	
	inYUV.close();
	outYUV.close();

	return 0;
}

YHistogram.h

#pragma once
void EXTRACT(unsigned char* in_YUV, unsigned char* in_Y, int width, int height);
void Histogram(unsigned char* in_Y, unsigned char* out_Y, int width, int height);
void Y2YUV(unsigned char* out_Y, unsigned char* out_YUV, int width, int height);

YHistogram.cpp

#include <iostream>
#include <cstdio>
#include "YHistogram.h"

using namespace std;

void EXTRACT(unsigned char* in_YUV, unsigned char* in_Y, int width, int height)
{
	int TotalSize = width * height * 3;
	int YSize = width * height;

	for (int i = 0; i < YSize; i++)
	{
		*(in_Y + i) = *(in_YUV + i);
	}
}

void Histogram(unsigned char* in_Y, unsigned char* out_Y, int width, int height)
{
	int TotalSize = width * height * 3;
	int YSize = width * height;
	double* count_p = new double[256];
	int* CCOUNT = new int[256];
	int* COUNT = new int[256];

	//初始化为0
	for (int i = 0; i < 256; i++)
	{
		COUNT[i] = 0;
	}

	//进行直方图均衡化
	
	//先计算总数
	for (int i = 0; i < YSize; i++)
	{
		int y = in_Y[i];
		COUNT[y]++;
	}

	//计算频数
	for (int i = 0; i < 256; i++)
	{
		count_p[i] = (double)COUNT[i] / YSize;
	}

	CCOUNT[0] = int(255 * count_p[0] + 0.5);
	
	//计算累加和
	for (int i = 1; i < 256; i++)
	{
		count_p[i] = count_p[i] + count_p[i - 1];
		CCOUNT[i] = int(255 * count_p[i] + 0.5);
	}

	//进行直方图均衡化
	for (int i = 0; i < YSize; i++)
	{
		out_Y[i] = CCOUNT[in_Y[i]];
	}
}

void Y2YUV(unsigned char* out_Y, unsigned char* out_YUV, int width, int height)
{
	int TotalSize = width * height * 3;
	int YSize = width * height;

	for (int i = 0; i < TotalSize; i++)
	{
		if (i < YSize)
		{
			out_YUV[i] = out_Y[i];
		}
		else
			out_YUV[i] = 128;
	}
}

实验结果

处理前	处理后