程序喷淋覆盖率

csharp · 发表于 2014-8-28 23:21:18

雪山飞狐(lzh) 发表于 2014-8-28 22:33
如果画一条开口spl 再强行将起点终点合在一起我这里测试会出问题改成这样正常

Join的Spline还没有试过

csharp · 发表于 2014-8-29 07:47:18

本帖最后由 csharp 于 2014-8-29 07:48 编辑

雪山飞狐(lzh) 发表于 2014-8-28 22:33
如果画一条开口spl 再强行将起点终点合在一起我这里测试会出问题改成这样正常

经测试 2014 上 Circle Ellipse 都可以了，不过有个问题是不是没有优化，对交点仅是起点、终点的返回 Clone，这样重新 ToCurve 就没有必要了

闭合 Spline 是 Curve3d 的 SplitCurves 所限只能期待 Autodesk 给出更好方法

csharp · 发表于 2014-8-29 10:32:56

在用 AcDbArc 进行 GetSplitCurves 时，有几个问题会导致程序异常

1 求得交点后再用 GetParameterAtPoint(pt) 求得的如果是 EndParam，有可能比源 Arc.EndParam 大一点点

2 如果有重复交点时就可能产生重复参数，这是 SplitCurves 就会异常，可能 CAD 允许有 0 长度 Line Polyline Spline ，不允许有 0 长度 Arc，Ellipse可能也是如此

雪山飞狐(lzh) · 发表于 2014-8-29 11:25:39

重复参数我的代码应该消除了交点在起点终点的问题我记得代码有去除的

csharp · 发表于 2014-8-29 11:36:06

本帖最后由 csharp 于 2014-8-29 11:37 编辑

对 Arc 去掉起点和终点，如果还剩一个点，这个点参数与 StartParam 或 EndParam 的差值小于误差时也会出错

雪山飞狐(lzh) · 发表于 2014-8-29 11:53:31

这个代码里面有去除吧? 不过圆可能有这个问题

雪山飞狐_lzh · 发表于 2014-8-29 16:35:50

本帖最后由雪山飞狐(lzh) 于 2014-8-29 19:02 编辑

应该是最终版了吧，，，基本修正了Bug 封闭曲线起点断开的问题也解决了

