OpenCV Mat矩阵乘法结果与手动计算不符，求原因及代码排查

给定信息 #

4x4变换矩阵 #

1.000, 0.000, -0.004, 0.784,
-0.000, 1.000, 0.001, -0.205,
0.004, -0.001, 1.000, -1.435,
0.000, 0.000, 0.000, 1.000

三维向量 #

507.347, 359.172, 229.584

执行的OpenCV代码 #

double data[] = { 507.347, 359.172, 229.584, 1 };

cv::Mat s(4, 1, CV_64F, data);
cv::Mat d = transformMatrix * s;

结果差异 #

• OpenCV计算结果：result = matrix * (507.347,359.172,229.584) = (507.394,359.158,229.605)
• 手动计算结果：result = 507.213, 359.197, 230.538

原因排查与分析 #

代码逻辑本身无问题 #

你写的矩阵乘法代码逻辑是对的：构造4x1的齐次列向量，执行矩阵×列向量的运算，符合齐次坐标下的变换规则，语法层面没有错误。

核心问题：变换矩阵的构造错误 #

从数值对比来看，OpenCV的计算结果和正确的矩阵乘列向量结果（即你的手动计算值）偏差很大，大概率是代码中transformMatrix的实际内容和你给出的矩阵不一致：

• 存储顺序错误：OpenCV的Mat默认是行优先存储，如果你在初始化transformMatrix时误将矩阵的列作为行来填充，或者反之，会导致矩阵被转置，运算结果自然完全不同。
• 元素输入错误：初始化transformMatrix时可能写错了元素值（比如符号、小数点位置），或者行列数设置错误（比如把4x4矩阵搞成了非预期尺寸）。

验证方法 #

直接在代码中输出transformMatrix的所有元素，和你给出的矩阵逐一对比，就能快速定位问题。比如可以用cv::cout << transformMatrix << std::endl;直接输出矩阵内容。

补充：手动计算的正确性验证 #

我们以第一行结果为例，用你给出的矩阵做正确的矩阵×列向量计算：
1.000*507.347 + 0.000*359.172 + (-0.004)*229.584 + 0.784*1 = 507.347 - 0.918336 + 0.784 ≈ 507.213，和你的手动计算结果一致，说明手动计算是正确的，问题确实出在代码中的transformMatrix构造环节。

内容的提问来源于stack exchange，提问作者EunJong Park