判断点与线段位置关系

白糖

通过简单的向量计算判断点与线段位置关系，没有设置点到线段距离容差，所以点击线段时不会判断在线段上

2. 
   
3. [CommandMethod("test")]
   
4. public void Test()
   
5. {
   
6. Database db = HostApplicationServices.WorkingDatabase;
   
7. Editor ed = Autodesk.AutoCAD.ApplicationServices.Application.DocumentManager.MdiActiveDocument.Editor;
   
8. 
   
9. using (Transaction trans = db.TransactionManager.StartTransaction())
   
10. {
    
11.   PromptPointOptions ptOpts = new PromptPointOptions("选择第一点A");
    
12.   PromptPointResult ptRes = ed.GetPoint(ptOpts);
    
13.   if (ptRes.Status != PromptStatus.OK)
    
14.   {
    
15.    throw new Autodesk.AutoCAD.Runtime.Exception();
    
16.   }
    
17.   Point3d pt1 = ptRes.Value;
    
18.   
    
19.   ptOpts = new PromptPointOptions("\n选择第二点B");
    
20.   ptRes = ed.GetPoint(ptOpts);
    
21.   if (ptRes.Status != PromptStatus.OK)
    
22.   {
    
23.    throw new Autodesk.AutoCAD.Runtime.Exception();
    
24.   }
    
25.   Point3d pt2 = ptRes.Value;
    
26.   
    
27.   //生成AB线段
    
28.   Line lineAB = new Line(pt1,pt2);
    
29.   AppendEntity(lineAB);
    
30.   trans.Commit();
    
31.   Vector3d vec1 = pt2 - pt1;
    
32.   ed.WriteMessage(vec1.ToString()+"\n");
    
33.   
    
34.   while (true)
    
35.   {
    
36.    ptOpts = new PromptPointOptions("\n选择第三点C");
    
37.    ptRes = ed.GetPoint(ptOpts);
    
38.    if (ptRes.Status != PromptStatus.OK)
    
39.    {
    
40.     throw new Autodesk.AutoCAD.Runtime.Exception();
    
41.    }
    
42.    Point3d pt3 = ptRes.Value;
    
43.    Vector3d vec2 = pt3 - pt1;
    
44.    
    
45.    //向量叉乘
    
46.    double result = vec2.X*vec1.Y-vec1.X*vec2.Y;
    
47.    
    
48.    if (result>0)
    
49.    {
    
50.     ed.WriteMessage("C点在AB线右侧");
    
51.    }else if (result<0)
    
52.    {
    
53.     ed.WriteMessage("C点在AB线左侧");
    
54.    }else if (result==0)
    
55.    {
    
56.     ed.WriteMessage("C点在AB线上");
    
57.    }
    
58.   }
    
59. }
    
60. }
    
61. public static ObjectId AppendEntity(Entity ent)
    
62. {
    
63. Database db = HostApplicationServices.WorkingDatabase;
    
64. ObjectId entid;
    
65. using (Transaction trans = db.TransactionManager.StartTransaction())
    
66. {
    
67.   BlockTable bt = (BlockTable)trans.GetObject(db.BlockTableId,OpenMode.ForRead);
    
68.   BlockTableRecord btr = (BlockTableRecord)trans.GetObject(bt[BlockTableRecord.ModelSpace],OpenMode.ForWrite);
    
69.   entid = btr.AppendEntity(ent);
    
70.   trans.AddNewlyCreatedDBObject(ent,true);
    
71.   trans.Commit();
    
72. }
    
73. return entid;
    
74. }

