Year: 2020

Abstract:
Group convolutional neural networks (G-CNNs) can be used to improve classical CNNs by equipping them with the geometric structure of groups. Central in the success of G-CNNs is the lifting of feature maps to higher dimensional disentangled representations in which data characteristics are effectively learned, geometric data-augmentations are made obsolete, and predictable behavior under geometric transformations (equivariance) is guaranteed via group theory. Currently, however, the practical implementations of G-CNNs are limited to either discrete groups (that leave the grid intact) or continuous compact groups such as rotations (that enable the use of Fourier theory). In this paper we lift these limitations and propose a modular framework for the design and implementation of G-CNNs for arbitrary Lie groups. In our approach the differential structure of Lie groups is used to expand convolution kernels in a generic basis of B-splines that is defined on the Lie algebra. This leads to a flexible framework that enables localized, atrous, and deformable convolutions in G-CNNs by means of respectively localized, sparse and non-uniform B-spline expansions. The impact and potential of our approach is studied on two benchmark datasets: cancer detection in histopathology slides in which rotation equivariance plays a key role and facial landmark localization in which scale equivariance is important. In both cases, GCNN architectures outperform their classical 2D counterparts and the added value of atrous and localized group convolutions is studied in detail.

heading: 1 INTRODUCTION
text: Group convolutional neural networks (G-CNNs) are a class of neural networks that are equipped with the geometry of groups. This enables them to profit from the structure and symmetries in signal data such as images (Cohen & Welling, 2016). A key feature of G-CNNs is that they are equivariant with respect to transformations described by the group, i.e., they guarantee predictable behavior under such transformations and are insensitive to both local and global transformations on the input data. Classical CNNs are a special case of G-CNNs that are equivariant to translations and, in contrast to unconstrained NNs, they make advantage of (and preserve) the basic structure of signal data throughout the network (LeCun et al., 1990). By considering larger groups (i.e. considering not just translation equivariance) additional geometric structure can be utilized in order to improve performance and data efficiency (see G-CNN literature in Sec. 2).
Part of the success of G-CNNs can be attributed to the lifting of feature maps to higher dimensional objects that are generated by matching kernels under a range of poses (transformations in the group). This leads to a disentanglement with respect to the pose and together with the group structure this enables a flexible way of learning high level representations in terms of low-level activated neurons observed in specific configurations, which we conceptually illustrate in Fig. 1. From a neuro-psychological viewpoint, this resembles a hierarchical composition from low- to high-level features akin to the recognition-by-components model by Biederman (1987), a viewpoint which is also adopted in work on capsule networks (Hinton et al., 2011; Sabour et al., 2017). In particular in (Lenssen et al., 2018) the group theoretical connection is made explicit with equivariant capsules that provide a sparse index/value representation of feature maps on groups (Gens & Domingos, 2014).
∗Work done at Eindhoven University of Techology
Representing low-level features via features maps on groups, as is done in G-CNNs, is also motivated by the findings of Hubel & Wiesel (1959) and Bosking et al. (1997) on the organization of orientation sensitive simple cells in the primary visual cortex V1. These findings are mathematically modeled by sub-Riemannian geometry on Lie groups (Petitot, 2003; Citti & Sarti, 2006; Duits et al., 2014) and led to effective algorithms in image analysis (Franken & Duits, 2009; Bekkers et al., 2015b; Favali et al., 2016; Duits et al., 2018; Baspinar, 2018). In recent work Montobbio et al. (2019) show that such advanced V1 modeling geometries emerge in specific CNN architectures and in Ecker et al. (2019) the relation between group structure and the organization of V1 is explicitly employed to effectively recover actual V1 neuronal activities from stimuli by means of G-CNNs.
G-CNNs are well motivated from both a mathematical point of view (Cohen et al., 2018a; Kondor & Trivedi, 2018) and neuro-psychological/neuro-mathematical point of view and their improvement over classical CNNs is convincingly demonstrated by the growing body of G-CNN literature (see Sec. 2). However, their practical implementations are limited to either discrete groups (that leave the grid intact) or continuous, (locally) compact, unimodular groups such as roto-translations (that enable the use of Fourier theory). In this paper we lift these limitations and propose a framework for the design and implementation of G-CNNs for arbitrary Lie groups.
The proposed approach for G-CNNs relies on a definition of B-splines on Lie groups which we use to expand and sample group convolution kernels. B-splines are piece-wise polynomials with local support and are classically defined on flat Euclidean spaces Rd. In this paper we generalize B-splines to Lie groups and formulate a definition using the differential structure of Lie groups in which B-splines are essentially defined on the (flat) vector space of the Lie algebra obtained by the logarithmic map, see Fig. 2. The result is a flexible framework for B-splines on arbitrary Lie groups and it enables the construction of G-CNNs with properties that cannot be achieved via traditional Fourier-type basis expansion methods. Such properties include localized, atrous, and deformable convolutions in G-CNNs by means of respectively localized, sparse and non-uniform B-splines.
Although concepts described in this paper apply to arbitrary Lie groups, we here concentrate on the analysis of input data that lives on Rd and consider G-CNNs for affine groups G = Rd o H
that are the semi-direct product of the translation group with a Lie group H . As such, only a few core definitions about the Lie group H (group product, inverse, Log, and action on Rd) need to be implemented in order to build full G-CNNs that are locally equivariant to the transformations in H .
The impact and potential of our approach is studied on two datasets in which respectively rotation and scale equivariance plays a key role: cancer detection in histopathology slides (PCam dataset) and facial landmark localization (CelebA dataset). In both cases G-CNNs out-perform their classical 2D counterparts and the added value of atrous and localized G-convolutions is studied in detail.

