CIFAR-10数据集应用：快速入门数据增强方法Mixup，显著提升图像识别准确度

格物钛公开数据集提供海量优质公开数据集搜索、数据托管、一站式搜索，与全球 AI 开发者共创公开数据集社区。本文介绍的是以格物钛公开数据集平台中的 CIFAR-10 数据集为基础，通过数据增强方法 Mixup，显著提升图像识别准确度。

关于作者：Ta-Ying Cheng，牛津大学博士研究生，Medium 技术博主，多篇文章均被平台官方刊物 Towards Data Science 收录（翻译：颂贤）。

深度学习蓬勃发展的这几年来，图像分类一直是最为火热的领域之一。传统上的图像识别严重依赖像是扩张/侵蚀或者是频域变换这样的处理方法，但特征提取的困难性限制了这些方法的进步空间。现如今的神经网络则显著提高了图像识别的准确率，因为神经网络能够寻找输入图像和输出标签之间的关系，并以此不断地调整它的识别策略。

然而，神经网络往往需要大量的数据进行训练，而优质的训练数据并不是唾手可得的。因此现在许多人都在研究如何能够实现所谓的数据增强（Data augmentation），即在一个已有的小数据集中凭空增加数据量，来达到以一敌百的效果。本文就将带大家认识一种简单而有效的数据增强策略 Mixup，并介绍直接在 PyTorch 中实现 Mixup 的方法。

为什么需要数据增强？

神经网络架构内的参数是根据给定的数据进行训练和更新的。但由于训练数据只覆盖了某一部分可能数据的分布情况，网络很可能就会在分布的“能见”部分过度拟合。因此，我们拥有的训练数据越多，理论上就越能覆盖整个分布的情况（这也正是为什么以数据为中心的 AI（data-centric AI）非常重要）。当然，在数据量有限的情况下，我们也并不是没有办法。通过数据增强，我们就可以尝试通过微调原有数据的方式产生新数据，并将其作为“新”样本送入网络进行训练。

什么是Mixup？

▲ 图1：Mixup的简易演示图

假设我们现在要做的事情是给猫和狗的图片做分类，并且我们已经有了一组标注好了是猫是狗的数据（例如[1, 0] -> 狗, [0, 1] -> 猫），那么 Mixup 简单来说就是将两张图像及其标签平均化为一个新数据。

具体而言，我们可以用数学公式写出 Mixup 的概念：

其中， 和 分别是混合 （标签为 ）和 （标签为 ）后的图像和标签，而 λ 则是从给定的贝塔分布中取得的随机数。

由此，Mixup 能够为我们提供不同数据类别之间的连续数据样本，并因此直接扩大了给定训练集的分布，从而使网络在测试阶段更加强大。

Mixup的万用性

Mixup 其实只是一种数据增强方法，它和任何用于分类的网络架构都是正交的。也就是说，我们可以在任何要进行分类任务的网络中对相应的数据集使用 Mixup 方法。Mixup 的提出者张宏毅等人基于其最初发表的论文《mixup: Beyond Empirical Risk Minimization》对多个数据集和架构进行了实验，发现了 Mixup 在神经网络之外的应用中也能体现其强大能力。

计算环境

库

我们将通过 PyTorch（包括 torchvision）来构建整个程序。Mixup 需要的从 beta 分布中生成的样本，我们可以从 NumPy 库中获得。我们还将使用 random 来为 Mixup 寻找随机图像。下面的代码能够导入我们需要的所有库：

"""
Import necessary libraries to train a network using mixup
The code is mainly developed using the PyTorch library
"""
import numpy as np
import pickle
import random
import torch
import torch.nn as nn
import torch.nn.functional as F
import torchvision
import torchvision.transforms as transforms
from torch.utils.data import Dataset, DataLoader

数据集

为了演示，我们将用传统的图像分类任务来说明 Mixup 的强大，那么这种情况下 CIFAR-10 则会是非常理想的数据集。CIFAR-10 包含 10 个类别的 60000 张彩色图像（每类 6000 张），按 5:1 的比例分为训练和测试集。这些图像分类起来相当简单，但比最基本的数字识别数据集 MNIST 要难一些。

