5

我有一个问题，我找不到任何解决方案。

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

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 中检查了结果。两个矩阵都显示和读取为浮点数，它们的深度是 a 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);
}

功能打印矩阵：

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，因为如果我分别打印元素，我会在逆数中获得相同的奇怪结果。

4

1 回答 1

12

正如彼得指出的那样，问题出在我的代码中。我仍然不明白其中的原因，但它就像：

我将数据CV_64F视为浮点数，这是一个错误，需要将其视为双精度值，用于写入和读取值。（<double>）

但是CV_32F可以被视为浮动，访问将是<float>.

于 2012-06-14T12:50:57.847 回答