heading: 2 RELATED WORK
text: G-CNNs The introduction of G-CNNs to the machine learning community by Cohen & Welling (2016) led to a growing body of G-CNN literature that consistently demonstrates an improvement of G-CNNs over classical CNNs. It can be roughly divided into work on discrete G-CNNs (Cohen & Welling, 2016; Dieleman et al., 2016; Winkels & Cohen, 2018; Worrall & Brostow, 2018; Hoogeboom et al., 2018), regular continuous G-CNNs (Oyallon & Mallat, 2015; Bekkers et al., 2015a; 2018b; Weiler et al., 2018a; Zhou et al., 2017; Marcos et al., 2017) and steerable continuous G-CNNs (Cohen et al., 2018b; Worrall et al., 2017; Kondor & Trivedi, 2018; Thomas et al., 2018; Weiler et al., 2018a; Esteves et al., 2018a; Andrearczyk et al., 2019). Since 3D rotations can only be sampled in very restrictive ways (without destroying the group structure) the construction of 3D roto-translation G-CNNs is limited. In order to avoid having to sample all together, steerable (G-)CNNs can be used. These are specialized G-CNNs in which the kernels are expanded in circlar/spherical harmonics and computations take place using the basis coefficients only (Chirikjian & Kyatkin, 2000; Franken, 2008; Almsick, 2007; Skibbe & Reisert, 2017). The latter approach is however only possible for unimodular groups such as roto-translations.
Scale equivariance In this paper we experiment with scale-translation G-CNNs, which is the first direct application of G-CNNs to achieve equivariance beyond roto-translations. Scale equivariance is however addressed in several settings (Henriques & Vedaldi, 2017; Esteves et al., 2018b; Marcos et al., 2018; Tai et al., 2019; Worrall & Welling, 2019; Jaderberg et al., 2015; Li et al., 2019), of which (Worrall & Welling, 2019) is most related. There, scale-space theory and semi-group theory is used to construct scale equivariant layers that elegantly take care of moving band-limits due to rescaling. Although our work differs in several ways (e.g. non-learned lifting layer, discrete group convolutions via atrous kernels, semi-group theory), the first two layers of deep scale-space networks relate to our lifting layer by treating our B-splines as a superposition of dirac deltas transformed under the semi-group action of (Worrall & Welling, 2019), as we show in App. C.1. Li et al. (2019) achieve scale invariance by sharing weights among kernels with different dilation rates. Their approach can be considered a special case of our proposed B-spline G-CNNs with kernels that do not encode scale interactions. Related work by Tai et al. (2019) and Henriques & Vedaldi (2017) relies on the same Lie group principles as we do in this paper (the Log map) to construct convenient coordinate systems, such as log-polar coordinates Esteves et al. (2018b), to handle equivariance. Such methods are however generally not translation equivariant and do not deal with local symmetries as they act globally on feature maps, much like spatial transformer networks (Jaderberg et al., 2015).
B-splines and vector fields in deep learning The current work can be seen as a generalization of the B-spline based SE(2) CNNs of Bekkers et al. (2015a; 2018b), see Sec. 3.3. Closely related is also the work of Fey et al. (2018) in which B-splines are used to generalize CNNs to non-Euclidean
data (graphs). There it is proposed to perform convolution via B-spline kernels on Rd that take as inputs vectors u(i, j) ∈ Rd that relate any two points i, j ∈ G in the graph to each other. How u(i, j) is constructed is left as a design choice, however, in (Fey et al., 2018) this is typically done by embedding the graph in an Euclidean space where points relate via offset vectors. In our work on Lie G-CNNs, two points g, g′ ∈ G in the Lie group G relate via the logarithmic map u(g, g′) = Log g−1g′. Another related approach in which convolutions take place on manifolds in terms of ”offset vectors” is the work by Cohen et al. (2019). There, points relate via the exponential map with respect to gauge frames rather than the left-invariant vector fields as in this paper, see App. C.2.

heading: 3 LIE GROUP CNNS

heading: 3.1 PRELIMINARIES AND NOTATION
text: The following describes the essential tools required for deriving a generic framework G-CNNs. Although we treat groups mostly in an abstract setting, we here provide examples for the rototranslation group and refer to App. B for more details, explicit formula’s, and figures for several groups. Implementations are available at https://github.com/ebekkers/gsplinets.
Group A group is defined by a set G together with a binary operator ·, the group product, that satisfies the following axioms: Closure: For all h, g ∈ G we have h · g ∈ G; Identiy: There exists an identity element e; Inverse: for each g ∈ G there exists an inverse element g−1 ∈ G such that g−1 · g = g · g−1 = e; and Associativity: For each g, h, i ∈ G we have (g · h) · i = g · (h · i). Lie group and Lie algebra If furthermore the group has the structure of a differential manifold and the group product and inverse are smooth, it is called a Lie group. The differentiability of the group induces a notion of infinitesimal generators (see also the exponential map below), which are elements of the Lie algebra g. The Lie algebra consists of a vector space (of generators), that is typically identified with the tangent space Te(G) at the identity e, together with a bilinear operator called the Lie bracket. In this work the Lie bracket is not of interest and we simply say g = Te(G). Exponential and logarithmic map We expand vectors in g in a basis {Ai}ni=1 and write A =∑n i a
iAi ∈ g, with components a = (a1, . . . , an) ∈ Rn. This allows us to identify the Lie algebra with Rn. Lie groups come with a logarithmic map that maps elements from the typically non-flat manifolds G to the flat Euclidean vector space g via A = Log g (conversely g = ExpA), see Fig. 2.
Semi-direct product groups In this paper we specifically consider (affine) Lie groups of type G = Rd oH that are the semi-direct product of the translation group Rd with a Lie group H that acts on Rd. Let h x denote the action of h ∈ H on x ∈ Rd; it describes how a vector in Rd is transformed by h. Then the group product of Rd oH is given by
g1 · g2 = (x1, h1) · (x2, h2) = (x1 + h1 x2, h1 · h2), (1) with g1 = (x1, h1), g2 = (x2, h2) ∈ G, x1,x2 ∈ Rd and h1, h2 ∈ H . For example the special Euclidean motion group SE(2) is constructed by choosing H = SO(2), the group of 2× 2 rotation matrices with matrix multiplication as the group product. The group product of G is then given by
(x1,Rθ1) · (x2,Rθ2) = (x1 + Rθ1 .x2,Rθ1 .Rθ2), with x1,x2 ∈ R2 and Rθ1 ,Rθ2 ∈ SO(2) rotation matrices parameterized by a rotation angle θi, and in which rotations act on vectors in R2 simply by matrix vector multiplication.
Group representations We consider linear transformations LG→L2(X)g : L2(X) → L2(X) that transform functions (or feature maps) f ∈ L2(X) on some space X as representations of a group G if they share the group structure via
(LG→L2(X)g1 ◦ L G→L2(X) g2 f)(x) = (L G→L2(X) g1·g2 f)(x),
with ◦ denoting function composition. Thus, a concatenation of two such transformations, parameterized by g1 and g2, can be described by a single transformation parameterized by g1 · g2. For semi-direct product groups G = Rd oH such a representation can be split into
LG→L2(X)g = LR d→L2(X) x ◦ L H→L2(X) h . (2)
Equivariance An operator Φ : L2(X)→ L2(Y ) is said to be equivariant to G when it satisfies ∀g∈G : LG→L2(Y )g ◦ Φ = Φ ◦ LG→L2(X)g . (3)

