边缘方向直方图

//
//输入：
//	image: 灰度图像
//	dst：边缘图像
//	edge_hist_feature：边缘直方图特征
// 	#define EDGE_HIST_SIZE 16
//	double edge_feature[EDGE_HIST_SIZE];
int calEdgeHistogram(IplImage *image,IplImage * dst,double * edge_hist_feature)
{
	if(NULL ==  image || NULL == edge_hist_feature)
		return 1;

	CvHistogram *hist = 0; // 直方图
	IplImage* canny;//边缘图像
	IplImage* gradient_im;
	IplImage* dx; // x方向的sobel差分
	IplImage* dy; // y方向的sobel差分 

	CvMat* canny_m;
	CvMat* gradient; // 梯度值
	CvMat* gradient_dir; //梯度的方向
	CvMat* dx_m; // 
	CvMat* dy_m;
	CvMat* mask;
	CvSize  size;

	int i,j;
	float theta;
	float max_val;

	int hdims = EDGE_HIST_SIZE;     // 划分HIST的个数，越高越精确
	float hranges_arr[] = {-PI/2,PI/2}; // 直方图的上界和下界
	float* hranges = hranges_arr;

	size = cvGetSize(image);
	canny = cvCreateImage(cvGetSize(image),8,1);//边缘图像
	dx = cvCreateImage(cvGetSize(image),32,1);//x方向上的差分