public static bool IsEqual(this Tolerance tol,double d1, double d2)
{
return Math.Abs(d1 - d2) < tol.EqualPoint;
}
public static List<CompositeCurve3d> GetSplitCurves(this CompositeCurve3d c3d, List<double> pars)
{
Interval inter = c3d.GetInterval();
Curve3d[] c3ds = c3d.GetCurves();
pars.Sort();
for (int i = pars.Count - 1; i > 0; i--)
{
if (Tolerance.Global.IsEqual(pars[i],pars[i - 1]))
pars.RemoveAt(i);
}
if (pars.Count == 0)
return new List<CompositeCurve3d>();
if (c3ds.Length == 1 && c3ds[0].IsClosed())
{
//闭合曲线不允许打断于一点
if (pars.Count > 1)
{
//如果包含起点
if (Tolerance.Global.IsEqual(pars[0], inter.LowerBound))
{
pars[0] = inter.LowerBound;
//又包含终点，去除终点
if (Tolerance.Global.IsEqual(pars[pars.Count - 1], inter.UpperBound))
{
pars.RemoveAt(pars.Count - 1);
if (pars.Count == 1)
return new List<CompositeCurve3d>();
}
}
else if (Tolerance.Global.IsEqual(pars[pars.Count - 1], inter.UpperBound))
{
pars[pars.Count - 1] = inter.UpperBound;
}
//加入第一点以支持反向打断
pars.Add(pars[0]);
}
else
{
return new List<CompositeCurve3d>();
}
}
else
{
//非闭合曲线加入起点和终点
if (Tolerance.Global.IsEqual(pars[0], inter.LowerBound))
pars[0] = inter.LowerBound;
else
pars.Insert(0, inter.LowerBound);
if (Tolerance.Global.IsEqual(pars[pars.Count - 1], inter.UpperBound))
pars[pars.Count - 1] = inter.UpperBound;
else
pars.Add(inter.UpperBound);
}
List<CompositeCurve3d> curves = new List<CompositeCurve3d>();
for (int i = 0; i < pars.Count - 1; i++)
{
List<Curve3d> cc3ds = new List<Curve3d>();
//复合曲线参数转换到包含曲线参数
CompositeParameter cp1 = c3d.GlobalToLocalParameter(pars[i]);
CompositeParameter cp2 = c3d.GlobalToLocalParameter(pars[i + 1]);
if (cp1.SegmentIndex == cp2.SegmentIndex)
{
cc3ds.Add(
c3ds[cp1.SegmentIndex].GetSubCurve(
cp1.LocalParameter,
cp2.LocalParameter));
}
else
{
inter = c3ds[cp1.SegmentIndex].GetInterval();
cc3ds.Add(
c3ds[cp1.SegmentIndex].GetSubCurve(
cp1.LocalParameter,
inter.UpperBound));
for (int j = cp1.SegmentIndex + 1; j < cp2.SegmentIndex; j++)
{
cc3ds.Add((Curve3d)c3ds[j].Clone());
}
inter = c3ds[cp2.SegmentIndex].GetInterval();
cc3ds.Add(
c3ds[cp2.SegmentIndex].GetSubCurve(
inter.LowerBound,
cp2.LocalParameter));
}
curves.Add(new CompositeCurve3d(cc3ds.ToArray()));
}
if (c3d.IsClosed() && c3ds.Length > 1)
{
var cs = curves[curves.Count - 1].GetCurves().ToList();
cs.AddRange(curves[0].GetCurves());
curves[curves.Count - 1] = new CompositeCurve3d(cs.ToArray());
curves.RemoveAt(0);
}
return curves;
}

