clc;clear;close all;
data= xlsread('Brent等厚度.xlsx'); 
x = data(:,1)               
xmax = max(x);                                         
xmin = min(x);
y = data(:,2)
ymax = max(y);
ymin = min(y);
z = data(:,3);
z1 = data(:,4);
N = 45;                                              
[X,Y] = meshgrid(linspace(xmin,xmax,N),linspace(ymin,ymax,N));
Z = griddata(x,y,z,X,Y,'v4');
Z1 = griddata(x,y,z1,X,Y,'v4');
figure(1)
contourf(X,Y,Z,N, 'LineColor','none');  
title('Brent等厚度(2D)');                      
colormap('JET');                                     
xlabel('LINE')
ylabel('CDP')
colorbar;
%saveas(gcf,'data1','png')                     
figure(2)S1=size(Z1);%曲面1大小
C1=ones(S1);
mesh(X,Y,Z,C1)
hold on;
mesh(X,Y,Z1)
colormap([1,0,1;0,1,1])
title('Brent层(3D)')                            
xlabel('LINE')
ylabel('CDP')
zlabel('T0')
set (gca,'ZDir','reverse')