def discretize_edge(edge, deflection=0.1):
curve_adaptator = BRepAdaptor_Curve(edge)
discretizer = GCPnts_QuasiUniformDeflection()
discretizer.Initialize(
curve_adaptator, deflection, curve_adaptator.FirstParameter(), curve_adaptator.LastParameter()
)
if not discretizer.IsDone():
raise AssertionError("Discretizer not done.")
points = [curve_adaptator.Value(discretizer.Parameter(i)).Coord() for i in range(1, discretizer.NbPoints() + 1)]
# return tuples representing the single lines of the egde
edges = []
for i in range(len(points) - 1):
edges.append((points[i], points[i + 1]))
return np.asarray(edges, dtype=np.float32)
def _transform_curve(edge, tolerance: float = 1e-3):
"""Converts a curved edge into a series of straight lines (facetets) with
the provided tolerance.
Args:
edge (cadquery.Wire): The CadQuery wire to redraw as a series of
straight lines (facet)
tolerance: faceting toleranceto use when faceting cirles and
splines. Defaults to 1e-3.
Returns:
cadquery.Wire
"""
curve = edge._geomAdaptor() # adapt the edge into curve
start = curve.FirstParameter()
end = curve.LastParameter()
points = GCPnts_QuasiUniformDeflection(curve, tolerance, start, end)
verts = (cq.Vector(points.Value(i + 1)) for i in range(points.NbPoints()))
return cq.Wire.makePolygon(verts)
def makeSVGedge(e):
"""
"""
cs = StringIO.StringIO()
curve = e._geomAdaptor() # adapt the edge into curve
start = curve.FirstParameter()
end = curve.LastParameter()
points = GCPnts_QuasiUniformDeflection(curve, DISCRETIZATION_TOLERANCE,
start, end)
if points.IsDone():
point_it = (points.Value(i + 1) for i in range(points.NbPoints()))
p = next(point_it)
cs.write("M{},{} ".format(p.X(), p.Y()))
for p in point_it:
cs.write("L{},{} ".format(p.X(), p.Y()))
return cs.getvalue()
def discretize_edge(a_topods_edge,
deflection=0.1,
algorithm="QuasiUniformDeflection"):
"""Take a TopoDS_Edge and returns a list of points
The more deflection is small, the more the discretization is precise,
i.e. the more points you get in the returned points
algorithm: to choose in ["UniformAbscissa", "QuasiUniformDeflection"]
"""
if not is_edge(a_topods_edge):
raise AssertionError(
"You must provide a TopoDS_Edge to the discretize_edge function.")
if a_topods_edge.IsNull():
print(
"Warning : TopoDS_Edge is null. discretize_edge will return an empty list of points."
)
return []
curve_adaptator = BRepAdaptor_Curve(a_topods_edge)
first = curve_adaptator.FirstParameter()
last = curve_adaptator.LastParameter()
if algorithm == "QuasiUniformDeflection":
discretizer = GCPnts_QuasiUniformDeflection()
elif algorithm == "UniformAbscissa":
discretizer = GCPnts_UniformAbscissa()
elif algorithm == "UniformDeflection":
discretizer = GCPnts_UniformDeflection()
else:
raise AssertionError("Unknown algorithm")
discretizer.Initialize(curve_adaptator, deflection, first, last)
if not discretizer.IsDone():
raise AssertionError("Discretizer not done.")
if not discretizer.NbPoints() > 0:
raise AssertionError("Discretizer nb points not > 0.")
points = []
for i in range(1, discretizer.NbPoints() + 1):
p = curve_adaptator.Value(discretizer.Parameter(i))
points.append(p.Coord())
return points