目标检测的图像特征提取之(一)HOG特征

1HOG特征:html

       方向梯度直方图(Histogram of Oriented Gradient, HOG)特征是一种在计算机视觉和图像处理中用来进行物体检测的特征描述子。它经过计算和统计图像局部区域的梯度方向直方图来构成特征。Hog特征结合SVM分类器已经被普遍应用于图像识别中,尤为在行人检测中得到了极大的成功。须要提醒的是,HOG+SVM进行行人检测的方法是法国研究人员Dalal2005CVPR上提出的,而现在虽然有不少行人检测算法不断提出,但基本都是以HOG+SVM的思路为主。算法

1)主要思想:app

       在一副图像中,局部目标的表象和形状(appearance and shape)可以被梯度或边缘的方向密度分布很好地描述。(本质:梯度的统计信息,而梯度主要存在于边缘的地方)。ide

2)具体的实现方法是:函数

       首先将图像分红小的连通区域,咱们把它叫细胞单元。而后采集细胞单元中各像素点的梯度的或边缘的方向直方图。最后把这些直方图组合起来就能够构成特征描述器。性能

3)提升性能:学习

       把这些局部直方图在图像的更大的范围内(咱们把它叫区间或block)进行对比度归一化(contrast-normalized),所采用的方法是:先计算各直方图在这个区间(block)中的密度,而后根据这个密度对区间中的各个细胞单元作归一化。经过这个归一化后,能对光照变化和阴影得到更好的效果。编码

4)优势:spa

       与其余的特征描述方法相比,HOG有不少优势。首先,因为HOG是在图像的局部方格单元上操做,因此它对图像几何的和光学的形变都能保持很好的不变性,这两种形变只会出如今更大的空间领域上。其次,在粗的空域抽样、精细的方向抽样以及较强的局部光学归一化等条件下,只要行人大致上可以保持直立的姿式,能够允许行人有一些细微的肢体动做,这些细微的动做能够被忽略而不影响检测效果。所以HOG特征是特别适合于作图像中的人体检测的。orm

 

2HOG特征提取算法的实现过程:

大概过程:

HOG特征提取方法就是将一个image(你要检测的目标或者扫描窗口):

1)灰度化(将图像看作一个x,y,z(灰度)的三维图像);

2)采用Gamma校订法对输入图像进行颜色空间的标准化(归一化);目的是调节图像的对比度,下降图像局部的阴影和光照变化所形成的影响,同时能够抑制噪音的干扰;

3)计算图像每一个像素的梯度(包括大小和方向);主要是为了捕获轮廓信息,同时进一步弱化光照的干扰。

4)将图像划分红小cells(例如6*6像素/cell);

5)统计每一个cell的梯度直方图(不一样梯度的个数),便可造成每一个celldescriptor

6)将每几个cell组成一个block(例如3*3cell/block),一个block内全部cell的特征descriptor串联起来便获得该blockHOG特征descriptor

7)将图像image内的全部blockHOG特征descriptor串联起来就能够获得该image(你要检测的目标)的HOG特征descriptor了。这个就是最终的可供分类使用的特征向量了。

 

 

具体每一步的详细过程以下:

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

     为了减小光照因素的影响,首先须要将整个图像进行规范化(归一化)。在图像的纹理强度中,局部的表层曝光贡献的比重较大,因此,这种压缩处理可以有效地下降图像局部的阴影和光照变化。由于颜色信息做用不大,一般先转化为灰度图;

     Gamma压缩公式:

     好比能够取Gamma=1/2

 

2)计算图像梯度

        计算图像横坐标和纵坐标方向的梯度,并据此计算每一个像素位置的梯度方向值;求导操做不只可以捕获轮廓,人影和一些纹理信息,还能进一步弱化光照的影响。

