Wie berechne ich den Schnittpunkt zwischen zwei Linien

Ich bin versucht zu berechnen, den Schnittpunkt zwischen den Linien für eine Optical-Flow-Algorithmus mit Hilfe einer Hough-Transformation. Allerdings bin ich nicht immer die Punkte, die ich sein sollte, wenn ich mit meinem Algorithmus zur Berechnung der Schnittpunkte.

Speichere ich die Zeilen, auf die eine Instanz einer Klasse, die ich erstellt namens ImageLine. Hier ist der code für mein Schnittpunkt-Methode.

Point ImageLine::intersectionWith(ImageLine other)
{
    float A2 = other.Y2() - other.Y1();
    float B2 = other.X2() - other.X1();
    float C2 = A2*other.X1() + B2*other.Y1();

    float A1 = y2 - y1;
    float B1 = x2 - x1;
    float C1 = A1 * x1 + B1 * y1;

   float det = A1*B2 - A2*B1;
   if (det == 0)
   {
        return Point(-1,-1);
   }
   Point d = Point((B2 * C1 - B1 * C2) / det, -(A1 * C2 - A2 * C1) / det);
   return d;
}

Ist diese Methode richtig, oder habe ich etwas falsch gemacht? Soweit ich das beurteilen kann, sollte es funktionieren, da es für einen einzigen Punkt, den ich hart codiert durch, allerdings habe ich nicht in der Lage gewesen, um eine gute Kreuzung bei der Verwendung von echten Daten.

InformationsquelleAutor AndrewSchade | 2013-05-13

  1. 6

    Angesichts der Mathe-Seite: wenn wir zwei line-Gleichungen:

    y = m1 * x + c1

    y = m2 * x + c2

    Den Schnittpunkt: (X , Y), der zwei Zeilen beschrieben, die durch die folgenden Gleichungen:

    Y = m1 * X + c1

    Y = m2 * X + c2

    Punkt erfüllt sowohl Gleichung, d.h.:

    m1 * X + c1 = m2 * X + c2

    (Y - c1) /m1 = (Y - c2) /m2

    damit die Schnittpunkt-Koordinaten sind:

    intersectionX = (c2 - c1) /(m1 - m2)

    intersectionY = (m1*c1 - c2*m2) /m1-m2 oder intersectionY = m1 * intersectionX + c1

    Hinweis: c1, m1, c2, m2 berechnet werden, indem man 2 beliebige Punkte einer Zeile und setzen Sie in die Zeile Gleichungen.

    Y = m1c2 - c1m2 / m2-m1 = m1c2 - c1-m1 ???

    InformationsquelleAutor sethi

  2. 2

    (det == 0) ist unwahrscheinlich, um wahr zu sein, wenn Sie mit floating-point-Arithmetik, weil es nicht präzise.

    Etwas wie (fabs(det) < epsilon) wird Häufig verwendet, für einige geeignete Wert von epsilon (sagen wir 1e-6).

    Wenn das nicht beheben, zeigen einige aktuelle zahlen, zusammen mit dem erwarteten Ergebnis und dem tatsächlichen Ergebnis.

    Beim Umgang mit echten Daten, fast parallele Linien sind auch eher unwahrscheinlich... das ist aber nicht das eigentliche Problem, was passiert, ist, dass dort, wo echte Linien sollen sich überschneiden ist nicht zurückgekehrt. (Ich visuelle Darstellung der Daten für debugging-Zwecke.
    Also, was wird zurückgegeben? Was sind Ihre Beiträge? Was sind die Zwischenwerte?

    InformationsquelleAutor Useless

  3. 0

    Wenn die Formeln korrekt sind, versuchen Sie deklarieren alle Ihre intermediate-Argumente, die als 'double'. Wobei der Unterschied von fast parallelen Linien führen könnte, die in Ihren Produkten sehr nah an einander, so dass die 'float' kann nicht erhalten genug Präzision.

    InformationsquelleAutor NPAssoc

  4. 0

    Ausführliche Formel, gehen Sie bitte zu auf dieser Seite.

    Aber ich Liebe code so, hier, probiert diesen code (ich bekomme es aus github, also alle Kredit geht an den Autor des Codes):

    ///Calculate intersection of two lines.
///\return true if found, false if not found or error
bool LineLineIntersect(double x1, double y1, //Line 1 start
    double x2, double y2, //Line 1 end
    double x3, double y3, //Line 2 start
    double x4, double y4, //Line 2 end
    double &ixOut, double &iyOut) //Output 
{
    //http://mathworld.wolfram.com/Line-LineIntersection.html

    double detL1 = Det(x1, y1, x2, y2);
    double detL2 = Det(x3, y3, x4, y4);
    double x1mx2 = x1 - x2;
    double x3mx4 = x3 - x4;
    double y1my2 = y1 - y2;
    double y3my4 = y3 - y4;

    double xnom = Det(detL1, x1mx2, detL2, x3mx4);
    double ynom = Det(detL1, y1my2, detL2, y3my4);
    double denom = Det(x1mx2, y1my2, x3mx4, y3my4);
    if(denom == 0.0)//Lines don't seem to cross
    {
        ixOut = NAN;
        iyOut = NAN;
        return false;
    }

    ixOut = xnom / denom;   
    iyOut = ynom / denom;
    if(!isfinite(ixOut) || !isfinite(iyOut)) //Probably a numerical issue
        return false;

    return true; //All OK
}

    InformationsquelleAutor 123iamking

Schreibe einen Kommentar