直方图均衡化

作用：这种方法通常用来增加许多图像的全局对比度，尤其是当图像的有用数据的对比度相当接近的时候。通过这种方法，亮度可以更好地在直方图上分布。这样就可以用于增强局部的对比度而不影响整体的对比度，直方图均衡化通过有效地扩展常用的亮度来实现这种功能。

优点：这种方法对于背景和前景都太亮或者太暗的图像非常有用，这种方法尤其是可以带来X光图像中更好的骨骼结构显示以及曝光过度或者曝光不足照片中更好的细节。这种方法的一个主要优势是它是一个相当直观的技术并且是可逆操作，如果已知均衡化函数，那么就可以恢复原始的直方图，并且计算量也不大。

缺点：它对处理的数据不加选择，它可能会增加背景噪声的对比度并且降低有用信号的对比度。

具体理论

代码实现

 

#include "stdafx.h"
#include "cv.h"
#include "highgui.h"
#include <memory.h>
#include <stdio.h>
#include <iostream>
using namespace std;
int img[300][300];
int times[256];
int imgChanged[300][300];
void getPointVlaue(IplImage *pImage)
{
for(int i = 0;i<pImage->width;i++)
{
for(int j = 0;j<pImage->height;j++)
{
img[i][j] = (int)((uchar*)(pImage->imageData+i*pImage->widthStep))[3*j];
}
}
}
void getEveryPiontTimes(IplImage *pImage)
{
int i;
for(i = 0;i<pImage->width;i++)
{
for(int j = 0;j<pImage->height;j++)
{
times[ (int)((uchar*)(pImage->imageData+i*pImage->widthStep))[3*j] ]++;
}
}

/* 
**得到该灰度图像的灰度值出现次数
int sum = 0;
for(i=0;i<255;i++)
{
if(times[i]!=0)
{
sum=sum+times[i];
}
printf("%d %d| ",i,times[i]);
printf("%d ",sum);
}
*/
}
int Round(float CDF,float minCDF,float width,float height)
{
float round = (CDF - minCDF)/(width*height - minCDF)*255;
return round;
}

 



void getCDF(IplImage *pImage)
{
int CDF[256],i,sum = 0,flag = 0,minCDF = 0;
for(i=0;i<255;i++)
{
if(times[i]!=0)
{
sum = times[i] + sum;
CDF[i] = sum;
if(flag == 0)
{
minCDF = sum;
flag = 1;
}
}
else
{
CDF[i] = 0;
}
}

for( i = 0;i<pImage->width;i++)
{
for(int j = 0;j<pImage->height;j++)
{
img[i][j] = Round(CDF[(int)((uchar*)
(pImage->imageData+i*pImage->widthStep))[3*j]],minCDF,pImage->width,pImage->height);
((uchar*)(pImage->imageData+i*pImage->widthStep))[3*j] = img[i][j];
}
}
}
int main(int argc, char* argv[])
{
const char* imagepath = "C:\\Users\\微安\\Pictures\\test.jpg";
IplImage *pImage;
pImage=cvLoadImage(imagepath,0);
cvNamedWindow("AtBegining");
cvShowImage( "AtBegining", pImage);
// getPointVlaue(pImage);
getEveryPiontTimes(pImage);
getCDF(pImage);
cvNamedWindow("AfterChanged");
cvShowImage( "AfterChanged", pImage);
cvWaitKey();
cvReleaseImage(&pImage);
cvDestroyWindow("AtBegining");
cvDestroyWindow("AfterChanged");

return 0;
}