1 # -*- coding: utf-8 -*-
5 from .geomsmesh import geompy
6 from .geomsmesh import geomPublish
7 from .geomsmesh import geomPublishInFather
10 from .sortFaces import sortFaces
11 from .extractionOrientee import extractionOrientee
13 def identifieFacesPeau(ifil, verticesPipePeau, facesOnside, wireFondFiss,
14 verticesEdgesFondIn, pipexts, cercles,
15 fillingFaceExterne, centreFondFiss):
17 inventaire des faces de peau : face de peau percee du pipe, extremites du pipe
18 La partition avec le pipe peut creer un vertex (et un edge) de trop sur le cercle projete,
19 quand le cercle est tres proche de la face.
20 dans ce cas, la projection du cercle sur la face suivie d'une partition permet
21 d'eviter le point en trop
25 facesAndFond = facesOnside
26 facesAndFond.append(wireFondFiss)
28 partitionPeauByPipe = geompy.MakePartition(facesAndFond, pipexts, [], [], geompy.ShapeType["FACE"], 0, [], 1)
30 logging.debug("probleme partition face pipe, contournement avec MakeSection")
32 for pipext in pipexts:
33 sections.append(geompy.MakeSection(facesOnside[0], pipext))
34 partitionPeauByPipe = geompy.MakePartition(facesAndFond, sections, [], [], geompy.ShapeType["FACE"], 0, [], 1)
36 # contrle edge en trop sur edges circulaires
37 if len(verticesPipePeau) > 0: # --- au moins une extremite du pipe sur cette face de peau
40 facesPeau = geompy.ExtractShapes(partitionPeauByPipe, geompy.ShapeType["FACE"], False)
41 facesPeauSorted, minsur, maxsurf = sortFaces(facesPeau)
42 for i, face in enumerate(facesPeauSorted[:-1]): # on ne teste que la ou les petites faces "circulaires"
43 nbv = geompy.NumberOfEdges(face)
44 logging.debug("nombre d'edges sur face circulaire: %s", nbv)
46 edgeEnTrop.append(True) # TODO : distinguer les cas avec deux faces circulaires dont l'une est correcte
48 edgeEnTrop.append(False)
49 refaire = sum(edgeEnTrop)
51 dc = [(geompy.MinDistance(verticesEdgesFondIn[0], fac), i) for i, fac in enumerate(facesPeauSorted[:-1])]
53 logging.debug("dc sorted: %s", dc)
54 i0 = dc[0][1] # indice de facesPeauSorted qui correspond a verticesEdgesFondIn[0], donc 0 pour cercles
56 for i, bad in enumerate(edgeEnTrop):
62 outilPart[j] = geompy.MakeProjection(cercles[j],facesOnside[0])
64 partitionPeauByPipe = geompy.MakePartition(facesAndFond, outilPart, [], [], geompy.ShapeType["FACE"], 0, [], 1)
67 name="partitionPeauByPipe%d"%ifil
68 geomPublish(initLog.debug, partitionPeauByPipe, name)
69 [edgesPeauFondIn, edgesPeauFondOut, edgesPeauFondOn] = extractionOrientee(fillingFaceExterne, partitionPeauByPipe, centreFondFiss, "EDGE", 1.e-3)
70 [facesPeauFondIn, facesPeauFondOut, facesPeauFondOn] = extractionOrientee(fillingFaceExterne, partitionPeauByPipe, centreFondFiss, "FACE", 1.e-3)
72 if len(verticesPipePeau) > 0: # --- au moins une extremite du pipe sur cette face de peau
73 facesPeauSorted, minsur, maxsurf = sortFaces(facesPeauFondOn)
74 facePeau = facesPeauSorted[-1] # la plus grande face
76 facePeau =geompy.MakePartition(facesPeauFondOn, [], [], [], geompy.ShapeType["FACE"], 0, [], 1)
77 facesPeauSorted = [facePeau]
78 name="facePeau%d"%ifil
79 geomPublish(initLog.debug, facePeau, name)
81 return (facePeau, facesPeauSorted, edgesPeauFondIn)