3D Human Pose Estimation from Deep Multi-View 2D Pose

2019·Arxiv

Abstract

Abstract

Human pose estimation - the process of recognizing a human’s limb positions and orientations in a video - has many important applications including surveillance, diagnosis of movement disorders, and computer animation. While deep learning has lead to great advances in 2D and 3D pose estimation from single video sources, the problem of estimating 3D human pose from multiple video sensors with overlapping fields of view has received less attention. When the application allows use of multiple cameras, 3D human pose estimates may be greatly improved through fusion of multi-view pose estimates and observation of limbs that are fully or partially occluded in some views. Past approaches to multi-view 3D pose estimation have used probabilistic graphical models to reason over constraints, including per-image pose estimates, temporal smoothness, and limb length. In this paper, we present a pipeline for multi-view 3D pose estimation of multiple individuals which combines a state-of-art 2D pose detector with a factor graph of 3D limb constraints optimized with belief propagation. We evaluate our results on the TUM-Campus and Shelf datasets for multi-person 3D pose estimation and show that our system significantly out-performs the previous state-of-the-art with a simpler model of limb dependency.

I. INTRODUCTION

Articulated human pose estimation is a long-studied problem in the field of computer vision [1], [2] which involves estimating a parameterized 2D or 3D human body model from video or still imagery. Pose estimation has many direct applications including motion capture for film and game production and may also be done as an pre-precessing step for higher-level vision problems such as action recognition. Most of the recent progress made in this area has leveraged advances in deep neural networks, particularly convolutional networks [3], to achieve impressive results in the areas of 2D and 3D pose estimation from single images [3], [4].

While significant breakthroughs, the utility of these approaches are limited by the single-view nature of the source data. 2D pose estimation techniques only detect joint and limb locations in image space, and cannot reveal anything directly about the positions and orientations of limbs in 3D. Approaches exist which perform 3D pose estimation on a single actor in a single image [5], but the range of 3D poses is heavily limited by the training dataset and limbs which are occluded must be inferred rather than directly measured. For applications which demand precise measurement of all 3D limbs, a synchronized network of overlapping calibrated cameras can be used to resolve single-view ambiguities and improve 3D triangulation accuracy.

To accelerate research in multi-person, multi-camera 3D pose estimation, Amin, Andriluka, Rohrbach, Schiele introduced the TUM Campus and Shelf datasets [6], [7], each of which contain video footage of several actors taken from 2-3 calibrated cameras. This dataset presents an opportunity to face challenges not present in single image pose estimation. For example, actors in the scene are frequently fully occluded by other actors in one view, and in this situation it becomes necessary to keep track of the actors through the occlusion event.

To address these challenges [6], and later [7], use the sum-product belief propagation algorithm to optimize over the space of 3D pose configurations, which they discretize by triangulating the results of a pose detector. The factor graph they use features a large number of parameters for factors such as joint collision and temporal smoothing. For their system to achieve optimal results, however, they need to use a structured SVM solver to optimize over a series of hyper parameters [7].

The system we propose, by contrast, achieves significantly better results with fewer constraints on limb motion. First, we use a modern 2D pose detector [4] which leverages advances in deep learning to generate better 2D limb hypotheses. Next we use a factor graph optimization, where factors are all constructed to be readily interpretable, and as such it is very easy to identify reasonable values for each factor without the need to resort to a more complex hyperparameter optimization scheme. Our final factor graph uses only three factors: limb location priors, collision terms, and temporal smoothing.

In this paper, we present our approach: a full pipeline for multiple person, multiple camera 3D pose estimation. After presenting our technique, we empirically show the effectiveness of our technique on the Campus and Shelf datasets, where we in some cases significantly out-perform the previous state-of-art.

II. RELATED WORK

A. 2D Human Pose Estimation

Approaches to human pose estimation can be categorized broadly into two separate tasks based on goals: single person pose estimation and multiple person pose estimation. In single person pose detection, a tight bounding box is provided around the person in question. Recently, stacked hourglass networks [3] and their variants have emerged as the leading deep learning methods for single person pose estimation. A stacked hourglass network is a fully-convolutional CNN architecture with multiple ”hourglass modules” which squeeze the image representation to a very small size, allowing for each unit in the CNN to have a very large receptive field.