heading: 3.2 GROUP CONVOLUTIONAL NEURAL NETWORKS
text: Equivariance of artificial neural networks (NNs) with respect to a groupG is a desirable property as it guarantees that no information is lost when applying transformations to the input, the information is merely shifted to different locations in the network. It turns out that if we want NNs to be equivariant, then our only option is to use layers whose linear operator is defined by group convolutions. We summarize this in Thm. 1. To get there, we start with the traditional definition of NN layers via
y = φ(Kwx+ b), with x ∈ X the input vector, Kw : X → Y a linear map parameterized by a weight vector w, with b ∈ Y a bias term, and φ a point-wise non-linearity. In classical neural networks X = RNx and Y = RNy are Euclidean vector spaces and the linear map Kw = RNy×Nx is a weight matrix. In this work we instead focus on structured data and consider feature maps on some domainX as functions f : X → RN , the space of which we denote with (L2(X))N . In this caseX = (L2(X))Nx and Y = (L2(Y ))Ny are the spaces of multi-channel feature maps, b ∈ RNy , and Kw : X → Y is a kernel operator. Now when we constrain the linear operator Kw to be equivariant under transformations in some group G we arrive at group convolutional neural networks. This is formalized in the following theorem on equivariant maps between homogeneous spaces (see (Duits & Burgeth, 2007; Kondor & Trivedi, 2018; Cohen et al., 2018a) for related statements). Theorem 1. Let operator K : L2(X) → L2(Y ) be linear and bounded, let X,Y be homogeneous spaces on which Lie group G act transitively, and dµX a Radon measure on X , then
1. K is a kernel operator, i.e., ∃k̃∈L1(Y×X) : (Kf)(y) = ∫ X k̃(y, x)f(x)dµX ,
2. under the equivariance constraint of Eq. (3) the map is defined by a one-argument kernel
k̃(y, x) = dµX(g
−1 y x)
dµX(x) k(g−1y x) (4)
for any gy ∈ G such that y = gy y0 for some fixed origin y0 ∈ Y ,
3. if Y ≡ G/H is the quotient of G with H = StabG(y0) = {g ∈ G|g y0 = y0} then the kernel is constrained via
∀h∈H ,∀x∈X : k(x) = dµX(g
−1 y x)
dµX(x) k(h−1 x), (5)
Proof. See App. A
Corollary 1. IfX = Rd is a homogeneous space of an affine Lie groupG = RdoH and dµX(x) = dx is the Lebesgue measure on Rd then the kernel front-factor simplifies to dµX(g
−1 x) dµX(x)
= 1| deth| with |deth| denoting the determinant of the matrix representation of h, for any g = (x, h) ∈ G. If X = G and dµX(x) is a Haar measure on G then dµX(g −1 x)
dµX(x) = 1.
In view of Thm. 1 we see that standard CNNs are a special case of G-CNNs that are equivariant to translations. In this case the domain of the feature maps X = Rd coincides with the space of translation vectors in the translation group G = (Rd,+). It is well known that if we want the networks to be translation and rotation equivariant (G = Rd o SO(d)), but stick to planar feature maps, then the kernels should be rotation invariant, which of course limits representation power. This constraint is due Thm. 1 item 3, since the domain of such features maps is a quotient group Rd ≡ G/{0} × SO(d) in which rotations (SO(d)) are factored out of the roto-translation group G = Rd o SO(d). Thus, in order to maximize representation power (without constrains on k) the feature maps should be lifted to the higher dimensional domain of the group itself (i.e. Y = G). We therefore propose to build G-CNNs with the following 3 types of layers (illustrated in Fig. 1):
• Lifting layer (X = Rd, Y = G): In this layer K is defined by lifting correlations (Kf)(g) = (k?̃f)(g) := 1| deth| ( LG→L2(R d) g k, f ) L2(Rd,dx) ,
with g = (x, h), which by splitting of the representation (Eq. (2)) can be written as
(Kf)(g) = (k?̃f)(g) = (kh ?Rd f)(x), (6)
with kh(x) = 1| deth| ( LH→L2(R d) h k ) (x) the transformed kernel. Lifting correlations thus match a kernel with the input feature map under all possible transformations in G.
• Group correlation layer (X = G, Y = G): In this caseK is defined by group correlations
(KF )(g) = (K ? F )(g) := ( LG→L2(G)g K,F ) L2(G,dµ) = ∫ G K(g−1g̃)F (g̃)dµ(g̃),
with dµ(g) a Haar measure on G. We can again split this cross-correlation into a transformation of K followed by a spatial cross-correlation via
(KF )(g) = (K ? F )(g) = (Kh ?Rd F )(x), (7)
with Kh(x̃, h̃) = K(h−1 x̃ , h−1 · h̃ ) the convolution kernel transformed by h ∈ H and in which we overload ?Rd to indicate cross-correlation on the Rd part of G = Rd oH . • Projection layer (X = G, Y = Rd): In this case K is a linear projection defined by
(KF )(x) = ∫ H F (x, h̃)dµ(h̃), (8)
where we simply integrate over H instead of using a kernel that would otherwise be constant over H and spatially isotropic with respect to H .
App. B provides explicit formula’s for these layers for several groups. E.g., in the SE(2) we get: SE(2)-lifting: (k?̃f)(x, θ) = ∫ R2 k(R−1θ .(x ′ − x))f(x′)dx′,
SE(2)-correlation: (K?F )(x, θ) = ∫ R2 ∫ S1 K(R−1θ .(x ′ − x), θ′ − θmod 2π)F (x′, θ′)dx′dθ′.

