X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FTools%2FblocFissure%2Fgmu%2FprolongeWire.py;h=52b3991b4540059b9314027dc2f47b9deec52f28;hb=c98d9fcd7f02c1f1f5c24dd3e709ed75228d66c4;hp=72e766dffec6d199cfe0c2e67701314331107445;hpb=072a73120b6db7bba2389aa7ada0cde20e22ee57;p=modules%2Fsmesh.git diff --git a/src/Tools/blocFissure/gmu/prolongeWire.py b/src/Tools/blocFissure/gmu/prolongeWire.py index 72e766dff..52b3991b4 100644 --- a/src/Tools/blocFissure/gmu/prolongeWire.py +++ b/src/Tools/blocFissure/gmu/prolongeWire.py @@ -2,6 +2,10 @@ import logging from geomsmesh import geompy +from geomsmesh import geomPublish +from geomsmesh import geomPublishInFather +import initLog +from orderEdgesFromWire import orderEdgesFromWire # ----------------------------------------------------------------------------- # --- prolongation d'un wire par deux segments tangents @@ -12,18 +16,43 @@ def prolongeWire(aWire, extrem, norms, long): logging.info("start") if geompy.NumberOfEdges(aWire) > 1: edges = geompy.ExtractShapes(aWire, geompy.ShapeType["EDGE"]) + uneSeuleEdge = False else: edges = [aWire] + uneSeuleEdge = True + edgesBout = [] for i, v1 in enumerate(extrem): exts = [geompy.MakeTranslationVectorDistance(v1, norms[i], l) for l in (-long, long)] dists = [(geompy.MinDistance(v, aWire), i , v) for i, v in enumerate(exts)] dists.sort() v2 = dists[-1][-1] - #v2 = geompy.MakeTranslationVectorDistance(v1, norms[i], long) edge = geompy.MakeEdge(v1, v2) edges.append(edge) + edgesBout.append(edge) name = "extrem%d"%i - #geompy.addToStudy(edge,name) - wireProlonge = geompy.MakeWire(edges) - geompy.addToStudy(wireProlonge, "wireProlonge") + geomPublish(initLog.debug, edge, name) + try: + wireProlonge = geompy.MakeWire(edges) + geomPublish(initLog.debug, wireProlonge, "wireProlonge") + except: + logging.warning("probleme MakeWire, approche pas a pas") + if uneSeuleEdge: + edgelist = [aWire] + accessList = [0] + else: + edgelist, accessList = orderEdgesFromWire(aWire) + edge1 = edgelist[accessList[0]] + if geompy.MinDistance(edgesBout[0], edge1) < 1.e-4 : + i0 = 0 + i1 = 1 + else: + i0 = 1 + i1 = 0 + wireProlonge = edgesBout[i0] + for i in range(len(edgelist)): + wireProlonge = geompy.MakeWire([wireProlonge, edgelist[accessList[i]]]) + geomPublish(initLog.debug, wireProlonge, "wireProlonge_%d"%i) + wireProlonge = geompy.MakeWire([wireProlonge,edgesBout[i1]]) + geomPublish(initLog.debug, wireProlonge, "wireProlonge") + logging.warning("prolongation wire pas a pas OK") return wireProlonge