¿Cómo detecta dónde se cruzan dos segmentos de línea? [cerrado]

518

¿Cómo determino si dos líneas se cruzan o no, y si lo hacen, en qué punto x, y?

ReyNestor
fuente
Puede ser útil pensar en los bordes del rectángulo como líneas separadas en lugar del polígono completo.
Ryan Graham
Nota para el moderador : la discusión sobre si esta publicación está o no en el tema pertenece a Meta Stack Overflow. Otros comentarios sobre esto aquí serán eliminados.
Martijn Pieters

Respuestas:

659

Hay un buen enfoque para este problema que utiliza productos cruzados de vectores. Defina el producto cruzado vectorial bidimensional v  ×  w como v x  w y  -  v y  w x .

Suponga que los dos segmentos de línea van de p a p  +  r y de q a q  +  s . Entonces, cualquier punto en la primera línea es representable como p  +  t  r (para un parámetro escalar  t ) y cualquier punto en la segunda línea como q  +  u  s (para un parámetro escalar  u ).

Dos segmentos de línea que se cruzan

Las dos líneas se cruzan si podemos encontrar t y u tal que:

p + t  r = q + u  s

Fórmulas para el punto de intersección.

Cruza ambos lados con s , obteniendo

( p + t  r ) × s = ( q + u  s ) × s

Y como s  ×  s = 0, esto significa

t  ( r × s ) = ( q - p ) × s

Y por lo tanto, resolviendo para t :

t = ( q - p ) × s / ( r × s )

Del mismo modo, podemos resolver por usted :

( p + t  r ) × r = ( q + u  s ) × r

u  ( s × r ) = ( p - q ) × r

u = ( p - q ) × r / ( s × r )

Para reducir el número de pasos de cálculo, es conveniente reescribir esto de la siguiente manera (recordando que s  ×  r = -  r  ×  s ):

u = ( q - p ) × r / ( r × s )

Ahora hay cuatro casos:

  1. Si r  ×  s  = 0 y ( q  -  p ) ×  r  = 0, entonces las dos líneas son colineales.

    En este caso, exprese los puntos finales del segundo segmento ( q y q  +  s ) en términos de la ecuación del primer segmento de línea ( p + t r ):

    t 0 = ( q - p ) ·  r / ( r  ·  r )

    t 1 = ( q + s - p ) ·  r / ( r  ·  r ) = t 0 + s  ·  r / ( r  ·  r )

    Si el intervalo entre t 0 y t 1 se cruza con el intervalo [0, 1], entonces los segmentos de línea son colineales y se superponen; de lo contrario son colineales y disjuntos.

    Observe que si s y r punto en direcciones opuestas, a continuación, s  ·  r <0 y así el intervalo a ser comprobados es [ t 1 , t 0 ] en lugar de [ t 0 , t 1 ].

  2. Si r  ×  s  = 0 y ( q  -  p ) ×  r  ≠ 0, entonces las dos líneas son paralelas y no se cruzan.

  3. Si r  ×  s  ≠ 0 y 0 ≤  t  ≤ 1 y 0 ≤  u  ≤ 1, los dos segmentos de línea se encuentran en el punto p + t  r = q + u  s .

  4. De lo contrario, los dos segmentos de línea no son paralelos pero no se intersecan.

Crédito: este método es la especialización bidimensional del algoritmo de intersección de líneas 3D del artículo "Intersección de dos líneas en tres espacios" de Ronald Goldman, publicado en Graphics Gems , página 304. En tres dimensiones, el caso habitual es que las líneas son asimétricas (ni paralelas ni intersectadas) en cuyo caso el método proporciona los puntos de aproximación más cercana de las dos líneas.

Gareth Rees
fuente
55
@myrkos: No. El primer segmento de línea se ejecuta "de p a p + r", de modo que cuando se representa en términos paramétricos como "p + tr", el segmento corresponde a 0 ≤ t ≤ 1. De manera similar para el otro segmento.
Gareth Rees
77
Gareth, siento que me falta algo, pero ¿cómo se divide (un vector) por un vector? Sus soluciones para t y u terminan con / (r × s), pero (r × s)es un vector, ¿verdad? Un vector (0, 0, rx * sy - ry * sx). Y el lado izquierdo es similarmente un vector paralelo al eje z. Entonces ... ¿acabo de dividir el componente z por el otro componente z? ¿Es la fórmula para t realmente |(q − p) × s| / |(r × s)|?
LarsH
77
@LarsH: ver el primer párrafo.
Gareth Rees
35
Para aquellos interesados, aquí hay una implementación simple de C #, tomando las coordenadas de inicio y fin de PointF para las líneas, que parece estar funcionando: ideone.com/PnPJgb
Matt
24
Arme una implementación de JavaScript después de @Matt. Hice correcciones para los errores señalados por Tekito.
pgkelley
230

FWIW, la siguiente función (en C) detecta las intersecciones de línea y determina el punto de intersección. Se basa en un algoritmo en "Los trucos de los gurús de programación de juegos de Windows de Andre LeMothe ". No es diferente a algunos de los algoritmos en otras respuestas (por ejemplo, Gareth). LeMothe luego usa la Regla de Cramer (no me pregunte) para resolver las ecuaciones por sí mismas.

Puedo dar fe de que funciona en mi débil clon de asteroides, y parece tratar correctamente los casos límite descritos en otras respuestas por Elemental, Dan y Wodzu. También es probablemente más rápido que el código publicado por KingNestor porque es todo multiplicación y división, ¡sin raíces cuadradas!

Supongo que hay algún potencial para dividir por cero allí, aunque no ha sido un problema en mi caso. Suficientemente fácil de modificar para evitar el choque de todos modos.

// Returns 1 if the lines intersect, otherwise 0. In addition, if the lines 
// intersect the intersection point may be stored in the floats i_x and i_y.
char get_line_intersection(float p0_x, float p0_y, float p1_x, float p1_y, 
    float p2_x, float p2_y, float p3_x, float p3_y, float *i_x, float *i_y)
{
    float s1_x, s1_y, s2_x, s2_y;
    s1_x = p1_x - p0_x;     s1_y = p1_y - p0_y;
    s2_x = p3_x - p2_x;     s2_y = p3_y - p2_y;

    float s, t;
    s = (-s1_y * (p0_x - p2_x) + s1_x * (p0_y - p2_y)) / (-s2_x * s1_y + s1_x * s2_y);
    t = ( s2_x * (p0_y - p2_y) - s2_y * (p0_x - p2_x)) / (-s2_x * s1_y + s1_x * s2_y);

    if (s >= 0 && s <= 1 && t >= 0 && t <= 1)
    {
        // Collision detected
        if (i_x != NULL)
            *i_x = p0_x + (t * s1_x);
        if (i_y != NULL)
            *i_y = p0_y + (t * s1_y);
        return 1;
    }

    return 0; // No collision
}

Por cierto, debo decir que en el libro de LeMothe, aunque aparentemente él tiene el algoritmo correcto, el ejemplo concreto muestra enchufes en los números incorrectos y hace cálculos incorrectos. Por ejemplo:

(4 * (4 - 1) + 12 * (7 - 1)) / (17 * 4 + 12 * 10)

= 844 / 0,88

= 0,44