Alternatively, multi-person pose estimation seeks a single architecture which, given a single image with multiple people present, can identify the individual people before or during the actual pose estimation process. These techniques can be ”topdown”, where people are first identified and then their poses are detected on an individual basis [8], or ”bottom-up”, where joints or limbs are detected first, and then aggregated into human detections only afterwards. To this end, Cao, Simon Wei, and Sheikh present Part Affinity Fields [4], which, in addition to detecting individual joints, also detects a vector indicating the orientation of each joint. These joints are then aggregated into separate skeletons, allowing for simultaneous person and part detection.

B. 3D Human Pose Estimation

Exploiting the success of 2D human pose estimation techniques, there has been a lot of work done in the area of 3D pose estimation, as inferred from a single image [9]. For example, [5] uses a modified stacked hourglass to progressively produce 3D pose detections with more and more depth.

3D pose estimation results can be improved when multiple cameras are available. [10] use a graphical model for 3D pose estimation and perform inference with a particle message passing scheme. [11] use a Sum of Gaussian based appearance model, which is later expanded in [12] to make use of deep learned features.

There also exist multiple datasets for 3D human pose estimation in controlled indoor environments [13]–[15], where multiple cameras observe a single actor in a studio. Many alternate datasets address the problem of 3D human pose detection with multiple humans, or in outdoor environments [6]. In particular, [6] introduces the TUM Campus and Shelf datasets, each of which feature three actors in a single video taken from at least three synchronized cameras. Following up on this, [16] and [7] use factor graphs optimized with belief propagation to perform inference on poses detected with 2D pose estimators.

III. APPROACH

The human body is represented as a graph with N = 14 nodes, where each node represents a joint in the human body (e.g. left hand, right shoulder, top of head, etc.), and edges represent adjacent joints. Our objective is to estimate, for each individual and timestep t = 0...T, the 3D joint locations where each . The input data comes in the form of video frames from each camera c, which have been pre-calibrated such that projection matrices have been computed and video frames have been aligned temporally.

Our processing pipeline for estimating the poses for each person from the video has multiple stages. First, 2D detections of joints are computed for each individual in the image using a state-of-art detector. When multiple people are present in the scene, these detections are used to identify specific individuals by triangulating 2D joint detections from different views. For each individual, a conditional random field (CRF) is constructed over 3D limb location variables to model the dependency on information sources including 2D detection strength, temporal smoothness, and collision terms. Inference on each factor graph is performed using the sum-product belief propagation algorithm to obtain posterior probabilities on limb location, from which the maximum likelihood 3D location is selected. These steps are described in detail in the following subsections.

A. 2D Pose Detection and Cross-frame Association

The first step in our pipeline is to obtain 2D limb detections in each frame of video from each camera. Many published approaches to this problem have made their software available and we choose to use the open source OpenPose library [4], which uses part-affinity fields to simultaneously extract human poses and detect individuals. Along with offering state-of-art 2D detection performance, this software can expose raw limb likelihood heat maps (see Figure 1a) used to compute the final 2D limb position. These heat maps allow propagation of uncertainty in the 2D detection step to the 3D optimization step, where a final decision can be made with all multi-view evidence. Each of these heat maps assigns a value in [0, 1] to each pixel in indicating the likelihood that the given pixel is displaying the projection of joint i.

Given a collection of 2D poses within each frame, we must match skeletons belonging to the same person both across time and multiple views. Consider an individual detected in 2D at time t for camera c. Here, each represent the pixel locations of a 2D joint, i.e. a local maximum of the corresponding . We determine that two detected individuals and from two consecutive frames are actually the same person if their bounding boxes have an intersection over union (IoU) score of at least 0.7. If a detection does not have an IoU match with any bounding box in the previous frame, then we determine it to be the first occurrence of that individual.

To determine matches between different cameras at the same point in time, i.e. between two poses and , we cast a ray into the scene for each pair of joints and in the pose, and measure the distance between the closest points of each ray (e.g., the distance between the wrist detection in and the same detection in ). If the average distance between these points is less than 2cm for all joints that are visible in both skeletons, we determine that the two skeletons belong to the same person.

When we have completed the matchings, we remove from each frame any identities that only appear in a single camera,

Fig. 1: (a) OpenPose outputs a heat map for each image, where intensity corresponds to the likelihood that a particular joint exists in that location. (b) In order to generate a set of possible states in for some joint , we convert into a probability mass function and sample from it.

leaving behind only detections that can be fed directly into our multi-view system.

B. Conditional Random Field for 3D Pose

At this stage, the joint locations from matched detections in multiple views can be back-projected into 3D and triangulated to obtain a 3D pose estimate. The result is noisy, lacks temporal smoothness frame-to-frame, and may occasionally contain gross errors for certain limbs if 2D detections are ambiguous. To improve the quality of the pose estimates, a CRF on 3D joint positions is constructed to allow the introduction of additional consistency constraints. The CRF is represented by a factor graph, where the joint positions are the variables and constraints on these joints become factors:

where Z is the partition function. The factors f are not necessarily true probability distributions although we do individually normalize all of our factors into probabilities, as we found that doing so leads to improved numerical stability.

The first of these factors, , represents the likelihood of a joint given evidence aggregated from the detections in the source images . The second factor is , which encourages the movement of joints between time steps to be temporally smooth. constrains certain symmetric pairs of joints not to collide. The set of edges represents the joints that cannot collide, which in our implementation represent symmetric body parts, e.g. right and left wrist, right and left

Fig. 2: A factor graph diagram for a simplified network involving two symmetric joints and and three timesteps. The joint variables are denoted with circles and the data, temporal, and collision factors are denoted with green, purple, and gray rectangles, respectively.

ankle, etc. An illustration of the dependencies of these factors on variables in Equation 1 is shown in Figure 2 for a subset of the full body network. Additional constraints, such as constant limb lengths, may be easily appended to the factorization if desired.

To make inference on the network tractable, the joint positions are limited to take on a discrete number of possible states. To discretize the state space for Y, we use a method similar to that of [7]. Let be the discrete set of M possible states that can take on. Our goal is to a select a small set of that well sample the 2D pose heat maps when projected into each image.

We iteratively select a subset of cameras of size with (a given joint i will generally be visible from either 2 or 3 cameras in the data we have). We then randomly sample from the 2D heat map in each camera and back- project the pixel coordinate into a ray. A 3D joint hypothesis is formed by computing the point in closest to all rays. For our implementation, we draw 16 samples this way when C = 3 and 64 samples when C = 2. In both cases, the final result is M = 64 samples, because when C = 3, there are three subsets of cameras with and one with . See Figure1b.

In the subsections that follow, we provide detailed descriptions of the potential functions in (1) and our CRF optimization procedure using the sum-product belief propagation algorithm [17] to obtain posterior estimates on the joint positions.

1) Data Term: The data term in the factor graph measures 2D pose confidence of the projected 3D joint hypothesis. This term is straight-forward to calculate: for each possible state , we calculate by projecting from world coordinates onto each and averaging the results:

where is the projection of s in the image. Referring back to the 2D pose heat maps, as opposed to the single optimal 2D locations, allows this data term to fuse weak evidence in the case of ambiguity, such as bimodal distributions around left/right joints in close proximity to each other.

2) Temporal Term: The purpose of the temporal term is to bias joints to move smoothly frame-to-frame. Without some form of temporal smoothing, we found that animations of the estimated 3D body models contained significant jitter as a result of each frame being calculated separately. Having information flow between time steps also helps the system to recover from short-term confusion events such as occlusion of limbs by the torso during walking.

We model temporal smoothing with a ternary term, involving 3D joint estimates from three consecutive time steps: , and . Specifically:

where is the expected position of the joint using a constant velocity assumption and samples immediately before and immediately after the current time step. A state hypothesis is penalized if it veers too far from the constant-velocity prediction. We experimentally set cm, though we found that the results were not very sensitive to the exact choice of so long as the choice was reasonably small.

3) Collision Term: The collision term prevents symmetric left/right joints from colliding with each other in 3D. This common issue in pose detection is generally caused by the fact that symmetric joints, like the right and left feet, look very similar to each other and are often in close proximity. This tends to result in a bimodal distribution in the 2D pose heat maps, which in turn can cause both symmetric joints in a pair to be assigned to the same location in 3D space.

In contrast to [7], who handle collision with a trained Gaussian on the distance between limbs, we instead use a modified sigmoid function that evaluates to 1 when the two points are reasonably far apart, and 0 when they are too close:

where is a hyperparameter that offsets the center of the sigmoid, and is a hyperparameter that tightens the sigmoid. We empirically set cm and cm. As with , we find that as long as the values chosen are reasonable, they don’t have a huge effect on the results.

C. Belief Propagation Optimization

The factor graph defines the dependencies between joint variables and input video frames across time and space in the total probability distribution. With the graph constructed, the next objective is to perform inference to obtain posterior estimates of the 3D joint positions . Belief propagation is a message-passing algorithm which estimates the posterior distribution in CRFs and has been used in many computer vision applications [18]. In the case where the graphs have loops, such as ours, belief propagation is not guaranteed to converge to the true posterior but may be used as an approximation regardless.

As belief propagation will iteratively optimize all frames of a video, an intelligent schedule for organizing the message passing between factors and variables is needed. Our implementation updates all factor-to-variable messages of the same type before moving on to the next type, i.e. all messages are updated, then all messages, etc. The temporal factor update passes messages forward in time through the whole video clip, then back again, so that temporal information is rapidly propagated to distant frames. Future work could easily adapt this schedule to process live video by temporally optimizing over a sliding window of frames rather than the entire video clip. We run all of our tests with five iterations of belief propagation across the full network. We experimentally found that running more than five iterations did not change the values because by the fifth iteration the values had effectively converged.

IV. EXPERIMENTS

We evaluate our 3D pose estimation pipeline on the TUM Campus and TUM Shelf [6] datasets, which have labeled 3D ground truth for multiple people moving in the scene. The inclusion of multiple people adds an extra dimension of challenge not present in many other 3D pose estimation datasets [14], [15], because as the people move and interact it is often the case that they obscure each other.

The TUM Campus and Shelf datasets both feature scenes shot from 3 calibrated cameras containing 2-3 people at any given time. We evaluated the Campus dataset on 220 frames of video with resolution and the Shelf dataset on 300 frames with resolution . During these test portions of the sequence, actors frequently come in and out of view, occlude each other, and stand stationary for long periods of time.

To measure the performance of our algorithm, we use the PCP (percentage of correctly estimated parts) metric, following the example set by [6]. Using the PCP metric, originally described by [19], two adjacent joints are considered to be correctly estimated if their respective distances from their ground truth locations are less than of the length of the limb in the ground truth.

A. Comparison to Other Methods

In this section, we compare the performance of our system to other state-of-the-art methods on the Shelf and Campus

TABLE I: Performance breakdown by part for on the Campus and Shelf datasets. Best results are in bold. The “Average” column is weighted by how often each actor appears in the dataset.

Fig. 3: A superimposed comparison of our method’s output (blue) with the ground truth pose (red) on a single frame of the Campus dataset. Note that the head and neck joints of our output are systematically lower than on the ground truth, negatively affecting our accuracy.

datasets, as broken down by individual joints. The results of these experiments are shown in Table I.

With the exception of Actor 1 in the Campus dataset, our method out-performs all other methods on the six actors of the Campus and Shelf datasets. Within the Campus dataset, we find that the method from [7] performs slightly better on each actor’s arms, but that our method makes up the difference by performing better on the legs and torso. We attribute this to differences in the output of the 2D part detectors.

On the other hand, our method performs almost unilaterally better on the Shelf dataset. This discrepancy is likely due to the resolution difference between the two datasets; at a native resolution of , images from the Campus dataset must be scaled up to be used by the OpenPose 2D pose detector, which processes input at a resolution of . The Shelf dataset’s images, on the other hand, can make use of the full discriminative power of OpenPose’s part detector.

During computation of these metrics, the 3D head pose output by our system was manually adjusted by a fixed amount to account for a discrepancy between the OpenPose

Fig. 4: The benefits of temporal smoothing. (a) Using only the data term, it’s possible for certain joints to be incorrectly placed on isolated frames. (b) The temporal term can propagate information from nearby frames to prevent these mistakes.

and ground truth body models. OpenPose is trained on the Microsoft COCO 2016 Dataset [20] to detect the head joint around the level of the nose, whereas the ground truth labels the head joint as the top of head. Similarly, OpenPose detects the neck as lying directly between the shoulder blades, as opposed to just above them. An example is shown in Figure 3. To perform a fair comparison with ground truth, the head and neck detections output by our system were translated upward in the z direction by 10cm for evaluation. This approximates the correct offset in these datasets since the range of actor pose is limited to upright poses. In the future we may extend the systems to learn the correct offsets, but longer term we expect that a community standardized human body model will emerge to eliminate these types of discrepancies in future datasets.

We also note that while Belagiannis et al. do out-perform our method on one actor out of six, our method always achieves accuracy of at least 82%, whereas the other methods we compare against score below 80% on at least half of the actors, making our method more consistent across sequences.

TABLE II: The effects of each term on the results for the Campus and Shelf datasets.

B. Analysis of Individual Factors

Table I breaks down the effect of each term in the factor graph on our final result accuracy. When the data term alone is used, belief propagation is not run and the 3D pose hypothesis with the highest average 2D pose likelihood is selected for the answer. As such, the data-only numbers reflect the performance of a trivial triangulation of the 2D pose output without any additional reasoning. We note that performance with the data term is competitive with other state-of-the-art methods, and in the case of the Shelf dataset actually out-performs all other methods. This strong performance can be attributed to the highly effective 2D part detection of OpenPose.

The main contribution of the temporal term on performance is its ability to correct errors that exist only in an isolated frame. Figure 4 illustrates this; errors are corrected by information flow from adjacent frames which do not exhibit the same flaws. Additionally, when viewing animated videos of the body-models estimated with and without temporal terms side-by-side, we observe that the temporal result is significantly less jittery. This smoothness is critical to applications such as motion capture, but is not well-captured by the PCP metric.

The collision term has a very minor effect in these results and corrects a small number of minor errors caused by joints being improperly disambiguated. Its effect may be more pronounced in datasets with more complex human motions not captured in the Campus and Shelf datasets. However, we conclude that for simple applications it could likely be removed for a modest improvement in runtime at a negligible accuracy penalty.

V. CONCLUSION

We advance the previous state-of-art in multi-person multi-camera pose estimation by combining a modern 2D pose detector based on deep networks with a factor graph optimization to reason over multiple views. We place fewer constraints on limb interactions than previous techniques, yet achieve higher metrics in the Campus and Shelf evaluation datasets. Future work will include adapting the factor graph to allow processing of live video by iterating over short sliding temporal windows.

ACKNOWLEDGEMENTS

The research was supported in part by the Office of Naval Research (N00014-16-P-2040).

REFERENCES

[1] P. F. Felzenszwalb and D. P. Huttenlocher, “Pictorial structures for object recognition,” Int. J. Comput. Vision, vol. 61, no. 1, pp. 55–79, Jan. 2005.

[2] M. Andriluka, S. Roth, and B. Schiele, “Pictorial structures revisited: People detection and articulated pose estimation,” in 2009 IEEE Conference on Computer Vision and Pattern Recognition, June 2009, pp. 1014–1021.

[3] A. Newell, K. Yang, and J. Deng, “Stacked Hourglass Networks for Human Pose Estimation,” mar 2016.

[4] Z. Cao, T. Simon, S.-E. Wei, and Y. Sheikh, “Realtime Multi-Person 2D Pose Estimation using Part Affinity Fields,” nov 2016.

[5] G. Pavlakos, X. Zhou, K. G. Derpanis, and K. Daniilidis, “Coarse-to- fine volumetric prediction for single-image 3d human pose,” CoRR, vol. abs/1611.07828, 2016.

