Skip to content

Commit

Permalink
revise pooling
Browse files Browse the repository at this point in the history
  • Loading branch information
astonzhang committed Aug 27, 2018
1 parent 0980d9d commit 514d7bc
Show file tree
Hide file tree
Showing 3 changed files with 30 additions and 22 deletions.
3 changes: 1 addition & 2 deletions chapter_convolutional-neural-networks/conv-layer.md
Expand Up @@ -107,8 +107,7 @@ Y
conv2d = nn.Conv2D(1, kernel_size=(1, 2))
conv2d.initialize()
# 二维卷积层使用 4 维输入输出,格式为(批量大小,通道数,高,宽),这里批量大小和通道数均为
# 1。
# 二维卷积层使用 4 维输入输出,格式为(批量,通道,高,宽),这里批量大小和通道数均为 1。
X = X.reshape((1, 1, 6, 8))
Y = Y.reshape((1, 1, 6, 7))
Expand Down
Expand Up @@ -34,7 +34,7 @@ def comp_conv2d(conv2d, X):
# (1,1)代表批量大小和通道数(后面章节将介绍)均为 1。
X = X.reshape((1, 1) + X.shape)
Y = conv2d(X)
# 我们不关心前两维:批量大小和通道数
# 我们不关心前两维:批量和通道
return Y.reshape(Y.shape[2:])
X = nd.random.uniform(shape=(8, 8))
Expand Down
47 changes: 28 additions & 19 deletions chapter_convolutional-neural-networks/pooling.md
@@ -1,18 +1,27 @@
# 池化层

回忆在[卷积层](./conv-layer.md)小节里介绍的图像物体边缘检测应用中,我们构造了卷积核来精确的找到像素变化的位置。例如如果输出`Y[i, j]=1`那么表示`X[i, j]``X[i, j+1]`数值不一样这可能意味着物体边缘通过这两个元素之间。但实际中图像里我们感兴趣的物体不会总出现在固定位置即使我们连续拍摄同一个物体也极有可能出现像素上的偏移。这样导致同一个边缘对应的输出可能出现在`Y`中不同位置,进而对后面的模式识别造成不便。
回忆在[二维卷积层](conv-layer.md)小节里介绍的图像物体边缘检测应用中,我们构造卷积核从而精确地找到了像素变化的位置。设任意二维数组`X``i``j`列的元素为`X[i, j]`。如果我们构造的卷积核输出`Y[i, j]=1`那么说明输入中`X[i, j]``X[i, j+1]`数值不一样这可能意味着物体边缘通过这两个元素之间。但实际中图像里我们感兴趣的物体不会总出现在固定位置即使我们连续拍摄同一个物体也极有可能出现像素上的偏移。这样导致同一个边缘对应的输出可能出现在卷积输出`Y`中的不同位置,进而对后面的模式识别造成不便。

这一节我们介绍池化层(pooling layer),它的提出是为了缓解卷积层对位置的过度敏感性。
这一节我们介绍池化(pooling)层,它的提出是为了缓解卷积层对位置的过度敏感性。

## 二维最大、平均池化层
## 二维最大池化层和平均池化层

池化层同卷积层一样每次对输入数据的一个固定形状窗口元素计算输出。不同于卷积层里计算输入和核的互相关性,池化层直接计算窗口内元素的最大值或者平均值。图5.6展示了$2\times 2$最大池化层,其输出的第一个元素是输入的左上$2\times 2$窗口里的四个元素的最大值。然后同卷积层一样依次向左或向下移动窗口来计算其余的输出。
同卷积层一样,池化层每次对输入数据的一个固定形状窗口(又称池化窗口)中的元素计算输出。不同于卷积层里计算输入和核的互相关性,池化层直接计算池化窗口内元素的最大值或者平均值。该运算也称最大池化或平均池化。在最大池化中,池化窗口从输入数组的最左和最上方开始,按从左往右、从上往下的顺序,依次在输入数组上滑动。当池化窗口滑动到某一位置时,窗口中的输入子数组的最大值即输出数组中相应位置的元素。图5.6中的输出数组高和宽分别为2,其中的四个元素由取最大值运算得出:

![$2\times 2$最大池化层。](../img/pooling.svg)
$$
\max(0,1,3,4)=4,\\
\max(1,2,4,5)=5,\\
\max(3,4,6,7)=7,\\
\max(4,5,7,8)=8.\\
$$

如果这个池化层的输入是来自于前面我们构造了卷积核的卷积层的输出。假设此卷积层输入是`X`,那么不管是`X[i, j]``X[i, j+1]`值不同,还是`X[i, j+1]``X[i, j+2]`不同,池化层输出均有`Y[i, j]=1`。换句话说,使用$2\times 2$最大池化层,只要卷积层识别的模式在高和宽上移动不超过一个元素,我们均可以将其检测出来。
![池化窗口形状为$2\times 2$的最大池化层。阴影部分为第一个输出元素及其计算所使用的输入元素:$\max(0,1,3,4)=4$](../img/pooling.svg)

池化层的前向计算实现在`pool2d`函数里。它跟[“卷积层”](conv-layer.md)一节里`corr2d`函数非常类似,唯一的区别是在计算`Y[h, w]`上。
平均池化的工作原理与最大池化类似。池化窗口形状为$p \times q$的池化层称为$p \times q$池化层。其中的池化计算叫做$p \times q$池化。

