四维图 matlab.zip

%qq为pm2.5矩阵，第一列为pm2.5值，第二列为纬度，第三列为经度，第四列为高度 F1=qq(:,1)'.*10000; x=qq(:,2)'; y=qq(:,3)'; z=qq(:,4)'; n=length(F1); mx1=(min(x));mx2=(max(x)); my1=(min(y));my2=(max(y)); for i=1:n F(i)=round((F1(i)-min(F1))/(max(F1)-min(F1))*100)+1; end f=max(F); H=colormap(jet(f+2)); fc=255/10*F/255; for i=1:n plot3(x(i),y(i),z(i),'linestyle','.','color',H(F(i),:)) hold on end colorbar xlabel('x'),ylabel('y'),zlabel('z') axis([mx1 mx2 my1 my2]) %box on