练习查找自相交Pline

csharp · 发表于 2014-6-27 11:48:35

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？立即注册

×

[CommandMethod("Test")]
public static void IsPolylineSeftInters()
{
Document doc = Application.DocumentManager.MdiActiveDocument;
Editor ed = doc.Editor;
Database db = doc.Database;
Transaction tr = db.TransactionManager.StartTransaction();
TypedValue typevalue = new TypedValue((int) DxfCode.Start, "LWPOLYLINE");
SelectionFilter filter = new SelectionFilter(new TypedValue[1]{typevalue });
PromptSelectionResult psr = ed.SelectAll( filter );
if (psr.Status != PromptStatus.OK) return;
using (tr)
{
try
{
SelectionSet ss = psr.Value;
ObjectId[] ids = ss.GetObjectIds();
List<Entity> entLst = new List<Entity>();
for (int i = 0; i < ids.Length; i++)
{
Entity ent = (Entity)ids[i].GetObject(OpenMode.ForRead);
if (IsSelfIntersect(ent))
{
entLst.Add(ent);
}
}
if (entLst.Count > 0)
{
for (int i = 0; i < entLst.Count; i++)
{
entLst[i].UpgradeOpen();
entLst[i].ColorIndex = 1;
entLst[i].Highlight();
}
}
tr.Commit();
}
catch (Exception)
{
throw;
}
}
}
/// <summary>
/// 判断Polyline是否自相交
/// </summary>
/// <param name="ent">Entity(Pline)</param>
/// <returns>true or false</returns>
private static bool IsSelfIntersect(Entity ent)
{
Point3dCollection pts = new Point3dCollection();
ent.IntersectWith(ent, Intersect.OnBothOperands, pts, IntPtr.Zero, IntPtr.Zero);
Curve pl = ent as Curve;
Point3d ep = pl.EndPoint;
double[] pams = new double[pts.Count];
for (int i = 0; i < pams.Length ; i++)
{
pams[i] = pl.GetParameterAtPoint(pts[i]);
}
return (IsGrade(pams) == 1 | EndPointIsOn(ent, ep) | NumberInters( ent,pams)) ? true : false;
}
/// <summary>
/// 交点数是否等于总顶点数 - 2
/// </summary>
/// <param name="ent"></param>
/// <param name="pams"></param>
/// <returns></returns>
private static bool NumberInters(Entity ent, double[] pams)
{
Curve pl = ent as Curve;
return Convert .ToInt32( pl.EndParam - 1) != pams.Length;
}
/// <summary>
/// 判断交点数是否等于顶点
/// </summary>
/// <param name="pams"></param>
/// <returns></returns>
private static int IsGrade(double[] pams)
{
var numPlace = pams.SkipWhile( (num, index) => num -1 == index);
return numPlace.Count() == 0 ? 0:1;
}
/// <summary>
/// 判断端点是否和其它顶点重合
/// </summary>
/// <param name="ent">Entity(多义线)</param>
/// <param name="p">端点</param>
/// <returns>重合 true</returns>
private static bool EndPointIsOn(Entity ent, Point3d p)
{
Curve pl = ent as Curve;
int n = Convert.ToInt32(pl.EndParam);
bool ret = false ;
for (int i = 0; i < n-1; i++)
{
if (pl.GetPointAtParameter(Convert .ToDouble( i)).Equals(p))
{
ret = true;
break;
}
}
return ret;
}

