文章同步发于知乎。最新版以
github
为主。如果看完文章有所收获,一定要先点赞后收藏。毕竟,赠人玫瑰,手有余香。
MobileNet
论文的主要贡献在于提出了一种深度可分离卷积架构(DW+PW 卷积),先通过理论证明这种架构比常规的卷积计算成本(Mult-Adds
)更小,然后通过分类、检测等多种实验证明模型的有效性。
一个大小为 2
维卷积核),沿着 feature map
的左上角移动到右下角,过滤器每移动一次,将过滤器参数矩阵和对应特征图 feature map
的数量(通道数)为 shape
为 $ (c_1, h_1, w_1)$ 的滤波器( 3
维卷积核),将每个输入 featue map 对应输出像素点位置计算和的值相加,即得到输出 feature map 对应像素点的值。又因为输出 feature map
的数量为
2D
卷积计算过程动态图如下,通过这张图能够更直观理解卷积核如何执行滑窗操作,又如何相加并输出 feature map
,动态图来源 这里。
Input Shape : (3, 7, 7) — Output Shape : (2, 3, 3) — K : (3, 3) — P : (1, 1) — S : (2, 2) — D : (2, 2) — G : 1
上述动态图侧重于滤波器的划窗计算过程,而下图则侧重于解释标准卷积层五个概念的关系,图片来源 这里。
- 输入 Input Channel
- 卷积核组 WeightsBias
- 过滤器 Filter
- 卷积核 kernel
- 输出 Feature Map
在此例中,输入是三维数据
Group Convolution
分组卷积,最早见于 AlexNet
。常规卷积与分组卷积的输入 feature map 与输出 feature map 的连接方式如下图所示,图片来自CondenseNet。
分组卷积的定义:对输入 feature map
进行分组,然后分组分别进行卷积。假设输入 feature map 的尺寸为 shape
为 concatenate
) 得到最终的得到最终尺寸为
分组卷积的意义:分组卷积是现在网络结构设计的核心,它通过通道之间的稀疏连接(也就是只和同一个组内的特征连接)来降低计算复杂度。一方面,它允许我们使用更多的通道数来增加网络容量进而提升准确率,但另一方面随着通道数的增多也对带来更多的
从以上公式可以得出
分组卷积的深入理解:对于
当分组卷积的分组数量 = 输入 feature map 数量 = 输出 feature map 数量,即 DW
卷积的卷积核权重尺寸为
常规卷积的卷积核权重 shape 都为(
C_out, C_in, kernel_height, kernel_width
),分组卷积的卷积核权重shape
为(C_out, C_in/g, kernel_height, kernel_width
),DW
卷积的卷积核权重shape
为(C_in, 1, kernel_height, kernel_width
)。
对于分组卷积,存在一定的限制:卷积层的输出通道数必须是分组数的整数倍,即
假设,$c_2$ 是输出通道数,$c_1$ 是输入通道数,$g$ 是分组数。
和标准卷积每个 Filter 内都有一个或多个卷积核 Kernal,对应每个输入通道(Input Channel)的特性不同,分组卷积和 DW 卷积的特点如下:
- 分组卷积:分组卷积是将输入通道分成若干组,每组的滤波器只与其同组的输入 feature map 进行卷积,最终将每组的输出通道拼接在一起得到最终输出。
- DW 卷积:每个 Filter 内只有一个卷积核 Kernal,对应每个输入通道(Input Channel),即对于每个输入通道分别使用一个固定大小的卷积核进行卷积操作。
分组卷积的极致是分组数数等于输入通道数,这其实就是 DW
卷积,可视化如下:
另外,对于 pytorch
的卷积层 api 是同时支持普通卷积、分组卷积/DW 卷积的。但值得注意的是,对于分组卷积,卷积层的输出通道数必须是分组数的整数倍,否则代码会报错!
import torch
input = torch.randn([20, 10, 224, 224]) # input_channels = 10
conv3x3 = torch.nn.Conv2d(in_channels = 10, output_channels = 5, kernel_size=3, groups=5)
output = conv3x3(input)
print(conv3x3.weight.shape)
print(output.shape)
如果将 groups=5
改为 groups=6
或者将 output_channels = 5
改为 6
,则会报错:
ValueError: in_channels must be divisible by groups
ValueError: out_channels must be divisible by groups
深度可分离卷积(depthwise separable convolutions)的提出最早来源于 Xception
论文,Xception 的论文中提到,对于卷积来说,卷积核可以看做一个三维的滤波器:通道维+空间维(Feature Map 的宽和高),常规的卷积操作其实就是实现通道相关性和空间相关性的联合映射。Inception 模块的背后存在这样的一种假设:卷积层通道间的相关性和空间相关性是可以退耦合(完全可分)的,将它们分开映射,能达到更好的效果(the fundamental hypothesis behind Inception is that cross-channel correlations and spatial correlations are sufficiently decoupled that it is preferable not to map them jointly.)。
引入深度可分离卷积的 Inception,称之为 Xception,其作为 Inception v3 的改进版,在 ImageNet 和 JFT 数据集上有一定的性能提升,但是参数量和速度并没有太大的变化,因为 Xception 的目的也不在于模型的压缩。深度可分离卷积的 Inception 模块如图 Figure 4 所示。
Figure 4 中的“极限” Inception 模块与本文的主角-深度可分离卷积模块相似,区别在于:深度可分离卷积先进行 channel-wise
的空间卷积,再进行
MobileNets
是谷歌 2017 年提出的一种高效的移动端轻量化网络,其核心是深度可分离卷积(depthwise separable convolutions
),深度可分离卷积的核心思想是将一个完整的卷积运算分解为两步进行,分别为 Depthwise Convolution(DW
卷积) 与 Pointwise Convolution(PW
卷积)。深度可分离卷积的计算步骤和滤波器尺寸如下所示。
图片来源论文 MobileNets: Efficient Convolutional Neural Networks for Mobile Vision Applications。
注意本文 DW 和 PW 卷积计算量的计算与论文有所区别,本文的输出 Feature map 大小是
$D_G \times D_G$ , 论文公式是$D_F \times D_F$。
不同于常规卷积操作, Depthwise Convolution 的一个卷积核只负责一个通道,一个通道只能被一个卷积核卷积(不同的通道采用不同的卷积核卷积),也就是输入通道、输出通道和分组数相同的特殊分组卷积,因此 Depthwise(DW
)卷积不会改变输入特征图的通道数目。深度可分离卷积的 DW
卷积步骤如下图:
DW
卷积的计算量
上述 Depthwise 卷积的问题在于它让每个卷积核单独对一个通道进行计算,但是各个通道的信息没有达到交换,从而在网络后续信息流动中会损失通道之间的信息,因此论文中就加入了 Pointwise 卷积操作,来进一步融合通道之间的信息。PW 卷积是一种特殊的常规卷积,卷积核的尺寸为 PW
卷积的过程如下图:
假设输入特征图大小为 PW
卷积的计算量
综上:Depthwise
和 Pointwise
卷积这两部分的计算量相加为
可以看出 Depthwise + Pointwise
卷积的计算量相较于标准卷积近乎减少了 backbone
都采用了这种卷积结构,再加上模型蒸馏,剪枝,能让移动端更高效的推理。
Block
结构如下图所示:
左边是带 bn
和 relu
的标准卷积层,右边是带 bn 和 relu 的深度可分离卷积层。
Block
网络的 pytorch 代码如下:
class MobilnetV1Block(nn.Module):
"""Depthwise conv + Pointwise conv"""
def __init__(self, in_channels, out_channels, stride=1):
super(MobilnetV1Block, self).__init__()
# dw conv kernel shape is (in_channels, 1, ksize, ksize)
self.dw = nn.Conv2d(in_channels, in_channels, kernel_size=3,stride=stride,padding=1, groups=in_channels, bias=False)
self.bn1 = nn.BatchNorm2d(in_channels)
self.pw = nn.Conv2d(in_channels, out_channels, kernel_size=1, stride=1, padding=0, bias=False)
self.bn2 = nn.BatchNorm2d(out_channels)
def forward(self, x):
out1 = F.relu(self.bn1(self.dw(x)))
out2 = F.relu(self.bn2(self.pw(out1)))
return out2
MobileNet v1
的 pytorch
模型导出为 onnx
模型后,深度可分离卷积 block
结构图如下图所示。
仅用 MobileNets 的 Mult-Adds
(乘加操作)次数更少来定义高性能网络是不够的,确保这些操作能够有效实施也很重要。例如非结构化稀疏矩阵运算(unstructured sparse matrix operations)通常并不会比密集矩阵运算(dense matrix operations)快,除非是非常高的稀疏度。
这句话是不是和
shufflenet v2
的观点一致,即不能仅仅以 FLOPs 计算量来表现网络的运行速度,除非是同一种网络架构。
MobileNet 模型结构将几乎所有计算都放入密集的 1×1 卷积中(dense 1 × 1 convolutions),卷积计算可以通过高度优化的通用矩阵乘法(GEMM
)函数来实现。 卷积通常由 GEMM
实现,但需要在内存中进行名为 im2col
的初始重新排序,然后才映射到 GEMM
。 caffe 框架就是使用这种方法实现卷积计算。 1×1
卷积不需要在内存中进行重新排序,可以直接使用 GEMM
(最优化的数值线性代数算法之一)来实现。
如表 2 所示,MobileNet 模型的 1x1
卷积占据了 95%
的计算量和 75%
的参数,剩下的参数几乎都在全连接层中, 3x3 的 DW 卷积核常规卷积占据了很少的计算量(Mult-Adds)和参数。
尽管基本的 MobileNet
体系结构已经很小且网络延迟 latency
很低,但很多情况下特定用例或应用可能要求模型变得更小,更快。为了构建这些更小且计算成本更低的模型,我们引入了一个非常简单的参数 width 乘数
。宽度乘数
其中 1、0.75、0.5
和 0.25
。$\alpha = 1$ 是基准 MobileNet 模型,$\alpha < 1$ 是缩小版的 MobileNets
。宽度乘数的作用是将计算量和参数数量大约减少 latency
和模型大小之间的权衡。 它用于定义新的精简结构,需要从头开始进行训练模型。基准 MobileNet
模型的整体结构定义如表 1 所示。
减少模型计算成本的的第二个超参数(hyper-parameter)是分辨率因子 224、192、160
或 128
。$\rho = 1$ 时是基准(baseline
) MobilNet,$\rho < 1$ 时缩小版 MobileNets
。分辨率乘数的作用是将计算量减少
- 整个网络不算平均池化层与
softmax
层,且将DW
卷积和PW
卷积计为单独的一层,则MobileNet
有28
层网络。+ 在整个网络结构中步长为2的卷积较有特点,卷积的同时充当下采样的功能; - 第一层之后的
26
层都为深度可分离卷积的重复卷积操作,分为4
个卷积stage
; - 每一个卷积层(含常规卷积、深度卷积、逐点卷积)之后都紧跟着批规范化和
ReLU
激活函数; - 最后一层全连接层不使用激活函数。
作者分别进行了 Stanford Dogs dataset 数据集上的细粒度识别、大规模地理分类、人脸属性分类、COCO 数据集上目标检测的实验,来证明与 Inception V3
、GoogleNet
、VGG16
等 backbone
相比,MobilNets
模型可以在计算量(Mult-Adds
)数 10 被下降的情况下,但是精度却几乎不变。
论文提出了一种基于深度可分离卷积的新模型架构,称为 MobileNets
。 在相关工作章节中,作者首先调查了一些让模型更有效的重要设计原则,然后,演示了如何通过宽度乘数和分辨率乘数来构建更小,更快的 MobileNet,通过权衡合理的精度以减少模型大小和延迟。 然后,我们将不同的 MobileNets
与流行的模型进行了比较,这些模型展示了出色的尺寸,速度和准确性特性。 最后,论文演示了 MobileNet 在应用于各种任务时的有效性。
自己复现的基准 MobileNet v1 代模型 pytorch 代码如下:
import torch
import torch.nn as nn
import torch.nn.functional as F
import torchvision.models as models
from torch import flatten
class MobilnetV1Block(nn.Module):
"""Depthwise conv + Pointwise conv"""
def __init__(self, in_channels, out_channels, stride=1):
super(MobilnetV1Block, self).__init__()
# dw conv kernel shape is (in_channels, 1, ksize, ksize)
self.dw = nn.Sequential(
nn.Conv2d(in_channels, 64, kernel_size=3,
stride=stride, padding=1, groups=4, bias=False),
nn.BatchNorm2d(in_channels),
nn.ReLU(inplace=True)
)
# print(self.dw[0].weight.shape) # print dw conv kernel shape
self.pw = nn.Sequential(
nn.Conv2d(in_channels, out_channels, kernel_size=1,
stride=1, padding=0, bias=False),
nn.BatchNorm2d(out_channels),
nn.ReLU(inplace=True)
)
def forward(self, x):
x = self.dw(x)
x = self.pw(x)
return x
def convbn_relu(in_channels, out_channels, stride=2):
return nn.Sequential(nn.Conv2d(in_channels, out_channels, kernel_size=3, stride=stride,
padding=1, bias=False),
nn.BatchNorm2d(out_channels),
nn.ReLU(inplace=True))
class MobileNetV1(nn.Module):
# (32, 64, 1) means MobilnetV1Block in_channnels is 32, out_channels is 64, no change in map size.
stage_cfg = [(32, 64, 1),
(64, 128, 2), (128, 128, 1), # stage1 conv
(128, 256, 2), (256, 256, 1), # stage2 conv
(256, 512, 2), (512, 512, 1), (512, 512, 1), (512, 512, 1), (512, 512, 1), (512, 512, 1), # stage3 conv
(512, 1024, 2), (1024, 1024, 1) # stage4 conv
]
def __init__(self, num_classes=1000):
super(MobileNetV1, self).__init__()
self.first_conv = convbn_relu(3, 32, 2) # Input image size reduced by half
self.stage_layers = self._make_layers(in_channels=32)
self.linear = nn.Linear(1024, num_classes) # 全连接层
def _make_layers(self, in_channels):
layers = []
for x in self.stage_cfg:
in_channels = x[0]
out_channels = x[1]
stride = x[2]
layers.append(MobilnetV1Block(in_channels, out_channels, stride))
in_channels = out_channels
return nn.Sequential(*layers)
def forward(self, x):
"""Feature map shape(h、w) is 224 -> 112 -> 56 -> 28 -> 14 -> 7 -> 1"""
x = self.first_conv(x)
x = self.stage_layers(x)
x = F.avg_pool2d(x, 7) # x shape is 7*7
x = flatten(x, 1) # x = x.view(x.size(0), -1)
x = self.linear(x)
return x
if __name__ == "__main__":
model = MobileNetV1()
model.eval() # set the model to inference mode
input_data = torch.rand(1, 3, 224, 224)
outputs = model(input_data)
print("Model output size is", outputs.size())
程序运行结果如下:
Model output size is torch.Size([1, 1000])
在降低 FLOPs
计算量上,MobileNet
的网络架构设计确实很好,但是 MobileNet
模型在 GPU
、DSP
和 TPU
硬件上却不一定性能好,原因是不同硬件进行运算时的行为不同,从而导致了 FLOPs
少不等于 latency
低的问题。
如果要实际解释 TPU
与 DSP
的运作原理,可能有点麻烦,可以参考下图,从结果直观地理解他们行为上的差异。考虑一个简单的 convolution
,在 CPU
上 latency
随着 input
与 output
的channel
上升正相关的增加。然而在 DSP
上却是阶梯型,甚至在更高的 channel
数下存在特别低latency
的甜蜜点。
在一定的程度上,网络越深越宽,性能越好。宽度,即通道(channel
)的数量,网络深度,即 layer
的层数,如 resnet18
有 18
个卷积层。注意我们这里说的和宽度学习一类的模型没有关系,而是特指深度卷积神经网络的(通道)宽度。
- 网络深度的意义:
CNN
的网络层能够对输入图像数据进行逐层抽象,比如第一层学习到了图像边缘特征,第二层学习到了简单形状特征,第三层学习到了目标形状的特征,网络深度增加也提高了模型的抽象能力。 - 网络宽度的意义:网络的宽度(通道数)代表了滤波器(
3
维)的数量,滤波器越多,对目标特征的提取能力越强,即让每一层网络学习到更加丰富的特征,比如不同方向、不同频率的纹理特征等。
FLOPs
低不等于latency
低,尤其是在有加速功能的硬体 (GPU
、DSP
与TPU
)上不成立。MobileNet conv block
(depthwise separable convolution
)在有加速功能的硬件(专用硬件设计-NPU
芯片)上比较没有效率。- 低
channel
数的情况下 (如网路的前几层),在有加速功能的硬件使用普通convolution
通常会比separable convolution
有效率。 - 在大多数的硬件上,
channel
数为8
的倍数比较有利计算。 DSP
与TPU
上,一般我们需要运算为uint8
形式,quantization
(低精度量化)是常见的技巧。DSP
与TPU
上,h-Swish
与squeeze-and-excitation
效果不明显 (因为硬体设计与uint8
压缩的关系)。DSP
与TPU
上,5x5
convolution
比较没效率。
- Group Convolution分组卷积,以及Depthwise Convolution和Global Depthwise Convolution
- 理解分组卷积和深度可分离卷积如何降低参数量
- 深度可分离卷积(Xception 与 MobileNet 的点滴)
- MobileNetV1代码实现
- Depthwise卷积与Pointwise卷积
- 【CNN结构设计】深入理解深度可分离卷积
- FLOPs与模型推理速度
- MobileDets: FLOPs不等于Latency,考量不同硬体的高效架构
- MobileNets: Efficient Convolutional Neural Networks for Mobile Vision Applications