让我们再次回到本节开始提到的图像物体边缘检测例子。现在我们将卷积层的输出作为$2\times 2$最大池化的输入。设该卷积层输入是`X`、池化层输出为`Y`。无论是`X[i, j]``X[i, j+1]`值不同,还是`X[i, j+1]``X[i, j+2]`不同,池化层输出均有`Y[i, j]=1`。也就是说,使用$2\times 2$最大池化层,只要卷积层识别的模式在高和宽上移动不超过一个元素,我们依然可以将它检测出来。

下面我们把池化层的前向计算实现在`pool2d`函数里。它跟[“二维卷积层”](conv-layer.md)一节里`corr2d`函数非常类似,唯一的区别是在计算输出`Y`上。

```{.python .input n=11}
from mxnet import nd
Expand All @@ -30,29 +39,29 @@ def pool2d(X, pool_size, mode='max'):
return Y
```

构造图5.6中的数据来验证实现的正确性
我们可以构造图5.6中的输入数组`X`来验证二维最大池化层的输出

```{.python .input n=13}
X = nd.array([[0, 1, 2], [3, 4, 5], [6, 7, 8]])
pool2d(X, (2, 2))
```

同时我们试一试平均池化层
同时我们实验一下平均池化层

```{.python .input n=14}
pool2d(X, (2, 2), 'avg')
```

## 填充和步幅

同卷积层一样,池化层也可以填充输入高宽两侧的数据和调整窗口的移动步幅来改变输入大小。我们将通过`nn`模块里的二维最大池化层MaxPool2D来演示它的工作机制。我们先构造一个`(1, 1, 4, 4)`形状的输入数据,前两个维度分别是批量和通道。
同卷积层一样,池化层也可以在输入的高和宽两侧的填充并调整窗口的移动步幅来改变输出形状。池化层填充和步幅与卷积层填充和步幅的工作机制一样。我们将通过`nn`模块里的二维最大池化层MaxPool2D来演示池化层填充和步幅的工作机制。我们先构造一个形状为`(1, 1, 4, 4)`的输入数据,前两个维度分别是批量和通道。

```{.python .input n=15}
X = nd.arange(16).reshape((1, 1, 4, 4))
X
```

MaxPool2D类里默认步幅设置成跟池化窗大小一样。下面使用`(3, 3)`窗口,默认获得`(3, 3)`步幅
`MaxPool2D`类里默认步幅设置成跟池化窗口大小一样。下面使用形状为`(3, 3)`的池化窗口,默认获得形状为`(3, 3)`的步幅

```{.python .input n=16}
pool2d = nn.MaxPool2D(3)
Expand All @@ -67,7 +76,7 @@ pool2d = nn.MaxPool2D(3, padding=1, strides=2)
pool2d(X)
```

当然,我们也可以是非方形的窗口,并且指定各个方向上的填充和步幅
当然,我们也可以指定非正方形的池化窗口,并分别指定高和宽上的填充和步幅

```{.python .input n=8}
pool2d = nn.MaxPool2D((2, 3), padding=(1, 2), strides=(2, 3))
Expand All @@ -76,14 +85,14 @@ pool2d(X)

## 多通道

在处理多通道输入数据时,池化层对每个输入通道分别池化,而不是像卷积层那样混合输入通道。这个意味着池化层的输出通道跟输入通道数相同。下面我们将`X``X+1`在通道维度上合并来构造通道数为2输入
在处理多通道输入数据时,池化层对每个输入通道分别池化,而不是像卷积层那样对各通道的输入按通道相加。这意味着池化层的输出通道数跟输入通道数相同。下面我们将数组`X``X+1`在通道维上连结来构造通道数为2的输入

```{.python .input n=9}
X = nd.concat(X, X + 1, dim=1)
X
```

做池化后我们发现输出通道仍然是2,而且通道0的结果跟之前一致
池化后,我们发现输出通道数仍然是2

```{.python .input n=10}
pool2d = nn.MaxPool2D(3, padding=1, strides=2)
Expand All @@ -92,14 +101,14 @@ pool2d(X)

## 小结

- 池化层通过滑动窗口计算结果,其通常直接取输入窗口内元素的最大值或者平均值作为输出
- 池化层的一个主要作用是缓解卷积层对位置的敏感性。
* 池化层取池化窗口中输入元素的最大值或者平均值作为输出
* 池化层的一个主要作用是缓解卷积层对位置的敏感性。

## 练习

- 分析池化层的计算复杂度。假设输入大小为$c\times h\times w$,我们使用$p_h\times p_w$的池化窗,而且使用$(p_h, p_w)$填充和$(s_h, s_w)$步幅,那么这个池化层的前向计算需要多少操作
- 想一想最大池化层和平均池化层的区别主要在哪里
- 你觉得最小池化层这个想法怎么样
* 分析池化层的计算复杂度。假设输入形状为$c\times h\times w$,我们使用形状为$p_h\times p_w$的池化窗口,而且使用$(p_h, p_w)$填充和$(s_h, s_w)$步幅。这个池化层的前向计算复杂度有多大
* 想一想,最大池化层和平均池化层在作用上可能有哪些区别
* 你觉得最小池化层这个想法有没有意义

## 扫码直达[讨论区](https://discuss.gluon.ai/t/topic/6406)

Expand Down

0 comments on commit 514d7bc

Please sign in to comment.