8. 圖像直方圖

image.png

image.png
#include <opencv.hpp>
#include <iostream>

using namespace cv;
using namespace std;

int main(int argc, char** argv) {

    Mat src = imread("D:/HUANGHAI_WORK/tools/opencv4/opencv/sources/samples/data/home.jpg");


    namedWindow("input", WINDOW_AUTOSIZE);
    imshow("input", src);
    

    vector<Mat> mv;
    split(src, mv);

    //計(jì)算直方圖
    int histSize = 256;
    float range[] = { 0,255 };
    const float* histRange = { range };
    Mat b_hist, g_hist, r_hist;
    calcHist(&mv[0], 1, 0, Mat(), b_hist, 1, &histSize, &histRange, true, false);
    calcHist(&mv[1], 1, 0, Mat(), g_hist, 1, &histSize, &histRange, true, false);
    calcHist(&mv[2], 1, 0, Mat(), r_hist, 1, &histSize, &histRange, true, false);

    Mat result = Mat::zeros(Size(600, 400), CV_8UC3);
    int margin = 50;
    int nm = result.rows;
    normalize(b_hist, b_hist, 0, nm,NORM_MINMAX,-1,Mat());
    normalize(g_hist, g_hist, 0, nm, NORM_MINMAX, -1, Mat());
    normalize(r_hist, r_hist, 0, nm, NORM_MINMAX, -1, Mat());
    float step = 500.0 / 256.0;
    for (int i = 0; i < 255; i++)
    {
        float f1 = b_hist.at<float>(i, 0);
        float f2 = b_hist.at<float>(i+1, 0);
        line(result, Point(step*i+50, 50 + nm - b_hist.at<float>(i, 0)),Point(step*i+50+1,50+nm- b_hist.at<float>(i + 1, 0)),Scalar(255,0,0),1,8,0);
        line(result, Point(step*i+50, 50 + nm - g_hist.at<float>(i, 0)), Point(step*i+50 + 1, 50 + nm - g_hist.at<float>(i + 1, 0)), Scalar(0, 255, 0), 1, 8, 0);
        line(result, Point(step*i+50, 50 + nm - r_hist.at<float>(i, 0)), Point(step*i+50 + 1, 50 + nm - r_hist.at<float>(i + 1, 0)), Scalar(0, 0, 255), 1, 8, 0);

    }
    imshow("result", result);


    waitKey(0);
    destroyAllWindows();
    return 0;
}
image.png

image.png
//直方圖均衡化
void test2() {
    Mat src = imread("D:/HUANGHAI_WORK/tools/opencv4/opencv/sources/samples/data/home.jpg");


    namedWindow("input", WINDOW_AUTOSIZE);
    imshow("input", src);

    Mat gray, dst;
    cvtColor(src, gray, COLOR_BGR2GRAY);
    imshow("gray", gray);
    equalizeHist(gray, dst);
    imshow("eq", dst);

    //計(jì)算直方圖
    int histSize = 256;
    float range[] = { 0,255 };
    const float* histRange = { range };
    Mat b_hist, g_hist;
    calcHist(&gray, 1, 0, Mat(), b_hist, 1, &histSize, &histRange, true, false);
    calcHist(&dst, 1, 0, Mat(), g_hist, 1, &histSize, &histRange, true, false);


    Mat result = Mat::zeros(Size(600, 400), CV_8UC3);
    int margin = 50;
    int nm = result.rows;
    normalize(b_hist, b_hist, 0, nm, NORM_MINMAX, -1, Mat());
    normalize(g_hist, g_hist, 0, nm, NORM_MINMAX, -1, Mat());
    float step = 500.0 / 256.0;
    for (int i = 0; i < 255; i++)
    {
        float f1 = b_hist.at<float>(i, 0);
        float f2 = b_hist.at<float>(i + 1, 0);
        line(result, Point(step*i + 50, 50 + nm - b_hist.at<float>(i, 0)), Point(step*i + 50 + 1, 50 + nm - b_hist.at<float>(i + 1, 0)), Scalar(0, 0, 255), 1, 8, 0);
        line(result, Point(step*i + 50, 50 + nm - g_hist.at<float>(i, 0)), Point(step*i + 50 + 1, 50 + nm - g_hist.at<float>(i + 1, 0)), Scalar(0, 255, 255), 1, 8, 0);

    }
    imshow("result", result);
}
image.png

