【学习笔记】计算机视觉与深度学习(1.线性分类器)

学习视频：
鲁鹏-计算机视觉与深度学习

1 图像分类任务

图像分类任务是计算机视觉的核心任务，其目标是根据图像信息中所反映的不同特征，把不同类别的图像区分开来。

比如这只狗，我们的目的是给计算机一些选项，比如说猫、狗、飞机、车，让计算机看到这张图片，并得出“这是一只狗”的结论。给出一些猫、飞机、车的图片，机器也能够分类正确。能做到这一点的机器，可以被称作是一种“分类器”。

2 数据驱动的图像分类方法

分为以下三步：

1. 数据集构建
2. 分类器设计与学习
3. 分类器决策

分类器的学习过程图：

分类器的决策过程图：

3 图像表示方法

3.1 像素表示

根据RGB的原理，图片中的每一个像素的颜色都可以用一个三元组来表示。假设我们要表示的图像为$n\times m$的大小的图像，那么就可以用$3\times n\times m$维向量来表示这张图像。

3.2 全局特征表示(如GIST)

通过提取整张图像中的全部像素体现出来的特征来表征图像，目的是降低维度。全局提取特征的方法适用于景色类图像的分类，细节较多的分类效果不好。

3.3 局部特征表示(如SIFT)

通过提取整张图像中的部分特征点来表示图像，如人脸识别，识别出两个眼睛一个嘴巴一个鼻子作为特征，这时挡住一个眼睛，也能够做一个较好的判断。同样的场景用于全局特征表示方法，由于其使用所有的像素作为特征提取，去掉一个眼睛会对数据造成很大的变化，导致失败率很高。

4 图像分类任务的评价指标

4.1 正确率与错误率

正确率(accuracy) = 分对的样本数 / 全部样本数
错误率(error rate) = 1 - 正确率

4.2 Top1指标与Top5指标

对于同一张图片，分类器给出其得出的五个最有可能正确的答案，并将他们按照正确的可能性从大到小排序。

Top1指标：可能性最大的答案正确，才被认为分类正确。
Top5指标：可能性最大的五个答案中有正确答案，即被认为分类正确。

5 线性分类器

5.1 图像表示

课程中采用的图像表示方法为像素表示法，将图像的每一个像素的颜色用一个三元组$(x_i,y_i,z_i)$表示，其中$0\le x_i,y_i,z_i\le 255$。然后将所有点的$x_i,y_i,z_i$依次排开，假设图像有$n$个点，那么构造一个向量$\mathbf{x}$(本系列笔记中会将代表向量的数学符号像这样加粗)：

$$\mathbf{x}=\left[\begin{array}{c}{x}_0\\ y_0\\ z_0\\ x_1\\ y_1\\ z_1\\ ⋮\end{array}\right]$$

当所选用的图像是$32\times 32$时，所得到的$\mathbf{x}$就是一个$32\times 32\times 3=3072$维的向量。

5.2 简介

线性分类器是一种线性映射，将输入的图像特征映射为类别函数。

设$\mathbf{x}$为输入的$d$维图像向量，$c$为类别个数。对于第$i$个类别的线性分类器，我们有：

$$f_i(\mathbf{x},{\mathbf{w}}_i)={\mathbf{w}}_i^T\mathbf{x}+b_i,i=1\cdots c$$

$f_i$为第$i$个类的线性分类器的打分。${\mathbf{w}}_i={\left[\begin{array}{cccc}{w}_{i1}& w_{i2}& \cdots & w_{id}\end{array}\right]}^T$是一个列向量表示第$i$个类别的权值向量，$b_i$为偏置。

得到的$f_i$是一个常数（不难看出，${\mathbf{w}}_i^T$是$1\times d$的，$\mathbf{x}$是$d\times 1$的）。

对于每一个类别，都有属于自己的$\mathbf{x}$和$b$。通过公式可以看出，求解分数的这个变换是线性的，所以这个分类器叫做线性分类器。

5.3 决策规则

如果$f_i(\mathbf{x})>f_j(\mathbf{x}),\forall j\ne i$，则决策输入图像$\mathbf{x}$输入第$i$类。

其实就是通过线性分类器对这个图像打完分以后，哪一个类别的分数最高，就认为这张图像属于哪一个类别。

如图所示，将$i$个类别的权值向量并到一起可以形成一个权值矩阵$\mathbf{W}$，将每一个类的偏移构成一个便宜向量$\mathbf{b}$，将每一个类最后得到的得分构成一个得分向量$\mathbf{f}$。于是不同于刚刚针对一个类的转移式，我们可以写出一个对于一整个线性分类器的向量转移式：
$$\mathbf{f}(\mathbf{x},\mathbf{W})=\mathbf{Wx}+\mathbf{b}$$

5.4 权值

线性分类器的权值可以看做是一个模板，每一个类别通过学习给定的图像，不断调整其中的数值，使其数值更加贴合对应类别的特征。比如说一个具有10个类别的线性分类器，训练后的模板长这个样子：

其模板若隐若现能看出对应类别的图像，但由于是对很多图片的统合学习，所以会出现${\mathbf{w}}_8$这样子的两个头的马（有的图片马头朝左，有的图片马头朝右导致的）。

输入图像与评估模板的匹配程度越高，分类器输出的分数就越高。

5.5 决策边界

假设我们只用$2$个特征来表征一个图像，通过画分界面将图片分类开来，大概是这样：

分数为$0$的线就是分界面，${\mathbf{w}}_i\mathbf{x}+b_i=0,i=1\cdots c$。
因为只有分数大于$0$，这个图像才有可能为这个类别，假设汽车类分界面，相当于一侧包含所有可能为汽车的图像，而另一侧没有。

我们对线性分类器学习，主要是学习权值矩阵$\mathbf{W}$和偏置向量$\mathbf{b}$，换句话说就是在学习这个分界面。

5.6 损失函数

找到最优的分类模型，还需要损失函数与优化算法的帮忙。我们通过损失函数来判断当前模型的预测情况。具体来说，就是使用测试集对分类器进行测试，观察分类器得出的结果和实际结果的差距，并用损失函数量化表达出来。

损失函数搭建了模型性能和模型参数之间的桥梁，指导模型进行参数优化。具体地说，损失函数是一个函数，用于度量给定分类器的预测值与真实值的不一致程度，其输出通常是一个非负实值。

其输出的非负实值可以作为反馈信号来对分类器参数进行调整，以降低当前示例对应的损失值，提升分类器的分类效果。

损失函数的一般定义

$$L=\frac{1}{N}\sum _i{L}_i(f({\mathbf{x}}_i,\mathbf{W}),y_i)$$

其中$x_i$表示数据集里面第$i$张图片；$f(x_i,W)$为分类器对$x_i$的类别预测；$y_i$是该图片的真实类别；$L_i$为第$i$个样本的预测损失值；$L$为数据集的损失，它是数据集中所有样本损失的平均。

多类支撑向量机损失

$$s_{ij}=f_i({\mathbf{x}}_i,{\mathbf{w}}_j,b_j)={\mathbf{w}}_j^T{\mathbf{x}}_i+b_j$$
这个式子中的部分字母和前面所述不同：
$j$：类别标签，取值范围为 { 1 , 2 , ⋯   , c } \{1,2,\cdots,c\} {1,2,⋯,c}
${\mathbf{w}}_j,b_j$：第$j$个类别的权值向量与偏置参数
${\mathbf{x}}_i$：表示数据集里第$i$个样本
$s_{ij}$：第$i$个样本第$j$类别的预测分数

第$i$个样本的多类支撑向量机损失定义如下：

