利用HOG特征和SVM算法对图片进行二分类实现图像目标的检测

一、HOG方向梯度直方图

1.HOG简介

    方向梯度直方图（Histogram of Oriented Gradient, HOG）特征是一种在计算机视觉和图像处理中用来进行物体检测的特征描述子。它通过计算和统计图像局部区域的梯度方向直方图来构成特征。Hog特征结合SVM分类器已经被广泛应用于图像识别中。HOG的思想就是在一副图像中，局部目标的表象和形状（appearance and shape）能够被梯度或边缘的方向密度分布很好地描述。

2.HOG基本步骤简述

(1)标准化gamma空间和颜色空间

    第一步，为了减少光照因素的影响，首先需要将整个图像进行规范化（归一化）。在图像的纹理强度中，局部的表层曝光贡献的比重较大，所以，这种压缩处理能够有效地降低图像局部的阴影和光照变化。因为颜色信息作用不大，通常先转化为灰度图；

(2)计算图像梯度

    第二步，计算图像横坐标和纵坐标方向的梯度，并据此计算每个像素位置的梯度方向值；求导操作不仅能够捕获轮廓，人影和一些纹理信息，还能进一步弱化光照的影响。

(3)为每个细胞单元构建梯度方向直方图

    第三步，为每个细胞单元构建梯度方向直方图的目的是为局部图像区域提供一个编码，同时能够保持对图像中人体对象的姿势和外观的弱敏感性。我们将图像分成若干个“单元格cell”，例如每个cell为66个像素。假设我们采用9个bin的直方图来统计这66个像素的梯度信息。但bin并不是连续的，而是20度为一个bin，然后相同bin是关于180度对称。

(4)把细胞单元组合成大的块（block），块内归一化梯度直方图

    第四步，由于局部光照的变化以及前景-背景对比度的变化，使得梯度强度的变化范围非常大。这就需要对梯度强度做归一化。归一化能够进一步地对光照、阴影和边缘进行压缩,方法是把各个细胞单元组合成大的、空间上连通的区间（blocks）。例如本次实验所使用的block是(16, 16)，窗口（即图片）大小是(64, 128)，block的滑动步长是(8,8)，共有105个block，block内的cell是(8, 8)，共有4个cell，而cell内的bin是9个，那么总共有36715=3780个特征。

(5)收集HOG特征

    最后一步就是将检测窗口中所有重叠的块进行HOG特征的收集，并将它们结合成最终的特征向量供分类使用，在opencv中使用相应的方法。图1 HOG检测步骤流程。

图1 HOG检测步骤流程

二、SVM支持向量机

    支持向量机（支持向量机（Support Vector Machine, SVM）是一类按监督学习(supervised learning)方式对数据进行二元分类的广义线性分类器(generalized linear classifier)，将样本分为正样本和负样本，其决策边界是对学习样本求解的最大边距超平面(maximum-margin hyperplane)。SVM一般用于解决模式识别领域中的数据分类问题，属于有监督学习算法的一种。

三、实验流程

1.数据集准备

    数据集分为正负样本数据集，2个文件夹，分别存储在pos（pos1为作者的图片数据集）和neg文件夹下，正样本有820张图片，负样本有1931张图片，均是64*128的格式。

2.图像预处理

    对图像进行3以下3步的预处理：首先将视频分解为图片，其次是图片缩放，最后是将图片裁剪为64 * 128。

3.实验结果

    图2和图3分别是我仿造源作者代码的实验结果以及作者的实验结果，明显的看出我的实验结果并不准确。

图2 自己的实验结果

图3 原作的实验结果

四、程序代码

1.视频图像的预处理

(1)视频分解

'''
视频分解图片 video_Decomposition(video_name)
1 load 2 info 3 parse 4 imshow imwrite
负样本 64*128 正样本：64*128
'''
import cv2
import numpy as np
def video_Decomposition(video_name, begin):
    cap = cv2.VideoCapture(video_name)# 获取一个视频打开cap 1 file name
    isOpened = cap.isOpened # 判断是否打开
    print(isOpened)
    fps = cap.get(cv2.CAP_PROP_FPS)#帧率
    width = int(cap.get(cv2.CAP_PROP_FRAME_WIDTH)) #width
    height = int(cap.get(cv2.CAP_PROP_FRAME_HEIGHT))#height
    print(fps, width, height)
    i = 0 + 205 * (begin - 1) 
    while(isOpened):
        if(i == 205 * (begin - 1) + 205):
            break
        else:
            i = i + 1
        (flag, frame) = cap.read() #读取每一章 flag frame
        fileName = 'pos_Original//'  + str(i) + '.jpg'
        print(fileName)
        if flag == True:
            cv2.imwrite(fileName, frame, [cv2.IMWRITE_JPEG_QUALITY,100])
    print(i)
    print("video_Decomposition OK!")
for i in range(1, 5):
    video_Decomposition(str(i) + '.mp4', i)

(2)图片裁剪缩放