直方图均衡化

#include <opencv2\core\core.hpp>
#include <opencv2\imgproc\imgproc.hpp>
#include<opencv2\opencv.hpp>
#include<iostream>
using namespace cv;

//计算单通道灰度图直方图
Mat Hist(Mat& image){
Mat hist;
float rang[2]={0.0,255.0};
int channels[1]={0};
int dims=1;
int histSize[1]={256};
const float* ranges[1]={rang};
calcHist(&image,1,channels,Mat(),hist,dims,histSize,ranges);
return hist;
}

//画出直方图
Mat GetHistimage(Mat& image){
Mat hist=Hist(image);
Mat histimage(256,256,CV_8U,Scalar(255));
double maxVal=0;
double minVal=0;
cv::minMaxLoc(hist, &minVal, &maxVal, 0, 0);
int hpt = static_cast<int>(0.9*256);
for(int i=0;i<256;i++){
float binVal = hist.at<float>(i);
int intensity = static_cast<int>(binVal*hpt/maxVal);
rectangle(histimage,Point(i,256),Point(i,256-intensity),Scalar::all(0));
}
return histimage;
}

Mat Equalize(Mat& image){
Mat result;
equalizeHist(image,result);
return result;
}

int main(){
Mat image=imread("F:\\opencv_test\\6.tif",0);
Mat originalHist=GetHistimage(image);
Mat eImage=Equalize(image);
Mat eImageHist=GetHistimage(eImage);

namedWindow("eImage",0);
imshow("eImage",eImage);

namedWindow("eImageHist",0);
imshow("eImageHist",eImageHist);

namedWindow("original",0);
imshow("original",image);

namedWindow("originalHist",0);
imshow("originalHist",originalHist);

waitKey(0);
return 0;
}