Swin Transformer: Hierarchical Vision Transformer using Shifted Windows
论文名称:Swin Transformer: Hierarchical Vision Transformer using Shifted Windows
作者:Ze Liu ,Yutong Lin,Yue Cao,Han Hu,Yixuan Wei,Zheng Zhang,Stephen Lin,Baining Guo
动机与介绍
动机:
正如无数物理学家为之奋斗的大一统理论,物理学和人类文明的进步都能从中受益很多;那么是否存在一个统一的基础模型在 NLP 与 CV 之间呢?本文正是为此做出了探索。
自 AlexNet 以来,CNN 作为骨干(backbone)在计算机视觉中得到了广泛应用;另一方面,自然语言处理中的网络结构的演变则走了一条不同的道路,自 2017 年以来,Transformer 便一直占据统治地位。
Transformer 是为序列建模和转换任务而设计的,它以关注数据中的长期依赖关系而著称。其在 NLP 领域的巨大成功吸引了人们研究它对 CV 的适应性,最近的实验显示其在图像分类和联合视觉语言建模方面有所成效。
同样,也有很多 NLP 的工作尝试使用卷积建模,为实现 NLP/CV 共享基础模型而做出努力。
NLP 与 CV 的主要区别:

- multi-scale:在 CV 中,多尺度是十分重要的,而在 NLP 中,每个词都被表示为相同大小的 token,表现为尺度不变;
- locality:在一张图片中,相近的区域往往会很相似,表现为空间平滑性,而在一个句子中,相邻的词语却很少拥有相似的意思;
- translation invariance:众所周知,卷积的平移不变性带来了很好的泛化能力,不论猫在图片的哪一个位置,都会被识别出来,而在 NLP 中,相同词语所在句子的不同位置中,则表示完全不 同的意思。
利用上述性质,将 Transformer 和 CNN 中优秀的先验相结合,提出了 Swin Transformer.
本文的主要贡献有:
- 提出了一种分层 Transformer,其可以作为计算机视觉的通用主干网络,并且在各类下游任务上取得 SOTA;
- 通过 Shift Windows 实现了对输入图像尺寸的线性时间复杂度。

Method
整体结构

上图是 Swin Transformer 中最小版本的可视化结构图,其主要流程如下:
- 通过 Patch Partition 将输入的 RGB 图像分割成不重叠的 Patch,堆叠进 B 维度;
- 使用 Linear Embedding 将通道映射至 C;
- 紧接着使用两个连续的 Swin Transformer Block,将上述组合称为 Stage 1;
- 为了获得分层表示,通过 Patch Merging 对 Stage 1 的输出继续进行分块,并且同样会使用 Linear Layer 进行通道降维,再使用几个连续的 Swin Transformer Block;如此,便能构成更多的 Stage.
Shifted Window based Self-Attention

连续的 Swin Transformer 如上图所示,其主要流程如下:
-
对于第一个 Swin Transformer Block,会先对输入 使用 LayerNorm,然后使用 W-MSA(Window based Self-Attention),并且使用残差连接得到 ,可以写成如下形式:
-
接下来使用 LN、MLP(两层、GELU 激活函数)和残差连接的得到最终输出 ,可以写成如下形式:
-
对于接下来的 Swin Transformer Block,会将其 W-MSA 替换成 SW-MSA(Shifted Window based Self-Attention),可写成如下形式:
至此便完成了连续的 Swin Transformer Block 的构建,由于需要将这两种组合起来达到信息交换的目的,因此层数的设置应为偶数。
Self-attention in non-overlapped windows
为了实现线性的时间复杂度,提出在 Window(窗口)中进行建模,窗口以非重叠的方式均匀地划分图像,这种方式在局部窗口中进行 Patch 的关系建模,计算注意力时,会将 Patch 展品与标准多头自注意力的时间复杂度对比如下:
其中输入包含 个 Patch.
由于 M 是固定的,所有 W-MSA 对输入图像尺寸的复杂度呈线性。
Shifted window partitioning in successive blocks
虽然W-MSA 解决了 MSA 时间复杂度随输入二次增长的问题,但是不同窗口间没有信息交流,这显然会限制模型的建模能力。
为了保持高效的同时进行有效建模,提出了 Shifted Window:

通过控制不同框的大小,实现上一层不同 Window 之间的信息交流,但是这样较难实现,并且 Window 的数量会从 增加到 ,并且某些 Window 的大小会小于 ,因此提出了一种更简单的方法来实现这个功能:

