拾取一个点坐标，找到包围这个坐标的最小封闭区域问题

kuaigoumanzhu

开始没一点思路，后来在论坛看了一天的帖子。现在是先拾取一个点坐标，然后找到所有封闭区域的顶点坐标然后再计算这个点是否在封闭区域内。

1. BlockTable bt = trans.GetObject(db.BlockTableId, OpenMode.ForRead) as BlockTable;
   
2.                 BlockTableRecord btr = trans.GetObject(bt[BlockTableRecord.ModelSpace], OpenMode.ForWrite) as BlockTableRecord;
   
3.                 Polyline p = new Polyline();
   
4.                 Point3d a = new Point3d();
   
5.                 PromptPointOptions opt = new PromptPointOptions("选择点");
   
6. PromptPointResult rest = ed.GetPoint(opt);
   
7. if (rest.Status == PromptStatus.OK)
   
8.                 {
   
9.                     a = rest.Value;
   
10.                 }
    
11. TypedValue[] pline ={ new TypedValue((int)DxfCode.Start, "LWPOLYLINE") };
    
12.                 SelectionFilter sf = new SelectionFilter(pline);
    
13.                 PromptSelectionResult res = ed.SelectAll(sf);
    
14.                 SelectionSet st = res.Value;
    
15.                 if (st != null)
    
16.                 {
    
17.                     ObjectId[] ids = st.GetObjectIds();
    
18.                     Point3dCollection p3dc=new Point3dCollection();
    
19.                     foreach (ObjectId id in ids)
    
20.                     {
    
21.                         p = trans.GetObject(id, OpenMode.ForRead) as Polyline;
    
22.                         //Point3d c = new Point3d((p.c));
    
23.                         int j=p.NumberOfVertices-1;
    
24.                         bool result = false;
    
25.                         for (int i = 0; i < p.NumberOfVertices; i++)
    
26.                         {
    
27.                             p3dc.Add(p.GetPoint3dAt(i));
    
28.                             //nodes.Add(p.GetPoint3dAt(i));
    
29.                             ed.WriteMessage(p.GetPoint3dAt(i).ToString());
    
30.                            
    
31.                             if (((p.GetPoint3dAt(i).Y)<=a.Y)&&(a.Y<p.GetPoint3dAt(j).Y)||((p.GetPoint3dAt(j).Y<a.Y))&&a.Y<(p.GetPoint3dAt(i).Y))
    
32.                             {
    
33.                                 if ((a.X - p.GetPoint3dAt(i).X) < (p.GetPoint3dAt(j).X - p.GetPoint3dAt(i).X) * (a.Y - p.GetPoint3dAt(i).Y) / (p.GetPoint3dAt(j).Y - p.GetPoint3dAt(i).Y))
    
34.                                 {
    
35.                                     result = !result;
    
36.                                 }
    
37.                             }
    
38.                             j = i;
    
39.                         }
    
40.                         if (result)
    
41.                         {
    
42.                             MessageBox.Show("在封闭区域中");
    
43.                         }
    
44.                     }

BlockTable bt = trans.GetObject(db.BlockTableId, OpenMode.ForRead) as BlockTable;<br />                 BlockTableRecord btr = trans.GetObject(bt[BlockTableRecord.ModelSpace], OpenMode.ForWrite) as BlockTableRecord;<br />                 Polyline p = new Polyline();<br />                 Point3d a = new Point3d();<br />                 PromptPointOptions opt = new PromptPointOptions("选择点");<br /> PromptPointResult rest = ed.GetPoint(opt);<br /> if (rest.Status == PromptStatus.OK)<br />                 {<br />                     a = rest.Value;<br />                 }<br /> TypedValue[] pline ={ new TypedValue((int)DxfCode.Start, "LWPOLYLINE") };<br />                 SelectionFilter sf = new SelectionFilter(pline);<br />                 PromptSelectionResult res = ed.SelectAll(sf);<br />                 SelectionSet st = res.Value;<br />                 if (st != null)<br />                 {<br />                     ObjectId[] ids = st.GetObjectIds();<br />                     Point3dCollection p3dc=new Point3dCollection();<br />                     foreach (ObjectId id in ids)<br />                     {<br />                         p = trans.GetObject(id, OpenMode.ForRead) as Polyline;<br />                         //Point3d c = new Point3d((p.c));<br />                         int j=p.NumberOfVertices-1;<br />                         bool result = false;<br />                         for (int i = 0; i < p.NumberOfVertices; i++)<br />                         {<br />                             p3dc.Add(p.GetPoint3dAt(i));<br />                             //nodes.Add(p.GetPoint3dAt(i));<br />                             ed.WriteMessage(p.GetPoint3dAt(i).ToString());<br />                             <br />                             if (((p.GetPoint3dAt(i).Y)<=a.Y)&&(a.Y<p.GetPoint3dAt(j).Y)||((p.GetPoint3dAt(j).Y<a.Y))&&a.Y<(p.GetPoint3dAt(i).Y))<br />                             {<br />                                 if ((a.X - p.GetPoint3dAt(i).X) < (p.GetPoint3dAt(j).X - p.GetPoint3dAt(i).X) * (a.Y - p.GetPoint3dAt(i).Y) / (p.GetPoint3dAt(j).Y - p.GetPoint3dAt(i).Y))<br />                                 {<br />                                     result = !result;<br />                                 }<br />                             }<br />                             j = i;<br />                         }<br />                         if (result)<br />                         {<br />                             MessageBox.Show("在封闭区域中");<br />                         }<br />                     }

这个是测试时候的代码。

kuaigoumanzhu

这种方法需要比较的太多，还有那个算点是否在封闭区域内的算法，也不明白。有什么更好的没有大家多多指点啊。

雪山飞狐_lzh

http://bbs.mjtd.com/thread-84681-1-1.html

dafeilang

(bpoly (getpoint "ss"))  找输入点周围的封闭图形并创建多段线

然后可以求创建的面域的面积 （也是在网上学的，这个命令）

不好意思，搜求坐标点的方法的时候进错版了，我以为是Lisp版

yanglin112

新手学习！