Rethinking the Inception Architecture for Computer Vision

Christian Szegedy, Vincent Vanhoucke, Sergey Ioffe, Jonathon Shlens, Zbigniew Wojna

Introduction

Since the 2012 ImageNet competition winning entry by Krizhevsky et al , their network “AlexNet” has been successfully applied to a larger variety of computer vision tasks, for example to object-detection , segmentation , human pose estimation , video classification , object tracking , and superresolution .

These successes spurred a new line of research that focused on finding higher performing convolutional neural networks. Starting in 2014, the quality of network architectures significantly improved by utilizing deeper and wider networks. VGGNet and GoogLeNet yielded similarly high performance in the 2014 ILSVRC classification challenge. One interesting observation was that gains in the classification performance tend to transfer to significant quality gains in a wide variety of application domains. This means that architectural improvements in deep convolutional architecture can be utilized for improving performance for most other computer vision tasks that are increasingly reliant on high quality, learned visual features. Also, improvements in the network quality resulted in new application domains for convolutional networks in cases where AlexNet features could not compete with hand engineered, crafted solutions, e.g. proposal generation in detection.

Although VGGNet has the compelling feature of architectural simplicity, this comes at a high cost: evaluating the network requires a lot of computation. On the other hand, the Inception architecture of GoogLeNet was also designed to perform well even under strict constraints on memory and computational budget. For example, GoogleNet employed only 5 million parameters, which represented a 12×12\times reduction with respect to its predecessor AlexNet, which used 6060 million parameters. Furthermore, VGGNet employed about 3x more parameters than AlexNet.

The computational cost of Inception is also much lower than VGGNet or its higher performing successors . This has made it feasible to utilize Inception networks in big-data scenarios, , where huge amount of data needed to be processed at reasonable cost or scenarios where memory or computational capacity is inherently limited, for example in mobile vision settings. It is certainly possible to mitigate parts of these issues by applying specialized solutions to target memory use , or by optimizing the execution of certain operations via computational tricks . However, these methods add extra complexity. Furthermore, these methods could be applied to optimize the Inception architecture as well, widening the efficiency gap again.

Still, the complexity of the Inception architecture makes it more difficult to make changes to the network. If the architecture is scaled up naively, large parts of the computational gains can be immediately lost. Also, does not provide a clear description about the contributing factors that lead to the various design decisions of the GoogLeNet architecture. This makes it much harder to adapt it to new use-cases while maintaining its efficiency. For example, if it is deemed necessary to increase the capacity of some Inception-style model, the simple transformation of just doubling the number of all filter bank sizes will lead to a 4x increase in both computational cost and number of parameters. This might prove prohibitive or unreasonable in a lot of practical scenarios, especially if the associated gains are modest. In this paper, we start with describing a few general principles and optimization ideas that that proved to be useful for scaling up convolution networks in efficient ways. Although our principles are not limited to Inception-type networks, they are easier to observe in that context as the generic structure of the Inception style building blocks is flexible enough to incorporate those constraints naturally. This is enabled by the generous use of dimensional reduction and parallel structures of the Inception modules which allows for mitigating the impact of structural changes on nearby components. Still, one needs to be cautious about doing so, as some guiding principles should be observed to maintain high quality of the models.

General Design Principles

Here we will describe a few design principles based on large-scale experimentation with various architectural choices with convolutional networks. At this point, the utility of the principles below are speculative and additional future experimental evidence will be necessary to assess their accuracy and domain of validity. Still, grave deviations from these principles tended to result in deterioration in the quality of the networks and fixing situations where those deviations were detected resulted in improved architectures in general.

Avoid representational bottlenecks, especially early in the network. Feed-forward networks can be represented by an acyclic graph from the input layer(s) to the classifier or regressor. This defines a clear direction for the information flow. For any cut separating the inputs from the outputs, one can access the amount of information passing though the cut. One should avoid bottlenecks with extreme compression. In general the representation size should gently decrease from the inputs to the outputs before reaching the final representation used for the task at hand. Theoretically, information content can not be assessed merely by the dimensionality of the representation as it discards important factors like correlation structure; the dimensionality merely provides a rough estimate of information content.

Higher dimensional representations are easier to process locally within a network. Increasing the activations per tile in a convolutional network allows for more disentangled features. The resulting networks will train faster.

