【Matlab】多年NDVI数据的sen趋势分析

具体内容及详细教程请关注微信公众号：夫也的笔记

最近老师在课上安排了一个汇报作业，要求对某一地区的多年NDVI数据进行趋势分析，对于小白的我来说，真是！唉！但是功夫不负有心人，在我充分利用百度搜索、谷歌搜索和Bing搜索后，终于让我找到了相关教程和知识。只是需要用Matlab软件进行一系列操作，废话不多说上代码。
趋势分析有很多种，今天介绍的是Theil-Sen median趋势分析，一般还要结合Mann-Kendall建议，这个MK检验后续文章会写。
Theil-Sen median趋势分析是一种稳健的非参数计算方法，计算式如下

公式中SET指计算n(n-1)/2 个数据组合的斜率的中位数，ETi和ETj代表i和j年的ET值。如果SET＞0，则ET呈上升趋势，否则，ET呈下降趋势。
Matlab具体代码如下：

% @author 1154318421@qq.com
[a,R]=geotiffread('E:\GIS\NDVI\2000.tif');%先导入投影信息
info=geotiffinfo('E:\GIS\NDVI\2000.tif');
[m,n]=size(a);
cd=2020-2000+1;%时间跨度，根据需要自行修改
datasum=zeros(m*n,cd)+NaN;
k=1;
for year=2000:2020 %起始年份filename=['E:\GIS\ENVI\xinjiang',int2str(year),'2000.tif'];data=importdata(filename);data=reshape(data,m*n,1);datasum(:,k)=data;k=k+1;
end
result=zeros(m,n)+NaN;
for i=1:size(datasum,1)data=datasum(i,:);if min(data)>0 %判断是否是有效值,我这里的有效值必须大于0valuesum=[];for k1=2:cdfor k2=1:(k1-1)cz=data(k1)-data(k2);jl=k1-k2;value=cz./jl;valuesum=[valuesum;value];endendvalue=median(valuesum);result(i)=value;end
end
filename=['E:\GIS\基于SEN的NDVI变化趋势.tif'];
geotiffwrite(filename,result,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag)在这里插入代码片

【Matlab】多年NDVI数据的sen趋势分析相关推荐

matlab中NDVI趋势分析,基于matlab 的长时间栅格数据的sen趋势分析
基于matlab 的长时间栅格数据的sen趋势分析 sen趋势分析是进行趋势分析的方法之一,常配合MK检验来使用,在学术论文中也非常常见,本文主要介绍如果通过matlab来实现对长时间序列的栅格数据进 ...
长时间序列遥感数据植被物候提取/遥感数据产品分析暨MODIS NDVILAI多年产品数据批处理分析/Python长时间序列遥感数据处理及在全球变化、物候提取、植被变绿与固碳分析、生物量估算与趋势分析
基于MATLAB长时间序列遥感数据植被物候提取与分析 1.本课程基于matlab语言 2.提供所有代码 3.以实践案例为课程内容主线,原理与操作相结合 4.根据讲解内容,布置作业,巩固所学内容及拓展在 ...
Matlab处理气象数据（九）逐点变化的显著性分析
利用上面M-K趋势分析的结果,只让置信度大于95%的数据在地图上显示. 这里,可以先自定义一个colorbar. mycolor = colormap; %mycolor名称任意 save mycol ...
将MATLAB的变量数据导入到C/C++程序中的方法
2016-5-21日:我发现自己真是傻,既然对象是图像,那何必把像素值导出来呢?直接将图像输出,然后VC再读入图像不就可以了? 最开始,我想在VC中直接调用MATLAB-Workspace中的变量,即 ...
利用matlab对rosbag数据,通过matlab提取rosbag数据
通过matlab提取rosbag数据代码如下: 最近在做AGV导航实验,导航小车是基于ROS实现的,判断小车导航是否运行稳定,可以直接通过rviz直观看出来.但是如何通过数据定性地分析数据,这就需要 ...
R:matlab交互，数据调用
目录: 前言首先安装R.matlab包第二步,加载R.matlab包前言最近在学习使用R语言,但是数据源来自于之前用matlab保存的数据文件(.mat格式),因此有了使用R读取mat文件的需 ...
分酒问题matlab代码,matlab葡萄酒分类数据归一化问题
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 matlab葡萄酒分类数据归一化问题% 选定训练集和测试集 % 将第一类的1-30,第二类的60-95,第三类的131-153做为训练集 train_wi ...
matlab读入txt数据_教程合集 | MATLAB文件读写（以nc与txt为例）
周一到!我们又和大家见面啦!最近几周的主题是,有了数据,那我们怎么读写数据文件呢?气象中的特殊资料又有哪些处理方法呢?以txt和nc文件做引子,欢迎各位拍砖~ No.1 基础通用知识读写文件时有一些 ...
最新Modis影像数据下载完整流程---以MOD13Q1 NDVI数据下载为例（2020年12月2日更新）
本文详细讲解Modis数据的下载流程,以MOD13Q1影像为例,该数据为16天合成的空间分辨率为250m的NDVI产品,原始数据的投影为正弦曲线投影. Modis数据MRT和ArcGIS处理方法集锦: ...

【Matlab】多年NDVI数据的sen趋势分析

【Matlab】多年NDVI数据的sen趋势分析相关推荐

最新文章

热门文章

【Matlab】 多年NDVI数据的sen趋势分析

【Matlab】 多年NDVI数据的sen趋势分析相关推荐

最新文章

热门文章

【Matlab】多年NDVI数据的sen趋势分析

【Matlab】多年NDVI数据的sen趋势分析相关推荐