三维线面交点

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

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

 

    原文作者:年少的我
    原文地址: https://blog.csdn.net/DDremote/article/details/85634002
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