Building Statistical Shape Spaces
for 3D Human Modeling
Abstract
Statistical models of 3D human shape and pose learned from scan databases have developed into valuable tools to solve a variety of vision and graphics problems. Unfortunately, most publicly available models are of limited expressiveness as they were learned on very small databases that hardly reflect the true variety in human body shapes. In this paper, we contribute by rebuilding a widely used statistical body representation from the largest commercially available scan database, and making the resulting model available to the community (visit http://humanshape.mpiinf.mpg.de). As preprocessing several thousand scans for learning the model is a challenge in itself, we contribute by developing robust best practice solutions for scan alignment that quantitatively lead to the best learned models. We make implementations of these preprocessing steps also publicly available. We extensively evaluate the improved accuracy and generality of our new model, and show its improved performance for human body reconstruction from sparse input data.
keywords:
statistical human body model, nonrigid template fitting1 Introduction
Statistical human shape models represent variations in human physique and pose using lowdimensional parameter spaces, and are valuable tools to solve difficult vision and graphics problems, e.g. in pose tracking or animation. Despite significant progress in modeling the statistics of the complete 3D human shape and pose AllenTSO03 (); AnguelovSCA05 (); Guan2012 (); ChenTBH13 (); Neophytou2013 (); HaslerBM2009 (); Jain:2010:MovieReshape () only few publicly available statistical 3D body shape spaces exist HaslerBM2009 (); Jain:2010:MovieReshape (). Further on, the public models are often learned on only small datasets with limited shape variations HaslerBM2009 (). The reason is a lack of large representative public datasets and the significant effort required to process and align raw laser scans prior to learning a statistical shape space.
This paper contributes by systematically constructing a model of 3D human shape and pose from the largest commercially available dataset of 3D laser scans RobinetteTCP99 () and making it publicly available to the research community (Section 2). Our model is based on a simplified and efficient variant of the SCAPE model AnguelovSCA05 () (henceforth termed SSCAPE space) that was described by Jain et al. Jain:2010:MovieReshape () and used for different applications in computer vision and graphics Jain:2010:MovieReshape (); pishchulin11bmvc (); pishchulin12cvpr (); HeltenPAE13 (); conf/cvpr/MundermannCA07 (), but was never learned from such a complete dataset. This compact shape space learns a probability distribution from a dataset of 3D human laser scans. It models variations due to changes in identity using a principal component analysis (PCA) space, and variations due to pose using a skeletonbased surface skinning approach. This representation makes the model versatile and computationally efficient.
Prior to statistical analysis, the human scans have to be processed and aligned to establish correspondence. We contribute by evaluating different variants of the stateoftheart techniques for nonrigid template fitting and posture normalization to process the raw data AllenTSO03 (); HaslerBM2009 (); WuhrerPIS12 (); Neophytou2013 (). Our findings are not entirely new methods, but best practices and specific solutions for automatic preprocessing of large scan databases for learning the SSCAPE model in the best way (Section 3). First, shape and posture fitting of an initial shape model to a raw scan prior to nonrigid deformation considerably improves the results. Second, multiple passes over the dataset improve initialization and thus increase the overall fitting accuracy and statistical model qualities. Third, posture normalization prior to shape space learning leads to much better generalization and specificity.
The main contribution of our work is a set of SSCAPE spaces learned from the largest database that is currently commercially available RobinetteTCP99 (). The differences in our SSCAPE spaces stem from differences in the registration and prealignment of the human body scans. We evaluate different data processing techniques in Section 4 and the resulting shape spaces in Section 5. Finally, in Section 6 we compare our SSCAPE spaces to the state of the art SSCAPE space learned from a publicly available database HaslerBM2009 () for the application of reconstructing full 3D body models from partial depth data. Experimental evaluation clearly demonstrates the advantages of our more expressive shape models in terms of shape space quality and performance on the task of reconstructing 3D human body shapes from partial depth observations.
We release the new shape spaces with code to (1) preprocess raw scans and (2) fit a shape space to a raw scan for public usage. We believe this contribution is required for future development in human body modeling. Visit http://humanshape.mpiinf.mpg.de to download code and models.
1.1 Related work
Datasets
Several datasets have been collected to analyze populations of 3D human bodies. Many publicly available research datasets allow for the analysis of shape and posture variations jointly; unfortunately they feature data of at most 100 individuals AnguelovSCA05 (); HaslerBM2009 (); Bogo:CVPR:2014 (), which limits the range of shape variations. We therefore use CAESAR database RobinetteTCP99 (), the largest commercially available dataset to date that contains 3D scans of over American and European subjects in a standard pose, as it represents much richer sample of the human physique.
Statistical shape spaces of 3D human bodies
Building statistical shape spaces of human bodies is challenging, as there is strong and intertwined 3D shape and posture variability yielding a complex function of multiple correlated shape and posture parameters. Methods to learn this shape space usually follow one of two routes. The first group of methods learn shape and posturerelated deformations separately and combine them afterwards AnguelovSCA05 (); Guan2012 (); ChenTBH13 (); Jain:2010:MovieReshape (); Neophytou2013 (); Loper2015 (). These methods are inspired by SCAPE model AnguelovSCA05 () that couples a shape space learned from variations in body shape with a posture space learned from deformations of a single subject. This method has recently been enhanced to capture deformations related to breathing Tsoli2014 () and dynamic motions PonsMoll2015 (). Most SCAPElike models use a set of transformations per triangle to encode shape variations in a shape space. Hence, to convert between the vertex coordinates of a processed scan and its representation in shape space, a computationally demanding optimization problem needs to be solved. To overcome this difficulty, a simplified version of SCAPE model (SSCAPE) was proposed Jain:2010:MovieReshape (). SSCAPE operates on vertex coordinates directly and models pose variation with an efficient skeletonbased surface skinning approach Jain:2010:MovieReshape (); HeltenPAE13 (); conf/cvpr/MundermannCA07 (). Recently, two alternative multilinear shape spaces have been proposed that also operate directly on vertex coordinates Neophytou2013 (); Loper2015 ().
Another group of methods intends to perform simultaneous analysis of shape and posture variations AllenLAC06 (); HaslerBM2009 (). These methods learn skinning weights for corrective enveloping of posturerelated shape variations, which allows to explore both shape and posture variations using a single shape space. Furthermore, it allows for realistic muscle bulging as shape and posture are correlated NeumannEG2013 (). It has been shown, however, that for many applications in computer vision and graphics this level of detail is not required and simpler and computationally more efficient shape spaces can be used Jain:2010:MovieReshape (); HeltenPAE13 (); pishchulin11bmvc (); pishchulin12cvpr ().
Mesh registration
Mesh registration is performed on the scans to bring them in correspondence for statistical analysis. Two surveys vanKaickASO11 (); TamRO313 () review such techniques, and a full review is beyond the scope of this paper. Allen et al. AllenTSO03 () use nonrigid template fitting to compute correspondences between human body shapes in similar posture. This technique has been extended to work for varying postures AnguelovSCA05 (); AllenLAC06 (); HaslerBM2009 () and in scenarios where no landmarks are available WuhrerLFP11 (). In this work, we evaluate a nonrigid template fitting approach inspired by Allen et al. AllenTSO03 ().
Applications
Statistical spaces of human body shape and posture are applicable in many areas including computer vision, computer graphics, and ergonomic design; our new model that was learned on a large commercially available dataset is beneficial in each of these applications. Statistical shape spaces have been used to predict body shapes from partial data, such as image sequences and depth images Seo2006 (); balan07CVPR (); Sigal2007 (); Xi2007a (); GuanEHS09 (); HaslerMPA10 (); WeissH3D11 (); Boisvert2013 (); HeltenPAE13 () and semantic parameters Seo2003 (); AllenTSO03 (); Chu2010 (); Baek2012 (); WuhrerE3D13 (); Rupprecht3DS13 (). Furthermore, they have been used to estimate body shapes from images BalanTNT08 () and 3D scans HaslerEBS09 (); Wuhrer14CVIU () of dressed subjects. Given a 3D body shape, statistical shape spaces can be used to modify input images ZhouPRO10 () or videos Jain:2010:MovieReshape (), to automatically generate training sets for people detection pishchulin11bmvc (); pishchulin12cvpr (), or to simulate clothing on people Guan2012 ().
2 Statistical modeling with SCAPE
We briefly recap the efficient version of the SCAPE model Jain:2010:MovieReshape () we build on and discuss its differences to the original SCAPE model AnguelovSCA05 () in more detail. For learning the model, both methods assume that a template mesh T containing vertices has been deformed to each raw scan in a database. All scans of the database are assumed to be rigidly aligned, e.g. by Procrustes Analysis goodall91jrss ().
2.1 Original SCAPE model
In the original SCAPE model, the transformation of each triangle of T is modeled as combination of three linear transformations and controlling posture, and controlling body shape. Index indicates one particular scan T is fitted to. Fitting result after rigid alignment with T is denoted as instance mesh .
Shape deformations encode pertriangle deformations that can be applied to change the body shape of the person in the same standard posture. A lowdimensional space of plausible shape deformations is computed by performing PCA on the training dataset captured in standard posture.
To represent posture changes, two transformations are used: represents the posture of the person as rotation induced by the deformation of an underlying rigid skeleton, and encodes the individual deformations of each triangle that originates from varying body shape or nonrigid posture dependent surface deformations such as muscle bulging. Computing for each triangle separately is an underconstrained problem. Therefore, smoothing is applied such that of neighboring triangles become dependent. Finally, the dimensionality of the transformations and is reduced with the help of a kinematic chain model.
In this way, SCAPE obtains a flexible model that covers a wide range of possible shape and posture deformations. However, as the model does not explicitly encode vertex positions, a computationally expensive optimization problem needs to be solved in order to reconstruct the mesh surface.
2.2 Simplified SCAPE (SScape) space
The aforementioned computational overhead is often prohibitive in applications where speed is more important than the overall reconstruction quality, or when many samples need to be drawn from the shape space. SSCAPE space Jain:2010:MovieReshape () reconstructs vertex positions in a given posture and shape without need of solving a Poisson system. To learn the model, only laser scans in a standard posture are used. Meshes are used to learn a PCA model that represents each shape using a parameter vector and can generate new models (represented in homogeneous coordinates) with body shape in posture as
(1) 
with the dimension of the PCA space is the matrix computed using PCA and is the mean body shape of the training database.
This shape space only covers variations in body shape but not in posture. To enable the latter an articulated skeleton is fitted to the average human shape and linear blend skinning weights are used to attach surface to bones. This allows to deform a body with fixed shape into an arbitrary posture as
(2) 
where is the homogeneous coordinate of the th vertex of , is the number of bones used for the rigging, is the transformation of the th bone, and are the rigging weights. We use the rigging and skeleton consisting of bones proposed by Jain et al. Jain:2010:MovieReshape (). The skeleton is controlled by pose parameters corresponding to rigid transformations, joint angles and scale.
For reconstructing a model of shape in skeleton posture , the method first calculates a personalized mesh using , and subsequently applies linear blend skinning to the personalized mesh to obtain the final mesh . This can be expressed in matrix notation as
(3) 
where is a blockdiagonal matrix containing the pervertex transformations. While decoupling of shape and posture modeling by SSCAPE results in lower level of details (e.g. posturespecific deformations such as muscle bulging may be missing), it leads much faster reconstruction speed, especially when the personalized mesh and skeleton can be precomputed. We argue that in many applications speed may be much more important than the overall reconstruction quality and build on this simple and efficient shape space in this work.
3 Data processing
This section describes our preprocessing procedure that allows to establish correspondences between raw laser scans. We demonstrate bestpractice solutions for nonrigid template fitting, effective initialization strategies, introduce novel humanintheloop bootstrapping approach that allows to improve the correspondences, and finally explore postures normalization strategies. Tools to reproduce these steps are made publicly available.
3.1 Nonrigid template fitting
Our method to fit a human shape template T to a human scan is inspired by Allen et al. AllenTSO03 (). In nonrigid template fitting (henceforth abbreviated NRD), each vertex of T is transformed by a affine matrix , which allows for twelve degrees of freedom during the transformation. The aim is to find a set of matrices that align vertices of the deformed template M to the corresponding points of in the best possible way. The fitting is done by minimizing a combination of data, smoothness and landmark errors.
Data term
The data term requires each vertex of the transformed template to be as close as possible to its corresponding vertex of , and takes the form
(4) 
where weights the error contribution of each vertex, denotes the Frobenius norm, and is a closest compatible point in . If surface normals of closest points are less than apart and the distance between the points is less than mm, we set to 1, otherwise to 0.
Smoothness term
Fitting using only may lead to situations where neighboring vertices of T match to disparate vertices in . To enforce smooth surface deformations we use a smoothness term that requires affine transformations applied to connected vertices to be similar, i.e.
(5) 
Landmark term
Although using and would suffice to fit two surfaces that are close to each other, the optimization may get stuck in a local minimum when T and are far apart. A remedy is to identify a set of points on T corresponding to known anthropometric landmarks on . In each CAESAR scan these are obtained by placing markers on each subject prior to scanning. Our landmark term penalizes misalignments between landmark locations
(6) 
where is the landmark index on T, and is the landmark point on . Although there are only 64 landmarks compared to the total number of vertices, good landmark fitting is enough to get the deformed surface of T close to and avoid local convergence.
Combined energy
The three terms are combined into a single objective
(7) 
For optimization we use LBFGSB zbMATH01235219 (). We vary the weights , and according to the following empirically found schedule. We first perform a single iteration of optimization without data term by setting , , , which allows to bring the surfaces into a rough correspondence. We then allow the data term to contribute by setting , , . In addition, we relax smoothness and landmark weights after each iteration of fitting to and , thus allowing the data term to dominate. This is repeated until . Reducing increases the flexibility of deformation and allows T to better reproduce fine details, while reducing is necessary due to unreliable placement of landmarks in some scans.
3.2 Initialization
For nonrigid template fitting to succeed, T should be prealigned to . We explore two initialization strategies.
A first standard way to initialize NRD is to use a static template with annotated landmarks. Corresponding landmarks are then used to rigidly align to T.
A second way to initialize the fitting is to start with a SSCAPE space that was learned from a small registered dataset. Fitting the shape space to a scan is achieved by finding shape and posture parameters and such that (see Eq. 3) is close to . To this end, (Eq. 7) is minimized with respect to and . To minimize depending on and , we use the vertices of , and set the pervertex deformations to the identity. That is, the deformation of the body shape is exclusively controlled by the parameters and . As in this case neighboring vertices do not move independently, the term is not required, and we set .
To find a good local minimum, good initialization is required. We found a twostep optimization approach to work well in practice. First, we set and and optimize with respect to while fixing , which fits the posture of to with the help of landmarks. Second, we set and and optimize with respect to and iteratively. For increased efficiency, each iteration optimizes with respect to in a first step and with respect to in a second step. After each iteration, the set is recomputed. This iterative procedure is repeated until does not change significantly. Iterative interior point method is used for optimization.
3.3 Bootstrapping
In many cases, even after nonrigid template fitting (NRD), fitted mesh M is far from the target human scan . Learning from registered scans with a high fitting error may capture unrealistic shape deformations. We thus propose the following humanintheloop bootstrapping learning process: after each fitting iteration we visually examine each registered scan, discard registered scans of low quality, and learn a SSCAPE space using the registered scans that passed the visual inspection; learned SSCAPE space is then used during initialization of the next fitting pass and the process is repeated. This bootstrapping process is performed for multiple iterations until nearly all registered scans pass the visual inspection. Note that visual inspection is required, as low average fitting errors do not always correspond to good results, since the fitting of localized areas may be inaccurate.
3.4 Posture normalization
The SSCAPE space used in this study decouples learning of shape and posture variations and learns shape variations via PCA on the registered scans captured in a standard posture. However even standard postures may still contain slight posture variations, mostly due to movements of arms. Thus PCA may learn global shape variations due to variation in posture. In order to address this issue we perform posture normalization of the registered scans based on two approaches WuhrerPIS12 (); Neophytou2013 (), as explained in the following.
Wuhrer et al. WuhrerPIS12 () factor out variations due to posture changes by performing PCA on localized Laplacian coordinates. While this approach leads to better shape spaces, it is difficult to directly apply this approach to the SSCAPE spaces learned using Cartesian coordinates. We therefore compute a posturenormalized version of each fitted mesh in the following way: we start with a mean shape computed over all and use WuhrerPIS12 () to optimize the localized Laplacian coordinates of to be as close as possible to . This leads to fittings that have the body shape of in the common normalized posture of .
Neophytou and Hilton Neophytou2013 () normalize the posture of each processed scan using a skeleton model and Laplacian surface deformation. While such normalization may introduce artifacts around joints when the posture is changed significantly, this approach is suitable to account for minor posture variations of CAESAR scans. We use this method to modify the posture of each fitted mesh .
4 Evaluation of template fitting
We now evaluate different components of our registration procedure on CAESAR dataset RobinetteTCP99 (). Each CAESAR scan contains manually placed landmarks. We exclude several landmarks located on open hands, as those are missing for our template, resulting in landmarks used for registration. Furthermore, we remove all laser scans without landmarks and corrupted scans, resulting in 4308 scans.
4.1 Implementation details
Nonrigid template fitting requires a human shape template as input, and the initialization procedure requires an initial shape space. We use registered scans of 111 individuals in neutral posture of the MPI Human Shape dataset HaslerBM2009 () to compute these initializations.
However, MPI scans have artifacts such as spiky nonsmooth surfaces in the areas of head and neck. We smooth these areas by identifying problematic vertices and by iteratively recomputing their positions as an average position of direct neighbors. Furthermore, due to privacy reasons, head vertices of each human scan were replaced by the same dummy head, which is not representative and of low quality at the backside. We adjust the vertex compatibility criteria to compute nearest neighbors during NRD by allowing deviation of the head face normals while increasing the distance threshold to mm.
We employ the algorithm from Section 3.1 to compute correspondences for the CAESAR dataset. One inconsistency between the datasets is that the hands in the MPI Human Shape dataset are closed, while they are open in the CAESAR dataset. As remedy, we set and to zero for hand vertices in Eq. 7, thus only allowing to contribute. Prior to fitting, we subsample each CAESAR scan to have a total number of vertices that exceeds the number of vertices of T by a factor of three ( vertices in T vs. vertices in ). This provides a good tradeoff between fitting quality and computational efficiency.
4.2 Quality measure
Measuring the accuracy of surface fitting is not straightforward, as no ground truth correspondence between and T is available. We evaluate the fitting accuracy by finding the nearest neighbor in for each fitted template vertex. If this neighbor is not further than mm from its correspondence in T and its face normals do not deviate more than , the Euclidean vertextovertex distance is computed. In our experiments we report both the proportion of vertices falling below a certain threshold and the distance per vertex averaged over all fitted templates. In the following, we first show the effects of various types of initialization and weighting schemes in the NRD procedure on the fitting error. Second, we show the effect of performing multiple bootstrapping rounds.


4.3 Initialization
First, we evaluate two different initialization strategies used in our fitting procedure. We compare the results when using an average human template (NRD) to the case when additionally using the SSCAPE space learned on the MPI Human Shape dataset (SSCAPE + NRD) for initialization. We also demonstrate the effects of both nonrigid deformation schemes on the fitting accuracy and finally compare to the results when using the publicly available SSCAPE space by Jain et al. Jain:2010:MovieReshape () alone (SSCAPE).
The results are shown in Fig. 1. The total fitting error in Fig. 1(a) shows that NRD achieves good fitting results in the low error range of 0 – 10 mm, as it can produce good template fits for the areas where T is close to . However, as NRD is a modelfree method, the smooth topology of T may not be preserved during the deformation, e.g. convex surfaces of T may be deformed into nonconvex surfaces after NRD. This leads to large fitting errors for areas of T that are far from . SSCAPE + NRD uses a shape space fitting prior to NRD, which allows for a better initial alignment of T to . Note that SSCAPE + NRD results in a better fitting accuracy in the high error range of . The fitting result by SSCAPE + NRD favorably compares against using SSCAPE alone. Although SSCAPE results into deformations preserving the human body shape topology, the shape space is learned from the relatively specialized MPI Human Shape dataset containing mostly young adults and thus cannot represent all shape variations.
We also analyze the differences in the mean fitting errors per vertex in Fig. 1(b). NRD achieves good fitting results for most of the vertices. However, the arms are not fitted well due to differences in body posture of T and . Furthermore, the average fitting error is not smooth, which shows that despite using , NRD may produce nonsmooth deformations. In contrast, the result of SSCAPE + NRD is smoother and has a lower fitting error for the arms. Clearly, the average fitting error of SSCAPE is much higher, with notably worse fitting results for arms, belly and chest.
4.4 Nrd parameters
Second, we evaluate the influence of the weight relaxation during NRD on the fitting accuracy. Specifically, we compare the standard weighting scheme where weights are relaxed in each iteration (SSCAPE + NRD) to the case where the weights stay constant (SSCAPE + NRD CW). Fig. 1(a) shows that the total fitting error of SSCAPE + NRD is lower than SSCAPE + NRD CW. This is because SSCAPE + NRD CW enforces higher localized rigidity by keeping weights constantly high, while SSCAPE + NRD relaxes the weights so that T can fit more accurately to . This explanation is supported by consistently higher pervertex mean fitting errors in case of SSCAPE + NRD CW compared to SSCAPE + NRD, as shown in Fig. 1(b). The highest differences are in the areas of high body shape variability, such as belly and chest. Different weight reduction schemes such as , and , lead to better fitting accuracy compared to constant weights, with the latter scheme achieving slightly better results and faster convergence rates. We thus use the proposed weight reduction scheme in the following.


4.5 Bootstrapping
Third, we evaluate the fitting accuracy before and after performing multiple rounds of bootstrapping. To that end, we use the output of SSCAPE + NRD (iteration 0) to learn a new statistical shape space, which is in turn used to initialize NRD during the second pass over the data (iteration 1). This process is repeated for five passes. The number of registered scans that passed the visual inspection after each round is , , , and , respectively. This results show that bootstrapping allows to register and thus to learn from an increasing number of scans. Fitting results are shown in Fig. 2. The closeup shows that although the overall fitting accuracy before and after bootstrapping is similar, bootstrapping allows to slightly improve the fitting accuracy in the range of . Fitting results after three passes over the dataset (iteration 2) are slightly better compared to the initial fitting (iteration 0), and the accuracy is further increased after five passes (iteration 4). Fig. 2 (b) shows sample fitting results before and after several bootstrapping rounds. Largest improvements are achieved for the belly and chest  areas with large shape variability. The fitting improves with an increasing number of bootstrapping rounds. We use the fitting results after five passes (iteration 4) to learn the SSCAPE space used in the following.
5 Evaluation of statistical shape space
In this section, we evaluate the SSCAPE space using the statistical quality measures of generalization and specificity Styner2003 ().
5.1 Quality measure
We use two complementary measures of shape statistics. Generalization evaluates the ability of a shape space to represent unseen instances of the object class. Good generalization means the shape space is capable of learning the characteristics of an object class from a limited number of training samples, poor generalization indicates overfitting of the training set. Generalization is measured using leaveoneout cross reconstruction of training samples, i.e. the shape space is learned using all but one training sample and the resulting shape space is fitted to the excluded sample. The fitting error is measured using the mean vertextovertex Euclidean distance. Generalization is reported as mean fitting error averaged over the complete set of trials, and plotted as a function of the number of shape space parameters. It is expected that the mean error decreases until convergence as the number of shape space parameters increases.
Specificity measures the ability of a shape space to generate instances of the object class that are similar to the training samples. The specificity test is performed by generating a set of instances randomly drawn from the learned shape space and by comparing them to the training samples. The error is measured as average distance of the generated instances to their nearest neighbors in the training set. It is expected that the mean distance increases until convergence with increasing number of shape space parameters. We follow Styner et al. Styner2003 () and generate 10,000 random samples.
Generalization 


Specificity 

(a) bootstrapping  (b) # training samples  (c) posture normalization 
5.2 Bootstrapping
We evaluate the influence of bootstrapping on the quality of the statistical shape space by comparing models obtained after zero, one, two and four iterations of bootstrapping. The geometry of the training samples changes in each bootstrapping round, which makes the generalization and specificity results incomparable across different shape spaces. We thus use the training samples obtained after four iterations of bootstrapping as “ground truth”, i.e., the reconstruction error of generalization and the nearest neighbor distance of specificity for each shape space is computed w.r.t. fitting results after four bootstrapping rounds. This allows for a fair comparison across different statistical shape spaces.
The results are shown in Fig. 3(a). Generalization error is already low after a single iteration of bootstrapping because after one iteration, the shape space is learned from a significantly larger number of training samples, thereby using samples with higher shape variation that were discarded in the iteration. The following rounds of bootstrapping have little influence on generalization and specificity, with the shape space after four iterations resulting in a slightly lower specificity error than for previous iterations for a small number of shape parameters.
5.3 Number of training samples
To evaluate the influence of the number of training samples, we vary the number of samples obtained after four bootstrapping iterations. Specifically, we consider subsets of , , and () training samples. To compute a shape space, the desired number of training shapes are sampled from the entire set of training samples. For generalization, we crossevaluate on all training samples by leaving one sample out and by sampling the desired number of training shapes from the remaining samples. For specificity, we compute the nearestneighbor distances to all training samples to find the closest sample.
The results are shown in Fig. 3(b). The shape space learned from the smallest number of samples performs worst. Increasing the number of samples consistently improves the performance with the best results achieved when using the maximum number. Both generalization and specificity error reduction is most pronounced when increasing the number of samples from to . Further increasing the number of samples to affects specificity much stronger than generalization. This shows that the shape space learned from only samples generalizes well, while its generative qualities are poor. Increasing the number of samples from to only slightly reduces both generalization and specificity errors, which shows that a highquality statistical shape space can be learned from samples.
5.4 Posture normalization
Finally, we evaluate the generalization and specificity of the shape space obtained when performing posture normalization using the methods of Wuhrer et al. WuhrerPIS12 () (WSX) and Neophytou and Hilton Neophytou2013 () (NH). The results are shown in Fig. 3 (c). Posture normalization significantly improves generalization and specificity, with WSX achieving the best result. The reduction of the average fitting error in case of generalization is highest for a low number of shape parameters. This is because both WSX and NH lead to shape spaces that are more compact compared to the shape space obtained with unnormalized training shapes. Additionally, both posturenormalized shape spaces exhibit much better specificity. Compared to the shape space trained before posture normalization, randomly generated samples from both shape spaces trained after WSX and NH exhibit less variation in posture and are thus more similar to their corresponding posturenormalized training samples.
Finally, we qualitatively examine the first five PCA components learned by the following SSCAPE spaces: the current stateoftheart shape space SSCAPE Jain:2010:MovieReshape (), our shape space without posture normalization and with posture normalization using WSX and NH. The results are shown in Fig. 4. Major modes of shape variation by SSCAPE (row 1) are affected by global and local posturerelated deformations, such as moving of arms or tilting the body. In contrast, the principal modes of variation by our shape space (row 2) are mostly due to shape changes, which is achieved due to better template fitting procedure and a more representative training set. However, small posture variations are still part of the learned shape space. Performing posture normalization of the training samples prior to learning the shape space completely factors out changes due to posture, as can be seen in the major principal components of WSX (row 3) and NH (row 4).
SSCAPE Jain:2010:MovieReshape () 


ours 

ours+ WSX 

ours+ NH 

st. dev.  
PCA id  I  II  III  IV  V 
6 Human body reconstruction
Finally, we evaluate our improved SSCAPE spaces on the task of estimating human body shape from sparse and noisy visual input. We follow the approach by Helten et al. HeltenPAE13 () to estimate the body shape of a person from two sequentially taken front and back depth images. First, body shape and posture are fitted independently to each depth image. Second, the obtained results are used as initialization of a method that jointly optimizes over shape and independently optimizes over posture parameters. This optimization strategy is used because the shape in both depth scans is of the same person, but the pose may differ.
6.1 Dataset and experimental setup
We use a publicly available dataset HeltenPAE13 () containing Kinect body scans of three males and three females. Examples of the Kinect scans are shown in Fig. 6(a). For each subject, a highresolution laser scan was captured to determine “ground truth” body shape. Following the evaluation protocol of Helten et al. HeltenPAE13 () we first fit a shape space to the depth data, then fit shape space to the ground truth scan, and finally compute the fitting error as a vertextovertex Euclidean distance between the vertices of the depthfitted mesh and the ground truthfitted mesh. As the required landmarks are not available for this dataset, we manually placed landmarks on each depth and laser scan.
6.2 Quantitative evaluation
For quantitative evaluation, we compare the following four shape models presented above: the current stateoftheart shape space Jain:2010:MovieReshape (), our shape space without posture normalization and with posture normalization using WSX and NH. In our experiments, we vary the number of shape space parameters and the number of training samples. To evaluate the fitting accuracy, we report the proportion of vertices below a certain threshold.
# train  # PCA coefficients  

sampl.  
20  30  40  50  
all 

1000 

100 

50 
The results are shown in Fig. 5, where the number of shape space parameters varies in the columns and the number of training samples varies in the rows. In all cases our SSCAPE spaces learned from the CAESAR dataset significantly outperform the shape space by Jain et al., which is learned from the far less representative MPI Human Shape dataset. Our models achieve good fitting accuracy when using as few as shape parameters, and the performance stays stable when increasing the number of shape parameters up to (first row). In contrast, the performance of the shape space by Jain et al. drops, possibly due to overfitting to unrealistic shape deformations in noisy depth data. Interestingly, better performance by our models is evident even in the case when all models are learned from the same number of training samples (third and fourth rows). This shows that the CAESAR data has higher shape variability than the MPI Human Shape data. In the majority of cases, the shape space learned from the posturenormalized samples with NH outperforms the shape space learned from samples without posture normalization. This shows that the posture normalization method of Neophytou and Hilton Neophytou2013 () helps to improve the accuracy of fitting to noisy depth data. Surprisingly, the shape space learned from samples without posture normalization outperforms the shape space learned from the posturenormalized samples with WSX in most cases. Overall, the quantitative results show the advantages of our approach of building SSCAPE spaces learned from a large representative set of training samples with additional posture normalization.
Female 1 


Male 1 

front/back scans  SSCAPE Jain:2010:MovieReshape ()  ours  ours  ours  
+WSX  +NH 
6.3 Qualitative evaluation
To qualitatively evaluate the fitting, we visualize the pervertex fitting errors. We consider the SSCAPE spaces learned from all available training samples and use shape space parameters. For visualization we choose two subjects, male and female, where the differences among the shape spaces are most pronounced.
Results are shown in Fig. 6. Our shape spaces better fit the data, in particular in the areas of belly and chest. This is to be expected, as we learn from the larger and more representative CAESAR dataset. Both shape spaces trained from posture normalized models can better fit the arms compared to nonnormalized models.
7 Conclusion
In this work we address the challenging problem of building an efficient and expressive 3D body shape space from the largest commercially available 3D body scan dataset RobinetteTCP99 (). We carefully design and evaluate different data preprocessing steps required to obtain highquality body shape models. To that end we evaluate different template fitting procedures. We observe that shape and posture fitting of an initial shape space to a scan prior to nonrigid deformation considerably improves the fitting results. Our findings indicate that multiple passes over the dataset improve initialization and thus increase the overall fitting accuracy and statistical shape space qualities. Furthermore, we show that posture normalization prior to learning a shape space leads to significantly better generalization and specificity of the SSCAPE spaces. Finally, we demonstrate the advantages of our learned shape spaces over the stateoftheart shape space of Jain et al. Jain:2010:MovieReshape () learned on largest publicly available dataset HaslerBM2009 () on the task of human body shape reconstruction from noisy depth data.
We release our SSCAPE spaces, registered CAESAR scans, raw scan preprocessing code, code to fit a SSCAPE space to a raw scan and evaluation code for public usage^{1}^{1}1Available at http://humanshape.mpiinf.mpg.de. We believe this contribution is required for future development in human body modeling.
Acknowledgements
We thank Alexandros Neophytou and Adrian Hilton for sharing their posture normalization code. We also thank Mónica Vidriales and Gautham Adithya for their contributions to model fitting and evaluation code. This work was partially funded by the Cluster of Excellence MMCI.
References
 [1] B. Allen, B. Curless, and Z. Popović. The space of human body shapes: reconstruction and parameterization from range scans. TG’03.
 [2] B. Allen, B. Curless, Z. Popović, and A. Hertzmann. Learning a correlated model of identity and posedependent body shape variation for realtime synthesis. In SCA’06.
 [3] D. Anguelov, P. Srinivasan, D. Koller, S. Thrun, J. Rodgers, and J. Davis. SCAPE: shape completion and animation of people. TG’05.
 [4] S.Y. Baek and K. Lee. Parametric human body shape modeling framework for humancentered product design. CAD’12.
 [5] A. Balan and M. Black. The naked truth: Estimating body shape under clothing. In ECCV’08.
 [6] A. Balan, L. Sigal, M. Black, J. Davis, and H. Haussecker. Detailed human shape and pose from images. In CVPR’07.
 [7] F. Bogo, J. Romero, M. Loper, and M. Black. FAUST: Dataset and evaluation for 3D mesh registration. In CVPR’14.
 [8] J. Boisvert, C. Shu, S. Wuhrer, and P. Xi. Threedimensional human shape inference from silhouettes: Reconstruction and validation. MVAP’13.
 [9] Y. Chen, Z. Liu, and Z. Zhang. Tensorbased human body modeling. In CVPR’13.
 [10] C.H. Chu, Y.T. Tsai, C. Wang, and T.H. Kwok. Exemplarbased statistical model for semantic parametric design of human body. Comp. in Ind.’10.
 [11] C. Goodall. Procrustes Methods in the Statistical Analysis of Shape. J. R. Stat. Soc. Ser. B Stat. Methodol.’91.
 [12] P. Guan, L. Reiss, D. Hirshberg, A. Weiss, and M. Black. DRAPE: DRessing Any PErson. TG’12.
 [13] P. Guan, A. Weiss, A. Bălan, and M. Black. Estimating human shape and pose from a single image. In ICCV’09.
 [14] N. Hasler, H. Ackermann, B. Rosenhahn, T. Thormählen, and H.P. Seidel. Multilinear pose and body shape estimation of dressed subjects from image sets. In CVPR’10.
 [15] N. Hasler, C. Stoll, B. Rosenhahn, T. Thormählen, and H.P. Seidel. Estimating body shape of dressed humans. Comput. & Graph.’09.
 [16] N. Hasler, C. Stoll, M. Sunkel, B. Rosenhahn, and H.P. Seidel. A statistical model of human pose and body shape. CGF’09.
 [17] T. Helten, A. Baak, G. Bharai, M. Müller, H.P. Seidel, and C. Theobalt. Personalization and evaluation of a realtime depthbased full body scanner. In 3DV’13.
 [18] A. Jain, T. Thormählen, H.P. Seidel, and C. Theobalt. MovieReshape: tracking and reshaping of humans in videos. TG’10.
 [19] Matthew Loper, Naureen Mahmood, Javier Romero, Gerard PonsMoll, and Michael Black. Smpl: A skinned multiperson linear model. TG’15.
 [20] L. Mündermann, S. Corazza, and T. Andriacchi. Accurately measuring human movement using articulated icp with softjoint constraints and a repository of articulated models. In CVPR’07.
 [21] A. Neophytou and A. Hilton. Shape and pose space deformation for subject specific animation. In 3DV’13.
 [22] T. Neumann, K. Varanasi, N. Hasler, M. Wacker, M. Magnor, and C. Theobalt. Capture and statistical modeling of armmuscle deformations. CGF’13.
 [23] L. Pishchulin, A. Jain, M. Andriluka, T. Thormählen, and B. Schiele. Articulated people detection and pose estimation: Reshaping the future. In CVPR’12.
 [24] L. Pishchulin, A. Jain, C. Wojek, T. Thormaehlen, and B. Schiele. In good shape: Robust people detection based on appearance and shape. In BMVC’11.
 [25] Gerard PonsMoll, Javier Romero, Naureen Mahmood, and Michael Black. DYNA: a model of dynamic human shape in motion. TG’15.
 [26] K. Robinette, H. Daanen, and E. Paquet. The CAESAR project: A 3D surface anthropometry survey. In 3DIM’99.
 [27] C. Rupprecht, O. Pauly, C. Theobalt, and S. Ilic. 3d semantic parameterization for human shape modeling: Application to 3d animation. In 3DV’13.
 [28] H. Seo and N. MagnenatThalmann. An automatic modeling of human bodies from sizing parameters. In I3D ’03.
 [29] H. Seo, Y. In Yeo, and K. Wohn. 3D body reconstruction from photos based on range scan. In Edutainment’06.
 [30] L. Sigal, A. Balan, and M. Black. Combined discriminative and generative articulated pose and nonrigid shape estimation. In NIPS’07.
 [31] M. Styner, K. Rajamani, L.P. Nolte, G. Zsemlye, G. Székely, C. Taylor, and R. Davies. Evaluation of 3D correspondence methods for model building. In IPMI’03.
 [32] G. Tam, Z.Q. Cheng, Y.K. Lai, F. Langbein, Y. Liu, D. Marshall, R. Martin, X.F. Sun, and P. Rosin. Registration of 3D point clouds and meshes: A survey from rigid to nonrigid. TVCG’13.
 [33] Aggeliki Tsoli, Naureen Mahmood, and Michael Black. Breathing life into shape: Capturing, modeling and animating 3d human breathing. TG’14.
 [34] O. van Kaick, H. Zhang, G. Hamarneh, and D. CohenOr. A survey on shape correspondence. CGF’11.
 [35] A. Weiss, D. Hirshberg, and M. Black. Home 3D body scans from noisy image and range data. In ICCV’11.
 [36] S. Wuhrer and C. Shu. Estimating 3d human shapes from measurements. MVAP’13.
 [37] S. Wuhrer, C. Shu, and P. Xi. Landmarkfree posture invariant human shape correspondence. The Vis. Comput.’11.
 [38] S. Wuhrer, C. Shu, and P. Xi. Postureinvariant statistical shape analysis using Laplace operator. Comput. & Graph.’12.
 [39] Stefanie Wuhrer, Leonid Pishchulin, Alan Brunton, Chang Shu, and Jochen Lang. Estimation of human body shape and posture under clothing. CVIU’14.
 [40] P. Xi, W.S. Lee, and C. Shu. A datadriven approach to humanbody cloning using a segmented body database. In PG’07.
 [41] S. Zhou, H. Fu, L. Liu, D. CohenOr, and X. Han. Parametric reshaping of human bodies in images. TG’10.
 [42] C. Zhu, R. Byrd, P. Lu, and J. Nocedal. Algorithm 778: LBFGSB Fortran subroutines for largescale boundconstrained optimization. TOMS’97.