当前位置：首页 > 编程日记 > 正文

利用OpenCV求取图像的重心

编程日记 2024-12-30 05:20:00

转自：http://blog.csdn.net/lxiaoxiaot/article/details/6539834

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

这个随机变量的属性可以用统计特征--矩（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获得指定维的矩信息。

例如：

// 计算二值化图像imgYellowThresh中物体的重心坐标

CvMoments *moments = (CvMoments*)malloc(sizeof(CvMoments));

cvMoments(imgYellowThresh, moments, 1); // The actual moment values double moment10 = cvGetSpatialMoment(moments, 1, 0); double moment01 = cvGetSpatialMoment(moments, 0, 1); double area = cvGetSpatialMoment(moments, 0, 0);

// 计算重心

posX = moment10 / area;

posY = moment01 / area;

转自：http://blog.163.com/forever_871226/blog/static/34424308201141851736984/

/** 计算二值图像的重心
* @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 &center)
{
//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;
}

转自：http://www.opencv.org.cn/forum/viewtopic.php?p=1016

#include "cv.h"
#include "highgui.h"

void main( int argc, char** argv )
{
IplImage* src;
CvMoments moments;
CvMat *region;
CvPoint pt1,pt2;
double m00 = 0, m10, m01, mu20, mu11, mu02, inv_m00;
double a, b, c;
int xc, yc;
// 第一条命令行参数确定了图像的文件名。
if( (src=cvLoadImage("two.bmp", 0))!= 0)
//if( (src=cvLoadImage("fbb.jpg", 0))!= 0)
//if( argc == 2 && (src=cvLoadImage(argv[1], 0))!= 0)
{
IplImage* dst = cvCreateImage( cvGetSize(src), 8, 3 );
CvMemStorage* storage = cvCreateMemStorage(0);
CvSeq* contour = 0;

cvThreshold( src, src, 100, 255, CV_THRESH_BINARY );//100 is the thredhold
cvNot( src, src );
cvNamedWindow( "Source", 1 );
cvShowImage( "Source", src );

//cvFindContours( src, storage, &contour, sizeof(CvContour), CV_RETR_CCOMP, CV_CHAIN_APPROX_SIMPLE,cvPoint(0,0) );

cvFindContours( src, storage, &contour, sizeof(CvContour),CV_RETR_TREE, CV_CHAIN_APPROX_SIMPLE, cvPoint(0,0) );
cvZero( dst );

for( ; contour != 0; contour = contour->h_next )
{
//CvScalar color = CV_RGB( rand()&255, rand()&255, rand()&255 );
CvScalar color = CV_RGB( 255, 0,0 );
/* 用1替代 CV_FILLED 所指示的轮廓外形 */
cvDrawContours( dst, contour, color, color, -1, CV_FILLED, 8,cvPoint(0,0) );//you can change 1 to CV_FILLED
contour = cvApproxPoly( contour, sizeof(CvContour), storage, CV_POLY_APPROX_DP, 3, 1 );
//CvRect* r = (CvRect*)cvGetSeqElem( contour,1);
region=(CvMat*)contour;
cvMoments( region, &moments,0 );
//cvMoments( &contour, &moments,0 );
// cvDrawContours( cnt_img, _contours, CV_RGB(255,0,0), CV_RGB(0,255,0), _levels, 3, CV_AA, cvPoint(0,0) ); CV_FILLED

/
m00 = moments.m00;
m10 = moments.m10;
m01 = moments.m01;
mu11 = moments.mu11;
mu20 = moments.mu20;
mu02 = moments.mu02;

//if( fabs(m00) < DBL_EPSILON )break;

inv_m00 = 1. / m00;
xc = cvRound( m10 * inv_m00 );
yc = cvRound( m01 * inv_m00 );
a = mu20 * inv_m00;
b = mu11 * inv_m00;
c = mu02 * inv_m00;
/
pt1.x=xc-1;pt1.y=yc;
pt2.x=xc+1;pt2.y=yc;
cvLine( dst, pt1, pt2, CV_RGB(0,255,0), 2, CV_AA, 0 );
pt1.x=xc;pt1.y=yc-1;
pt2.x=xc;pt2.y=yc+1;
cvLine( dst, pt1, pt2, CV_RGB(0,255,0), 2, CV_AA, 0 );
}

cvNamedWindow( "Components", 1 );
cvShowImage( "Components", dst );
cvWaitKey(0);
}
}

https://www.dkcj.cn/info/34048.html

利用OpenCV求取图像的重心

相关文章：

C++中Ansi、Unicode、UTF8字符串之间的转换和写入

MyBatis的扩展点(plugins)

linux中使用CST时间

到「黄埔学院」去：打造AI首席架构师，第二期限量招募！

linux守护进程的创建

tesseract3.01的训练和使用

旷视推出鼻纹识别，用AI寻找丢失宠物

Qt浅谈之一：内存泄露（总结）

LINUX新手入门-1.装系统

Log4cplus1.04的使用

FRVT赛程全纪录：格灵深瞳全球排名前五

反转比特位(文章最后有干货)【转】

过关斩将打进Kaggle竞赛Top 0.3%，我是这样做的

JavaBean规范

vigra1.8.0的使用

17个Python小窍门

用产品思路建设中台，这走得通吗？| 白话中台

利用cvMinAreaRect2求取轮廓最小外接矩形

电脑开机显示Invalidsystemdisk

Windows7 64位下vs2008配置OpenCV2.3.1

腾讯拥抱开源：首次公布开源路线图，技术研发向共享、复用和开源迈进

Babylon.js 3.3发布：更强大的粒子系统和WebVR支持

基于Erlang语言的视频相似推荐系统 | 深度

MinGW简介

Confluence 6 创建小组的公众空间

windows 7 可以清除的文件

阿里最新论文解读：考虑时空域影响的点击率预估模型DSTN

windows7 64位机上配置MinGW+Codeblocks+ wxWidgets

jQuery带动画的弹出对话框

陶哲轩实分析习题 13.4.6