heading: 3.3 B-SPLINES ON LIE GROUPS
text: Central in our formulation of G-CNNs is the transformation of convolution kernels under the action of H as described above in Eqs. (6) and (7) in the continuous setting. For the implementation of G-CNNs the kernels and their transformations need to be sampled on a discrete grid. We expand on the idea’s in (Bekkers et al., 2015a; 2018b; Weiler et al., 2018b) to express the kernels in an analytic form which we can then sample under arbitrary transformations in G to perform the actual computations. In particular we generalize the approach of Bekkers et al. (2015a; 2018b) to expand group correlation kernels in a basis of shifted cardinal B-splines, which are localized polynomial functions on Rd with finite support. In (Bekkers et al., 2015a; 2018b), B-splines on Rd could be used to construct kernels on SE(2) by identifying the group with the space of positions and orientations and simply using periodic splines on the 1D orientation axis S1. However, in order to construct B-splines on arbitrary Lie groups, we need a generalization. In the following we propose a new definition of B-splines on Lie groups H which enables us to construct the kernels on Rd oH that are required in the G-correlations (Eq. (7)).
Definition 1 (Cardinal B-spline on Rn). The 1D cardinal B-Spline of degree n be is defined as Bn(x) := (
1[− 12 , 1 2 ] ∗(n) 1[− 12 , 12 ]
) (x), (9)
where ∗(n) denotes n-fold convolution of the indicator function 1[− 12 , 12 ]. The multi-variate cardinal B-spline on Rd, with coordinates x = (x0, . . . , xd)T ∈ Rd, is defined via the tensor product
BR d,n(x) := (Bn ⊗ · · · ⊗Bn)(x)︸ ︷︷ ︸
d times
= Bn(x0)B n(x1) . . . B n(xd). (10)
Cardinal B-splines are piece-wise polynomials and are localized on support [−n+12 , n+1 2 ]. Functions can be expanded in a basis of shifted cardinal B-splines, which we simply refer to as B-splines.
Definition 2 (B-splines on Rn). A B-spline is a function f : Rd → R expanded in a basis that consists of shifted and scaled copies of the cardinal B-spline
f(x) := N∑ i=1 ciB Rd,n ( x− xi sx ) , (11)
and is fully characterized by spline degree n, scale sx, set of centers {xi}Ni=1 with xi ∈ Rd and corresponding coefficients c = (c1, c2, . . . , cN )T ∈ RN . The B-spline is called uniform if the set of centers {xi}Ni=1 forms a uniform grid on Rd, in which the distance ‖xj−xi‖ between neighbouring centers xi,xj ∈ Rd is constant along each axis and equal to sx. Definition 3 (B-splines on Lie group H). A B-spline on H is a function f : H → R expanded in a basis that consists of shifted (by left multiplication) and scaled copies of the cardinal B-spline
f(h) := N∑ i=1 ciB Rd,n ( Log h−1i h sh ) , (12)
with h ∈ H and Log : H → h the logarithmic map on H . The B-spline is fully characterized by the spline degree n, scale sh, set of centers {hi}Ni=1 with hi ∈ H and corresponding coefficients c = (c1, c2, . . . , cN )
T ∈ RN . The spline is called uniform if the distance ‖Log h−1i hj‖ between neighbouring centers hi, hj ∈ H is constant.
Examples of B-splines on Lie groups H are given in Fig. 3. In this paper we choose to expand convolution kernels on G = Rd o H as the tensor product of B-splines on Rd and H respectively and obtain functions f : Rd ×H via
f(x, h) := N∑ i=1 ciB Rd,n ( x− xi sx ) BR d,n ( Log h−1i h sh ) . (13)
Note that that one could also directly define B-splines on G via (12), however, this splitting ensures we can use a regular Cartesion grid on the Rd part. In our experiments we use B-splines as in (13) and consider the coefficients c as trainable parameters and the centers (xi and/or hi) and scales (sx and/or sh) are fixed by design. Some design choices are the following (and illustrated in Fig. 4).
Global vs localized uniform B-splines The notion of a uniform B-spline globally coveringH exists only for a small set of Lie groups, e.g. for any 1D group and abelian groups, and it is not possible to construct uniform B-splines on Lie groups in general due to non-zero commutators. Nevertheless, we find that it is possible to construct approximately uniform B-splines on compact groups either by constructing a grid of centers {hi}Ni=1 on H that approximately uniformly covers H , e.g. by using a repulsion model in which ‖Log h−1i · hj‖ between any two grid points hi, hj ∈ H is maximized (as is done in Fig. 3), or by specifying a uniform localized grid on the lie algebra h and obtaining the centers via the exponential map. The latter approach is in fact possible for any Lie group and leads to a notion of localized convolution kernels that have a finite support on H , see Fig. 4.
Atrous B-splines Atrous convolutions, i.e. convolutions with sparse kernels defined by weights interleaved with zeros (Holschneider et al., 1990), are commonly used to increase the effective receptive field size and add a notion of scale to deep CNNs (Yu & Koltun, 2016; Chen et al., 2018). Atrous convolution kernels can be constructed with B-splines by fixing the scale factors sx and sh, e.g. to the grid size, and increasing the distance between the center points xi and hi.
Non-uniform/deformable B-splines In non-uniform B-splines the centers xi and hi do not necessarily need to lie on a regular grid. Then, deformable CNNs, first proposed by Dai et al. (2017), are obtained by treating the centers as trainable parameters. For B-spline CNNs on Rd of order n = 1 this in fact leads to the deformable convolution layers as defined in (Dai et al., 2017).
Modular design The design of G-correlation layers (Eqs. (6-8)) using B-spline kernels (Eqs. (11- 13)) results in a generic and modular construction of G-CNNs that are equivariant to Lie groups G and only requires a few group specific definitions (see examples in App. B): The group structure of H (group product · and inverse), the action of H on Rd (together with the group structure of H this automatically defines the structure of G = Rd oH), and the logarithmic map Log : H → h.

heading: 4 EXPERIMENTS

heading: 4.1 ROTO-TRANSLATION CNNS
text: Data The PatchCamelyon (PCam) dataset (Veeling et al., 2018) consists of 327,680 RGB patches taken from histopathologic scans of lymph node sections and is derived from Camelyon16 (Ehteshami Bejnordi et al., 2017). The patches are binary labeled for the presence of metastasis. The classification problem is truly rotation invariant as image features appear under arbitrary rotations at all levels of abstraction, e.g. from edges (low-level) to individual cells to the tissue (high-level).
Experiments G-CNNs ensure roto-translation equivariance both locally (low-level) and globally (high-level) and invariance is achieved by means of pooling. In our experiments we test the performance of roto-translation G-CNNs (with G = R2 o SO(2)) against a 2D baseline and investigate the effect of different choices (local, global, atrous) for defining the kernels on the SO(2)-part of the network, cf. Eq. (13) and Fig. 4. Each network has the same architecture (detailed in App. D) but the kernels are sampled with varying resolution on H = SO(2), denoted with Nh, and with varying resolution of the B-splines, which is achieved by varying sh and the number of basis functions onH , denoted withNk. Each network has approximately the same number of trainable weights. Each network is trained with data-augmentation (see App. D for details), but for reference we also compare to a 2D and a SE(2) model (Nk = Nh = 12) which are trained without 90◦ rotation augmentation.
The results are summarized in Fig. 5. Here Nh = 1 means the kernels are transformed for only one rotation, which coincides with standard 2D convolutions (our baseline). A result labeled ”dense” with Nh = 16 and Nk = 8 means the convolution kernels are rotated 16 times and the kernels are expanded in a B-spline basis with 8 basis functions to fully cover H . The label ”local” means the basis is localized with Nk basis functions with a spacing of sh = 2π16 between them, with sh equal to the grid resolution. Atrous kernels are spaced equidistantly on H and have sh 2πNh .
Results We generally observe that a finer sampling of SO(2) leads to better results up untilNh = 12 after which results slightly degrade. This is line with findings in (Bekkers et al., 2018a). The degradation after this point could be explained by overfitting; there is a limit on the resolution of the signal generated by rotating 5x5 convolution kernels; at some point the splines are described in more detail than the data and thus an unnecessary amount of coefficients are trained. One could still benefit from sampling coarse kernels (low Nk) on a fine grid (high Nh), e.g. compare the cases Nh > Nk for fixed Nk. This is in line with findings in (Weiler et al., 2018b) where a fixed circular harmonic basis is used. Generally, atrous kernels tend to outperform dense kernels as do the localized kernels in the low Nk regime. Finally, comparing the models with and without 90◦ augmentation show that such augmentations are crucial for the 2D model but hardly affect the SE(2) model. Moreover, the SE(2) model without outperforms the 2D model with augmentation. This confirms the theory: G-CNNs guarantee both local and global equivariance by construction, whereas with augmentations valuable network capacity is spend on learning (only) global invariance. The very modest drop in the SE(2) case may be due to the discretization of the network on a grid after which it is no longer purely equivariant but rather approximately, which may be compensated for via augmentations.

