emgucv 图像重心

原創

碎片球球

2020-02-21 11:30

不规则区域的矩，表示把一个归一化的灰度级图像函数理解为一个二维随机变量的概率密度。

这个随机变量的属性可以用统计特征--矩（Moments）来描述。通过假设非零的像素值表示区域，矩可以用于二值或灰度级的区域描述。

M_pq = sigma(i)sigma(j) i^pj^qf(i,j)

其中x,y,i,j是区域点的座标（在数字图像中的像素座标）。

令Xc,Yc表示区域重心的座标，则：

Xc = M₁₀/M₀₀;

Yc = M₀₁/M₀₀;

在二值图像的情况下，M₀₀表示区域的面积。

OpenCV中可以使用函数cvMoments来计算二值图像的矩信息。

使用函数cvGetSpatialMoment获得指定维的矩信息。

以上介绍转自http://blog.csdn.net/lxiaoxiaot/article/details/6539834

以下代码转自博客http://blog.163.com/forever_871226/blog/static/34424308201141851736984/ 使用OpenCV计算图像重心

/** 计算二值图像的重心
* @param[in] src  输入的待处理图像
* @param[out] center 重心座标
* @retval 0  操作成功
* @retval -1 操作失败
* @note 输入图像是二值化图像
* @note xc=M10/M00, yc=M01/M00, 其中 Mx_order,y_order=SUMx,y(I(x,y)*x^x_order*y^y_order)
*/
static int aoiGravityCenter(IplImage *src, CvPoint ¢er)
{
 //if(!src)
 // return GRAVITYCENTER__SRC_IS_NULL;
 double m00, m10, m01;
 CvMoments moment;
 cvMoments( src, &moment, 1);
 m00 = cvGetSpatialMoment( &moment, 0, 0 );
 if( m00 == 0) 
  return 1;
 m10 = cvGetSpatialMoment( &moment, 1, 0 );
 m01 = cvGetSpatialMoment( &moment, 0, 1 );
 center.x = (int) (m10/m00);
 center.y = (int) (m01/m00);
 return 0;
}

以下代码是按照此方式进行转换的EmguCV图像重心求法，并进行了绘制。

                        Image<Gray, Byte> imageCenter = imageProcessG.GetSubRect(new System.Drawing.Rectangle(item.X, item.Y, item.Width, item.Height));
                        double m00, m01, m10;
                        MCvMoments moment = new MCvMoments();
                        CvInvoke.cvMoments(imageCenter.Ptr, ref moment, 1);
                        m00 = CvInvoke.cvGetSpatialMoment(ref moment, 0, 0);

                        m10 = CvInvoke.cvGetSpatialMoment(ref moment, 1, 0);
                        m01 = CvInvoke.cvGetSpatialMoment(ref moment, 0, 1);

                        //System.Drawing.PointF center = new System.Drawing.PointF(item.Center.X, item.Center.Y);
                        System.Drawing.PointF center = new System.Drawing.PointF(item.X+(float)(m10 / m00), item.Y+(float)(m01 / m00));
                        Ellipse tempE = new Ellipse(center, new System.Drawing.SizeF(2, 2), (float)Math.PI);
                        imageResults.Draw(tempE, new Bgr(System.Drawing.Color.Red), 2);