Question

我有一个可能适合初学者的问题。

我想在c ++中将20x2矩阵乘以2x2矩阵。

我尝试使用openCV，但是我收到错误

cvarrToMat中的错误参数（未知数组类型）

这是我在openCV中使用的代码，以便检查问题，如果它是我的代码或openCV中的问题，但它仍然无法正常工作，我可以编译出错，但是当我测试代码时得到问题“cvarrToMat中的错误参数（未知数组类型）”

#include <stdio.h>
#include <stdlib.h>
//#include "/usr/include/opencv/cv.h"
#include <cv.h>
#include <cvaux.h>
#include <highgui.h>
#include <math.h>
#include <iostream>

  int main()
{


double a[] = {1, 2, 3, 4};
CvMat Ma;
cvInitMatHeader(&Ma, 2, 2, CV_32FC1, a);


double b[] ={0, -1, 1, 0};

CvMat Mb;
cvInitMatHeader(&Mb, 2, 2, CV_32FC1, b);

CvMat Mc;
CvMat Mc1;
cvMatMul(&Ma, &Mb, &Mc);

return 0;
}

Answer 1

将代码与example in the OpenCV docs进行比较，似乎您忘记初始化输出矩阵Mc：

double a[] = { 1, 2, 3, 4,
               5, 6, 7, 8,
               9, 10, 11, 12 };

double b[] = { 1, 5, 9,
               2, 6, 10,
               3, 7, 11,
               4, 8, 12 };

double c[9];
CvMat Ma, Mb, Mc ;

cvInitMatHeader(&Ma, 3, 4, CV_64FC1, a);
cvInitMatHeader(&Mb, 4, 3, CV_64FC1, b);
cvInitMatHeader(&Mc, 3, 3, CV_64FC1, c);

cvMatMulAdd(&Ma, &Mb, 0, &Mc);
// the c array now contains the product of a (3x4) and b (4x3)

根据文档，cvMatMul(&Ma, &Mb, &Mc)与cvMatMulAdd(&Ma, &Mb, 0, &Mc)相同。

Answer 2

好。这个问题的答案实际上取决于一些事情。您已经说过，您知道如何手动执行此操作，因此在代码中执行此操作将取决于您如何表示矩阵。现在，如果这是一次性的事情，你只需要答案，我会建议像MATLAB这样的语言。如果它是较大程序的一部分，并且您将进行大量需要高效的矩阵乘法，我建议使用高质量的优化库，如boost::ublas。

如果它是一次性的事情并且你真的想用C ++做它并且你真的不想/知道如何使用像ublas这样的第三方库，那么（非优化的）矩阵乘法就像以下内容：

template<typename T>
struct matrix2d
{
private:
    std::vector<std::vector<T>> data;
    size_t _rows, _columns;
public:
    matrix2d(size_t rows, size_t columns)
        :_rows(rows)
        ,_columns(columns)
    {
        data.resize(_rows, std::vector<T>(_columns));
    }

    size_t rows() const { return _rows; } 
    size_t columns() const { return _columns; } 

    T& operator()(size_t row, size_t column)
    {
        return data[row][column];
    }

    const T& operator()(size_t row, size_t column) const
    {
        return data[row][column];
    }
};

template<typename T>
void mmult(const matrix2d<T>& m1, const matrix2d<T>&m2, matrix2d<T>& result)
{
    for (size_t r = 0 ; r<m1.rows() ; ++r)
        for (size_t c = 0; c<m2.columns() ; ++c)
            for (size_t n = 0; n<m1.columns() ; ++n)
                result(r, c) = m1(r, n) * m2(n, c);
}


int main()
{

    matrix2d<double> m1(20, 2);
    matrix2d<double> m2(2, 2);
    matrix2d<double> result(m1.rows(), m2.columns());
    mmult(m1, m2, result);
}

Answer 3

可能您应该发布您正在呼叫的功能的原型，以及您的矩阵声明和您的通话。我不认为每个人都熟悉openCV。

矩阵乘法

3 个答案: