7个hu不变矩c++程序小试

 我使用OpenCV中的函数cvGetHumoments（）计算7个hu不变矩时始终有错误，为了得到结果我自己用C++编写了一段计算7个hu不变矩的程序。

说明：图片必须是等边长的，即宽和高相等。****或许图片宽不能大于40像素，如果是边缘检测后的线条图像，图像宽高的限制会降低。如果你有什么好的想法，不妨留言交流一下。

#include "cv.h"

#include "cxcore.h"

#include "highgui.h"

#include <iostream.h>

float MOMENTS_M(IplImage* MAT,int p,int q)        //空间距

{

       CvScalar scl;

       float M_sum_y;

       float M_sum_xy;

 

       int H,W;

       H=MAT->height;

       W=MAT->width;

 

       for(int i=0;i<H;i++)

       {     for(int j=0;j<W;j++)

              {

                     scl=cvGet2D(MAT,i,j);

                     M_sum_y+=scl.val[0]*pow(j,q);

              }

              M_sum_xy+=pow(i,p)*M_sum_y;

       }

       return M_sum_xy;

}         

 

float MOMENTS_U(IplImage* MAT,int p,int q,float X_,float Y_)   //中心距

{

       CvScalar scl;

       float N_sum_y;

       float N_sum_xy;

 

       int H,W;

       H=MAT->height;

       W=MAT->width;

 

       for(int i=0;i<H;i++)

       {     for(int j=0;j<W;j++)

              {

                     scl=cvGet2D(MAT,i,j);

                     N_sum_y+=scl.val[0]*pow(j-Y_,q);

              }

              N_sum_xy+=pow(i-X_,p)*N_sum_y;

       }

       return N_sum_xy;

}          

//归一化的中心距

 

 

void main()

{

       int M00,M01,M10,M11,M02,M20,M12,M21,M03,M30;

       float X_,Y_;

       float U00,U02,U20,U11,U03,U30,U21,U12;

       float N02,N20,N11,N03,N30,N12,N21;

       float H1,H2,H3,H4,H5,H6,H7;

 

//     IplImage* img;

       IplImage* mat;

 

       mat=cvLoadImage("pic5_00.jpg",0);

//     mat=cvCreateImage(img->height,img->width,8);

 

//     cvCopy(img,mat);

 

       M00=MOMENTS_M(mat,0,0);

       M01=MOMENTS_M(mat,0,1);

       M10=MOMENTS_M(mat,1,0);

       M11=MOMENTS_M(mat,1,1);

       M02=MOMENTS_M(mat,0,2);

       M20=MOMENTS_M(mat,2,0);

       M12=MOMENTS_M(mat,1,2);

       M21=MOMENTS_M(mat,2,1);

       M03=MOMENTS_M(mat,0,3);

       M30=MOMENTS_M(mat,3,0);

       X_=M10/M00;

       Y_=M01/M00;

 

       U00=M00;

       U02=M02-M01*Y_;

       U20=M20-M10*X_;

       U11=M11-M10*Y_;

       U03=M03-3*M02*X_+2*M10*X_*X_;

       U30=M30-3*M02*Y_+2*M01*Y_*Y_;

       U12=M12-2*M11*Y_-M02*X_+2*M10*Y_*Y_;

       U21=M21-2*M11*X_-M20*Y_+2*M01*X_*X_;

       cout<<U00<<"XXX"<<U02<<"XXX"<<U20<<"XXX"<<U11<<"XXX"<<U03<<"XXX"<<U30<<"XXX"<<U12<<"XXX"<<U21<<endl;

 

       N02=U02/U00/U00;

       N20=U20/U00/U00;

       N11=U11/U00/U00;

       N03=U03/U00/U00/sqrt(U00);

       N30=U30/U00/U00/sqrt(U00);

       N12=U12/U00/U00/sqrt(U00);

       N21=U21/U00/U00/sqrt(U00);

 

       cout<<N02<<"XXX"<<N20<<"XXX"<<N11<<"XXX"<<N03<<"XXX"<<N30<<"XXX"<<N12<<"XXX"<<N21<<endl;

 

//hu7个不变距

       H1=N20+N02;

       H2=(N20-N02)* (N20-N02)+4*N11*N11;

       H3=(N30-3*N12)* (N30-3*N12)+(3*N21-N03)* (3*N21-N03);

       H4=(N30+N12)* (N30+N12)+(N21+N03)* (N21+N03);

       H5=((N30-3*N12)*(N30+N12)*((N30+N12)*(N30+N12)-3*(N21+N03)*(N21+N03))

              +(3*N21-N03)*(N21+N03)*(3*(N30+N12)* (N30+N12)-(N21+N03)* (N21+N03)))*100000;

       H6=(N20-N02)*((N30+N12)* (N30+N12)-(N21+N03)* (N21+N03))

              +4*N11*(N30+N12)*(N21+N03);

       H7=((3*N21-N03)*(N30+N12)*((N30+N12)*(N30+N12)-3*(N21+N03)*(N21+N03))

              +(3*N12-N30)*(N21+N03)*(3*(N30+N12)* (N30+N12)- (N21+N03)* (N21+N03)))*100000;

 

       cout<<H1<<"XXX"<<H2<<"XXX"<<H3<<"XXX"<<H4<<"XXX"<<H5<<"XXX"<<H6<<"XXX"<<H7<<endl;

}