<br /> [CommandMethod("test")]<br /> public void Test()<br /> {<br /> Database db = HostApplicationServices.WorkingDatabase;<br /> Editor ed = Autodesk.<a href="http://bbs.mjtd.com/forum-41-1.html" target="_blank" class="relatedlink">AutoCAD</a>.ApplicationServices.Application.DocumentManager.MdiActiveDocument.Editor;<br /> <br /> using (Transaction trans = db.TransactionManager.StartTransaction())<br /> {<br />   PromptPointOptions ptOpts = new PromptPointOptions("选择第一点A");<br />   PromptPointResult ptRes = ed.GetPoint(ptOpts);<br />   if (ptRes.Status != PromptStatus.OK)<br />   {<br />    throw new Autodesk.AutoCAD.Runtime.Exception();<br />   }<br />   Point3d pt1 = ptRes.Value;<br />   <br />   ptOpts = new PromptPointOptions("\n选择第二点B");<br />   ptRes = ed.GetPoint(ptOpts);<br />   if (ptRes.Status != PromptStatus.OK)<br />   {<br />    throw new Autodesk.AutoCAD.Runtime.Exception();<br />   }<br />   Point3d pt2 = ptRes.Value;<br />   <br />   //生成AB线段<br />   Line lineAB = new Line(pt1,pt2);<br />   AppendEntity(lineAB);<br />   trans.Commit();<br />   Vector3d vec1 = pt2 - pt1;<br />   ed.WriteMessage(vec1.ToString()+"\n");<br />   <br />   while (true)<br />   {<br />    ptOpts = new PromptPointOptions("\n选择第三点C");<br />    ptRes = ed.GetPoint(ptOpts);<br />    if (ptRes.Status != PromptStatus.OK)<br />    {<br />     throw new Autodesk.AutoCAD.Runtime.Exception();<br />    }<br />    Point3d pt3 = ptRes.Value;<br />    Vector3d vec2 = pt3 - pt1;<br />    <br />    //向量叉乘<br />    double result = vec2.X*vec1.Y-vec1.X*vec2.Y;<br />    <br />    if (result>0)<br />    {<br />     ed.WriteMessage("C点在AB线右侧");<br />    }else if (result<0)<br />    {<br />     ed.WriteMessage("C点在AB线左侧");<br />    }else if (result==0)<br />    {<br />     ed.WriteMessage("C点在AB线上");<br />    }<br />   }<br /> }<br /> }<br /> public static ObjectId AppendEntity(Entity ent)<br /> {<br /> Database db = HostApplicationServices.WorkingDatabase;<br /> ObjectId entid;<br /> using (Transaction trans = db.TransactionManager.StartTransaction())<br /> {<br />   BlockTable bt = (BlockTable)trans.GetObject(db.BlockTableId,OpenMode.ForRead);<br />   BlockTableRecord btr = (BlockTableRecord)trans.GetObject(bt[BlockTableRecord.ModelSpace],OpenMode.ForWrite);<br />   entid = btr.AppendEntity(ent);<br />   trans.AddNewlyCreatedDBObject(ent,true);<br />   trans.Commit();<br /> }<br /> return entid;<br /> }

cdinten

语句：
double result = vec2.X*vec1.Y-vec1.X*vec2.Y;
考虑Vector3d结构体的方法CrossProduct 之后向量的Length==0，
另外，还可以使用IsParallelTo方法，看向量AC和向量CB是不是平行的，一旦平行必定共线。
思想不错，但是你抛出异常不处理让我感觉不太好，而且一旦半路取消就会有异常抛出。

cdinten

另外，还要啰嗦一下，你的代码和我说的方法都对于点在线段的延长线和在线段上是区分不出来的，呵呵
你可以使用Length属性

cdinten

看样子你是新学习这个的吧，觉得有时间可以考虑一下我说的意见，当个练习，祝你成功

白糖

cdinten 发表于 2012-12-7 22:42
语句：
double result = vec2.X*vec1.Y-vec1.X*vec2.Y;
考虑Vector3d结构体的方法CrossProduct 之后向量的 ...

谢谢提醒，还是不够严谨~
但是我装的是AutoCAD 2007，没有叉乘方法~
共线判断确实可行啊！

cdinten

白糖 发表于 2012-12-7 22:52
谢谢提醒，还是不够严谨~
但是我装的是AutoCAD 2007，没有叉乘方法~
共线判断确实可行啊！

哦，对的，我没看仔细，不过水平线怎么分左右侧？

白糖

cdinten 发表于 2012-12-7 22:55
哦，对的，我没看仔细，不过水平线怎么分左右侧？

我个人理解的左右，其实是指的顺时针方向的点在右侧，逆时针方向的在左侧
这样说不知道你是否认同~

cdinten

白糖 发表于 2012-12-7 23:07
我个人理解的左右，其实是指的顺时针方向的点在右侧，逆时针方向的在左侧
这样说不知道你是否认同~

嗯，相当于角度的正负方向

mkhsj928

左右之分不是以顺时针和逆时针来分，而是以曲线的走向来分。