三维线面交点

在三维空间中，如果直线不与平面平行，将存在交点。线面交点有两种方法可以求出，我所用的是第二种方法：向量法。具体代码如下：

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