如何在OpenCV中计算矩阵的直方图?

时间:2015-10-07 16:51:34

标签: c++ opencv matrix histogram

我有一个由T的整数组成的矩阵,我想获得一个大小为0 to N的数据矩阵,其中每个元素都有每个数字的出现次数。

所以如果输入矩阵是

N+1

然后我的输出应该是

0 0 0
1 1 0
0 0 2

是否有一种简单的(内置)方法可以在C ++ OpenCV中实现这一点,比如MATLAB中的// occurences of 0 1 2 hist = [6, 2, 1] 函数?

谢谢!

2 个答案:

答案 0 :(得分:3)

只是提供OpenCV calcHist的替代方法,您可以使用std::vector,并在灰度值位置递增计数器:

#include <opencv2\opencv.hpp>
#include <iostream>
#include <algorithm>
#include <numeric>
using namespace std;
using namespace cv;

int main(void)
{
    // Grayscale image
    Mat1b img = imread("path_to_image", IMREAD_GRAYSCALE);

    // Compute histogram
    vector<int> hist(256, 0);
    for (int r = 0; r < img.rows; ++r)
        for (int c = 0; c < img.cols; ++c)
            ++hist[img(r, c)];

    cout << "number of 0: " << hist[0] << endl;
    cout << "number of 1: " << hist[1] << endl;
    // etc...


    // Sort according to frequency

    vector<int> indices(256);
    iota(indices.begin(), indices.end(), 0); // 0, 1, ... 255
    sort(indices.begin(), indices.end(), [&hist](int lhs, int rhs) { return hist[lhs] > hist[rhs]; });

    cout << "1st frequent number: " << indices[0] << " with N: " << hist[indices[0]] << endl;
    cout << "2nd frequent number: " << indices[1] << " with N: " << hist[indices[1]] << endl;
    // etc

    return 0;
}

更新代码段,根据频率降低对结果进行排序。

答案 1 :(得分:1)

此代码读取图像中的所有灰度值,并导致频繁出现的值(值出现的次数)。即

像素值&#39; 0&#39;如图所示,

像素值&#39; 1&#39;看来,...&amp;所以一直到256岁。

#include <opencv\cv.h>
#include <highgui\highgui.hpp>
using namespace std;
using namespace cv;
    int main()
{
cv::Mat img = cv::imread("5.jpg",0);

//for(int j=0;j<img.rows;j++) 
//{
//  for (int i=0;i<img.cols;i++)
//  {
//    int a;
//  a=img.at<uchar>(j,i);
//     cout<<a<<endl; 
//  }
//}
vector<int> values_rgb;
for(int i=0; i<20; i++)
{
    for(int j=0; j<20; j++)
    {
        int value_rgb = img.at<uchar>(i,j);
        values_rgb.push_back(value_rgb);
        //cout << "(" << i << "," << j << ") : " << value_rgb <<endl;
    }
}
// Sorting of values in ascending order 
vector<int> counter_rg_values;
for(int l=0; l<256; l++)

{
    for(int k=0; k<values_rgb.size(); k++)
    {
        if(values_rgb.at(k) == l)
        {
            counter_rg_values.push_back(l);
        }
    }
}
//for(int m=0;m<counter_rg_values.size();m++)
//cout<<m<<" "<< counter_rg_values[m] <<endl;
int m=0;
for(int n=0;n<256;n++)
{
    int c=0;
    for(int q=0;q<counter_rg_values.size();q++)
    {
        if(n==counter_rg_values[q])
        {
            //int c;
        c++;
        m++;
        }
    }

    cout<<n<<"= "<< c<<endl;
}
cout<<"Total number of elements "<< m<<endl;

cv::imshow("After",img);
waitKey(0);
}

关于https://stackoverflow.com/a/32902450/3853072

的上述讨论的链接