3DShapeNets
for cvpr2015 project 3DShapeNets
view repo
3D shape is a crucial but heavily underutilized cue in today's computer vision systems, mostly due to the lack of a good generic shape representation. With the recent availability of inexpensive 2.5D depth sensors (e.g. Microsoft Kinect), it is becoming increasingly important to have a powerful 3D shape representation in the loop. Apart from category recognition, recovering full 3D shapes from viewbased 2.5D depth maps is also a critical part of visual understanding. To this end, we propose to represent a geometric 3D shape as a probability distribution of binary variables on a 3D voxel grid, using a Convolutional Deep Belief Network. Our model, 3D ShapeNets, learns the distribution of complex 3D shapes across different object categories and arbitrary poses from raw CAD data, and discovers hierarchical compositional part representations automatically. It naturally supports joint object recognition and shape completion from 2.5D depth maps, and it enables active object recognition through view planning. To train our 3D deep learning model, we construct ModelNet  a largescale 3D CAD model dataset. Extensive experiments show that our 3D deep representation enables significant performance improvement over thestateofthearts in a variety of tasks.
READ FULL TEXT VIEW PDF
With the advent of affordable depth sensors, 3D capture becomes more and...
read it
We propose a novel approach for 3D shape completion by synthesizing
mult...
read it
Mesh is an important and powerful type of data for 3D shapes and widely
...
read it
Recovering textured 3D models of nonrigid human body shapes is challeng...
read it
In viewbased 3D shape recognition, extracting discriminative visual
rep...
read it
3D convolutional neural networks (3DCNN) have been used for object
reco...
read it
The goal of this paper is to compare surfacebased and volumetric 3D obj...
read it
for cvpr2015 project 3DShapeNets
for cvpr2015 project 3DShapeNets
Since the establishment of computer vision as a field five decades ago, 3D geometric shape has been considered to be one of the most important cues in object recognition. Even though there are many theories about 3D representation (e.g. [5, 22]), the success of 3Dbased methods has largely been limited to instance recognition (e.g. modelbased keypoint matching to nearest neighbors [24, 31]). For object category recognition, 3D shape is not used in any stateoftheart recognition methods (e.g. [11, 19]), mostly due to the lack of a good generic representation for 3D geometric shapes. Furthermore, the recent availability of inexpensive 2.5D depth sensors, such as the Microsoft Kinect, Intel RealSense, Google Project Tango, and Apple PrimeSense, has led to a renewed interest in 2.5D object recognition from depth maps (e.g. Sliding Shapes [30]). Because the depth from these sensors is very reliable, 3D shape can play a more important role in a recognition pipeline. As a result, it is becoming increasingly important to have a strong 3D shape representation in modern computer vision systems.
Apart from category recognition, another natural and challenging task for recognition is shape completion: given a 2.5D depth map of an object from one view, what are the possible 3D structures behind it? For example, humans do not need to see the legs of a table to know that they are there and potentially what they might look like behind the visible surface. Similarly, even though we may see a coffee mug from its side, we know that it would have empty space in the middle, and a handle on the side.
Datadriven visualization: For each neuron, we average the top 100 training examples with highest responses ( 0.99) and crop the volume inside the receptive field. The averaged result is visualized by transparency in 3D (Gray) and by the average surface obtained from zerocrossing (Red). 3D ShapeNets are able to capture complex structures in 3D space, from lowlevel surfaces and corners at L1, to objects parts at L2 and L3, and whole objects at L4 and above. 
In this paper, we study generic shape representation for both object category recognition and shape completion. While there has been significant progress on shape synthesis [7, 17] and recovery [27], they are mostly limited to partbased assembly and heavily rely on expensive part annotations. Instead of handcoding shapes by parts, we desire a datadriven way to learn the complex shape distributions from raw 3D data across object categories and poses, and automatically discover a hierarchical compositional part representation. As shown in Figure 1, this would allow us to infer the full 3D volume from a depth map without the knowledge of object category and pose a priori. Beyond the ability to jointly hallucinate missing structures and predict categories, we also desire the ability to compute the potential information gain for recognition with regard to missing parts. This would allow an active recognition system to choose an optimal subsequent view for observation, when the category recognition from the first view is not sufficiently confident.
To this end, we propose 3D ShapeNets to represent a geometric 3D shape as a probabilistic distribution of binary variables on a 3D voxel grid. Our model uses a powerful Convolutional Deep Belief Network (Figure 2
) to learn the complex joint distribution of all 3D voxels in a datadriven manner. To train this 3D deep learning model, we construct ModelNet, a largescale object dataset of 3D computer graphics CAD models. We demonstrate the strength of our model at capturing complex object shapes by drawing samples from the model. We show that our model can recognize objects in singleview 2.5D depth images and hallucinate the missing parts of depth maps. Extensive experiments suggest that our model also generalizes well to real world data from the NYU depth dataset
[23], significantly outperforming existing approaches on singleview 2.5D object recognition. Further it is also effective for nextbestview prediction in view planning for active object recognition [25].There has been a large body of insightful research on analyzing 3D CAD model collections. Most of the works [7, 12, 17] use an assemblybased approach to build deformable partbased models. These methods are limited to a specific class of shapes with small variations, with surface correspondence being one of the key problems in such approaches. Since we are interested in shapes across a variety of objects with large variations and part annotation is tedious and expensive, assemblybased modeling can be rather cumbersome. For surface reconstruction of corrupted scanning input, most related works [3, 26]
are largely based on smooth interpolation or extrapolation. These approaches can only tackle small missing holes or deficiencies. Templatebased methods
[27] are able to deal with large space corruption but are mostly limited by the quality of available templates and often do not provide different semantic interpretations of reconstructions.The great generative power of deep learning models has allowed researchers to build deep generative models for 2D shapes: most notably the DBN [15] to generate handwritten digits and ShapeBM [10] to generate horses, etc. These models are able to effectively capture intraclass variations. We also desire this generative ability for shape reconstruction but we focus on more complex real world object shapes in 3D. For 2.5D deep learning, [29] and [13] build discriminative convolutional neural nets to model images and depth maps. Although their algorithms are applied to depth maps, they use depth as an extra 2D channel instead of modeling full 3D. Unlike [29], our model learns a shape distribution over a voxel grid. To the best of our knowledge, we are the first work to build 3D deep learning models. To deal with the dimensionality of high resolution voxels, inspired by [21]^{1}^{1}1The model is precisely a convolutional DBM where all the connections are undirected, while ours is a convolutional DBN., we apply the same convolution technique in our model.
Unlike static object recognition in a single image, the sensor in active object recognition [6] can move to new view points to gain more information about the object. Therefore, the NextBestView problem [25] of doing view planning based on current observation arises. Most previous works in active object recognition [9, 16] build their view planning strategy using 2D color information. However this multiview problem is intrinsically 3D in nature. Atanasov et al, [1, 2]
implement the idea in real world robots, but they assume that there is only one object associated with each class reducing their problem to instancelevel recognition with no intraclass variance. Similar to
[9], we use mutual information to decide the NBV. However, we consider this problem at the precise voxel level allowing us to infer how voxels in a 3D region would contribute to the reduction of recognition uncertainty.To study 3D shape representation, we propose to represent a geometric 3D shape as a probability distribution of binary variables on a 3D voxel grid. Each 3D mesh is represented as a binary tensor: 1 indicates the voxel is inside the mesh surface, and 0 indicates the voxel is outside the mesh (i.e., it is empty space). The grid size in our experiments is
.To represent the probability distribution of these binary variables for 3D shapes, we design a Convolutional Deep Belief Network (CDBN). Deep Belief Networks (DBN) [15] are a powerful class of probabilistic models often used to model the joint probabilistic distribution over pixels and labels in 2D images. Here, we adapt the model from 2D pixel data to 3D voxel data, which imposes some unique challenges. A 3D voxel volume with reasonable resolution (say ) would have the same dimensions as a highresolution image (). A fully connected DBN on such an image would result in a huge number of parameters making the model intractable to train effectively. Therefore, we propose to use convolution to reduce model parameters by weight sharing. However, different from typical convolutional deep learning models (e.g. [21]), we do not use any form of pooling in the hidden layers – while pooling may enhance the invariance properties for recognition, in our case, it would also lead to greater uncertainty for shape reconstruction.
The energy, , of a convolutional layer in our model can be computed as:
(1) 
where denotes each visible unit, denotes each hidden unit in a feature channel , and denotes the convolutional filter. The “” sign represents the convolution operation. In this energy definition, each visible unit is associated with a unique bias term to facilitate reconstruction, and all hidden units in the same convolution channel share the same bias term . Similar to [19]
, we also allow for a convolution stride.
A 3D shape is represented as a
voxel grid with 3 extra cells of padding in both directions to reduce the convolution border artifacts. The labels are presented as standard one of
softmax variables. The final architecture of our model is illustrated in Figure 2(a). The first layer has 48 filters of size 6 and stride 2; the second layer has 160 filters of size 5 and stride 2 (i.e., each filter has parameters); the third layer has 512 filters of size 4; each convolution filter is connected to all the feature channels in the previous layer; the fourth layer is a standard fully connected RBM with 1200 hidden units; and the fifth and final layer with 4000 hidden units takes as input a combination of multinomial label variables and Bernoulli feature variables. The top layer forms an associative memory DBN as indicated by the bidirectional arrows, while all the other layer connections are directed topdown.We first pretrain the model in a layerwise fashion followed by a generative finetuning procedure. During pretraining, the first four layers are trained using standard Contrastive Divergence
[14], while the top layer is trained more carefully using Fast Persistent Contrastive Divergence (FPCD) [32]. Once the lower layer is learned, the weights are fixed and the hidden activations are fed into the next layer as input. Our finetuning procedure is similar to wake sleep algorithm [15] except that we keep the weights tied. In the wake phase, we propagate the data bottomup and use the activations to collect the positive learning signal. In the sleep phase, we maintain a persistent chain on the topmost layer and propagate the data topdown to collect the negative learning signal. This finetuning procedure mimics the recognition and generation behavior of the model and works well in practice. We visualize some of the learned filters in Figure 2(b).During pretraining of the first layer, we collect learning signal only in receptive fields which are nonempty. Because of the nature of the data, empty spaces occupy a large proportion of the whole volume, which have no information for the RBM and would distract the learning. Our experiment shows that ignoring those learning signals during gradient computation results in our model learning more meaningful filters. In addition, for the first layer, we also add sparsity regularization to restrict the mean activation of the hidden units to be a small constant (following the method of [20]). During pretraining of the topmost RBM where the joint distribution of labels and highlevel abstractions are learned, we duplicate the label units 10 times to increase their significance.
After training the CDBN, the model learns the joint distribution of voxel data and object category label . Although the model is trained on complete 3D shapes, it is able to recognize objects in singleview 2.5D depth maps (e.g., from RGBD sensors). As shown in Figure 3, the 2.5D depth map is first converted into a volumetric representation where we categorize each voxel as free space, surface or occluded, depending on whether it is in front of, on, or behind the visible surface (i.e., the depth value) from the depth map. The free space and surface voxels are considered to be observed, and the occluded voxels are regarded as missing data. The test data is represented by , where refers to the observed free space and surface voxels, while
refers to the unknown voxels. Recognizing the object category involves estimating
.We approximate the posterior distribution by Gibbs sampling. The sampling procedure is as follows. We first initialize to a random value and propagate the data bottom up to sample for a label from . Then the high level signal is propagated down to sample for voxels . We clamp the observed voxels in this sample and do another bottom up pass. 50 iterations of updown sampling are sufficient to get a shape completion , and its corresponding label . The above procedure is run in parallel for a large number of particles resulting in a variety of completion results corresponding to potentially different classes. The final category label corresponds to the most frequently sampled class.
Object recognition from a singleview can sometimes be challenging, both for humans and computers. However, if an observer is allowed to view the object from another view point when recognition fails from the first view point, we may be able to significantly reduce the recognition uncertainty. Given the current view, our model is able to predict which next view would be optimal for discriminating the object category.
The inputs to our nextbestview system are observed voxels of an unknown object captured by a depth camera from a single view, and a finite list of nextview candidates representing the camera rotation and translation in 3D. An algorithm chooses the nextview from the list that has the highest potential to reduce the recognition uncertainty. Note that during this view planning process, we do not observe any new data, and hence there is no improvement on the confidence of .
The original recognition uncertainty, , is given by the entropy of conditioned on the observed :
(2) 
where the conditional probability can be approximated as before by sampling from and marginalizing .
When the camera is moved to another view , some of the previously unobserved voxels may become observed based on its actual shape. Different views will result in different visibility of these unobserved voxels . A view with the potential to see distinctive parts of objects (e.g. arms of chairs) may be a better next view. However, since the actual shape is partially unknown^{2}^{2}2If the 3D shape is fully observed, adding more views will not help to reduce the recognition uncertainty in any algorithm purely based on 3D shapes, including our 3D ShapeNets., we will hallucinate that region from our model. As shown in Figure 4, conditioning on , we can sample many shapes to generate hypotheses of the actual shape, and then render each hypothesis to obtain the depth maps observed from different views, . In this way, we can simulate the new depth maps for different views on different samples and compute the potential reduction in recognition uncertainty.
Mathematically, let denote the new observed voxels (both free space and surface) in the next view . We have , and they are unknown variables that will be marginalized in the following equation. Then the potential recognition uncertainty for is measured by this conditional entropy,
(3) 
The above conditional entropy could be calculated by first sampling enough from , doing the 3D rendering to obtain 2.5D depth map in order to get from , and then taking each to calculate as before.
According to information theory, the reduction of entropy is the mutual information between and conditioned on . This meets our intuition that observing more data will always potentially reduce the uncertainty. With this definition, our view planning algorithm is to simply choose the view that maximizes this mutual information,
(4) 
Our view planning scheme can naturally be extended to a sequence of view planning steps. After deciding the best candidate to move for the first frame, we physically move the camera there and capture the other object surface from that view. The object surfaces from all previous views are merged together as our new observation , allowing us to run our view planning scheme again.
Training a deep 3D shape representation that captures intraclass variance requires a large collection of 3D shapes. Previous CAD datasets (e.g., [28]) are limited both in the variety of categories and the number of examples per category. Therefore, we construct ModelNet, a largescale 3D CAD model dataset.
To construct ModelNet, we downloaded 3D CAD models from 3D Warehouse, and Yobi3D search engine indexing 261 CAD model websites. We query common object categories from the SUN database [33] that contain no less than 20 object instances per category, removing those with too few search results, resulting in a total of 660 categories. We also include models from the Princeton Shape Benchmark [28]. After downloading, we remove miscategorized models using Amazon Mechanical Turk. Turkers are shown a sequence of thumbnails of the models and answer “Yes” or “No” as to whether the category label matches the model. The authors then manually checked each 3D model and removed irrelevant objects from each CAD model (e.g, floor, thumbnail image, person standing next to the object, etc) so that each mesh model contains only one object belonging to the labeled category. We also discarded unrealistic (overly simplified models or those only containing images of the object) and duplicate models. Compared to [28], which consists of 6670 models in 161 categories, our new dataset is 22 times larger containing 151,128 3D CAD models belonging to 660 unique object categories. Examples of major categories and dataset statistics are shown in Figure 5.
We choose 40 common object categories from ModelNet with 100 unique CAD models per category. We then augment the data by rotating each model every 30 degrees along the gravity direction (i.e., 12 poses per model) resulting in models in arbitrary poses. Pretraining and finetuning each took about two days on a desktop with one Intel XEON E52690 CPU and one NVIDIA K40c GPU. Figure 6 shows some shapes sampled from our trained model.
Deep learning has been widely used as a feature extraction technique. Here, we are also interested in how well the features learned from 3D ShapeNets compare with other stateoftheart 3D mesh features. We discriminatively finetune 3D ShapeNets by replacing the top layer with class labels and use the 5th layer as features. For comparison, we choose Light Field descriptor
[8] (LFD, 4,700 dimensions) and Spherical Harmonic descriptor [18] (SPH, 544 dimensions), which performed best among all descriptors [28].


