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);                                    
N = 45;                                              
[X,Y] = meshgrid(linspace(xmin,xmax,N),linspace(ymin,ymax,N));
Z = griddata(x,y,z,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)
mesh(X,Y,Z)
title('BRENT层等厚度构造(3D)')
colormap('JET');                                    
xlabel('LINE')
ylabel('CDP')
zlabel('T0')
colorbar;