load ('abu_jianbeijing.mat');
Z=abu_jianbeijing;
M=Z;
imsequence=Z(1:256,1:320,1:100);%取到了数据
imsequence1=imsequence(:,:,1:100);
imsequence2=imsequence1(:,:,1:100);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%多项式拟合
for j=1:1:256
for k=1:1:320
T=imsequence2(j,k,1:100);
for i=1:1:100
TT(i)=T(i);%%%%%减去初始温度!
end
t=1:1:100;
% p=polyfit(t,TT,5);
TTT=reshape(TT,1,100); %返回一个1*580的矩阵TT,TTT中数据从TT中得到
p=polyfit(log(t),log(TTT),2); %%%%%%重要!由于原来的TT维数为1*1*580,而t为1*580,因此需要将TT维数变为与t一样,才可以进行拟合!
for w=1:1:3
yy(j,k,w)=p(w); %%%%此处注意,不要再使用y,因为前面用过了!!!!
end
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%以上程序提取了每个像素点对应的多项式系数
for i=1:1:3
for j=1:1:256
for k=1:1:320
zz(j,k)=yy(j,k,i); %%%%此处注意,不要再使用z,y,因为前面用过了!!!!
end
end
m=1:1:100;
n=1:1:100;
figure(i)
imagesc(real(zz))
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%上面的程序已经完全成功!
评论0