ArcMap对多年NDVI进⾏线性趋势和F显著性检验 1⽅法
1.1 趋势分析
采⽤⼀元线性回归分析和最⼩⼆乘法，逐像元拟合近n年间遥感影像每个栅格NDVI的斜率Slope，得到多年NDVI的变化趋势，综合分析流域多年植被覆盖变化额⽅向和速率。 式中，slope为变化趋势；yi为第xi年的值；当slope > 0，呈增加趋势；当slope < 0，呈下降趋势；
1.2 F检验麦迪霉素
为了进⼀步评价植被覆盖变化状况，采⽤F检验法对NDVI变化趋势进⾏显著性分析，⽤于表⽰趋势变化置信度的⾼低。声带小节
式中，n为研究时序；U为误差平⽅和；Q为回归平⽅和；~yi为拟合回归值；⼀yi为n年的平均值；yi为第xi年的值。
2 数据处理
2.1 时间上
1、打开Arcgis，导⼊2000-2015年延河流域的NDVI数据，将各个年份的NDVI的数值乘以10000
【Spatial Analyst】-- 【地图代数】 -- 【栅格计算器】（批处理）
地图代数表达式为： ⽬标年份的NDVI数据 * 10000
选择输出路径和输出⽂件名
2、提取 NDVI 年平均值
【System Toolboxes】 ->【 Spatial Analyst Tools.tbx】 -> 【Multivariate】 ->【Band Collection Statistics】（波段集统计）
按照图中⽅式选择 2000-2015 年 NDVI 数据，设置⽂件保存路径和名称，OK。
在设置的路径下，得到 txt ⽂件“ndvi 平均值”
3、得到⼀元线性回归⽅程
在Excel中打开“ndvi 平均值”⽂件，选定⽬标区域，
硫氰酸钠【插⼊】 -- 【图表】（散点图）克里雅人
血凝素⿏标右击⼀个散点，添加【趋势线】
选择 【线性】、【显⽰公式】、【显⽰R平⽅】
输出结果
4、G 显著性检验
利⽤ SPSS 软件，计算 NDVI 与年份之间的⼀元线性回归⽅程的显著性
【分析】->【回归】->【线性】，选择因变量和⾃变量，计算结果。
模型汇总中 R 2为 0.610，⽅差分析（Anova）中，F 值为 21.940，查 F 分布临界值表，当 α=0.05 时，F ~ F(1, n-2)，即f1=1，f2=n-2=16-2=14，F0.05临界值为 4.60， F =21.940 > F0.05 = 4.60，可知拟合的线性⽅程是显著的。系数中，常量为-37204.226，变量系数为 21.157，Sig.为 0.000b < 0.05，表明系数也是显著的。拟合⽅程显著情况下，才可以进⾏预测分析。
5、H 预测分析
利⽤拟合得到的显著线性⽅程 y = 21.157x – 37204 (R2 = 0.610， F=21.940 )。其中 21.157 为系数，-37204 为常数，R2为曲线和拟合度。可以预测未来的 NDVI 变化。
2.2 空间上
论表见代理
2.2.1 系数 b 的计算
根据公式，b 值可划分为以下⼏个部分分别计算：
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议