有许多方法可以下载 CIFAR-10 数据集，比如多伦多大学网站里就包含了相关数据集。在这里，我推荐大家使用格物钛的公开数据集平台（graviti.cn/open-datasets），因为在这个平台上，如果使用他们的 SDK，不用下载也可以获取免费的数据集资源。事实上，这个公开数据集平台包含了行业内数百个知名的优质数据集，每个数据集都有相关的作者说明，以及不同训练任务的标签，例如分类或目标检测。当然，大家也可以在这个平台下载其他分类数据集，如 CompCars 或 SVHN，来测试 Mixup 在不同场景下的性能。

▲ 数据集广告图（CIFAR-10）

硬件要求

一般来说，我们最好用 GPU（显卡）来训练神经网络，因为它能显著提高训练速度。不过如果只有 CPU 可用，我们还是可以对程序进行简单测试的。如果你想让程序能够自行确定所需硬件，使用以下代码即可：

"""
Determine if any GPUs are available
"""
device = torch.device('cuda' if torch.cuda.is_available() else 'cpu')

实现

网络

这里，我们的目标是要测试 Mixup 的性能，而不是调试网络本身，所以我们只需要简单实现一个 4 层卷积层和 2 层全连接层的卷积神经网络（CNN）即可。为了比较使用和不使用 Mixup 的区别，我们将应用同一个网络来确保比较的准确性。

我们可以使用下列代码来搭建上面所说的简单网络：

"""
Create a simple CNN
"""
class CNN(nn.Module):
    def __init__(self):
        super(CNN, self).__init__()

        # Network consists of 4 convolutional layers followed by 2 fully-connected layers
        self.conv11 = nn.Conv2d(3, 64, 3)
        self.conv12 = nn.Conv2d(64, 64, 3)
        self.conv21 = nn.Conv2d(64, 128, 3)
        self.conv22 = nn.Conv2d(128, 128, 3)
        self.fc1 = nn.Linear(128 * 5 * 5, 256)
        self.fc2 = nn.Linear(256, 10)
    def forward(self, x):
       x = F.relu(self.conv11(x))
       x = F.relu(self.conv12(x))
       x = F.max_pool2d(x, (2,2))
       x = F.relu(self.conv21(x))
       x = F.relu(self.conv22(x))
       x = F.max_pool2d(x, (2,2))

       # Size is calculated based on kernel size 3 and padding 0
       x = x.view(-1, 128 * 5 * 5)
       x = F.relu(self.fc1(x))
       x = self.fc2(x)

       return nn.Sigmoid()(x)

Mixup

Mixup 阶段是在数据集加载过程中完成的，所以我们必须写入我们自己的数据集，而不是使用 torchvision.datasets 所提供的默认数据集。

下面的代码简单地实现了 Mixup，并结合使用了 NumPy 的贝塔函数。

"""
Dataset and Dataloader creation
All data are downloaded found via Graviti Open Dataset which links to CIFAR-10 official page
The dataset implementation is where mixup take place
"""

class CIFAR_Dataset(Dataset):
    def __init__(self, data_dir, train, transform):
        self.data_dir = data_dir
        self.train = train
        self.transform = transform
        self.data = []
        self.targets = []

        # Loading all the data depending on whether the dataset is training or testing
        if self.train:
            for i in range(5):
                with open(data_dir + 'data_batch_' + str(i+1), 'rb') as f:
                    entry = pickle.load(f, encoding='latin1')
                    self.data.append(entry['data'])
                    self.targets.extend(entry['labels'])
        else:
            with open(data_dir + 'test_batch', 'rb') as f:
                entry = pickle.load(f, encoding='latin1')
                self.data.append(entry['data'])
                self.targets.extend(entry['labels'])

        # Reshape it and turn it into the HWC format which PyTorch takes in the images
        # Original CIFAR format can be seen via its official page
        self.data = np.vstack(self.data).reshape(-1, 3, 32, 32)
        self.data = self.data.transpose((0, 2, 3, 1))

    def __len__(self):
        return len(self.data)

    def __getitem__(self, idx):

        # Create a one hot label
        label = torch.zeros(10)
        label[self.targets[idx]] = 1.

        # Transform the image by converting to tensor and normalizing it
        if self.transform:
            image = transform(self.data[idx])

        # If data is for training, perform mixup, only perform mixup roughly on 1 for every 5 images
        if self.train and idx > 0 and idx%5 == 0:

            # Choose another image/label randomly
            mixup_idx = random.randint(0, len(self.data)-1)
            mixup_label = torch.zeros(10)
            label[self.targets[mixup_idx]] = 1.
            if self.transform:
                mixup_image = transform(self.data[mixup_idx])

            # Select a random number from the given beta distribution
            # Mixup the images accordingly
            alpha = 0.2
            lam = np.random.beta(alpha, alpha)
            image = lam * image + (1 - lam) * mixup_image
            label = lam * label + (1 - lam) * mixup_label

        return image, label