heading: 4.2 SCALE-TRANSLATION CNNS
text: Data The CelebA dataset (Liu et al., 2015) contains 202,599 RGB images of varying size of celebrities together with labels for attributes (hair color, glasses, hat, etc) and 5 annotated facial landmarks (2 eyes, 1 nose, 2 corners of the mouth). We reformatted the data as follows. All images are isotropically scaled to a maximum width or height of 128 and if necessary padded in the other dimension with zeros to obtain a size of 128x128. For each image we took the distance between the eyes as a reference for the size of the face and categorized each image into above and below average size. For each unique celebrity with at least 1 image per class, we randomly sampled 1 image per class. The final dataset consists of 17,548 images of 128x128 of 8,774 celebrities with faces at varying scales. Each image is labeled with 5 heatmaps constructed by sampling a Gaussian with standard deviation 1.5 centered around each landmark.
Experiments We train a scale-translation G-CNN (with G = R2 o R+) with different choices for kernels. The ”dense” networks have kernels defined over the whole discretization of H = R+ and thus consider interactions between features at all scales. The ”local” networks consider only
interaction between neighbouring scales via localized kernels (Nk = 3) or no scale interaction at all (Nk = 1). Either way, each G-CNN is a multi-scale network in which kernels are applied at a range of scales. We compared against a 2D baseline with fixed-scale kernels which we tested for several scales separately. In the G-CNNs,H is uniformly sampled (w.r.t. to the metric onH) on a fixed scale range, generating the discrete sets Hd = {e(i−1)sh}Nhi=1 with sh = 12 ln 2. Each G-CNN is sampled with the same resolution in H with sh, and each B-spline basis function is centered on the discrete grid (i.e. hi ∈ Hd). We note that the discretization of H is formally no longer a group as it is not closed, however, the group structure still applies locally. The result is that information may leak out of the domain in a similar way as happens spatially in standard zero-padded 2D CNNs (translational G-CNNs), in which the discretized domain of translations is also no longer (locally) compact. This information leak can be avoided by using localized kernels of size Nk = 1 along the H axis, as is also done in (Worrall & Welling, 2019; Li et al., 2019). The networks are trained without dataaugmentation, except for the 2D network, which for reference we train with and without random scale augmentations at train-time. Augmentations were done with a random scale-factor between 1 and 1.4. We found that scale augmentations beyond 1.4 did not improve results.
Results Fig. 5 summarizes the results. By testing our 2D baseline at several scales we observe that there is an optimal scale (h = 2) that gives a best trade off between the scale variations in the data. This set of experiments is also used to rule out the idea that G-CNNs outperform the 2D baseline simply because they have a larger effective receptive field size. For large scale ranges the G-CNNs start to outperform 2D CNNs as these networks consider both small and large scale features simultaneously (multi-scale behavior). Comparing different G-CNN kernel specifications we observe that enabling interactions between neighbouring scales, via localized kernels on H (”local”, Nh = 3), outperforms both all-scale interactions (”dense”) and no-scale interaction at all (Nh = 1). This finding is in line with those in (Worrall & Welling, 2019). Finally, although 2D CNNs moderately benefit from scale augmentations, it is not possible to achieve the performance of G-CNNs. It seems that a multi-scale approach (G-CNNs) is essential.

heading: 5 CONCLUSION
text: This paper presents a flexible framework for building G-CNNs for arbitrary Lie groups. The proposed B-spline basis functions, which are used to represent convolution kernels, have unique properties that cannot be achieved by classical Fourier based basis functions. Such properties include the construction of localized, atrous, and deformable convolution kernels. We experimentally demonstrated the added value of localized and atrous group convolutions on two different applications, considering two different groups. In particular in experiments with scale-translation G-CNNs, kernel localization was important. The B-spline basis functions can be considered as smooth pixels on Lie groups and they enable us to design G-CNNs using familiar notions from classical CNN design (localized, atrous, and deformable convolutions). Future work will focus on exploring these options further in new applications that could benefit from equivariance constraints, for which the tools now are available for a large class of transformation groups via the proposed Lie group B-splines.

heading: ACKNOWLEDGMENTS
text: Remco Duits (Eindhoven University of Technology (TU/e)) is gratefully acknowledged for his contributions to the formulation and proof of Thm. 1 and for helpful discussions on Lie groups. I thank Maxime Lafarge (TU/e) for advice on setting up the PCam experiments. This work is part of the research programme VENI with project number 17290, which is (partly) financed by the Dutch Research Council (NWO).

heading: A PROOF OF THEOREM 1
text: The following proofs the three sub-items of Thm. 1.
1. It follows from Dunford-Pettis Theorem, see e.g. (Arendt & Bukhvalov, 1994, Thm 1.3), (Kantorovich & Akilov, 1982, Ch 9, Thm 5), or (Duits, 2005, Thm 1), that if K is linear and bounded it is an integral operator.
2. The left-equivariance constraint then imposes bi-left-invariance of the kernel k̃ as follows, where ∀g∈G and ∀f∈L2(X):
(K ◦ LG→L2(X)g )(f) = (LG→L2(Y )g ◦ K)(f) ⇔∫ X k̃(y, x)f(g−1x)dx = ∫ X k̃(g−1y, x)f(x)dx in r.h.s. integral x←g−1x⇔∫
X
k̃(y, x)f(g−1x)dx = ∫ X
k̃(g−1y, g−1x)f(g−1x)d(g−1x)⇔∫ X k̃(y, x)f(g−1x)dx = ∫ X k̃(g−1y, g−1x)f(g−1x) 1| det g|dx. (14)
Since (14) should hold for all f ∈ L2(X) we obtain
∀g∈G : k̃(y, x) = 1| det g| k̃(g −1y, g−1x). (15)
Furthermore, sinceG acts transitively on Y we have that ∀y,y0∈Y ∃gy∈G such that y = gyy0 and thus
k̃(y, x) = k̃(gy y0, x) (15) = 1| det gy| k̃(y0, g −1 y x) =: 1 | det gy|k(g −1 y x)
for every gy ∈ G such that y = gy y0 with arbitrary fixed origin y0 ∈ Y . 3. Every homogeneous space Y of G can be identified with a quotient group G/H . Choose
an origin y0 ∈ Y s.t. ∀h∈H : h y0 = y0, i.e., H = StabG y0, then
k̃(y0, x) = k̃(h y0, x)⇔ k(x) = 1| deth|k(h −1 x).
We further remark that When Y ≡ G = G/{e}, with e ∈ G the identify element ofG, the symmetry constraint of Eq. (5) vanishes. Thus, in order to construct equivariant maps without constraints on the kernel the functions should be lifted to the group G.

heading: B EXAMPLES OF LIE GROUPS
text: In the following sub-sections some explicit examples of Lie groups H are given, together with their actions on Rd and the Log operators. The required tools for building B-spline based G-CNNs for Lie groups of the form G = Rd o H are then automatically derived from these core definitions. E.g., the action of H on a space X defines a left-regular representation on functions on X via
(LH→L2(X)g f)(x) = f(g−1 x).
When X = G is the group itself, the action equals the group product. The group structure of semidirect product groups G = RdoH is automatically derived from the action of H on Rd, see Eq. (1) and is in turn used to define the representations (see Eq. (2)). Some examples are given below.
B.1 TRANSLATION GROUP G = (Rd,+)
The group of translations is given by the space of translation vectors Rd with the group product and inverse given by
g · g′ = (x + x′) g−1 = (−x),
with g = (x), g′ = (x′) ∈ G with x,x′ ∈ Rd. The identity element is e = (0). The left-regular representation on d-dimensional functions f ∈ L2(Rd) produces translations of f via
(LG→L2(R d) g f)(g ′) = f(g−1 · g′) = f(x′ − x).
Now, since the space of translations can be identified with Rd, the lifting (Eq. (6)) and group correlations (Eq. (7)) coincide with the standard definition of cross-correlation. I.e.,
(k?̃f)(x) = (k?f)(x) = ∫ R2 k(x′ − x)f(x′)dx′.
The logarithmic map is simply given by

