Python代码实现:
import numpy as np p1= np.array([x1,y1,z1])# p2 = np.array([x2,y2,z2]) plane_normal= np.array([a,b,c])#a,b,c,d平面方程系数 P1D= (np.vdot(p1,plane_normal)+d)/np.sqrt(np.vdot(plane_normal,plane_normal)) P1D2= (np.vdot(p2-P1,plane_normal))/np.sqrt(np.vdot(plane_normal,plane_normal)) n= P1D2/P1D p= p1 + n*(p2- p1)#所求交点
原文链接:http://blog.csdn.net/smallflyingpig/article/details/51234711?locationNum=8&fps=1