学习笔记之 基于matlab的栅格数据相关分析及显著性检验

相关分析在科研过程中常常会碰到，表征了一种两个因素间的相关程度，其值的大小只能说明相关性的强弱。
比如当我们有2000-2015年的NPP数据和2000-2015年的降水数据时，我们想查看两者在空间上随时间变化的
相关性。本文以产水和NPP为例进行说明
具体代码如下：
%authour@[email protected]

%将两者多年的数据放在三个不同的矩阵中
nppsum=zeros(3587*4642,16);
for year=2000:2015
    filename=strcat('F:\课题项目\data\',int2str(year),'npp.tif');
    data=importdata(filename);
    data=reshape(data,3587*4642,1);
    nppsum(:,year-1999)=data;
end

wcsum=zeros(3587*4642,16);
for year=2000:2015
    filename=strcat('F:\课题项目\data\',int2str(year),'water_yield.tif');
    data=importdata(filename);
    data=reshape(data,3587*4642,1);
    wcsum(:,year-1999)=data;
end
%相关性和显著性
npp_wc_xgx=zeros(3587,4642);
npp_wc_p=zeros(3587,4642);
for i=1:length(nppsum)
    npp=nppsum(i,:);
    if min(npp)>