5.5. 卷积神经网络(LeNet)

    • 图像在同一列邻近的像素在这个向量中可能相距较远。它们构成的模式可能难以被模型识别。
    • 对于大尺寸的输入图像,使用全连接层容易造成模型过大。假设输入是高和宽均为 像素的彩色照片(含3个通道)。即使全连接层输出个数仍是256,该层权重参数的形状是 5.5. 卷积神经网络(LeNet) - 图2 :它占用了大约3GB的内存或显存。这带来过复杂的模型和过高的存储开销。 卷积层尝试解决这两个问题。一方面,卷积层保留输入形状,使图像的像素在高和宽两个方向上的相关性均可能被有效识别;另一方面,卷积层通过滑动窗口将同一卷积核与不同位置的输入重复计算,从而避免参数尺寸过大。

    卷积神经网络就是含卷积层的网络。本节里我们将介绍一个早期用来识别手写数字图像的卷积神经网络:LeNet[1]。这个名字来源于LeNet论文的第一作者YannLeCun。LeNet展示了通过梯度下降训练卷积神经网络可以达到手写数字识别在当时最先进的结果。这个奠基性的工作第一次将卷积神经网络推上舞台,为世人所知。

    LeNet分为卷积层块和全连接层块两个部分。下面我们分别介绍这两个模块。

    卷积层块里的基本单位是卷积层后接最大池化层:卷积层用来识别图像里的空间模式,如线条和物体局部,之后的最大池化层则用来降低卷积层对位置的敏感性。卷积层块由两个这样的基本单位重复堆叠构成。在卷积层块中,每个卷积层都使用

    的窗口,并在输出上使用sigmoid激活函数。第一个卷积层输出通道数为6,第二个卷积层输出通道数则增加到16。这是因为第二个卷积层比第一个卷积层的输入的高和宽要小,所以增加输出通道使两个卷积层的参数尺寸类似。卷积层块的两个最大池化层的窗口形状均为 5.5. 卷积神经网络(LeNet) - 图4 ,且步幅为2。由于池化窗口与步幅形状相同,池化窗口在输入上每次滑动所覆盖的区域互不重叠。

    下面我们通过类来实现LeNet模型。

    1. import d2lzh as d2l
    2. import mxnet as mx
    3. from mxnet import autograd, gluon, init, nd
    4. from mxnet.gluon import loss as gloss, nn
    5. import time
    6.  
    7. net = nn.Sequential()
    8. net.add(nn.Conv2D(channels=6, kernel_size=5, activation='sigmoid'),
    9. nn.MaxPool2D(pool_size=2, strides=2),
    10. nn.Conv2D(channels=16, kernel_size=5, activation='sigmoid'),
    11. nn.MaxPool2D(pool_size=2, strides=2),
    12. # Dense会默认将(批量大小, 通道, 高, 宽)形状的输入转换成
    13. # (批量大小, 通道 * 高 * 宽)形状的输入
    14. nn.Dense(120, activation='sigmoid'),
    15. nn.Dense(84, activation='sigmoid'),
    16. nn.Dense(10))

    接下来我们构造一个高和宽均为28的单通道数据样本,并逐层进行前向计算来查看每个层的输出形状。

    1. In [2]:
    1. X = nd.random.uniform(shape=(1, 1, 28, 28))
    2. net.initialize()
    3. for layer in net:
    4. X = layer(X)
    5. print(layer.name, 'output shape:\t', X.shape)
    1. pool0 output shape: (1, 6, 12, 12)
    2. conv1 output shape: (1, 16, 8, 8)
    3. pool1 output shape: (1, 16, 4, 4)
    4. dense0 output shape: (1, 120)
    5. dense1 output shape: (1, 84)
    6. dense2 output shape: (1, 10)

    可以看到,在卷积层块中输入的高和宽在逐层减小。卷积层由于使用高和宽均为5的卷积核,从而将高和宽分别减小4,而池化层则将高和宽减半,但通道数则从1增加到16。全连接层则逐层减少输出个数,直到变成图像的类别数10。

    下面我们来实验LeNet模型。实验中,我们仍然使用Fashion-MNIST作为训练数据集。

    1. In [3]:
    1. In [4]:
    1. def try_gpu(): # 本函数已保存在d2lzh包中方便以后使用
    2. try:
    3. ctx = mx.gpu()
    4. _ = nd.zeros((1,), ctx=ctx)
    5. except mx.base.MXNetError:
    6. ctx = mx.cpu()
    7. return ctx
    8.  
    9. ctx = try_gpu()
    10. ctx
    1. Out[4]:
    1. gpu(0)

    相应地,我们对“softmax回归的从零开始实现”一节中描述的evaluate_accuracy函数略作修改。由于数据刚开始存在CPU使用的内存上,当ctx变量代表GPU及相应的显存时,我们通过一节中介绍的函数将数据复制到显存上,例如gpu(0)

    1. In [5]:

    我们同样对“softmax回归的从零开始实现”一节中定义的train_ch3函数略作修改,确保计算使用的数据和模型同在内存或显存上。

    1. In [6]:
    1. # 本函数已保存在d2lzh包中方便以后使用
    2. def train_ch5(net, train_iter, test_iter, batch_size, trainer, ctx,
    3. num_epochs):
    4. print('training on', ctx)
    5. loss = gloss.SoftmaxCrossEntropyLoss()
    6. for epoch in range(num_epochs):
    7. train_l_sum, train_acc_sum, n, start = 0.0, 0.0, 0, time.time()
    8. X, y = X.as_in_context(ctx), y.as_in_context(ctx)
    9. with autograd.record():
    10. y_hat = net(X)
    11. l = loss(y_hat, y).sum()
    12. l.backward()
    13. trainer.step(batch_size)
    14. y = y.astype('float32')
    15. train_l_sum += l.asscalar()
    16. train_acc_sum += (y_hat.argmax(axis=1) == y).sum().asscalar()
    17. n += y.size
    18. test_acc = evaluate_accuracy(test_iter, net, ctx)
    19. print('epoch %d, loss %.4f, train acc %.3f, test acc %.3f, '
    20. 'time %.1f sec'
    21. % (epoch + 1, train_l_sum / n, train_acc_sum / n, test_acc,
    22. time.time() - start))

    我们重新将模型参数初始化到设备变量之上,并使用Xavier随机初始化。损失函数和训练算法则依然使用交叉熵损失函数和小批量随机梯度下降。

    1. In [7]:
    1. lr, num_epochs = 0.9, 5
    2. net.initialize(force_reinit=True, ctx=ctx, init=init.Xavier())
    3. trainer = gluon.Trainer(net.collect_params(), 'sgd', {'learning_rate': lr})
    4. train_ch5(net, train_iter, test_iter, batch_size, trainer, ctx, num_epochs)
    1. training on gpu(0)
    2. epoch 1, loss 2.3197, train acc 0.102, test acc 0.100, time 1.9 sec
    3. epoch 2, loss 2.1476, train acc 0.173, test acc 0.525, time 1.6 sec
    4. epoch 3, loss 1.0506, train acc 0.578, test acc 0.657, time 1.6 sec
    5. epoch 4, loss 0.7780, train acc 0.698, test acc 0.728, time 1.6 sec
    6. epoch 5, loss 0.6767, train acc 0.732, test acc 0.758, time 1.6 sec
    • 卷积神经网络就是含卷积层的网络。
    • LeNet交替使用卷积层和最大池化层后接全连接层来进行图像分类。
    • 尝试基于LeNet构造更复杂的网络来提高分类准确率。例如,调整卷积窗口大小、输出通道数、激活函数和全连接层输出个数。在优化方面,可以尝试使用不同的学习率、初始化方法以及增加迭代周期。

    [1] LeCun, Y., Bottou, L., Bengio, Y., & Haffner, P. (1998).Gradient-based learning applied to document recognition. Proceedings ofthe IEEE, 86(11), 2278-2324.