MATLAB: Deformed surf using griddata

griddatameshgridsurf

I have a set of data as attached. but when I surf them after using griddata. the result is deformed surf AS SHOWN IN THE PICTURE with red squares. Z values on the surf were supposed to follow the edges of the shape. I think its because there is something wrong with (z_surf) matrix in my code. But I don't know how to solve it.
Hope you can help me to solve this problem.
Thank you.
 x_limit=xyz_board_ordered_dvd(:,1);
 y_limit=xyz_board_ordered_dvd(:,2);
 z_limit=xyz_board_ordered_dvd(:,3);
 dx=1;
 dy=1;
x_edge=floor(min(x_limit)):dx:ceil(max(x_limit));
y_edge=floor(min(y_limit)):dy:ceil(max(y_limit));
[x_surf,y_surf]=meshgrid(x_edge,y_edge);
z_surf=griddata(x_limit,y_limit,z_limit,x_surf,y_surf);
surf(x_surf,y_surf,z_surf)
xlabel('X')
ylabel('Y')
xlim([0 120])
ylim([0 120])

Best Answer

Easy, peasy.
Griddata will fail here. In fact, my own gridfit will produce some artifacts.
But essentially, you have a simple ribbon of the form z(y), which seems to be independent of x.
x = xyzSTL_board_ordered_dvd(:,1);
y = xyzSTL_board_ordered_dvd(:,2);
z = xyzSTL_board_ordered_dvd(:,3);
plot(y,z,'o')
There is NO x in that relationship. So all you need to do is fit the curve z(y). Use a spline. Start out by removing any replicates. Using my consolidator tool, as found on the file exchange :
[yhat,zhat] = consolidator(y,z,@mean);
plot(yhat,zhat,'-')
spl = pchip(yhat,zhat);
Use pchip here. This is fairly important.
[xgrid,ygrid] = ndgrid(linspace(min(x),max(x),10),linspace(min(y),max(y),100));
zgrid = repmat(ppval(spl,linspace(min(y),max(y),100)),10,1);
surf(xgrid,ygrid,zgrid)
Related Question