Question

我正在尝试使用opencv EM算法库。

我首先尝试了一些关于EM算法的例子，比如this和this，但我在这些代码中遇到了同样的问题。在EM培训步骤中花了很多时间。

这是我从this website进行图像分割的测试代码：

#include <opencv2/opencv.hpp>
#include <opencv2/legacy/legacy.hpp>
using namespace std;
using namespace cv;
Mat asSamplesVectors(Mat& img);
/** @function main */
int main(int argc, char** argv)
{
    Mat image; int no_of_clusters = 2;
    image = imread("images/flower.jpg");
    Mat samples = asSamplesVectors(image);

    cout << "Starting EM training" << endl;
    EM em(no_of_clusters);
    em.train(samples);
    cout << "Finished training EM" << endl;

    vector<Mat> segmented;
    for (int i = 0; i < no_of_clusters; i++)
        segmented.push_back(Mat::zeros(image.rows, image.cols, CV_8UC3));

    int index = 0;
    for (int y = 0; y < image.rows; y++) {
        for (int x = 0; x < image.cols; x++) {
            int result = em.predict(samples.row(index++))[1];
            segmented[result].at<Point3i>(y, x, 0) = image.at<Point3i>(y, x, 0);
        }
    }
    imshow("result", samples);
    return(0);
}

Mat asSamplesVectors(Mat& img) {
    Mat float_img;
    img.convertTo(float_img, CV_32F);

    Mat samples(img.rows * img.cols, 3, CV_32FC1);

    /* Flatten  */
    int index = 0;
    for (int y = 0; y < img.rows; y++) {
        Vec3f* row = float_img.ptr<Vec3f>(y);
        for (int x = 0; x < img.cols; x++)
            samples.at<Vec3f>(index++, 0) = row[x];
    }
    return samples;
}

我在VS2013中测试，opencv 2.4.9和opencv 2.4.11。

这是我的测试image。

在此代码中，它显示错误message。

错误步骤为segmented[result].at<Point3i>(y, x, 0) = image.at<Point3i>(y, x, 0);

我的第一个问题是：为什么EM训练如此缓慢。通常是吗？

第二个问题是：我的代码错误有什么问题，是否有更好的例子让我知道EM的用法？

Answer 1

我找到了第一个问题的答案，我的测试图像太大而无法训练。训练图像的大小需要缩小。

我使用此article来修改像素，它可以正常工作。

Opencv：EM算法训练缓慢

1 个答案: