%根据读入的散点生成三角网和三维模型
clear;
%[x,y,z]=textread('mm.xyz','%f%f%f');
[x,y,z]=textread('gridpoint.txt','%f%f%f');
TRI = delaunay(x,y,{'Qt','Qbb','Qc','Qz'});
%triplot(TRI,x,y)%画出三角形
figure(1);
trisurf(TRI,x,y,z)
%trimesh(TRI,x,y,z)
%显示指定范围的等高线：
figure(2);
%tix = 679000:500:683000;
%tiy = 234000:500: 238000;
%[XI,YI] = meshgrid(tix,tiy);%产生方格交叉点的坐标
%ZI = griddata(x,y,z,XI,YI);%根据散点x,y,z拟合出曲面，内插产生各方格点的Z坐标
%利用方格点数据内插出等高线：
tix = 681500:50:683000;
tiy = 234000:50:235500;
[XI,YI] = meshgrid(tix,tiy);%产生方格交叉点的坐标
ZI = griddata(x,y,z,XI,YI,'cubic',{'Qt','Qbb','Qc','Qz'});%根据散点x,y,z拟合出曲面，内插产生各方格点的Z坐标
contour(XI,YI, ZI );
%Plot the gridded data along with the nonuniform data points used to generate it:
%mesh(XI,YI,ZI);
%hold
%plot3(XI,YI,ZI,'o')
%plot3(x,y,z,'o')
%hold off
%hold on;%两图叠加，保留上一个图
%subplot(2,1,2);
%plot(x,y,'or');
%hold off
%Plot the gridded data along with the nonuniform data points used to generate it:
%trisurf(TRI,x,y,z)
%trimesh(TRI,x,y,z)
%mesh(x,y,z);
% ZI=zeros(length(XI) );plot3(XI,YI,ZI,'o')
% plot3(XI,YI,ZI,'o')
% plot(XI,YI,'o')
%[C,h] = contour(X,Y,Z);
%set(h,'ShowText','on','TextStep',get(h,'LevelStep')*2)
%colormap cool