heading: Log g = x,
text: so the convolution kernels can be defined by regular splines on Rd, cf. Eq. (11).
B.2 THE 2D ROTATION GROUP H = SO(2)
The special orthogonal group SO(2) consists of all orthogonal 2×2 matrix with determinant 1, i.e., rotation matrices of the form
Rθ = ( cos θ − sin θ sin θ cos θ ) ,
and the group product and inverse is given by the matrix product and matrix inverse:
h · h′ = (Rθ.Rθ′) = (Rθ+θ′) h−1 = ( R−1θ ) ,
with h = (Rθ), h′ = (Rθ′) ∈ SO(2) with θ, θ′ ∈ S1. The identity element is e = (R0) = (I). The action of H on R2 is given by matrix vector multiplication:
h x = Rθ.x, with x ∈ R2. Combining the group structure of 2D translations with rotations in SO(2) as a semi-direct product group gives us the roto-translation group SE(2) = R2oSO(2), also known as the special Euclidean motion group. The group structure of SE(2) is automatically derived from that of SO(2), see Sec. 3.1 for details.
The left-regular representation of SO(2) are given by
(LSO(2)→L2(R 2) h f)(x ′) = f(h−1 x′) = f(R−1θ .x ′),
(LSO(2)→L2(SO(2))h F )(h ′) = F (h−1 · h′) = F (Rθ′−θ),
(LSO(2)→L2(SE(2))h F )(x ′, h′) = F (h−1 x′, h−1 · h′) = F (R−1θ .x ′,Rθ′−θ).
Note that the latter two representations in terms of the rotation parameters θ, θ′ ∈ S1 represents the periodic shift θ′ − θmod 2π along the rotation axis. See Fig. 6 for an illustration for the transformation of SE(2) convolution kernels via the representation LSO(2)→L2(SE(2))h . The determinant of the Jacobian of the action of H on Rd, see corollary 1, is |deth| = 1. Using the above group structure we can write out the explicit forms for the lifting (Eq. (6)) and group correlations (Eq. (7)) as follows:
SE(2)-lifting: (k?̃f)(x, θ) = ∫ R2 k(R−1θ .(x ′ − x))f(x′)dx′,
SE(2)-correlation: (K?F )(x, θ) = ∫ R2 ∫ S1 K(R−1θ .(x ′ − x), θ′ − θmod 2π)F (x′, θ′)dx′dθ′.
The logarithmic map on H is given by the matrix logarithm
LogRθ =
( 0 −θmod 2π
θmod 2π 0
) ∈ Te(SO(2)) = span {A1}
Which in terms of the Lie algebra basis {A1}withA1 = (
0 −1 1 0
) gives a vector with coefficient
a1 = θmod 2π. The B-spline basis, centered around each hi = (Rθi) ∈ H with scale sh ∈ R, as depicted in Fig. 3, is thus computed via
BR,n (
Log h−1i · h sh
) = BR,n ( θ − θi mod 2π
sh
) .

heading: B.3 SCALING GROUP H = (R+,×)
text: We call the positive real line R+, together with multiplication, the scaling group. The group product and inverse are given by
h · h′ = (ss′) h−1 = ( 1 s ) ,
with h = (s), h′ = (s′) ∈ H with s, s′ ∈ R+. The identity element is e = (1). The action of H on Rd is given by scalar multiplication
h x = sx, with x ∈ Rd. By combining 2D translations with scaling as a semi-direct product group we obtain the scale-translation group, which we denote with R2 o R+. The left-regular representation of the scaling group are given by
(LR +→L2(Rd) h f)(x ′) = f(h−1 x′) = f( 1sx ′),
(LR +→L2(R+) h F )(h ′) = F (h−1 · h′) = F ( s ′ s ),
(LR +→L2(R2×R+) h F )(x ′, h′) = F (h−1 x′, h−1 · h′) = F ( 1sx ′, s ′ s ).
The scaling of a R2 o R+ group convolution kernel is thus achieved by a planar scaling and a logarithmic shift along the scale axis. This is illustrated in Fig. 7. The determinant of the Jacobian of this action is |deth| = sd. Using the above group structure we can write out the explicit forms for the lifting (Eq. (6)) and group correlations (Eq. (7)) as follows:
R2 oR+-lifting: (k?̃f)(x, s) = ∫ R2 1 sd k( 1s′ (x ′ − x))f(x′)dx′,
R2 oR+-correlation: (K?F )(x, s) = ∫ R2 ∫ R+ 1 sd K( 1s′ (x ′ − x), ss′ )F (x ′, θ′)dx′ds′.
The logarithmic map on H is provided by the natural logarithm as follows

