如何在OpenCV 2.4.2中使用CV :: EM类?

时间:2012-07-07 02:09:49

标签: opencv

有关如何在最新的OpenCV 2.4.2中使用EM算法的问题。 我曾经使用过C版本并且工作得很好,但是自从系统升级以来,似乎已经在OpenCV 2.4.2中删除了C API。

这是EM的一个简单情况:假设我有一个半径列表,被认为来自两种球,我想估计这两种类型的均值/方差。换句话说,这是一维问题。

我尝试使用新的API编写EM的C ++版本,但尚未使用它。

int nsamples = radius_list.size();
int ncluster = 2;                 //we assume a bimodal model
Mat samples = Mat::zeros(nsamples, 1, CV_32FC1);

// init data
for (int i = 0; i < radius_list.size(); ++i) {
    int value = radius_list[i];
    samples.at<float>(i, 0) = value;
}
EM em_model = EM(ncluster, EM::COV_MAT_SPHERICAL);

if (!em_model.train(samples)) {
    cerr << "error training the EM model" << endl;
    exit(-1);
}

const Mat& means = em_model.get<Mat>("means");
int mean1 = means.at<float>(0, 0);
int mean2 = means.at<float>(1, 0);
cout << "mean1 = " << mean1 << ", mean2 = " << mean2 << endl;

const vector<Mat>& covs  = em_model.get<vector<Mat> >("covs");
int scale1 = covs[0].at<float>(0, 0);
int scale2 = covs[1].at<float>(0, 0);
cout << "scale1 = " << scale1 << ", scale2 = " << scale2 << endl;

问题是:虽然if()没有抱怨,但检索到的均值和比例值 是我的机器上的垃圾值,-2147483648。

请告知如何修改代码以使其正常工作。我还在学习OpenCV中的各种C ++ API。

谢谢大家!

1 个答案:

答案 0 :(得分:5)

您正在进行隐式类型转换,这会分散编译器的注意力。均值,权重和协方差矩阵不是整数但是加倍(你可以通过将Mat.depth()结果打印到屏幕来检查它)所以改变所有的行:

int mean1 = means.at<float>(0, 0);

代码如下:

double mean1 = means.at<double>(0, 0);

此致 拉法尔

相关问题