图像中像素点(x,y)的梯度为:

       最经常使用的方法是:首先用[-1,0,1]梯度算子对原图像作卷积运算,获得x方向(水平方向,以向右为正方向)的梯度份量gradscalx,而后用[1,0,-1]T梯度算子对原图像作卷积运算,获得y方向(竖直方向,以向上为正方向)的梯度份量gradscaly。而后再用以上公式计算该像素点的梯度大小和方向。

 

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

        第三步的目的是为局部图像区域提供一个编码,同时可以保持对图像中人体对象的姿式和外观的弱敏感性。

咱们将图像分红若干个“单元格cell”,例如每一个cell6*6个像素。假设咱们采用9bin的直方图来统计这6*6个像素的梯度信息。也就是将cell的梯度方向360度分红9个方向块,如图所示:例如:若是这个像素的梯度方向是20-40度,直方图第2bin的计数就加一,这样,对cell内每一个像素用梯度方向在直方图中进行加权投影(映射到固定的角度范围),就能够获得这个cell的梯度方向直方图了,就是该cell对应的9维特征向量(由于有9bin)。

        像素梯度方向用到了,那么梯度大小呢?梯度大小就是做为投影的权值的。例如说:这个像素的梯度方向是20-40度,而后它的梯度大小是2(假设啊),那么直方图第2bin的计数就不是加一了,而是加二(假设啊)。

         细胞单元能够是矩形的(rectangular),也能够是星形的(radial)。

 

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

       因为局部光照的变化以及前景-背景对比度的变化,使得梯度强度的变化范围很是大。这就须要对梯度强度作归一化。归一化可以进一步地对光照、阴影和边缘进行压缩。

        做者采起的办法是:把各个细胞单元组合成大的、空间上连通的区间(blocks)。这样,一个block内全部cell的特征向量串联起来便获得该blockHOG特征。这些区间是互有重叠的,这就意味着:每个单元格的特征会以不一样的结果屡次出如今最后的特征向量中。咱们将归一化以后的块描述符(向量)就称之为HOG描述符。

        区间有两个主要的几何形状——矩形区间(R-HOG)和环形区间(C-HOG)。R-HOG区间大致上是一些方形的格子,它能够有三个参数来表征:每一个区间中细胞单元的数目、每一个细胞单元中像素点的数目、每一个细胞的直方图通道数目。

       例如:行人检测的最佳参数设置是:3×3细胞/区间、6×6像素/细胞、9个直方图通道。则一块的特征数为:3*3*9

 

5)收集HOG特征

      最后一步就是将检测窗口中全部重叠的块进行HOG特征的收集,并将它们结合成最终的特征向量供分类使用。

    

6)那么一个图像的HOG特征维数是多少呢?

        顺便作个总结:Dalal提出的Hog特征提取的过程:把样本图像分割为若干个像素的单元(cell),把梯度方向平均划分为9个区间(bin),在每一个单元里面对全部像素的梯度方向在各个方向区间进行直方图统计,获得一个9维的特征向量,每相邻的4个单元构成一个块(block),把一个块内的特征向量联起来获得36维的特征向量,用块对样本图像进行扫描,扫描步长为一个单元。最后将全部块的特征串联起来,就获得了人体的特征。例如,对于64*128的图像而言,8*8的像素组成一个cell,每2*2cell组成一个块,由于每一个cell9个特征,因此每一个块内有4*9=36个特征,以8个像素为步长,那么,水平方向将有7个扫描窗口,垂直方向将有15个扫描窗口。也就是说,64*128的图片,总共有36*7*15=3780个特征。

HOG维数,16×16像素组成的block,8x8像素的cell

 

注释:

行人检测HOG+SVM

整体思路:
一、提取正负样本hog特征
二、投入svm分类器训练,获得model
三、由model生成检测子
四、利用检测子检测负样本,获得hardexample
五、提取hardexample的hog特征并结合第一步中的特征一块儿投入训练,获得最终检测子。

深刻研究hog算法原理:
1、hog概述

Histograms of Oriented Gradients,顾名思义,方向梯度直方图,是目标的一种描述的方式,既是描述子。
2、hog提出
hog是05年一位nb的博士提出来的,论文连接 http://wenku.baidu.com/view/676f2351f01dc281e53af0b2.html
3、算法理解
        终于到10月了,终于能够松一口气了,整理一下hog的算法流程。
