相关分析、偏相关+t值检验、复相关+F值检验

# -*- coding: utf-8 -*-     
import arcpy
from arcpy import env
from arcpy.sa import *

#读取数据
path=“K:/Data3“
workspace1=path+“/growth“
arcpy.env.workspace=unicode（workspace1“utf8“）
YS=arcpy.ListRasters（“2*““tif“）
XS=arcpy.ListRasters（“3*““tif“）
ZS=arcpy.ListRasters（“2*““tif“）
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension（“Spatial“）

#计算Y值平均值
yMean=CellStatistics（YS “MEAN“ “DATA“）
#计算X值平均值
xMean=CellStatistics（XS “MEAN“ “DATA“）
#计算Z值平均值
zMean=CellStatistics（ZS “MEAN“ “DATA“）
#计算r1r2r3分别代表xyxzyz相关系数
#数据初始化
n=len（XS）
r1xyr1xr1yr2xzr2xr2z=0.00.00.00.00.00.0
for i in range（0n）:
    #计算xy参数
    r1xy=（XS[i]-xMean）*（YS[i]-