将原有的窗口以 M/2 的大小进行偏移,将多出的部分移动到相对的位置,这样就实现了不同 Window 之间的信息交流,偏移 M/2 是为了保证足够大范围进行信息交流。
需要注意的一点是,实际计算的过程中会使用 Mask,将上图右侧移动过来的位置给盖住,原因是这部分计算注意力没有意义,因为我们只需要计算邻近部分的注意力,cyclic shift 会造成图片对侧的部分移动到另一侧,显然我们是不需要计算对侧部分的
**虽然使用了 cyclic shift 依旧保持 window 的数量不增加,实际上我们仍然需要在 9 个 window 里计算自注意力 **

如上图所示,注意 mask 的尺寸,其直接作用在计算好的注意力矩阵上
-
window0 所包含的“像素”在空间位置上都是相邻的,因此 mask 为 0;
-
而 window1 所包含的 1 和 2 显然不是相邻的,而我们只需要他们计算各自区域内的注意力(相当于只在 1 和 2 内进行计算分别计算),window1 的 mask 如此之混乱,原因就在于注意力的计算过程,涉及到矩阵相乘,从 知乎 上偷了张图:
举个例子,window1 中的如上图最左边的所示,大小为 7×7,在映射到 KQV 之前,我们需要先将其 flatten,再经过 kqv 的映射,就成了上图中的 Q 和 K,当 Q 和 K 矩阵相乘时,我们实际上只需要黄色的地方与黄色的地方计算,粉色的地方与粉色的地方进行计算;矩阵相乘之后我们会得到右边的结果,橙色的地方即是我们不需要,因为我们使用 mask 与其相加,消除影响。
Relative position bias
添加了相对位置偏置 ,其描述每个 Window 相对于其它 Window 的相对位置,注意力公式可以写成:
该相对位置偏置可以学习,相较于绝对位置嵌入拥有更好的性能,主要原因是使用相同的绝对位置权重时,对每个 Query 来说,其拥有不同的窗口范围,这些不同的窗口相对于实际的共享窗口是有偏差的,消融实验中发现,某些任务绝对位置嵌入会甚至会造成精度降低。

由于每个轴上的相对位置的取值范围都是 ,于是生成一个小的偏置矩阵 ,相对位置偏置 从 中采样而来。
Patch merging
由于 transformer 中并没有分辨率的说法,作者设计了一个“下采样”的方法——Patch merging,具体实现方式是 CNN 中空间到深度的变换(Pixel Shuffle),将空间信息堆叠进通道中,相当于变相扩大了 Window 的大小

