1 Introduction
Understanding 3D shape semantics from a large collection of 3D geometries has been a popular research direction over the past few years in both the graphics and vision communities. Many applications such as autonomous driving, robotics, and biostructure analysis depend on the ability to analyze 3D shape collections and the information associated with them.
Background
It is common practice to encode 3D shape information such as segmentation masks, geometric features, keypoints, reflectance, materials, etc. as perpoint functions defined on the shape surface, known as probe functions. We are interested, in a joint analysis setting, in discovering common latent structures among such probe functions defined on a collection of related 3D shapes. With the emergence of large 3D shape databases Chang et al. (2015), a variety of datadriven approaches, such as cycleconsistencybased optimization Huang et al. (2014)
and spectral convolutional neural networks
Bruna et al. (2014), have been applied to a range of tasks including semisupervised part cosegmentation Huang et al. (2011, 2014)and supervised keypoint/region correspondence estimation
Yi et al. (2017).However, one major obstacle in joint analysis is that each 3D shape has its own individual functional space, and linking related functions across shapes is challenging. To clarify this point, we contrast 3D shape analysis with 2D image processing. Under the functional point of view, each 2D image is a function defined on the regular 2D lattice, so all images are functions over a common underlying parameterizing domain. In contrast, with discretized 3D shapes, the probe functions are generally defined on heterogeneous shape graphs/meshes, whose nodes are points on each individual shape and edges link adjacent points. Therefore, the functional spaces on different 3D shapes are independent and not naturally aligned, making joint analysis over the probe functions nontrivial.
To cope with this problem, in the classical framework, ideas from manifold harmonics and linear algebra have been introduced. To analyze meaningful functions that are often smooth, a compact set of basis functions are computed by the eigendecomposition of the shape graph/mesh Laplacian matrix. Then, to relate basis functions across shapes, additional tools such as functional maps must be introduced Ovsjanikov et al. (2012) to handle the conversions among functional bases. This, however, raises further challenges since functional map estimation can be challenging for nonisometric shapes, and errors are often introduced in this step. In fact, functional maps are computed from corresponding sets of probe functions on the two shapes, something which we neither assume nor need.
Approach
Instead of a twostage procedure to first build independent functional spaces and then relate them through correspondences (functional or traditional), we propose a novel correspondencefree framework that directly learns consistent bases across a shape collection that reflect the shared structure of the set of probe functions. We produce a compact encoding for meaningful functions over a collection of related 3D shapes by learning a small functional basis for each shape using neural networks. The set of functional bases of each shape, a.k.a a shapedependent dictionary, is computed as a set of functions on a point cloud representing the underlying geometry — a functional set whose span will include probe functions on that shape. The training is accomplished in a very simple manner by giving the network sequences of pairs consisting of a shape geometry (as point clouds) and a semantic probe function on that geometry (that should be in the associated basis span). Our shapes are correlated, and thus the semantic functions we train on reflect the consistent structure of the shapes. The neural network will maximize its representational capacity by learning consistent bases that reflect this shared functional structure, leading in turn to consistent sparse function encodings. Thus, in our setting, consistent functional bases emerge from the network without explicit supervision.
We also demonstrate how to impose different constraints to the network optimization problem so that atoms in the dictionary exhibit desired properties adaptive to application scenarios. For instance, we can encourage the atoms to indicate smallest parts in the segmentation, or single points in keypoint detection. This implies that our model can serve as a collaborative filter that takes any mixture of semantic functions as inputs, and find the finest granularity that is the shared latent structure. Such a possibility can particularly be useful when the annotations in the training data are incomplete and corrupted. For examples, users may desire to decompose shapes into specific parts, but all shapes in the training data have only partial decomposition data without labels on parts. Our model can aggregate the partial information across the shapes and learn the full decomposition.
We remark that our network can be viewed as a function autoencoder, where the decoding is required to be in a particular format (a basis selection in which our function is compactly expressible). The resulting
canonicalization of the basis (the consistency we have described above) is something also recently seen in other autoencoders, for example in the quotientspace autoencoder of E. Mehr (2018) that generates shape geometry into a canonical pose.In experiments, we test our model with existing neural network architectures, and demonstrate the performance on labeled/unlabeled segmentation and keypoint correspondence problem in various datasets. In addition, we show how our framework can be utilized in learning synchronized basis functions with random continuous functions.
Contribution
Though simple, our model has advantages over the previous bases synchronization works Wang et al. (2013, 2014); Yi et al. (2017)
in several aspects. First, our model does not require precomputed basis functions. Typical bases such as Laplacian (on graphs) or LaplaceBeltrami (on mesh surfaces) eigenfunctions need extra preprocessing time for computation. Also, small perturbation or corruption in the shapes can lead to big differences. We can avoid the overhead of such preprocesssing by
predicting dictionaries while also synchronizing them simultaneously. Second, our dictionaries are applicationdriven, so each atom of the dictionary itself can attain a semantic meaning associated with smallscale geometry, such as a small part or a keypoint, while LB eigenfunctions are only suitable for approximating continuous and smooth functions (due to basis truncation). Third, the previous works define canonical bases, and the synchronization is achieved from the mapping between each individual set of bases and the canonical bases. In our model, the neural network becomes the synchronizer, without any explicit canonical bases. Lastly, compared with classical dictionary learning works that assume a universal dictionary for all data instances, we obtain a datadependent dictionary that allows nonlinear distortion of atoms but still preserves consistency. This endows us additional modeling power without sacrificing model interpretability.1.1 Related Work
Since much has already been discussed above, we only cover important missing ones here.
Learning compact representations of signals has been widely studied in many forms such as factor analysis and sparse dictionaries. Sparse dictionary methods learn an overcomplete basis of a collection of data that is as succinct as possible and have been studied in natural language processing
Deerwester et al. (1990); Hofmann (1999), timefrequency analysis Chen et al. (2001); Lewicki and Sejnowski (2000), video Olshausen (2002); Alfaro et al. (2016), and images Lee et al. (2007); Zeiler et al. (2010); Bristow et al. (2013). Encoding sparse and succinct representations of signals has also been observed in biological neurons
Olshausen and Field (1997, 1996, 2004).Since the introduction of functional maps Ovsjanikov et al. (2012), shape analysis on functional spaces has also been further developed in a variety of settings Pokrass et al. (2013); Kovnatsky et al. (2015); Huang et al. (2014); Eynard et al. (2016); Rodolà et al. (2016); Nogneng and Ovsjanikov (2017)
, and mappings between precomputed functional spaces have been studied in a deep learning context as well
Litany et al. (2017). In addition to our work, deep learning on point clouds has also been done on shape classification Qi et al. (2017a, b); Klokov and Lempitsky (2017); Wang et al. (2018b), semantic scene segmentation Huang et al. (2018), instance segmentation Wang et al. (2018a), and 3D amodal object detection Qi et al. (2018). We bridge these areas of research in a novel framework that learns, in a datadriven endtoend manner, dataadaptive dictionaries on the functional space of 3D shapes.2 Problem Statement
Given a collection of shapes , each of which has a sample function of specific semantic meaning (e.g. indicator of a subset of semantic parts or keypoints), we consider the problem of sharing the semantic information across the shapes, and predicting a functional dictionary for each shape that linearly spans all plausible semantic functions on the shape ( denotes the neural network weights). We assume that a shape is given as points sampled on its surface, a function
is represented with a vector in
(a scalar per point), and the atoms of the dictionary are represented as columns of a matrix , where is a sufficiently large number for the size of the dictionary. Note that the column space of can include any function if it has all Dirac delta functions of all points as columns. We aim at finding a much lowerdimensional vector space that also contains all plausible semantic functions. We also force the columns of to encode atomic semantics in applications, such as atomic instances in segmentation, by adding appropriate constraints.3 Deep Functional Dictionary Learning Framework
General Framework
We propose a simple yet effective loss function, which can be applied to any neural network architecture processing a 3D geometry as inputs. The neural network takes pairs of a shape
including points and a function as inputs in training, and outputs a matrix as a dictionary of functions on the shape. The loss function needs to be designed for minimizing both 1) the projection error from the input function to the vector space , and 2) the number of atoms in the dictionary matrix. This gives us the following loss function:(1) 
where is a linear combination weight vector, is a weight for regularization. is a function that measures the projection error, and the norm is a regularizer inducing structured sparsity, encouraging more columns to be zero vectors. We may have a set of constraints on both and depending on the applications. For example, when the input function is an indicator (binary) function, we constrain all elements in both and to be in range. Other constraints for specific applications are also introduced at the end of this section.
Note that our loss minimization is a minmin optimization problem; the inner minimization, which is embedded in our loss function in Equation 1, optimizes the reconstruction coefficients based on the shape dependent dictionary predicted by the network, and the outer minimization, which minimizes our loss function, updates the neural network weights to predict a best shape dependent dictionary. The nested minimization generally does not have an analytic solution due to the constraint on . Thus, it is not possible to directly compute the gradient of without . We solve this by an alternating minimization scheme as described in Algorithm 1. In a single gradient descent step, we first minimize over with the current , and then compute the gradient of while fixing . The minimization over is a convex quadratic programming, and the scale is very small since is a very thin matrix (). Hence, a simplex method can very quickly solve the problem in every gradient iteration.
Adaptation in Weaklysupervised Cosegmentation
Some constraints for both and can be induced from the assumptions of the input function and the properties of the dictionary atoms. In the segmentation problem, we take an indicator function of a set of segments as an input, and we desire that each atom in the output dictionary indicates an atomic part (Figure 1 (a)). Thus, we restrict both and to have values in the range. Also, the atomic parts in the dictionary must partition the shape, meaning that each point must be assigned to one and only one atom. Thus, we add sumtoone constraint for every row of . The set of constraints for the segmentation problem is defined as follows:
(2) 
where is the th element of matrix , and and are vectors/matrices with an appropriate size. The first constraint on is incorporated in solving the inner minimization problem, and the second and third constraints on can simply be implemented by using softmax activation at the last layer of the network.
Adaptation in Weaklysupervised Keypoint Correspondence Estimation
Similarly with the segmentation problem, the input function in the keypoint correspondence problem is also an indicator function of a set of points (Figure 1 (b)). Thus, we use the same range constraint for both and . Also, each atom needs to represent a single point, and thus we add sumtoone constraint for every column of :
(3) 
For robustness, a distance function from the keypoints can be used as input instead of the binary indicator function. Particularly some neural network architectures such as PointNet Qi et al. (2017a) do not exploit local geometry context. Hence, a spatially localized distance function can avoid overfitting to the Dirac delta function. We use a normalized Gaussianweighed distance function in our experiment: , where is th element of the distance function from the keypoint , is th point coordinates, is Euclidean distance, and is Gaussianweighting parameter (0.001 in our experiment). The distance function is normalized to sum to one, which is consistent with our constraints in Equation 3. The sum of any subset of the keypoint distance functions becomes an input function in our training.
Adaptation in Smooth Function Approximation and Mapping
For predicting atomic functions whose linear combination can approximate any smooth function, we generate the input function by taking a random linear combination of LB bases functions (Figure 1 (c)). We also use a unit vector constraint for each atom of the dictionary:
(4) 
4 Experiments
We demonstrate the performance of our model in keypoint correspondence and segmentation problems with different datasets. We also provide qualitative results of synchronizing atomic functions on nonrigid shapes. While any neural network architecture processing 3D geometry can be employed in our model (e.g. PointNet Qi et al. (2017a), PointNet++ Qi et al. (2017b), KDNET Klokov and Lempitsky (2017), DGCNN Wang et al. (2018b), ShapePFCN Kalogerakis et al. (2017)), we use PointNet Qi et al. (2017a) architecture in the experiments due to its simplicity. Note that our output is a set of dimensional row vectors for all points. Thus, we can use the PointNet segmentation architecture without any modification. Code for all experiments below is available in https://github.com/mhsung/deepfunctionaldictionaries.
4.1 ShapeNet Keypoint Correspondence
Yi et al. (2017) provide keypoint annotations on 6,243 chair models in ShapeNet Chang et al. (2015). The keypoints are manually annotated by experts, and all of them are matched and aligned across the shapes. Each shape has up to 10 keypoints, while most of the shapes have missing keypoints. In the training, we take a random subset of keypoints of each shape to feed an input function, and predict a function dictionary in which atoms indicate every single keypoint. In the experiment, we use a 8020 random split for training/test sets ^{1}^{1}1Yi et al. (2017) use a select subset of models in their experiment, but this subset is not provided by the authors. Thus, we use the entire dataset and make our own train/test split., train the network with point clouds as provided by Yi et al. (2017), and set and .
Figure 2 (at the top) illustrates examples of predicted keypoints when picking the points having the maximum value in each atom. The colors denote the order of atoms in dictionaries, which is consistent across all shapes despite their different geometries. The outputs are also evaluated by the percentage of correct keypoints (PCK) metric as done in Yi et al. (2017) while varying the Euclidean distance threshold (Figure 2 at the bottom). We report the results for both when finding the best onetoone correspondences between the ground truth and predicted keypoints for each shape (red line) and when finding the correspondences between ground truth labels and atom indices for all shapes (green line). These two plots are identical, meaning that the order of predicted keypoints is rarely changed in different shapes. Our results also outperform the previous works Huang et al. (2013); Yi et al. (2017) by a big margin.
4.2 ShapeNet Semantic Part Segmentation
ShapeNet Chang et al. (2015) contains 16,881 shapes in 16 categories, and each shape has semantic part annotations Yi et al. (2016) for up to six segments. Qi et al. (2017a) train PointNet segmentation using shapes in all categories, and the loss function is defined as the cross entropy per point with all labels. We follow their experiment setup by using the same split of training/validation/test sets and the same sampled point cloud as inputs. The difference is that we do not leverage the labels of segments in training, and consider the parts as unlabeled segments. We also deal with the more general situation that each shape may have incomplete segmentation by taking an indicator function of a random subset of segments as an input.
mean  airplane  bag  cap  car  chair  earphone  guitar  knife  lamp  laptop  motorbike  mug  pistol  rocket  skateboard  table  

