ArcMap对多年NDVI进行线性趋势和F显著性检验.pdf

ArcMap对多年NDVI进行线性趋势和F显著性检验.pdf

《ArcMap对多年NDVI进行线性趋势和F显著性检验.pdf》由会员分享，可在线阅读，更多相关《ArcMap对多年NDVI进行线性趋势和F显著性检验.pdf（13页珍藏版）》请在冰点文库上搜索。

ArcMap对多年NDVI进线性趋势和F显著性检验1法1.1趋势分析采元线性回归分析和最乘法，逐像元拟合近n年间遥感影像每个栅格NDVI的斜率Slope，得到多年NDVI的变化趋势，综合分析流域多年植被覆盖变化额向和速率。

式中，slope为变化趋势；yi为第xi年的值；当slope0,呈增加趋势；当slope【SpatialAnalystTools.tbx】-【Multivariate】-【BandCollectionStatistics】（波段集统计）按照图中式选择2000-2015年NDVI数据，设置件保存路径和名称，OK。

在设置的路径下，得到txt件“ndvi平均值”3、得到元线性回归程在Excel中打开“ndvi平均值”件，选定标区域，【插】-【图表】（散点图）标右击个散点，添加【趋势线】选择【线性】、【显公式】、【显R平】输出结果4、G显著性检验利SPSS软件，计算NDVI与年份之间的元线性回归程的显著性【分析】-【回归】-【线性】，选择因变量和变量，计算结果。

模型汇总中R2为0.610，差分析（Anova）中，F值为21.940，查F分布临界值表，当=0.05时，FF（1,n-2），即f1=1，f2=n-2=16-2=14，F0.05临界值为4.60，F=21.940F0.05=4.60，可知拟合的线性程是显著的。

系数中，常量为-37204.226，变量系数为21.157，Sig.为0.000b【SpatialAnalystTools.tbx】-【MapAlgebra】-【RasterCalculator】按照下图的步骤分别计算b1、b2、b值。

2.2.2常数a的计算a值计算按照公式计算，分别计算x和y平均值，即2000-2015年年份和NDVI的平均值，年份的平均值为2007.5。

打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【Local】-【CellStatistics】计算NDVI平均值计算a值打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【MapAlgebra】-【RasterCalculator】2.2.3线性回归程的显著性检验计算F值计算2000-2015年拟合NDVI值打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【MapAlgebra】-【Raster】（批处理）命名为Fit_ndvi2000、Fit_ndvi2001、Fit_ndvi20计算Q值打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【MapAlgebra】-【Raster】计算U值打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【MapAlgebra】-【Raster】计算F值，F值计算公式如下：

打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【MapAlgebra】-【Raster】2.2.4成果输出线性变化速率（系数b值）和显著性检验结果（F值）。

线性变化速率（b值）：

是以0为界值，将线性变化速率划分为正值和负值，正值表线性增加趋势、负值表线性减少趋势将b重分类：

打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【Reclass】-【Reclassify】显著性检验结果（F值）：

根据显著性临界值（=0.05时，n=16是，F表临界值为4.60），将F值划分为显著减少趋势、不显著减少趋势、不显著增加趋势、显著增加趋势4类。

将F重分类：

打开ArcGIS，【SystemToolboxes】-【SpatialAnalystTools.tbx】-【Reclass】-【Reclassify】趋势及显著性叠加检验：

利RasterCalculator按照图计算F检验结果2.2.5统计积和例导出属性表，在Excel中打开计算各个类型的积和例