逆矩阵OpenCV。 Matrix.inv()无法正常工作

时间:2012-06-13 12:45:45

标签: c++ opencv matrix-inverse

我有一个问题,我无法找到任何解决方案。

我必须用一个已知矩阵的逆矩阵进行一些计算。

Matrix homography=

1.1688, 0.23, 62.2,

-0.013,1.225, -6.29,

0, 0, 1, 

然后:

Mat homoInv=homography.inv();

矩阵的内容是:

1.81381e-29, 15.1628, -7.57361e+17, 

0, -0, 0, 

5.4561e-33, -2.40123e+34, -1.38198e-05

那当然是错的,因为我已经在Matlab中检查了结果。两个矩阵都显示并作为浮点读取,其深度为64FC1

有没有人知道可以做些什么?

全部谢谢

更多代码:

int main(int argc, char ** argv )
{ 
  Mat homogra(3,3,CV_64FC1);
  Mat coord(3,1,CV_64FC1);
  Mat result(target.size(),CV_8UC1);
  homogra.at<float>(0,0)=1.1688;
  homogra.at<float>(0,1)=0.23;
  homogra.at<float>(0,2)=(-62.20);
  homogra.at<float>(1,0)=(-0.013);
  homogra.at<float>(1,1)=1.225;
  homogra.at<float>(1,2)=-6.29;
  homogra.at<float>(2,0)=0;
  homogra.at<float>(2,1)=0;
  homogra.at<float>(2,2)=1;
  printMatrix(homogra);

  Mat inverse=homogra.inv();
  printMatrix(inverse);
}

function printMatrix:

void printMatrix(Mat M){
cout<<"Tipo de matriz:"<<M.type()<<endl;
 // dont print empty matrices
  if (M.empty()){
    cout << "---" << endl;
    return;
  }
  // loop through columns and rows of the matrix
  for(int i=0; i < M.rows; i++){
      for(int j=0; j < M.cols ; j++){
      cout << M.at<float>(i,j) << ", "<<endl;
      }
    cout<<"Change\n"<<endl;
}
  }

但错误不在printMatrix中,因为如果我单独打印元素,我会在逆数中获得相同的奇怪结果。

1 个答案:

答案 0 :(得分:10)

解决: 正如彼得在我的代码中指出的那样,问题就在于此。我仍然不明白这个原因,但它就像:

我将数据CV_64F视为浮点数,这是一个错误,需要将其视为双精度,用于写入值并读取它们。(<double>

但是CV_32F可以视为float,访问权限为<float>

谢谢大家