图像区域的最小二乘拟合椭圆与偏心率

禹晶、肖创柏、廖庆敏《数字图像处理（面向新工科的电工电子信息基础课程系列教材）》P312

对区域最小二乘拟合椭圆，它的形状和方向能够很好地代表该区域的形状和方向。

计算$x$轴与椭圆长轴之间的角度（以度为单位），

1. 计算椭圆的长轴和短轴：

   • 椭圆的长轴和短轴可以通过计算图像区域的二阶矩（协方差矩阵）得到。
   • 协方差矩阵的特征向量和特征值分别对应于椭圆的主轴方向和轴长。

2. 确定长轴的方向：

   • 协方差矩阵的最大特征值对应的特征向量表示椭圆的长轴方向。

3. 计算角度：

   • 假设 $x$ 轴正方向为水平向右（0 度），长轴方向向量与 $x$ 轴之间的夹角即为所求角度。
   • 可以使用反正切函数（arctan2）来计算两个向量之间的夹角。

具体步骤

根据图像区域计算二阶矩（$m_{00},m_{10},m_{01},m_{11},m_{20},m_{02}$）

1. 计算图像区域的协方差矩阵：
   假设我们有一个图像区域 $I(x,y)$，其中 $x$ 和 $y$ 是像素的位置坐标。我们可以通过以下步骤计算该区域的协方差矩阵：

(1) 计算质心：首先计算图像区域的质心 $(\bar{x},\bar{y})$。
$$\bar{x}=\frac{{\sum }_{x,y}x\cdot I(x,y)}{{\sum }_{x,y}I(x,y)}=m_{10}$$
$$\bar{y}=\frac{{\sum }_{x,y}y\cdot I(x,y)}{{\sum }_{x,y}I(x,y)}=m_{01}$$

(2) 计算协方差矩阵：协方差矩阵 $C$ 是一个 2x2 的矩阵，其元素如下：
$$C=\left(\begin{array}{cc}{\sigma }_{xx}& {\sigma }_{xy}\\ {\sigma }_{yx}& {\sigma }_{yy}\end{array}\right)$$
其中：
$${\sigma }_{xx}=\frac{{\sum }_{x,y}(x-\bar{x}{)}^2\cdot I(x,y)}{{\sum }_{x,y}I(x,y)}=\frac{m_{20}}{m_{00}}$$
$${\sigma }_{yy}=\frac{{\sum }_{x,y}(y-\bar{y}{)}^2\cdot I(x,y)}{{\sum }_{x,y}I(x,y)}=\frac{m_{02}}{m_{00}}$$
$${\sigma }_{xy}={\sigma }_{yx}=\frac{{\sum }_{x,y}(x-\bar{x})(y-\bar{y})\cdot I(x,y)}{{\sum }_{x,y}I(x,y)}=\frac{m_{11}}{m_{00}}$$

2. 求解协方差矩阵的特征值和特征向量：

   • 求解特征方程 $|C-\lambda I|=0$ 得到特征值 ${\lambda }_1,{\lambda }_2$（${\lambda }_1\ge {\lambda }_2$）。
   • 对应的特征向量 ${\mathbf{v}}_1,{\mathbf{v}}_2$，其中 ${\mathbf{v}}_1$ 对应长轴方向。

3. 计算角度：

   • 特征向量 ${\mathbf{v}}_1=(v_{1x},v_{1y})$。
   • 角度 $\theta$ 为：
     $$\theta =\arctan 2(v_{1y},v_{1x})$$
   • 注意，arctan2 返回的角度范围是 $[-\pi ,\pi ]$，需要将其转换为 $[-90^{\circ },90^{\circ }]$：
     $$\theta =\{\begin{array}{ll}\theta \times \frac{180}{\pi }& \text{if }\theta \ge 0\\ (\theta +2\pi )\times \frac{180}{\pi }& \text{if }\theta <0\end{array}$$
     • 若 $\theta >90^{\circ }$，则 $\theta =180^{\circ }-\theta$。

一个例子

禹晶、肖创柏、廖庆敏《数字图像处理》资源二维码