gdal计算相关系数

参考：http://baike.baidu.com/view/172091.htm　

一、相关系数定义

      相关系数又称线性相关系数.它是衡量变量之间线性相关程度的指标。样本相关系数用r表示,总体相关系数用ρ表示,相关系数的取值范围为[-1,1]。|r|值越大，误差Q越小，变量之间的线性相关程度越高；|r|值越接近0，Q越大，变量之间的线性相关程度越低。

　　相关系数又称皮（尔生）氏积矩相关系数，说明两个现象之间相关关系密切程度的统计分析指标。

 

　　相关系数用希腊字母γ表示，γ值的范围在-1和+1之间。γ>0为 正相关，γ<0为 负相关。γ=0表示不相关；

　　γ的绝对值越大，相关程度越高。

　　两个现象之间的相关程度，一般划分为四级：

　　如两者呈正相关，r呈正值，r=1时为完全正相关；如两者呈负相关则r呈负值，而r=-1时为完全负相关。完全正相关或负相关时，所有图点都在直线回归线上；点子的分布在直线回归线上下越离散，r的绝对值越小。当例数相等时，相关系数的绝对值越接近1，相关越密切；越接近于0，相关越不密切。当r=0时，说明X和Y两个变量之间无直线关系。通常|r|大于0.8时，认为两个变量有很强的线性相关性。

 

2、gdal计算不同波段之间相关系数

（1）

利用这个公式计算相关系数时，我首先想到的是利用

CPLErr GDALRasterBand::ComputeStatistics  ( int  bApproxOK,  
  double *  pdfMin,  
  double *  pdfMax,  
  double *  pdfMean,  
  double *  pdfStdDev,  
  GDALProgressFunc  pfnProgress,  
  void *  pProgressData   
 )

这个函数计算出波段的平均值，标准差，最好只需计算一下E（XY）就好了。但是这个方法计算出来的结果不是很精确，在计算波段1与波段1之间的相关系数时都有一些误差，1.0000000000000002，可以忽略啊。

（2）

使用这个公式时，通过读取图像按部就班计算出公式中的各个参数，波段1与波段1的相关系数为1。