Set vector direction to point away from a closed surface (or a point coordinate)?
8 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Dear all,
I am plotting a closed potato-shaped surface using "patch", then adding some tangent planes using
plane=surf(x,y,z)
But when I am plotting the normals of those planes, using
normals=plane.FaceNormals
some of them are pointing "inwards" towards the potato, while others point outward away from it. How can I make them all point away from the surface (or the center) of the potato?
Thank you kindly!
0 Commenti
Risposte (1)
David Goodmanson
il 19 Dic 2017
Hi LMS,
Say you set up a point B in the middle of the volume. Assuming the surface is not excessively curved, from B you should be able to reach every point on the skin with a straight line without leaving the potato. Suppose dface is the displacement vector from B to the middle of a face. The the dot product of dface with an outward face normal should be positive. The code below checks for that and and changes the sign of the face normal accordingly.
I'm assuming you have m x n matrices of the x,y and z coordinates. If you fed x and y vectors to surf, then you will have to use meshgrid to get the appropriate x and y matrices.
nvecs is a (m-1) x (n-1) x 3 matrix of the FaceNormal vectors, which you got in the way you described.
v = cat(3,x,y,z); % construct vertices of faces
B = mean(mean(v)); % center, approximately
v = v - B; % translate to put center point at the origin
% find centers of faces
vcf = (1/4)*(v(1:end-1,1:end-1,:)+v(1:end-1,2:end,:) ...
+ v( 2:end,1:end-1,:)+v( 2:end,2:end,:));
sdot = sign(dot(vcf,nvecs,3)) % want these all to be positive
nvecs_new = nvecs.*repmat(sdot,1,1,3);
sdot = sign(dot(vcf,nvecs_new,3)) % check, all positive
I think the third line will have to be done with repmat if you don't have a fairly new version of Matlab.
You can eyeball a value for B, but it would have to be entered as B(1,1,:) = [bx by bz]
1 Commento
Esteban Garcia
il 29 Set 2022
Is there any other way ?, because my surface can be excessively curved...
Thank you
Vedere anche
Categorie
Scopri di più su Surface and Mesh Plots in Help Center e File Exchange
Prodotti
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!