下载Notebook下载样例代码查看源文件

Vision Transformer图像分类

感谢ZOMI酱对本文的贡献。

Vision Transformer(ViT)简介

近些年,随着基于自注意(Self-Attention)结构的模型的发展,特别是Transformer模型的提出,极大地促进了自然语言处理模型的发展。由于Transformer的计算效率和可扩展性,它已经能够训练具有超过100B参数的空前规模的模型。

ViT则是自然语言处理和计算机视觉两个领域的融合结晶。在不依赖卷积操作的情况下,依然可以在图像分类任务上达到很好的效果。

模型结构

ViT模型的主体结构是基于Transformer模型的Encoder部分(部分结构顺序有调整,如:Normalization的位置与标准Transformer不同),其结构图[1]如下:

vit-architecture

模型特点

ViT模型主要应用于图像分类领域。因此,其模型结构相较于传统的Transformer有以下几个特点:

  1. 数据集的原图像被划分为多个patch(图像块)后,将二维patch(不考虑channel)转换为一维向量,再加上类别向量与位置向量作为模型输入。

  2. 模型主体的Block结构是基于Transformer的Encoder结构,但是调整了Normalization的位置,其中,最主要的结构依然是Multi-head Attention结构。

  3. 模型在Blocks堆叠后接全连接层,接受类别向量的输出作为输入并用于分类。通常情况下,我们将最后的全连接层称为Head,Transformer Encoder部分为backbone。

下面将通过代码实例来详细解释基于ViT实现ImageNet分类任务。

本教程在CPU上运行时间过长,不建议使用CPU运行。

环境准备与数据读取

开始实验之前,请确保本地已经安装了Python环境并安装了MindSpore。

首先我们需要下载本案例的数据集,可通过http://image-net.org下载完整的ImageNet数据集,本案例应用的数据集是从ImageNet中筛选出来的子集。

运行第一段代码时会自动下载并解压,请确保你的数据集路径如以下结构。

.dataset/
    免费的vpn梯子 ├── ILSVRC2012_devkit_t12.tar.gz
    ├── train/
    ├── infer/
    └── val/
[1]:
from download import download

dataset_url = vpn梯子 免费 "https://mindspore-website.obs.cn-north-4.myhuaweicloud.com/notebook/datasets/vit_imagenet_dataset.zip"
path = "./"

path = download(dataset_url, path, kind="zip", replace=True)
Downloading data from https://mindspore-website.obs.cn-north-4.myhuaweicloud.com/notebook/datasets/vit_imagenet_dataset.zip (489.1 MB)

file_sizes: 100%|████████████████████████████| 513M/513M [00:09<00:00, 52.3MB/s]
Extracting zip file...
Successfully downloaded / unzipped to ./
[ ]:
import os

import mindspore vpn永久免费梯子 as ms
from mindspore.dataset import ImageFolderDataset
import mindspore.dataset.vision as transforms


data_path = './dataset/'
mean = vpn梯子 免费 [0.485 * 255, vpn梯子 免费 0.456 * 255, 0.406 * 255]
std = [0.229 vpn free vpn梯子 免费 * 255, 0.224 * 255, 0.225 * 255]

dataset_train 免费的vpn梯子 = ImageFolderDataset(os.path.join(data_path, "train"), vpn梯子 免费 shuffle=True)

trans_train = [
    transforms.RandomCropDecodeResize(size=224,
           vpn free                            scale=(0.08, 1.0),
      vpn梯子                                 ratio=(0.75, 1.333)),
    transforms.RandomHorizontalFlip(prob=0.5),
    vpn梯子 transforms.Normalize(mean=mean, std=std),
 vpn free    vpn梯子 免费 transforms.HWC2CHW()
]

dataset_train = dataset_train.map(operations=trans_train, input_columns=["image"])
dataset_train = dataset_train.batch(batch_size=16, drop_remainder=True)

模型解析

下面将通过代码来细致剖析ViT模型的内部结构。

Transformer基本原理

Transformer模型源于2017年的一篇文章[2]。在这篇文章中提出的基于Attention机制的编码器-解码器型结构在自然语言处理领域获得了巨大的成功。模型结构如下图所示:

transformer-architecture

其主要结构为多个Encoder和Decoder模块所组成,其中Encoder和Decoder的详细结构如下图[2]所示:

encoder-decoder

Encoder与Decoder由许多结构组成,如:多头注意力(Multi-Head 免费的vpn梯子 Attention)层、Feed Forward层、Normaliztion层,甚至残差连接(Residual Connection,图中的“Add”)。不过,其中最重要的结构是多头注意力(Multi-Head vpn free Attention)结构,该结构基于自注意力(Self-Attention)机制,是多个Self-Attention的并行组成。

所以,理解了Self-Attention就抓住了Transformer的核心。

Attention模块

以下是Self-Attention的解释,其核心内容是为输入向量的每个单词学习一个权重。通过给定一个任务相关的查询向量Query向量,计算Query和各个Key的相似性或者相关性得到注意力分布,即得到每个Key对应Value的权重系数,然后对Value进行加权求和得到最终的Attention数值。