PointNet Qi et al. (2017a)  82.4  81.4  81.1  59.0  75.6  87.6  69.7  90.3  83.9  74.6  94.2  65.5  93.2  79.3  53.2  74.5  81.3 
Ours  84.6  81.2  72.7  79.9  76.5  88.3  70.4  90.0  80.5  76.1  95.1  60.5  89.8  80.8  57.1  78.3  88.1 
mean  airplane  bag  cap  car  chair  earphone  guitar  knife  lamp  laptop  motorbike  mug  pistol  rocket  skateboard  table  

Ours (per shape)  84.6  81.2  72.7  79.9  76.5  88.3  70.4  90.0  80.5  76.1  95.1  60.5  89.8  80.8  57.1  78.3  88.1 
Ours (per cat.)  77.3  79.0  67.5  66.9  75.4  87.8  58.7  90.0  79.7  37.1  95.0  57.1  88.8  78.4  46.0  75.8  78.4 
Evaluation
For evaluation, we binarize
by finding the maximum value in each row, and consider each column as an indicator of a segment. The accuracy is measured based on the average of each shape mean IoU similarly with Qi et al. (2017a), but we make a difference since our method does not exploit labels. In ShapeNet, some categories have optional labels, and shapes may or may not have a part with these optional labels (e.g. armrests of chairs). Qi et al. (2017a) take into account the optional labels even when the segment does not exist in a shape ^{2}^{2}2IoU becomes zero if the label is assigned to any point in prediction, and one otherwise.. But we do not predict labels of points, and thus such cases are ignored in our evaluation.We first measure the performance of segmentation by finding the correspondences between ground truth and predicted segments for each shape. The best onetoone correspondences are found by running the Hungarian algorithm on mean IoU values. Table 1 shows the results of our method when using and , and the results of the labelbased PointNet segmentation Qi et al. (2017a). When only considering the segmentation accuracy, our approach outperforms the original PointNet segmentation trained with labels.
We also report the average mean IoUs when finding the best correspondences between part labels and the indices of dictionary atoms per category. As shown in Table 2, the accuracy is still comparable in most categories, indicating that the order of column vectors in are mostly consistent with the semantic labels. There are a few exceptions; for example, lamps are composed of shades, base, and tube, and half of lamps are ceiling lamps while the others are standing lamps. Since PointNet learns perpoint features from the global coordinates of the points, shades and bases are easily confused when their locations are switched (Figure 3). Such problem can be resolved when using a different neural network architecture learning more from the local geometric contexts. For more analytic experiments, refer to the supplementary material.
mean  ceiling  floor  wall  beam  column  window  door  table  chair  sofa  bookcase  board  

SGPN Wang et al. (2018a)  64.7  67.0  71.4  66.8  54.5  45.4  51.2  69.9  63.1  67.6  64.0  54.4  60.5 
Ours  69.1  95.4  99.2  77.3  48.0  39.2  68.2  49.2  56.0  53.2  35.3  31.6  42.2 
4.3 S3DIS Instance Segmentation
Stanford 3D Indoor Semantic Dataset (S3DIS) Armeni et al. (2016)
is a collection of real scan data of indoor scenes with annotations of instance segments and their semantic labels. When segmenting instances in such data, the main difference with the semantic segmentation of ShapeNet is that there can exist multiple instances of the same semantic label. Thus, the approach of classifying points with labels is not applicable. Recently,
Wang et al. (2018a)tried to solve this problem by leveraging the PointNet architecture. Their framework named SGPN learns a similarity metric among points, enabling every point to generate a instance proposal based on proximity in the learned feature space. The perpoint proposals are further merged in a heuristic post processing step. We compare the performance of our method with the same experiment setup with SGPN. The input is a
point cloud of a floor block in the scenes, and each block contains up to 150 instances. Thus, we use and . Refer to Wang et al. (2018a) for the details of the data preparation. In the experiments of both methods, all 6 areas of scenes except area 5 are used as a training set, and the area 5 is used as a test set.Evaluation
We evaluate the performance of instance proposal prediction in each block of the scenes. ^{3}^{3}3Wang et al. (2018a) propose a heuristic process of merging prediction results of each block and generating instance proposals in a scene, but we measure the performance for each block in order to factor out the effect of this postprocessing step.
As an evaluation metric, we use proposal recall
Hosang et al. (2016), which measures the percentage of ground truth instances covered by any prediction within a given IoU threshold. In both SGPN and our model, the outputs are nonoverlapped segments, thus we do not consider the number of proposals in the evaluation. Figure 6 depicts the proposal recall of both methods when varying the IoU threshold from 0.5 to 1.0. The recall of our method is greater than the baseline throughout all threshold levels. The recalls for each semantic part label with IoU threshold 0.5 are reported in Table 3. Our method performs well specifically for large objects such as ceilings, floors, walls, and windows. Note that Wang et al. (2018a) start their training from a pretrained model for semantic label prediction, and also their framework consumes point labels as supervision in the training to jointly predict labels and segments. Our model is trained from scratch and labelfree.Consistency with semantic labels
Although it is hard to expect strong correlations among semantic part labels and the indices of dictionary atoms in this experiment due to the large variation of scene data, we still observe weak consistency between them. Figure 6 illustrates confusion among semantic part labels. This confusion is calculated by first creating a vector for each label in which the th element indicates the count of the label in the th atom, normalizing this vector, and taking a dot product for every pair of labels. Ceilings and floors are clearly distinguished from the others due to their unique positions and scales. Some groups of objects having similar heights (e.g. doors, bookcases, and boards; chairs and sofas) are confused with each other frequently, but objects in different groups are discriminated well.
4.4 MPIFAUST Human Shape Bases Synchronization
In this experiment, we aim at finding synchronized atomic functions in a collection of shapes for which linear combination can approximate any continuous function. Such synchronized atomic functions can be utilized in transferring any information on one shape to the other without having pointwise correspondences. Here, we test with 100 nonrigid human body shapes in MPIFAUST dataset Bogo et al. (2014). Since the shapes are deformable, it is not appropriate to process Euclidean coordinates of a point cloud as inputs. Hence, instead of a point cloud and the PointNet, we use HKS Sun et al. (2009) and WKS Aubry et al. (2011) point descriptors for every vertex, and process them using 7 residual layers shared for all points as proposed in Litany et al. (2017). The point descriptors cannot clearly distinguish symmetric parts in a shape, so the output atomic functions also become symmetric. To break the ambiguity, we sample four points using farthest point sampling in each shape, find their onetoone correspondences in other shapes using the same point descriptor, and use the geodesic distances from these points as additional point features. As input functions, we compute LaplaceBeltrami operators on shapes, and take a random linear combination of the first ten bases.
Figure 8 visualizes the output atomic function when we train the network with and . The order of atomic functions are consistent in all shapes. In Figure 8, we show how the information in one shape is transferred to the other shapes using our atomic functions. We project the indicator function of each segment (at left in figure) to the function dictionary space of the base shape, and unproject them in the function dictionary space of the other shapes. The transferred segment functions are blurry since the network is trained with only continuous functions, but still indicate proper areas of the segments.
5 Conclusion
We have investigated a problem of jointly analyzing probe functions defined on different shapes, and finding a common latent space through a neural network. The learning framework we proposed predicts a function dictionary of each shape that spans input semantic functions, and finds the atomic functions in a consistent order without any correspondence information. Our framework is very general, enabling easy adaption to any neural network architecture and any application scenario. We have shown some examples of constraints in the loss function that can allow the atomic functions to have desired properties in specific applications: the smallest parts in segmentation, and single points in keypoint correspondence.
In the future, we will further explore the potential of our framework to be applied to various applications and even in different data domains. Also, we will investigate how the power of neural network decomposing a function space to atoms can be enhanced through different architectures and a hierarchical basis structure.
Acknowledgments
We thank the anonymous reviewers for their comments and suggestions. This project was supported by a DoD Vannevar Bush Faculty Fellowship, NSF grants CHS1528025 and IIS1763268, and an Amazon AWS AI Research gift.
References
 Alfaro et al. [2016] Anali Alfaro, Domingo Mery, and Alvaro Soto. Action recognition in video using sparse coding and relative features. In CVPR, 2016.
 Armeni et al. [2016] I. Armeni, O. Sener, A. R. Zamir, H. Jiang, I. Brilakis, M. Fischer, and S. Savarese. 3d semantic parsing of largescale indoor spaces. In CVPR, 2016.
 Aubry et al. [2011] M. Aubry, U. Schlickewei, and D. Cremers. The wave kernel signature: A quantum mechanical approach to shape analysis. In ICCV Workshops, 2011.
 Bogo et al. [2014] Federica Bogo, Javier Romero, Matthew Loper, and Michael J. Black. FAUST: Dataset and evaluation for 3D mesh registration. In CVPR, 2014.
 Bristow et al. [2013] Hilton Bristow, Anders Eriksson, and Simon Lucey. Fast convolutional sparse coding. In CVPR, 2013.
 Bruna et al. [2014] Joan Bruna, Wojciech Zaremba, Arthur Szlam, and Yann Lecun. Spectral networks and locally connected networks on graphs. In ICLR, 2014.
 Chang et al. [2015] Angel X. Chang, Thomas A. Funkhouser, Leonidas J. Guibas, Pat Hanrahan, QiXing Huang, Zimo Li, Silvio Savarese, Manolis Savva, Shuran Song, Hao Su, Jianxiong Xiao, Li Yi, and Fisher Yu. Shapenet: An informationrich 3d model repository. CoRR, abs/1512.03012, 2015.
 Chen et al. [2001] Scott Shaobing Chen, David L Donoho, and Michael A Saunders. Atomic decomposition by basis pursuit. SIAM review, 2001.
 Deerwester et al. [1990] Scott Deerwester, Susan T Dumais, George W Furnas, Thomas K Landauer, and Richard Harshman. Indexing by latent semantic analysis. Journal of the American society for information science, 1990.
 E. Mehr [2018] V. Guitteny M. Cord E. Mehr, N. Thome. Manifold learning in quotient spaces. In CVPR, 2018.
 Eynard et al. [2016] Davide Eynard, Emanuele Rodola, Klaus Glashoff, and Michael M Bronstein. Coupled functional maps. In 3DV, pages 399–407, 2016.

