使用opencv作物件识别（一） —— 积分直方图加速HOG特征计算

博客分类：
图像识别、机器学习、数据挖掘

CC++C#

方向梯度直方图（Histograms of Oriented Gradients，简称HOG特征）结合支持向量机（ support vector machine，简称SVM），被广泛应用于图像识别中，尤其在行人检测中获得了极大的成功。

积分直方图可以用于快速计算原始图像矩形区域内的HOG特征。积分直方图的概念类似与viola和jones在脸部识别中所用的积分图像。

下面的代码给出了，对于一幅给定的图像，如何快速计算积分直方图，以及如何使用其进行HOG特征的演算（关键处以给出注释）：

C代码

/*Function to calculate the integral histogram*/
IplImage** calculateIntegralHOG(IplImage* in)
{
/*Convert the input image to grayscale*/
IplImage* img_gray = cvCreateImage(cvGetSize(in), IPL_DEPTH_8U,1);
cvCvtColor(in, img_gray, CV_BGR2GRAY);
cvEqualizeHist(img_gray,img_gray);
/*Calculate the derivates of the grayscale image in the x and y directions using a sobel operator and obtain 2 gradient images for the x and y directions*/
IplImage *xsobel, *ysobel;
xsobel = doSobel(img_gray, 1, 0, 3);
ysobel = doSobel(img_gray, 0, 1, 3);
cvReleaseImage(&img_gray);
/* Create an array of 9 images (9 because I assume bin size 20 degrees and unsigned gradient ( 180/20 = 9), one for each bin which will have zeroes for all pixels, except for the pixels in the original image for which the gradient values correspond to the particular bin. These will be referred to as bin images. These bin images will be then used to calculate the integral histogram, which will quicken the calculation of HOG descriptors */
IplImage** bins = (IplImage**) malloc(9 * sizeof(IplImage*));
for (int i = 0; i < 9 ; i++) {
bins[i] = cvCreateImage(cvGetSize(in), IPL_DEPTH_32F,1);
cvSetZero(bins);
}
/* Create an array of 9 images ( note the dimensions of the image, the cvIntegral() function requires the size to be that), to store the integral images calculated from the above bin images. These 9 integral images together constitute the integral histogram */
IplImage** integrals = (IplImage**) malloc(9 * sizeof(IplImage*));
for (int i = 0; i < 9 ; i++) {
integrals[i] = cvCreateImage(cvSize(in->width + 1, in->height + 1),
IPL_DEPTH_64F,1);
}
/* Calculate the bin images. The magnitude and orientation of the gradient at each pixel is calculated using the xsobel and ysobel images.{Magnitude = sqrt(sq(xsobel) + sq(ysobel) ), gradient = itan (ysobel/xsobel) }. Then according to the orientation of the gradient, the value of the corresponding pixel in the corresponding image is set */
int x, y;
float temp_gradient, temp_magnitude;
for (y = 0; y < in->height; y++) {
/* ptr1 and ptr2 point to beginning of the current row in the xsobel and ysobel images respectively. ptrs point to the beginning of the current rows in the bin images */
float* ptr1 = (float*) (xsobel->imageData + y * (xsobel->widthStep));
float* ptr2 = (float*) (ysobel->imageData + y * (ysobel->widthStep));
float** ptrs = (float**) malloc(9 * sizeof(float*));
for (int i = 0; i < 9 ;i++){
ptrs[i] = (float*) (bins[i]->imageData + y * (bins->widthStep));
}
/*For every pixel in a row gradient orientation and magnitude are calculated and corresponding values set for the bin images. */
for (x = 0; x <in->width; x++) {
/* if the xsobel derivative is zero for a pixel, a small value is added to it, to avoid division by zero. atan returns values in radians, which on being converted to degrees, correspond to values between -90 and 90 degrees. 90 is added to each orientation, to shift the orientation values range from {-90-90} to {0-180}. This is just a matter of convention. {-90-90} values can also be used for the calculation. */
if (ptr1[x] == 0){
temp_gradient = ((atan(ptr2[x] / (ptr1[x] + 0.00001))) * (180/ PI)) + 90;
}
else{
temp_gradient = ((atan(ptr2[x] / ptr1[x])) * (180 / PI)) + 90;
}
temp_magnitude = sqrt((ptr1[x] * ptr1[x]) + (ptr2[x] * ptr2[x]));
/*The bin image is selected according to the gradient values. The corresponding pixel value is made equal to the gradient magnitude at that pixel in the corresponding bin image */
if (temp_gradient <= 20) {
ptrs[0][x] = temp_magnitude;
}
else if (temp_gradient <= 40) {
ptrs[1][x] = temp_magnitude;
}
else if (temp_gradient <= 60) {
ptrs[2][x] = temp_magnitude;
}
else if (temp_gradient <= 80) {
ptrs[3][x] = temp_magnitude;
}
else if (temp_gradient <= 100) {
ptrs[4][x] = temp_magnitude;
}
else if (temp_gradient <= 120) {
ptrs[5][x] = temp_magnitude;
}
else if (temp_gradient <= 140) {
ptrs[6][x] = temp_magnitude;
}
else if (temp_gradient <= 160) {
ptrs[7][x] = temp_magnitude;
}
else {
ptrs[8][x] = temp_magnitude;
}
}
}
cvReleaseImage(&xsobel);
cvReleaseImage(&ysobel);
/*Integral images for each of the bin images are calculated*/
for (int i = 0; i <9 ; i++){
cvIntegral(bins[i], integrals[i]);
}
for (int i = 0; i <9 ; i++){
cvReleaseImage(&bins[i]);
}
/*The function returns an array of 9 images which consitute the integral histogram*/
return (integrals);
}