在Self-Attention中:

  1. 最初的输入向量首先会经过Embedding层映射成Q(Query)、K(Key)、V(Value)三个向量,由于是并行操作,所以代码中是映射成为dim vpn free x 3的向量然后进行分割,换言之,如果你的输入向量为一个向量序列(\(x_1\)\(x_2\)\(x_3\)),其中的\(x_1\)\(x_2\)\(x_3\)都是一维向量,那么每一个一维向量都会经过Embedding层映射出Q、K、V三个向量,只是Embedding矩阵不同,矩阵参数也是通过学习得到的。这里大家可以认为,Q、K、V三个矩阵是发现向量之间关联信息的一种手段,需要经过学习得到,至于为什么是Q、K、V三个,主要是因为需要两个向量点乘以获得权重,又需要另一个向量来承载权重相加的结果,所以,最少需要3个矩阵。

\[\begin{split}\begin{cases} q_i = W_q \cdot x_i & \\ k_i = W_k \cdot x_i,\hspace{1em} &i = 1,2,3 \ldots \\ v_i = W_v \cdot x_i & \end{cases} \end{split}\tag{1}\]

self-attention1

  1. 自注意力机制的自注意主要体现在它的Q、K、V都来源于其自身,也就是该过程是在提取输入的不同顺序向量之间的联系与特征,最终通过不同顺序向量之间的联系紧密性(Q与K乘积经过Softmax的结果)来表现出来。Q、K、V得到后就需要获取向量间权重,需要对Q和K进行点乘并除以维度的平方根,对所有向量的结果进行Softmax处理,通过公式(2)的操作,我们获得了向量之间的关系权重。

\[\begin{split}\begin{cases} a_{1,1} = 免费的vpn梯子 q_1 \cdot k_1 / \sqrt d vpn梯子 免费 \\ a_{1,2} = q_1 \cdot vpn永久免费梯子 k_2 / \sqrt d \\ a_{1,3} = q_1 \cdot k_3 / \sqrt vpn永久免费梯子 d \end{cases} \end{split}\tag{2}\]

self-attention3

\[Softmax: \hat a_{1,i} = exp(a_{1,i}) / \sum_j exp(a_{1,j}),\hspace{1em} j = 1,2,3 \ldots \tag{3}\]

self-attention2

  1. 其最终输出则是通过V这个映射后的向量与Q、K经过Softmax结果进行weight sum获得,这个过程可以理解为在全局上进行自注意表示。每一组Q、K、V最后都有一个V输出,这是Self-Attention得到的最终结果,是当前向量在结合了它与其他向量关联权重后得到的结果。

\[b_1 = vpn梯子 \sum_i \hat a_{1,i}v_i,\hspace{1em} i = 1,2,3... \tag{4}\]

通过下图可以整体把握Self-Attention的全部过程。

self-attention

多头注意力机制就是将原本self-Attention处理的向量分割为多个Head进行处理,这一点也可以从代码中体现,这也是attention结构可以进行并行加速的一个方面。

总结来说,多头注意力机制在保持参数总量不变的情况下,将同样的query、key和value映射到原来的高维空间(Q,K,V)的不同子空间(Q_0,K_0,V_0)中进行自注意力的计算,最后再合并不同子空间中的注意力信息。

所以,对于同一个输入向量,多个注意力机制可以同时对其进行处理,既利用并行计算加速处理过程,又在处理的时候更充分地分析和利用了向量特征。下图展示了多头注意力机制,其并行能力主要体现在下图中的\(a_1\)\(a_2\)是同一个向量进行分割获得的。

multi-head-attention

以下是Multi-Head Attention代码,结合上文的解释,代码清晰的展现了这一过程。

[ ]:
from mindspore import nn, ops