heading: Log h = ln s.
text: The B-spline basis, centered around each hi = (si) ∈ H with scale si ∈ R, as depicted in Fig. 3, is thus computed via
BR,n (
Log h−1i · h sh
) = BR,n ( ln si −1s
sh
) .
B.4 THE 3D ROTATION GROUP H = SO(3)
The 3D rotation group is given by space of 3 × 3 orthogonal matrices with determinant 1, with the group product and inverse given by matrix product and matrix inverse:
h · h′ = (R.R′) h−1 = ( R−1 ) .
The action of SO(3) on R3 is given by matrix-vector multiplication h x = R.x,
with x ∈ R3. Combining the group structure of 3D translations with rotations in SO(3) as a semidirect product group gives us the roto-translation group SE(3) = R3 o SO(3), also known as the 3D special Euclidean motion group. The left-regular representation are then
(LSO(3)→L2(R 3) h f)(x ′) = f(h−1 x′) = f(R−1.x′), (LSO(3)→L2(SO(3))h F )(h ′) = F (h−1 · h′) = F (R−1.R′),
(LSO(3)→L2(SE(3))h F )(x ′, h′) = F (h−1 x′, h−1 · h′) = F (R−1.x′,R−1.R′).
The determinant of the Jacobian of the action of SO(3) on R3 is |deth| = 1. Using the above group structure we can write out the explicit forms for the lifting (Eq. (6)) and group correlations (Eq. (7)) as follows:
SE(3)-lifting: (k?̃f)(x,R) = ∫ R2 k(R−1.(x′ − x))f(x′)dx′,
SE(3)-correlation: (K?F )(x,R) = ∫ R2 ∫ SO(3) K(R−1.(x′ − x),R−1.R′)F (x′,R′)dx′dR′,
with dR′ denoting the Haar measure on SO(3), which depends on the parameterization, see details below.
The logarithmic map from the group SO(3) to the Lie algebra so(3) is given by the matrix logarithm and the resulting matrix can be expanded in a basis {A1, A2, A3} for the Lie algebra
LogR = 3∑ i=1 aiAi,
with
A1 = ( 0 0 0 0 0 −1 0 1 0 ) , A2 = ( 0 0 1 0 0 0 −1 0 0 ) , A3 = ( 0 −1 0 1 0 0 0 0 0 ) ,
which at the origin represent an infinitisimal rotation around the x, y, and z axis respectively. A cardinal B-spline centered at some Ri ∈ SO(3) with scale sh can then be computed in terms of these coefficients via BR,n ( Log R−1i .R
sh
) .
In practice it is often convenient to rely on a parameterization of the group and define the group structure in terms of these parameters. A common choice is to do this via ZYZ Euler angles via
Rα,β,γ = Rez,γ .Rey,β .Rez,α,
with Re,θ a rotation of θ around a reference axis e, and α ∈ [0, 2π], β ∈ [0, π], γ ∈ [0, 2π]. A Haar measure in terms of this parameterization is then given by dµ(R) = sinβdαdβdγ. We will use this parameterization in construction of the quotient group SO(3)/SO(2) next.
B.5 THE 2-SPHERE H = S2 ≡ SO(3)/SO(2)
Here we define the 2-sphere as a group quotient of SO(3) and remark that the same equations as in the SO(3) case (App. B.4) are used for the lifting and group correlations. This section describes how to construct convolution kernels on the sphere using the logarithmic map of SO(3), which can then be used to build G-CNNs on Rd o S2.
The 2-sphere is defined as S2 = {x ∈ R3 | ‖x‖ = 1}. Any point on the sphere can be obtained by rotating a reference vector z = (0, 0, 1)T with elements of SO(3), i.e., ∀n∈S2 ,∃R∈SO(3) : n = R.z. In other words, the group SO(3) acts transitively on S2. In ZYZ Euler angle parameterization of SO(3) all angles α leave the reference vector z in place, meaning that for each n ∈ S2 we have several R ∈ SO(3) that map z to the same n. As such, we can treat S2 as the quotient group SO(3)/SO(2), where SO(2) refers to the sub-group of rotations around the z-axis.
In order to define B-splines on the 2-sphere we need a logarithmic map from a point in S2 to the (Euclidean) tangent vector space Te(S2) at the origin. We will construct this logarithmic map using the Log define for SO(3). Let us parameterize the sphere with
n(β, γ) = Rez,γ .Rey,β .z.
Any rotation Rα,β,γ with arbitrary α maps to the same n(β, γ) ∈ S2. As such, there are also many vectorsA = LogRα,β,γ ∈ Te(SO(3)) that map to a suitable rotation matrix via the exponential map R = expA. We aim to find the vector in Te(SO(3)) for which c3 = 0, which via the exponential map generate torsion free exponential curves. The Log of any Rα,β,γ with α = −γ results in such a vector (Portegies et al., 2015). As such we define
LogS2 n(β, γ) := LogSO(3) R−γ,β,γ ,
which maps any point in S2 to a 2-dimensional vector space Te(S2) ⊂ Te(SO(3)). A B-spline on S2 can then be defined via
f(β, γ) = N∑ i=1 ciBR,n
( LogS2 R −1 0,βi,γi .R0,β,γ .z
sh
) , (16)
in which individual splines basis functions are centered around points n(βi, γi).
We remark that the group product R−1αi,βi,γi .R0,β,γ generates different rotations when varying αi, that however still map to the same n. The vectors obtained by taking LogS2 of the rotation matrices rotate with the choice for αi. Since the B-splines are approximately isotropic we neglect this effect and simply set αi = 0 in Eq. (16). Finally, we remark that the superposition of shifted splines (as in Eq. (16)) is not isotropic by construction, which is desirable when using the spline as a convolution kernel to lift functions to SO(3). When constraining G-CNNs to generate feature maps on S2, the kernels are constrained to be isotropic. Alternatively on could stay on S2 entirely and resort to gauge-equivariant networks (Cohen et al., 2019), for which the proposed splines are highly suited to move from the discrete setting (as in (Cohen et al., 2019)) to the continuous setting, see also App. C.2. For examples of splines on S2 see Figs. 3 and 4.

heading: C RELATED WORK

heading: C.1 DEEP SCALE-SPACS

heading: C.1.1 SCALE SPACE LIFTING AND CORRELATIONS
text: In (Worrall & Welling, 2019) images f ∈ L2(Rd) are lifted to a space of positions and scale parameters by constructing a scale space via
f↑(x, s) := fs(x) := (Gs ? f)(x)
with Gs(x) = (4πs)−d/2e− ‖x‖2
4s . The kernels and images are sampled on a discrete grid. Let Ω ⊂ Zd be the support of the kernel. Then the discrete scale space correlation is given by (Worrall & Welling, 2019, Eq. (19))
(K ?S f ↑)(x, s) = ∑ x̃∈Ω ∑ s̃∈Hd K(x̃, s̃)fs̃s(sx̃+ x),
with Hd the discretized set of scales, e.g., Hd = { 2i−1 }Nh i=1
, where we remark that here we use the convention of scaling of a function by s ≥ 1 instead of using the dilation parameter a = 1s in (Worrall & Welling, 2019). Next we remark that the scale space correlation without any scale interaction (Hd = {1}) is defined by a 2D correlation kernel via
(k ?S f ↑)(x, s) = ∑ x̃∈Zd k(x̃)fs(sx̃ + x),
which can be regarded as a discrete atrous/dilated correlation on each of the scale slices of f↑(x, s) with kernels dilated by a factor s.

heading: C.1.2 RELATION TO LIFTING CORRELATIONS (EQ. (6)) WITH B-SPLINES
text: Let our lifting correlation kernel k be given in a B-spline basis via Eq. (11) and let c : Ω ⊂ Zd → R be the map that assigns the weights to each B-spline center xi ∈ Ω with Ω the set of spline centers (i.e. c(xi) = ci). Let the Gaussian kernel Gs(x) be approximated by a scaled B-spline (up to a factor) and define BR
d,n s := 1 sd BR d,n( 1sx). With such an approximation (see also (Bouma et al., 2007)) our lifting correlation via Eq. (6) coincides with the lifting of (Worrall & Welling, 2019) followed by their non-scale interacting scale-space correlation, i.e.,
(k?̃f)(x, h) = (c ?S f ↑)(x, s),
with c : Ω→ R the spline coefficients. We show this by rewriting
(k?̃f)(x, h) = ∫ Rd 1 | deth|k(h −1 (x̃− x))f(x̃)dx̃
= ∫ Rd 1 sd ∑ xi∈Ω c(xi)B Rd,n( 1s (x̃− x)− xi)f(x̃)dx̃
= ∑ xi∈Ω c(xi) ∫ Rd BR d,n s (x̃− (x + sxi))f(x̃)dx̃
= ∑ xi∈Ω c(xi)(B Rd,n s ? f)(x + sxi)
≈ ∑ xi∈Ω c(xi)fs(x + sxi).

heading: C.2 GAUGE EQUIVARIANT NETWORKS

