求解物体重心

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


using namespace std;
int main( int argc, char** argv )
{
IplImage* src;  
src=cvLoadImage("222.png",CV_LOAD_IMAGE_GRAYSCALE); 
    // IplImage* result = cvCreateImage( cvGetSize(src), 8, 3 ); 
IplImage* dst = cvCreateImage( cvGetSize(src), 8, 3 );  
CvMemStorage* storage = cvCreateMemStorage(0);  
CvSeq* contour = 0; 
IplImage* result;
result=cvCloneImage(src);
cvThreshold( src, src,120, 255, CV_THRESH_BINARY );//二值化   
cvNamedWindow( "Source", 1 );  
cvShowImage( "Source", src );  
IplImage*pContourImg= cvCreateImage( cvGetSize(src), 8, 1 ); 
cvCopy(src,pContourImg);

 // 提取轮廓  
    int contour_num = cvFindContours(pContourImg, storage, &contour, sizeof(CvContour), CV_RETR_CCOMP, CV_CHAIN_APPROX_SIMPLE);  
    cvZero(dst);        // 清空数组  
    CvSeq *_contour = contour;   
    double maxarea = 10000;  
    double minarea = 10;  
    int m = 0;  
    for( ; contour != 0; contour = contour->h_next )    
    {    
  
        double tmparea = fabs(cvContourArea(contour));  
        if(tmparea < minarea)     
        {    
            cvSeqRemove(contour, 0); // 删除面积小于设定值的轮廓  
            continue;  
        }    
        CvRect aRect = cvBoundingRect( contour, 0 );   
        if ((aRect.width/aRect.height)<1)    
        {    
            cvSeqRemove(contour, 0); //删除宽高比例小于设定值的轮廓  
            continue;  
        }    
        if(tmparea > maxarea)    
        {    
            maxarea = tmparea;  
        }    
        m++;  
        // 创建一个色彩值  
    //    CvScalar color = CV_RGB( 0, 0, 255 );  
  
     /*   max_level 绘制轮廓的最大等级。如果等级为0，绘制单独的轮廓。如果为1，绘制轮廓及在其后的相同的级别下轮廓  
        如果值为2，所有的轮廓。如果等级为2，绘制所有同级轮廓及所有低一级轮廓，诸此种种  
        如果值为负数，函数不绘制同级轮廓，但会升序绘制直到级别为abs(max_level)-1的子轮廓 */ 
     //   cvDrawContours(dst, contour, color, color, 0, 1, 8);   //绘制外部和内部的轮廓  
    }    
    contour = _contour;  
    int count = 0;  double tmparea=0;
    for(; contour != 0; contour = contour->h_next)  
    {    
        count++;  
         tmparea = fabs(cvContourArea(contour));  
        if (tmparea >= maxarea)    
        {    
            CvScalar color = CV_RGB( 0, 255, 0);  
            cvDrawContours(dst, contour, color, color, -1, 1, 8);  
cout<<"222"<<endl;
cout<<"面积为"<<tmparea<<endl;
cout<<endl;


CvRect aRect = cvBoundingRect( contour, 0 ); 
//找重心
{
CvPoint2D32f center = cvPoint2D32f(0, 0);
int countOfPoint = 0;
for(int i = aRect.x; i < aRect.x + aRect.width; ++i){
for(int j = aRect.y; j < aRect.y + aRect.height; ++j){
if(*(src->imageData + src->widthStep * j + i) != 0){
center.x += i;
center.y += j;
countOfPoint++;
}
}
}


center.x /= countOfPoint;
center.y /= countOfPoint;
cvCircle(dst, cvPoint(center.x, center.y), 5, cvScalar(0, 255), 2);
}













        }    
    }    
    printf("%f\r\n", tmparea );
//CvRect aRect = cvBoundingRect( contmax, 0 );
//cvSetImageROI( pContourImg,aRect);
 //显示图像
 cvShowImage( "contour", dst );
 cvShowImage( "result",result);
 cvSaveImage("contour.bmp",pContourImg);
 cvWaitKey(0);


 //销毁窗口
 cvDestroyWindow( "src" );
 cvDestroyWindow( "result" );
 cvDestroyWindow( "contour" );
 //释放图像
 cvReleaseImage( &src );
 cvReleaseImage( &result );
 cvReleaseImage( &pContourImg );
 cvReleaseMemStorage(&storage);


 return 0;





 return 0;
}