resqpy.olio.vector_utilities.points_in_triangles¶
- resqpy.olio.vector_utilities.points_in_triangles(p, t, da, projection='xy', edged=False)[source]¶
Returns 2D numpy bool array indicating which of points da are within which triangles.
- Parameters
p (numpy float array of shape (N, 2 or 3)) – points in use as vertices of triangles
t (numpy int array of shape (M, 3)) – indices into first axis of p defining the triangles
da (numpy float array of shape (D, 2 or 3)) – points to test for
projection (string, default 'xy') – one of ‘xy’, ‘xz’ or ‘yz’ being the direction of projection, ie. which elements of the second axis of p and da to use; must be ‘xy’ if p and da have shape (N, 2)
edged (bool, default False) – if True, points lying exactly on the edge of a triangle are included as being in the triangle, otherwise they are excluded
- Returns
numpy bool array of shape (M, D) indicating which points are within which triangles
Note
the triangles do not need to be in a consistent clockwise or anti-clockwise order