for loop
U1=(U0.*d.A0)./((d+h).*(A0+A1));U1=(U0.*d.*A0)./((d+h).*(A0+A1));inmask = inpolygon(xg(:),yg(:), bx , by );%define which pixel exist in the polygon
inmask = inpolygon(xg(:),yg(:), x(b) , y(b) );%define which pixel exist in the polygondata=xlsread('ctd_data.xlsx');x=data(:,5); %first colum in excel file, as distance
y=data(:,3)*-1;%second colum in excel file, as depth
c=data(:,4); %c is the concentration of Temperature
figure(1)scatter(x,y,30,c,'filled')colorbaraxis tight; hold onshading interpylabel('Depth (m)'); title('Temperature distribution');figure(2)[xg,yg] = meshgrid(min(x):0.1:max(x),min(y):0.1:max(y)); %make a grid to contain x,y
f = scatteredInterpolant(x,y,c,'linear'); %assign c values to x y
vg = f(xg,yg);% assign c values to finer pixel xg yg
%find boundary
[g,xd] = findgroups(x)yd_top = splitapply(@min,y,g)yd_bot = splitapply(@max,y,g)yd = [yd_top;flip(yd_bot)];xd = [xd;flip(xd)]inmask = inpolygon(xg(:),yg(:), xd , yd );%define which pixel exist in the polygonvg(~inmask) = nan; %NAN values for the pixel not in the polygon
h = pcolor(xg,yg,vg); %ploting in 2D of vg by color
h.EdgeColor = 'none'; %to get rid off the edge line
h.FaceColor = 'interp';colorbarxlabel('Distance (km)'); ylabel('Depth (m)');
Best Answer