opencv利用calcCovarMatrix函数求协方差矩阵存在的问题

来源：互联网发布：c语言怎么判断奇偶编辑：程序博客网时间：2024/04/28 23:57

今天用到opencv求解矩阵的协方差矩阵，无意中发现了一个问题，写出来供大家参考。

#include <opencv2/opencv.hpp>#include <iostream>using namespace std;using namespace cv;void displayMat(Mat& mat){for (int i = 0;i<mat.rows;i++)  {  for (int j = 0;j<mat.cols ;j++)  {  printf("%f  ",mat.at<double>(i,j));  }  printf("\n");  }}int main(){double data[5][4] = { 1, 2, 3 , 4 ,        5, 6, 7 , 8 ,        9, 0, 1 , 3 ,        2, 4, 6 , 8 ,        1, 3, 5 , 9      };  Mat mat(5,4,CV_64FC1) ;  for (int i =0;i<mat.rows;i++)  {  for (int j = 0;j<mat.cols ;j++)  {  mat.at<double>(i,j) = data[i][j];  }}Mat covMat;  Mat meanMat;  calcCovarMatrix(mat,covMat,meanMat,CV_COVAR_NORMAL|CV_COVAR_ROWS);cout<<"cov:"<<endl;//打印矩阵displayMat(covMat);cout<<"mean:"<<endl;displayMat(meanMat);return 0;}

实验结果：

cov:47.200000  -11.000000  -15.200000  -19.200000-11.000000  20.000000  21.000000  19.000000-15.200000  21.000000  23.200000  23.200000-19.200000  19.000000  23.200000  29.200000mean:3.600000  3.000000  4.400000  6.400000

只看这个结果你肯定不会发现任何问题，下面让我们看看matlab中的结果

 data=[1, 2, 3 , 4;       5, 6, 7 , 8;       9, 0, 1 , 3;        2, 4, 6 , 8;       1, 3, 5 , 9 ];>> cov(data)ans =   11.8000   -2.7500   -3.8000   -4.8000   -2.7500    5.0000    5.2500    4.7500   -3.8000    5.2500    5.8000    5.8000   -4.8000    4.7500    5.8000    7.3000>> cov(data)*4ans =   47.2000  -11.0000  -15.2000  -19.2000  -11.0000   20.0000   21.0000   19.0000  -15.2000   21.0000   23.2000   23.2000  -19.2000   19.0000   23.2000   29.2000

这次你肯定会发现问题了，是的，opencv中求得的协方差矩阵要比matlab中的大4倍！

经过我多次实验的出结论：opencv中求得的协方差矩阵比matlab中的大（n-1）倍！（n是矩阵的行数）

不知道是什么原因，有知道的大神可以说一下。

0 0