重要代码分析
Window operation
window_partition:
将输入图像分割成 大小的 patch,并堆叠进 Batch 维度。
def window_partition(x, window_size):
"""
Args:
x: (B, H, W, C)
window_size (int): window size
Returns:
windows: (num_windows*B, window_size, window_size, C)
"""
B, H, W, C = x.shape
x = x.view(B, H // window_size, window_size,
W // window_size, window_size, C)
windows = x.permute(0, 1, 3, 2, 4, 5).contiguous(
).view(-1, window_size, window_size, C)
return windows
window_reverse:
恢复,用于残差连接之前。
def window_reverse(windows, window_size, H, W):
"""
Args:
windows: (num_windows*B, window_size, window_size, C)
window_size (int): Window size
H (int): Height of image
W (int): Width of image
Returns:
x: (B, H, W, C)
"""
B = int(windows.shape[0] / (H * W / window_size / window_size))
x = windows.view(B, H // window_size, W // window_size,
window_size, window_size, -1)
x = x.permute(0, 1, 3, 2, 4, 5).contiguous().view(B, H, W, -1)
return x
WindowAttention
该部分代码为 W-MSA 和 SW-MSA,具体切换依赖于输入数据和 mask,该模块只计算负责 Window 内的自注意力。
class WindowAttention(nn.Module):
r""" Window based multi-head self attention (W-MSA) module with relative position bias.
It supports both of shifted and non-shifted window.
Args:
dim (int): Number of input channels.
window_size (tuple[int]): The height and width of the window.
num_heads (int): Number of attention heads.
qkv_bias (bool, optional): If True, add a learnable bias to query, key, value. Default: True
qk_scale (float | None, optional): Override default qk scale of head_dim ** -0.5 if set
attn_drop (float, optional): Dropout ratio of attention weight. Default: 0.0
proj_drop (float, optional): Dropout ratio of output. Default: 0.0
"""
def __init__(self, dim, window_size, num_heads, qkv_bias=True, qk_scale=None, attn_drop=0., proj_drop=0.):
super().__init__()
self.dim = dim
self.window_size = window_size # Wh, Ww
self.num_heads = num_heads
head_dim = dim // num_heads
self.scale = qk_scale or head_dim ** -0.5
# define a parameter table of relative position bias
self.relative_position_bias_table = nn.Parameter(
torch.zeros((2 * window_size[0] - 1) * (2 * window_size[1] - 1), num_heads)) # 2*Wh-1 * 2*Ww-1, nH
# get pair-wise relative position index for each token inside the window
coords_h = torch.arange(self.window_size[0])
coords_w = torch.arange(self.window_size[1])
coords = torch.stack(torch.meshgrid([coords_h, coords_w])) # 2, Wh, Ww
coords_flatten = torch.flatten(coords, 1) # 2, Wh*Ww
relative_coords = coords_flatten[:, :, None] - \
coords_flatten[:, None, :] # 2, Wh*Ww, Wh*Ww
relative_coords = relative_coords.permute(
1, 2, 0).contiguous() # Wh*Ww, Wh*Ww, 2
relative_coords[:, :, 0] += self.window_size[0] - \
1 # shift to start from 0
relative_coords[:, :, 1] += self.window_size[1] - 1
relative_coords[:, :, 0] *= 2 * self.window_size[1] - 1
relative_position_index = relative_coords.sum(-1) # Wh*Ww, Wh*Ww
self.register_buffer("relative_position_index",
relative_position_index)
self.qkv = nn.Linear(dim, dim * 3, bias=qkv_bias)#使用一个线性层生成QKV,使用切片分开
self.attn_drop = nn.Dropout(attn_drop)
self.proj = nn.Linear(dim, dim)#对输出进行映射
self.proj_drop = nn.Dropout(proj_drop)
trunc_normal_(self.relative_position_bias_table, std=.02)
self.softmax = nn.Softmax(dim=-1)
def forward(self, x, mask=None):
"""
Args:
x: input features with shape of (num_windows*B, N, C)
mask: (0/-inf) mask with shape of (num_windows, Wh*Ww, Wh*Ww) or None
"""
B_, N, C = x.shape
qkv = self.qkv(x).reshape(B_, N, 3, self.num_heads, C //
self.num_heads).permute(2, 0, 3, 1, 4)
# make torchscript happy (cannot use tensor as tuple)
q, k, v = qkv[0], qkv[1], qkv[2]
q = q * self.scale #意义不明
attn = (q @ k.transpose(-2, -1))
relative_position_bias = self.relative_position_bias_table[self.relative_position_index.view(-1)].view(
self.window_size[0] * self.window_size[1], self.window_size[0] * self.window_size[1], -1) # Wh*Ww,Wh*Ww,nH
relative_position_bias = relative_position_bias.permute(
2, 0, 1).contiguous() # nH, Wh*Ww, Wh*Ww
attn = attn + relative_position_bias.unsqueeze(0)
if mask is not None: #当使用SW-MSA时,会使用mask
nW = mask.shape[0]
attn = attn.view(B_ // nW, nW, self.num_heads, N,
N) + mask.unsqueeze(1).unsqueeze(0)
attn = attn.view(-1, self.num_heads, N, N)
attn = self.softmax(attn)
else:
attn = self.softmax(attn)
attn = self.attn_drop(attn)
x = (attn @ v).transpose(1, 2).reshape(B_, N, C)
x = self.proj(x)
x = self.proj_drop(x)
return x
def extra_repr(self) -> str:
return f'dim={self.dim}, window_size={self.window_size}, num_heads={self.num_heads}'
def flops(self, N):
# calculate flops for 1 window with token length of N
flops = 0
# qkv = self.qkv(x)
flops += N * self.dim * 3 * self.dim
# attn = (q @ k.transpose(-2, -1))
flops += self.num_heads * N * (self.dim // self.num_heads) * N
# x = (attn @ v)
flops += self.num_heads * N * N * (self.dim // self.num_heads)
# x = self.proj(x)
flops += N * self.dim * self.dim
return flops
SwinTransformer
SwinTransformerBlock:
class SwinTransformerBlock(nn.Module):
r""" Swin Transformer Block.
Args:
dim (int): Number of input channels.
input_resolution (tuple[int]): Input resulotion.
num_heads (int): Number of attention heads.
window_size (int): Window size.
shift_size (int): Shift size for SW-MSA.
mlp_ratio (float): Ratio of mlp hidden dim to embedding dim.
qkv_bias (bool, optional): If True, add a learnable bias to query, key, value. Default: True
qk_scale (float | None, optional): Override default qk scale of head_dim ** -0.5 if set.
drop (float, optional): Dropout rate. Default: 0.0
attn_drop (float, optional): Attention dropout rate. Default: 0.0
drop_path (float, optional): Stochastic depth rate. Default: 0.0
act_layer (nn.Module, optional): Activation layer. Default: nn.GELU
norm_layer (nn.Module, optional): Normalization layer. Default: nn.LayerNorm
"""
def __init__(self, dim, input_resolution, num_heads, window_size=7, shift_size=0,
mlp_ratio=4., qkv_bias=True, qk_scale=None, drop=0., attn_drop=0., drop_path=0.,
act_layer=nn.GELU, norm_layer=nn.LayerNorm):
super().__init__()
self.dim = dim
self.input_resolution = input_resolution
self.num_heads = num_heads
self.window_size = window_size
self.shift_size = shift_size
self.mlp_ratio = mlp_ratio
if min(self.input_resolution) <= self.window_size:
# if window size is larger than input resolution, we don't partition windows
self.shift_size = 0
self.window_size = min(self.input_resolution)
assert 0 <= self.shift_size < self.window_size, "shift_size must in 0-window_size"
self.norm1 = norm_layer(dim)
self.attn = WindowAttention(
dim, window_size=to_2tuple(self.window_size), num_heads=num_heads,
qkv_bias=qkv_bias, qk_scale=qk_scale, attn_drop=attn_drop, proj_drop=drop)
self.drop_path = DropPath(
drop_path) if drop_path > 0. else nn.Identity()
self.norm2 = norm_layer(dim)
mlp_hidden_dim = int(dim * mlp_ratio)
self.mlp = Mlp(in_features=dim, hidden_features=mlp_hidden_dim,
act_layer=act_layer, drop=drop)
if self.shift_size > 0: # shift_size表示使用SW-MSA
# calculate attention mask for SW-MSA
H, W = self.input_resolution
img_mask = torch.zeros((1, H, W, 1)) # 1 H W 1
h_slices = (slice(0, -self.window_size),
slice(-self.window_size, -self.shift_size),
slice(-self.shift_size, None))
w_slices = (slice(0, -self.window_size),
slice(-self.window_size, -self.shift_size),
slice(-self.shift_size, None))
cnt = 0
for h in h_slices:
for w in w_slices:
img_mask[:, h, w, :] = cnt
cnt += 1
# nW, window_size, window_size, 1
mask_windows = window_partition(img_mask, self.window_size)
mask_windows = mask_windows.view(-1,
self.window_size * self.window_size)
attn_mask = mask_windows.unsqueeze(1) - mask_windows.unsqueeze(2)
attn_mask = attn_mask.masked_fill(
attn_mask != 0, float(-100.0)).masked_fill(attn_mask == 0, float(0.0))
else:
attn_mask = None
self.register_buffer("attn_mask", attn_mask)
def forward(self, x):
H, W = self.input_resolution
B, L, C = x.shape
assert L == H * W, "input feature has wrong size"
shortcut = x
x = self.norm1(x)
x = x.view(B, H, W, C)
# cyclic shift
if self.shift_size > 0:
shifted_x = torch.roll(
x, shifts=(-self.shift_size, -self.shift_size), dims=(1, 2))
else:
shifted_x = x
# partition windows
# nW*B, window_size, window_size, C
x_windows = window_partition(shifted_x, self.window_size)
# nW*B, window_size*window_size, C
x_windows = x_windows.view(-1, self.window_size * self.window_size, C)
# W-MSA/SW-MSA
# nW*B, window_size*window_size, C
attn_windows = self.attn(x_windows, mask=self.attn_mask)
# merge windows
attn_windows = attn_windows.view(-1,
self.window_size, self.window_size, C)
shifted_x = window_reverse(
attn_windows, self.window_size, H, W) # B H' W' C ,还原
# reverse cyclic shift
if self.shift_size > 0:
x = torch.roll(shifted_x, shifts=(
self.shift_size, self.shift_size), dims=(1, 2)) # 使用torch.roll实现shift
else:
x = shifted_x
x = x.view(B, H * W, C)
# FFN
x = shortcut + self.drop_path(x)
x = x + self.drop_path(self.mlp(self.norm2(x)))
return x
stage: