matlab中偏相关,基于matlab的逐像元偏相关分析

本文介绍基于matlab的逐像元的三者间的偏相关，以产水，NPP和土壤保持为例进行说明

具体代码如下：

%将三者多年的数据放在三个不同的矩阵中

[a,R]=geotiffread('F:\校级课题项目\data\屏障带\2002water_yield.tif');%先导入投影信息

info=geotiffinfo('F:\校级课题项目\data\屏障带\2002water_yield.tif');

nppsum=zeros(size(a,1)*size(a,2),16);

for year=2000:2015

filename=strcat('F:\校级课题项目\data\屏障带\',int2str(year),'npp.tif');

data=importdata(filename);

data=reshape(data,size(a,1)*size(a,2),1);

nppsum(:,year-1999)=data;

end

scsum=zeros(size(a,1)*size(a,2),16);

for year=2000:2015

filename=strcat('F:\校级课题项目\data\屏障带\',int2str(year),'sc.tif');

data=importdata(filename);

data=reshape(data,size(a,1)*size(a,2),1);

scsum(:,year-1999)=data;

end

wcsum=zeros(size(a,1)*size(