[CommandMethod("Test")] 
        public static void IsPolylineSeftInters() 
        { 
            Document doc = Application.DocumentManager.MdiActiveDocument; 
            Editor ed = doc.Editor; 
            Database db = doc.Database; 
            Transaction tr = db.TransactionManager.StartTransaction(); 
 
            TypedValue typevalue = new TypedValue((int) DxfCode.Start, "LWPOLYLINE"); 
            SelectionFilter filter = new SelectionFilter(new TypedValue[1]{typevalue }); 
 
            PromptSelectionResult psr = ed.SelectAll( filter ); 
 
            if (psr.Status != PromptStatus.OK) return; 
            using (tr) 
            { 
                try 
                { 
                    SelectionSet ss = psr.Value; 
                    ObjectId[] ids = ss.GetObjectIds(); 
                    List<Entity> entLst = new List<Entity>(); 
                    for (int i = 0; i < ids.Length; i++) 
                    { 
                        Entity ent = (Entity)ids[i].GetObject(OpenMode.ForRead); 
                        if (IsSelfIntersect(ent)) 
                        { 
                            entLst.Add(ent); 
                        } 
                    } 
                    if (entLst.Count > 0) 
                    { 
                        for (int i = 0; i < entLst.Count; i++) 
                        { 
                            entLst[i].UpgradeOpen(); 
                            entLst[i].ColorIndex = 1; 
                            entLst[i].Highlight(); 
                        } 
                    } 
                    tr.Commit(); 
                } 
                catch (Exception) 
                { 
                     
                    throw; 
                } 
            } 
        } 
        /// <summary> 
        /// 判断Polyline是否自相交 
        /// </summary> 
        /// <param name="ent">Entity(Pline)</param> 
        /// <returns>true or false</returns> 
        private static bool IsSelfIntersect(Entity ent) 
        { 
            Point3dCollection pts = new Point3dCollection(); 
            ent.IntersectWith(ent, Intersect.OnBothOperands, pts, IntPtr.Zero, IntPtr.Zero); 
            Curve pl = ent as Curve; 
            Point3d ep = pl.EndPoint; 
            double[] pams = new double[pts.Count]; 
            for (int i = 0; i < pams.Length ; i++) 
            { 
                pams[i] = pl.GetParameterAtPoint(pts[i]); 
            } 
            return (IsGrade(pams) == 1 | EndPointIsOn(ent, ep) | NumberInters( ent,pams)) ? true : false; 
        } 
        /// <summary> 
        /// 交点数是否等于总顶点数 - 2 
        /// </summary> 
        /// <param name="ent"></param> 
        /// <param name="pams"></param> 
        /// <returns></returns> 
        private static bool NumberInters(Entity ent, double[] pams) 
        { 
            Curve pl = ent as Curve; 
            return Convert .ToInt32( pl.EndParam - 1) != pams.Length; 
        } 
        /// <summary> 
        /// 判断交点数是否等于顶点 
        /// </summary> 
        /// <param name="pams"></param> 
        /// <returns></returns> 
        private static int IsGrade(double[] pams) 
        { 
            var numPlace = pams.SkipWhile( (num, index) => num -1 == index); 
            return numPlace.Count() == 0 ? 0:1; 
        } 
        /// <summary> 
        /// 判断端点是否和其它顶点重合 
        /// </summary> 
        /// <param name="ent">Entity(多义线)</param> 
        /// <param name="p">端点</param> 
        /// <returns>重合 true</returns> 
        private static bool EndPointIsOn(Entity ent, Point3d p) 
        { 
            Curve pl = ent as Curve; 
            int n = Convert.ToInt32(pl.EndParam); 
            bool ret = false ; 
            for (int i = 0; i < n-1; i++) 
            { 
                if (pl.GetPointAtParameter(Convert .ToDouble( i)).Equals(p)) 
                { 
                    ret = true; 
                    break; 
                } 
            } 
            return ret; 
        }

csharp · 发表于 2014-6-27 12:51:12

本帖最后由 csharp 于 2014-6-27 13:02 编辑

三个条件
1 交点数量等于顶点数 - 2
2 交点的参数等于索引值
3 Endpoint单独判断不和任意顶点重合

这里面 using try … Highlight tr.Commit(); 无法完成亮显

csharp · 发表于 2014-6-27 15:21:04

