Neural Network Encapsulation
A capsule is a collection of neurons which represents different variants of a pattern in the network. The routing scheme ensures only certain capsules which resemble lower counterparts in the higher layer should be activated. However, the computational complexity becomes a bottleneck for scaling up to larger networks, as lower capsules need to correspond to each and every higher capsule. To resolve this limitation, we approximate the routing process with two branches: a master branch which collects primary information from its direct contact in the lower layer and an aide branch that replenishes master based on pattern variants encoded in other lower capsules. Compared with previous iterative and unsupervised routing scheme, these two branches are communicated in a fast, supervised and one-time pass fashion. The complexity and runtime of the model are therefore decreased by a large margin. Motivated by the routing to make higher capsule have agreement with lower capsule, we extend the mechanism as a compensation for the rapid loss of information in nearby layers. We devise a feedback agreement unit to send back higher capsules as feedback. It could be regarded as an additional regularization to the network. The feedback agreement is achieved by comparing the optimal transport divergence between two distributions (lower and higher capsules). Such an add-on witnesses a unanimous gain in both capsule and vanilla networks. Our proposed EncapNet performs favorably better against previous state-of-the-arts on CIFAR10/100, SVHN and a subset of ImageNet.
Keywords:Network architecture design; capsule feature learning.
Convolutional neural networks (CNNs)  have been proved to be quite successful in modern deep learning architectures [2, 3, 4, 5] and achieved better performance in various computer vision tasks [6, 7, 8]. By tying the kernel weights in convolution, CNNs have the translation invariance property that can identify the same pattern irrespective of the spatial location. Each neuron in CNNs is a scalar and can detect different (low-level details or high-level regional semantics) patterns layer by layer. However, in order to detect the same pattern with various variants in viewpoint, rotation, shape, etc., we need to stack more layers, which tends to “memorize the dataset rather than generalize a solution” .
A capsule [10, 11] is a group of neurons whose output, in form of a vector instead of a scalar, represents various perspectives of an entity, such as pose, deformation, velocity, texture, object parts or regions, etc. It captures the existence of a feature and its variant. Not only does a capsule detect a pattern but also it is trained to learn the many variants of the pattern. This is what CNNs are incapable of. The concept of capsule provides a new perspective on feature learning via instance parameterization of entities (known as capsules) to encode different variants within a capsule structure, thus achieving the feature equivariance property111Equivariance is the detection of feature patterns that can transform to each other. and being robust to adversaries. Intuitively, the capsule detects a pattern (say a face) with a certain variant (it rotates 20 degree clockwise) rather than realizes that the pattern matches a variant in the higher layer.
One basic capsule layer consists of two steps: capsule mapping and agreement routing, which is depicted in Fig. 1(a). The input capsules are first mapped into the space of their higher counterparts via a transform matrix. Then the routing process involves all capsules between adjacent layers to communicate by the routing co-efficients; it ensures only certain lower capsules which resemble higher ones (in terms of cosine similarity) can pass on information and activate the higher counterparts. Such a scheme can be seen as a feature clustering and is optimized by coordinate descent through several iterations. However, the computational complexity in the first mapping step is the main bottleneck to apply the capsule idea in CNNs; lower capsules have to generate correspondence for every higher capsule (e.g., a typical choice  is 2048 capsules with 16 dimension, resulting in 8 million parameters in the transform matrix).
To tackle this drawback, we propose an alternative to estimate the original routing summation by introducing two branches: one is the master branch that serves as the primary source from the direct contact capsule in the lower layer; another is the aide branch that strives for searching other pattern variants along the channel and replenishes side information to master. These two branches are intertwined by their co-efficients so that feature patterns encoded in lower capsules could be fully leveraged and exchanged. Such a one-pass approximation is fast, light-weight and supervised, compared to the current iterative, short-lived and unsupervised routing scheme.
Furthermore, the routing effect in making higher capsule have agreement with lower capsule can be extended as a direct loss function. In deep neural networks, information is inevitably lost through stack of layers. To reduce the rapid loss of information in nearby layers, a loss function can be included to enforce that neurons or capsules in the higher layer can be used for reconstructing the counterparts in lower layers. Based on this motivation, we devise an agreement feedback unit which sends back higher capsules as a feedback signal to better supervise feature learning. This could be deemed as a regularization on network. Such a feedback agreement is achieved by measuring the distance between the two distributions using optimal transport (OT) divergence, namely the Sinkhorn loss. The OT metric (e.g., Wasserstein loss) is promised to be superior than other options to modeling data on general space. This add-on regularization is inserted during training and disposed of for inference. The agreement enforcement has witnessed a unanimous gain in both capsule and vanilla neural networks.
Altogether, bundled with the two mechanisms aforementioned, we (i) encapsulate the neural network in an approximate routing scheme with master/aide interaction, (ii) enforce the network’s regularization by an agreement feedback unit via optimal transport divergence. The proposed capsule network is denoted as EncapNet and performs superior against previous state-of-the-arts for image recognition tasks on CIFAR10/100, SVHN and a subset of ImageNet. The code and dataset are available https://github.com/hli2020/nn_capsulation.
2 CapNet: Agreement Routing Analysis
2.1 Preliminary: capsule formulation
Let denote the input and output capsules in a layer, where indicates the index of capsules. The dimension and the number of capsules at input and output are , respectively, i.e., . The first step is a mapping from lower capsules to higher counterparts: where is a transform matrix and we define the intermediate output as mapped activation (called prediction vector in ) from to . The second step is an agreement routing process to aggregate all lower capsules into higher ones. The mapped activation is multiplied by a routing coefficient through several iterations in an unsupervised manner: .
This is where the highlight of capsule idea resides in. It could be deemed as a voting process: the activation of higher capsules should be entirely dependent on the resemblance from the lower entities. Prevalent routing algorithms include the coordinate descent optimization  and the Gaussian mixture clustering via Expectation-Maximum (EM) , to which we refer as dynamic and EM routing, respectively. For dynamic routing, given , we have:
where is the softmax input to obtain ; is computed from via , i.e., . The update of the routing co-efficient is conducted in a coordinate descent manner which optimizes and alternatively. For EM routing, given , and the activation response of input capsules , we iteratively aggregate input capsules into Gaussian clusters:
where the mean of cluster is deemed as the output capsule . M-step generates the activation alongside the mean and std w.r.t. higher capsules; these variables are further fed into E-step to update the routing co-efficients . The output from a capsule layer is thereby obtained after iterating times.
2.2 Agreement routing analysis in CapNet
Effectiveness of the agreement routing. Figure 2 illustrates the training dynamics on routing between adjacent capsules as the network evolves. In essence, the routing process is a weighted average from all lower capsules to the higher entity (Eqn.(4)). Intuitively, given a sample which belongs to the -th class, the network tries to optimize capsule learning such that the length (existence probability) of in the final capsule layer should be the largest. This requires the magnitude of its lower counterparts who resemble capsule should occupy a majority and have a higher length compared to others that are dissimilar to . Take the top row of Dynamic case for instance. At the first epoch, the kernel weights are initialized with Gaussian hence most capsules are orthogonal to each other and have the same length. As training goes (epoch 20 and 80), the percentage and length of “blurring” capsules, whose cosine similarity is around zero, goes down and the distribution evolves into a polarization: the most similar and dissimilar capsules gradually take the majority and hold a higher length than other ’s. As training approaches termination (epoch 200), such a phenomenon is further polarized and the network is at a stable state where the most resembled and non-resembled capsules have a higher percentage and length than the others. The role of agreement routing is to adjust the magnitude and relevance from lower capsules to higher capsules, such that the activation of relevant higher counterparts could be appropriately turned on and the pattern information from lower capsules be passed on.
The analysis for EM routing draws a unanimous conclusion. The polarization phenomenon is further intensified (c.f. (h) vs (d) in Fig. (2)). The percentage of dissimilar capsules is lower (20% vs 37%) whilst the length of similar capsules is higher (0.02 vs 0.01): implying that EM is potentially a better routing solution than dynamic, which is also verified by (a) vs (b) in Table 4.
Moreover, it is observed that replacing scalar neurons in traditional CNNs with vector capsules and routing is effective, c.f. (a-b) vs (c) in Table 4. We adopt the same blob shape for each layer in vanilla CNNs for fair comparison. However, when we increase the parameters of CNNs to the same amount as that of CapNet, the former performs better in (d). Due to the inherent design, CapNet requires more parameters than the traditional CNNs, c.f. (a) vs (c) in Table 4 with around 152 Mb for CapNet vs 24 Mb for vanilla CNNs.
The capsule network is implemented in a group convolution fashion by [10, 11], which is depicted in Fig. 1(b). It is assumed that the vector capsules are placed in the same way as the scalar neurons in vanilla CNNs. The spatial capsules in a channel share the same transform kernel since they search for the same patterns in different locations. The channel capsules own different kernels as they represent various patterns encapsulated in a group of neurons.
Computational complexity in CapNet. From an engineering perspective, the original design for capsules in CNN structure (see Fig. 1(b)) is to save computation cost in the capsule mapping step; otherwise it would take more kernel parameters (assuming spatial size is 8) to fulfill the mapping step. However, the burden is not eased effectively since step one has to generate a mapping for each and every capsule in the subsequent layer. The output channel size of the transform kernel in Tab. 4(a-b) is 1,048,576 (). If we feed the network with a batch size of 128 (even smaller option, e.g., 32), OOM (out-of-memory) occurs due to the super-huge volume of the transform kernel. The subtle difference of parameter size between dynamic and EM is that additionally the latter has a larger convolutional output before the first capsule operation to generate activations; and it has a set of trainable parameters in the EM routing. Another impact to consider is the routing co-efficient matrix of size , the computation cost for this part is lightweight though and yet it takes longer runtime than traditional CNNs due to the routing iteration times to update , especially for EM method that involves two update alternations.
Inspired by the routing-by-agreement scheme to aggregate feature patterns in the network and bearing in mind that the current solution has a large computation complexity, we resort to some alternative that both inherits the spirit of routing-to-interact among capsules and implements such a process in a fast and accurate fashion. This leads to the proposed scheme stated below.
|method||param #||mem. size||runtime||CIFAR-10||MNIST|
|(a) CapNet, dynamic||151.24||3,961 (8)||0.444||14.28||0.37|
|(b) CapNet, EM||152.44||10,078 (4)||0.957||12.66||0.31|
|(c) vanilla CNN, same shape||24.44||1,652 (128)||0.026||14.43||0.38|
|(d) vanilla CNN, similar param||146.88||2,420 (128)||0.146||12.09||0.33|
|(e) EncapNet, master||25.76||1,433 (128)||0.039||13.87||0.31|
|(f) EncapNet, master/aide||60.68||1,755 (128)||0.061||11.93||0.25|
3 EncapNet: Neural Network Encapsulation
3.1 Approximate routing with master/aide interaction
Recall that higher capsules are generated according to the voting co-efficient across all entities (capsules) in the lower layer:
Eqn. (4) can be grouped into two parts: one is a main mapping that directly receives knowledge from its lower counterpart , whose spatial location is the same as ’s; another is a side mapping which sums up all the remaining lower capsules, whose spatial location is different from ’s. Hence the original unsupervised and short-lived routing process can be approximated in a supervised manner (see Fig. 1(c)):
where is a location set along the channel dimension that directly maps lower capsules (there might be more than one) to higher ; is the complimentary set of that contains the remaining locations along the channel; is the spatial kernel size; altogether indicates the location set of all contributing lower capsules to create a higher capsule. Formally, we define and in Eqn. (6) as the master and aide activation, respectively, with their co-efficients denoted as and .
The master branch looks for the same pattern in two consecutive layers and thus only sees a window from its direct lower capsule. The aide branch, on the other hand, serves as a side unit to replenish information from capsules located in other channels. The convolution kernels in both branches use the spatial locality: kernels only attend to a small neighborhood of size and on the input capsule to generate the intermediate activation and . The master and aide activations in these two branches are communicated by their co-efficients in an interactive manner. Co-efficient is the output of group convolution; the input source is from both and , leveraging information encoded in capsules from both the master and aide branches.
After the interaction as shown in Fig. 1(c), we append the batch normalization , rectified non-linearity unit  and squash operations at the end. These connectivities are not shown in the figure for brevity. To this end, we have encapsulated one layer of the neural network with each neuron being replaced by a capsule, where interaction among them is achieved by the master/aide scheme, and denote the whole pipeline as the capConv layer. An encapsulated module is illustrated in Fig. 3(a), where several capConv layers are cascaded with a skip connection. There are two types of capConv. Type I is to increase the dimension of capsules across modules and merge spatially-distributed capsules. The kernel size in the master branch is set to be 3 in this type. Type II is to increase the depth of the module for a length of ; the dimension of capsule is unchanged; nor does the number of spatial capsules. The kernel size for the master branch in this type is set to be 1. The capFC block is a per-capsule-dimension operation of the fully-connected layer as in standard neural network. Table 2 gives an example of the proposed network, called EncapNet.
Comparison to CapNet. Compared to the heavy computation of generating a huge number of mappings for each higher capsule in CapNet, our design only requires two mappings in the master and aide branch. The computational complexity is reduced by a large margin: the kernel size in the transform matrix in the first step is times fewer and the routing scheme in the second step is times fewer ( being the spatial size of feature map). Take the previous setting in Table 4 for instance, our design leads to 1024 and 256 times fewer parameters than the original 8,388,608 and 4,194,304 parameters in these two steps. To this end, we replace the unsupervised, iterative routing process [10, 11] with a supervised, one-pass master/aide scheme. Compared with (a-b) in Table 4, our proposed method (e-f) has fewer parameters, less runtime, and better performance. It is also observed that the side information from the aide branch is a necessity to replenish the master branch, with baseline error 13.87% decreasing to 11.93% on CIFAR-10, c.f. (e) vs. (f) in Table 4.
3.2 Network regularization by feedback agreement
Motivated by the agreement routing where higher capsules should be activated if there is a good ‘agreement’ with lower counterparts, we include a loss that requires the higher layer to be able to recover the lower layer. The influence of such a constraint (loss) is used during training and removed during inference.
To put the intuition aforementioned in math notation, let and be a sample in space and , respectively, where are sample indices. Consider a set of observations, e.g. capsules at lower layer, , we design a loss which enforces samples on space as input (e.g., capsules at higher layer) can be mapped to on space through a differentiable function , i.e., . The data distribution, denoted as , for the generated set of samples should be as much close as the distribution for . In summary, our goal is to find that minimizes a certain loss or distance between two distributions 333In some literature, i.e., [14, 15], it is called the probability measure and commonly denoted as or ; a coupling is the joint distribution (measure). We use distribution or measure interchangeably in the following context. is the set of probability distributions over a metric space .: .
In this paper, we opt for an optimal transport (OT) metric to measure the distance. The OT metric between two joint probability distributions supported on two metric spaces is defined as the solution of the linear program :
where is a coupling; is the set of couplings that consists of joint distributions over the product space with marginals . Our formulation skips some mathematic notations; details are provided in [15, 16]. Intuitively, implies how much “mass” must be transported from to in order to transform the distribution into ; is the “ground cost” to move a unit mass from to . As is well known, Eqn. (7) becomes the p-Wasserstein distance (or loss, divergence) between probability measures when is equipped with a distance and , for some exponent .
Note that the expectation in Eqn. (7) is used for mini-batches of size . In our case, and are equal to the training batch size. Since both input measures are discrete for the indices and (capsules in the network), the coupling can be treated as a non-negative matrix , namely , where represents the Dirac unit mass distribution at point . Rephrasing the continuous case of Eqn. (7) into a discrete version, we have the desired OT loss:
where satisfies . indicates the Frobenius dot-product for two matrices and . Now the problem boils down to computing given some ground cost . We adopt the Sinkhorn algorithm  in an iterative manner, which is promised to have a differentiable loss function . Starting with , Sinkhorn iterates read :
where the Gibbs kernel is defined as ; is a control factor. For a given budget of iterations, we have:
which serves as a proxy for the OT coupling. Equipped with the computation of and having some form of cost in hand, we can minimize the optimal transport divergence along with other loss in the network.
In practice, we introduce a bias fix to the original OT distance in Eqn. (8), namely the Sinkhorn divergence . Given two sets of samples and accordingly distributions , the revision is defined as:
where is the module index. By tuning in from to , the Sinkhorn divergence has the property of taking the best of both OT (non-flat geometry) and MMD  (high-dimensional rigidity) loss, which we find in experiments improves performance.
The overall workflow to calculate a Sinkhorn divergence444The term Sinkhorn used in this paper is two-folds: one is to indicate the computation of via a Sinkhorn iterates; another is to imply the revised OT divergence. is depicted in Fig. 3(b). Note that our ultimate goal of applying OT loss is to make feature learning in the mainstream (blue blocks) better aligned across capsules in the network. It is added during training and abominated for inference. Therefore the design for Sinkhorn divergence has two principles: light-weighted and capsule-minded. Sub-networks and should increase as minimal parameters to the model as possible; the generator should be encapsulated to match the data structure. Note that the Sinkhorn divergence is optimized to minimize loss w.r.t. both , instead of the practice in [15, 19, 14] via an adversarial manner.
Discussions. (i) There are alternatives besides the OT metric for , e.g., the Kullback-Leibler (KL) divergence, which is defined as or Jenson-Shannon (JS) divergence. In , it is observed that these distances are not sensible when learning distributions supported by low dimensional manifolds on . Often the model manifold and the “true” distribution’s support often have a non-negligible intersection, implying that and are non-existent or infinite in some cases. In comparison, the optimal transport loss is continuous and differentiable on under mild assumptions nonetheless. (ii) Our design of feedback agreement unit is not limited to the capsule framework. Its effectiveness on vanilla CNNs is also verified by experimental results in Section 5.1.
Design choices in OT divergence. We use a deconvolutional version of the capConv block as the mapping function for reconstructing lower layer neurons from higher layer neurons. Before feeding into the cost function , samples from two distributions are passed into a feature extractor . The extractor is modeled by a vanilla neural network and can be regarded as a dimensionality reduction of onto a lower-dimension space. There are many options to design the cost function , such as cosine distance or norm. Moreover, it is found in experiments that if the gradient flow in the Sinkhorn iterates process is ignored as does in , the result gets slightly better. Remind that is dependent on (so does ); hence the whole OT unit can be trained in the standard optimizers (such as Adam ).
Overall loss function. The final loss of EncapNet is a weighted combination from both the Sinkhorn divergence across modules and the marginal loss  for capsule in the classification task: , where is the ground truth and class capsule outputs of the capFC layer, respectively; is a hyper-parameter to negotiate between these two losses (set to be 10).
4 Related Work
Capsule network. Wang et al.  formulated the routing process as an optimization problem that minimizes the clustering-like loss and a KL regularization term. They proposed a more general way to regularize the objective function, which shares similar spirit as the agglomerative fuzzy k-means algorithm . Shahroudnejad et al.  explained the capsule network inherently constructs a relevance path, by way of dynamic routing in an unsupervised way, to eliminate the need for a backward process. When a group of capsules agree for a parent one, they construct a part-whole relationship which can be considered as a relevance path. A variant capsule network  is proposed where capsule’s activation is obtained based on the eigenvalue of a decomposed voting matrix. Such a spectral perspective witnesses a faster convergence and a better result than the EM routing  on a learning-to-diagnose problem.
Attention vs. routing. In , Mnih et al. proposed a recurrent module to extract information by adaptively selecting a sequence of regions and to only attend the selected locations. DasNet  allows the network to iteratively focus attention on convolution filters via feedback connections from higher layer to lower ones. The network generates an observation vector, which is used by a deterministic policy to choose an action, and accordingly changes the weights of feature maps for better classifying objects. Vaswani et al.  formulated a multi-head attention for machine translation task where attention coefficients are calculated and parameterized by a compatibility function. Attention models aforementioned tries to learn the attended weights from lower neurons to higher ones. The lower activations are weighted by the learned parameters in attention module to generate higher activations. However, the agreement routing scheme [10, 11] is a top-down solution: higher capsules should be activated if and only if the most similar lower counterparts have a large response. The routing co-efficients is obtained by recursively looking back at lower capsules and updated based on the resemblance. Our approximate routing can be deemed as a bottom-up approach which shares similar spirit as attention models.
The experiments are conducted on CIFAR-10/100 , SVHN  and a large-scale dataset called “h-ImageNet”. We construct the fourth one as a subset of the ILSVRC 2012 classification database . It consists of 200 hard classes whose top-1 accuracy, based on the prediction output of the ResNet-18  model is lower than other classes. The ResNet-18 baseline model on h-ImageNet has a 41.83% top-1 accuracy. The dataset has a collection of 255725, 17101 images for training and validation, compared with CIFAR’s 50000 for training and 10000 for test. We manually crop the object with some padding for each image (if the bounding box is not provided) since the original image has too much background and might be too large (over 1500 pixels); after the pre-processing, each image size is around 50 to 500, compared with CIFAR’s input. “h-ImageNet” is proposed for fast verifying ML algorithms on a large-scale dataset which shares similar distribution as ImageNet.
Implementation details. The general settings are the same across datasets if not specified afterwards. Initial learning rate is set to 0.0001 and reduced by 90% with a schedule in epoch unit. Maximum epoch is 600. Adam  is used with momentum 0.9 and weight decay . Batch size is 128.
5.1 Ablative analysis
In this subsection we analyze the connectivity design in the encapsulated module and the many choices in the OT divergence unit. The depth of EncapNet and ResNet are the same 18 layers () for fair comparison. Their structures are depicted in Table 2. Remind that the comparison of capConv block with CapNet is reported in Table 4 and analyzed in Section 3.1.
|module||output size||cap dim.||EncapNet_v1||ResNet|
|-||, conv||, conv|
Design in capConv block. Table 3 (1-4) reports the different incoming sources of the co-efficients in the master and aide branches. Without using aide, case (1) serves as baseline where higher capsules are only generated from the master activation. Note that the 9.83% result is already superior than all cases in Table 4, due to the increase of network depth. Result show that obtaining from the activation in its own branch is better than obtaining from the other activations, c.f., cases (2) and (3). When the incoming source of co-efficient is from both branches, denoted as “maser/aide_v3” in (4), the pattern information from lower capsules is fully interacted by both master and aide branches; hence we achieve the best result of 7.41% when compared with cases (2) and (3). Table 3 (5-7) reports the result of adding skip connection based on case (4). It is observed that the skip connection used in both types of the capConv block make the network converge faster and get better performance (5.82%). Our final candidate model employs an additional OT unit with two Sinkhorn losses imposed on each module. One is the connectivity as shown in Fig. 3(a) where is half the size of ; another connectivity is the same as the skip connection path shown in the figure, where shares the same size with ; the “deconvolutional” generator in this connectivity has a stride of 1. It performances better (4.55%) than using one OT divergence alone (4.58%).
Network regularization design. Fig. 4 illustrates the training loss curve with and without OT (Sinkhorn) divergence. It is found that the performance gain is more evident for EncapNet than ResNet (21% vs 4% increase on two networks, respectively). Moreover, we testify the KL divergence option as a distance measurement to substitute the Sinkhorn divergence, shown as case (b) in Table 3. The error rate decreases for both model, suggesting that the idea of imposing regularization on the network training is effective; such an add-on is to keep feature patterns better aligned across layers. The subtlety is that the gain clearly differs when we replace Sinkhorn with KL in EncapNet while these two options barely matter in ResNet.
|capConv Design||error (%)|
|(1) master (baseline)||9.83|
|(8) two OTs||4.55|
|(a) capConv block (baseline)||5.82||8.03555ResNet-20 reported in  has a 8.75% error rate; some online third party implementation (link anonymised for submission) obtains 6.98%; we run the 18-layer model in PyTorch with settings stated in the context.|
|(d1) remove bias fix||4.71||-|
|(d2) do BP in||4.77||-|
|(d3) no extractor||5.79||-|
|(d4) use vanilla||5.01||-|
|(d5) use in||4.90||-|
Furthermore, we conduct a series of experiments (d-*) to prove the rationale of the Sinkhorn divergence design in Section 3.2. Without the bias fix, the result is inferior since it does not leverage both OT and MMD divergences (case d1); if we back-propagate the gradient in the path, the error rate slightly increases; the role of feature extractor is to down-sample both inputs to the same shape on a lower dimension for the subsequent pipeline to process. If we remove this functionality and directly compare the raw inputs using cosine distance, the error increases by a large margin to 5.79%, compared with baseline 5.82%; if we adopt norm to measure the distance between raw inputs, loss will not converge (not shown in Table 3). This verifies the urgent necessity of having a feature extractor; if the generator recovering from employs a standard CNN, the performance is inferior (5.01%) than the capsule version of the generator since data flows in form of capsules in the network; finally if we adopt norm to calculate after the feature extractor, the performance degrades as well.
5.2 Comparison to state-of-the-arts
As shown in Table 4, (a) on CIFAR-10/100 and SVHN, we achieve a better performance of 3.10%, 24.01% and 1.52 % compared to previous entires. The multi-crop test is a key factor to further enhance the result, which is widely used by other methods as well. (b) on h-ImageNet, v1 is the 18-layer structure and has a reasonable top-1 accuracy of 51.77%. We further increase the depth of EncapNet (known as v2) by stacking more capConv blocks, making a depth of 101 to compare with the ResNet-101 model. To ease runtime complexity due to the master/aide intertwined communication, we replace some blocks in the shallow layers with master alone. v3 has a larger input size. Moreover, we have the ultimate version of EncapNet with data augmentation (v3) and obtain an error rate of 40.05%, compared with the runner-up WRN  42.51%. Training on h-ImageNet roughly takes 2.9 days with 8 GPUs and batch size 256. (c) we have some preliminary results on the ILSVRC-CLS (complete-ImageNet) dataset, which are reported in terms of the top-5 error in Table 4.
|EncapNet||3.10 (3.56 0.12)||24.18||1.52 (1.87 0.11)||EncapNet_v3||42.76|
|ResNet ||6.43||-||-||WRN ||42.51|
|ELU ||6.55||24.28||-||ResNet-101 ||44.13|
|Batch NIN ||6.75||28.86||1.81||VGG ||55.76|
|Rec-CNN ||7.09||31.75||1.77||GoogleNet ||60.18|
|dasNet ||9.22||33.78||-||GoogleNet ||7.89|
|Maxout ||9.35||38.57||2.47||VGG ||8.43|
|AlexNet ||11.00||-||-||ResNet-101 ||6.21|
In this paper, we analyze the role of routing-by-agreement to aggregate feature clusters in the capsule network. To lighten the computational load in the original framework, we devise an approximate routing scheme with master-aide interaction. The proposed alternative is light-weight, supervised and one-time pass during training. The twisted interaction ensures that the approximation can make best out of lower capsules to activate higher capsules. Motivated by the routing process to make capsules better aligned across layers, we send back higher capsules as feedback signal to better supervise the learning across capsules. Such a network regularization is achieved by minimizing the distance of two distributions using optimal transport divergence during training. This regularization is also found to be effective for vanilla CNNs.
Acknowledgment. We would like to thank Jonathan Hui for the wonderful blog on capsule research, Gabriel PeyrÃ© and Yu Liu for helpful discussions in early stage. H. Li and X. Guo are funded by Hong Kong Ph.D. Fellowship Scheme.
-  LeCun, Y., Bottou, L., Bengio, Y., Haffner, P.: Gradient-based learning applied to document recognition. In: Proceedings of the IEEE. Volume 86. (1998) 2278–2324
-  Krizhevsky, A., Sutskever, I., Hinton, G.E.: Imagenet classification with deep convolutional neural networks. In: NIPS. (2012) 1106–1114
-  Simonyan, K., Zisserman, A.: Very deep convolutional networks for large-scale image recognition. In: ICLR. (2015)
-  Szegedy, C., Liu, W., Jia, Y., Sermanet, P., Reed, S., Anguelov, D., Erhan, D., Vanhoucke, V., Rabinovich, A.: Going deeper with convolutions. In: CVPR. (2015)
-  He, K., Zhang, X., Ren, S., Sun, J.: Deep residual learning for image recognition. In: CVPR. (2016)
-  Li, H., Liu, Y., Ouyang, W., Wang, X.: Zoom out-and-in network with map attention decision for region proposal and object detection. In: International Journal of Computer Vision (IJCV). (2018)
-  Li, H., Liu, Y., Zhang, X., An, Z., Wang, J., Chen, Y., Tong, J.: Do we really need more training data for object localization. In: IEEE International Conference on Image Processing. (2017)
-  Chi, Z., Li, H., Huchuan, Yang, M.H.: Dual deep network for visual tracking. IEEE Trans. on Image Processing (2017)
-  Hui, J.: Understanding Matrix capsules with EM Routing. https://jhui.github.io/2017/11/14/Matrix-Capsules-with-EM-routing-Capsule-Network (2017) Accessed: 2018-03-10.
-  Sabour, S., Frosst, N., Hinton, G.: Dynamic routing between capsules. In: NIPS. (2017)
-  Hinton, G.E., Sabour, S., Frosst, N.: Matrix capsules with EM routing. In: ICLR. (2018)
-  Ioffe, S., Szegedy, C.: Batch normalization: Accelerating deep network training by reducing internal covariate shift. In: ICML. (2015)
-  Nair, V., Hinton, G.E.: Rectified linear units improve restricted boltzmann machines. In: ICML. (2010) 807–814
-  Arjovsky, M., Chintala, S., Bottou, L.: Wasserstein GAN. arXiv preprint: 1701.07875 (2017)
-  Genevay, A., PeyrÃ©, G., Cuturi, M.: Learning generative models with sinkhorn divergences. arXiv preprint: 1706.00292 (2017)
-  Cuturi, M.: Sinkhorn distances: Lightspeed computation of optimal transport. NIPS (2013)
-  Sinkhorn, R.: A relationship between arbitrary positive matrices and doubly stochastic matrices. Ann. Math. Statist. (2) (06) 876–879
-  Gretton, A., Borgwardt, K., Rasch, M.J., Scholkopf, B., Smola, A.J.: A kernel method for the two-sample problem. NIPS (2007)
-  Salimans, T., Zhang, H., Radford, A., Metaxas, D.: Improving GANs using optimal transport. In: International Conference on Learning Representations. (2018)
-  Kingma, D.P., Ba, J.: Adam: A method for stochastic optimization. In: ICLR. (2015)
-  Wang, D., Liu, Q.: An optimization view on dynamic routing between capsules. In: Submit to ICLR workshop. (2018)
-  Li, M.J., Ng, M.K., ming Cheung, Y., Huang, J.Z.: Agglomerative fuzzy k-means clustering algorithm with selection of number of clusters. IEEE Transactions on Knowledge and Data Engineering 20 (2008) 1519–1534
-  Shahroudnejad, A., Mohammadi, A., Plataniotis, K.N.: Improved explainability of capsule networks: Relevance path by agreement. In: arXiv preprint:1802.10204. (2018)
-  Bahadori, M.T.: Spectral capsule networks. In: ICLR workshop. (2018)
-  Mnih, V., Heess, N., Graves, A., Kavukcuoglu, K.: Recurrent models of visual attention. In: NIPS. (2014)
-  Stollenga, M.F., Masci, J., Gomez, F., Schmidhuber, J.: Deep networks with internal selective attention through feedback connections. In: NIPS. (2014) 3545–3553
-  Vaswani, A., Shazeer, N., Parmar, N., Uszkoreit, J., Jones, L., Gomez, A.N., Kaiser, L., Polosukhin, I.: Attention is all you need. arXiv preprint: 1706.03762 (2017)
-  Krizhevsky, A., Hinton, G.: Learning multiple layers of features from tiny images. In: Technical Report. (2009)
-  Netzer, Y., Wang, T., Coates, A., Bissacco, A., Wu, B., Ng, A.Y.: Reading digits in natural images with unsupervised feature learning. In: NIPS workshop. (2011)
-  Russakovsky, O., Deng, J., Su, H., Krause, J., Satheesh, S., Ma, S., Huang, Z., Karpathy, A., Khosla, A., Bernstein, M., Berg, A.C., Fei-Fei, L.: ImageNet Large Scale Visual Recognition Challenge. International Journal of Computer Vision (IJCV) 115(3) (2015) 211–252
-  Zagoruyko, S., Komodakis, N.: Wide residual networks. In: BMVC. (2016)
-  Mishkin, D., Matas, J.: All you need is a good init. arXiv preprint:1511.06422 (2015)
-  Snoek, J., Rippel, O., Swersky, K., Kiros, R., Satish, N., Sundaram, N., Patwary, M., Prabhat, M., Adams, R.: Scalable bayesian optimization using deep neural networks. In: ICML. (2015)
-  Clevert, D.A., Unterthiner, T., Hochreiter, S.: Fast and accurate deep network learning by exponential linear units. arXiv preprint: 1511.07289 (2015)
-  Chang, J.R., Chen, Y.S.: Batch-normalized maxout network in network. In: arXiv preprint:1511.02583. (2015)
-  Liang, M., Hu, X.: Recurrent convolutional neural network for object recognition. In: The IEEE Conference on Computer Vision and Pattern Recognition (CVPR). (June 2015)
-  Agostinelli, F., Hoffman, M., Sadowski, P., Baldi, P.: Learning activation functions to improve deep neural networks. In: ICLR workshop. (2015)
-  Lee, C.Y., Xie, S., Gallagher, P., Zhang, Z., Tu, Z.: Deeply-supervised nets. arXiv preprint: 1409.5185 (2014)
-  Lin, M., Chen, Q., Yan, S.: Network in network. In: ICLR. (2014)
-  Goodfellow, I.J., Warde-farley, D., Mirza, M., Courville, A., Bengio, Y.: Maxout networks. In: ICML. (2013)
-  Frogner, C., Zhang, C., Mobahi, H., Araya-Polo, M., Poggio, T.: Learning with a wasserstein loss. NIPS (2015)
Appendix A More Details and Results in CapNet
A.i Configurations in Section 2.2
All the experiments in Table 4 and Figure 2 employ a shallow network of six layers. The first four layers across vanilla CNNs, CapNet and EncapNet use the same recipe: one standard convolution plus a BN-ReLU follow-up. The stride and padding in convolution is 3 and 1, respectively; these four modules are for increasing the channel number and down-sampling the feature map. After these modules, the activation map is of size .
For CapNet (case (a) and (b)), both the fifth and sixth layer are capsule-minded and implemented by group convolution, as shown in Fig. 1(b). The transform matrix (kernel in convolution) in the fifth layer is super huge (of size ) since the mapping step involves transforming capsule to each and every in the subsequent layer. After routing (dynamic or EM), the output has a shape of and is rearranged into the blob shape in a spatial-channel distributed manner () for the next layer. The budget of routing iteration is set to .
There are two subtleties when we implement the routing scheme. For dynamic routing, we find the loss does not converge if we conduct the softmax operation of co-efficients along the dimension (as suggested in ); instead we operate softmax along the dimension. For EM routing, it is observed that adding a normalization operation before data is fed into the sigmoid function to generate activation is crucial to make the network stable and faster towards convergence, since the std acquired after the M-step is small. Having them through a log function results in most values being into the range, say over 100; this will further generate all 1’s after sigmoid. If we add normalization to set the range of input data before sigmoid around zero, the output will be diversified in .
For the compared vanilla CNNs (case (c) and (d)), the fifth module is a conv-BN-ReLU cascade with average pooling to make spatial size of feature maps to be ; the last module is a fully connected layer. The network for case (d), “similar param”, is achieved by increasing the channel number such that the model size amounts to a similar level as case (a) and (b).
For the compared EncapNet (case (e) and (f)), the fifth layer is a capConv layer with approximate routing and master/aide interaction, which is proposed in the main paper. The sixth layer is a capFC layer with outputs being the same as CapNet’s ().
A.ii Ablative study on CapNet
In preliminary experiments, we conduct some ablative analysis on the routing scheme and different loss choices in CapNet, which is shown in Table 5.
|routing iter.||CIFAR-10||routing co-efficient||CIFAR-10|
First, the number of routing iteration. We find the error would get lower if we conduct the routing algorithm via coordinate descent () than a simple averaging across all lower capsules (, implying there is no routing). As the number of iteration goes up, the performance gets slightly better; since it would cost more runtime, we fix to be 3 considering both efficiency and accuracy. Second, the initialization of the routing co-efficients. The zero initialization is suggested in ; we also investigate other options. The performance difference between random and learnable initialization of the co-efficients does not vary much. Hence we adopt the original zero manner of initialization. At last, we investigate different loss choices under various optimizers. In general, RMSProp and Adam are better in sense of training capsule networks than does SGD. The popular cross entropy loss is inferior than the marginal  and spread  loss.
Appendix B More Details in EncapNet
The capsule generator is composed of a convTranspose2d operation with grouping, followed by a BN-ReLU unit; the value of grouping is the dimension of input capsules. The kernel size in the convolution is 3; the stride could be 2 or 1, depending on whether the size of output activations matches the size of input. The generator appends a squash operation at the end. The feature extractor (also known as “critic”) consists of two consecutive convolutions with each followed by a BN-ReLU sequel. The kernel size, padding and stride of these convolutions are (3,1,2); the number of output channel in the first convolution is one forth of the input’s while the second convolution only has one output channel.
The control factor to compute the kernel is set to . As suggested in , a stronger regularization (larger ) leads to stable result and faster convergence; hence can be set small. The length of Sinkhorn iterates is set to .
Appendix C Comparison to State-of-the-arts
The mild data augmentation includes the multi-crop test and random cropping; and the strong augmentation additionally includes random horizontal flipping, color jittering and longer training. For random cropping, a patch of size 32 is randomly cropped out of a resized image of size 34 on CIFAR and SVHN; a patch of size 128 (224) is randomly cropped out of a resized image of size 156 (256) on h-ImageNet across all versions. The number in parentheses indicates the setting in v3 version. For color jittering, the input color channels are randomly disturbed in brightness, contrast, saturation and hue with .
The h-ImageNet database666For detailed statistics, please refer to the GitHub repo mentioned earlier. is proposed to train an ImageNet-like dataset in a manageable time. EncapNet_v3 takes around 2.5 days on 4 GPUs, c.f. the full ImageNet database trainig  being one week on 8 GPUs. However, the correlation between h-ImageNet and the original database is not investigated by time of submission; whether the effectiveness on h-ImageNet still holds on ILSVRC for most algorithms will be left as future work.