Opencv计算一幅图像的hu矩

Opencv计算一幅图像的hu矩

首先Hu矩是干嘛得呢？用hu矩来匹配图像

下面是计算hu矩的过程  概念先了解一下吧

普通矩 -----》
0阶矩（m00） ：目标区域的质量
1阶矩（m01，m10） ：目标区域的质心
2阶矩（m02，m20，m11） ：旋转半径
3阶矩（m03，m30，m12，m21） ：目标的方位和斜度，反应目标

中心矩：构造平移不变性 ---------》 
质心坐标：x = m10/m00 y = m01/m00；

归一化中心矩:构造尺度不变性 ----------》（mu02，mu20，mu11，mu03，mu30，mu12，mu21）

hu矩：构造旋转不变性 （h1，h2，。。。。，h7）

代码如下：

#include<iostream>
#include<opencv2/core/core.hpp>
#include<opencv2/highgui/highgui.hpp>
//计算一幅图像的hu矩
int main()
{
	IplImage* src;
	src = cvCreateImage(cvSize(10,10),8,1); //这里是创建一幅图像 也可以自己加载一幅图像计算
	cvZero(src); //把所以值都赋值为0
	int xx,yy;
	for( yy = 0;yy<5;yy++){
		for( xx = 0;xx<5;xx++){
			cvSetReal2D(src,yy,xx,255); //把前5行5列的像素变为255 白色
		}
	}
	/*cvNamedWindow("1");
	cvShowImage("1",src);*/
	double m00,m10,m01; //计算0阶矩和1阶矩
	CvMoments moment;
	cvMoments(src,&moment,2);  //2： 只要比0大的数都可以 表示像素255的点都用1来表示
	m00 = cvGetSpatialMoment(&moment,0,0);
	m10 = cvGetSpatialMoment(&moment,1,0);
	m01 = cvGetSpatialMoment(&moment,0,1);

	double x = m10/m00;
	double y = m01/m00;

	CvHuMoments humoment;
	cvGetHuMoments(&moment,&humoment); //计算出hu矩的值

	cvWaitKey(0);
	return 0;
}