$$\begin{array}{rl}{L}_i& =\sum _{j\ne y_i}\{\begin{array}{ll}0& if{s}_{y_i}\ge s_{ij}+1\\ s_{ij}-s_{y_i}+1& otherwise\end{array}\\ & =\sum _{j\ne y_i}\max (0,s_{ij}-s_{y_i}+1)\end{array}$$

其中$s_{y_i}$为第$i$个样本真实类别的预测分数。
正确类别的得分比不正确类别的得分高出1分，就没有损失；否则，就会产生损失。

思考：

问 多类支撑向量机损失$L_i$的最大/最小值会是多少？
答 预测非常离谱（分值打到了无穷大的极端情况）时预测值可以达到无穷大；完全预测正确时达到最小值$0$。

问 如果初始化时$\mathbf{w}$和$b$很小，损失$L$会是多少？
答 有$c$个类别时，令$\mathbf{w}$和$b$中的全部权值取$0$，得到$s_{ij}$和$s_{y_i}$为$0$，代入式中得到$L=c-1$。该特性可以用来检测代码编写是否正确。

问 考虑所有类别（包括$j=y_i$），损失$L_i$会有什么变化？
答 最终结果会比原来增加$1$，不会对结果又太大影响，也因此不采用$j=y_i$。

问 在总损失$L$计算时，如果用求和代替平均？
答 根据公式，相当于结果乘上$N$，对结果本身不产生影响。

问 如果使用$L_i=\sum _{j\ne y_i}\max (0,s_j-s_{y_i}+1{)}^2$？
答 大于$1$的损失被放大，小于$1$的损失被缩小，本质上的损失函数会产生一些差异，影响模型优化，可能会出现不同的结果。

问 假设存在一个$\mathbf{W}$使得损失函数$L=0$，这个$\mathbf{W}$是唯一的吗？
答 不是。

5.7 正则项与超参数

通过5.6节的思考中最后一个问题可以得知，可能会存在多种情况，使得损失函数$L=0$，这时我们需要用一种方法来判断哪一种更好，或者说找到一种方法去选择。

做法：在损失函数公式后面加入一个正则项。
$$L=\frac{1}{N}\sum _i{L}_i(f({\mathbf{x}}_i,\mathbf{W}),y_i)+\lambda R(\mathbf{W})$$

式中$\lambda R(\mathbf{W})$为正则项。

其中我们将$\frac{1}{N}\sum _i{L}_i(f({\mathbf{x}}_i,\mathbf{W}),y_i)$称作数据损失，模型预测需要和训练集相匹配；将$\lambda R(\mathbf{W})$称作正则损失，防止模型在训练集上学习得“太好”。

$R(\mathbf{W})$是一个与权值有关，跟图像无关的函数。
$\lambda$是一个超参数，控制着正则损失在总损失中所占的比重。

超参数

超参数是学习过程之前设置的参数，而不是学习中得到的。超参数一般都会对模型的性能有重要的影响。

观察
$$L=\frac{1}{N}\sum _i{L}_i(f({\mathbf{x}}_i,\mathbf{W}),y_i)+\lambda R(\mathbf{W})$$

当$\lambda =0$时，优化结果仅与数据损失相关。
当$\lambda =\infty$时，优化结果与数据损失无关，仅考虑权重损失。此时，系统最优解为$\mathbf{W}=0$。

L2正则项

$$R(\mathbf{W})=\sum _k\sum _l{W}_{k,l}^2$$

举例说明

L2正则损失对大数值权值进行惩罚，喜欢分散数值，鼓励分类器将所有维度的特征都用起来，而不是强烈的依赖其中少数几维特征。

正则项让模型有了偏好！

常用的正则项

L1正则项：$R(\mathbf{W})=\sum _k\sum _l|W_{k,l}|$
L2正则项：$R(\mathbf{W})=\sum _k\sum _l{W}_{k,l}^2$
Elastic net(L1+L2)：$R(\mathbf{W})=\sum _k\sum _l\beta W_{k,l}^2+|W_{k,l}|$

5.8 参数优化

参数优化是机器学习的核心步骤之一，它利用损失函数的输出值作为反馈信号来调整分类器参数，以提升分类器对训练样本的预测性能。

损失函数$L=\frac{1}{N}\sum _{i=1}^N{L}_i+\lambda R(\mathbf{W})$，是一个与参数$\mathbf{W}$有关的函数，优化的目标是找到使损失函数$L$达到最优的那组参数$\mathbf{W}$。

有直接方法：$$\frac{\partial L}{\partial W}=0$$
但由于$L$形式一般很复杂，很难求解，因此提出梯度下降算法。

梯度下降算法

梯度下降算法是一种简单而高效的迭代优化算法。

如图所示，红色点左上角的蒙上眼睛的笑脸是我们当前的位置。假设当前我们在$L$函数的图像上遇到了一个极值点，我们通过一点一点向着往下的坡的方向去移动。于是有两个问题：

1. 往哪儿走？——负梯度方法
2. 要走多远？——步长来决定

梯度下降：利用所有样本计算损失并更新梯度。

while true
		权值的梯度 ← 计算梯度(损失、训练样本、权值)
		权值 ← 权值 - 学习率 * 权值的梯度

数值法

一维变量，函数求导：
$$\frac{dL(w)}{dw}=\underset{h\to 0}{\lim }\frac{L(w+h)-L(w)}{h}$$

解析法

示例：损失函数$L(w)=w^2$，求$w=1$点处的梯度。
$\nabla L(w)=2w$
${\nabla }_{w=1}L(w)=2$
解析法精确、快速，但是导数函数推导易错。

一般求解梯度使用解析法，通过数值法验证（梯度检查）。

如何计算多类支撑向量机损失的导数函数？

$$L_i=\sum _{j\ne y_i}\max (0,s_{ij}-s_{y_i}+1)$$

将$s_{ij}={\mathbf{w}}_j^T{\mathbf{x}}_i+b_j$代入上式，得到：

$$L_i=\sum _{j\ne y_i}\max (0,{\mathbf{w}}_j^T{\mathbf{x}}_i+b_j-{\mathbf{w}}_{y_i}^T{\mathbf{x}}_i+b_{y_i}+1)$$

现尝试通过解析法求解，我们先以简单例子入手：假设$y=\max (0,w^2-1)$，那么有
$$\frac{\partial y}{\partial w}=\{\begin{array}{ll}2w& if{w}^2-1\ge 0\\ 0& otherwise\end{array}$$

同理，对于我们要求的$\frac{\partial L_i}{\partial {\mathbf{w}}_j}$，也可以分段：
当${\mathbf{w}}_j^T{\mathbf{x}}_i+b_j-{\mathbf{w}}_{y_i}^T{\mathbf{x}}_i+b_{y_i}+1<0$时，有$\frac{\partial L_i}{\partial {\mathbf{w}}_j}=0$。
当${\mathbf{w}}_j^T{\mathbf{x}}_i+b_j-{\mathbf{w}}_{y_i}^T{\mathbf{x}}_i+b_{y_i}+1\ge 0$时，有$\frac{\partial L_i}{\partial {\mathbf{w}}_j}={\mathbf{x}}_i$。这个导数涉及到了矩阵求导，推导过程如下：

观察$L_i$式中只有${\mathbf{w}}_j^T{\mathbf{x}}_i$这一项与${\mathbf{w}}_j$有关，其他的都可以忽略不计，于是我们可以把$L_i$的式子（仅在求导时）看作：
$$L_i={\mathbf{w}}_j^T{\mathbf{x}}_i=x_{i1}{w}_{j1}+x_{i2}{w}_{j2}+\cdots +x_{id}{w}_{jd}$$

根据矩阵求导，我们有
$$\frac{\partial L_i}{\partial {\mathbf{w}}_j}=\left[\begin{array}{c}\frac{\partial L_i}{w_{j1}}\\ \frac{\partial L_i}{w_{j2}}\\ ⋮\\ \frac{\partial L_i}{w_{jd}}\end{array}\right]=\left[\begin{array}{c}{x}_{i1}\\ x_{i2}\\ ⋮\\ x_{id}\end{array}\right]={\mathbf{x}}_i$$

所以我们有：

$$\frac{\partial L_i}{\partial {\mathbf{w}}_j}=\{\begin{array}{ll}{\mathbf{x}}_i& if{\mathbf{w}}_j^T{\mathbf{x}}_i+b_j-{\mathbf{w}}_{y_i}^T{\mathbf{x}}_i+b_{y_i}+1\ge 0\\ 0& otherwise\end{array}$$

而$\frac{\partial L_i}{\partial b_j}$较好求，不多说明：

$$\frac{\partial L_i}{\partial b_j}=\{\begin{array}{ll}1& if{\mathbf{w}}_j^T{\mathbf{x}}_i+b_j-{\mathbf{w}}_{y_i}^T{\mathbf{x}}_i+b_{y_i}+1\ge 0\\ 0& otherwise\end{array}$$

利用这些已求导数，我们可以求出$L(\mathbf{W})$的梯度。

$$L(\mathbf{W})=\frac{1}{N}\sum _{i=1}^N{L}_i(x_i,y_i,\mathbf{W})+\lambda R(\mathbf{W})$$

$${\nabla }_{\mathbf{W}}L(\mathbf{W})=\frac{1}{N}\sum _{i=1}^N{\nabla }_{\mathbf{W}}{L}_i(x_i,y_i,\mathbf{W})+\lambda {\nabla }_{\mathbf{W}}R(\mathbf{W})$$

while true
		权值的梯度 ← 计算梯度(损失、训练样本、权值)
		权值 ← 权值 - 学习率 * 权值的梯度

不难发现，当$N$很大时，权值的梯度所需的计算量很大。

随机梯度下降算法

每次随机选择一个样本$x_i$，计算损失并更新梯度。

$$L(\mathbf{W})=L_i(x_i,y_i,\mathbf{W})+\lambda R(\mathbf{W})$$

$${\nabla }_{\mathbf{W}}L(\mathbf{W})={\nabla }_{\mathbf{W}}{L}_i(x_i,y_i,\mathbf{W})+\lambda {\nabla }_{\mathbf{W}}R(\mathbf{W})$$

while true
		数据 ← 从训练数据采样(训练数据, 1)
		权值的梯度 ← 计算梯度(损失, 训练样本, 权值)
		权值 ← 权值 - 学习率 * 权值的梯度

但是单个样本的训练可能会带来很多噪声，不是每次迭代都向着整体最优化的方向。

小批量梯度下降算法

每次随机选取$m$个（批量的大小，也是一个超参数）个样本，计算损失并更新梯度。

$$L(\mathbf{W})=\frac{1}{m}\sum _{i=1}^m{L}_i(x_i,y_i,\mathbf{W})+\lambda R(\mathbf{W})$$

$${\nabla }_{\mathbf{W}}L(\mathbf{W})=\frac{1}{m}\sum _{i=1}^N{\nabla }_{\mathbf{W}}{L}_i(x_i,y_i,\mathbf{W})+\lambda {\nabla }_{\mathbf{W}}R(\mathbf{W})$$

一些名词的介绍（会实现的代码中）
iteration：表示1次迭代，每次迭代更新1次网络结构的参数；
batch_size：1次迭代所使用的样本量；
epoch：1个epoch表示过了1遍训练集中的所有样本。

通常使用2的幂数作为批量大小，比如说32或64或128个样本。

while true
		小批量数据 ← 从训练数据采样(训练数据, 批量大小)
		权值的梯度 ← 计算梯度(损失, 小批量数据, 权值)
		权值 ← 权值 - 学习率 * 权值的梯度

6 数据集

6.1 数据集划分

训练集→训练模型，寻找最优的分类器；
测试集→评估模型，评测泛化能力。

问 如果模型含有超参数（比如正则化强度），如何找到泛化能力最好的超参数？
答 再设立一个验证集，用于选择超参数。

训练集用于给定的超参数时分类器参数的学习；
验证集用于选择超参数；
测试集评估泛化能力。

（这样可以避免训练过程中看到测试集，使得出现过拟合现象）

问 如果数据很少，那么可能验证集包含的样本就很少，从而无法在统计上代表数据。
这个问题很容易发现，如果在划分数据前进行不同的随机打乱，最终得到的模型性能差别很大，那么就存在这个问题。

K折交叉验证

带打乱顺序的K折交叉验证

6.2 数据预处理

数据能否直接使用？有哪些处理方式？

通常不太会直接使用原始数据进行学习。

假设原始数据中的数据分别为$x_i$，均值为$\bar{x}$，最小值为$\min (x)$，最大值为$\max (x)$，标准差为$\sigma (x)$。

去均值的方法：$x_i\leftarrow x_i-\bar{x}$。
归一化的方法有两种：

1. min-max归一化：$x_i\leftarrow \frac{x_i-\min (x)}{\max (x)-\min (x)}$
2. mean归一化：$x_i\leftarrow \frac{x_i-\bar{x}}{\max (x)-\min (x)}$

标准化的方法：$x_i\leftarrow \frac{x_i-\bar{x}}{\sigma (x)}$

事实上，归一化和标准化都是把数据图像中心平移到原点。但是归一化没有改变数据分布的形状，而标准化使样本数据的分布近似为某种分布（通常为正态分布）。

原始数据中可以看出，$x$增大时，$y$也有总体增大的趋势，所以我们可以认为此时$x$和$y$是相关的。但学习的过程中，我们不希望把$x$和$y$两种东西放在一起去考虑，这时我们采取去相关的操作。去相关后，可以看出$y$的变化方式与$x$没有关系了。如此一来，我们可以把$x$和$y$分开单独考虑，实现“降维”的效果。

白化为去相关后进行标准化。

神经网络中用到更多的是去均值和标准化、归一化。