Added vertex filleting example

Author: tpaviot

examples/core_topology_vertex_filleting.py

@@ -0,0 +1 @@
+#!/usr/bin/env python##Copyright 2008-2015 Jelle Feringa ([email protected])####This file is part of pythonOCC.####pythonOCC is free software: you can redistribute it and/or modify##it under the terms of the GNU Lesser General Public License as published by##the Free Software Foundation, either version 3 of the License, or##(at your option) any later version.####pythonOCC is distributed in the hope that it will be useful,##but WITHOUT ANY WARRANTY; without even the implied warranty of##MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the##GNU Lesser General Public License for more details.####You should have received a copy of the GNU Lesser General Public License##along with pythonOCC.  If not, see <http://www.gnu.org/licenses/>.# A sample that shows how to generate the gear geometry according# to knowledgefrom OCC.BRepFilletAPI import BRepFilletAPI_MakeFilletfrom OCC.BRep import BRep_Toolfrom OCC.BRepPrimAPI import BRepPrimAPI_MakeBoxfrom OCC.TopExp import (TopExp_Explorer,                        topexp_MapShapesAndAncestors,                        topexp_FirstVertex,                        topexp_LastVertex)from OCC.TopAbs import TopAbs_VERTEX, TopAbs_EDGEfrom OCC.TopTools import (TopTools_IndexedDataMapOfShapeListOfShape,                          TopTools_ListIteratorOfListOfShape)from OCC.TopoDS import topods_Vertex, topods_Edgefrom OCC.Display.SimpleGui import init_displaydisplay, start_display, add_menu, add_function_to_menu = init_display()# create shapecube = BRepPrimAPI_MakeBox(100, 100, 100).Shape()topExp = TopExp_Explorer()topExp.Init(cube, TopAbs_VERTEX)# get two verticesvertA = topods_Vertex(topExp.Current())topExp.Next()vertB = topods_Vertex(topExp.Current())def vertex_fillet(cube, vert):    # apply a fillet on incident edges on a vertex    afillet = BRepFilletAPI_MakeFillet(cube)    cnt = 0    # find edges from vertex    _map = TopTools_IndexedDataMapOfShapeListOfShape()    topexp_MapShapesAndAncestors(cube, TopAbs_VERTEX, TopAbs_EDGE, _map)    results = _map.FindFromKey(vert)    topology_iterator = TopTools_ListIteratorOfListOfShape(results)    while topology_iterator.More():        edge = topods_Edge(topology_iterator.Value())        topology_iterator.Next()        first, last = topexp_FirstVertex(edge), topexp_LastVertex(edge)        vertex, first_vert, last_vert = BRep_Tool().Pnt(vert), BRep_Tool().Pnt(first), BRep_Tool().Pnt(last)        if edge.Orientation():            if not vertex.IsEqual(first_vert, 0.001):                afillet.Add(0, 20., edge)            else:                afillet.Add(20, 0, edge)        cnt += 1    afillet.Build()    if afillet.IsDone():        return afillet.Shape()    else:        raise AssertionError('you failed on me you fool!')filleted_vertA = vertex_fillet(cube, vertA)display.DisplayShape(filleted_vertA, update=True)start_display()