ff=zeros(40,40,40);
[x,y,z]=distancep(10,500,-150,-10,40,‘1.csv’);
ff(:,:,10)=z;
[x,y,z]=distancep(10,500,-150,-10,40,‘2.csv’);
ff(:,:,20)=z;
x=linspace(10,500,40);
y=linspace(-150,-10,40);
z=linspace(1,40,40);
[xx,yy,zz]=meshgrid(x,y,z);
hold on;
%h=slice(xx,yy,zz,ff,[],[],[10,20],30);%有30个等级
%h=contourslice(xx,yy,zz,ff,[],[],[10,20],30);%有30个等级,好用
%alpha(0.8) % 透明控制
%h=contourslice(xx,yy,zz,ff,[],[],[10,20],30);%有30个等级,好用
surf(xx(:,:,1),yy(:,:,1),zz(:,:,10),ff(:,:,10));
surf(xx(:,:,1),yy(:,:,1),zz(:,:,20),ff(:,:,20));
shading interp %插值