基于Matlab的通径分析法代码

基本原理

解方程即可求出通径系数ai， ai 表示自变量 xi 对因变量 y 的直接通径系数，为 xi 对 y 的直接作用效应； r xi,xj* ai 表示自变量 xi 通过 xj 对因变量 y 的间接通径系数，为xi 通过 xj 对因变量 y 的间接作用效应。
以一个例子说明：气温 （降水） 对植被绿度的直接影响是在不考虑降水 （气温） 变化的影响时，由于气温 （降水） 的变化对植被绿度产生的影响；此外，气温 （降水） 的变化会引起降水（气温） 变化，进而影响植被绿度，即气温 （降水） 对植被绿度的间接影响；气温 （降水）对植被绿度的综合影响包括直接影响和间接影响。即植被的变化不能忽视气候因子间的相互作用。

---

部分源代码

这段MATLAB代码用于计算气温和降水对植被指数（NDVI）的直接和间接影响，以及它们的综合影响，并将结果保存为GeoTIFF文件。以下是对代码的详细注释：

% 清空命令行窗口和清除所有变量
clc
clear
% 读取NDVI数据和获取文件信息
[a,R]=geotiffread('2020_ndvi.tif');  
info=geotiffinfo('2020_ndvi.tif');
[m,n]=size(a); % 获取NDVI图像的尺寸

% 设置气温和降水数据文件的路径
Tpath='气温'; % 气温路径
Ppath='降水'; % 降水路径

% 初始化存储气温和降水数据的矩阵
Tsum=zeros(m*n,21); % 气温数据矩阵
Psum=zeros(m*n,21); % 降水数据矩阵

% 循环读取2000年到2020年的气温和降水数据，并转换为列向量
for year=2000:2020
    temp=importdata([Tpath,int2str(year),'temp.tif']);  
    Tsum(:,year-1981)=reshape(temp,m*n,1);
    prec=importdata([Ppath,int2str(year),'prec.tif']);
    Psum(:,year-1981)=reshape(prec,m*n,1);   
end

% 初始化存储气温和降水相关性的向量
xgx=zeros(m*n,1)+nan;

% 计算每个像素点气温和降水的相关性
for i=1:m*n
    if min(Tsum(i,:))~=-999   % 排除背景值
        r=corrcoef(Tsum(i,:),Psum(i,:)); % 计算相关系数
        xgx(i)=r(2); % 存储相关系数
    end
end
xgx=reshape(xgx,m,n); % 将相关系数向量重新整形为图像尺寸

% 读取气温和降水对NDVI直接影响的GeoTIFF文件
At=importdata(('C:\Users\DELL\Desktop\通径分析\结果\气温对NDVI直接影响.tif'));
Ap=importdata(('C:\Users\DELL\Desktop\通径分析\结果\降水对NDVI直接影响.tif'));

% 计算气温和降水对NDVI的间接影响
Atp=xgx.*Ap;   % 气温对NDVI间接影响
Apt=xgx.*At;   % 降水对NDVI间接影响

% 保存间接影响结果为GeoTIFF文件
out1=strcat('气温对NDVI间接影响.tif');
geotiffwrite(out1,Atp,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
out2=strcat('降水对NDVI间接影响.tif');
geotiffwrite(out2,Apt,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);

% 计算气温和降水对NDVI的综合影响
rtn=At+Atp; % 气温对NDVI综合影响
rpn=Ap+Apt; % 降水对NDVI综合影响

% 保存综合影响结果为GeoTIFF文件
out3=strcat('结果\气温对NDVI综合影响.tif');
geotiffwrite(out3,rtn,R,'GeoKeyDirectoryTag',info.Geo

欢迎关注gongzhonghao 趣品科研，获取更多代码和前沿论文资讯等相关内容

出图效果

参考文献

[1] 蔡甲冰, 许迪, 刘钰, 等. 冬小麦返青后腾发量时空尺度效应的通径分析. 农业工程学报, 2011, 27(8): 69-76. [Cai Jiabing, Xu Di, Liu Yu, et al. Path analysis on spatio-temporal scaling effect of crop evapotranspiration in growing seasons for winter wheat after reviving. Transactions of the Chinese Society of Agricultural Engineering, 2011, 27(8): 69-
76.]
[2] 徐春华, 张华, 张兰, 等. 基于通径分析的兰州北山三种典型植物光合作用影响因子. 生态学杂志, 2015, 34(5): 1289-
1294. [Xu Chunhua, Zhang Hua, Zhang Lan, et al. Factors influencing photosynthesis of three typical species in Beishan
Mountain of Lanzhou based on path analysis. Chinese Journal of Ecology, 2015, 34(5): 1289-1294. ]
[3] 谭丽萍, 刘苏峡, 莫兴国, 等. 华北人工林水热碳通量环境影响因子分析. 植物生态学报, 2015, 39(8): 773-784. [Tan
Liping, Liu Suxia, Mo Xingguo, et al. Environmental controls over energy, water and carbon fluxes in a plantation in
Northern China. Chinese Journal of Plant Ecology, 2015, 39(8): 773-784. ]
[4]谢舒笛,莫兴国,胡实等.三北防护林工程区植被绿度对温度和降水的响应[J].地理研究,2020,39(01):152-165.