Hofmann [1999]
Thomas Hofmann.
Probabilistic latent semantic analysis.
In
Proceedings of the Fifteenth conference on Uncertainty in artificial intelligence
, 1999.  Hosang et al. [2016] Jan Hosang, Rodrigo Benenson, Piotr Dollar, and Bernt Schiele. What makes for effective detection proposals? IEEE TPAMI, 2016.
 Huang et al. [2013] QiXing Huang, Hao Su, and Leonidas Guibas. Finegrained semisupervised labeling of large shape collections. In SIGGRAPH Asia, 2013.
 Huang et al. [2018] Qiangui Huang, Weiyue Wang, and Ulrich Neumann. Recurrent slice networks for 3d segmentation on point clouds. In CVPR, 2018.

Huang et al. [2011]
Qixing Huang, Vladlen Koltun, and Leonidas Guibas.
Joint shape segmentation with linear programming.
In SIGGRAPH Asia, 2011.  Huang et al. [2014] Qixing Huang, Fan Wang, and Leonidas Guibas. Functional map networks for analyzing and exploring large shape collections. In SIGGRAPH, 2014.
 Kalogerakis et al. [2017] Evangelos Kalogerakis, Melinos Averkiou, Subhransu Maji, and Siddhartha Chaudhuri. 3D shape segmentation with projective convolutional networks. In CVPR, 2017.
 Klokov and Lempitsky [2017] Roman Klokov and Victor S. Lempitsky. Escape from cells: Deep kdnetworks for the recognition of 3d point cloud models. In ICCV, 2017.
 Kovnatsky et al. [2015] Artiom Kovnatsky, Michael M Bronstein, Xavier Bresson, and Pierre Vandergheynst. Functional correspondence by matrix completion. In CVPR, 2015.
 Lee et al. [2007] Honglak Lee, Alexis Battle, Rajat Raina, and Andrew Y Ng. Efficient sparse coding algorithms. In NIPS, 2007.
 Lewicki and Sejnowski [2000] Michael S Lewicki and Terrence J Sejnowski. Learning overcomplete representations. Neural computation, 2000.
 Litany et al. [2017] Or Litany, Tal Remez, Emanuele Rodola, Alex Bronstein, and Michael Bronstein. Deep functional maps: Structured prediction for dense shape correspondence. In CVPR, 2017.
 Nogneng and Ovsjanikov [2017] Dorian Nogneng and Maks Ovsjanikov. Informative descriptor preservation via commutativity for shape matching. In Eurographics, 2017.
 Olshausen [2002] Bruno A Olshausen. Sparse coding of timevarying natural images. Journal of Vision, 2002.
 Olshausen and Field [1996] Bruno A Olshausen and David J Field. Emergence of simplecell receptive field properties by learning a sparse code for natural images. Nature, 1996.
 Olshausen and Field [1997] Bruno A Olshausen and David J Field. Sparse coding with an overcomplete basis set: A strategy employed by v1? Vision research, 1997.
 Olshausen and Field [2004] Bruno A Olshausen and David J Field. Sparse coding of sensory inputs. Current opinion in neurobiology, 2004.
 Ovsjanikov et al. [2012] Maks Ovsjanikov, Mirela BenChen, Justin Solomon, Adrian Butscher, and Leonidas Guibas. Functional maps: a flexible representation of maps between shapes. In SIGGRAPH, 2012.
 Pokrass et al. [2013] Jonathan Pokrass, Alexander M Bronstein, Michael M Bronstein, Pablo Sprechmann, and Guillermo Sapiro. Sparse modeling of intrinsic correspondences. In Eurographics, 2013.
 Qi et al. [2018] Charles R Qi, Wei Liu, Chenxia Wu, Hao Su, and Leonidas J Guibas. Frustum pointnets for 3d object detection from rgbd data. In CVPR, 2018.
 Qi et al. [2017a] Charles Ruizhongtai Qi, Hao Su, Kaichun Mo, and Leonidas J. Guibas. Pointnet: Deep learning on point sets for 3d classification and segmentation. In CVPR, 2017a.
 Qi et al. [2017b] Charles Ruizhongtai Qi, Li Yi, Hao Su, and Leonidas J. Guibas. Pointnet++: Deep hierarchical feature learning on point sets in a metric space. In NIPS, 2017b.
 Rodolà et al. [2016] Emanuele Rodolà, Luca Cosmo, Michael M Bronstein, Andrea Torsello, and Daniel Cremers. Partial functional correspondence. In SGP, 2016.
 Sun et al. [2009] Jian Sun, Maks Ovsjanikov, and Leonidas Guibas. A concise and provably informative multiscale signature based on heat diffusion. In SGP, 2009.
 Wang et al. [2014] F. Wang, Q. Huang, M. Ovsjanikov, and L. J. Guibas. Unsupervised multiclass joint image segmentation. In CVPR, 2014.
 Wang et al. [2013] Fan Wang, Qixing Huang, and Leonidas J. Guibas. Image cosegmentation via consistent functional maps. In ICCV, 2013.
 Wang et al. [2018a] Weiyue Wang, Ronald Yu, Qiangui Huang, and Ulrich Neumann. Sgpn: Similarity group proposal network for 3d point cloud instance segmentation. In CVPR, 2018a.
 Wang et al. [2018b] Yue Wang, Yongbin Sun, Ziwei Liu, Sanjay E. Sarma, Michael M. Bronstein, and Justin M. Solomon. Dynamic graph cnn for learning on point clouds. arXiv, 2018b.
 Yi et al. [2016] Li Yi, Vladimir G. Kim, Duygu Ceylan, IChao Shen, Mengyan Yan, Hao Su, Cewu Lu, Qixing Huang, Alla Sheffer, and Leonidas Guibas. A scalable active framework for region annotation in 3d shape collections. In SIGGRAPH Asia, 2016.
 Yi et al. [2017] Li Yi, Hao Su, Xingwen Guo, and Leonidas J Guibas. Syncspeccnn: Synchronized spectral cnn for 3d shape segmentation. In CVPR, 2017.
 Zeiler et al. [2010] Matthew D Zeiler, Dilip Krishnan, Graham W Taylor, and Rob Fergus. Deconvolutional networks. In CVPR, 2010.