Spatial aggregation can be done over lower dimensional embeddings without much or any loss in representational power. For example, before performing a more spread out (e.g. 3×33\times 3) convolution, one can reduce the dimension of the input representation before the spatial aggregation without expecting serious adverse effects. We hypothesize that the reason for that is the strong correlation between adjacent unit results in much less loss of information during dimension reduction, if the outputs are used in a spatial aggregation context. Given that these signals should be easily compressible, the dimension reduction even promotes faster learning.

Balance the width and depth of the network. Optimal performance of the network can be reached by balancing the number of filters per stage and the depth of the network. Increasing both the width and the depth of the network can contribute to higher quality networks. However, the optimal improvement for a constant amount of computation can be reached if both are increased in parallel. The computational budget should therefore be distributed in a balanced way between the depth and width of the network.

Although these principles might make sense, it is not straightforward to use them to improve the quality of networks out of box. The idea is to use them judiciously in ambiguous situations only.

Factorizing Convolutions with Large Filter Size

Much of the original gains of the GoogLeNet network arise from a very generous use of dimension reduction. This can be viewed as a special case of factorizing convolutions in a computationally efficient manner. Consider for example the case of a 1×11\times 1 convolutional layer followed by a 3×33\times 3 convolutional layer. In a vision network, it is expected that the outputs of near-by activations are highly correlated. Therefore, we can expect that their activations can be reduced before aggregation and that this should result in similarly expressive local representations.

Here we explore other ways of factorizing convolutions in various settings, especially in order to increase the computational efficiency of the solution. Since Inception networks are fully convolutional, each weight corresponds to one multiplication per activation. Therefore, any reduction in computational cost results in reduced number of parameters. This means that with suitable factorization, we can end up with more disentangled parameters and therefore with faster training. Also, we can use the computational and memory savings to increase the filter-bank sizes of our network while maintaining our ability to train each model replica on a single computer.

Convolutions with larger spatial filters (e.g. 5×55\times 5 or 7×77\times 7) tend to be disproportionally expensive in terms of computation. For example, a 5×55\times 5 convolution with nn filters over a grid with mm filters is 25/9 = 2.78 times more computationally expensive than a 3×33\times 3 convolution with the same number of filters. Of course, a 5×55\times 5 filter can capture dependencies between signals between activations of units further away in the earlier layers, so a reduction of the geometric size of the filters comes at a large cost of expressiveness. However, we can ask whether a 5×55\times 5 convolution could be replaced by a multi-layer network with less parameters with the same input size and output depth. If we zoom into the computation graph of the 5×55\times 5 convolution, we see that each output looks like a small fully-connected network sliding over 5×55\times 5 tiles over its input (see Figure 1).

Since we are constructing a vision network, it seems natural to exploit translation invariance again and replace the fully connected component by a two layer convolutional architecture: the first layer is a 3×33\times 3 convolution, the second is a fully connected layer on top of the 3×33\times 3 output grid of the first layer (see Figure 1). Sliding this small network over the input activation grid boils down to replacing the 5×55\times 5 convolution with two layers of 3×33\times 3 convolution (compare Figure 4 with 5).

This setup clearly reduces the parameter count by sharing the weights between adjacent tiles. To analyze the expected computational cost savings, we will make a few simplifying assumptions that apply for the typical situations: We can assume that n=αmn=\alpha m, that is that we want to change the number of activations/unit by a constant alpha factor. Since the 5×55\times 5 convolution is aggregating, α\alpha is typically slightly larger than one (around 1.5 in the case of GoogLeNet). Having a two layer replacement for the 5×55\times 5 layer, it seems reasonable to reach this expansion in two steps: increasing the number of filters by α\sqrt{\alpha} in both steps. In order to simplify our estimate by choosing α=1\alpha=1 (no expansion), If we would naivly slide a network without reusing the computation between neighboring grid tiles, we would increase the computational cost. sliding this network can be represented by two 3×33\times 3 convolutional layers which reuses the activations between adjacent tiles. This way, we end up with a net 9+925×\frac{9+9}{25}\times reduction of computation, resulting in a relative gain of 28%28\% by this factorization. The exact same saving holds for the parameter count as each parameter is used exactly once in the computation of the activation of each unit. Still, this setup raises two general questions: Does this replacement result in any loss of expressiveness? If our main goal is to factorize the linear part of the computation, would it not suggest to keep linear activations in the first layer? We have ran several control experiments (for example see figure 2) and using linear activation was always inferior to using rectified linear units in all stages of the factorization. We attribute this gain to the enhanced space of variations that the network can learn especially if we batch-normalize the output activations. One can see similar effects when using linear activations for the dimension reduction components.