需要注意的是，我们并没有对所有的图像都进行 Mixup，而是大概每 5 张处理 1 张。我们还使用了一个 0.2 的贝塔分布。你可以自己为不同的实验改变分布以及被混合的图像的数量，或许你会取得更好的结果!

训练和评估

下面的代码展示的是训练过程。我们将批次大小设置为 128，学习率为 1e-3，总次数为 30 次。整个训练进行了两次，唯一区别是有没有使用 Mixup。需要注意的是， 损失函数需要由我们自己定义，因为目前 BCE 损失不允许使用带有小数的标签。

"""
Initialize the network, loss Adam optimizer
Torch BCE Loss does not support mixup labels (not 1 or 0), so we implement our own
"""
net = CNN().to(device)
optimizer = torch.optim.Adam(net.parameters(), lr=LEARNING_RATE)
def bceloss(x, y):
    eps = 1e-6
    return -torch.mean(y * torch.log(x + eps) + (1 - y) * torch.log(1 - x + eps))
best_Acc = 0


"""
Training Procedure
"""
for epoch in range(NUM_EPOCHS):
    net.train()
    # We train and visualize the loss every 100 iterations
    for idx, (imgs, labels) in enumerate(train_dataloader):
        imgs = imgs.to(device)
        labels = labels.to(device)
        preds = net(imgs)
        loss = bceloss(preds, labels)
        optimizer.zero_grad()
        loss.backward()
        optimizer.step()
        if idx%100 == 0:
            print("Epoch {} Iteration {}, Current Loss: {}".format(epoch, idx, loss))

    # We evaluate the network after every epoch based on test set accuracy
    net.eval()
    with torch.no_grad():
        total = 0
        numCorrect = 0
        for (imgs, labels) in test_dataloader:
            imgs = imgs.to(device)
            labels = labels.to(device)
            preds = net(imgs)
            numCorrect += (torch.argmax(preds, dim=1) == torch.argmax(labels, dim=1)).float().sum()
            total += len(imgs)
        acc = numCorrect/total
        print("Current image classification accuracy at epoch {}: {}".format(epoch, acc))
        if acc > best_Acc:
            best_Acc = acc

为了评估 Mixup 的效果，我们进行了三次对照试验来计算最终的准确性。在没有 Mixup 的情况下，该网络在测试集上的准确率约为 74.5%，而在使用了 Mixup 的情况下，准确率提高到了约 76.5%!

图像分类之外

Mixup 将图像分类的准确性带到了一个前所未有的高度，但研究表明，Mixup 的好处还能延伸到其他计算机视觉任务中，比如对抗性数据的生成和防御。另外也有相关文献在 Mixup 拓展到三维表示中，目前的结果表明 Mixup 在这一领域也十分有效的，例如 PointMixup（arxiv.org/abs/2008.06374）。

结语

由此，我们用 Mixup 做的小实验就大功告成啦！在这篇文章中，我们简单介绍了 Mixup 的概念并演示了如何在图像分类网络训练中应用 Mixup。完整的实现方式可以在这— GitHub 仓库（github.com/ttchengab/mixup.git）中找到。

关于格物钛

格物钛智能科技定位为面向机器学习的数据平台，致力于为 AI 开发者打造下一代新型基础设施，从根本上改变其与非结构化数据的交互方式。我们通过非结构化数据管理工具 TensorBay 和开源数据集社区 Open Datasets，帮助机器学习团队和个人降低数据获取、存储和处理成本，加速 AI 开发和产品创新，为人工智能赋能千行百业、驱动产业升级提供坚实基础。

格物钛｜公开数据集

graviti.cn/open-datasets

 Open Datasets ???? 

微信号｜Graviti_2019

微博｜@格物钛

https://www.graviti.cn

 订阅号：格物钛 ???? 

????

现在，在「知乎」也能找到我们了

进入知乎首页搜索「PaperWeekly」

点击「关注」订阅我们的专栏吧

·