多年NDVI数据时间序列数据的MK检验

最新推荐文章于 2024-08-21 16:26:07 发布

夫也

最新推荐文章于 2024-08-21 16:26:07 发布

阅读量1.1w

点赞数 25

分类专栏： MATLAB 文章标签： matlab

本文链接：https://blog.csdn.net/qq_44795205/article/details/112796872

版权

本文介绍了Mann-Kendall趋势检验法，这是一种非参数统计方法，适用于非正态分布的水文气象数据趋势检测，特别是降水、干旱频次等气候变化研究。在MATLAB中进行MK检验的代码示例被提供，该方法无需数据分布检验，可处理缺失值，关注数据的相对数量级而非数值本身。

摘要由CSDN通过智能技术生成

教程照片及其他详细信息请关注微信公众号：夫也的笔记
公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享

理论知识
Mann-Kendall趋势检验法是世界气象组织推荐并已广泛应用的非参数统计方法，能有效区分某一自然过程是处于自然波动还是存在确定的变化趋势。对于非正态分布的水文气象数据，Mann-Kendall秩次相关检验具有更加突出的适用性。Mann-Kendall也经常用于气候变化影响下的降水、干旱频次趋势检测。Mann-Kendall非参数秩次检验在数据趋势检测中极为有用，其特点表现为：（1）无需对数据系列进行特定的分布检验，对于极端值也可参与趋势检验；（2）允许系列有缺失值；（3）主要分析相对数量级而不是数字本身，这使得微量值或低于检测范围的值也可以参与分析；（4）在时间序列分析中，无需指定是否是线性趋势。

在前一篇文章中讲述了用sen法进行长时间的趋势分析，但并未对结果进行显著性检验，通常Sen与MK检验是结合在一起的，因此本文主要讲述如何进行MK检验。具体代码如下：

% @author 2287335652@qq.com
clear
[a,R]=geotiffread('D:\GIS\yearmax\1982.tif'); %先导入投影信息
info=geotiffinfo('D:\GIS\yearmax\1982.tif');%先导入投影信息
[m,n]=size(a);
cd=34;       %34年，时间跨度  
datasum=zeros(m*n,cd)&#