opencv中最小外接矩形2 的算法

Cvbox2D

本文来自CSDN博客，转载请标明出处：http://blog.csdn.net/mine1024/archive/2010/11/30/6044856.aspx

之前用到opencv最小外接矩形去表示一个类椭圆形的高度，特此记录备查。

对给定的 2D 点集，寻找最小面积的包围矩形，使用函数：

CvBox2D     cvMinAreaRect2(   const   CvArr *   points,   CvMemStorage *   storage = NULL   );

   points

   点序列或点集数组

   storage

   可选的临时存储仓

  函数 cvMinAreaRect2 通过建立凸外形并且旋转外形以寻找给定 2D 点集的最小面积的包围矩形。

其中返回的2D盒子定义如下：

1 typedef   struct   CvBox2D
2 {
3     CvPoint2D32f   center;   /*   盒子的中心   */
4     CvSize2D32f   size;   /*   盒子的长和宽   */
5     float   angle;   /*   水平轴与第一个边的夹角，用弧度表示 */
6 } CvBox2D;

注意夹角 angle 是水平轴逆时针旋转，与碰到的第一个边（不管是高还是宽）的夹角。 如下图

                                 

  可用函数 cvBoxPoints(box[count], point); 寻找盒子的顶点

1  void   cvBoxPoints(   CvBox2D   box,   CvPoint2D32f   pt[ 4 ]   )
2  {
3      double   angle   =   box . angle * CV_PI / 180 .
4      float   a   =   ( float )cos(angle) * 0 . 5f;
5      float   b   =   ( float )sin(angle) * 0 . 5f;
6 
7      pt[ 0 ] . x   =   box . center . x   -   a * box . size . height   -   b * box . size . width;
8      pt[ 0 ] . y   =   box . center . y   +   b * box . size . height   -   a * box . size . width;
9      pt[ 1 ] . x   =   box . center . x   +   a * box . size . height   -   b * box . size . width;
10     pt[ 1 ] . y   =   box . center . y   -   b * box . size . height   -   a * box . size . width;
11     pt[ 2 ] . x   =   2 * box . center . x   -   pt[ 0 ] . x;
12     pt[ 2 ] . y   =   2 * box . center . y   -   pt[ 0 ] . y;
13     pt[ 3 ] . x   =   2 * box . center . x   -   pt[ 1 ] . x;
14     pt[ 3 ] . y   =   2 * box . center . y   -   pt[ 1 ] . y;
15 }

简单证明此函数的计算公式:

 

   计算 x，由图可得到三个方程式： pt[ 1 ] . x   -   pt[ 0 ] . x   =   width * sin(angle)
                             pt[ 2 ] . x   -   pt[ 1 ] . x   =   height * cos(angle)
                             pt[ 2 ] . x   -   pt[ 0 ] . x   =   2 (box . center . x   -   pt[ 0 ] . x)

   联立方程可解得函数里的计算式，算 y 略。

写了个函数绘制CvBox2D

1  void   DrawBox(CvBox2D   box,IplImage *   img)
2  {
3      CvPoint2D32f   point[ 4 ];
4      int   i;
5      for   (   i = 0 ;   i < 4;   i + + )
6      {
7          point[i] . x   =   0 ;
8          point[i] . y   =   0 ;
9      }
10     cvBoxPoints(box,   point);   // 计算二维盒子顶点
11     CvPoint   pt[ 4 ];
12     for   (   i = 0 ;   i < 4;   i + + )
13     {
14         pt[i] . x   =   ( int )point[i] . x;
15         pt[i] . y   =   ( int )point[i] . y;
16     }
17     cvLine(   img,   pt[ 0 ],   pt[ 1 ],CV_RGB( 255 , 0 , 0 ),   2 ,   8 ,   0   );
18     cvLine(   img,   pt[ 1 ],   pt[ 2 ],CV_RGB( 255 , 0 , 0 ),   2 ,   8 ,   0   );
19     cvLine(   img,   pt[ 2 ],   pt[ 3 ],CV_RGB( 255 , 0 , 0 ),   2 ,   8 ,   0   );
20     cvLine(   img,   pt[ 3 ],   pt[ 0 ],CV_RGB( 255 , 0 , 0 ),   2 ,   8 ,   0   );
21 }