Supplementary Material
s.1 ShapeNet Semantic Part Segmentation – Analytic Experiments
Effect of and
In Table S3, we demonstrate the effect of changing parameter and . When the norm regularizer is not used (), the accuracy decreases as increases since parts can map to a number of smaller segments. After adding the regularizer with a weight , the accuracy becomes similar however we choose the number of columns . We found that the norm regularizer effectively forces the unnecessary columns to be close to a zero vector.
Training with partial segmentation
In the segmentation problem using unlabeled segments, learning from partial segmentation is a nontrivial task, while our method can easily learn segmentation from the partial information. To demonstrate this, we randomly select a set of parts in the entire training set with a fixed fraction, and ignore them when choosing a random subset of segments for input functions. The accuracy according to the fraction is shown in Table S3. Note that performance remains roughly the same even when we do not use 75% of segments in the training.
Training with noise
We test the robustness of our training system against noise in the input function . Table S3
describes the performance when switching each bit of the binary indicator function with a specific probability. The results show that our system is not affected by the small noise in the input functions.
0.0  0.5  1.0  

10  75.0  82.7  84.6 
25  71.2  83.8  85.2 
50  65.3  82.9  82.9 
Fraction  mIoU 

0.00  84.6 
0.25  86.1 
0.50  86.0 
0.75  84.5 
Probability  mIoU 

