@画长空_yin 噢，我突然明白了，是说这18副上个影像大小不一样，谢谢
基于Matlab的栅格数据一元线性回归及显著性检验（slope趋势分析）%by yinlichang3064@163.com在进行长时间序列的栅格数据分析时，如NDVI,fvc，LAI，NPP,GPP，需要知道每个格点的长期趋势。如果再arcgi...
发简信
IP属地：台湾
-
-
@画长空_yin 噢，我刚才看了一下ndvi是722*675，但是它成功了，这种情况正常么
基于Matlab的栅格数据一元线性回归及显著性检验（slope趋势分析）
%by yinlichang3064@163.com在进行长时间序列的栅格数据分析时，如NDVI,fvc，LAI，NPP,GPP，需要知道每个格点的长期趋势。如果再arcgi...
-
中科院在读博士，分享自己在生态遥感，水文方面的原创代码，已经转V X gong zhong号(地学分析与算法)了，请不要给我简书打赏！不要打赏！不要打赏！重要的事情说三遍！
-
你好，我用这段代码对ndvi实行趋势分析成功了，但是我用它对灯光进行分析时，reshape一直显示“请勿更改元素数目”，不知道为什么
[a,R]=geotiffread('F:\ZMJ2020\LY\2001-2018夜光遥感\2001至2018灯光数据全\灯光副本\2001.tif');%先导入某个图像的投影信息，为后续图像输出做准确
info=geotiffinfo('F:\ZMJ2020\LY\2001-2018夜光遥感\2001至2018灯光数据全\灯光副本\2001.tif');
[m,n]=size(a);
years=18; %表示有多少年份需要做回归
data=zeros(m*n,years);
k=1;
for year=2001:2018 %起始年份
file=['F:\ZMJ2020\LY\2001-2018夜光遥感\2001至2018灯光数据全\灯光副本\',int2str(year),'.tif'];%注意自己的名字形式，这里使用的名字是年prec2000.tif，根据这个可修改
bz=importdata(file);
bz=reshape(bz,m*n,1);
data(:,k)=bz;
k=k+1;
end基于Matlab的栅格数据一元线性回归及显著性检验（slope趋势分析）
%by yinlichang3064@163.com在进行长时间序列的栅格数据分析时，如NDVI,fvc，LAI，NPP,GPP，需要知道每个格点的长期趋势。如果再arcgi...