在量化人类活动对生态参量的影响时,常常引入残差分析法来进行量化,本文在matlab平台下对NDVI进行残差分析,选取的自变量是降水和温度。代码如下所示,通过下面代码,能够获取残差的空间分布及残差趋势以及统计得到的研究区每年的残差值。
%author yinlichang3064@163.com
[aa,R]=geotiffread('D:\年NDVI\2000anveage_ndvi.tif');%先投影信息
info=geotiffinfo('D:\年NDVI\2000anveage_ndvi.tif');
[m,n]=size(aa);
begin_year=2000;%用户修改
end_year=2017;%用户修改
long=end_year-begin_year+1;
temsum=zeros(m*n,long);
presum=zeros(m*n,long);
ndvisum=zeros(m*n,long);
k=1;
for year=begin_year:end_year
temp=importdata(['D:\年NDVI\',int2str(year),'TEMP.tif']) ; %根据自己名称修改,本数据名称为2000TEMP.tif
pre=importdata(['D:\年NDVI\',int2str(year),'PRE.tif']) ; %根据自己名称修改,本数据名称为2000PRE.tif
ndvi=importdata(['D:\年NDVI\',int2str(ye