python 检测直线 交点_在python中,找到直线和圆的交点的最有效方法是什么?

我有一个由许多点组成的多边形。我想找到多边形和圆的交集。提供[x0,y0]的圆心和r0的半径,我写了一个粗函数来简单求解圆和直线的二次方程。但是,一个接一个地找到多边形的每个线段的交集的效率如何呢?有更有效的方法吗?

我知道sympy已经提供了获取不同几何体交点的功能。但是如果我想处理大量的多边形,那么与用我自己的函数计算相比,sympy这样的外部库的效率如何?def LineIntersectCircle(p,lsp,lep):

# p is the circle parameter, lsp and lep is the two end of the line

x0,y0,r0 = p

x1,y1 = lsp

x2,y2 = esp

if x1 == x2:

if abs(r0) >= abs(x1 - x0):

p1 = x1, y0 - sqrt(r0**2 - (x1-x0)**2)

p2 = x1, y0 + sqrt(r0**2 - (x1-x0)**2)

inp = [p1,p2]

# select the points lie on the line segment

inp = [p for p in inp if p[1]>=min(y1,y2) and p[1]<=max(y1,y2)]

else:

inp = []

else:

k = (y1 - y2)/(x1 - x2)

b0 = y1 - k*x1

a = k**2 + 1

b = 2*k*(b0 - y0) - 2*x0

c = (b0 - y0)**2 + x0**2 - r0**2

delta = b**2 - 4*a*c

if delta >= 0:

p1x = (-b - sqrt(delta))/(2*a)

p2x = (-b + sqrt(delta))/(2*a)

p1y = k*x1 + b0

p2y = k*x2 + b0

inp = [[p1x,p1y],[p2x,p2y]]

# select the points lie on the line segment

inp = [p for p in inp if p[0]>=min(x1,x2) and p[0]<=max(x1,x2)]

else:

inp = []

return inp

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值