二维几何计算-直线与圆、二次曲线交点计算

---

1.直线与二次曲线交点

1.1 方程

二次曲线：$X^{T}AX+B^{T}X+c=0\cdots (1)$,
直线： $X(t)=P+t\vec{d}\cdots （2）$

1.2 交点

直接将(2)直接带入（1）可得：
$X(t{)}^{T}AX(t)+B^{T}X(t)+c=0$
$=>({\vec{d}}^{T}A\vec{d})t^2+{\vec{d}}^T(2AP+B)t+(P^{T}AP+B^T+c)=0$
$=>e_2{t}^2+e_{1}t+e_0=0$
此时，通过对二次根判别式进行求解，当然如果只是要做根判别计算的话，可以免去平方根计算，提高计算效率。

1.3 关于射线，线段

同理1.2 我们可以在考虑定以域后自行去推理。

2.直线与圆交点计算

2.1方程

Circle: ${\Vert X-C\Vert }^2=r^2\cdots (1)$
Line : $X(t)=P+t\vec{d}\cdots (2)$

2.2 交点

将（2） 带入(1) 可得：
$${\Vert \vec{d}\Vert }^2{t}^2+2\Vert \vec{d}\Vert \overrightarrow{△}t+{\Vert \overrightarrow{△}\Vert }^2-r^2=0$$
其中$\overrightarrow{△}=P-C$。
可以通过二次求根公式求出交点：
$$t=\frac{-\vec{d}\cdot \overrightarrow{△}\pm \sqrt{(\vec{d}\cdot \overrightarrow{△}{)}^2-{\Vert \overrightarrow{△}\Vert }^2({\Vert \overrightarrow{△}\Vert }^2-r^2)}}{{\Vert \vec{d}\Vert }^2}$$
令$\delta =(\vec{d}\cdot \overrightarrow{△}{)}^2-{\Vert \overrightarrow{△}\Vert }^2({\Vert \overrightarrow{△}\Vert }^2-r^2)$。若$\delta <0$,则不相交；若等于0，则有一个交点；若大于0，则有2个交点。

2.3 考虑直线为射线的情况

在我们求出t值后，需要做t>0判断即可。

2.4 考虑直线为线段情况

线段方程：$P_0+t\vec{d}=0,\vec{d}=P_1-P_0$,判断$t\in [0,1]$

2.5 考虑圆为圆弧的case

若有交点K，弧线AB应与K在直线AB的同一侧，如下：

判断依据为： $KR(K-A,B-A)\ge 0$,KR定义参照。