# -*- coding: utf-8 -*-
import arcpy
from arcpy import env
from arcpy.sa import *
# 读取数据 !!!
workspace1 = "F:/xiaolunwen_data/Partial correlation"
arcpy.env.workspace = unicode(workspace1, "utf8")
Rxy = Raster("F:/xiaolunwen_data/Partial correlation/PythonPPPXG/R/Rxy.tif")
Rxz = Raster("F:/xiaolunwen_data/Partial correlation/PythonPPPXG/R/Rxz.tif")
Ryz = Raster("F:/xiaolunwen_data/Partial correlation/PythonPPPXG/R/Ryz.tif")
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
print Rxy
print Rxz
print Ryz
# 计算偏相关系数
Rxyz = (Rxy - (Rxz * Ryz)) / SquareRoot((1 - Square(Rxz)) * (1 - Square(Ryz)))
print Rxyz
# 输出偏相关系数
out = "F:/xiaolunwen_data/Partial correlation/PythonPPPXG/偏相关系数/"
out1 = unicode(out, "utf8") + "Rxyz" + ".tif"
Rxyz.save(out1)
print ('xyz good luck!')
# 偏相关系数t值检验
#T = Rxyz * SquareRoot(21 - 2 - 1) / SquareRoot(1 - Square(Rxyz))
#print T
#out5 = unicode(out, "utf8") + "T" + ".tif"
#T.save(out5)
#print 5
# 计算复相关系数
#Fxyz = SquareRoot(1 - (1 - Square(Rxy)) * (1 - Square(Rxyz)))
#print Fxyz
#out6 = unicode(out, "utf8") + "Fxyz" + ".tif"
#Fxyz.save(out6)
#print 6
# 复相关系数做F值检验
#F = Square(Fxyz) / (1 - Square(Fxyz)) * ((16 - 2 - 1) / 2)
#print F
#out7 = unicode(out, "utf8") + "F" + ".tif"
#F.save(out7)
#print 7
偏相关分析（一阶偏相关、二阶偏相关）
需积分: 0 175 浏览量
2022-09-01
09:24:51
上传
评论 10
收藏 16KB ZIP 举报
代码.zip （21个子文件） 
二阶偏相关系数输出代码 
Rxyzw.py 1KB Rxzyw.py 1KB Rxwzy.py 1KB R输出代码 Rzw_out.py 1KB Ryw_out.py 1KB Rxz_out.py 1KB Rxw_out.py 1KB Rxy_out.py 1KB Ryz_out.py 1KB 一阶偏相关系数输出代码 Rwyz.py 885B Rzyw.py 885B Rxyw.py 883B Rxzy.py 1KB Rxwz.py 883B Ryzw.py 883B Rwzy.py 885B Rywz.py 883B Rxwy.py 883B Rxzw.py 883B Rxyz.py 1KB Rzwy.py 883B
评论0