2 Spatial Factorization into Asymmetric Convolutions

The above results suggest that convolutions with filters larger 3×33\times 3 a might not be generally useful as they can always be reduced into a sequence of 3×33\times 3 convolutional layers. Still we can ask the question whether one should factorize them into smaller, for example 2×22\times 2 convolutions. However, it turns out that one can do even better than 2×22\times 2 by using asymmetric convolutions, e.g. n×1n\times 1. For example using a 3×13\times 1 convolution followed by a 1×31\times 3 convolution is equivalent to sliding a two layer network with the same receptive field as in a 3×33\times 3 convolution (see figure 3). Still the two-layer solution is 33%33\% cheaper for the same number of output filters, if the number of input and output filters is equal.

By comparison, factorizing a 3×33\times 3 convolution into a two 2×22\times 2 convolution represents only a 11%11\% saving of computation.

In theory, we could go even further and argue that one can replace any n×nn\times n convolution by a 1×n1\times n convolution followed by a n×1n\times 1 convolution and the computational cost saving increases dramatically as nn grows (see figure 6). In practice, we have found that employing this factorization does not work well on early layers, but it gives very good results on medium grid-sizes (On m×mm\times m feature maps, where mm ranges between 1212 and 2020). On that level, very good results can be achieved by using 1×71\times 7 convolutions followed by 7×17\times 1 convolutions.

Utility of Auxiliary Classifiers

has introduced the notion of auxiliary classifiers to improve the convergence of very deep networks. The original motivation was to push useful gradients to the lower layers to make them immediately useful and improve the convergence during training by combating the vanishing gradient problem in very deep networks. Also Lee et al argues that auxiliary classifiers promote more stable learning and better convergence.

Interestingly, we found that auxiliary classifiers did not result in improved convergence early in the training: the training progression of network with and without side head looks virtually identical before both models reach high accuracy. Near the end of training, the network with the auxiliary branches starts to overtake the accuracy of the network without any auxiliary branch and reaches a slightly higher plateau.

Also used two side-heads at different stages in the network. The removal of the lower auxiliary branch did not have any adverse effect on the final quality of the network. Together with the earlier observation in the previous paragraph, this means that original the hypothesis of that these branches help evolving the low-level features is most likely misplaced. Instead, we argue that the auxiliary classifiers act as regularizer. This is supported by the fact that the main classifier of the network performs better if the side branch is batch-normalized or has a dropout layer. This also gives a weak supporting evidence for the conjecture that batch normalization acts as a regularizer.

Efficient Grid Size Reduction

Traditionally, convolutional networks used some pooling operation to decrease the grid size of the feature maps. In order to avoid a representational bottleneck, before applying maximum or average pooling the activation dimension of the network filters is expanded. For example, starting a d×dd\times d grid with kk filters, if we would like to arrive at a d2×d2\frac{d}{2}\times\frac{d}{2} grid with 2k2k filters, we first need to compute a stride-1 convolution with 2k2k filters and then apply an additional pooling step. This means that the overall computational cost is dominated by the expensive convolution on the larger grid using 2d2k22d^{2}k^{2} operations. One possibility would be to switch to pooling with convolution and therefore resulting in 2(d2)2k22(\frac{d}{2})^{2}k^{2} reducing the computational cost by a quarter. However, this creates a representational bottlenecks as the overall dimensionality of the representation drops to (d2)2k(\frac{d}{2})^{2}k resulting in less expressive networks (see Figure 9). Instead of doing so, we suggest another variant the reduces the computational cost even further while removing the representational bottleneck. (see Figure 10). We can use two parallel stride 2 blocks: PP and CC. PP is a pooling layer (either average or maximum pooling) the activation, both of them are stride 22 the filter banks of which are concatenated as in figure 10.

Inception-v2

