/// <summary> /// Overrides <see cref="CADability.Curve2D.GeneralCurve2D.Intersect (ICurve2D)"/> /// </summary> /// <param name="IntersectWith"></param> /// <returns></returns> public override GeoPoint2DWithParameter[] Intersect(ICurve2D IntersectWith) { Circle2D c2d = IntersectWith as Circle2D; if (c2d != null) { GeoPoint2D[] isp = Geometry.IntersectCC(center, radius, c2d.center, c2d.radius); GeoPoint2DWithParameter[] res = new GeoPoint2DWithParameter[isp.Length]; for (int i = 0; i < isp.Length; ++i) { res[i].p = isp[i]; res[i].par1 = this.PositionOf(isp[i]); res[i].par2 = c2d.PositionOf(isp[i]); } return(res); } Line2D l2d = IntersectWith as Line2D; if (l2d != null) { GeoPoint2D[] isp = Geometry.IntersectLC(l2d.StartPoint, l2d.EndPoint, center, radius); GeoPoint2DWithParameter[] res = new GeoPoint2DWithParameter[isp.Length]; for (int i = 0; i < isp.Length; ++i) { res[i].p = isp[i]; res[i].par1 = this.PositionOf(isp[i]); res[i].par2 = l2d.PositionOf(isp[i]); } return(res); } if (IntersectWith is Ellipse2D) // git auch für Arc, ist dort jeweils implementiert { GeoPoint2DWithParameter[] res = IntersectWith.Intersect(this); for (int i = 0; i < res.Length; ++i) { // Parameter im Ergebnis vertauschen double tmp = res[i].par1; res[i].par1 = res[i].par2; res[i].par2 = tmp; } return(res); } return(base.Intersect(IntersectWith)); // der allgemeine Fall }
/// <summary> /// Overrides <see cref="CADability.Curve2D.GeneralCurve2D.MinDistance (ICurve2D)"/> /// </summary> /// <param name="Other"></param> /// <returns></returns> public override double MinDistance(ICurve2D Other) { if (Other is Line2D) { return(Other.MinDistance(this)); // bei der Linie ist es definiert } else if (Other is Circle2D) { double minDist = Curves2D.SimpleMinimumDistance(this, Other); // dort werden die Schnitt- Start- Endpunkte miteinander verrechnet // jetzt nur noch gucken, ob sich die beiden Kreise oder Bögen annähern Circle2D c = Other as Circle2D; double cdist = Geometry.Dist(c.center, center); if (cdist < Precision.eps) { return(minDist); // mit gleichen Mittelpunkten funktioniert } // das folgende nicht, aber der Abstand ist schon von SimpleMinimumDistance bestimmt // da dort auch die Fußpunkte verwendung finden if (cdist > c.radius + radius) { // Annäherung zweier nicht schneidenden Kreise GeoVector2D dir = c.center - center; // von hier zum anderen dir.Norm(); GeoPoint2D p1 = center + radius * dir; GeoPoint2D p2 = c.center + c.radius * dir.Opposite(); double pos1 = this.PositionOf(p1); double pos2 = c.PositionOf(p2); if (pos1 >= 0.0 && pos1 <= 1.0 && pos2 >= 0.0 && pos2 <= 1.0) { minDist = Math.Min(minDist, cdist - c.radius - radius); } } else if (radius - c.radius > cdist) { // Kreis c liegt innerhalb von diesem Kreis GeoVector2D dir = c.center - center; // von hier zum anderen dir.Norm(); GeoPoint2D p1 = center + radius * dir; GeoPoint2D p2 = c.center + c.radius * dir; double pos1 = this.PositionOf(p1); double pos2 = c.PositionOf(p2); if (pos1 >= 0.0 && pos1 <= 1.0 && pos2 >= 0.0 && pos2 <= 1.0) { minDist = Math.Min(minDist, Geometry.Dist(p1, p2)); } } else if (c.radius - radius > cdist) { // dieser Kreis liegt innerhalb von c GeoVector2D dir = center - c.center; dir.Norm(); GeoPoint2D p1 = c.center + c.radius * dir; GeoPoint2D p2 = center + radius * dir; double pos1 = c.PositionOf(p1); double pos2 = PositionOf(p2); if (pos1 >= 0.0 && pos1 <= 1.0 && pos2 >= 0.0 && pos2 <= 1.0) { minDist = Math.Min(minDist, Geometry.Dist(p1, p2)); } } return(minDist); } else { return(base.MinDistance(Other)); } }
/// <summary> /// Overrides <see cref="CADability.Curve2D.GeneralCurve2D.Intersect (ICurve2D)"/> /// </summary> /// <param name="IntersectWith"></param> /// <returns></returns> public override GeoPoint2DWithParameter[] Intersect(ICurve2D IntersectWith) { // gesucht sind alle Schnittpunkte, also auch in der Verlängerung! Line2D l2d = IntersectWith as Line2D; if (l2d != null) { GeoPoint2D ip; if (Geometry.IntersectLL(startPoint, endPoint, l2d.StartPoint, l2d.EndPoint, out ip)) { double pos1 = this.PositionOf(ip); double pos2 = l2d.PositionOf(ip); GeoPoint2DWithParameter pwp = new GeoPoint2DWithParameter(); pwp.p = ip; pwp.par1 = pos1; pwp.par2 = pos2; return(new GeoPoint2DWithParameter[] { pwp }); } else { return(new GeoPoint2DWithParameter[0]); } } Circle2D c2d = IntersectWith as Circle2D; if (c2d != null) { GeoPoint2D[] isp = Geometry.IntersectLC(startPoint, endPoint, c2d.Center, c2d.Radius); GeoPoint2DWithParameter[] res = new GeoPoint2DWithParameter[isp.Length]; for (int i = 0; i < isp.Length; ++i) { res[i].p = isp[i]; res[i].par1 = this.PositionOf(isp[i]); res[i].par2 = c2d.PositionOf(isp[i]); } return(res); } Ellipse2D e2d = IntersectWith as Ellipse2D; if (e2d != null) { GeoPoint2D[] isp = Geometry.IntersectEL(e2d.center, e2d.majorAxis.Length, e2d.minorAxis.Length, e2d.majorAxis.Angle, startPoint, endPoint); GeoPoint2DWithParameter[] res = new GeoPoint2DWithParameter[isp.Length]; for (int i = 0; i < isp.Length; ++i) { res[i].p = isp[i]; res[i].par1 = this.PositionOf(isp[i]); res[i].par2 = e2d.PositionOf(isp[i]); } return(res); } Polyline2D p2d = IntersectWith as Polyline2D; if (p2d != null) { GeoPoint2DWithParameter[] res = p2d.Intersect(this); // sorum geht es for (int i = 0; i < res.Length; ++i) { double t = res[i].par1; res[i].par1 = res[i].par2; res[i].par2 = t; } return(res); } Path2D pa2d = IntersectWith as Path2D; if (pa2d != null) { GeoPoint2DWithParameter[] res = pa2d.Intersect(this); // sorum geht es for (int i = 0; i < res.Length; ++i) { double t = res[i].par1; res[i].par1 = res[i].par2; res[i].par2 = t; } return(res); } BSpline2D b2d = IntersectWith as BSpline2D; if (b2d != null) { GeoPoint2DWithParameter[] res = b2d.Intersect(this); // sorum geht es for (int i = 0; i < res.Length; ++i) { double t = res[i].par1; res[i].par1 = res[i].par2; res[i].par2 = t; } return(res); } return(base.Intersect(IntersectWith)); }