MATLAB计算DEM与NDVI之间的关系

    最近写论文的时候，需要做一个DEM数据与NDVI之间的空间关系，在网上也没有找到相关的资料，于是就自己写了一个程序，用来研究他们之间的关系。

    这个软件中分析的模型使用的是一元线性回归模型中的系数作为相关系数，每个栅格之间都有一一对应的关系，但是只是两个栅格之间一一对应，没有其他的值，所以没有办法用p值来进行检验，只能是作为一个定性的数据来使用，如果是要定量的话，可以私聊我，我把源码发给你，下面就是大家心心念念的程序了：

clc

clear

% 设置图像处理后导出坐标系和图像行列数值

[a,R]=geotiffread('G:\毕业论文材料\裁剪数据\年度数据\年NDVI值\00-20年平均值.tif');

info=geotiffinfo('G:\毕业论文材料\裁剪数据\年度数据\年NDVI值\00-20年平均值.tif');

[m,n]=size(a);

% 导入需要处理的NDVI数据

data1=importdata('G:\毕业论文材料\裁剪数据\年度数据\年NDVI值\00-20年平均值.tif');

data1=reshape(data1,m*n,1);

NDVIsum=data1;

% 导入需要处理的DEM数据

data2=importdata('G:\毕业论文材料\裁剪数据\DEM裁剪数据\DEM裁剪数据.tif');

data2=reshape(data2,m*n,1);

DEMsum=data2;

% 相关性

NDVI_wc_xgx=zeros(m,n);

NDVIsum=double(NDVIsum);

DEMsum=double(DEMsum);

DEMsum=[ones(size(DEMsum)) DEMsum];

for i=1:length(DEMsum)

   b=regress(NDVIsum(i),DEMsum(i));

   if b<=0

   b=NaN;

   end

   NDVI_wc_xgx(i)=b;

end

%输出语句

name1='G:\DEM与NDVI相关关系.tif';

geotiffwrite(name1,NDVI_wc_xgx,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag)

    结合两张影像可以看出，在总体上可以定性说明NDVI与DEM数据的相关关系，缺点是不能使用F检验或者T检验等方法进行验证。如果想要使用检验的话，计算时间太长 ，如果有需要我会挂在评论区供大家使用。如果对大家有用的话记得一键三联支持一下。