public static bool IsEqual(this Tolerance tol,double d1, double d2) 
        { 
            return Math.Abs(d1 - d2) < tol.EqualPoint; 
        } 
 
        public static List<CompositeCurve3d> GetSplitCurves(this CompositeCurve3d c3d, List<double> pars) 
        { 
            Interval inter = c3d.GetInterval(); 
            Curve3d[] c3ds = c3d.GetCurves(); 
 
            pars.Sort(); 
 
            for (int i = pars.Count - 1; i > 0; i--) 
            { 
                if (Tolerance.Global.IsEqual(pars[i],pars[i - 1])) 
                    pars.RemoveAt(i); 
            } 
 
 
            if (pars.Count == 0) 
                return new List<CompositeCurve3d>(); 
 
            if (c3ds.Length == 1 && c3ds[0].IsClosed()) 
            { 
                //闭合曲线不允许打断于一点 
                if (pars.Count > 1) 
                { 
                    //如果包含起点 
                    if (Tolerance.Global.IsEqual(pars[0], inter.LowerBound)) 
                    { 
                        pars[0] = inter.LowerBound; 
                        //又包含终点，去除终点 
                        if (Tolerance.Global.IsEqual(pars[pars.Count - 1], inter.UpperBound)) 
                        { 
                            pars.RemoveAt(pars.Count - 1); 
                            if (pars.Count == 1) 
                                return new List<CompositeCurve3d>(); 
                        } 
                    } 
                    else if (Tolerance.Global.IsEqual(pars[pars.Count - 1], inter.UpperBound)) 
                    { 
                        pars[pars.Count - 1] = inter.UpperBound; 
                    } 
 
                    //加入第一点以支持反向打断 
                    pars.Add(pars[0]); 
 
                } 
                else 
                { 
                    return new List<CompositeCurve3d>(); 
                } 
            } 
            else 
            { 
                //非闭合曲线加入起点和终点 
                if (Tolerance.Global.IsEqual(pars[0], inter.LowerBound)) 
                    pars[0] = inter.LowerBound; 
                else 
                    pars.Insert(0, inter.LowerBound); 
 
                if (Tolerance.Global.IsEqual(pars[pars.Count - 1], inter.UpperBound)) 
                    pars[pars.Count - 1] = inter.UpperBound; 
                else 
                    pars.Add(inter.UpperBound); 
            } 
          
          
            List<CompositeCurve3d> curves = new List<CompositeCurve3d>(); 
            for (int i = 0; i < pars.Count - 1; i++) 
            { 
                List<Curve3d> cc3ds = new List<Curve3d>(); 
 
                //复合曲线参数转换到包含曲线参数 
                CompositeParameter cp1 = c3d.GlobalToLocalParameter(pars[i]); 
                CompositeParameter cp2 = c3d.GlobalToLocalParameter(pars[i + 1]); 
                if (cp1.SegmentIndex == cp2.SegmentIndex) 
                { 
                    cc3ds.Add( 
                        c3ds[cp1.SegmentIndex].GetSubCurve( 
                            cp1.LocalParameter, 
                            cp2.LocalParameter)); 
                } 
                else 
                { 
                    inter = c3ds[cp1.SegmentIndex].GetInterval(); 
                    cc3ds.Add( 
                        c3ds[cp1.SegmentIndex].GetSubCurve( 
                            cp1.LocalParameter, 
                            inter.UpperBound)); 
 
                    for (int j = cp1.SegmentIndex + 1; j < cp2.SegmentIndex; j++) 
                    { 
                        cc3ds.Add((Curve3d)c3ds[j].Clone()); 
                    } 
 
                    inter = c3ds[cp2.SegmentIndex].GetInterval(); 
                    cc3ds.Add( 
                        c3ds[cp2.SegmentIndex].GetSubCurve( 
                            inter.LowerBound, 
                            cp2.LocalParameter)); 
 
                } 
                curves.Add(new CompositeCurve3d(cc3ds.ToArray())); 
 
            } 
 
            if (c3d.IsClosed() && c3ds.Length > 1) 
            { 
                var cs = curves[curves.Count - 1].GetCurves().ToList(); 
                cs.AddRange(curves[0].GetCurves()); 
                curves[curves.Count - 1] = new CompositeCurve3d(cs.ToArray()); 
                curves.RemoveAt(0); 
            } 
 
            return curves; 
 
        }

public static Curve3d GetSubCurve(this Curve3d curve, double from, double to)
{
Interval inter = curve.GetInterval();
bool atStart = Tolerance.Global.IsEqual(inter.LowerBound, from);
bool atEnd = Tolerance.Global.IsEqual(inter.UpperBound, to);
if (atStart && atEnd)
return (Curve3d)curve.Clone();
if (curve is NurbCurve3d)
{
if (from < to)
{
NurbCurve3d clone = (NurbCurve3d)curve.Clone();
if (atStart || atEnd)
{
clone.HardTrimByParams(from, to);
return clone;
}
else
{
clone.HardTrimByParams(inter.LowerBound, to);
clone.HardTrimByParams(from, to);
return clone;
}
}
else
{
NurbCurve3d clone1 = (NurbCurve3d)curve.Clone();
clone1.HardTrimByParams(from, inter.UpperBound);
NurbCurve3d clone2 = (NurbCurve3d)curve.Clone();
clone2.HardTrimByParams(inter.LowerBound, to);
clone1.JoinWith(clone2);
return clone1;
}
}
else
{
Curve3d clone = (Curve3d)curve.Clone();
clone.SetInterval(new Interval(from, to, 0));
return clone;
}
}

