我有一个问题,我无法找到任何解决方案。
我必须用一个已知矩阵的逆矩阵进行一些计算。
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
中,因为如果我单独打印元素,我会在逆数中获得相同的奇怪结果。
答案 0 :(得分:10)
解决: 正如彼得在我的代码中指出的那样,问题就在于此。我仍然不明白这个原因,但它就像:
我将数据CV_64F
视为浮点数,这是一个错误,需要将其视为双精度,用于写入值并读取它们。(<double>
)
但是CV_32F
可以视为float,访问权限为<float>
。
谢谢大家