参考博客
invert()介绍
invert是opencv自带的求解逆矩阵或者伪逆矩阵的函数,具体形式如下:
double invert(InputArray src, OutputArraydst, int flags=DECOMP_LU);
- src: 输入,浮点型(32位或者64位)的M×N的矩阵,当参数3的使用方法为DECOMP_CHOLESKY DECOMP_LU DECOMP_EIG时函数功能为求逆,此时需保证M=N(参见参数flag)
- dst: 输出,与输入矩阵类型一致的N×M的矩阵
- flag:求逆方法,提供4种可选择的方法:DECOMP_CHOLESKY(基于CHOLESKY分解的方法), DECOMP_LU(基于LU分解的方法), DECOMP_EIG(基于特征值分解的方法), DECOMP_SVD(基于奇异值分解的方法)。其中,前三种方法要求输入的矩阵必须为方阵,此时计算结果为矩阵的逆;最后一种方法为对非方阵的伪逆计算,对矩阵的形状没有要求。函数接口的默认参数为DECOMP_LU方法(应该是效率较高的一种方法)
使用时出现的问题
刚开始的代码:
Mat temp(1, 2, CV_8UC1, Jw);
Mat temp1(2, 1, CV_8UC1,Jw);
cout<<temp.at<double>(0,0)<<" "<<temp.at<double>(0,1)<<endl;
invert(temp, temp1, DECOMP_SVD);
出现的问题:
发现是矩阵的类型发生了错误
于是更改代码如下:
Mat temp(1, 2, CV_64FC1, Jw);
Mat temp1(2, 1, CV_64FC1,Jw);
cout<<temp.at<double>(0,0)<<" "<<temp.at<double>(0,1)<<endl;
invert(temp, temp1, DECOMP_SVD);
问题解决。