OpenCV：获取图像重心

//-----------------------------------【程序说明】----------------------------------------------
//        程序名称:：《获取图像重心》
//        开发所用系统：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);
}