知道了自交点的判别和交点参数就可以写在自交点处断开了

sryf · 发表于 2014-6-27 15:44:24

路过，看看。。。。。。。。。。。。。。。。。。。。。。。。

muyu301 · 发表于 2014-6-27 21:33:02

多谢分享，学习了:)

PLgis · 发表于 2014-6-29 08:27:37

这里面 using try … Highlight tr.Commit(); 无法完成亮显

把entLst.Highlight();
改为：
Highlight(FullSubentityPath subId, bool highlightAll);
是否可以实现亮显？

csharp · 发表于 2014-7-3 09:22:40

判断端点和其它顶点重合循环应该是 i=1, i < n -2, i++

csharp · 发表于 2014-7-10 16:47:11

本帖最后由 csharp 于 2014-7-10 18:51 编辑

用 CurveCurveIntersector3d 计算更简单快捷，GetIntersectionParameters 返回的不是 Pamram ，是距离起点的距离值(GetDistAtPoint)

[CommandMethod("MyGroup", "MyCommand", "MyCommandLocal", CommandFlags.Modal)]
public void MyCommand() // This method can have any name
{
Document doc = Application.DocumentManager.MdiActiveDocument;
Editor ed = doc.Editor;
Database db = doc.Database;
Transaction tr = db.TransactionManager.StartTransaction();
PromptEntityResult per = ed.GetEntity("\nPick Polyline");
if (per.Status != PromptStatus.OK) return;
ObjectId plid = per.ObjectId;
using (tr)
{
try
{
Curve cv = (Curve) plid.GetObject(OpenMode.ForRead);
Plane plan = cv.GetPlane();
Vector3d vv = plan.Normal;
Curve3d c3d = cv.GetGeCurve();
CurveCurveIntersector3d ccid = new CurveCurveIntersector3d(c3d,c3d,vv);
for (int i = 0; i < ccid.NumberOfIntersectionPoints; i++)
{
double[] pams = ccid.GetIntersectionParameters(i);
for (int j = 0; j < pams.Length ; j++)
{
ed.WriteMessage("\n" + pams[j].ToString());
}
}
}
catch (Exception)
{
throw;
}
}
}

[CommandMethod("MyGroup", "MyCommand", "MyCommandLocal", CommandFlags.Modal)] 
        public void MyCommand() // This method can have any name 
        { 
            Document doc = Application.DocumentManager.MdiActiveDocument; 
            Editor ed = doc.Editor; 
            Database db = doc.Database; 
            Transaction tr = db.TransactionManager.StartTransaction(); 
            PromptEntityResult per = ed.GetEntity("\nPick Polyline"); 
            if (per.Status != PromptStatus.OK) return; 
            ObjectId plid = per.ObjectId; 
            using (tr) 
            { 
                try 
                { 
                    Curve cv = (Curve) plid.GetObject(OpenMode.ForRead); 
                    Plane plan = cv.GetPlane(); 
                    Vector3d vv = plan.Normal; 
                    Curve3d c3d = cv.GetGeCurve(); 
                    CurveCurveIntersector3d ccid = new CurveCurveIntersector3d(c3d,c3d,vv); 
                    for (int i = 0; i < ccid.NumberOfIntersectionPoints; i++) 
                    { 
                        double[] pams = ccid.GetIntersectionParameters(i); 
                        for (int j = 0; j < pams.Length ; j++) 
                        { 
                            ed.WriteMessage("\n" + pams[j].ToString()); 
                        } 
                         
                    } 
                } 
                catch (Exception) 
                { 
                     
                    throw; 
                } 
            } 
        }

没有自交点时 ccid.NumberOfIntersectionPoints 为 0 ！

一楼用AcDb 绕弯了，这个才是正道，有了参数在用GetSplitCurves就简单

账号		自动登录	找回密码
密码			立即注册

[原创] 练习 查找自相交Pline

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

[原创] 练习查找自相交Pline