//-----------------------------------【程序说明】----------------------------------------------
// 程序名称::《获取图像重心》
// 开发所用系统:Windows 8.1 Pro
// 开发所用IDE版本:Visual Studio 2013 Ultimate
// 开发所用OpenCV版本: 2.4.8
// 2014年5月6日 Integtated by Thunthersine
//------------------------------------------------------------------------------------------------
#include "stdafx.h"
#include "cv.h"
#include "highgui.h"
#include "cvaux.h"
#include "cxcore.h"
#include <stdio.h>
#include <vector>
using namespace std;
using namespace cv;
void main()
{
double m00, m10, m01, x, y; //定义x,y为重心坐标。
IplImage* img1 = cvLoadImage("tyre.jpg", 1);
IplImage* img1_gray = cvCreateImage(cvGetSize(img1), 8, 1);
cvCvtColor(img1, img1_gray, CV_BGR2GRAY);
IplImage* img_edge1 = cvCreateImage(cvGetSize(img1), 8, 1);
cvThreshold(img1_gray, img_edge1, 240, 255, CV_THRESH_BINARY);
CvMemStorage* storage1 = cvCreateMemStorage();
CvSeq* contour1 = NULL;
int outline = cvFindContours(img_edge1, storage1, &contour1, sizeof(CvContour), CV_RETR_LIST); //获取轮廓数
printf("轮廓数:%d\n", outline);
//轮廓的矩
CvMoments *moments1 = new CvMoments();
cvMoments(contour1, moments1, 0);
//Hu矩
CvHuMoments *huMonents1 = new CvHuMoments();
cvGetHuMoments(moments1, huMonents1);
printf("遍历三阶矩:\n");
for (int xOrder = 0; xOrder <= 3; xOrder++)
for (int yOrder = 0; yOrder <= 3; yOrder++)
{
if (xOrder + yOrder <= 3)
{
double spatialMoment = cvGetSpatialMoment(moments1, xOrder, yOrder);
printf(" %d,%d : 空间距-%.2f\n", xOrder, yOrder, spatialMoment);
}
}
m00 = cvGetSpatialMoment(moments1, 0, 0);
m10 = cvGetSpatialMoment(moments1, 1, 0);
m01 = cvGetSpatialMoment(moments1, 0, 1);
x = (int)(m10 / m00);
y = (int)(m01 / m00);
printf("轮胎中心坐标为:%f,%f\n", x,y);
}