如何使用上面的函数来计算9维的方向梯度直方图呢？如下：

C代码

/* The following function takes as input the rectangular cell for which the histogram of oriented gradients has to be calculated, a matrix hog_cell of dimensions 1x9 to store the bin values for the histogram, the integral histogram, and the normalization scheme to be used. No normalization is done if normalization = -1 */
void calculateHOG_rect(CvRect cell, CvMat* hog_cell,
IplImage** integrals, int normalization) {
/* Calculate the bin values for each of the bin of the histogram one by one */
for (int i = 0; i < 9 ; i++){
float a =((double*)(integrals[i]->imageData + (cell.y) * (integrals->
widthStep)))[cell.x];
float b = ((double*) (integrals->imageData + (cell.y + cell.height) *
(integrals->widthStep)))[cell.x + cell.width];
float c = ((double*) (integrals->imageData + (cell.y) * (integrals-
>widthStep)))[cell.x + cell.width];
float d = ((double*) (integrals->imageData + (cell.y + cell.height) *
(integrals->widthStep)))[cell.x];
((float*) hog_cell->data.fl) = (a + b) - (c + d);
}
/*Normalize the matrix*/
if (normalization != -1){
cvNormalize(hog_cell, hog_cell, 1, 0, normalization);
}
}

使用opencv作物件识别（一） —— 积分直方图加速HOG特征计算相关推荐

OpenCV For Python识别图片中的米粒个数，并计算米粒的面积平均面积
精确分割请看分水岭算法实现https://blog.csdn.net/qq_36623595/article/details/109273629 使用OpenCV构建图像识别算法,识别图片中的米粒个数 ...
python opencv 人体/人脸识别简易demo
前言参考文章: pythonopencv检测行人_[图像处理]使用OpenCV实现人脸和行人检测 Python如何实现行人识别-人体识别本文主要讲述关于opencv官方提供的现成模型,采取Hog特 ...
OpenCV：使用python-cv2+Hog特征+SVM实现狮子识别
文章目录 SVM Hog特征 Hog特征+SVM实现狮子识别 SVM 支持向量机:寻求一个最优的超平面,实现样本的分类下面我们用SVM实现一个根据身高体重对男女生分类的问题 import cv2 i ...
svm性别识别及HOG特征原理浅析
HOG特征原理主要参考http://blog.csdn.net/abcjennifer/article/details/7365651 HOG即histogram of oriented gra ...
基于LBP纹理特征计算GLCM的纹理特征统计量+SVM/RF识别纹理图片
作者 | Raini 出品 | 北京图特摩斯科技有限公司 (thutmose.cn) 说在前头: 局部特征检测方法斑点Blob检测,LoG检测 , DoG,DoH检测,SIFT算法,SUFT算法边 ...
youcans 的 OpenCV 学习课—6.灰度变换与直方图处理
youcans 的 OpenCV 学习课-6.灰度变换与直方图处理本系列面向 Python 小白,从零开始实战解说 OpenCV 项目实战. 空间域的图像处理方法直接对图像的像素点进行处理,空间域图 ...
OpenCV 识别图片中的米粒个数，并计算米粒的平均面积和长度（转）
介绍 OpenCV+Python 使用OpenCV构建图像识别算法,识别图片中的米粒个数,并计算米粒的平均面积和长度软件架构模块:OpenCV 4.0.0.21 编程语言:Python 3.7.2 ...
OpenCV 识别图片中的米粒个数，并计算米粒的平均面积和长度
介绍 OpenCV+Python 使用OpenCV构建图像识别算法,识别图片中的米粒个数,并计算米粒的平均面积和长度软件架构模块:OpenCV 4.0.0.21 编程语言:Python 3.7.2 ...
【OpenCV 4开发详解】直方图应用
本文首发于"小白学视觉"微信公众号,欢迎关注公众号本文作者为小白,版权归人民邮电出版社发行所有,禁止转载,侵权必究! 经过几个月的努力,小白终于完成了市面上第一本OpenCV 4 ...

使用opencv作物件识别（一） —— 积分直方图加速HOG特征计算

使用opencv作物件识别（一） —— 积分直方图加速HOG特征计算

使用opencv作物件识别（一） —— 积分直方图加速HOG特征计算相关推荐

最新文章

热门文章