2

有一个关于如何在最新的 OpenCV 2.4.2 中使用 EM 算法的问题。我曾经使用 C 版本，它工作得非常好，但自从系统升级以来，似乎 C API 已在 OpenCV 2.4.2 中删除。

这是 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。

谢谢你们！

4

1 回答 1

5

您进行隐式类型转换会分散编译器的注意力。均值、权重和协方差矩阵不是整数，而是双精度数（您可以通过将 Mat.depth() 结果打印到屏幕上来检查），因此将所有行从：

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

喜欢代码：

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

问候， 拉法尔

于 2012-11-19T11:51:30.930 回答