0.00  84.6 
0.05  85.8 
0.10  85.9 
0.20  85.1 
s.2 Siamese Structure for Correspondence Supervision
While our framework empirically performs well on generating consistent function dictionaries even without correspondences, we further investigate about how the correspondence supervision can be incorporated in our framework when it is provided. We consider the case when the correspondence information is given as a pair or functions in different shapes. Note that this setup does not require to have full correspondence information for all pairs. The correspondence of functions means that two functions are represented with the same linear combination weight when the order of dictionary atoms are consistent. Thus, we build a Siamese neural network structure processing two corresponding functions, and minimize the inner problem in the loss function Equation 1 jointly with the shared variable .
We test this approach with the ShapeNet part segmentation problem. Every time when feeding the input function in the training, we find the other shape that have a corresponding function, and randomly choose one of them. The comparison with the vanilla framework is shown in Table S4 and S5. and are used in both experiments. When finding the best onetoone correspondences between ground truth part labels and atom indices in each category, the Siamese structure shows 3.0% improvement in average mean IoU, meaning that the output dictionaries make less confusion when distinguishing semantic parts with the indices of atoms. It also gives better accuracy when finding the correspondences in each shape.
mean  airplane  bag  cap  car  chair  earphone  guitar  knife  lamp  laptop  motorbike  mug  pistol  rocket  skateboard  table  

Vanilla  77.3  79.0  67.5  66.9  75.4  87.8  58.7  90.0  79.7  37.1  95.0  57.1  88.8  78.4  46.0  75.8  78.4 
Siamese  80.3  78.6  73.7  44.8  76.9  87.7  65.0  90.6  85.2  60.4  94.7  60.5  93.6  78.5  55.8  76.1  80.1 
mean  airplane  bag  cap  car  chair  earphone  guitar  knife  lamp  laptop  motorbike  mug  pistol  rocket  skateboard  table  

Vanilla  84.6  81.2  72.7  79.9  76.5  88.3  70.4  90.0  80.5  76.1  95.1  60.5  89.8  80.8  57.1  78.3  88.1 
Siamese  85.6  82.2  75.7  74.5  77.5  88.4  73.5  91.0  85.2  77.9  95.9  63.4  93.6  80.7  62.4  80.7  88.9 
Comments
There are no comments yet.