MixoSrSSKrSSKJrJrJr SSKJr SSKJ r SSK J r SSK J r SS jrS rSS jrSS jrSS jrSSjrg)zT intersections.py ------------------ Primarily mesh-plane intersections (slicing). N)geometrygroupingutil)transformations) triangles)tol)points_to_barycentricc X^^SnUU4SjnSnUU4Sjn [R"T[RS9m[R"T[RS9mTRS:wdTRS:wa [ S5eUc UR n O2[R"U[R S9nUR Un UbUn O$[R"URT- T5n [R"[UR5[RS9n SX[R*:'S X[R:'Xn U"U 5n XU4n[R"[X5VVs/sHunnU"XXUR5PM snn5nU(ah[R "U Vs/sHn[R""U5S PM sn5nUR$R&S :XdeUcUU4$UUU4$U$s snnfs snf) a Find a the intersections between a mesh and a plane, returning a set of line segments on that plane. Parameters --------- mesh : Trimesh object Source mesh to slice plane_normal : (3,) float Normal vector of plane to intersect with mesh plane_origin : (3,) float Point on plane to intersect with mesh return_faces : bool If True return face index each line is from local_faces : None or (m,) int Limit section to just these faces. cached_dots : (n, 3) float If an external function has stored dot products pass them here to avoid recomputing. Returns ---------- lines : (m, 2, 3) float List of 3D line segments in space. face_index : (m,) int Index of mesh.faces for each line Only returned if return_faces was True cX[R"USS9n[R"[U5[RS9S-n[ S5HnX!SS2U4SU- -- nM [R"S[ S9nSUS 'XBnS USS&SUS 'XBnS USS&SUS S /'XBnXvU4$)a Figure out which faces correspond to which intersection case from the signs of the dot product of each vertex. Does this by bitbang each row of signs into an 8 bit integer. code : signs : intersects 0 : [-1 -1 -1] : No 2 : [-1 -1 0] : No 4 : [-1 -1 1] : Yes; 2 on one side, 1 on the other 6 : [-1 0 0] : Yes; one edge fully on plane 8 : [-1 0 1] : Yes; one vertex on plane 2 on different sides 12 : [-1 1 1] : Yes; 2 on one side, 1 on the other 14 : [0 0 0] : No (on plane fully) 16 : [0 0 1] : Yes; one edge fully on plane 20 : [0 1 1] : No 28 : [1 1 1] : No Parameters ---------- signs: (n,3) int, all values are -1,0, or 1 Each row contains the dot product of all three vertices in a face with respect to the plane Returns --------- basic : (n,) bool Which faces are in the basic intersection case one_vertex : (n,) bool Which faces are in the one vertex case one_edge : (n,) bool Which faces are in the one edge case raxisdtypeNTF )npsortzeroslenint8rangebool)signs signs_sortedcodedikeyone_edge one_vertexbasics o/var/www/eduai.edurigo.com/storigo/production/storigo_env/lib/python3.13/site-packages/trimesh/intersections.pytriangle_cases"mesh_plane..triangle_cases6sFwwu1- \*"'':R?qA !Q$'1q50 0Ehhr&B:AAZ AQG  (**c>XS:HnXS:gRS5n[T TX$RSS9upV[R"X#UU45RS5nU$)NrF line_segmentsr-r.r)reshape plane_linesTr column_stack) rfacesvertices vertex_plane edge_thrupoint_intersectvalidlines plane_normal plane_origins r'handle_on_vertex$mesh_plane..handle_on_vertextswaZ( 1*%--g6 !, ,(=U" u*=!> PQYY   r*c:XS:HRS5nX#nU$)Nrr,)r2)rr6r7edgespointss r'handle_on_edge"mesh_plane..handle_on_edges%qj!))'2 r*c l>[R"USS/S9n[R"XU[R"USSS9XU[R"USSS945R S5n[ TTX$RSS9upVUR5(deUR S 5$) Nr-r)uniquer r.r,Fr/r1) runique_value_in_rowrr5rollr2r3r4all) rr6r7unique_elementrB intersectionsr;r=r>s r' handle_basic mesh_plane..handle_basics!55eRGL%bggnaa89%bggnaa89    ''   + ,(9  yy{{{$$Z00r*rrz%Plane origin and normal must be (3,)!r-rrr")r asanyarrayfloat64shape ValueErrorr6int64dotr7rrrr mergevstackziphstacknonzerorkind)meshr=r> return_faces local_faces cached_dotsr(r?rDrMr6dotsrcaseshandlerschr<indexs `` r' mesh_planerfsJ<+|  1(==RZZ@L==RZZ@LT!\%7%74%?@AA mmKrxx@  ;' vvdmml2LA HHS'rww 7E!E#)) E  LE 5 !E?H II:=e:NO:N$!Q58UXt}} -:NO E U;U2::a=+U;<{{3&&&  %< k%((( L P r=r]r_)r.rr-rNr.)r-r.r.)runitizerrPrQrUr7r4rplane_transformlinalginveyerfappendtftransform_pointsr2array)r\r>r=heights vertex_dotsbase_transform translation transforms face_indexsegmentsheight new_originnew_dotsr<reto_3Dto_2Dlines_2Ds r'mesh_multiplaners6<< -L==RZZ@LmmG2::6G&& (D'G'GHK--\WNYY]]>2N&&)KJJH!F%:; '!#%  #D~3 e$% &&u}}W'=uEArrE?**:6!% AF*BJJ7J  ++r*cn[R"U5n[R"U5RS5n[R"USUS- 5n[R"[R"U5RS55n[R "XUS- R 5n[R "XR 5n[R"U5[R:nU(a[R "U[R"XS- 55n[R"U5[R"U5:gn [R"[R"U5[R:[R"U5[R:5n [R"Xy5n[R"Xz5n[R"XWXg5n USUn U [R"U S5XG--n X4$)a! Calculate plane-line intersections Parameters --------- plane_origin : (3,) float Point on plane plane_normal : (3,) float Plane normal vector endpoints : (2, n, 3) float Points defining lines to be tested line_segments : bool If True, only returns intersections as valid if vertices from endpoints are on different sides of the plane. Returns --------- intersections : (m, 3) float Cartesian intersection points valid : (n, 3) bool Indicate whether a valid intersection exists for each input line segment rrrr-r)rrPr2rrlrUr4absr zero transposesign logical_or logical_anddivide) r>r= endpointsr0line_dirtbr;testdifferent_sidesrZd intersections r'r3r3%s2 i(I==.66q9L||IaL9Q<78H<< l ; C CA FGL |Yq\9<<=A |ZZ(A FF1I EvvlBLL! 1L$MN''!* 5--q CHH 4bffTlSXX6MNu6u. !(AH%AQ<&L"**Q"88?"JJL  r*c[R"U[RS9n[R"U[RS9n[R"U[RS9n[R"U[RS9nX- n[R"Xa5n[R"X15n[R "U5S:n [R "XyX5n X9U RS5-n XU - n X/n U(aU RU 5 U(aU RU5 U $)a Given one line per plane find the intersection points. Parameters ----------- plane_origins : (n,3) float Point on each plane plane_normals : (n,3) float Normal vector of each plane line_origins : (n,3) float Point at origin of each line line_directions : (n,3) float Direction vector of each line return_distance : bool Return distance from origin to point also return_denom : bool Return denominator, so you can check for small values Returns ---------- on_plane : (n,3) float Points on specified planes valid : (n,) bool Did plane intersect line or not distance : (n,) float [OPTIONAL] Distance from point denom : (n,) float [OPTIONAL] Denominator rgh㈵>r) rrPrQr diagonal_dotrrr2rq) plane_origins plane_normals line_originsline_directionsreturn_distance return_denomorigin_vectorsprojection_oriprojection_dirr;distanceon_planeresults r' planes_linesrXsNMM-rzzBMMM-rzzBM==RZZ@LmmO2::FO#1N&&~EN&&FN FF> "T )Eyy.0EFH%(8(8(AAH U##H F h n% Mr*c 2[U5S:XaXU4$USLnUbXnUbUnO[R"X- U5n[R"[U5[RS9n SX[ R *:'SX[ R :'Xn U RS[RS9n [R"U 5RS[RS9n [R"U S:[R"U 5S:*5n X*:Hn U S:HnUR5(aC[R"XU5unn[R"X5nUX'UX'US:X'XnXnUU nU[R"XS:5nU[R"XS:5nU [R"XS:5nU [R"XS:5n[U5[U5-S:Xa[U5S:Xav[R"S [RS9[R"S [RS9U(a&[R"S [RS94nU$S4nU$[R "UR#S5[U5S S 9unnUUnUR#S 5nU(aUUOSnUUU4$Un[R$"USSS9U- n UU- RU5n![R"U U5n"SU"U"S:H'[R&"U!U"5n#U#SS2SS2S4U -U-n$Un%[R"S 5n&[R"S 5n'U$U S:U SS2SS24n([U(5n)U)S:Gas[R("US:H5Sn*U[R*"[-U)5[-U)54SS9[R*"U*S-S-U*S-S-4SS94n+[R."U+[R0"[U%5[U%5SU)--5R#U)S5SS9n,U([R*"[-U)5[-U)54SS9[R*"U*S-S-R2U*R24SS9SS24R#SU)-S5n&[R."U%U&SS9n%[4R6"U,5n-[R."UU-SS9nU$U S:U SS2SS24n.[U.5n/U/S:Ga.[R("US:H5Sn0U[-U/5U04R#U/S5n1[R."U1[R0"[U%5[U%5SU/--5R#U/S5SS9n2U.[R*"[-U/5[-U/54SS9[R*"U0R2U0S-S-R24SS9SS24R#SU/-S5n'[R."U%U'SS9n%[R."UU2SS9n[R "UR#S5[U%5S S 9unnU%UnUR#S 5nSnU(a[9[R:"UUSSS9U&5n3[9[R:"UUSSS9U'5n4[R<"U3U4/5n5[R<"UUUU/5n6[R>"S[R:"U6SSS9U55n7[R<"UU7/5UnUUU4$)a Slice a mesh (given as a set of faces and vertices) with a plane, returning a new mesh (again as a set of faces and vertices) that is the portion of the original mesh to the positive normal side of the plane. Parameters --------- vertices : (n, 3) float Vertices of source mesh to slice faces : (n, 3) int Faces of source mesh to slice plane_normal : (3,) float Normal vector of plane to intersect with mesh plane_origin : (3,) float Point on plane to intersect with mesh uv : (n, 2) float, optional UV coordinates of source mesh to slice face_index : ((m,) int) Indexes of faces to slice. When no mask is provided, the default is to slice all faces. cached_dots : (n, 3) float If an external function has stored dot products pass them here to avoid recomputing Returns ---------- new_vertices : (n, 3) float Vertices of sliced mesh new_faces : (n, 3) int Faces of sliced mesh new_uv : (n, 2) int or None UV coordinates of sliced mesh rNrrr-)rrr.g)rr)rr.T) minlengthreturn_inverserkr g-q=rz ijk,ij->ik) rrrUrrr rVsumrranytmnormalsrQrTrunique_bincountr2rIrwherestackrrqaranger4rtriangulate_quadsr repeat concatenateeinsum)8r7r6r=r>uvrzr_have_uvr`r signs_sum signs_asumonedgeinsidercheckr; dot_check new_facesr cut_trianglescut_faces_quad cut_faces_tricut_signs_quad cut_signs_triemptyrGinverse final_vert final_facefinal_uvornumdenomdist int_points new_verticesnew_quad_verticesnew_tri_verticesquad_int_points num_quads quad_int_indsquad_int_vertsnew_quad_facesnew_tri_faces_from_quadstri_int_pointsnum_tris tri_int_inds tri_int_verts new_tri_facesquad_barycentricstri_barycentricsall_barycentricscut_uvnew_uvs8 r'slice_faces_planers5V 8}""nG! vvh-|< HHS]"'' 2E E#))  E  LE  q 0I"""9J ^^J!ORVVI->!-C DF + %F QH||~~zz(?";< uFF5/  "$s? I If%M2>>&a-@AN"..a@AM2>>&a-@AN"..a@AM >S//14 y>Q rzz2rxx06=rzz2E LDHE L#22   b !S]4  f% __W- !(2f:d:x// A 2A"A !   .C FF1l #EE%3, 99S% DaDj!A%)JL(xx'!)a-!8!Q!>?OO$I1}1!45a8 ' HHeI&i(89 B HH}q(A- 0AQ/FGa P R    IIc,'\):Q])J K S S1   , HHeI&i(89 B HH )Q.11=??C! L    '!i- # yy/@qI #+#=#=n#M IIi)AJ  a 8!Q >?N>"H!|xx  34Q7 %eHo|&CDLLXWXY   IIc,'\):Q\)I J R R!   * HHeHouX7a @ HHlnn q(8A'=&@&@A J    '!h, " yy/?aH IIiQ? .."\):4OFG f%J)JH1 IIh~. :>+<>N*OPN!3R 5F GH<611)EGWX>>2v,/7 z8 ++r*c j UcgSSKJn SSKJn SSKJn SSKJn SSKJ n [R"U[RS 9n[R"U[RS 9nURS :H=(d [R"US 5=(aM URS :H=(d [R"US 5=(a URUR:Hn U (d [!S 5eUR"R%5n UR&R%5n[)URS 5=(aB [R"URR*5[-UR"5S4:H=(a U(+nU(a$URR*R%5OSnSU;aSUS'[/UR1S 5UR1S 55GHtunn[3U UUUUUS9upnU(dM!U(a [5S5e[6R8"U 5unnU Un UUnUUSS2SS24USS2SS24:gR;SS9n[<R>"UU*S9n[R@RCU5n[DRF"U U5n[<RH"U5nURKSS9 [RL"USS2S45S:nUUUR;SS9nUUSS2S4USS2S4:gn[6RN"USS9n[-U5S:aGM`U"U 5nU/nU RQUUUSS2SS245HnU "UUSS9unn[DRF"[RR"U5U5n URUU 5un!n"U!RW5S:a[RXR[S5 U"Un#U#SS2SS24U#SS2SS24:gR;SS9U#SS2S4U#SS2S4:g-n$UR]U#U$5 M [R^"U5nGMw U(a)U "UURR`R%5S9OSn%U"SXU%S.UD6$)a@ Slice a mesh with a plane returning a new mesh that is the portion of the original mesh to the positive normal side of the plane. Parameters --------- mesh : Trimesh object Source mesh to slice plane_normal : (3,) float Normal vector of plane to intersect with mesh plane_origin : (3,) float Point on plane to intersect with mesh cap : bool If True, cap the result with a triangulated polygon face_index : ((m,) int) Indexes of mesh.faces to slice. When no mask is provided, the default is to slice all faces. cached_dots : (n, 3) float If an external function has stored dot products pass them here to avoid recomputing engine : None or str Triangulation engine passed to `triangulate_polygon` kwargs : dict Passed to the newly created sliced mesh Returns ---------- new_mesh : Trimesh object Sliced mesh Nr)cKDTreer)Trimesh)triangulate_polygon)polygons)TextureVisualsrrOrkz)plane origins and normals must be (n, 3)!rr.processF)r7r6rr=r>rzz)face_index and cap can't be used togetherr rhg:0yE>) require_countrT)engineforce_verticesz%triangulate may have inserted vertex!)rmaterial)r7r6visual)1 scipy.spatialrbasercreationrpathrrrrrPrQrRris_shaperSr7copyr6hasattrrrrXr2rNotImplementedErrorr unique_rowsrJrrmrnrorrrsfaces_to_edgesrr group_rowsedges_to_polygonsstack_3DquerymaxlogdebugrqrWr)&r\r=r>rzcaprkwargsrrrrrshape_okr7r6has_uvrrirjrGrfrr vertices_2DrBr unique_edgetreepvnfnvn3rvidnfnf_okrs& r'slice_mesh_planer s<R |&-&==RZZ@L==RZZ@L   t # Kt}}\7'K 5   4 ' O4==w+O 5   ,"4"4 4 DEE}}!!#H JJOO E  T"Zrxx '?C DVXYCZ'Z' #)   dB!yW%|';';G'D0!   3)*UVV'228K++A.E JJAJ vvk!Q$/047H(5/--1-56E%1+q!t45E"--e1EK6{Q8$DCE//k0BKPQSUTUSUPUDVW,QvdSB))$--*;UC $ 3 #<<>D(HHNN#JKWAqrEbBQBi/44!4<1a4BqRStH@TU RY'XIIe$Ew|HN"t{{';';'@'@'BCSW  KH& KF KKr*)FNN)T)FF)NNN)NFN)__doc__numpyrrrrrrrrr constantsr r rfrr3rrr rr*r'rss&&#, zzT,n0p @P v,z  QLr*