以延河流域2000-2015年NDVI数据为例
1 数据准备
(此步骤为可选)打开Arcgis，导入2000-2015年延河流域的NDVI数据，并将10000作为每年的NDVI数字
(【Spatial Analyst】-----【地图代数】-----【网格计算机】(批量处理) ) ) ) )。
代数式是【目标年的NDVI数据* 10000】
选择输出路径和输出文件名
提取NDVI的年平均值
【System Toolboxes】--【Spatial Analyst Tools.tbx】-- 【Multivariate】((多变量分析) ) (bandcollectionstatistics ) (带宽集中统计)
如图所示选择2000-2015年NDVI数据，设置文件保存路径和名称，确定。
通过设定的路径，获取txt文件“ndvi平均值”。
2 数据处理
2.1 Excel
选择并插入目标区域----图表(散点图) )。
用鼠标右键单击任意散布点，选择【添加趋势线】
选择【线性】、【显示式】、【显示R2】
结果如下
2.2 SPSS
利用SPSS软件计算NDVI与年之间的一元线性回归方程及其显著性
点击数据，粘贴时间、NDVI数据
选择变量并尝试更改变量名称和变量类型(数字)。
【分析】-【回归】-【线性】
选择因变量【NDVI_MEAN】和自变量【Year】
计算结果
在模型的统计中，r侧为0.610，方差分析(Anova )中f值为21.940，如果观察F分布阈值表，发现=0.05，f(f(1，n-2 )，即f1=1，F2=n-2 得到显著的线性方程y=21.157 x37204 (R2=0. 610，F=21.940 )。
其中21.157是系数，-37204是常数。
2.3 Matlab
2.3.1源代码
clcclearalldata=xls read (c : (users ) administrator )桌面) mean.xlsx ) %数据存储路径x=data ) 3360，1 ) %读取事务处理1 ) p=polyfit(x ) x，y，n )，将%多项式曲线拟合到x参数，y是变量，n是系数yi=polyval(a ) a，x ) %得到相应位置' o ' %散布点hold on yi ) )并继续绘制一维线性拟合函数R2=(corrcoef ) y，yi ).^2) %并绘制R2 )2*2)矩阵R2=R2 ) 1， 要计算1 )的' bold ' %水平轴标签ylabel('ndvi )、' FontWeight '、' bold ' %垂直轴标签text(min(x )、max(y )-100、strcat ) ) forI=1:length(x ) q=q ) y ) I )-yi ) I ) ^2； u=u(yi(I ) -均值(y ) y ) ) 2； ENDF=u/(q/) ) 16-2 )； %计算显著性检查FP_value=1-fcdf(f，1，14 ) %p值text(min(x ) max (y )，strcat('p='，string ) p_value ) ) )