10 classes  SPH [18]  LFD [8]  Ours 


classification  79.79 %  79.87 %  83.54% 
retrieval AUC  45.97%  51.70%  69.28% 
retrieval MAP  44.05%  49.82%  68.26% 


40 classes  SPH [18]  LFD [8]  Ours 


classification  68.23%  75.47%  77.32% 
retrieval AUC  34.47%  42.04%  49.94% 
retrieval MAP  33.26%  40.91%  49.23% 

We conduct 3D classification and retrieval experiments to evaluate our features. Of the 48,000 CAD models (with rotation enlargement), 38,400 are used for training and 9,600 for testing. We also report a smaller scale result on a 10category subset (corresponding to NYU RGBD dataset [23]
) of the 40category data. For classification, we train a linear SVM to classify meshes using each of the features mentioned above, and use average category accuracy to evaluate the performance.
For retrieval, we use distance to measure the similarity of the shapes between each pair of testing samples. Given a query from the test set, a ranked list of the remaining test data is returned according to the similarity measure^{3}^{3}3For our feature and SPH we use the norm, and for LFD we use the distance measure from [8].. We evaluate retrieval algorithms using two metrics: (1) mean area under precisionrecall curve (AUC) for all the testing queries^{4}^{4}4We interpolate each precisionrecall curve.; (2) mean average precision (MAP) where AP is defined as the average precision each time a positive sample is returned.
We summarize the results in Table 1 and Figure 7. Since both of the baseline mesh features (LFD and SPH) are rotation invariant, from the performance we have achieved, we believe 3D ShapeNets must have learned this invariance during feature learning. Despite using a significantly lower resolution mesh as compared to the baseline descriptors, 3D ShapeNets outperforms them by a large margin. This demonstrates that our 3D deep learning model can learn better features from 3D data automatically.
To evaluate 3D ShapeNets for 2.5D depthbased object recognition task, we set up an experiment on the NYU RGBD dataset with Kinect depth maps [23]. We select 10 object categories from ModelNet that overlap with the NYU dataset. This results in 4,899 unique CAD models for training 3D ShapeNets.
We create each testing example by cropping the 3D point cloud from the 3D bounding boxes. The segmentation mask is used to remove outlier depth in the bounding box. Then we directly apply our model trained on CAD models to the NYU dataset. This is absolutely nontrivial because the statistics of real world depth are significantly different from the synthetic CAD models used for training. In Figure
9, we visualize the successful recognitions and reconstructions. Note that 3D ShapeNets is even able to partially reconstruct the “monitor” despite the bad scanning caused by the reflection problem. To further boost recognition performance, we discriminatively finetune our model on the NYU dataset using back propagation. By simply assigning invisible voxels as 0 (i.e. considering occluded voxels as free space and only representing the shape as the voxels on the 3D surface) and rotating training examples every 30 degrees, finetuning works reasonably well in practice.As a baseline approach, we use nearestneighbor matching in our low resolution voxel space. Testing depth maps are converted to voxel representation and compared with each of the training samples. As a more sophisticated high resolution baseline, we match the testing point cloud to each of our 3D mesh models using Iterated Closest Point method [4] and use the top 10 matches to vote for the labels. We also compare our result with [29] which is the stateoftheart deep learning model applied to RGBD data. To train and test their model, 2D bounding boxes are obtained by projecting the 3D bounding box to the image plane, and object segmentations are also used to extract features. 1,390 instances are used to train the algorithm of [29] and perform our discriminative finetuning, while the remaining 495 instances are used for testing all five methods. Table 3 summarizes the recognition results. Using only depth without color, our finetuned 3D ShapeNets outperforms all other approaches with or without color by a significant margin.


