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

#讀取數(shù)據(jù)
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相關(guān)系數(shù)
#數(shù)據(jù)初始化
n=len（XS）
r1xyr1xr1yr2xzr2xr2z=0.00.00.00.00.00.0
for?i?in?range（0n）:
????#計算xy參數(shù)
????r1xy=（XS[i]-xMean）*（YS[i]-