直线与平面交点

eachy

写了个 LispFunction ，很简单，
(linetoplane p1 p2 p3 p4 p5）
p1 p2 为线上两点， p3 p4 p5 为i构成平面的非共线三点

2. 
   
3. // (C) Copyright 2014 by  
   
4. using System;
   
5. using Autodesk.AutoCAD.Runtime;
   
6. using Autodesk.AutoCAD.ApplicationServices;
   
7. using Autodesk.AutoCAD.DatabaseServices;
   
8. using Autodesk.AutoCAD.Geometry;
   
9. using Autodesk.AutoCAD.EditorInput;
   
10. using Exception = System.Exception;
    
11. 
    
12. [assembly: CommandClass(typeof(AutoCAD_CSharp_plug_in34.MyCommands))]
    
13. 
    
14. namespace AutoCAD_CSharp_plug_in34
    
15. {
    
16.     public class MyCommands
    
17.     {
    
18. 
    
19.         [LispFunction("LineToPlane")]
    
20.         public object MyLispFunction(ResultBuffer args) // This method can have any name
    
21.         {
    
22.             TypedValue[] vars = args.AsArray();
    
23.             if (vars != null && vars.Length == 5
    
24.                 && vars[0].TypeCode == (int)LispDataType.Point3d
    
25.                 && vars[1].TypeCode == (int)LispDataType.Point3d
    
26.                 && vars[2].TypeCode == (int)LispDataType.Point3d
    
27.                 && vars[3].TypeCode == (int)LispDataType.Point3d
    
28.                 && vars[4].TypeCode == (int)LispDataType.Point3d)
    
29.             {
    
30.                 Document doc = Application.DocumentManager.MdiActiveDocument;
    
31.                 Database db = doc.Database;
    
32.                 Transaction tr = db.TransactionManager.StartTransaction();
    
33.                 using (tr)
    
34.                 {
    
35.                     try
    
36.                     {
    
37.                         Point3d p1 = (Point3d)vars[0].Value;
    
38.                         Point3d p2 = (Point3d)vars[1].Value;
    
39.                         Point3d p3 = (Point3d)vars[2].Value;
    
40.                         Point3d p4 = (Point3d)vars[3].Value;
    
41.                         Point3d p5 = (Point3d)vars[4].Value;
    
42.                         Line3d ln = new Line3d(p1, p2);
    
43.                         Vector3d vv = p4 - p3;
    
44.                         Vector3d vu = p5 - p3;
    
45.                         Plane plan = new Plane(p3, vv.CrossProduct(vu));
    
46.                         if (ln.IsParallelTo(plan))
    
47.                             return null;
    
48.                         else
    
49.                         {
    
50.                             Point3d[] pnt = ln.IntersectWith(plan);
    
51.                             if (pnt != null)
    
52.                             {
    
53.                                 return pnt[0];
    
54.                             }
    
55.                             else
    
56.                             {
    
57.                                 return null;
    
58.                             }
    
59.                         }
    
60.                         tr.Commit();
    
61.                     }
    
62.                     catch (Exception)
    
63.                     {
    
64.                         throw;
    
65.                     }
    
66.                     
    
67.                 }
    
68.             }
    
69.             else
    
70.             {
    
71.                 return null;
    
72.             }
    
73.         }
    
74.     }
    
75. }
    

<br /> // (C) Copyright 2014 by  <br /> using System;<br /> using Autodesk.AutoCAD.Runtime;<br /> using Autodesk.AutoCAD.ApplicationServices;<br /> using Autodesk.AutoCAD.DatabaseServices;<br /> using Autodesk.AutoCAD.Geometry;<br /> using Autodesk.AutoCAD.EditorInput;<br /> using Exception = System.Exception;<br /> <br /> [assembly: CommandClass(typeof(AutoCAD_CSharp_plug_in34.MyCommands))]<br /> <br /> namespace AutoCAD_CSharp_plug_in34<br /> {<br />     public class MyCommands<br />     {<br /> <br />         [LispFunction("LineToPlane")]<br />         public object MyLispFunction(ResultBuffer args) // This method can have any name<br />         {<br />             TypedValue[] vars = args.AsArray();<br />             if (vars != null && vars.Length == 5<br />                 && vars[0].TypeCode == (int)LispDataType.Point3d<br />                 && vars[1].TypeCode == (int)LispDataType.Point3d<br />                 && vars[2].TypeCode == (int)LispDataType.Point3d<br />                 && vars[3].TypeCode == (int)LispDataType.Point3d<br />                 && vars[4].TypeCode == (int)LispDataType.Point3d)<br />             {<br />                 Document doc = Application.DocumentManager.MdiActiveDocument;<br />                 Database db = doc.Database;<br />                 Transaction tr = db.TransactionManager.StartTransaction();<br />                 using (tr)<br />                 {<br />                     try<br />                     {<br />                         Point3d p1 = (Point3d)vars[0].Value;<br />                         Point3d p2 = (Point3d)vars[1].Value;<br />                         Point3d p3 = (Point3d)vars[2].Value;<br />                         Point3d p4 = (Point3d)vars[3].Value;<br />                         Point3d p5 = (Point3d)vars[4].Value;<br />                         Line3d ln = new Line3d(p1, p2);<br />                         Vector3d vv = p4 - p3;<br />                         Vector3d vu = p5 - p3;<br />                         Plane plan = new Plane(p3, vv.CrossProduct(vu));<br />                         if (ln.IsParallelTo(plan))<br />                             return null;<br />                         else<br />                         {<br />                             Point3d[] pnt = ln.IntersectWith(plan);<br />                             if (pnt != null)<br />                             {<br />                                 return pnt[0];<br />                             }<br />                             else<br />                             {<br />                                 return null;<br />                             }<br />                         }<br />                         tr.Commit();<br />                     }<br />                     catch (Exception)<br />                     {<br />                         throw;<br />                     }<br />                     <br />                 }<br />             }<br />             else<br />             {<br />                 return null;<br />             }<br />         }<br />     }<br /> }<br />

相同功能的 XdGe 函数实现如下

1. (defun linetoPlane (p1 p2 p3 p4 p5 / line3d plane pnt)
   
2.   (and (setq line3d (xdge::constructor "kLine3d" p1 p2))
   
3.        (setq plane (xdge::constructor "kPlane" p3 p4 p5))
   
4.        (setq pnt (xdge::getpropertyvalue line3d "intersectWith" plane))
   
5.        (xdge::free)
   
6.   )
   
7.   pnt
   
8. )

(defun linetoPlane (p1 p2 p3 p4 p5 / line3d plane pnt)<br />   (and (setq line3d (xdge::constructor "kLine3d" p1 p2))<br />        (setq plane (xdge::constructor "kPlane" p3 p4 p5))<br />        (setq pnt (xdge::getpropertyvalue line3d "intersectWith" plane))<br />        (xdge::free)<br />   )<br />   pnt<br /> )

测试代码，修改下点样式，不然可能不容易发现

1. (defun c:tt (/ p1 p2 p3 p4 p5 p)
   
2.   (if (and (setq p1 (getpoint "\n线上地一点: "))
   
3.            (setq p2 (getpoint p1 "\n线上第二点: "))
   
4.            (setq p3 (getpoint "\n平面上第一点: "))
   
5.            (setq p4 (getpoint p3 "\n平面上第二点: "))
   
6.            (setq p5 (getpoint p3 "\n平面上第三点: "))
   
7.            (setq p (linetoplane p1 p2 p3 p4 p5))
   
8.       )
   
9.     (entmake (list '(0 . "point") (cons 10 p)))
   
10.   )
    
11.   (princ)
    
12. )

(defun c:tt (/ p1 p2 p3 p4 p5 p)<br />   (if (and (setq p1 (getpoint "\n线上地一点: "))<br />            (setq p2 (getpoint p1 "\n线上第二点: "))<br />            (setq p3 (getpoint "\n平面上第一点: "))<br />            (setq p4 (getpoint p3 "\n平面上第二点: "))<br />            (setq p5 (getpoint p3 "\n平面上第三点: "))<br />            (setq p (linetoplane p1 p2 p3 p4 p5))<br />       )<br />     (entmake (list '(0 . "point") (cons 10 p)))<br />   )<br />   (princ)<br /> )

csharp

平面三点的合法性应判断

st788796

还可以利用 Calc 计算器计算

命令: (setq p1 (getpoint))
(3546.4 79.8683 823.661)

命令: (setq p2 (getpoint))
(4140.34 561.425 0.0)

命令: (setq p3 (getpoint))
(3546.4 561.425 0.0)

命令: (setq p4 (getpoint))
(4140.34 79.8683 0.0)

命令: (setq p5 (getpoint))
(4140.34 561.425 411.831)

命令: cal
>> 表达式: ilp(p1,p2,p3,p4,p5)
3991.85795,441.035465,205.915305

命令: (c:cal "ilp(p1,p2,p3,p4,p5)")
(3991.86 441.035 205.915)

ill 和 ilp 函数可以确定交点。

ill(p1,p2,p3,p4)确定两条直线 (p1、p2) 和 (p3、p4) 的交点。所有点都被认为是三维点。

ilp(p1,p2,p3,p4,p5)确定直线（通过 p1、p2）和平面（通过三点 p3、p4、p5）的交点。