bathtub  bed  chair  desk  dresser  monitor  nightstand  sofa  table  toilet  all  


[29] Depth  0.000  0.729  0.806  0.100  0.466  0.222  0.343  0.481  0.415  0.200  0.376 
NN  0.429  0.446  0.395  0.176  0.467  0.333  0.188  0.458  0.455  0.400  0.374 
ICP  0.571  0.608  0.194  0.375  0.733  0.389  0.438  0.349  0.052  1.000  0.471 
3D ShapeNets  0.142  0.500  0.685  0.100  0.366  0.500  0.719  0.277  0.377  0.700  0.437 
3D ShapeNets finetuned  0.857  0.703  0.919  0.300  0.500  0.500  0.625  0.735  0.247  0.400  0.579 


[29] RGB  0.142  0.743  0.766  0.150  0.266  0.166  0.218  0.313  0.376  0.200  0.334 
[29] RGBD  0.000  0.743  0.693  0.175  0.466  0.388  0.468  0.602  0.441  0.500  0.448 



bathtub  bed  chair  desk  dresser  monitor  nightstand  sofa  table  toilet  all  


Ours  0.80  1.00  0.85  0.50  0.45  0.85  0.75  0.85  0.95  1.00  0.80 
Max Visibility  0.85  0.85  0.85  0.50  0.45  0.85  0.75  0.85  0.90  0.95  0.78 
Furthest Away  0.65  0.85  0.75  0.55  0.25  0.85  0.65  0.50  1.00  0.85  0.69 
Random Selection  0.60  0.80  0.75  0.50  0.45  0.90  0.70  0.65  0.90  0.90  0.72 

