任意椭圆的外凸包矩形算法

How to get the limits of rotated ellipse?

任意椭圆的外凸包矩形算法

问题来源：人脸检测数据集FDDB，涉及到将原标签中的椭圆形状转化为矩形框用于人脸检测，其中就有生成斜椭圆的外接包围矩形的方法，结果如下图。

注意：不是下图的最小面积矩形，最小面积矩形指两对边分别与主轴和副轴平行的斜矩形，参见下图。

二次曲线的方程表达式

首先从二次曲线的方程式入手，它可以写成以下形式：
$$A{x}^2+2Bxy+C{y}^2+2Dx+2Ey+F=0$$

其中针对椭圆的标准方程：
$$\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$$

如果 $a=b$ 上述方程就是半径为 $a$ 的圆，所以椭圆有时候称为变形的圆。

标准椭圆通过平移和旋转等几何变换后，

$$\begin{gathered} x=x^{\prime }\cos \theta -y^{\prime }\sin \theta +h, \\ y=x^{\prime }\sin \theta +y^{\prime }\cos \theta +k \end{gathered}$$

此标准方程变形为:
$$A^{\prime }{x}^{\prime 2}+2B^{\prime }{x}^{\prime }{y}^{\prime }+C^{\prime }{y}^{\prime 2}+2D^{\prime }{x}^{\prime }+2E^{\prime }{y}^{\prime }+F^{\prime }=0$$

其中

$A^{\prime }=A{\cos }^2\theta +2B\cos \theta \sin \theta +C{\sin }^2\theta$
$2B^{\prime }=2B\cos 2\theta -(A-C)\sin 2\theta$
$C^{\prime }=A{\sin }^2\theta -2B\cos \theta \sin \theta +C{\cos }^2\theta$
$2D^{\prime }=2Ah\cos \theta +2B(k\cos \theta +h\sin \theta )+2Ck\sin \theta +2D\cos \theta +2E\sin \theta$
$2E^{\prime }=-2Ah\sin \theta +2B(h\cos \theta -k\sin \theta )+2Ck\cos \theta -2D\sin \theta +2E\cos \theta$
$F^{\prime }=A{h}^2+2Bhk+C{k}^2+2Dh+2Ek+F$

二次曲线方程的系数关系

记录

$$\begin{gathered} A+C=A^{\prime }+C^{\prime }=H \\ {B}^2-AC=B^{\prime 2}-A^{\prime }{C}^{\prime }=\delta \\ \left|\begin{array}{ccc}A& B& D\\ B& C& E\\ D& E& F\end{array}\right|=\left|\begin{array}{ccc}{A}^{\prime }& B^{\prime }& D^{\prime }\\ B^{\prime }& C^{\prime }& E^{\prime }\\ D^{\prime }& E^{\prime }& F^{\prime }\end{array}\right|=\Delta \end{gathered}$$

是二次曲线方程系数之间的三个基本不变量。

二次曲线分类定理

1. 当 $\Delta =0$ 时, 二次曲线蜕化
   1. $\delta <0$ 时，点椭圆
   2. $\delta >0$ 时，相交两直线
   3. $\delta =0$ 时，平行两直线（包括相重）
2. 当 $\Delta \ne 0$ 时
   1. 当 $\delta <0$ 时，二次曲线为椭圆（有可能为虚椭圆）。
   2. 当 $\delta >0$ 时，二次曲线为双曲线。
   3. 当 $\delta =0$ 时，二次曲线为抛物线。

在Geogebra中，如何获取二次曲线方程或多项式系数

Coefficients(conic)

Coefficients command, that displays the coefficients of the conic, using the form

$$a{x}^2+b{y}^2+c+dxy+ex+fy=0$$

$A=a,B=d/2.0,C=b,D=e/2.0,E=f/2.0,F=c$

针对一般方程计算椭圆的边界

假设 二次曲线的方程为