Here we are connecting the dots from above and propose a new architecture with improved performance on the ILSVRC 2012 classification benchmark. The layout of our network is given in table 1. Note that we have factorized the traditional 7×77\times 7 convolution into three 3×33\times 3 convolutions based on the same ideas as described in section 3.1. For the Inception part of the network, we have 33 traditional inception modules at the 35×3535\times 35 with 288288 filters each. This is reduced to a 17×1717\times 17 grid with 768768 filters using the grid reduction technique described in section 5. This is is followed by 55 instances of the factorized inception modules as depicted in figure 5. This is reduced to a 8×8×12808\times 8\times 1280 grid with the grid reduction technique depicted in figure 10. At the coarsest 8×88\times 8 level, we have two Inception modules as depicted in figure 6, with a concatenated output filter bank size of 2048 for each tile. The detailed structure of the network, including the sizes of filter banks inside the Inception modules, is given in the supplementary material, given in the model.txt that is in the tar-file of this submission. However, we have observed that the quality of the network is relatively stable to variations as long as the principles from Section 2 are observed. Although our network is 4242 layers deep, our computation cost is only about 2.52.5 higher than that of GoogLeNet and it is still much more efficient than VGGNet.

Model Regularization via Label Smoothing

Here we propose a mechanism to regularize the classifier layer by estimating the marginalized effect of label-dropout during training.

We propose a mechanism for encouraging the model to be less confident. While this may not be desired if the goal is to maximize the log-likelihood of training labels, it does regularize the model and makes it more adaptable. The method is very simple. Consider a distribution over labels u(k)u(k), independent of the training example xx, and a smoothing parameter ϵ\epsilon. For a training example with ground-truth label yy, we replace the label distribution q(kx)=δk,yq(k|x)=\delta_{k,y} with

which is a mixture of the original ground-truth distribution q(kx)q(k|x) and the fixed distribution u(k)u(k), with weights 1ϵ1-\epsilon and ϵ\epsilon, respectively. This can be seen as the distribution of the label kk obtained as follows: first, set it to the ground-truth label k=yk=y; then, with probability ϵ\epsilon, replace kk with a sample drawn from the distribution u(k)u(k). We propose to use the prior distribution over labels as u(k)u(k). In our experiments, we used the uniform distribution u(k)=1/Ku(k)=1/K, so that

We refer to this change in ground-truth label distribution as label-smoothing regularization, or LSR.

Note that LSR achieves the desired goal of preventing the largest logit from becoming much larger than all others. Indeed, if this were to happen, then a single q(k)q(k) would approach 11 while all others would approach . This would result in a large cross-entropy with q(k)q^{\prime}(k) because, unlike q(k)=δk,yq(k)=\delta_{k,y}, all q(k)q^{\prime}(k) have a positive lower bound.

Another interpretation of LSR can be obtained by considering the cross entropy:

Thus, LSR is equivalent to replacing a single cross-entropy loss H(q,p)H(q,p) with a pair of such losses H(q,p)H(q,p) and H(u,p)H(u,p). The second loss penalizes the deviation of predicted label distribution pp from the prior uu, with the relative weight ϵ1ϵ\frac{\epsilon}{1-\epsilon}. Note that this deviation could be equivalently captured by the KL divergence, since H(u,p)=DKL(up)+H(u)H(u,p)=D_{KL}(u\|p)+H(u) and H(u)H(u) is fixed. When uu is the uniform distribution, H(u,p)H(u,p) is a measure of how dissimilar the predicted distribution pp is to uniform, which could also be measured (but not equivalently) by negative entropy H(p)-H(p); we have not experimented with this approach.

In our ImageNet experiments with K=1000K=1000 classes, we used u(k)=1/1000u(k)=1/1000 and ϵ=0.1\epsilon=0.1. For ILSVRC 2012, we have found a consistent improvement of about 0.2%0.2\% absolute both for top-11 error and the top-55 error (cf. Table 3).

Training Methodology

We have trained our networks with stochastic gradient utilizing the TensorFlow distributed machine learning system using 5050 replicas running each on a NVidia Kepler GPU with batch size 3232 for 100100 epochs. Our earlier experiments used momentum with a decay of 0.90.9, while our best models were achieved using RMSProp with decay of 0.90.9 and ϵ=1.0\epsilon=1.0. We used a learning rate of 0.0450.045, decayed every two epoch using an exponential rate of 0.940.94. In addition, gradient clipping with threshold 2.02.0 was found to be useful to stabilize the training. Model evaluations are performed using a running average of the parameters computed over time.

Performance on Lower Resolution Input

A typical use-case of vision networks is for the the post-classification of detection, for example in the Multibox context. This includes the analysis of a relative small patch of the image containing a single object with some context. The tasks is to decide whether the center part of the patch corresponds to some object and determine the class of the object if it does. The challenge is that objects tend to be relatively small and low-resolution. This raises the question of how to properly deal with lower resolution input.

The common wisdom is that models employing higher resolution receptive fields tend to result in significantly improved recognition performance. However it is important to distinguish between the effect of the increased resolution of the first layer receptive field and the effects of larger model capacitance and computation. If we just change the resolution of the input without further adjustment to the model, then we end up using computationally much cheaper models to solve more difficult tasks. Of course, it is natural, that these solutions loose out already because of the reduced computational effort. In order to make an accurate assessment, the model needs to analyze vague hints in order to be able to “hallucinate” the fine details. This is computationally costly. The question remains therefore: how much does higher input resolution helps if the computational effort is kept constant. One simple way to ensure constant effort is to reduce the strides of the first two layer in the case of lower resolution input, or by simply removing the first pooling layer of the network.

For this purpose we have performed the following three experiments:

299×299299\times 299 receptive field with stride 22 and maximum pooling after the first layer.

151×151151\times 151 receptive field with stride 11 and maximum pooling after the first layer.

79×7979\times 79 receptive field with stride 11 and without pooling after the first layer.

All three networks have almost identical computational cost. Although the third network is slightly cheaper, the cost of the pooling layer is marginal and (within 1%1\% of the total cost of the)network. In each case, the networks were trained until convergence and their quality was measured on the validation set of the ImageNet ILSVRC 2012 classification benchmark. The results can be seen in table 2. Although the lower-resolution networks take longer to train, the quality of the final result is quite close to that of their higher resolution counterparts.

However, if one would just naively reduce the network size according to the input resolution, then network would perform much more poorly. However this would an unfair comparison as we would are comparing a 16 times cheaper model on a more difficult task.

Also these results of table 2 suggest, one might consider using dedicated high-cost low resolution networks for smaller objects in the R-CNN context.

Experimental Results and Comparisons

Table 3 shows the experimental results about the recognition performance of our proposed architecture (Inception-v2) as described in Section 6. Each Inception-v2 line shows the result of the cumulative changes including the highlighted new modification plus all the earlier ones. Label Smoothing refers to method described in Section 7. Factorized 7×77\times 7 includes a change that factorizes the first 7×77\times 7 convolutional layer into a sequence of 3×33\times 3 convolutional layers. BN-auxiliary refers to the version in which the fully connected layer of the auxiliary classifier is also batch-normalized, not just the convolutions. We are referring to the model in last row of Table 3 as Inception-v3 and evaluate its performance in the multi-crop and ensemble settings.

All our evaluations are done on the 48238 non-blacklisted examples on the ILSVRC-2012 validation set, as suggested by . We have evaluated all the 50000 examples as well and the results were roughly 0.1% worse in top-5 error and around 0.2% in top-1 error. In the upcoming version of this paper, we will verify our ensemble result on the test set, but at the time of our last evaluation of BN-Inception in spring indicates that the test and validation set error tends to correlate very well.

Conclusions

We have provided several design principles to scale up convolutional networks and studied them in the context of the Inception architecture. This guidance can lead to high performance vision networks that have a relatively modest computation cost compared to simpler, more monolithic architectures. Our highest quality version of Inception-v3 reaches 21.2%21.2\%, top-11 and 5.6%5.6\% top-5 error for single crop evaluation on the ILSVR 2012 classification, setting a new state of the art. This is achieved with relatively modest (2.5×2.5\times) increase in computational cost compared to the network described in Ioffe et al . Still our solution uses much less computation than the best published results based on denser networks: our model outperforms the results of He et al – cutting the top-55 (top-11) error by 25%25\% (14%14\%) relative, respectively – while being six times cheaper computationally and using at least five times less parameters (estimated). Our ensemble of four Inception-v3 models reaches 3.5%3.5\% with multi-crop evaluation reaches 3.5%3.5\% top-55 error which represents an over 25%25\% reduction to the best published results and is almost half of the error of ILSVRC 2014 winining GoogLeNet ensemble.

We have also demonstrated that high quality results can be reached with receptive field resolution as low as 79×7979\times 79. This might prove to be helpful in systems for detecting relatively small objects. We have studied how factorizing convolutions and aggressive dimension reductions inside neural network can result in networks with relatively low computational cost while maintaining high quality. The combination of lower parameter count and additional regularization with batch-normalized auxiliary classifiers and label-smoothing allows for training high quality networks on relatively modest sized training sets.

References