For our view planning strategy, computation of the term is critical. When the observation is ambiguous, samples drawn from should come from a variety of different categories. When the observation is rich, samples should be limited to very few categories. Since is the surface of the completions, we could just test the shape completion performance . In Figure 8, our results give reasonable shapes across different categories. We also match the nearest neighbor in the training set to show that our algorithm is not just memorizing the shape and it can generalize well.
To evaluate our view planning strategy, we use CAD models from the test set to create synthetic renderings of depth maps. We evaluate the accuracy by running our 3D ShapeNets model on the integration depth maps of both the first view and the selected second view. A good viewplanning strategy should result in a better recognition accuracy. Note that nextbestview selection is always coupled with the recognition algorithm. We prepare three baseline methods for comparison : (1) random selection among the candidate views; (2) choose the view with the highest new visibility (yellow voxels, NBV for reconstruction); (3) choose the view which is farthest away from the previous view (based on camera center distance). In our experiment, we generate 8 view candidates randomly distributed on the sphere of the object, pointing to the region near the object center and, we randomly choose 200 test examples (20 per category) from our testing set. Table 3 reports the recognition accuracy of different view planning strategies with the same recognition 3D ShapeNets. We observe that our entropy based method outperforms all other strategies.
To study 3D shape representation for objects, we propose a convolutional deep belief network to represent a geometric 3D shape as a probability distribution of binary variables on a 3D voxel grid. Our model can jointly recognize and reconstruct objects from a singleview 2.5D depth map (e.g. from popular RGBD sensors). To train this 3D deep learning model, we construct ModelNet, a largescale 3D CAD model dataset. Our model significantly outperforms existing approaches on a variety of recognition tasks, and it is also a promising approach for nextbestview planning. All source code and data set are available at our project website.
This work is supported by gift funds from Intel Corporation and Project X grant to the Princeton Vision Group, and a hardware donation from NVIDIA Corporation. Z.W. is also partially supported by Hong Kong RGC Fellowship. We thank Thomas Funkhouser, Derek Hoiem, Alexei A. Efros, Andrew Owens, Antonio Torralba, Siddhartha Chaudhuri, and Szymon Rusinkiewicz for valuable discussion.
On visual similarity based 3d model retrieval.
In Computer graphics forum, 2003.The shape boltzmann machine: a strong model of object shape.
In CVPR, 2012.SUN database: Largescale scene recognition from abbey to zoo.
In CVPR, 2010.
Comments
There are no comments yet.