在三维空间中,如果直线不与平面平行,将存在交点。线面交点有两种方法可以求出,我所用的是第二种方法:向量法。具体代码如下:
def line_intersection(P1, P2, coeffi):
#P1、P2为存贮平面外两点xyz坐标的列表,coeffi为存贮平面法向量a、b、c、d(常数项)的列表ax+by+cZ+d=0
P12 = [P2[0]-P1[0],P2[1]-P1[1],P2[2]-P1[2]]
i = 0
num,den,n = 0,0,0
num = coeffi[0]*P1[0]+coeffi[1]*P1[1]+coeffi[2]*P1[2] + coeffi[3]
den = coeffi[0]*P1P2[0]+coeffi[1]*P1P2[1]+coeffi[2]*P1P2[2]
if(abs(den)<1e-4):
#线面平行
flage = False
n=num/den
P=[]
for i in range(3):
P.append(P1[i]+n*P1P2[i])
return P