直方图均衡化

我们可以认为一副高质量的图像应该平均使用所有像素强度，直方图均衡化，即使得直方图尽可能的平坦；equalizeHist可以用来得到均衡化后的直方图，该函数返回类型是void;

#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;
}