对tif数据进行最小二乘计算拟合时出现散点向上,拟合线向下,不知道是哪里出了问题
tifdiry = 'E:\2023上半年\20230215\b81.tif'; %文件夹下tif数据的名称?
tifdirx = 'E:\2023上半年\20230215\b41.tif'; %文件夹下tif数据的名称?
savdir = [ 'E:\2023上半年\20230215\处理数据\'];
[tm1,R] = geotiffread(tifdiry);
[tm2,R] = geotiffread(tifdirx);
e1=reshape(tm1,size(tm1,1)*size(tm1,2),1);
e2=reshape(tm2,size(tm2,1)*size(tm2,2),1);
x=e2;
y=e1;
p = polyfit(x,y,1);
y1 = polyval(p,x);
% mdl = fitlm(x,y);%求一元线性拟合的参数
% r2 = num2str(mdl.Rsquared.Ordinary);%即一元线性拟合的R平方
figure; % 绘图
scatter(x,y,'.');
hold on;
plot(x,y1,'Color','r');
x=e1?;y=e2?