July 27, 2022, 10:01am
What’s the logic behind letting lines intersect what seems like each and every circle in a scene?
For instance, the line in green intersects both circles 0 and 1, although it should only intersect circle 0 when interpreted as an infinite line.
However, it finds an intersection with circle 1 as well (cf. GH panel)?
The same goes for “physical” line circle intersections?
Here’s my script:
import Rhino.Geometry as rg
for i in range(len(C)):
rc, t1, _, t2, _ = rg.Intersect.Intersection.LineCircle(L, C[i])
print "Circle", i, ":", t1, t2, rc
Is this a bug?
line-circle-intersection-ghpyhton.gh (6.8 KB)
July 27, 2022, 11:22am
looks like a bug to me.
the Intersection returns the closest point from the circle to the line.
import Rhino.Geometry as rg
pts = 
rc, t1, p1, t2, p2 = rg.Intersect.Intersection.LineCircle(L, C)
if int(rc) > 0 and (L.DistanceTo(p1,False) < 0.001):
if int(rc) > 1 and (L.DistanceTo(p2,False) < 0.001):
print "Circle", ":", t1, t2, rc, bool(p1), bool(p2)
( i changed to item-Access and added output pts of the scripting component)
you should always check rc:
Single is returned, only t1 and point1 will have valid values. If Multiple is returned, t2 and point2 will also be filled out.
check the enum “rc” and the distance to the line as workarround
line-circle-intersection-ghpyhton_tp02.gh (10.7 KB)
July 27, 2022, 12:21pm
Yes, I made a little deep dive into the topic, and it seems that to find an intersection between a circle and line, you can construct a line from the center of the circle perpendicular to the line. If the intersection point of the original and the new line lies within the circle, there’s a two point intersection between the circle and the line and if it’s on the circle there’s a single point intersection.
This should be
pts.append(pi) in Python.
int(rc) > 0
This is shorter but I like the more explicit version better.
I’m aware of that, but it seems that t2 simply inherits the value of t1 when the intersection is single, and is not not filled out.
July 27, 2022, 12:41pm
dive into the topic
the planar / 2d case is described here:
An (infinite) line determined by two points (x_1,y_1) and (x_2,y_2) may intersect a circle of radius r and center (0, 0) in two imaginary points (left figure), a degenerate single point (corresponding to the line being tangent to the circle; middle...
(edit my initial point add → append… sorry i checked if the error is the same (yes) in a c# script)
July 27, 2022, 1:12pm
Here’s my workaround function, inspired by your suggestion.
It’s only meant to check for intersections, nothing more.
MTOL = sc.doc.ModelAbsoluteTolerance
def intersecting(line, circle):
relation, t1, pt1, t2, _ = rg.Intersect.Intersection.LineCircle(line, circle)
if relation == rg.Intersect.LineCircleIntersection.Multiple:
if (t1 >= 0.0 and t1 <= 1.0) and (t2 >= 0.0 and t2 <= 1.0):
elif relation == rg.Intersect.LineCircleIntersection.Single:
if line.DistanceTo(pt1, False) < MTOL:
@dale, is this indeed a bug?