在 C++ 中实现 matlab hist()答案

【问题标题】：implementing matlab hist() in c++在 C++ 中实现 matlab hist()
【发布时间】：2015-06-02 11:08:12
【问题描述】：

这是我要实现的代码行

kb = [];
    for k = 1:length(nRef)
        for n=1:length(dCmpInd{k})
            x = [centroid(nRef{k}, 1), centroid(dCmpInd{k}(n),1)];
            y = [centroid(nRef{k}, 2), centroid(dCmpInd{k}(n),2)];

            [x,ind] = sort(x);
            y = y(ind);
            kb = [kb diff(y) / diff(x)];
        end
    end

    theta = (atan(kb));

    [N, X] = hist(abs(theta),2);

这是我的 C++ 代码：

 std::vector<double> kb;
    std::vector<double> theta;
    for (int k = 0; k < nRef.size(); k++)
    {
        for (int n = 0; n < dCmpInd[k].size(); n++)
        {
            double x1 = centroids[nRef[k]].m_X; double x2 = centroids[dCmpInd[k][n]].m_X;
            double y1 = centroids[nRef[k]].m_Y; double y2 = centroids[dCmpInd[k][n]].m_Y;
            if (x1 <x2)
            {
                double tempdiff = (y2-y1)/(x2-x1);              
                kb.push_back(tempdiff);
                theta.push_back(abs(atan(tempdiff)));
            }
            else
            {
                double tempdiff = (y1-y2)/(x1-x2);
                kb.push_back(tempdiff);
                theta.push_back(abs(atan(tempdiff)));
            }
        }
    }

有没有快速实现的方法：

[N,X] = hist(theta,2);

我也可以使用 openCV 2.4.10，但 calcHist() 并不完全相同，我需要创建 2 个 bin。

我的输入是一维数组：

0.00598881028540019 1.56120677124307    0.00598881028540019 0.00669537049433832 1.37723800334516    1.37723800334516    1.36424594043624    1.56120677124307    0.0152220988707370

输出是：

X= 0.394793300524817    1.17240228100365
 N = 4 5

【问题讨论】：

使用 OpenCV 的功能有什么问题？无论如何，实现自定义直方图类/函数应该是不必要的（除非你想要一些练习），有足够的示例可用。
@Baiz 有一个用于 matlab 的功能，您可以在其中将一维数组排序到所需的箱数中。在我的情况下= 2。我无法在 openCV 中实现这一点。也许我只是不知道怎么做，但我已经尝试过了。
你是如何调用 OpenCV 的 calcHist 函数的？查看documentation，参数 histSize 是您需要指定的。因为你真正想要的只是一个带有两个 bin 的直方图。

标签： c++ matlab opencv image-processing

【解决方案1】：

calcHist 工作正常。尝试以下方式选择 histSize 和 range。希望能帮助到你！。

float ary[9] = { 0.00598881028540019, 1.56120677124307,0.00598881028540019, 0.00669537049433832,1.37723800334516, 1.37723800334516, 1.36424594043624,1.56120677124307, 0.0152220988707370 };
cv::Mat srcMat = cv::Mat(1, 9, CV_32FC1, ary);
int histSize = 2;
float range[] = { 0.0, 2.0 };
const float* histRange = { range };
bool uniform = true;
bool accumulate = false;
cv::Mat histOut;
cv::calcHist(&srcMat, 1, 0, Mat(),histOut, 1, &histSize, &histRange, uniform, accumulate);
cout << histOut << endl;

【讨论】：

好的，我看到你指出了垃圾箱的大小，但是我从哪里得到 X= 0.394793300524817 1.17240228100365？
我不认为你可以直接使用 calcHist 或任何函数来获得它。这些不过是垃圾箱中心。您可以使用这个 matlab 脚本来重现结果。 binwidth = (maxy - miny) ./ x; xx = miny + binwidth*(0:x); xx(length(xx)) = maxy; x = xx(1:length(xx)-1) + binwidth/2;

【解决方案2】：

虽然这不是一个完整的实现，但这是我所做的

HistogramResults Utilities::Hist(std::vector<double> vector, std::vector<double> x)
{
    if (vector.empty())
    {
        return HistogramResults();
    }
    if(x.empty())
    {
        x.push_back(10);
    }
    std::vector<double> N(x);
    std::vector<double>::iterator iter = std::min_element(std::begin(vector),std::end(vector));
    double minY = *iter;
    iter = std::max_element(std::begin(vector),std::end(vector));
    double maxY = *iter;
    std::vector<double> xx;
    if (x.size() == 1)
    {
        if (minY == maxY)
        {
            minY = minY - floor(x[0] /2) - 0.5;
            maxY = maxY - ceil(x[0] /2) - 0.5;
        }
        double binwidth = (maxY - minY)/x[0];

        xx.push_back(-1*std::numeric_limits<double>::infinity());
        for (int i = 0; i <= x[0]; i++)
        {
            xx.push_back( minY +  binwidth * i);
        }
        xx[x[0]+1] = maxY;
        xx.push_back(std::numeric_limits<double>::infinity());
        x.clear();

        for (int i = 1; i < xx.size()-2; i++)
        {
            x.push_back(xx[i] + binwidth/2);
        }       
    }
    else
    {
    }

    std::vector<int> nn =HistC(vector,xx);
    nn[1] = nn[1] +nn[0];
    nn[nn.size()-1] == nn[nn.size()-1] + nn[nn.size()-2];
    nn.pop_back();
    nn.erase (nn.begin());
    return HistogramResults(nn,x);
}

std::vector<int> Utilities::HistC(std::vector<double> y, std::vector<double> edges)
{
    std::vector<int> bins(edges.size()-1);
    for (int i = 0; i < y.size(); i++)
    {
        for (int j = 0; j < edges.size()-1; j++)
        {
            if( (edges[j] < y[i]) && (y[i] <= edges[j+1]))
            {
                bins[j]= bins[j]+1;
            }
        }
    }
    return bins;
}

这里是 HistorgramResults 类

#include "HistogramResults.h"


HistogramResults::HistogramResults(std::vector<int> numberOfElementN, std::vector<double> averageX)
{
    m_numberOfElementN = numberOfElementN;
    m_averageX = averageX;
}

HistogramResults::HistogramResults(void)
{
}

HistogramResults::~HistogramResults(void)
{
}

【讨论】：