有关如何在最新的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。
谢谢大家!
答案 0 :(得分:5)
您正在进行隐式类型转换,这会分散编译器的注意力。均值,权重和协方差矩阵不是整数但是加倍(你可以通过将Mat.depth()结果打印到屏幕来检查它)所以改变所有的行:
int mean1 = means.at<float>(0, 0);
代码如下:
double mean1 = means.at<double>(0, 0);
此致 拉法尔