Deformed surf using griddata ?

3 ビュー (過去 30 日間)
Faez Alkadi
Faez Alkadi 2017 年 8 月 24 日
コメント済み: Faez Alkadi 2017 年 8 月 25 日
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])
  1 件のコメント
Walter Roberson
Walter Roberson 2017 年 8 月 24 日
I notice that only 554 of the 604 points are unique.

サインインしてコメントする。

採用された回答

John D'Errico
John D'Errico 2017 年 8 月 24 日
編集済み: John D'Errico 2017 年 8 月 24 日
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)
  8 件のコメント
Stephen23
Stephen23 2017 年 8 月 25 日
Another option would be to fill "missing data" areas with NaNs, which will not be plotted.
Faez Alkadi
Faez Alkadi 2017 年 8 月 25 日
編集済み: Faez Alkadi 2017 年 8 月 25 日
John and Stephen,
Thank you so much for the explanation,
The data I got is from an .stl surface file. I'm wondering if there is any easier way to meshgrid and surf the a face of an .stl file.

サインインしてコメントする。

その他の回答 (1 件)

Walter Roberson
Walter Roberson 2017 年 8 月 24 日
You can get closer to what you want if you add the 'nearest' option to the griddata() call.
You will get an odd artifact, which I figure is due to peculiarities about which is the "closest" point in 3-space.
There is no reason for the surface to follow the edge using griddata(). Gridding takes place in 3-space, so when you get sufficiently far from the edges, the closest points might be from a different z than you were picturing.
If you have a surface that you know should be a ribbon, then there are ways to decompose that.
  2 件のコメント
Faez Alkadi
Faez Alkadi 2017 年 8 月 24 日
編集済み: Faez Alkadi 2017 年 8 月 24 日
Is there any other way to surf it and make it follow the shape edges other than using griddata?
Thank you so much
Faez Alkadi
Faez Alkadi 2017 年 8 月 25 日
Walter,
The data I got is from an .stl surface file. I'm wondering if there is any easier way to meshgrid and surf the a face of an .stl file.
Thank you

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeSurfaces and Volumes についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by