lonx=min_x:20:max_x;
laty=min_y:20:max_y;
[Xi,Yi]=meshgrid(lonx,laty);
bathymetry=interp2(XX,YY,ZZ,Xi,Yi);
figure(1)
pcolor(Xi,Yi,bathymetry),shading flat ,colorbar;
saveas(gcf,'bat_fine','eps');
Also you can striaght away try this:
clear;
format long e
T=readtable('points.txt');
A=T{:,:};
X=A(:,1);
Y=A(:,2);
Z=A(:,3);
min_x=min(X);
max_x=max(X);
min_y=min(Y);
max_y=max(Y);
n=length(X);
lonx=min_x:20:max_x;
laty=min_y:20:max_y;
[Xi,Yi]=meshgrid(lonx,laty);
bathymetry=interp2(X,Y,Z,Xi,Yi);
ZZ = griddata(X,Y,Z,Xi,Yi);
figure(1)
pcolor(lonx,laty,bathymetry),shading flat ,colorbar;
saveas(gcf,'bat_fine','eps');
Best Answer