基于matlab-arcgis求栅格数据的相关性和显著性分析

clc;
clear;
%栅格相关性分析
[a,R]=geotiffread(‘C:\Users\smv16\Desktop\栅格数据\ndvi\2001_modis1.tif’);
info=geotiffinfo(‘C:\Users\smv16\Desktop\栅格数据\ndvi\2001_modis1.tif’);
[m,n]=size(a);
vpdsum=zeros(mn,15);
ndvisum=zeros(mn,15);
k=1;
for year=2000:2014
fname1=strcat(‘C:\Users\smv16\Desktop\栅格数据\Wet-VPD’,int2str(year),‘Q1.tif’);
data1=importdata(fname1);
data1=reshape(data1,mn,1);
vpdsum(:,k)=data1;
fname2=strcat(‘C:\Users\smv16\Desktop\栅格数据\ndvi’,int2str(year),’_modis1.tif’);
data2=importdata(fname2);
data2=reshape(data2,mn,1);
ndvisum(:,k)=data2;
k=k+1;
end

%相关性和显著性
r=zeros(m,n)+nan; % 相关性
p=zeros(m,n)+nan; % 显著性
for i=1:m*n
ndvidata=ndvisum(i,:)’;
vpddata=vpdsum(i,:)’;
[r1,p1]=corrcoef(vpddata,ndvidata);
r(i)=r1(2);
p(i)=p1(2);
end
outname1=‘C:\Users\smv16\Desktop\栅格数据\相关性\相关性.tif’;
outname2=‘C:\Users\smv16\Desktop\栅格数据\相关性\显著性.tif’;
outname3=‘C:\Users\smv16\Desktop\栅格数据\相关性\通过0.05显著性检验的相关性.tif’;
r0=r;
r0(p>0.05)=NaN;
%%输出图像
geotiffwrite(outname1,r,R,‘GeoKeyDirectoryTag’,info.GeoTIFFTags.GeoKeyDirectoryTag);
geotiffwrite(outname2,p,R,‘GeoKeyDirectoryTag’,info.GeoTIFFTags.GeoKeyDirectoryTag);
geotiffwrite(outname3,r0,R,‘GeoKeyDirectoryTag’,info.GeoTIFFTags.GeoKeyDirectoryTag);