$$A{x}^2+2Bxy+C{y}^2+2Dx+2Ey+F=0$$

题目要求：该椭圆的最小包围矩形（凸包矩形）

思考路径：求出斜椭圆的 $x_{min},x_{max},y_{min},y_{max}$

(1) 用直线 $y=k$ 去截椭圆, 相交就是两个交点，相切就是一个交点，不相交就是无解。
故可以通过一元二次方程判别式来求解。
显然，我们要求相切的情形。
椭圆方程可以看成关于$x$的一元二次方程:
$$A{x}^2+2(By+D)x+(C{y}^2+2Ey+F)=0$$

由判别式 $\Delta =0$, 我们得到：
$$(By+D{)}^2-A(C{y}^2+2Ey+F)=0$$

整理得到：
$$(B^2-AC)y^2+2(BD-AE)y+(D^2-AF)=0$$

记 $\delta =B^2-AC,\beta =BD-AE,\gamma =D^2-AF$

则可以求得到：
$$y_{max}=\frac{-\beta +\sqrt{{\beta }^2-\delta \gamma }}{\delta }$$
$$y_{min}=\frac{-\beta -\sqrt{{\beta }^2-\delta \gamma }}{\delta }$$

(2) 同理，用直线 $x=h$ 去截椭圆, 相交就是两个交点，相切就是一个交点，不相交就是无解。
椭圆方程可以看成关于$y$的一元二次方程:
$$C{y}^2+2(Bx+E)x+(A{x}^2+2Dy+F)=0$$

由判别式 $\Delta =0$, 我们得到：
$$(Bx+E{)}^2-C(A{x}^2+2Dx+F)=0$$

整理得到：
$$(B^2-AC)x^2+2(BE-CD)x+(E^2-CF)=0$$

记 $\delta =B^2-AC,\alpha =BE-CD,\lambda =E^2-CF$

则可以求得：
$$x_{max}=\frac{-\alpha +\sqrt{{\alpha }^2-\delta \lambda }}{\delta }$$
$$x_{min}=\frac{-\alpha -\sqrt{{\alpha }^2-\delta \lambda }}{\delta }$$

(3) 最后得到矩形的四个顶点为 $(x_{min},y_{min}),(x_{max},y_{min}),(x_{max},y_{max}),(x_{min},y_{max})$

用 Geogebra 得到的结果

利用函数 Coefficients(conic) 得到系数
转换后，直接计算即可得到包围倾斜椭圆的正矩形。

命令 Coefficients(椭圆ellipse)：得到列表 { a , b , c , d , e , f } \{a,b,c,d,e,f\} {a,b,c,d,e,f}, 代表椭圆方程为 $a{x}^2+b{y}^2+c+dxy+ex+fy=0$

离心率

到定点（焦点）的距离与到定直线（准线）的距离的商是常数 e（离心率）的点的轨迹。

• 当e>1时，为双曲线的一支，
• 当e=1时，为抛物线，
• 当 0< e <1 时，为椭圆，
• 当e=0时，为一点或一个圆。

离心率求法 $e=\frac{c}{a}$

在Geogebra软件中，还可以用命令 e=Eccentricity( 圆锥曲线 ) 得到圆锥曲线的离心率。
从而得到：圆锥曲线上一点的焦半径长度等于该点到相应准线的距离乘以离心率， $c=ae$。
在椭圆中有 $c^2=a^2-b^2,a^2+b^2=2a^2-c^2=a^2(2-e^2)$

蒙日圆

蒙日圆上任一点到圆锥曲线的两条切线是相互垂直的。

椭圆的蒙日圆的圆心就是椭圆的中心，半径就是 $r=\sqrt{a^2+b^2}=a\sqrt{2-e^2}$

利用蒙日圆得到椭圆外凸包矩形

参见 GeoGebra网站

面积差值对比

参考文献

1. 如何计算包围倾斜椭圆的最小矩形

2. How to get the limits of rotated ellipse?