MATLAB sen趋势分析+MK显著性检验植被覆盖度 FVC NDVI

原创

已于 2025-01-07 13:13:35 修改 · 2k 阅读

62 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #开发语言

于 2025-01-07 12:40:25 首次发布

以下代码可直接复制使用！！！

示例数据：2000-2020年年际FVC数据，一共21张tif图

NDVI等数据也同样可以使用该方法，数据来源：立方数据学社，国家地球系统科学数据中心等都可以找到！

1、sen趋势分析

% @author geo_data_analysis@163.com
% 基于Sen的趋势值
clc
clear
[a,R]=geotiffread('E:\data\fvc\NDVI2000_VFC.tif');  
info=geotiffinfo('E:\data\fvc\NDVI2000_VFC.tif');
[m,n]=size(a);
datasum=zeros(m*n,21)+NaN; 
k=1;
for year=2000:2020
    filename=['E:\data\fvc\','NDVI',int2str(year),'_VFC.tif'];
    data=importdata(filename);
    data=reshape(data,m*n,1);
    datasum(:,k)=data;
    k=k+1;
end
trend=zeros(m,n)+NaN;
for i=1:size(datasum,1)
    data=datasum(i,:);
    if min(data)>0 
        valuesum=[];
        for k1=2:21
            for k2=1:(k1-1)
                cz=data(k1)-data(k2);
                jl=k1-k2;
                value=cz./jl;
                valuesum=[valuesum;value];
            end
        en