Hi Uwe,
In pythonOCC an edge has an underlying geometry: straight line, parabola,
circle, ellipse, or spline, to name a few. By a segment, I meant a straight
line geometry bounded by two vertices. BRepExtrema_DistShapeShape is
written for arbitrary shapes. In the case of the distance between a vertex
and an edge, it works for edges of any geometry. If all you need is the
distance between a straight line and a vertex, BRepExtrema_DistShapeShape is
overkill, and it will be slower than a more specialized approach. So, are
your edges straight lines bounded by two vertices or something else (like a
spline)? If they're straight lines, I currently lack the time to create an
example. However, I'd simply be putting the formula at:
http://mathworld.wolfram.com/Point-LineDistance3-Dimensional.html into
python numpy. If you need help coding it up, let me know, and I can try to
help later next week.
If I understand correctly, your current approach already pre-filters distant
edges with the midpoint/KDTree search, so my bounding box advice is
superfluous.
If you're crashing for a small set of edges (<100k), I think
BRepExtrema_DistShapeShape has a memory leak. If that's the case, we need
to find whether it's a pythonOCC bug or an OCC bug and notify the
developers.
If it is a leak, have you tried not destroying the
BRepExtrema_DistShapeShape instance every time in the loop as a workaround?
The pertinent lines in your example would look like this:
b1 = BRepExtrema_DistShapeShape()
for k in pntDic:
for l in curveDic:
v = BRepBuilderAPI_MakeVertex(pntDic[k])
b1.LoadS1(v.Vertex())
b1.LoadS2(curveDic[l])
b1.Perform()
dist = b1.Value()
All that said, I don't understand Lie Pablo's response. Perhaps I'm missing
something, and I'm not really helping.
- Charles
_______________________________________________
Pythonocc-users mailing list
Pythonocc-users@gna.org
https://mail.gna.org/listinfo/pythonocc-users
