三维立体图matlab代码.doc

clc

clear

close all

D=load('ZuoBiao.dat');

M=load('JinShu.dat');

x=D(:,2);

y=D(:,3);

z=D(:,4);

c=D(:,5);

xi=linspace(min(x),max(x),100);

yi=linspace(min(y),max(y),100);

[xi,yi]=meshgrid(xi,yi);

zi=griddata(x,y,z,xi,yi);

ci=griddata(x,y,c,xi,yi);

marker={'*','o','s','^','p'};

color={'k','r','y','c','b'};

mat={'As','Cd','Cr','Cu','Hg','Ni','Pb','Zn'};

str={'等高线','生活区','工业区','山林区','交通区','绿地区'};

for j=1:8

% 等高线图形

figure

contourf(xi,yi,zi,0:10:500);

% set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2)

% clabel(C,h,[0:10:50,50:50:300,300:100:500]);

title(['金属 ',mat{j}, ' 二维等高线分布图'])

xlabel('X')

ylabel('Y')

colormap summer

colorbar

grid on

hold on

for i=1:5

loc=c==i;

plot(x(loc),y(loc),marker{i},'markerfacecolor',color{i},'MarkerEdgeColor',color{i});

end

legend(str,'location','best')

for k=1:length(x)

text(x(k)-200,y(k)+200,num2str(M(k,j+1)),'fontsize',8);

end

end

% 三维体图

figure

h=surf(xi,yi,zi);

set(h,'cdata',ci);

colormap hsv

title('三维图立体（颜色条表示分类）')

xlabel('X')

ylabel('Y')

colorbar

hidden off

hold on

for i=1:5

loc=c==i;

plot3(x(loc),y(loc),z(loc),marker{i},'markerfacecolor',color{i});

end

str{1}='三维图';

legend(str,'location','best')