opencv3计算图像轮廓的中心矩

#include<iostream>
#include<opencv2/opencv.hpp>
#include<vector>

using namespace cv;
using namespace std;

int main()
{
Mat srcImage = imread("group.jpg");
imshow("【原图】", srcImage);

//首先对图像进行空间的转换
Mat grayImage;
cvtColor(srcImage, grayImage, CV_BGR2GRAY);
//对灰度图进行滤波
GaussianBlur(grayImage, grayImage, Size(3, 3), 0, 0);
imshow("【滤波后的图像】", grayImage);

//为了得到二值图像，对灰度图进行边缘检测
Mat cannyImage;
Canny(grayImage, cannyImage, 128, 255, 3);
//在得到的二值图像中寻找轮廓
vector<vector<Point>> contours;
vector<Vec4i> hierarchy;
findContours(cannyImage, contours, hierarchy, RETR_TREE, CHAIN_APPROX_SIMPLE, Point(0, 0));

//绘制轮廓
for (int i = 0; i < (int)contours.size(); i++)
{
drawContours(cannyImage, contours, i, Scalar(255), 1, 8);
}

//计算轮廓的矩
vector<Moments> g_vMoments(contours.size());
for (int i = 0; i < (int)contours.size(); i++)
{
g_vMoments[i] = moments(contours[i], true);
}
//输出所有轮廓的矩
for (int i = 0; i < (int)contours.size(); i++)
{
cout << "【用矩计算出来的第" << i << "个轮廓的面积为：】" << g_vMoments[i].m00 << endl;
}

/**************************************************中心矩计算公式**************************************************/
//利用计算得到的 矩 计算中心矩
vector<Point2f> centerMoments(contours.size());
for (int i = 0; i < (int)contours.size(); i++)
{
centerMoments[i] = Point2f(float(g_vMoments[i].m10 / g_vMoments[i].m00), float(g_vMoments[i].m01 / g_vMoments[i].m00));
}
//将得到的中心矩显示出来
for (int i = 0; i < (int)contours.size(); i++)
{
//参数中的 -1 表示填充
circle(cannyImage, (Point)centerMoments[i], 5, Scalar(255), -1, 8);
}

imshow("【处理后的图像】", cannyImage);

waitKey(0);

return 0;
}