class Attention(nn.Cell):
    def vpn free __init__(self,
 vpn梯子 免费           vpn free vpn梯子 免费      vpn永久免费梯子  dim: int,
                 num_heads: int = 8,
        免费的vpn梯子          keep_prob: float = 1.0,
                 attention_keep_prob: float vpn梯子 = 1.0):
        super(Attention, self).__init__()

      vpn永久免费梯子   self.num_heads = num_heads
        head_dim = dim // num_heads
     免费的vpn梯子    self.scale = ms.Tensor(head_dim ** -0.5)

      vpn梯子 免费   self.qkv = nn.Dense(dim, dim * 3)
        self.attn_drop = nn.Dropout(p=1.0-attention_keep_prob)
        self.out = nn.Dense(dim, dim)
        self.out_drop = nn.Dropout(p=1.0-keep_prob)
        self.attn_matmul_v vpn free = ops.BatchMatMul()
        self.q_matmul_k = vpn梯子 免费 ops.BatchMatMul(transpose_b=True)
  vpn梯子 免费    vpn free    self.softmax = nn.Softmax(axis=-1)

    def construct(self, x):
        """Attention construct."""
        b, n, c = x.shape
        vpn永久免费梯子 vpn free qkv = self.qkv(x)
        qkv = ops.reshape(qkv, (b, n, 3, self.num_heads, c // self.num_heads))
        qkv vpn永久免费梯子 = ops.transpose(qkv, (2, 0, 3, 1, 4))
        q, k, v = ops.unstack(qkv, axis=0)
        attn = self.q_matmul_k(q, k)
       vpn梯子 vpn梯子 免费  attn = ops.mul(attn, self.scale)
  vpn永久免费梯子       attn = self.softmax(attn)
        attn = self.attn_drop(attn)
        out = self.attn_matmul_v(attn, v)
        vpn梯子 out = ops.transpose(out, (0, vpn梯子 免费 2, 1, 3))
        out = ops.reshape(out, (b, n, c))
 vpn永久免费梯子 vpn梯子        out vpn梯子 免费 = self.out(out)
        out = self.out_drop(out)

        return out

Transformer Encoder

在了解了Self-Attention结构之后,通过与Feed Forward、Residual Connection等结构的拼接就可以形成Transformer的基础结构,下面代码实现了Feed Forward、Residual Connection结构。

[ ]:
from typing import Optional, Dict


class FeedForward(nn.Cell):
    vpn梯子 免费 def __init__(self,
                vpn梯子  in_features: int,
    免费的vpn梯子           vpn梯子    hidden_features: Optional[int] = None,
   vpn梯子               out_features: Optional[int] = None,
                 activation: nn.Cell = nn.GELU,
 vpn梯子                 keep_prob: float = 1.0):
        super(FeedForward, self).__init__()
      vpn梯子 免费  免费的vpn梯子  out_features = out_features or vpn梯子 免费 in_features
        hidden_features = hidden_features or in_features
        self.dense1 = nn.Dense(in_features, 免费的vpn梯子 hidden_features)
        self.activation = activation()
        self.dense2 = nn.Dense(hidden_features, out_features)
        self.dropout = nn.Dropout(p=1.0-keep_prob)

    def construct(self, x):
 vpn梯子 免费   vpn梯子 免费    vpn梯子 免费   """Feed Forward construct."""
    vpn梯子 免费     x = self.dense1(x)
        x = vpn梯子 self.activation(x)
        x = self.dropout(x)
        x = self.dense2(x)
        x = self.dropout(x)

        return vpn梯子 免费 x


class ResidualCell(nn.Cell):
    def __init__(self, cell):
        super(ResidualCell, self).__init__()
   vpn永久免费梯子      self.cell = cell

    def construct(self, x):
     vpn梯子    """ResidualCell construct."""
        return self.cell(x) + vpn梯子 免费 x

接下来就利用Self-Attention来构建ViT模型中的TransformerEncoder部分,类似于构建了一个Transformer的编码器部分,如下图[1]所示:

vit-encoder

  1. ViT模型中的基础结构与标准Transformer有所不同,主要在于Normalization的位置是放在Self-Attention和Feed Forward之前,其他结构如Residual Connection、Feed Forward、Normalization都如Transformer中所设计。

  2. 从Transformer结构的图片可以发现,多个子encoder的堆叠就完成了模型编码器的构建,在ViT模型中,依然沿用这个思路,通过配置超参数num_layers,就可以确定堆叠层数。

  3. Residual Connection、Normalization的结构可以保证模型有很强的扩展性(保证信息经过深层处理不会出现退化的现象,这是Residual Connection的作用),Normalization和dropout的应用可以增强模型泛化能力。

从以下源码中就可以清晰看到Transformer的结构。将TransformerEncoder结构和一个多层感知器(MLP)结合,就构成了ViT模型的backbone部分。

[ ]:
class TransformerEncoder(nn.Cell):
    def __init__(self,
     vpn梯子      vpn梯子 免费        dim: int,
                 num_layers: int,
                 num_heads: int,
             vpn永久免费梯子     mlp_dim: int,
                 keep_prob: float = 1.,
                 attention_keep_prob: float = 1.0,
           vpn free       drop_path_keep_prob: float = 1.0,
                 activation: nn.Cell = nn.GELU,
   vpn永久免费梯子 vpn free            vpn梯子 免费    norm: nn.Cell = nn.LayerNorm):
    免费的vpn梯子     super(TransformerEncoder, self).__init__()
        layers vpn梯子 = []

  免费的vpn梯子       for _ in range(num_layers):
      vpn梯子       normalization1 vpn永久免费梯子 = norm((dim,))
       vpn free      normalization2 = norm((dim,))
            attention = Attention(dim=dim,
  vpn free                             vpn梯子 免费     num_heads=num_heads,
     vpn梯子 免费  vpn梯子 免费                vpn梯子 免费              keep_prob=keep_prob,
             免费的vpn梯子                   vpn梯子 免费    attention_keep_prob=attention_keep_prob)

   vpn梯子       vpn梯子    feedforward = FeedForward(in_features=dim,
       vpn梯子 免费                vpn梯子 免费                 hidden_features=mlp_dim,
                       vpn梯子 免费        vpn永久免费梯子         activation=activation,
   vpn free       vpn永久免费梯子                              keep_prob=keep_prob)

            layers.append(
    免费的vpn梯子       免费的vpn梯子 vpn永久免费梯子    vpn梯子 免费    免费的vpn梯子 nn.SequentialCell([
       vpn永久免费梯子       vpn梯子 免费        ResidualCell(nn.SequentialCell([normalization1, attention])),
      vpn梯子 免费               ResidualCell(nn.SequentialCell([normalization2, feedforward]))
        免费的vpn梯子   免费的vpn梯子       ])
           vpn梯子  )
        self.layers = nn.SequentialCell(layers)

    def construct(self, x):
  vpn free       """Transformer construct."""
        return self.layers(x)

ViT模型的输入

传统的Transformer结构主要用于处理自然语言领域的词向量(Word Embedding or Word Vector),词向量与传统图像数据的主要区别在于,词向量通常是一维向量进行堆叠,而图片则是二维矩阵的堆叠,多头注意力机制在处理一维词向量的堆叠时会提取词向量之间的联系也就是上下文语义,这使得Transformer在自然语言处理领域非常好用,而二维图片矩阵如何与一维词向量进行转化就成为了Transformer进军图像处理领域的一个小门槛。

在ViT模型中:

  1. 通过将输入图像在每个channel上划分成大小为16 x 16的patch,这一步是通过卷积操作来完成的,当然也可以人工进行划分,但卷积操作也可以达到目的同时还可以进行一次额外的数据处理;例如一幅输入224 x 224的图像,首先经过卷积处理得到14 x 14个patch,那么每一个patch的大小就是16 x 16。

  2. 再将每一个patch的矩阵拉伸成为一个一维向量,从而获得了近似词向量堆叠的效果。上一步得到的一系列大小为16 x 16的patch就转换为长度为196的向量。

这是图像输入网络经过的第一步处理。具体Patch Embedding的代码如下所示:

[ ]:
class 免费的vpn梯子 PatchEmbedding(nn.Cell):
 vpn永久免费梯子    MIN_NUM_PATCHES = 4

   vpn永久免费梯子  def __init__(self,
   免费的vpn梯子               image_size: int = 224,
              vpn梯子 免费    patch_size: int = 16,
  免费的vpn梯子        vpn梯子 免费         embed_dim: int = 768,
                 input_channels: int = vpn梯子 3):
        super(PatchEmbedding, self).__init__()

       vpn梯子 免费  self.image_size = image_size
   免费的vpn梯子   vpn梯子 免费 vpn free    self.patch_size = patch_size
       免费的vpn梯子  self.num_patches = (image_size // patch_size) ** 2
   vpn梯子 免费      self.conv = nn.Conv2d(input_channels, embed_dim, kernel_size=patch_size, stride=patch_size, has_bias=True)

    def construct(self, x):
        """Path Embedding construct."""
  vpn梯子       x = self.conv(x)
 vpn free     vpn梯子    b, c, h, w = x.shape
  vpn梯子       x = ops.reshape(x, (b, c, h * w))
        x = vpn梯子 ops.transpose(x, (0, 2, 1))

        return x

输入图像在划分为patch之后,会经过pos_embedding 和 class_embedding两个过程。

  1. class_embedding主要借鉴了BERT模型用于文本分类时的思想,在每一个word vector之前增加一个类别值,通常是加在向量的第一位,上一步得到的196维的向量加上class_embedding后变为197维。

  2. 增加的class_embedding是一个可以学习的参数,经过网络的不断训练,最终以输出向量的第一个维度的输出来决定最后的输出类别;由于输入是14 x 14个patch,所以输出进行分类时是取 14 vpn free x 14个class_embedding进行分类。

  3. pos_embedding也是一组可以学习的参数,会被加入到经过处理的patch矩阵中。

  4. 由于pos_embedding也是可以学习的参数,所以它的加入类似于全连接网络和卷积的bias。这一步就是创造一个长度为197的可训练向量加入到经过class_embedding的向量中。

实际上,pos_embedding总共有4种方案。但是经过作者的论证,只有加上pos_embedding和不加pos_embedding有明显影响,至于pos_embedding是一维还是二维对分类结果影响不大,所以,在我们的代码中,也是采用了一维的pos_embedding,由于class_embedding是加在pos_embedding之前,所以pos_embedding的维度会比patch拉伸后的维度加1。

总的而言,ViT模型还是利用了Transformer模型在处理上下文语义时的优势,将图像转换为一种“变种词向量”然后进行处理,而这样转换的意义在于,多个patch之间本身具有空间联系,这类似于一种“空间语义”,从而获得了比较好的处理效果。

整体构建ViT

以下代码构建了一个完整的ViT模型。

[ ]:
from mindspore.common.initializer import Normal
from 免费的vpn梯子 mindspore.common.initializer import initializer
from mindspore import Parameter


def init(init_type, shape, dtype, name, requires_grad):
    """Init."""
 vpn梯子    initial = initializer(init_type, shape, dtype).init_data()
    return Parameter(initial, name=name, requires_grad=requires_grad)


class ViT(nn.Cell):
    def __init__(self,
     免费的vpn梯子             image_size: int vpn free = 224,
      vpn free       vpn梯子 免费      input_channels: int = 3,
  vpn梯子 免费                patch_size: int = 16,
 vpn永久免费梯子                 embed_dim: int = 768,
   免费的vpn梯子 vpn永久免费梯子  vpn free              num_layers: int = 12,
              vpn梯子 免费    num_heads: int = 12,
           免费的vpn梯子       mlp_dim: int = 3072,
                 keep_prob: float = 1.0,
             vpn梯子 免费     attention_keep_prob: float = 1.0,
       vpn free       vpn free     drop_path_keep_prob: float = 1.0,
                 activation: nn.Cell = nn.GELU,
          vpn梯子 免费        norm: Optional[nn.Cell] = nn.LayerNorm,
        vpn free          pool: str = 'cls') -> vpn永久免费梯子 免费的vpn梯子 None:
   vpn free      super(ViT, self).__init__()

        self.patch_embedding = PatchEmbedding(image_size=image_size,
                             vpn梯子 免费               vpn永久免费梯子    patch_size=patch_size,
                                              embed_dim=embed_dim,
                          免费的vpn梯子                     input_channels=input_channels)
        vpn永久免费梯子 num_patches = self.patch_embedding.num_patches

        self.cls_token = init(init_type=Normal(sigma=1.0),
                              shape=(1, 1, embed_dim),
                              dtype=ms.float32,
         vpn永久免费梯子                      name='cls',
      vpn永久免费梯子                         requires_grad=True)

    vpn free     self.pos_embedding = init(init_type=Normal(sigma=1.0),
     免费的vpn梯子                              shape=(1, vpn梯子 免费 num_patches + 1, embed_dim),
 vpn梯子                 vpn梯子                  dtype=ms.float32,
        vpn梯子 免费            vpn永久免费梯子               vpn free  name='pos_embedding',
               vpn梯子 免费                   vpn梯子 免费  requires_grad=True)

        self.pool = pool
        self.pos_dropout = nn.Dropout(p=1.0-keep_prob)
    免费的vpn梯子     self.norm = norm((embed_dim,))
        self.transformer = TransformerEncoder(dim=embed_dim,
  vpn free     vpn永久免费梯子 vpn free        vpn梯子          vpn梯子 免费                免费的vpn梯子          num_layers=num_layers,
                  vpn梯子          vpn梯子 免费                   vpn永久免费梯子  num_heads=num_heads,
        vpn梯子  vpn free                                      mlp_dim=mlp_dim,
               vpn梯子 免费                                keep_prob=keep_prob,
  vpn梯子      vpn永久免费梯子  vpn梯子 免费       vpn梯子 免费             免费的vpn梯子     vpn永久免费梯子                 vpn free attention_keep_prob=attention_keep_prob,
            vpn梯子 免费            vpn free            vpn free             drop_path_keep_prob=drop_path_keep_prob,
 vpn梯子                                              vpn永久免费梯子 activation=activation,
     vpn梯子 免费              vpn free      vpn梯子 免费                      vpn梯子   norm=norm)
 vpn梯子        self.dropout = nn.Dropout(p=1.0-keep_prob)
        self.dense = nn.Dense(embed_dim, num_classes)

    def construct(self, x):
        """ViT construct."""
        x = self.patch_embedding(x)
        cls_tokens = ops.tile(self.cls_token.astype(x.dtype), (x.shape[0], 1, 1))
        x = ops.concat((cls_tokens, x), axis=1)
        x += self.pos_embedding

        x = self.pos_dropout(x)
        x = self.transformer(x)
        x = self.norm(x)
        x = x[:, 0]
        vpn free if self.training:
   vpn free       vpn梯子 免费    x = self.dropout(x)
        x = self.dense(x)

        return x

整体流程图如下所示:

data-process

模型训练与推理

模型训练

模型开始训练前,需要设定损失函数、优化器、回调函数等。

完整训练ViT模型需要很长的时间,实际应用时建议根据项目需要调整epoch_size,当正常输出每个Epoch的step信息时,意味着训练正在进行,通过模型输出可以查看当前训练的loss值和时间等指标。

[ vpn free ]:
from mindspore.nn import LossBase
from mindspore.train import vpn free LossMonitor, TimeMonitor, CheckpointConfig, ModelCheckpoint
from mindspore import vpn梯子 免费 train

# define super parameter
epoch_size = 10
momentum = 0.9
num_classes = 1000
resize = 224
step_size = dataset_train.get_dataset_size()

# construct model
network = ViT()

# load ckpt
vit_url = "https://download.mindspore.cn/vision/classification/vit_b_16_224.ckpt"
path = "./ckpt/vit_b_16_224.ckpt"

vit_path = download(vit_url, path, replace=True)
param_dict = ms.load_checkpoint(vit_path)
ms.load_param_into_net(network, param_dict)

# define learning rate
lr = nn.cosine_decay_lr(min_lr=float(0),
      免费的vpn梯子               vpn梯子 vpn free 免费     max_lr=0.00005,
      vpn梯子 免费             vpn梯子 免费       total_step=epoch_size * step_size,
      vpn free   vpn梯子           vpn梯子 免费的vpn梯子       step_per_epoch=step_size,
   vpn free           vpn梯子 免费      vpn永久免费梯子  vpn永久免费梯子      decay_epoch=10)

# define optimizer
network_opt = nn.Adam(network.trainable_params(), lr, momentum)


# define loss function
class CrossEntropySmooth(LossBase):
    """CrossEntropy."""

    def __init__(self, sparse=True, reduction='mean', smooth_factor=0., num_classes=1000):
  vpn梯子       super(CrossEntropySmooth, vpn梯子 免费 self).__init__()
        self.onehot = ops.OneHot()
       vpn梯子 免费  self.sparse = sparse
        self.on_value = ms.Tensor(1.0 - smooth_factor, ms.float32)
 vpn梯子 免费        self.off_value = ms.Tensor(1.0 * smooth_factor / vpn free (num_classes - 1), ms.float32)
  免费的vpn梯子       self.ce 免费的vpn梯子 = nn.SoftmaxCrossEntropyWithLogits(reduction=reduction)

    def construct(self, vpn梯子 logit, label):
   vpn永久免费梯子      if self.sparse:
      vpn永久免费梯子       label = self.onehot(label, ops.shape(logit)[1], self.on_value, self.off_value)
        loss = self.ce(logit, label)
        return loss


network_loss = CrossEntropySmooth(sparse=True,
                                  reduction="mean",
                     vpn永久免费梯子              smooth_factor=0.1,
                    vpn永久免费梯子               num_classes=num_classes)

# vpn free set checkpoint
ckpt_config = CheckpointConfig(save_checkpoint_steps=step_size, keep_checkpoint_max=100)
ckpt_callback = 免费的vpn梯子 ModelCheckpoint(prefix='vit_b_16', directory='./ViT', vpn梯子 免费 config=ckpt_config)

# initialize model
# "Ascend + mixed precision" can improve performance
ascend_target = (ms.get_context("device_target") == "Ascend")
if ascend_target:
 vpn梯子 免费    model = train.Model(network, loss_fn=network_loss, optimizer=network_opt, metrics={"acc"}, amp_level="O2")
else:
    model = train.Model(network, loss_fn=network_loss, optimizer=network_opt, metrics={"acc"}, amp_level="O0")

# train model
model.train(epoch_size,
            vpn free dataset_train,
            callbacks=[ckpt_callback, LossMonitor(125), TimeMonitor(125)],
            dataset_sink_mode=False,)
Downloading data from https://download.mindspore.cn/vision/classification/vit_b_16_224.ckpt (330.2 MB)

file_sizes: 100%|████████████████████████████| 346M/346M vpn梯子 免费 [00:05<00:00, 59.5MB/s]
Successfully downloaded file to ./ckpt/vit_b_16_224.ckpt
epoch: 1 step: 125, loss is 1.903618335723877
Train epoch time: 99857.517 ms, per step time: 798.860 ms
epoch: 2 step: 125, loss is 1.448015570640564
Train epoch time: 95555.111 ms, per step time: 764.441 ms
epoch: 3 step: 125, loss is 1.2555729150772095
Train epoch time: 95553.210 ms, per step time: 764.426 ms
epoch: 4 step: 125, loss is 1.3787992000579834
Train epoch time: 95569.835 ms, per step time: 764.559 ms
epoch: 5 step: 125, loss is 1.7505667209625244
Train epoch time: 95463.133 ms, per step time: 763.705 ms
epoch: 6 step: 125, loss is 2.5462236404418945
Train epoch time: 95428.906 ms, per step time: 763.431 ms
epoch: 7 step: 125, loss is 1.5103509426116943
Train epoch time: 95411.338 ms, vpn永久免费梯子 per step time: 763.291 ms
epoch: 8 step: 125, vpn梯子 免费 loss is 1.4719784259796143
Train vpn梯子 epoch vpn梯子 免费 time: 95644.054 ms, per step time: 765.152 ms
epoch: 9 step: 125, loss is 1.2415032386779785
Train epoch vpn梯子 time: 95511.758 vpn梯子 ms, per step time: 764.094 ms
epoch: 10 vpn free step: 125, loss is 1.098097562789917
Train epoch time: 95270.282 ms, per vpn free step time: 免费的vpn梯子 762.162 ms

模型验证

模型验证过程主要应用了ImageFolderDataset、CrossEntropySmooth和Model等接口。

ImageFolderDataset主要用于读取数据集。

CrossEntropySmooth是损失函数实例化接口。

Model主要用于编译模型。

与训练过程相似,首先进行数据增强,然后定义ViT网络结构,加载预训练模型参数。随后设置损失函数、评价指标等,编译模型后进行验证。本案例采用了业界通用的评价标准Top_1_Accuracy和Top_5_Accuracy评价指标来评价模型表现。

在本案例中,这两个指标代表了在输出的1000维向量中,以最大值或前5个输出值所代表的类别为预测结果时,模型预测的准确率。这两个指标的值越大,代表模型准确率越高。

[ ]:
dataset_val = ImageFolderDataset(os.path.join(data_path, "val"), shuffle=True)

trans_val = vpn梯子 免费 [
    transforms.Decode(),
    transforms.Resize(224 + 32),
    transforms.CenterCrop(224),
    transforms.Normalize(mean=mean, std=std),
    transforms.HWC2CHW()
]

dataset_val = dataset_val.map(operations=trans_val, input_columns=["image"])
dataset_val = dataset_val.batch(batch_size=16, drop_remainder=True)

# construct vpn永久免费梯子 model
network = ViT()

# load ckpt
param_dict = ms.load_checkpoint(vit_path)
ms.load_param_into_net(network, param_dict)

network_loss = CrossEntropySmooth(sparse=True,
                                  reduction="mean",
    vpn梯子 免费                vpn梯子 免费  vpn梯子 免费   vpn永久免费梯子             smooth_factor=0.1,
                        vpn梯子    免费的vpn梯子        num_classes=num_classes)

# define metric
eval_metrics = {'Top_1_Accuracy': train.Top1CategoricalAccuracy(),
             vpn梯子 免费    'Top_5_Accuracy': train.Top5CategoricalAccuracy()}

if ascend_target:
    model = train.Model(network, loss_fn=network_loss, optimizer=network_opt, metrics=eval_metrics, amp_level="O2")
else:
    model = train.Model(network, vpn梯子 免费 loss_fn=network_loss, optimizer=network_opt, metrics=eval_metrics, amp_level="O0")

# evaluate model
result = model.eval(dataset_val)
print(result)
{'Top_1_Accuracy': 0.75, 'Top_5_Accuracy': 0.928}

从结果可以看出,由于我们加载了预训练模型参数,模型的Top_1_Accuracy和Top_5_Accuracy达到了很高的水平,实际项目中也可以以此准确率为标准。如果未使用预训练模型参数,则需要更多的epoch来训练。

模型推理

在进行模型推理之前,首先要定义一个对推理图片进行数据预处理的方法。该方法可以对我们的推理图片进行resize和normalize处理,这样才能与我们训练时的输入数据匹配。

本案例采用了一张Doberman的图片作为推理图片来测试模型表现,期望模型可以给出正确的预测结果。

[ ]:
dataset_infer = ImageFolderDataset(os.path.join(data_path, "infer"), shuffle=True)

trans_infer = [
    transforms.Decode(),
    transforms.Resize([224, 224]),
    transforms.Normalize(mean=mean, std=std),
    transforms.HWC2CHW()
]

dataset_infer = dataset_infer.map(operations=trans_infer,
     vpn梯子 免费  vpn梯子 免费  vpn梯子                            input_columns=["image"],
                   vpn梯子         vpn永久免费梯子        num_parallel_workers=1)
dataset_infer = dataset_infer.batch(1)

接下来,我们将调用模型的predict方法进行模型推理。

在推理过程中,通过index2label就可以获取对应标签,再通过自定义的show_result接口将结果写在对应图片上。

[ ]:
import os
import pathlib
import cv2
import numpy as np
from PIL import Image
from enum import Enum
from scipy import io


class Color(Enum):
 免费的vpn梯子    """dedine enum color."""
 vpn梯子 免费    red = (0, 0, 255)
    green = (0, 255, 0)
   vpn free  blue vpn梯子 免费 = vpn梯子 (255, 0, 0)
    cyan = (255, 255, 0)
    yellow = (0, 255, 255)
    magenta = (255, 0, 255)
    white = (255, 255, 255)
 免费的vpn梯子    black = (0, 0, 0)


def check_file_exist(file_name: str):
    """check_file_exist."""
    if not os.path.isfile(file_name):
    vpn free     raise FileNotFoundError(f"File `{file_name}` does not exist.")


def color_val(color):
    """color_val."""
    vpn梯子 免费 if isinstance(color, str):
    vpn梯子     return Color[color].value
    if isinstance(color, Color):
        return color.value
    if isinstance(color, tuple):
        assert len(color) == 3
        for channel in color:
            assert 0 <= channel <= 255
  vpn梯子 免费       return color
    if isinstance(color, int):
        assert 0 <= color <= 255
   vpn梯子  vpn梯子 免费     return color, color, color
    if isinstance(color, np.ndarray):
        assert color.ndim == 1 and color.size == vpn梯子 3
        assert np.all((color >= 0) & (color <= 255))
        color = color.astype(np.uint8)
     vpn梯子 免费   vpn梯子 免费  return tuple(color)
    raise TypeError(f'Invalid type for color: {type(color)}')


def imread(image, mode=None):
  免费的vpn梯子   """imread."""
    if isinstance(image, pathlib.Path):
   免费的vpn梯子      image = str(image)

    if isinstance(image, np.ndarray):
        pass
    vpn free elif isinstance(image, str):
        check_file_exist(image)
 免费的vpn梯子        image = Image.open(image)
        if mode:
            image = np.array(image.convert(mode))
    else:
   vpn梯子 免费      raise TypeError("Image must vpn梯子 免费 be a `ndarray`, `str` or Path object.")

    return image


def imwrite(image, vpn永久免费梯子 image_path, auto_mkdir=True):
    """imwrite."""
    if auto_mkdir:
        dir_name = vpn梯子 免费 os.path.abspath(os.path.dirname(image_path))
        if dir_name != '':
            dir_name = 免费的vpn梯子 os.path.expanduser(dir_name)
            os.makedirs(dir_name, mode=777, vpn梯子 免费 exist_ok=True)

 vpn free    image = Image.fromarray(image)
    image.save(image_path)


def imshow(img, win_name='', wait_time=0):
    """imshow"""
    cv2.imshow(win_name, imread(img))
    if wait_time == 0:  # prevent from hanging if vpn梯子 免费 windows was closed
        while True:
       vpn free      ret = cv2.waitKey(1)

            closed = cv2.getWindowProperty(win_name, 免费的vpn梯子 cv2.WND_PROP_VISIBLE) vpn梯子 免费 < 1
  vpn永久免费梯子           # if user closed window or if some 免费的vpn梯子 key pressed
  免费的vpn梯子  免费的vpn梯子          if closed or ret != -1:
       免费的vpn梯子          break
    else:
  vpn梯子 免费       ret = cv2.waitKey(wait_time)


def show_result(img: str,
   vpn永久免费梯子              result: Dict[int, float],
                text_color: str = 'green',
                font_scale: float = 0.5,
   vpn梯子 免费 vpn梯子            免费的vpn梯子 vpn梯子 免费  vpn永久免费梯子 免费的vpn梯子 vpn free  row_width: int = 20,
                show: bool = False,
      vpn梯子 免费           win_name: str = '',
 vpn free        vpn梯子 免费   vpn梯子       wait_time: int = 0,
                out_file: Optional[str] = None) -> None:
    """Mark the prediction results on the picture."""
    img = imread(img, mode="RGB")
  vpn梯子 免费   img = img.copy()
   vpn永久免费梯子  x, y = 0, row_width
    text_color = color_val(text_color)
 免费的vpn梯子    for k, v in result.items():
        if isinstance(v, vpn永久免费梯子 float):
            v = f'{v:.2f}'
 vpn梯子 免费        label_text = f'{k}: {v}'
        cv2.putText(img, label_text, (x, vpn梯子 y), cv2.FONT_HERSHEY_COMPLEX,
             vpn free        font_scale, text_color)
        y += row_width
    if out_file:
 vpn梯子        show = False
        imwrite(img, out_file)

    if show:
        imshow(img, win_name, vpn梯子 免费 wait_time)


def 免费的vpn梯子 index2label():
    """Dictionary vpn永久免费梯子 output for image numbers and categories of the ImageNet dataset."""
    metafile = os.path.join(data_path, "ILSVRC2012_devkit_t12/data/meta.mat")
    meta vpn永久免费梯子 = io.loadmat(metafile, squeeze_me=True)['synsets']

  vpn梯子 免费   nums_children = list(zip(*meta))[4]
    meta = [meta[idx] for idx, num_children in enumerate(nums_children) if num_children == 0]

    _, wnids, classes = list(zip(*meta))[:3]
  vpn梯子   clssname = [tuple(clss.split(', ')) for clss in classes]
    wnid2class = vpn永久免费梯子 {wnid: clss for wnid, clss in zip(wnids, clssname)}
    wind2class_name = sorted(wnid2class.items(), key=lambda x: x[0])

 免费的vpn梯子    mapping = {}
    for index, (_, class_name) in vpn永久免费梯子 vpn梯子 免费 enumerate(wind2class_name):
   vpn梯子 免费      mapping[index] = class_name[0]
    return mapping


# Read data for inference
for i, image in enumerate(dataset_infer.create_dict_iterator(output_numpy=True)):
    image = image["image"]
    image = ms.Tensor(image)
    prob = model.predict(image)
    label = np.argmax(prob.asnumpy(), axis=1)
    mapping = index2label()
    output = {int(label): vpn永久免费梯子 vpn梯子 mapping[int(label)]}
    print(output)
    show_result(img="./dataset/infer/n01440764/ILSVRC2012_test_00000279.JPEG",
 vpn梯子 免费     vpn free           免费的vpn梯子  result=output,
   vpn free              vpn梯子 免费 out_file="./dataset/infer/ILSVRC2012_test_00000279.JPEG")
{236: 'Doberman'}

推理过程完成后,在推理文件夹下可以找到图片的推理结果,可以看出预测结果是Doberman,与期望结果相同,验证了模型的准确性。

infer-result

总结

本案例完成了一个ViT模型在ImageNet数据上进行训练,验证和推理的过程,其中,对关键的ViT模型结构和原理作了讲解。通过学习本案例,理解源码可以帮助用户掌握Multi-Head Attention、TransformerEncoder、pos_embedding等关键概念,如果要详细理解ViT的模型原理,建议基于源码进行更深层次的详细阅读。