public static Curve3d GetSubCurve(this Curve3d curve, double from, double to) 
        { 
 
            Interval inter = curve.GetInterval(); 
            bool atStart = Tolerance.Global.IsEqual(inter.LowerBound, from); 
            bool atEnd = Tolerance.Global.IsEqual(inter.UpperBound, to); 
 
            if (atStart && atEnd) 
                return (Curve3d)curve.Clone(); 
 
            if (curve is NurbCurve3d) 
            { 
                
                if (from < to) 
                { 
 
                    NurbCurve3d clone = (NurbCurve3d)curve.Clone(); 
                    if (atStart || atEnd) 
                    { 
                        clone.HardTrimByParams(from, to); 
                        return clone; 
                    } 
                    else 
                    { 
                        clone.HardTrimByParams(inter.LowerBound, to); 
                        clone.HardTrimByParams(from, to); 
                        return clone; 
                    } 
 
                } 
                else 
                { 
                    NurbCurve3d clone1 = (NurbCurve3d)curve.Clone(); 
                    clone1.HardTrimByParams(from, inter.UpperBound); 
                    NurbCurve3d clone2 = (NurbCurve3d)curve.Clone(); 
                    clone2.HardTrimByParams(inter.LowerBound, to); 
                    clone1.JoinWith(clone2); 
                    return clone1; 
                } 
 
            } 
            else 
            { 
                Curve3d clone = (Curve3d)curve.Clone(); 
                clone.SetInterval(new Interval(from, to, 0)); 
                return clone; 
            } 
        }

雪山飞狐_lzh · 发表于 2014-8-29 16:39:15

本帖最后由雪山飞狐(lzh) 于 2014-8-29 19:04 编辑

/// <summary>
/// 转换GeCompositeCurve3d为Db曲线
/// </summary>
/// <param name="curve">GeCompositeCurve3d</param>
/// <returns>Db曲线</returns>
public static Curve ToCurve(this CompositeCurve3d curve)
{
Curve3d[] cs = curve.GetCurves();
if (cs.Length == 0)
{
return null;
}
else if (cs.Length == 1)
{
return ToCurve(cs[0]);
}
else
{
if (cs[0] is NurbCurve3d)
{
NurbCurve3d nc3d = (NurbCurve3d)cs[0];
for (int i = 1; i < cs.Length; i++)
nc3d.JoinWith((NurbCurve3d)cs[i]);
return ToCurve(nc3d);
}
else
{
return ToPolyline(curve);
}
}
}
public static NurbCurve3d ToCurve3d(this Spline spl)
{
NurbCurve3d nc3d;
NurbsData ndata = spl.NurbsData;
KnotCollection knots = new KnotCollection();
foreach (Double knot in ndata.GetKnots())
knots.Add(knot);
if (ndata.Rational)
{
nc3d =
new NurbCurve3d(
ndata.Degree,
knots,
ndata.GetControlPoints(),
ndata.GetWeights(),
ndata.Periodic);
}
else
{
nc3d =
new NurbCurve3d(
ndata.Degree,
knots,
ndata.GetControlPoints(),
ndata.Periodic);
}
if (spl.HasFitData)
{
var fdata = spl.FitData;
var vec = new Vector3d();
if (fdata.TangentsExist && (fdata.StartTangent != vec || fdata.EndTangent != vec))
nc3d.SetFitData(fdata.GetFitPoints(), fdata.StartTangent, fdata.EndTangent);
}
return nc3d;
}

复制代码

csharp · 发表于 2014-8-29 16:45:20

一堆圆，断开，正常，再选再断开，有时能把 CAD 挂掉，原因查找中

我试试新代码

csharp · 发表于 2014-8-29 17:15:48

本帖最后由 csharp 于 2014-8-29 17:19 编辑

