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;
}