OpenCV：包含Vec3d的矩阵和包含双打的矩阵的矩阵乘法

Question

我正在使用OpenCV一段时间，现在我需要这种类型的乘法：

定义一个矩阵T，其中包含类型为Vec3d¹的元素。矩阵T的大小为：M X N。矩阵T必须与VectorPhi相乘，VectorPhi的大小为：N X 1，包含双s作为值。结果的每个元素必须是两个矩阵的矩阵相乘的结果。

我不想做组件式乘法，而是“实”矩阵乘法，例如，将t2的第一个元素与矩阵J的第一个元素相乘，然后将矩阵t3的第二个元素与矩阵J的第二个元素相乘。这样做直到你完成第一行的T，然后总结结果。结果是一个m1。

例如，如果T是一个3x2矩阵，而Phi是一个2x1矩阵，那么对于结果的第一个值，计算应该是T_11*Phi_11 T_12*Phi_21。目前，我使用两种速度较慢的for循环：

for (int i = 0; i<M; ++i){
    cv::Mat summedResult = cv::Mat(3, 1, CV_64F, double(0));
    for (uint32 j = 0; j<N; ++j){
        summedResult = summedResult + 
        (cv::Mat(mMatrixT.at<cv::Vec3d>(i, j)) * mMatrixPhi.at<double>(j));
    }
    // The result matrix contains values of type Vec3d again
    mResultMatrix.at<cv::Vec3d>(i) = cv::Vec3d(summedResult);
}

更一般地说：在OpenCV中，是否有可能有效地将包含Vec3ds和双代码的矩阵相乘？

^{1。包含双精度的三维向量}

^{2。坐标：1,1}

^{坐标： 1,2}




Answer 1

我仍然不知道你期望得到什么样的结果，但也许可以试试这个：

假设你有一个Vec3d的MxN矩阵和一个double类型的Nx1矩阵，你的结果将是一个Vec3d类型的Mx1矩阵：

for (int i = 0; i<M; ++i)
    {
        cv::Vec3d summedResult; // here this must be a Vec3d instead of a matrix, if I assume the right result expection
        for (uint32 j = 0; j<N; ++j)
        {
            summedResult = summedResult + (mMatrixT.at<cv::Vec3d>(i, j) * mMatrixPhi.at<double>(j));
        }
        // The result matrix contains values of type Vec3d again
        mResultMatrix.at<cv::Vec3d>(i) = summedResult;
    }

编辑：

啊，对不起，我没有读到最后，你提供的代码工作，但太慢了。。。嗯，我认为这是没有优化的，因为数学上这是没有定义的。你可以做的是将你的Vec3d矩阵转换成Mx（3*N）矩阵，将你的Nx1矩阵转换成（3*N）x1矩阵（在下一个值之前是相同值的3倍），并直接使用OpenCV矩阵产品。但这可能不会更快，因为两个矩阵的大小都是3*

编辑：将是不同的结果，因为每个元素将是Vec3d元素的总和。。。