OpenCV直方图均衡化

对于排除一个点的直方图均衡化，可以用以下代码实现

代码：

namespace cvWish
{

int stretch2(const cv::Mat& image,cv::Mat &result,int v) 
{
cv::Mat hist(1, 256, CV_32SC1);

int histSize = 256;  
float range[] = {0, 255};  
const float *ranges = { range };  
const int  channels = 0;  
cv::calcHist(&image, 1, 0,cv::Mat(), hist,1, &histSize, &ranges, true, false);

float scale[256];
float lookupF[256];
cv::Mat lookup(cv::Size(1, 256), CV_8U);

int pixNum = image.cols * image.rows;
for (int i =0; i <256; i++) {
scale[i] = hist.at<float>(i) / pixNum *255;
}

float  numMedian = scale[v];//
numMedian /= 256;

scale[v] = 0;
for (int i =0; i <256; i++) {
scale[i] += numMedian;
}

for (int i =0; i <256; i++) {
//scale[i] = hist.at<float>(i) / pixNum *255;
if (i ==0) 
{
lookupF[i] = scale[i];
} 
else 
{
lookupF[i] = lookupF[i -1] + scale[i];
}
}

for (int i =0; i <256; i++) 
{
lookup.at<uchar>(i) = static_cast<uchar>(lookupF[i]);
}
//lookup.at<uchar>(v) = static_cast<uchar>(lookupF[v]);//写法是错误的，不必修改了！

//cv::Mat result;
cv::LUT(image, lookup, result);

return 1;
}

void equalizeHistExcludeValue(const cv::Mat &mat_src, cv::Mat &mat_dst,int v)
{
stretch2(mat_src,mat_dst,v);
}
}