VS 2012  Net 4.0 (2014) 这句编译不过去
         else if (cs.Length == 1)
         {
            return ToCurve((cs[0]);
         }

雪山飞狐(lzh) · 发表于 2014-8-29 17:39:33

这个函数你直接用createfromgecurve代替吧估计还有点bug 晚上改

csharp · 发表于 2014-8-29 18:26:39

雪山飞狐(lzh) 发表于 2014-8-29 17:39
这个函数你直接用createfromgecurve代替吧估计还有点bug 晚上改

好，Db 那个断开也改差不多了，调好发上来

雪山飞狐(lzh) · 发表于 2014-8-29 19:05:41

GetSplitCurves的Bug 上面的代码已修正

csharp · 发表于 2014-8-29 20:40:21

雪山飞狐(lzh) 发表于 2014-8-29 19:05
GetSplitCurves的Bug 上面的代码已修正

改进的用 Db 曲线直接操作，目前测试还没有发现问题

/// <summary>
/// 曲线集交点断开
/// </summary>
/// <param name="curves">List<Curve></param>
/// <returns>DBObjectCollection</returns>
private static List<Curve> CurveBreakAtInters(Curve[] cvs)
{
int tol = cvs.Length;
List<Curve> ret = new List<Curve>();//新曲线
List<List<double>> pamss = new List<List<double>>(tol);//每个曲线上交点参数
//List<Extents3d> extss = new List<Extents3d>(tol);//包围盒
for (int i = tol - 1; i > -1; i--)
{
pamss.Add(new List<double>());
//extss.Add(cvs[i].GeometricExtents);
}
for (int cur = tol - 1; cur > -1; cur--)
{
Curve cv1 = cvs[cur];
List<double> cv1ps = pamss[cur];
//Extents3d ext1 = extss[cur];
for (int n = cur - 1; n > -1; n--)
{
Curve cv2 = cvs[n];
//Extents3d ext2 = extss[n];
List<double> cv2ps = pamss[n];
Point3dCollection points = new Point3dCollection();
//if (RectIntersects(ext1, ext2) || RectIntersects(ext2, ext1))//包围盒相交的才求交点
using (points)
{
cv1.IntersectWith(cv2, Intersect.OnBothOperands, points, IntPtr.Zero, IntPtr.Zero);
foreach (Point3d pt in points)
{
cv1ps.Add(cv1.GetParameterAtPoint(pt));
cv2ps.Add(cv2.GetParameterAtPoint(pt));
}
}
}
}
for (int i = 0; i < tol; i++)
{
Curve cv = cvs[i];
List<double> value = pamss[i];
if (value.Count == 0) //没有交点
{
continue;
}
else if (value.Count == 1 && cv.IsPeriodic && cv.IsPersistent)//闭合曲线只有一个交点
{
continue;
}
else
{
value.Sort();
int vl = value.Count - 1;
for (int j = vl; j > 0; j--)
{
if (IsEqual( value[j], value[j - 1]))
{
value.RemoveAt(j);
}
}
bool atStart = IsEqual( cvs[i].StartParam , value[0]);
bool atEnd = IsEqual( cvs[i].EndParam ,value[value.Count - 1]);
if (cvs[i] is Circle && atEnd && atStart)
{
value.RemoveAt(0);
}
else
{
if (atEnd)
{
value.RemoveAt(value.Count - 1);
}
if (atStart)
{
value.RemoveAt(0);
}
}
if (value == null || value.Count == 0)
{
continue;
}
else
{
double[] arrpt = value.ToArray();
DoubleCollection pts = new DoubleCollection(arrpt);
DBObjectCollection objs = cv.GetSplitCurves(pts);
var lst =
objs.Cast<Curve>().Where(c => c.GetDistanceAtParameter(c.EndParam) > 1e-6).ToList();
if (lst.Count > 0)
{
ret.AddRange(lst);
cvs[i].UpgradeOpen();
cvs[i].Erase();//断开的就删除
}
}
}
}
return ret;//仅包含新生成的曲线
}
private static bool IsEqual(double a, double b)
{
return Math.Abs(a - b) < Tolerance.Global.EqualPoint;
}

/// <summary> 
        /// 曲线集交点断开 
        /// </summary> 
        /// <param name="curves">List<Curve></param> 
        /// <returns>DBObjectCollection</returns> 
        private static List<Curve> CurveBreakAtInters(Curve[] cvs) 
        { 
            int tol = cvs.Length; 
            List<Curve> ret = new List<Curve>();//新曲线 
            List<List<double>> pamss = new List<List<double>>(tol);//每个曲线上交点参数 
            //List<Extents3d> extss = new List<Extents3d>(tol);//包围盒 
            for (int i = tol - 1; i > -1; i--) 
            { 
                pamss.Add(new List<double>()); 
                //extss.Add(cvs[i].GeometricExtents); 
            } 
            for (int cur = tol - 1; cur > -1; cur--) 
            { 
                Curve cv1 = cvs[cur]; 
                List<double> cv1ps = pamss[cur]; 
                //Extents3d ext1 = extss[cur]; 
                for (int n = cur - 1; n > -1; n--) 
                { 
                    Curve cv2 = cvs[n]; 
                    //Extents3d ext2 = extss[n]; 
                    List<double> cv2ps = pamss[n]; 
                    Point3dCollection points = new Point3dCollection(); 
                    //if (RectIntersects(ext1, ext2) || RectIntersects(ext2, ext1))//包围盒相交的才求交点 
                    using (points) 
                    { 
                        cv1.IntersectWith(cv2, Intersect.OnBothOperands, points, IntPtr.Zero, IntPtr.Zero); 
                        foreach (Point3d pt in points) 
                        { 
                            cv1ps.Add(cv1.GetParameterAtPoint(pt)); 
                            cv2ps.Add(cv2.GetParameterAtPoint(pt)); 
                        } 
                    } 
                } 
            } 
            for (int i = 0; i < tol; i++) 
            { 
                Curve cv = cvs[i]; 
                List<double> value = pamss[i]; 
                if (value.Count == 0) //没有交点 
                { 
                    continue; 
                } 
                else if (value.Count == 1 && cv.IsPeriodic && cv.IsPersistent)//闭合曲线只有一个交点 
                { 
                    continue; 
                } 
                else 
                { 
                    value.Sort(); 
                    int vl = value.Count - 1; 
                    for (int j = vl; j > 0; j--) 
                    { 
                        if (IsEqual( value[j], value[j - 1])) 
                        { 
                            value.RemoveAt(j); 
                        } 
                    } 
                    bool atStart = IsEqual( cvs[i].StartParam , value[0]); 
                    bool atEnd = IsEqual( cvs[i].EndParam ,value[value.Count - 1]); 
                    if (cvs[i] is Circle && atEnd && atStart) 
                    { 
                        value.RemoveAt(0); 
                    } 
                    else 
                    { 
                        if (atEnd) 
                        { 
                            value.RemoveAt(value.Count - 1); 
                        } 
                        if (atStart) 
                        { 
                            value.RemoveAt(0); 
                        } 
                    } 
                    if (value == null || value.Count == 0) 
                    { 
                        continue; 
                    } 
                    else 
                    { 
                        double[] arrpt = value.ToArray(); 
                        DoubleCollection pts = new DoubleCollection(arrpt); 
                        DBObjectCollection objs = cv.GetSplitCurves(pts); 
                        var lst = 
                            objs.Cast<Curve>().Where(c => c.GetDistanceAtParameter(c.EndParam) > 1e-6).ToList(); 
                        if (lst.Count > 0) 
                        { 
                            ret.AddRange(lst); 
                            cvs[i].UpgradeOpen(); 
                            cvs[i].Erase();//断开的就删除 
                        } 
                    } 
                } 
            } 
            return ret;//仅包含新生成的曲线 
        } 
 
        private static bool IsEqual(double a, double b) 
        { 
            return Math.Abs(a - b) < Tolerance.Global.EqualPoint; 
        }

雪山飞狐(lzh) · 发表于 2014-8-29 21:16:37

你的没有考虑自相交?另按参数断开的部分最好做成函数这应该是个通用的

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

[原创] 程序喷淋覆盖率

点评

点评

点评

[原创] 程序 喷淋覆盖率

点评

点评

点评

[原创] 程序喷淋覆盖率