首先要有一个总体的认识,每个目标都对应一个一维特征向量,这个向量一共有n维,这个n不是凭空瞎猜的,是有理有据, 打个比方,为何opencv自带的hog检测子是3781维的?这个问题在初期确实比较头疼,纠结了好长的时间,不过别着急,
咱们先来看一下opencv里的HOGDescriptor这个结构的构造函数HOGDescriptor(Size winSize,Size blocksize,Size blockStride,Size cellSize,...(后面的参数在这里用不到)),去查一下opencv默认的参数咱们能够看到,winSize(64,128),blockSize(16,16),blockStride(8,8),cellSize(8,8),很显然hog是将一个特征窗口win划分为不少的块block,在每个块里又划分为不少的细胞单元cell(即胞元),hog特征向量既是把这些全部的cell对应的小特征串起来获得一个高维的特征向量,那么这个窗口对应的一维特征向量维数n就等于窗口中的块数 x 块中的胞元数  x 每个胞元对应的特征向量数。
写到这里, 咱们计算一下3781如何获得的,窗口大小64x128,块大小16x16,块步长8x8, 那么窗口中块的数目是((64-16)/8+1)*((128-16)/8+1) = 7*15 =105个块,块大小为16x16,胞元大小为8x8,那么一个块中的胞元cell数目是 (16/8)*(16/8) = 4个胞元,到这里咱们能够看到要求最后须要的维数n,只须要计算每个胞元对应的向量,这个参数在哪呢?别急,咱们把每个胞元投影到9个bin(如何投影?这里卡了很长一段时间,后面会说),那么每个胞元对应的向量就是 9维,每一个bin对应该9维向量的一个数,如今看一下是否是计算窗口维数的三个需求量都知道了,n =  窗口中的块数 x 块中的胞元数  x 每个胞元对应的特征向量数,带入看一下n= 105x4x9 = 3780,这就是这个窗口对应的特征了。有人会说,为何opencv里的getDefaultPeopleDetector()获得的是3781维呢?这是由于另一维是一维偏移,(很崩溃是吧,我也崩溃好久。。。,下一段解释)。
咱们利用hog+svm检测行人,最终的检测方法是最基本的线性判别函数,wx + b = 0,刚才所求的3780维向量其实就是w,而加了一维的b就造成了opencv默认的3781维检测算子,而检测分为train和test两部分,在train期间咱们须要提取一些列训练样本的hog特征使用svm训练最终的目的是为了获得咱们检测的w以及b,在test期间提取待检测目标的hog特征x,带入方程是否是就能进行判别了呢?
**************************************************************************************************
                                           华丽的分割线
写到这里,至少对hog的运做流程有了一个大概的认识,在网上能看到不少的hog计算方法,神马归一化,计算梯度,对每一个胞元进行投影,千篇一概,对刚开始接触的人来讲,看无缺像懂了,但就是不知道怎么用,hog和svm如何配合,并且那些东西对咱们的初期的学期彻底没用,好处就是会用hog了,再回过头去看原理,才有收获,那些资料网上一堆,这里就不多此一举了。
另外值得一提的是在计算胞元特征的时候,须要向各个bin投影,这个投影里面大有文章,师兄毕业论文里就提到了,取名叫‘三维一次线性插值’,若是想深刻了解hog的能够仔细琢磨去。
**************************************************************************************************
                                         继续华丽的分割
下面说一下libsvm和CvSVM的使用,我以为libsvm更好用,不过cvsvm也是基于libsvm2.6(没记错的话)改写的,这两个的区别就是libsvm训练获得的是一个model,而cvsvm是xml文件,在计算最后的wx+b=0中的w向量的时候,对于libsvm直接处理model文件便可,可是对于cvsvm则能够跳过产生xml文件,直接使用cvsvm的对象中的属性便可(这里说的有点模糊,两者选一个便可,关系倒不是很大)
欢迎你们批评指正、交流学习