[6] S. Amin, M. Andriluka, M. Rohrbach, and B. Schiele, “Multi-view Pictorial Structures for 3D Human Pose Estimation,” in Procedings of the British Machine Vision Conference 2013, 2013.

[7] V. Belagiannis, S. Amin, M. Andriluka, B. Schiele, N. Navab, and S. Ilic, “IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE 1 3D Pictorial Structures Revisited: Multiple Human Pose Estimation.”

[8] G. Papandreou, T. Zhu, N. Kanazawa, A. Toshev, J. Tompson, C. Bregler, and K. P. Murphy, “Towards accurate multi-person pose estimation in the wild,” CoRR, vol. abs/1701.01779, 2017.

[9] X. Zhou, M. Zhu, G. Pavlakos, S. Leonardos, K. G. Derpanis, and K. Daniilidis, “MonoCap: Monocular Human Motion Capture using a CNN Coupled with a Geometric Prior,” jan 2017.

[10] L. Sigal, M. Isard, H. Haussecker, and M. J. Black, “Loose-limbed people: Estimating 3d human pose and motion using non-parametric belief propagation,” Int. J. Comput. Vision, vol. 98, no. 1, pp. 15–48, May 2012.

[11] A. Elhayek, E. de Aguiar, A. Jain, J. Tompson, L. Pishchulin, M. An- driluka, C. Bregler, B. Schiele, and C. Theobalt, “Efficient convnet-based marker-less motion capture in general scenes with a low number of cameras,” in 2015 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), June 2015, pp. 3810–3818.

[12] A. Elhayek, E. de Aguiar, A. Jain, J. Thompson, L. Pishchulin, M. Andriluka, C. Bregler, B. Schiele, and C. Theobalt, “Marconi—convnet-based marker-less motion capture in outdoor and indoor scenes,” IEEE Trans. Pattern Anal. Mach. Intell., vol. 39, no. 3, pp. 501–514, Mar. 2017.

[13] L. Sigal and M. J. Black, “Measure locally, reason globally: Occlusion- sensitive articulated pose estimation,” in Proceedings of the 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition - Volume 2, ser. CVPR ’06. Washington, DC, USA: IEEE Computer Society, 2006, pp. 2041–2048.

[14] L. Sigal, A. O. Balan, and M. J. Black, “Humaneva: Synchronized video and motion capture dataset and baseline algorithm for evaluation of articulated humanmotion,” International Journal of Computer Vision, vol. 87, no. 1, p. 4, Aug 2009.

[15] C. Ionescu, D. Papava, V. Olaru, and C. Sminchisescu, “Human3.6m: Large scale datasets and predictive methods for 3d human sensing in natural environments,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 36, no. 7, pp. 1325–1339, jul 2014.

[16] V. Belagiannis, S. Amin, M. Andriluka, B. Schiele, N. Navab, and S. Ilic, “3d pictorial structures for multiple human pose estimation,” in 2014 IEEE Conference on Computer Vision and Pattern Recognition, June 2014, pp. 1669–1676.

[17] F. R. Kschischang, B. J. Frey, and H. A. Loeliger, “Factor graphs and the sum-product algorithm,” IEEE Transactions on Information Theory, vol. 47, no. 2, pp. 498–519, Feb 2001.

[18] A. Quattoni, M. Collins, and T. Darrell, “Conditional random fields for object recognition,” in Advances in Neural Information Processing Systems 17, L. K. Saul, Y. Weiss, and L. Bottou, Eds. MIT Press, 2005, pp. 1097–1104.

[19] V. Ferrari, M. Marin-Jimenez, and A. Zisserman, “Progressive search space reduction for human pose estimation,” pp. 1–8, June 2008.

[20] T. Lin, M. Maire, S. J. Belongie, L. D. Bourdev, R. B. Girshick, J. Hays, P. Perona, D. Ramanan, P. Doll´ar, and C. L. Zitnick, “Microsoft COCO: common objects in context,” CoRR, vol. abs/1405.0312, 2014.

designed for accessibility and to further open science