网站建设猫腻,做网站时字幕怎么做,黑龙江省建设工程质量协会网站,成都网站建设全美由于图像的直方图表示图像像素灰度值的统计特性#xff0c;因此可以通过两幅图像的直方图特性比较 两幅图像的相似程度。从一定程度上来讲#xff0c;虽然两幅图像的直方图分布相似不代表两幅图像相似#xff0c;但是两幅图像相似则两幅图像的直方图分布一定相似。例如… 由于图像的直方图表示图像像素灰度值的统计特性因此可以通过两幅图像的直方图特性比较 两幅图像的相似程度。从一定程度上来讲虽然两幅图像的直方图分布相似不代表两幅图像相似但是两幅图像相似则两幅图像的直方图分布一定相似。例如在通过插值对图像进行缩放后虽然图像的直方图不会与之前完全一致但是两者之间一定具有很高的相似性因而可以通过比较两幅图像的直方图分布相似性对图像进行初步的筛选与识别。 OpenCV 4 中提供了用于比较两个图像直方图相似性的 compareHist() 函数。 1. compareHist()函数原型
double cv::compareHist(InputArray H1,
InputArray H2,
int method
) H1第一幅图像直方图。H2第二幅图像直方图与 H1 具有相同的尺寸。method比较方法标志。 该函数前两个参数为需要比较相似性的图像直方图由于不同尺寸的图像中像素数目可能不相同为了能够得到两个图像直方图正确的相似性需要输入同一种方式归一化后的图像直方图并且要求两个图像直方图具有相同的尺寸。该函数中第三个参数为比较相似性的方法选择不同的方法会得到不同的相似性系数会将计算得到的相似性系数以 double 类型返回。由于不同计算方法的规则不一因此相似性系数代表的含义也不相同。接下来介绍每种方法比较相似性的原理。 compareHist()函数比较直方图方法的可选择标志参数 1.1 HISTCMP_CORREL 该方法名为相关法 在该方法中如果两个图像直方图完全一致那么计算数值为 1 如果两个图像直方图完全不相关那么计算值为 0 。 其中N 是直方图的灰度值个数。
1.2 HISTCMP_CHISQR 该方法名为卡方法在该方法中如果两个图像直方图完全一致那么计算数值为 0两个图像的相似性越小计算数值越大。 1.3 HISTCMP_INTERSECT 该方法名为直方图相交法 该方法不会将计算结果归一化因此即使是两个完全一致的图像直方图来自于不同图像也会有不同的数值。例如由A 图像缩放后得到的两个完全一样的直方图相似性结果与由 B 图像缩放后得到的两个完全一样的直方图相似性结果可能不相同。但是当任意图像的直方图与 A 图像的直方图比较时数值越大 相似性越高数值越小相似性越低。 1.4 HISTCMP_BHATTACHARYYA 该方法名为巴塔恰里雅距离巴氏距离法 在该方法中如果两个图像直方图完全一致那么计算数值为 0 两个图像的相似性越小计算数值越大。 1.5 HISTCMP_CHISQR_ALT 该方法称为替代卡方法其判断两个直方图是否相似的方法与巴氏距离法相同常用于替代巴氏距离法用于纹理比较。 1.6 HISTCMP_KL_DIV 该方法名为相对熵法又名 Kullback-Leibler 散度法 在该方法中如果两个图像直方图完全一致那么计算数值为 0 两个图像的相似性越小计算数值越大。 2. 示例代码
#include opencv2/opencv.hpp
#include iostream using namespace cv;
using namespace std; void drawHist(Mat hist, int type, string name) //归一化并绘制直方图函数
{ int hist_w 512; int hist_h 400; int width 2; Mat histImage Mat::zeros(hist_h, hist_w, CV_8UC3); normalize(hist, hist, 1, 0, type, -1, Mat()); for (int i 1; i hist.rows; i) { rectangle(histImage, Point(width*(i - 1), hist_h - 1), Point(width*i - 1, hist_h - cvRound(hist_h*hist.atfloat(i - 1)) - 1), Scalar(255, 255, 255), -1); } imshow(name, histImage);
}
//主函数
int main()
{ //system(color F0); //更改输出界面颜色Mat img imread(../pic/gril_1.jpg); if (img.empty()) { cout 请确认图像文件名称是否正确 endl; return -1; } Mat gray, hist, gray2, hist2, gray3, hist3; cvtColor(img, gray, COLOR_BGR2GRAY); resize(gray, gray2, Size(), 0.5, 0.5); gray3 imread(../pic/hand.jpg, IMREAD_GRAYSCALE); const int channels[1] { 0 }; float inRanges[2] { 0,255 }; const float* ranges[1] { inRanges }; const int bins[1] { 256 }; calcHist(gray, 1, channels, Mat(), hist, 1, bins, ranges); calcHist(gray2, 1, channels, Mat(), hist2, 1, bins, ranges); calcHist(gray3, 1, channels, Mat(), hist3, 1, bins, ranges); drawHist(hist, NORM_INF, hist); drawHist(hist2, NORM_INF, hist2); drawHist(hist3, NORM_INF, hist3); //原图直方图与原图直方图的相关系数double hist_hist compareHist(hist, hist, HISTCMP_CORREL); cout apple_apple hist_hist endl; //原图直方图与缩小原图后的直方图的相关系数double hist_hist2 compareHist(hist, hist2, HISTCMP_CORREL); cout apple_apple256 hist_hist2 endl; //两幅不同图像直方图相关系数double hist_hist3 compareHist(hist, hist3, HISTCMP_CORREL); cout apple_lena hist_hist3 endl; waitKey(0); return 0;
}
3. 测试结果
