Sen趋势分析+MK显著性分析---matlab实现

%Sen趋势

clear
clc
%datasum2(:,a)=mean(data1(find(isnan(data1) == 0)));
[a,R]=geotiffread('C:\Users\。\Desktop\data\2000.tif');%先导入投影信息
info=geotiffinfo('C:\Users\。\Desktop\data\2000.tif');
[m,n]=size(a);
cd=2020-2000+1;%时间跨度，根据需要自行修改
datasum=zeros(m*n,cd)+NaN; 
k=1;
for year=2000:2020 %起始年份
%     for month=1:12
%          month=num2str(month,'%02d');
    filename=['C:\Users\。\Desktop\data\',int2str(year),'.tif'];
         data=importdata(filename);
         data(data<-1000)=nan;
         data=reshape(data,m*n,1);
         data=double(data);
         datasum(:,k)=data;
         %datasum(:,k)=data(find(isnan(data) == 0));
         k=k+1;
%     end
end
result=zeros(m,n)+NaN;
for i=1:size(datasum,1)
    data=datasum(i,:);
%      if min(data)>-10000 %判断是否是有效值,我这里的有效值必须大于0
        valuesum=[];
        for k1=2:cd
            for k2=1:(k1-1)
                cz=data(k1)-data(k2);
                jl=k1-k2;