opencv back projection

最新推荐文章于 2023-06-30 15:09:41 发布

littletomatodonkey

最新推荐文章于 2023-06-30 15:09:41 发布

阅读量700

点赞数

分类专栏： opencv 文章标签： opencv

本文链接：https://blog.csdn.net/u012526003/article/details/52781305

版权

opencv 专栏收录该内容

5 篇文章 0 订阅

订阅专栏

简介

反向投影是通过计算图像的直方图，得到图像的概率密度分布，并作变化，得到新的图像，在新的图像中，亮度越高，待变此区间内的像素值的概率密度最大，可以用于camshift算法等。

比如说：矩阵M = {1,2,3,4;4,6,7,8;9,10,11,12;9,10,14,15}，直方图的统计范围是[1,4],[5,8],[9,12],[13,16],则各个区间内的值的个数分别为{5,3,6,2}，得到最终的反向投影矩阵为{5,5,5,5;5,3,3,3;6,6,6,6;6,6,2,2}，可以通过归一化，对反向投影矩阵进行可视化，最终6所在的位置的灰度值最大，即最大的灰度值所在的地方对应的直方图的概率密度最大。

代码

#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/highgui/highgui.hpp"

#include <iostream>

using namespace cv;
using namespace std;

/// Global Variables
Mat src; Mat hsv; Mat hue;
int bins = 25;

/// Function Headers
void Hist_and_Backproj(int, void*);

/** @function main */
int main(int argc, char** argv)
{
    /// Read the image
    src = imread("Back_Projection_Theory2.jpg", 1);
    /// Transform it to HSV
    cvtColor(src, hsv, CV_BGR2HSV);

    /// Use only the Hue value
    hue.create(hsv.size(), hsv.depth());
    int ch[] = { 0, 0 };
    mixChannels(&hsv, 1, &hue, 1, ch, 1);

    /// Create Trackbar to enter the number of bins
    char* window_image = "Source image";
    namedWindow(window_image, CV_WINDOW_AUTOSIZE);
    createTrackbar("* Hue  bins: ", window_image, &bins, 180, Hist_and_Backproj);
    Hist_and_Backproj(0, 0);

    /// Show the image
    imshow(window_image, src);

    /// Wait until user exits the program
    waitKey(0);
    return 0;
}


/**
* @function Hist_and_Backproj
* @brief Callback to Trackbar
*/
void Hist_and_Backproj(int, void*)
{
    MatND hist;
    int histSize = MAX(bins, 2);
    float hue_range[] = { 0, 180 };
    const float* ranges = { hue_range };

    /// Get the Histogram and normalize it
    calcHist(&hue, 1, 0, Mat(), hist, 1, &histSize, &ranges, true, false);
    normalize(hist, hist, 0, 255, NORM_MINMAX, -1, Mat());

    /// Get Backprojection
    MatND backproj;
    calcBackProject(&hue, 1, 0, hist, backproj, &ranges, 1, true);

    /// Draw the backproj
    imshow("BackProj", backproj);

    /// Draw the histogram
    int w = 400; int h = 400;
    int bin_w = cvRound((double)w / histSize);
    Mat histImg = Mat::zeros(w, h, CV_8UC3);

    for (int i = 0; i < bins; i++)
    {
        rectangle(histImg, Point(i*bin_w, h), Point((i + 1)*bin_w, h - cvRound(hist.at<float>(i)*h / 255.0)), Scalar(0, 0, 255), -1);
    }

    imshow("Histogram", histImg);
}