You could try using the return matrices Xi,Yi,Zi to make an interpolating function. Then find z values corresponding to any x,y points you want.
An example, using scatteredInterpolant,
[t,r] = meshgrid(linspace(0,2*pi,361), ...
[x,y] = pol2cart(t,r);
P = peaks(x,y);
angR = [20 350]*pi/180;
radR = [.8 4];
[xi,yi,zi] = ...
zFun = ...
Excellent script; thank you very much for it.
I would like to interpolate across the center of the plot(constituted by the first row or the input vector), as well as between the first and last column of the input vector.
Could you advise me on how to proceed with this?