heading: C.2.1 GAUGE EQUIVARIANT CORRELATION
text: The following highlights commonalities between this paper and the work by Cohen et al. (2019) with respect to use of left-invariant vector fields in equivariant neural networks. Consider some Lie group G with Lie algebra g = Te(G), the exponential map Exp : g → G, and logarithmic map Log : G → g. Consider the group correlation between a kernel and function K,F : G → R, given in Eq. (7), which for B-spline kernels K with finite support Ω := supp(K) ⊂ G reads as
(K ? F )(g) = ∫ gΩ K(g−1 · g̃)F (g̃)µG(g̃), (17)
with µG(h) the Haar measure on G, and where write K for the convolution kernel on G, and K̃ for the corresponding kernel on g:
K(g) = ⇀ K(Log(g)), (18)
Let ⇀ Ω := supp( ⇀ K) ⊂ g be the support of ⇀ K. Finally let ⇀
Ω be localized such that Exp is a diffeomorphism (i.e., Exp( ⇀ Ω) = Ω and Log(Ω) = ⇀ Ω).
Now consider the definition of gauge equivariant correlation on manifolds as given by Eq. (3) of Cohen et al. (2019) for the case of scalar functions (in which case the trivial representation ρ(g) = 1 is to be used). In this case integration takes place over the Lie algebra, and gauge equivariant correlation is defined by
( ⇀ K ⇀ ?F )(g) := ∫ Ω ⇀ K(x)F (Expg(x))dx, (19)
with dx the Lebesgue measure on Rd, and where Expg denotes the exponential map from Tg(G)→ G. In our Lie group setting all tangent spaces can be identified with the tangent space at the origin (via the push-forward of left-multiplication) and we are able to write Expg := g · Expx. In the setting of gauge equivariant CNNs as in (Cohen et al., 2019) the exponential maps are generally dependent on g, using a separate reference/gauge frame (basis for the tangent space) at each g.

heading: C.2.2 RELATION TO G-CORRELATIONS (EQ. (7)) WITH B-SPLINES
text: For Lie groups the following identity holds between group correlations with localized B-splines on the one hand, in which integration takes place over the group G and elements are mapped to the algebra via Log, and gauge equivariant correlation on the other hand, in which integration takes place on the tangent spaces and vectors in these tangent spaces are mapped to the manifolds via Exp. In other words, given the definition of G cross-correlation in (7), denoted with ?, and gauge correlation in (19) or (Cohen et al., 2019, Eq. (1)), denoted with ⇀?, the two operators relate via
(K ? F )(g) = ( ⇀ K ⇀ ?F )(g). (20)
We show this by deriving
(K ? F )(g) = ∫ gΩ K(g−1 · g̃)F (g̃)dg̃
1 = ∫ Ω K(g̃)F (g · g̃)dg̃
2 = ∫ ⇀ Ω ⇀ K(x̃)F (g · Exp(x̃))dx̃,
= ∫ ⇀ Ω ⇀ K(x̃)F (Expg(x̃))dx̃ = ( ⇀ K?̃F )(g).
In the above dg̃ is a Haar measure on G. At 1= the substitution g̃ → g · g̃ is made and left-invariance of the Haar measure is used (d(g · g̃) = dg̃). At 2= we switch from integration over the region Ω in the Lie group to integration over region ⇀ Ω = Log(Ω) in the Lie algebra.

heading: D G-CNN ARCHITECTURES
text: This section describes the G-CNN architectures used in the experiments of Sec. 4 using the layers as defined in Sec. 3.2 and illustrated in Fig. 1. Two slightly different architectures are in the two different tasks (metastasis classifiation and landmark detection), but both are regular sequential GCNNs that start with a lifting layer (6), followed by a several group correlation layers (7), possibly alternated with spatial max-pooling, followed by a projection over H via (8), and end with a 1 × 1 convolution or fully connected layers. The architectures are summarized in Table. 1 and 2. Note that the output of the PCam architecture is two probabilities (1 for each class), and the output of the CelebA is five heatmaps (1 for each landmark).

heading: D.1 PCAM
text: The architecture for metastasis classification in the PCam dataset is given in Tab. 1. The input (64 × 64) is first cropped to 88 × 88 and is then used as input for the first layer (the lifting layer). None of the layers use spatial padding such that the image is eventually cropped to size 1× 1. Each layer is followed by batch normalization1 and a ReLU activation function, except for the last layer (layer 7) which is followed by adding a bias vector of length 2 and a softmax.
Note that the first five layers, including max pooling over rotations, encode the image into a 64- dimensional rotation invariant feature vector. The final two layers (6 and 7) can be regarded as a classical neural network classifier.
To reduce a possible orientation bias we aim to approximate the support of the kernels with a disk, rather than a rectangle. We do this by only considering splines with basis function centers {xi ∈ Zd | ‖xi‖ ≤ r, with radius r. For 5 × 5 kernels we set r = √ 5 by which we discard the basis functions at the corners of the 5 × 5 grid. The grid on H is uniformly sampled with Nh samples, giving the discretized gridHd = {(i−1)∗ 2πNh } Nh i=1. The global kernels (both dense and atrous) have
1We apply batch normalization over the domain of the feature maps, so over X = Rd or over X = Rd×H , as in (Cohen & Welling, 2016).
their centers also equidistant and globally cover S1, i.e., hi ∈ {(i− 1) ∗ 2πNk } Nk i=1, with the scales of the dense and atrous kernels respectively given by sh = 2πNk and sh = 2π Nh . The localized kernels have their centers on the grid with hi ∈ {i 2πNh } bNk/2c i=−bNk/2c and have scale sh = 2π Nh .
Finally, we follow the same data-augmentations at train time as proposed in (Liu et al., 2015). These include geometric augmentations such as 90◦ rotations and horizontal flips, as well as color augmentations such as brightness, saturation, hue and contrast variations.

heading: D.2 CELEBA
text: The architecture for landmark detecion in the CelebA dataset is biven in Tab. 2. The input is formatted according to the details in Sec. 4. In each layer zero padding is used in order to map the 128×128 input images to a 128 × 128 output heatmaps. Each layer is followed by batch normalization and a ReLU activation function, except for the last layer (layer 10) which is followed by adding a bias vector and a logistic sigmoid activation function.
Note that the result of the first 6 layers, including average pooling over scale, assign locally scaleinvariant feature vectors to each pixel. The final layers convert these feature maps into heatmaps via regular 2D convolutions.
Landmarks are localized via the argmax on each heatmap. The results in Fig. 5 show the success rate for localizing a landmark correctly. The success rate is computed as the average fraction of successful detections for all five landmarks in all images. A detection is considered successful if the distance to the actual landmark is less then 10 pixels.
The H axis is uniformly sampled (w.r.t. to the metric on H) on a fixed scale range, generating the discrete sets Hd = {e(i−1)sh}Nhi=1 with sh = 12 ln 2. The global kernels have their centers on this grid, i.e., hi ∈ Hd with the scale parameter the same as that of the grid. The local kernels also have their centers equidistant (with scale sh) to eachother, but are localized and given by hi ∈ {eish}bNk/2ci=−bNk/2c.

Review summaries:
In this paper, a framework for building group CNN with an arbitrary Lie group G is proposed. Generally, such a group CNN consists of 3 types of layers: a lifting layer which lifts a 2D image to a 3D data (G-image) whose domain is G; a group correlation layer which computes a 3D G-image from a 3D G-image; and a projection layer from a 3D G-image to a 2D image. To implement the convolutions in the lifting layer and group correlation layer which are defined in the continuous setting, the B-Spline basis functions are applied to expand the convolution kernels. Experimental results on tumor clarification and landmark localization show the superiority over CNN.