2014-04-11 8 views
5

Ho un segmento (x1, y1, x2, y2) che interseca un cerchio di raggio r. Come posso determinare quale punto di intersezione è più vicino a (x1, y1)?Segmento segmento e intersezione cerchio

circle-line

+1

Puoi spiegare cosa intendi per "il punto di intersezione più vicino" È un punto sul cerchio? Un punto sul segmento? Qualcosa non su entrambi? – andand

+4

Questa è una domanda di matematica di base e tradurre la soluzione in C# è davvero semplice e non rende la domanda più pertinente. –

+0

@e il punto di intersezione appartiene a entrambe le linee. – Dmitry

risposta

7

Per fare che prima trovare i punti di intersezione con il cerchio e poi prendere la più vicina al punto di partenza della linea

in modo da controllare questo codice

// cx, cy è il centro punto del cerchio

 public PointF ClosestIntersection(float cx, float cy, float radius, 
      PointF lineStart, PointF lineEnd) 
     { 
      PointF intersection1; 
      PointF intersection2; 
     int intersections = FindLineCircleIntersections(cx, cy, radius, lineStart, lineEnd, out intersection1, out intersection2); 

     if (intersections == 1) 
      return intersection1;//one intersection 

     if (intersections == 2) 
     { 
      double dist1 = Distance(intersection1, lineStart); 
      double dist2 = Distance(intersection2, lineStart); 

      if (dist1 < dist2) 
       return intersection1; 
      else 
       return intersection2; 
     } 

     return PointF.Empty;// no intersections at all 
    } 

    private double Distance(PointF p1, PointF p2) 
    { 
     return Math.Sqrt(Math.Pow(p2.X - p1.X, 2) + Math.Pow(p2.Y - p1.Y, 2)); 
    } 

    // Find the points of intersection. 
    private int FindLineCircleIntersections(float cx, float cy, float radius, 
     PointF point1, PointF point2, out PointF intersection1, out PointF intersection2) 
    { 
     float dx, dy, A, B, C, det, t; 

     dx = point2.X - point1.X; 
     dy = point2.Y - point1.Y; 

     A = dx * dx + dy * dy; 
     B = 2 * (dx * (point1.X - cx) + dy * (point1.Y - cy)); 
     C = (point1.X - cx) * (point1.X - cx) + (point1.Y - cy) * (point1.Y - cy) - radius * radius; 

     det = B * B - 4 * A * C; 
     if ((A <= 0.0000001) || (det < 0)) 
     { 
      // No real solutions. 
      intersection1 = new PointF(float.NaN, float.NaN); 
      intersection2 = new PointF(float.NaN, float.NaN); 
      return 0; 
     } 
     else if (det == 0) 
     { 
      // One solution. 
      t = -B/(2 * A); 
      intersection1 = new PointF(point1.X + t * dx, point1.Y + t * dy); 
      intersection2 = new PointF(float.NaN, float.NaN); 
      return 1; 
     } 
     else 
     { 
      // Two solutions. 
      t = (float)((-B + Math.Sqrt(det))/(2 * A)); 
      intersection1 = new PointF(point1.X + t * dx, point1.Y + t * dy); 
      t = (float)((-B - Math.Sqrt(det))/(2 * A)); 
      intersection2 = new PointF(point1.X + t * dx, point1.Y + t * dy); 
      return 2; 
     } 
    } 

Codice Intersezione modulo qui LINK

0.123.
+1

L'ordine delle distanze non cambia di sqrt, quindi è meglio lasciarle al quadrato per una migliore velocità. – Preza8