Eso me confundió por horas . :(

Gavin
fuente
99
función getLineIntersection (p0_x, p0_y, p1_x, p1_y, p2_x, p2_y, p3_x, p3_y) {var s1_x, s1_y, s2_x, s2_y; s1_x = p1_x - p0_x; s1_y = p1_y - p0_y; s2_x = p3_x - p2_x; s2_y = p3_y - p2_y; var s, t; s = (-s1_y * (p0_x - p2_x) + s1_x * (p0_y - p2_y)) / (-s2_x * s1_y + s1_x * s2_y); t = (s2_x * (p0_y - p2_y) - s2_y * (p0_x - p2_x)) / (-s2_x * s1_y + s1_x * s2_y);
cortijon

55
if (s> = 0 && s <= 1 && t> = 0 && t <= 1) {// Colisión detectada var intX = p0_x + (t * s1_x); var intY = p0_y + (t * s1_y); return [intX, intY]; } return null; // Sin colisión}
cortijon
13
buen algoritmo, sin embargo fyi no maneja casos donde el determinante es 0. (el -s2_x * s1_y + s1_x * s2_y arriba). Si es 0 (o cerca de 0) las líneas son paralelas o colineales. Si es colineal, entonces la intersección puede ser otro segmento de línea.
seand
16
Las dos operaciones de división se pueden evitar por velocidad (la división cuesta más que la multiplicación); si las líneas se cruzan, necesita una división, si no se cruzan, necesita cero. Primero se debe calcular el denominador y detenerse temprano si es cero (posiblemente agregando código para detectar colinealidad). Luego, en lugar de calcular sy tdirectamente, probar la relación entre los dos numeradores y el denominador. Solo si se confirma que las líneas se cruzan, realmente necesita calcular el valor de t(pero no s).
Qwertie
18
Hice pruebas de rendimiento en todos los algoritmos publicados aquí, y este es al menos dos veces más rápido que cualquiera de los otros. ¡Gracias por publicar!
lajos
63

El problema se reduce a esta pregunta: ¿se cruzan dos líneas de A a B y de C a D? Luego puedes preguntarlo cuatro veces (entre la línea y cada uno de los cuatro lados del rectángulo).

Aquí está la matemática vectorial para hacerlo. Supongo que la línea de A a B es la línea en cuestión y la línea de C a D es una de las líneas rectangulares. Mi notación es que Axes la "coordenada x de A" y Cyes la "coordenada y de C". Y " *" significa producto de punto, por ejemplo A*B = Ax*Bx + Ay*By.

E = B-A = ( Bx-Ax, By-Ay )
F = D-C = ( Dx-Cx, Dy-Cy ) 
P = ( -Ey, Ex )
h = ( (A-C) * P ) / ( F * P )

Este hnúmero es la clave. Si hestá entre 0y 1, las líneas se cruzan, de lo contrario no lo hacen. Si F*Pes cero, por supuesto, no puede hacer el cálculo, pero en este caso las líneas son paralelas y, por lo tanto, solo se cruzan en los casos obvios.

El punto exacto de intersección es C + F*h.

Más diversión:

Si hes exactamente 0 o 1las líneas se tocan en un punto final. Puede considerar esto como una "intersección" o no como mejor le parezca.

Específicamente, hes cuánto tienes que multiplicar la longitud de la línea para tocar exactamente la otra línea.

Por lo tanto, si h<0, significa que la línea del rectángulo está "detrás" de la línea dada (con "dirección" siendo "de A a B"), y si h>1la línea del rectángulo está "delante" de la línea dada.

Derivación:

A y C son vectores que apuntan al inicio de la línea; E y F son los vectores de los extremos de A y C que forman la línea.

Para cualquiera de las dos líneas no paralelas en el plano, debe haber exactamente un par de escalares gy hpara que esta ecuación sea válida:

A + E*g = C + F*h

¿Por qué? Debido a que dos líneas no paralelas deben cruzarse, lo que significa que puede escalar ambas líneas en cierta cantidad cada una y tocarse entre sí.

Al principio, esto parece una ecuación única con dos incógnitas! Pero no lo es cuando se considera que se trata de una ecuación vectorial 2D, lo que significa que esto es realmente un par de ecuaciones en xy y).

Tenemos que eliminar una de estas variables. Una manera fácil es hacer que el Etérmino sea cero. Para hacer eso, tome el producto de puntos de ambos lados de la ecuación usando un vector que punteará a cero con E. Ese vector que llamé Panteriormente, e hice la transformación obvia de E.

Ahora tienes:

A*P = C*P + F*P*h
(A-C)*P = (F*P)*h
( (A-C)*P ) / (F*P) = h
Jason Cohen
fuente
29
Este algoritmo es bueno. Pero hay un agujero en él como lo señaló Dan @ stackoverflow.com/questions/563198/… & Elemental @ stackoverflow.com/questions/563198/… Sería genial si actualizara su respuesta para referencia futura. Gracias.
Chantz el
2
¿Es este algoritmo numéricamente estable? He intentado un enfoque similar y resultó dar resultados extraños al trabajar en flotadores.
milosz
3
Parece haber otro problema con este algoritmo. Cuando se alimenta los puntos A = {1, 0} B = {2, 0} C = {0, 0} D = {1,0}, aunque los segmentos de línea se tocan claramente en un extremo, F P (y también E Q, en línea con la corrección del usuario a continuación) son ambos 0, lo que hace que la división entre 0 encuentre h y g. Todavía estoy trabajando en la solución para este, pero pensé que valía la pena señalar el problema.
candrews
12
Esta respuesta es simplemente incorrecta. Pruebe A = {0,0}, B = {0,1}, C = {0,2} D = {2,0}
Tim Cooper el
66
A + E*g = C + F*hLas dos líneas se cruzan si y solo si la solución a esa ecuación (suponiendo que no sean paralelas) tiene ambas, gyh entre 0 y 1 (in o exclusivo, dependiendo de si cuenta tocar en un punto final).
Daniel Fischer
46

He tratado de implementar el algoritmo descrito tan elegantemente por Jason arriba; Desafortunadamente, mientras trabajaba a través de las matemáticas en la depuración, encontré muchos casos para los que no funciona.

Por ejemplo, considere los puntos A (10,10) B (20,20) C (10,1) D (1,10) da h = .5 y, sin embargo, queda claro al examinar que estos segmentos no están cerca de cada uno otro.

Graficar esto deja en claro que el criterio 0 <h <1 solo indica que el punto de intercepción estaría en CD si existiera, pero no le dice a uno si ese punto está en AB. Para asegurarse de que haya un punto de cruce, debe hacer el cálculo simétrico para la variable gy el requisito de intercepción es: 0 <g <1 Y 0 <h <1

Elemental
fuente
2
Me he estado sacando el pelo tratando de entender por qué la respuesta aceptada no me funcionaba. ¡Muchas gracias!
Matt Bridges
1
También es notable que las condiciones de contorno funcionen en este caso (es decir, para h = 0 o h = 1 o g = 0 o g = 1 las líneas 'solo' se tocan
Elemental
Para las personas que tienen problemas para visualizar el resultado, hice una implementación de esto en Javascript: jsfiddle.net/ferrybig/eokwL9mp
Ferrybig
45

Aquí hay una mejora en la respuesta de Gavin. La solución de marcp es similar también, pero tampoco pospone la división.

Esto en realidad resulta ser una aplicación práctica de la respuesta de Gareth Rees también, porque el equivalente del producto cruzado en 2D es el producto perp-dot, que es de lo que este código usa tres. Cambiar a 3D y usar el producto cruzado, interpolando tanto syt al final, da como resultado los dos puntos más cercanos entre las líneas en 3D. De todos modos, la solución 2D:

int get_line_intersection(float p0_x, float p0_y, float p1_x, float p1_y, 
    float p2_x, float p2_y, float p3_x, float p3_y, float *i_x, float *i_y)
{
    float s02_x, s02_y, s10_x, s10_y, s32_x, s32_y, s_numer, t_numer, denom, t;
    s10_x = p1_x - p0_x;
    s10_y = p1_y - p0_y;
    s32_x = p3_x - p2_x;
    s32_y = p3_y - p2_y;

    denom = s10_x * s32_y - s32_x * s10_y;
    if (denom == 0)
        return 0; // Collinear
    bool denomPositive = denom > 0;

    s02_x = p0_x - p2_x;
    s02_y = p0_y - p2_y;
    s_numer = s10_x * s02_y - s10_y * s02_x;
    if ((s_numer < 0) == denomPositive)
        return 0; // No collision

    t_numer = s32_x * s02_y - s32_y * s02_x;
    if ((t_numer < 0) == denomPositive)
        return 0; // No collision

    if (((s_numer > denom) == denomPositive) || ((t_numer > denom) == denomPositive))
        return 0; // No collision
    // Collision detected
    t = t_numer / denom;
    if (i_x != NULL)
        *i_x = p0_x + (t * s10_x);
    if (i_y != NULL)
        *i_y = p0_y + (t * s10_y);

    return 1;
}

Básicamente, pospone la división hasta el último momento y mueve la mayoría de las pruebas hasta antes de que se realicen ciertos cálculos, agregando así salidas anticipadas. Finalmente, también evita la división por el caso cero que ocurre cuando las líneas son paralelas.

También es posible que desee considerar el uso de una prueba epsilon en lugar de la comparación contra cero. Las líneas que están extremadamente cercanas al paralelo pueden producir resultados que están ligeramente desviados. Esto no es un error, es una limitación con las matemáticas de coma flotante.

iMalc
fuente
1
Falla si algunos de los puntos tienen un valor de 0 ... eso no debería suceder, ¿verdad?
hfossli
1
He corregido un error introducido al diferir la división. t podría ser positivo cuando el número y el nombre eran negativos.
iMalc
2
No funciona si p0-p1 es vertical y p2-p3 es horizontal y los dos segmentos se cruzan. (se ejecuta el primer regreso)
Fabio Dalla Libera
El caso coolinear tiene dos posibilidades: no superposición y superposición. El primero debe devolver falso el segundo verdadero. En su código esto no está probado. siempre devuelve falso como la mayoría de las respuestas aquí. Es una pena que ninguna solución realmente funcione.
AlexWien
3
¿Puedes aclararme por qué todos estos usan nombres de variables tan vagos como en s32_ylugar de algo que describe cómo es point2YDifference?
Supuhstar el
40

Pregunta C: ¿Cómo detecta si dos segmentos de línea se cruzan o no?

He buscado el mismo tema y no estaba contento con las respuestas. Así que escribí un artículo que explica muy detalladamente cómo verificar si dos segmentos de línea se cruzan con muchas imágenes. Hay un código Java completo (y probado).

Aquí está el artículo, recortado en las partes más importantes:

El algoritmo, que verifica si el segmento de línea a se cruza con el segmento de línea b, se ve así:

Ingrese la descripción de la imagen aquí

¿Qué son los cuadros delimitadores? Aquí hay dos cuadros delimitadores de dos segmentos de línea:

ingrese la descripción de la imagen aquí

Si ambos cuadros delimitadores tienen una intersección, mueve el segmento de línea a para que un punto esté en (0 | 0). Ahora tiene una línea a través del origen definido por a. Ahora mueva el segmento de línea b de la misma manera y verifique si los nuevos puntos del segmento de línea b están en diferentes lados de la línea a. Si este es el caso, verifíquelo al revés. Si este es también el caso, los segmentos de línea se cruzan. Si no, no se cruzan.

Pregunta A: ¿Dónde se cruzan dos segmentos de línea?

Usted sabe que dos segmentos de línea a y b se cruzan. Si no lo sabe, verifíquelo con las herramientas que le di en la "Pregunta C".

Ahora puede pasar por algunos casos y obtener la solución con matemáticas de séptimo grado (vea el código y el ejemplo interactivo ).

Pregunta B: ¿Cómo detecta si dos líneas se cruzan o no?

Digamos que su punto A = (x1, y1), punto B = (x2, y2), C = (x_3, y_3), D = (x_4, y_4). Su primera línea está definida por AB (con A! = B), y su segunda línea por CD (con C! = D).

function doLinesIntersect(AB, CD) {
    if (x1 == x2) {
        return !(x3 == x4 && x1 != x3);
    } else if (x3 == x4) {
        return true;
    } else {
        // Both lines are not parallel to the y-axis
        m1 = (y1-y2)/(x1-x2);
        m2 = (y3-y4)/(x3-x4);
        return m1 != m2;
    }
}

Pregunta D: ¿Dónde se cruzan dos líneas?

Verifique con la pregunta B si se cruzan en absoluto.

Las líneas ayb están definidas por dos puntos para cada línea. Básicamente, puede aplicar la misma lógica que se utilizó en la pregunta A.

Martin Thoma
fuente
15
Para ser claros, la pregunta B en esta respuesta es realmente sobre dos líneas que se cruzan, no segmentos de línea. No me estoy quejando; No es incorrecto. Simplemente no quiero que nadie se engañe.
Phord
1
No hay "pregunta C". Y la pregunta D solo vuelve a la pregunta A.
Konrad Viltersten
21

La respuesta una vez aceptada aquí es incorrecta (desde entonces no se ha aceptado, ¡así que hurra!). No elimina correctamente todas las no intersecciones. Trivialmente puede parecer que funciona pero puede fallar, especialmente en el caso de que 0 y 1 se consideren válidos para h.

Considere el siguiente caso:

Líneas en (4,1) - (5,1) y (0,0) - (0,2)

Estas son líneas perpendiculares que claramente no se superponen.

A = (4,1)
B = (5,1)
C = (0,0)
D = (0,2)
E = (5,1) - (4,1) = (- 1,0)
F = (0,2) - (0,0) = (0, -2)
P = (0,1)
h = ((4,1) - (0,0)) punto (0,1) / ((0 , -2) punto (0,1)) = 0

De acuerdo con la respuesta anterior, estos dos segmentos de línea se encuentran en un punto final (valores de 0 y 1). Ese punto final sería:

(0,0) + (0, -2) * 0 = (0,0)

Entonces, aparentemente los dos segmentos de línea se encuentran en (0,0), que está en la línea CD, pero no en la línea AB. Entonces, ¿qué va mal? La respuesta es que los valores de 0 y 1 no son válidos y solo a veces SUCEDEN para predecir correctamente la intersección del punto final. Cuando la extensión de una línea (pero no la otra) se encuentra con el segmento de línea, el algoritmo predice una intersección de segmentos de línea, pero esto no es correcto. Me imagino que al probar comenzando con AB vs CD y luego también probando con CD vs AB, este problema se eliminaría. Solo si ambos caen entre 0 y 1 inclusive, se puede decir que se cruzan.

Recomiendo usar el método de producto cruzado vectorial si debe predecir puntos finales.

-Dan

Dan
fuente
44
La respuesta "aceptada" puede cambiar, por lo que debe llamarla de otra manera. (De hecho, creo que ha cambiado desde su comentario)
Johannes Hoff
14

Versión de Python de la respuesta de iMalc:

def find_intersection( p0, p1, p2, p3 ) :

    s10_x = p1[0] - p0[0]
    s10_y = p1[1] - p0[1]
    s32_x = p3[0] - p2[0]
    s32_y = p3[1] - p2[1]

    denom = s10_x * s32_y - s32_x * s10_y

    if denom == 0 : return None # collinear

    denom_is_positive = denom > 0

    s02_x = p0[0] - p2[0]
    s02_y = p0[1] - p2[1]

    s_numer = s10_x * s02_y - s10_y * s02_x

    if (s_numer < 0) == denom_is_positive : return None # no collision

    t_numer = s32_x * s02_y - s32_y * s02_x

    if (t_numer < 0) == denom_is_positive : return None # no collision

    if (s_numer > denom) == denom_is_positive or (t_numer > denom) == denom_is_positive : return None # no collision


    # collision detected

    t = t_numer / denom

    intersection_point = [ p0[0] + (t * s10_x), p0[1] + (t * s10_y) ]


    return intersection_point
Kris
fuente
Recuerde que debe hacer que sus números floten o cambiar la línea 8 para usardenom = float(...)
Jonno_FTW
11

Encontrar la intersección correcta de dos segmentos de línea es una tarea no trivial con muchos casos de borde. Aquí hay una solución bien documentada, funcional y probada en Java.

En esencia, hay tres cosas que pueden suceder al encontrar la intersección de dos segmentos de línea:

  1. Los segmentos no se cruzan

  2. Hay un punto de intersección único

  3. La intersección es otro segmento.

NOTA : En el código, supongo que un segmento de línea (x1, y1), (x2, y2) con x1 = x2 e y1 = y2 es un segmento de línea válido. Matemáticamente hablando, un segmento de línea consta de puntos distintos, pero estoy permitiendo que los segmentos sean puntos en esta implementación para completar.

El código está tomado de mi repositorio de github

/**
 * This snippet finds the intersection of two line segments.
 * The intersection may either be empty, a single point or the
 * intersection is a subsegment there's an overlap.
 */

import static java.lang.Math.abs;
import static java.lang.Math.max;
import static java.lang.Math.min;

import java.util.ArrayList;
import java.util.List;

public class LineSegmentLineSegmentIntersection {

  // Small epsilon used for double value comparison.
  private static final double EPS = 1e-5;

  // 2D Point class.
  public static class Pt {
    double x, y;
    public Pt(double x, double y) {
      this.x = x; 
      this.y = y;
    }
    public boolean equals(Pt pt) {
      return abs(x - pt.x) < EPS && abs(y - pt.y) < EPS;
    }
  }

  // Finds the orientation of point 'c' relative to the line segment (a, b)
  // Returns  0 if all three points are collinear.
  // Returns -1 if 'c' is clockwise to segment (a, b), i.e right of line formed by the segment.
  // Returns +1 if 'c' is counter clockwise to segment (a, b), i.e left of line
  // formed by the segment.
  public static int orientation(Pt a, Pt b, Pt c) {
    double value = (b.y - a.y) * (c.x - b.x) - 
                   (b.x - a.x) * (c.y - b.y);
    if (abs(value) < EPS) return 0;
    return (value > 0) ? -1 : +1;
  }

  // Tests whether point 'c' is on the line segment (a, b).
  // Ensure first that point c is collinear to segment (a, b) and
  // then check whether c is within the rectangle formed by (a, b)
  public static boolean pointOnLine(Pt a, Pt b, Pt c) {
    return orientation(a, b, c) == 0 && 
           min(a.x, b.x) <= c.x && c.x <= max(a.x, b.x) && 
           min(a.y, b.y) <= c.y && c.y <= max(a.y, b.y);
  }

  // Determines whether two segments intersect.
  public static boolean segmentsIntersect(Pt p1, Pt p2, Pt p3, Pt p4) {

    // Get the orientation of points p3 and p4 in relation
    // to the line segment (p1, p2)
    int o1 = orientation(p1, p2, p3);
    int o2 = orientation(p1, p2, p4);
    int o3 = orientation(p3, p4, p1);
    int o4 = orientation(p3, p4, p2);

    // If the points p1, p2 are on opposite sides of the infinite
    // line formed by (p3, p4) and conversly p3, p4 are on opposite
    // sides of the infinite line formed by (p1, p2) then there is
    // an intersection.
    if (o1 != o2 && o3 != o4) return true;

    // Collinear special cases (perhaps these if checks can be simplified?)
    if (o1 == 0 && pointOnLine(p1, p2, p3)) return true;
    if (o2 == 0 && pointOnLine(p1, p2, p4)) return true;
    if (o3 == 0 && pointOnLine(p3, p4, p1)) return true;
    if (o4 == 0 && pointOnLine(p3, p4, p2)) return true;

    return false;
  }

  public static List<Pt> getCommonEndpoints(Pt p1, Pt p2, Pt p3, Pt p4) {

    List<Pt> points = new ArrayList<>();

    if (p1.equals(p3)) {
      points.add(p1);
      if (p2.equals(p4)) points.add(p2);

    } else if (p1.equals(p4)) {
      points.add(p1);
      if (p2.equals(p3)) points.add(p2);

    } else if (p2.equals(p3)) {
      points.add(p2);
      if (p1.equals(p4)) points.add(p1);

    } else if (p2.equals(p4)) {
      points.add(p2);
      if (p1.equals(p3)) points.add(p1);
    }

    return points;
  }

  // Finds the intersection point(s) of two line segments. Unlike regular line 
  // segments, segments which are points (x1 = x2 and y1 = y2) are allowed.
  public static Pt[] lineSegmentLineSegmentIntersection(Pt p1, Pt p2, Pt p3, Pt p4) {

    // No intersection.
    if (!segmentsIntersect(p1, p2, p3, p4)) return new Pt[]{};

    // Both segments are a single point.
    if (p1.equals(p2) && p2.equals(p3) && p3.equals(p4))
      return new Pt[]{p1};

    List<Pt> endpoints = getCommonEndpoints(p1, p2, p3, p4);
    int n = endpoints.size();

    // One of the line segments is an intersecting single point.
    // NOTE: checking only n == 1 is insufficient to return early
    // because the solution might be a sub segment.
    boolean singleton = p1.equals(p2) || p3.equals(p4);
    if (n == 1 && singleton) return new Pt[]{endpoints.get(0)};

    // Segments are equal.
    if (n == 2) return new Pt[]{endpoints.get(0), endpoints.get(1)};

    boolean collinearSegments = (orientation(p1, p2, p3) == 0) && 
                                (orientation(p1, p2, p4) == 0);

    // The intersection will be a sub-segment of the two
    // segments since they overlap each other.
    if (collinearSegments) {

      // Segment #2 is enclosed in segment #1
      if (pointOnLine(p1, p2, p3) && pointOnLine(p1, p2, p4))
        return new Pt[]{p3, p4};

      // Segment #1 is enclosed in segment #2
      if (pointOnLine(p3, p4, p1) && pointOnLine(p3, p4, p2))
        return new Pt[]{p1, p2};

      // The subsegment is part of segment #1 and part of segment #2.
      // Find the middle points which correspond to this segment.
      Pt midPoint1 = pointOnLine(p1, p2, p3) ? p3 : p4;
      Pt midPoint2 = pointOnLine(p3, p4, p1) ? p1 : p2;

      // There is actually only one middle point!
      if (midPoint1.equals(midPoint2)) return new Pt[]{midPoint1};

      return new Pt[]{midPoint1, midPoint2};
    }

    /* Beyond this point there is a unique intersection point. */

    // Segment #1 is a vertical line.
    if (abs(p1.x - p2.x) < EPS) {
      double m = (p4.y - p3.y) / (p4.x - p3.x);
      double b = p3.y - m * p3.x;
      return new Pt[]{new Pt(p1.x, m * p1.x + b)};
    }

    // Segment #2 is a vertical line.
    if (abs(p3.x - p4.x) < EPS) {
      double m = (p2.y - p1.y) / (p2.x - p1.x);
      double b = p1.y - m * p1.x;
      return new Pt[]{new Pt(p3.x, m * p3.x + b)};
    }

    double m1 = (p2.y - p1.y) / (p2.x - p1.x);
    double m2 = (p4.y - p3.y) / (p4.x - p3.x);
    double b1 = p1.y - m1 * p1.x;
    double b2 = p3.y - m2 * p3.x;
    double x = (b2 - b1) / (m1 - m2);
    double y = (m1 * b2 - m2 * b1) / (m1 - m2);

    return new Pt[]{new Pt(x, y)};
  }

}

Aquí hay un ejemplo de uso simple:

  public static void main(String[] args) {

    // Segment #1 is (p1, p2), segment #2 is (p3, p4)
    Pt p1, p2, p3, p4;

    p1 = new Pt(-2, 4); p2 = new Pt(3, 3);
    p3 = new Pt(0, 0);  p4 = new Pt(2, 4);
    Pt[] points = lineSegmentLineSegmentIntersection(p1, p2, p3, p4);
    Pt point = points[0];

    // Prints: (1.636, 3.273)
    System.out.printf("(%.3f, %.3f)\n", point.x, point.y);

    p1 = new Pt(-10, 0); p2 = new Pt(+10, 0);
    p3 = new Pt(-5, 0);  p4 = new Pt(+5, 0);
    points = lineSegmentLineSegmentIntersection(p1, p2, p3, p4);
    Pt point1 = points[0], point2 = points[1];

    // Prints: (-5.000, 0.000) (5.000, 0.000)
    System.out.printf("(%.3f, %.3f) (%.3f, %.3f)\n", point1.x, point1.y, point2.x, point2.y);
  }
will.fiset
fuente
¡Funcionó para mi sistema de coordenadas geográficas! ¡Gracias! Pero es para la intersección de líneas infinitas, y estoy más buscando intersección de líneas finitas.
M. Usman Khan
8

Solo quería mencionar que se puede encontrar una buena explicación y una solución explícita en la serie Numeric Recipes. Tengo la tercera edición y la respuesta está en la página 1117, sección 21.4. Otra solución con una nomenclatura diferente se puede encontrar en un documento de Marina Gavrilova Prueba de intersección de sección de línea confiable . Su solución es, en mi opinión, un poco más simple.

Mi implementación está abajo:

bool NuGeometry::IsBetween(const double& x0, const double& x, const double& x1){
   return (x >= x0) && (x <= x1);
}

bool NuGeometry::FindIntersection(const double& x0, const double& y0, 
     const double& x1, const double& y1,
     const double& a0, const double& b0, 
     const double& a1, const double& b1, 
     double& xy, double& ab) {
   // four endpoints are x0, y0 & x1,y1 & a0,b0 & a1,b1
   // returned values xy and ab are the fractional distance along xy and ab
   // and are only defined when the result is true

   bool partial = false;
   double denom = (b0 - b1) * (x0 - x1) - (y0 - y1) * (a0 - a1);
   if (denom == 0) {
      xy = -1;
      ab = -1;
   } else {
      xy = (a0 * (y1 - b1) + a1 * (b0 - y1) + x1 * (b1 - b0)) / denom;
      partial = NuGeometry::IsBetween(0, xy, 1);
      if (partial) {
         // no point calculating this unless xy is between 0 & 1
         ab = (y1 * (x0 - a1) + b1 * (x1 - x0) + y0 * (a1 - x1)) / denom; 
      }
   }
   if ( partial && NuGeometry::IsBetween(0, ab, 1)) {
      ab = 1-ab;
      xy = 1-xy;
      return true;
   }  else return false;
}
marcp
fuente
No funciona para p1 = (0,0), p2 = (10,0), p3 = (9,0), p4 = (20,0)
padmalcom
Depende de su definición de "no funciona", supongo. Denom es 0, por lo que devolverá falso que me parece correcto ya que no se cruzan. Colinear no es lo mismo que intersectar.
marcp
8

Hay muchas soluciones disponibles anteriormente, pero creo que la solución a continuación es bastante simple y fácil de entender.

Dos segmentos Vector AB y Vector CD se cruzan si y solo si

  1. Los puntos finales ayb están en lados opuestos del segmento CD.
  2. Los puntos finales cyd están en el lado opuesto del segmento AB.

Más específicamente, ayb están en el lado opuesto del segmento CD si y solo si exactamente uno de los dos triples a, c, d y b, c, d está en el sentido contrario a las agujas del reloj.

Intersect(a, b, c, d)
 if CCW(a, c, d) == CCW(b, c, d)
    return false;
 else if CCW(a, b, c) == CCW(a, b, d)
    return false;
 else
    return true;

Aquí CCW representa en sentido antihorario que devuelve verdadero / falso según la orientación de los puntos.

Fuente: http://compgeom.cs.uiuc.edu/~jeffe/teaching/373/notes/x06-sweepline.pdf Página 2

zstring
fuente
2
Creo que debería ser un poco más específico: ¿cómo se CCWdefine la prueba? Con el signo del producto exterior?
ocramz
Gracias; este pseudocódigo permitió una implementación muy sencilla en Scratch; vea este proyecto: scratch.mit.edu/projects/129319027
Ruud Helderman el
8

C y Objetivo-C

Basado en la respuesta de Gareth Rees

const AGKLine AGKLineZero = (AGKLine){(CGPoint){0.0, 0.0}, (CGPoint){0.0, 0.0}};

AGKLine AGKLineMake(CGPoint start, CGPoint end)
{
    return (AGKLine){start, end};
}

double AGKLineLength(AGKLine l)
{
    return CGPointLengthBetween_AGK(l.start, l.end);
}

BOOL AGKLineIntersection(AGKLine l1, AGKLine l2, CGPoint *out_pointOfIntersection)
{
    // http://stackoverflow.com/a/565282/202451

    CGPoint p = l1.start;
    CGPoint q = l2.start;
    CGPoint r = CGPointSubtract_AGK(l1.end, l1.start);
    CGPoint s = CGPointSubtract_AGK(l2.end, l2.start);

    double s_r_crossProduct = CGPointCrossProductZComponent_AGK(r, s);
    double t = CGPointCrossProductZComponent_AGK(CGPointSubtract_AGK(q, p), s) / s_r_crossProduct;
    double u = CGPointCrossProductZComponent_AGK(CGPointSubtract_AGK(q, p), r) / s_r_crossProduct;

    if(t < 0 || t > 1.0 || u < 0 || u > 1.0)
    {
        if(out_pointOfIntersection != NULL)
        {
            *out_pointOfIntersection = CGPointZero;
        }
        return NO;
    }
    else
    {
        if(out_pointOfIntersection != NULL)
        {
            CGPoint i = CGPointAdd_AGK(p, CGPointMultiply_AGK(r, t));
            *out_pointOfIntersection = i;
        }
        return YES;
    }
}

CGFloat CGPointCrossProductZComponent_AGK(CGPoint v1, CGPoint v2)
{
    return v1.x * v2.y - v1.y * v2.x;
}

CGPoint CGPointSubtract_AGK(CGPoint p1, CGPoint p2)
{
    return (CGPoint){p1.x - p2.x, p1.y - p2.y};
}

CGPoint CGPointAdd_AGK(CGPoint p1, CGPoint p2)
{
    return (CGPoint){p1.x + p2.x, p1.y + p2.y};
}

CGFloat CGPointCrossProductZComponent_AGK(CGPoint v1, CGPoint v2)
{
    return v1.x * v2.y - v1.y * v2.x;
}

CGPoint CGPointMultiply_AGK(CGPoint p1, CGFloat factor)
{
    return (CGPoint){p1.x * factor, p1.y * factor};
}

Muchas de las funciones y estructuras son privadas, pero es bastante fácil saber qué está pasando. Esto es público en este repositorio https://github.com/hfossli/AGGeometryKit/

hfossli
fuente
¿De dónde viene AGPointZero en este código?
seanicus
1
@seanicus actualizó el ejemplo para usar CGPoint en su lugar
hfossli
6

Esto está funcionando bien para mí. Tomado de aquí .

 // calculates intersection and checks for parallel lines.  
 // also checks that the intersection point is actually on  
 // the line segment p1-p2  
 Point findIntersection(Point p1,Point p2,  
   Point p3,Point p4) {  
   float xD1,yD1,xD2,yD2,xD3,yD3;  
   float dot,deg,len1,len2;  
   float segmentLen1,segmentLen2;  
   float ua,ub,div;  

   // calculate differences  
   xD1=p2.x-p1.x;  
   xD2=p4.x-p3.x;  
   yD1=p2.y-p1.y;  
   yD2=p4.y-p3.y;  
   xD3=p1.x-p3.x;  
   yD3=p1.y-p3.y;    

   // calculate the lengths of the two lines  
   len1=sqrt(xD1*xD1+yD1*yD1);  
   len2=sqrt(xD2*xD2+yD2*yD2);  

   // calculate angle between the two lines.  
   dot=(xD1*xD2+yD1*yD2); // dot product  
   deg=dot/(len1*len2);  

   // if abs(angle)==1 then the lines are parallell,  
   // so no intersection is possible  
   if(abs(deg)==1) return null;  

   // find intersection Pt between two lines  
   Point pt=new Point(0,0);  
   div=yD2*xD1-xD2*yD1;  
   ua=(xD2*yD3-yD2*xD3)/div;  
   ub=(xD1*yD3-yD1*xD3)/div;  
   pt.x=p1.x+ua*xD1;  
   pt.y=p1.y+ua*yD1;  

   // calculate the combined length of the two segments  
   // between Pt-p1 and Pt-p2  
   xD1=pt.x-p1.x;  
   xD2=pt.x-p2.x;  
   yD1=pt.y-p1.y;  
   yD2=pt.y-p2.y;  
   segmentLen1=sqrt(xD1*xD1+yD1*yD1)+sqrt(xD2*xD2+yD2*yD2);  

   // calculate the combined length of the two segments  
   // between Pt-p3 and Pt-p4  
   xD1=pt.x-p3.x;  
   xD2=pt.x-p4.x;  
   yD1=pt.y-p3.y;  
   yD2=pt.y-p4.y;  
   segmentLen2=sqrt(xD1*xD1+yD1*yD1)+sqrt(xD2*xD2+yD2*yD2);  

   // if the lengths of both sets of segments are the same as  
   // the lenghts of the two lines the point is actually  
   // on the line segment.  

   // if the point isn’t on the line, return null  
   if(abs(len1-segmentLen1)>0.01 || abs(len2-segmentLen2)>0.01)  
     return null;  

   // return the valid intersection  
   return pt;  
 }  

 class Point{  
   float x,y;  
   Point(float x, float y){  
     this.x = x;  
     this.y = y;  
   }  

   void set(float x, float y){  
     this.x = x;  
     this.y = y;  
   }  
 }  
ReyNestor
fuente
8
Hay varios problemas con este código. Puede generar una excepción debido a la división por cero; es lento porque tiene raíces cuadradas; y a veces devuelve falsos positivos porque usa un factor fudge. ¡Puedes hacerlo mejor que esto!
Gareth Rees
Está bien como solución, pero la dada por Jason definitivamente es computacionalmente más rápida y evita muchos de los problemas con esta solución
Elemental el
6

Intenté algunas de estas respuestas, pero no me funcionaron (lo siento muchachos); Después de buscar más en la red, encontré esto .

Con una pequeña modificación a su código, ahora tengo esta función que devolverá el punto de intersección o, si no se encuentra ninguna intersección, devolverá -1, -1.

    Public Function intercetion(ByVal ax As Integer, ByVal ay As Integer, ByVal bx As Integer, ByVal by As Integer, ByVal cx As Integer, ByVal cy As Integer, ByVal dx As Integer, ByVal dy As Integer) As Point
    '//  Determines the intersection point of the line segment defined by points A and B
    '//  with the line segment defined by points C and D.
    '//
    '//  Returns YES if the intersection point was found, and stores that point in X,Y.
    '//  Returns NO if there is no determinable intersection point, in which case X,Y will
    '//  be unmodified.

    Dim distAB, theCos, theSin, newX, ABpos As Double

    '//  Fail if either line segment is zero-length.
    If ax = bx And ay = by Or cx = dx And cy = dy Then Return New Point(-1, -1)

    '//  Fail if the segments share an end-point.
    If ax = cx And ay = cy Or bx = cx And by = cy Or ax = dx And ay = dy Or bx = dx And by = dy Then Return New Point(-1, -1)

    '//  (1) Translate the system so that point A is on the origin.
    bx -= ax
    by -= ay
    cx -= ax
    cy -= ay
    dx -= ax
    dy -= ay

    '//  Discover the length of segment A-B.
    distAB = Math.Sqrt(bx * bx + by * by)

    '//  (2) Rotate the system so that point B is on the positive X axis.
    theCos = bx / distAB
    theSin = by / distAB
    newX = cx * theCos + cy * theSin
    cy = cy * theCos - cx * theSin
    cx = newX
    newX = dx * theCos + dy * theSin
    dy = dy * theCos - dx * theSin
    dx = newX

    '//  Fail if segment C-D doesn't cross line A-B.
    If cy < 0 And dy < 0 Or cy >= 0 And dy >= 0 Then Return New Point(-1, -1)

    '//  (3) Discover the position of the intersection point along line A-B.
    ABpos = dx + (cx - dx) * dy / (dy - cy)

    '//  Fail if segment C-D crosses line A-B outside of segment A-B.
    If ABpos < 0 Or ABpos > distAB Then Return New Point(-1, -1)

    '//  (4) Apply the discovered position to line A-B in the original coordinate system.
    '*X=Ax+ABpos*theCos
    '*Y=Ay+ABpos*theSin

    '//  Success.
    Return New Point(ax + ABpos * theCos, ay + ABpos * theSin)
End Function
Robert
fuente
6

Parece haber cierto interés en la respuesta de Gavin para la cual Cortijon propuso una versión de JavaScript en los comentarios e iMalc proporcionó una versión con un poco menos de cálculos . Algunos han señalado las deficiencias con varias propuestas de código y otros han comentado sobre la eficacia de algunas propuestas de código.

El algoritmo proporcionado por iMalc a través de la respuesta de Gavin es el que estoy usando actualmente en un proyecto de JavaScript y solo quería proporcionar una versión limpia aquí si puede ayudar a alguien.

// Some variables for reuse, others may do this differently
var p0x, p1x, p2x, p3x, ix,
    p0y, p1y, p2y, p3y, iy,
    collisionDetected;

// do stuff, call other functions, set endpoints...

// note: for my purpose I use |t| < |d| as opposed to
// |t| <= |d| which is equivalent to 0 <= t < 1 rather than
// 0 <= t <= 1 as in Gavin's answer - results may vary

var lineSegmentIntersection = function(){
    var d, dx1, dx2, dx3, dy1, dy2, dy3, s, t;

    dx1 = p1x - p0x;      dy1 = p1y - p0y;
    dx2 = p3x - p2x;      dy2 = p3y - p2y;
    dx3 = p0x - p2x;      dy3 = p0y - p2y;

    collisionDetected = 0;

    d = dx1 * dy2 - dx2 * dy1;

    if(d !== 0){
        s = dx1 * dy3 - dx3 * dy1;
        if((s <= 0 && d < 0 && s >= d) || (s >= 0 && d > 0 && s <= d)){
            t = dx2 * dy3 - dx3 * dy2;
            if((t <= 0 && d < 0 && t > d) || (t >= 0 && d > 0 && t < d)){
                t = t / d;
                collisionDetected = 1;
                ix = p0x + t * dx1;
                iy = p0y + t * dy1;
            }
        }
    }
};
Nolo
fuente
No entiendo cómo puedes entender lo que está sucediendo con líneas como t = dx2 * dy3 - dx3 * dy2;...
Supuhstar
@Supuhstar Tiene que ver con las matemáticas vectoriales y la definición de producto punto y producto cruzado. Por ejemplo, el código que publicó representa una operación de producto cruzado. Es una forma de proyectar un segmento de línea sobre otro para determinar dónde cae en el otro segmento de línea, antes del punto de partida en algún punto intermedio o después de la línea. Entonces t es un valor normalizado. Si está entre 0 y 1, entonces los dos segmentos se cruzan. Si es menor que 0 o mayor que uno, entonces no lo hacen.
Nolo
@Supuhstar También tenga en cuenta que para que la proyección encuentre el punto real, el resultado debe ser escalado. Ahí es donde t/dentra.
Nolo
1
Quiero decir, ¿cómo entiendes lo que está sucediendo de un vistazo con nombres de variables como ese? ¿Por qué no algo así crossProduct = (line1XDifference * line2YDifference) - (line2XDifference * line1YDifference)y scaledResult = crossProduct / dotProduct?
Supuhstar el
1
@Supuhstar Ah, entiendo lo que quieres decir. Erm, bueno, supongo que no hay una buena razón para hablar más allá de obsesionarse con la eficiencia, pero esa no es una muy buena razón en sí misma porque los compiladores hacen un muy buen trabajo al tomar la mayoría de los códigos que les das y hacerlo tan eficiente como sea posible. posible sin cambiar lo que se debe calcular. Por otro lado, los nombres, p1x, p1yetc. están destinados a describir puntos por sus valores x e y, por lo que p1xes una abreviatura de point1x, del mismo modo d1x, en mi mente es una abreviatura de la letra griega deltaXo podría decirse differenceInX. (más)
Nolo
5

Creo que hay una solución mucho más simple para este problema. Se me ocurrió otra idea hoy y parece funcionar bien (al menos en 2D por ahora). Todo lo que tiene que hacer es calcular la intersección entre dos líneas, luego verificar si el punto de intersección calculado está dentro de los cuadros de límite de ambos segmentos de línea. Si es así, los segmentos de línea se cruzan. Eso es.

EDITAR:

Así es como calculo la intersección (ya no sé dónde encontré este fragmento de código)

Point3D

viene de

System.Windows.Media.Media3D

public static Point3D? Intersection(Point3D start1, Point3D end1, Point3D start2, Point3D end2) {

        double a1 = end1.Y - start1.Y;
        double b1 = start1.X - end1.X;
        double c1 = a1 * start1.X + b1 * start1.Y;

        double a2 = end2.Y - start2.Y;
        double b2 = start2.X - end2.X;
        double c2 = a2 * start2.X + b2 * start2.Y;

        double det = a1 * b2 - a2 * b1;
        if (det == 0) { // lines are parallel
            return null;
        }

        double x = (b2 * c1 - b1 * c2) / det;
        double y = (a1 * c2 - a2 * c1) / det;

        return new Point3D(x, y, 0.0);
    }

y esta es mi clase (simplificada para el propósito de la respuesta) BoundingBox:

public class BoundingBox {
    private Point3D min = new Point3D();
    private Point3D max = new Point3D();

    public BoundingBox(Point3D point) {
        min = point;
        max = point;
    }

    public Point3D Min {
        get { return min; }
        set { min = value; }
    }

    public Point3D Max {
        get { return max; }
        set { max = value; }
    }

    public bool Contains(BoundingBox box) {
        bool contains =
            min.X <= box.min.X && max.X >= box.max.X &&
            min.Y <= box.min.Y && max.Y >= box.max.Y &&
            min.Z <= box.min.Z && max.Z >= box.max.Z;
        return contains;
    }

    public bool Contains(Point3D point) {
        return Contains(new BoundingBox(point));
    }

}
t3chb0t
fuente
3

Esta solución puede ayudar

public static float GetLineYIntesept(PointF p, float slope)
    {
        return p.Y - slope * p.X;
    }

    public static PointF FindIntersection(PointF line1Start, PointF line1End, PointF line2Start, PointF line2End)
    {

        float slope1 = (line1End.Y - line1Start.Y) / (line1End.X - line1Start.X);
        float slope2 = (line2End.Y - line2Start.Y) / (line2End.X - line2Start.X);

        float yinter1 = GetLineYIntesept(line1Start, slope1);
        float yinter2 = GetLineYIntesept(line2Start, slope2);

        if (slope1 == slope2 && yinter1 != yinter2)
            return PointF.Empty;

        float x = (yinter2 - yinter1) / (slope1 - slope2);

        float y = slope1 * x + yinter1;

        return new PointF(x, y);
    }
yazan
fuente
3

Porté la respuesta anterior de Kris a JavaScript. Después de probar numerosas respuestas diferentes, proporcionó los puntos correctos. Pensé que me estaba volviendo loco porque no estaba obteniendo los puntos que necesitaba.

function getLineLineCollision(p0, p1, p2, p3) {
    var s1, s2;
    s1 = {x: p1.x - p0.x, y: p1.y - p0.y};
    s2 = {x: p3.x - p2.x, y: p3.y - p2.y};

    var s10_x = p1.x - p0.x;
    var s10_y = p1.y - p0.y;
    var s32_x = p3.x - p2.x;
    var s32_y = p3.y - p2.y;

    var denom = s10_x * s32_y - s32_x * s10_y;

    if(denom == 0) {
        return false;
    }

    var denom_positive = denom > 0;

    var s02_x = p0.x - p2.x;
    var s02_y = p0.y - p2.y;

    var s_numer = s10_x * s02_y - s10_y * s02_x;

    if((s_numer < 0) == denom_positive) {
        return false;
    }

    var t_numer = s32_x * s02_y - s32_y * s02_x;

    if((t_numer < 0) == denom_positive) {
        return false;
    }

    if((s_numer > denom) == denom_positive || (t_numer > denom) == denom_positive) {
        return false;
    }

    var t = t_numer / denom;

    var p = {x: p0.x + (t * s10_x), y: p0.y + (t * s10_y)};
    return p;
}
Mono de código
fuente
2

Lo intenté de muchas maneras y luego decidí escribir el mío. Asi que aqui esta:

bool IsBetween (float x, float b1, float b2)
{
   return ( ((x >= (b1 - 0.1f)) && 
        (x <= (b2 + 0.1f))) || 
        ((x >= (b2 - 0.1f)) &&
        (x <= (b1 + 0.1f))));
}

bool IsSegmentsColliding(   POINTFLOAT lineA,
                POINTFLOAT lineB,
                POINTFLOAT line2A,
                POINTFLOAT line2B)
{
    float deltaX1 = lineB.x - lineA.x;
    float deltaX2 = line2B.x - line2A.x;
    float deltaY1 = lineB.y - lineA.y;
    float deltaY2 = line2B.y - line2A.y;

    if (abs(deltaX1) < 0.01f && 
        abs(deltaX2) < 0.01f) // Both are vertical lines
        return false;
    if (abs((deltaY1 / deltaX1) -
        (deltaY2 / deltaX2)) < 0.001f) // Two parallel line
        return false;

    float xCol = (  (   (deltaX1 * deltaX2) * 
                        (line2A.y - lineA.y)) - 
                    (line2A.x * deltaY2 * deltaX1) + 
                    (lineA.x * deltaY1 * deltaX2)) / 
                 ((deltaY1 * deltaX2) - (deltaY2 * deltaX1));
    float yCol = 0;
    if (deltaX1 < 0.01f) // L1 is a vertical line
        yCol = ((xCol * deltaY2) + 
                (line2A.y * deltaX2) - 
                (line2A.x * deltaY2)) / deltaX2;
    else // L1 is acceptable
        yCol = ((xCol * deltaY1) +
                (lineA.y * deltaX1) -
                (lineA.x * deltaY1)) / deltaX1;

    bool isCol =    IsBetween(xCol, lineA.x, lineB.x) &&
            IsBetween(yCol, lineA.y, lineB.y) &&
            IsBetween(xCol, line2A.x, line2B.x) &&
            IsBetween(yCol, line2A.y, line2B.y);
    return isCol;
}

Basado en estas dos fórmulas: (las simplifiqué de la ecuación de líneas y otras fórmulas)

fórmula para x

fórmula para y

Soroush Falahati
fuente
Funciona, pero intente ingresar esta coordenada (si es colineal / superpuesta, entonces devolverá un resultado falso): PointA1 = (0,0) PointA2 = (0,2) y PointB1 = (0,1) PointB2 = (0,5)
dns
@dns Bueno, eso se debe a que el código devuelve falso para líneas paralelas. Veo el problema, sin embargo, todavía no sé qué debería devolver la función, ya que hay un número infinito de respuestas.
Soroush Falahati
2

Esto se basa en la respuesta de Gareth Ree. También devuelve la superposición de los segmentos de línea si lo hacen. Codificado en C ++, V es una clase de vector simple. Donde el producto cruzado de dos vectores en 2D devuelve un solo escalar. Fue probado y aprobado por el sistema de prueba automática de mi escuela.

//Required input point must be colinear with the line
bool on_segment(const V& p, const LineSegment& l)
{
    //If a point is on the line, the sum of the vectors formed by the point to the line endpoints must be equal
    V va = p - l.pa;
    V vb = p - l.pb;
    R ma = va.magnitude();
    R mb = vb.magnitude();
    R ml = (l.pb - l.pa).magnitude();
    R s = ma + mb;
    bool r = s <= ml + epsilon;
    return r;
}

//Compute using vector math
// Returns 0 points if the lines do not intersect or overlap
// Returns 1 point if the lines intersect
//  Returns 2 points if the lines overlap, contain the points where overlapping start starts and stop
std::vector<V> intersect(const LineSegment& la, const LineSegment& lb)
{
    std::vector<V> r;

    //http://stackoverflow.com/questions/563198/how-do-you-detect-where-two-line-segments-intersect
    V oa, ob, da, db; //Origin and direction vectors
    R sa, sb; //Scalar values
    oa = la.pa;
    da = la.pb - la.pa;
    ob = lb.pa;
    db = lb.pb - lb.pa;

    if (da.cross(db) == 0 && (ob - oa).cross(da) == 0) //If colinear
    {
        if (on_segment(lb.pa, la) && on_segment(lb.pb, la))
        {
            r.push_back(lb.pa);
            r.push_back(lb.pb);
            dprintf("colinear, overlapping\n");
            return r;
        }

        if (on_segment(la.pa, lb) && on_segment(la.pb, lb))
        {
            r.push_back(la.pa);
            r.push_back(la.pb);
            dprintf("colinear, overlapping\n");
            return r;
        }

        if (on_segment(la.pa, lb))
            r.push_back(la.pa);

        if (on_segment(la.pb, lb))
            r.push_back(la.pb);

        if (on_segment(lb.pa, la))
            r.push_back(lb.pa);

        if (on_segment(lb.pb, la))
            r.push_back(lb.pb);

        if (r.size() == 0)
            dprintf("colinear, non-overlapping\n");
        else
            dprintf("colinear, overlapping\n");

        return r;
    }

    if (da.cross(db) == 0 && (ob - oa).cross(da) != 0)
    {
        dprintf("parallel non-intersecting\n");
        return r;
    }

    //Math trick db cross db == 0, which is a single scalar in 2D.
    //Crossing both sides with vector db gives:
    sa = (ob - oa).cross(db) / da.cross(db);

    //Crossing both sides with vector da gives
    sb = (oa - ob).cross(da) / db.cross(da);

    if (0 <= sa && sa <= 1 && 0 <= sb && sb <= 1)
    {
        dprintf("intersecting\n");
        r.push_back(oa + da * sa);
        return r;
    }

    dprintf("non-intersecting, non-parallel, non-colinear, non-overlapping\n");
    return r;
}
ColacX
fuente
2

Aquí hay una implementación básica de un segmento de línea en C #, con el código de detección de intersección correspondiente. Requiere una estructura de vector / punto 2D llamada Vector2f, aunque puede reemplazarla con cualquier otro tipo que tenga propiedades X / Y. También puede reemplazar floatcon doublesi eso se adapta mejor a sus necesidades.

Este código se usa en mi biblioteca de física .NET, Boing .

public struct LineSegment2f
{
    public Vector2f From { get; }
    public Vector2f To { get; }

    public LineSegment2f(Vector2f @from, Vector2f to)
    {
        From = @from;
        To = to;
    }

    public Vector2f Delta => new Vector2f(To.X - From.X, To.Y - From.Y);

    /// <summary>
    /// Attempt to intersect two line segments.
    /// </summary>
    /// <remarks>
    /// Even if the line segments do not intersect, <paramref name="t"/> and <paramref name="u"/> will be set.
    /// If the lines are parallel, <paramref name="t"/> and <paramref name="u"/> are set to <see cref="float.NaN"/>.
    /// </remarks>
    /// <param name="other">The line to attempt intersection of this line with.</param>
    /// <param name="intersectionPoint">The point of intersection if within the line segments, or empty..</param>
    /// <param name="t">The distance along this line at which intersection would occur, or NaN if lines are collinear/parallel.</param>
    /// <param name="u">The distance along the other line at which intersection would occur, or NaN if lines are collinear/parallel.</param>
    /// <returns><c>true</c> if the line segments intersect, otherwise <c>false</c>.</returns>
    public bool TryIntersect(LineSegment2f other, out Vector2f intersectionPoint, out float t, out float u)
    {
        var p = From;
        var q = other.From;
        var r = Delta;
        var s = other.Delta;

        // t = (q − p) × s / (r × s)
        // u = (q − p) × r / (r × s)

        var denom = Fake2DCross(r, s);

        if (denom == 0)
        {
            // lines are collinear or parallel
            t = float.NaN;
            u = float.NaN;
            intersectionPoint = default(Vector2f);
            return false;
        }

        var tNumer = Fake2DCross(q - p, s);
        var uNumer = Fake2DCross(q - p, r);

        t = tNumer / denom;
        u = uNumer / denom;

        if (t < 0 || t > 1 || u < 0 || u > 1)
        {
            // line segments do not intersect within their ranges
            intersectionPoint = default(Vector2f);
            return false;
        }

        intersectionPoint = p + r * t;
        return true;
    }

    private static float Fake2DCross(Vector2f a, Vector2f b)
    {
        return a.X * b.Y - a.Y * b.X;
    }
}
Drew Noakes
fuente
1

Un programa C ++ para verificar si dos segmentos de línea dados se cruzan

#include <iostream>
using namespace std;

struct Point
{
    int x;
    int y;
};

// Given three colinear points p, q, r, the function checks if
// point q lies on line segment 'pr'
bool onSegment(Point p, Point q, Point r)
{
    if (q.x <= max(p.x, r.x) && q.x >= min(p.x, r.x) &&
        q.y <= max(p.y, r.y) && q.y >= min(p.y, r.y))
       return true;

    return false;
}

// To find orientation of ordered triplet (p, q, r).
// The function returns following values
// 0 --> p, q and r are colinear
// 1 --> Clockwise
// 2 --> Counterclockwise
int orientation(Point p, Point q, Point r)
{
    // See 10th slides from following link for derivation of the formula
    // http://www.dcs.gla.ac.uk/~pat/52233/slides/Geometry1x1.pdf
    int val = (q.y - p.y) * (r.x - q.x) -
              (q.x - p.x) * (r.y - q.y);

    if (val == 0) return 0;  // colinear

    return (val > 0)? 1: 2; // clock or counterclock wise
}

// The main function that returns true if line segment 'p1q1'
// and 'p2q2' intersect.
bool doIntersect(Point p1, Point q1, Point p2, Point q2)
{
    // Find the four orientations needed for general and
    // special cases
    int o1 = orientation(p1, q1, p2);
    int o2 = orientation(p1, q1, q2);
    int o3 = orientation(p2, q2, p1);
    int o4 = orientation(p2, q2, q1);

    // General case
    if (o1 != o2 && o3 != o4)
        return true;

    // Special Cases
    // p1, q1 and p2 are colinear and p2 lies on segment p1q1
    if (o1 == 0 && onSegment(p1, p2, q1)) return true;

    // p1, q1 and p2 are colinear and q2 lies on segment p1q1
    if (o2 == 0 && onSegment(p1, q2, q1)) return true;

    // p2, q2 and p1 are colinear and p1 lies on segment p2q2
    if (o3 == 0 && onSegment(p2, p1, q2)) return true;

     // p2, q2 and q1 are colinear and q1 lies on segment p2q2
    if (o4 == 0 && onSegment(p2, q1, q2)) return true;

    return false; // Doesn't fall in any of the above cases
}

// Driver program to test above functions
int main()
{
    struct Point p1 = {1, 1}, q1 = {10, 1};
    struct Point p2 = {1, 2}, q2 = {10, 2};

    doIntersect(p1, q1, p2, q2)? cout << "Yes\n": cout << "No\n";

    p1 = {10, 0}, q1 = {0, 10};
    p2 = {0, 0}, q2 = {10, 10};
    doIntersect(p1, q1, p2, q2)? cout << "Yes\n": cout << "No\n";

    p1 = {-5, -5}, q1 = {0, 0};
    p2 = {1, 1}, q2 = {10, 10};
    doIntersect(p1, q1, p2, q2)? cout << "Yes\n": cout << "No\n";

    return 0;
}
Ayush Srivastava
fuente
1

Basado en la respuesta de @Gareth Rees, versión para Python:

import numpy as np

def np_perp( a ) :
    b = np.empty_like(a)
    b[0] = a[1]
    b[1] = -a[0]
    return b

def np_cross_product(a, b):
    return np.dot(a, np_perp(b))

def np_seg_intersect(a, b, considerCollinearOverlapAsIntersect = False):
    # /programming/563198/how-do-you-detect-where-two-line-segments-intersect/565282#565282
    # http://www.codeproject.com/Tips/862988/Find-the-intersection-point-of-two-line-segments
    r = a[1] - a[0]
    s = b[1] - b[0]
    v = b[0] - a[0]
    num = np_cross_product(v, r)
    denom = np_cross_product(r, s)
    # If r x s = 0 and (q - p) x r = 0, then the two lines are collinear.
    if np.isclose(denom, 0) and np.isclose(num, 0):
        # 1. If either  0 <= (q - p) * r <= r * r or 0 <= (p - q) * s <= * s
        # then the two lines are overlapping,
        if(considerCollinearOverlapAsIntersect):
            vDotR = np.dot(v, r)
            aDotS = np.dot(-v, s)
            if (0 <= vDotR  and vDotR <= np.dot(r,r)) or (0 <= aDotS  and aDotS <= np.dot(s,s)):
                return True
        # 2. If neither 0 <= (q - p) * r = r * r nor 0 <= (p - q) * s <= s * s
        # then the two lines are collinear but disjoint.
        # No need to implement this expression, as it follows from the expression above.
        return None
    if np.isclose(denom, 0) and not np.isclose(num, 0):
        # Parallel and non intersecting
        return None
    u = num / denom
    t = np_cross_product(v, s) / denom
    if u >= 0 and u <= 1 and t >= 0 and t <= 1:
        res = b[0] + (s*u)
        return res
    # Otherwise, the two line segments are not parallel but do not intersect.
    return None
Ibraim Ganiev
fuente
0

Si cada lado del rectángulo es un segmento de línea, y la porción dibujada por el usuario es un segmento de línea, entonces solo debe verificar la intersección del segmento dibujado por el usuario con los cuatro segmentos de línea lateral. Este debería ser un ejercicio bastante simple dados los puntos de inicio y finalización de cada segmento.

Harper Shelby
fuente
3
Tenga en cuenta que esta fue una respuesta razonable a la pregunta como se planteó originalmente, pero ahora que la pregunta ha sido editada en gran medida no tiene mucho sentido.
GS - Pídele disculpas a Monica el
0

Basado en la respuesta de t3chb0t:

int intersezione_linee(int x1, int y1, int x2, int y2, int x3, int y3, int x4, int y4, int& p_x, int& p_y)
{
   //L1: estremi (x1,y1)(x2,y2) L2: estremi (x3,y3)(x3,y3)
   int d;
   d = (x1-x2)*(y3-y4) - (y1-y2)*(x3-x4);
   if(!d)
       return 0;
   p_x = ((x1*y2-y1*x2)*(x3-x4) - (x1-x2)*(x3*y4-y3*x4))/d;
   p_y = ((x1*y2-y1*x2)*(y3-y4) - (y1-y2)*(x3*y4-y3*x4))/d;
   return 1;
}

int in_bounding_box(int x1, int y1, int x2, int y2, int p_x, int p_y)
{
    return p_x>=x1 && p_x<=x2 && p_y>=y1 && p_y<=y2;

}

int intersezione_segmenti(int x1, int y1, int x2, int y2, int x3, int y3, int x4, int y4, int& p_x, int& p_y)
{
    if (!intersezione_linee(x1,y1,x2,y2,x3,y3,x4,y4,p_x,p_y))
        return 0;

    return in_bounding_box(x1,y1,x2,y2,p_x,p_y) && in_bounding_box(x3,y3,x4,y4,p_x,p_y);
}
volperossa
fuente
0

Leí estos algoritmos del libro "geometría de vista múltiple"

siguiente texto usando

'como signo de transposición

* como producto punteado

x como producto cruzado, cuando se usa como operador

1. definición de línea

un punto x_vec = (x, y) 'se encuentra en la línea ax + by + c = 0

denotamos L = (a, b, c) ', el punto como (x, y, 1)' como coordenadas homogéneas

la ecuación lineal se puede escribir como

(x, y, 1) (a, b, c) '= 0 o x' * L = 0

2. intersección de líneas

tenemos dos líneas L1 = (a1, b1, c1) ', L2 = (a2, b2, c2)'

suponga que x es un punto, un vector yx = L1 x L2 (producto cruzado L1 L2).

tenga cuidado, x es siempre un punto 2D, lea las coordenadas homogéneas si está confundido acerca de (L1xL2) es un vector de tres elementos, y x es una coordenadas 2D.

según triple producto, sabemos que

L1 * (L1 x L2) = 0, y L2 * (L1 x L2) = 0, debido al coplano L1, L2

sustituimos (L1xL2) con el vector x, entonces tenemos L1 * x = 0, L2 * x = 0, lo que significa que x está en L1 y L2, x es el punto de intersección.

tenga cuidado, aquí x es coordenadas homogéneas, si el último elemento de x es cero, significa que L1 y L2 son paralelas.

Mass Zhou
fuente
0

Muchas respuestas han incluido todos los cálculos en una sola función. Si necesita calcular las pendientes de la línea, las intersecciones con el eje y las intersecciones con el eje x para usar en otra parte de su código, realizará esos cálculos de forma redundante. He separado las funciones respectivas, he usado nombres de variables obvios y he comentado mi código para que sea más fácil de seguir. Necesitaba saber si las líneas se cruzan infinitamente más allá de sus puntos finales, así que en JavaScript:

http://jsfiddle.net/skibulk/evmqq00u/

var point_a = {x:0, y:10},
    point_b = {x:12, y:12},
    point_c = {x:10, y:0},
    point_d = {x:0, y:0},
    slope_ab = slope(point_a, point_b),
    slope_bc = slope(point_b, point_c),
    slope_cd = slope(point_c, point_d),
    slope_da = slope(point_d, point_a),
    yint_ab = y_intercept(point_a, slope_ab),
    yint_bc = y_intercept(point_b, slope_bc),
    yint_cd = y_intercept(point_c, slope_cd),
    yint_da = y_intercept(point_d, slope_da),
    xint_ab = x_intercept(point_a, slope_ab, yint_ab),
    xint_bc = x_intercept(point_b, slope_bc, yint_bc),
    xint_cd = x_intercept(point_c, slope_cd, yint_cd),
    xint_da = x_intercept(point_d, slope_da, yint_da),
    point_aa = intersect(slope_da, yint_da, xint_da, slope_ab, yint_ab, xint_ab),
    point_bb = intersect(slope_ab, yint_ab, xint_ab, slope_bc, yint_bc, xint_bc),
    point_cc = intersect(slope_bc, yint_bc, xint_bc, slope_cd, yint_cd, xint_cd),
    point_dd = intersect(slope_cd, yint_cd, xint_cd, slope_da, yint_da, xint_da);

console.log(point_a, point_b, point_c, point_d);
console.log(slope_ab, slope_bc, slope_cd, slope_da);
console.log(yint_ab, yint_bc, yint_cd, yint_da);
console.log(xint_ab, xint_bc, xint_cd, xint_da);
console.log(point_aa, point_bb, point_cc, point_dd);

function slope(point_a, point_b) {
  var i = (point_b.y - point_a.y) / (point_b.x - point_a.x);
  if (i === -Infinity) return Infinity;
  if (i === -0) return 0;
  return i;
}

function y_intercept(point, slope) {
    // Horizontal Line
    if (slope == 0) return point.y;
  // Vertical Line
    if (slope == Infinity)
  {
    // THE Y-Axis
    if (point.x == 0) return Infinity;
    // No Intercept
    return null;
  }
  // Angled Line
  return point.y - (slope * point.x);
}

function x_intercept(point, slope, yint) {
    // Vertical Line
    if (slope == Infinity) return point.x;
  // Horizontal Line
    if (slope == 0)
  {
    // THE X-Axis
    if (point.y == 0) return Infinity;
    // No Intercept
    return null;
  }
  // Angled Line
  return -yint / slope;
}

// Intersection of two infinite lines
function intersect(slope_a, yint_a, xint_a, slope_b, yint_b, xint_b) {
  if (slope_a == slope_b)
  {
    // Equal Lines
    if (yint_a == yint_b && xint_a == xint_b) return Infinity;
    // Parallel Lines
    return null;
  }
  // First Line Vertical
    if (slope_a == Infinity)
  {
    return {
        x: xint_a,
      y: (slope_b * xint_a) + yint_b
    };
  }
  // Second Line Vertical
    if (slope_b == Infinity)
  {
    return {
        x: xint_b,
      y: (slope_a * xint_b) + yint_a
    };
  }
  // Not Equal, Not Parallel, Not Vertical
  var i = (yint_b - yint_a) / (slope_a - slope_b);
  return {
    x: i,
    y: (slope_a * i) + yint_a
  };
}
skibulk
fuente