image.png
void test3() {
    Mat src1 = imread("D:/HUANGHAI_WORK/tools/opencv4/opencv/sources/samples/data/aloeL.jpg");
    Mat src2 = imread("D:/HUANGHAI_WORK/tools/opencv4/opencv/sources/samples/data/aloeR.jpg");
    
    //計(jì)算直方圖
    int histSize[] = {256,256,256};
    float range[] = { 0,255 };
    float c1[] = { 0,255 };
    float c2[] = { 0,255 };
    float c3[] = { 0,255 };
    const float* histRange[] = { c1,c2,c3 };
    int channels[] = { 0,1,2 };
    Mat hist1, hist2;
    calcHist(&src1, 1, channels, Mat(), hist1, 3, histSize, histRange, true, false);
    calcHist(&src2, 1, channels, Mat(), hist2, 3, histSize, histRange, true, false);

    //歸一化
    normalize(hist1, hist1, 0, 1.0, NORM_MINMAX, -1, Mat());
    normalize(hist2, hist2, 0, 1.0, NORM_MINMAX, -1, Mat());

    //巴氏距離比較
    double b12 = compareHist(hist1, hist2, HISTCMP_BHATTACHARYYA);
    double b11 = compareHist(hist1, hist1, HISTCMP_BHATTACHARYYA);
    printf("巴氏距離:%.2f,%.2f \n", b12, b11);

    //相關(guān)性比較
    double c12 = compareHist(hist1, hist2, HISTCMP_CORREL);
    double c11 = compareHist(hist1, hist1, HISTCMP_CORREL);

    printf("相關(guān)性:%.2f,%.2f \n",c12,c11);
}
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末凉袱,一起剝皮案震驚了整個(gè)濱河市微王,隨后出現(xiàn)的幾起案子蠢熄,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 223,126評(píng)論 6 520
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡却妨,警方通過(guò)查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 95,421評(píng)論 3 400
  • 文/潘曉璐 我一進(jìn)店門(mén)括眠,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)彪标,“玉大人,你說(shuō)我怎么就攤上這事掷豺±萄蹋” “怎么了?”我有些...
    開(kāi)封第一講書(shū)人閱讀 169,941評(píng)論 0 366
  • 文/不壞的土叔 我叫張陵当船,是天一觀的道長(zhǎng)题画。 經(jīng)常有香客問(wèn)我,道長(zhǎng)德频,這世上最難降的妖魔是什么苍息? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 60,294評(píng)論 1 300
  • 正文 為了忘掉前任,我火速辦了婚禮壹置,結(jié)果婚禮上竞思,老公的妹妹穿的比我還像新娘。我一直安慰自己钞护,他們只是感情好盖喷,可當(dāng)我...
    茶點(diǎn)故事閱讀 69,295評(píng)論 6 398
  • 文/花漫 我一把揭開(kāi)白布。 她就那樣靜靜地躺著难咕,像睡著了一般课梳。 火紅的嫁衣襯著肌膚如雪距辆。 梳的紋絲不亂的頭發(fā)上,一...
    開(kāi)封第一講書(shū)人閱讀 52,874評(píng)論 1 314
  • 那天暮刃,我揣著相機(jī)與錄音跨算,去河邊找鬼。 笑死椭懊,一個(gè)胖子當(dāng)著我的面吹牛诸蚕,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播灾搏,決...
    沈念sama閱讀 41,285評(píng)論 3 424
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼挫望,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼立润!你這毒婦竟也來(lái)了狂窑?” 一聲冷哼從身側(cè)響起,我...
    開(kāi)封第一講書(shū)人閱讀 40,249評(píng)論 0 277
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤桑腮,失蹤者是張志新(化名)和其女友劉穎泉哈,沒(méi)想到半個(gè)月后,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體破讨,經(jīng)...
    沈念sama閱讀 46,760評(píng)論 1 321
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡丛晦,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,840評(píng)論 3 343
  • 正文 我和宋清朗相戀三年,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了提陶。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片烫沙。...
    茶點(diǎn)故事閱讀 40,973評(píng)論 1 354
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡,死狀恐怖隙笆,靈堂內(nèi)的尸體忽然破棺而出锌蓄,到底是詐尸還是另有隱情,我是刑警寧澤撑柔,帶...
    沈念sama閱讀 36,631評(píng)論 5 351
  • 正文 年R本政府宣布瘸爽,位于F島的核電站,受9級(jí)特大地震影響铅忿,放射性物質(zhì)發(fā)生泄漏剪决。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 42,315評(píng)論 3 336
  • 文/蒙蒙 一檀训、第九天 我趴在偏房一處隱蔽的房頂上張望柑潦。 院中可真熱鬧,春花似錦峻凫、人聲如沸妒茬。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 32,797評(píng)論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)乍钻。三九已至肛循,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間银择,已是汗流浹背多糠。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 33,926評(píng)論 1 275
  • 我被黑心中介騙來(lái)泰國(guó)打工, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留浩考,地道東北人夹孔。 一個(gè)月前我還...
    沈念sama閱讀 49,431評(píng)論 3 379
  • 正文 我出身青樓,卻偏偏與公主長(zhǎng)得像析孽,于是被迫代替她去往敵國(guó)和親搭伤。 傳聞我的和親對(duì)象是個(gè)殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,982評(píng)論 2 361

推薦閱讀更多精彩內(nèi)容