Intersection of 3D ray and surface of revolution
11 ビュー (過去 30 日間)
古いコメントを表示
Greetings. I'm looking for an efficient routine to find the 3D intersection point of a ray and a surface of revolution. Assume we have the generic 3D vase of height H, generated by revolving a curve of pre-stored data points. Sit the vase on the x-y plane with the center of the base over the origin. Now shoot a 3D ray from the origin in an az-el direction to an altitude of H, long enough to travel through the vase surface. The shape of the vase is such that the ray will only hit one point on the surface, probably between the predefined data points. Can the community weigh in on the most efficient way to find the intersection point? I have several brute-force solutions that look downright ugly but seem to produce reasonable results. Thanks.
4 件のコメント
採用された回答
Matt J
2013 年 4 月 25 日
編集済み: Matt J
2013 年 4 月 25 日
I'll assume, without loss of generality, that the CCW angle is zero which means we don't have to work with the 3D surface, only its pre-revolution 1D profile. Since the geometry is circularly symmetric, you can just rotate the solution by the actual, nonzero CCW angle at the end. Let the revolution profile be defined by (r,z) coordinates where vector R are the distance from the origin of your samples and vector Z are their z-coordinates. Let d be the length of your line projected onto the XY plane. Then, the solution Rmin,Zmin could be obtained by,
pp=interp1(Z,R,'linear','pp');
[~,idx]=min(abs(R-Z*d/H));
Zmin = Z(idx);
Zmin = fminsearch(@(z) abs( ppval(pp,z) - z*d/H ), Zmin);
Rmin=ppval(pp,Zmin);
2 件のコメント
Matt J
2013 年 4 月 25 日
編集済み: Matt J
2013 年 4 月 25 日
That's the drawback of restricting the scope of your question to simplistic examples, unfortunately. I suggest you Accept-click my answer, since it does address your posted question. Then pose your true question in a separate thread.
I need to interpolate between the pre-defined data points with whichever method gives accurate enough results.
When you re-post, you should expand on this. How will you decide which method "gives accurate enough results"?
その他の回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Surface and Mesh Plots についてさらに検索
製品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!