In this paper, we address the problem of autonomous cinematography using unmanned aerial vehicles (UAVs). Specifically, we focus on scenarios where an UAV must film an actor moving through an unknown environment at high speeds, in an unscripted manner. Filming dynamic actors among clutter is extremely challenging, even for experienced pilots. It takes high attention and effort to simultaneously predict how the scene is going to evolve, control the UAV, avoid obstacles and reach desired viewpoints. Towards solving this problem, we present a complete system that can autonomously handle the real-life constraints involved in aerial cinematography: tracking the actor, mapping out the surrounding terrain and planning maneuvers to capture high quality, artistic shots.
Consider the typical filming scenario in Fig 1. The UAV must accomplish a number of tasks. First, it must estimate the actor’s pose using an onboard camera and forecast their future motion. The pose estimation should be robust to changing viewpoints, backgrounds and lighting conditions. Accurate forecasting is key for anticipating events which require changing camera viewpoints. Secondly, the UAV must remain safe as it flies through through new environments. Safety requires explicit modelling of environmental uncertainty. Finally, the UAV must capture high quality videos which require maximizing a set of artistic guidelines. The key challenge is that all these tasks must be done in real-time under limited onboard computational resources.
There is a rich history of work in autonomous aerial filming that tackles these challenges in part. For instance, several works focus on artistic guidelines [1, 2, 3, 4] but often rely on perfect actor localization through a high-precision RTK GNSS or a motion-capture system. Additionally, while the majority of work in the area deals with collisions between UAV and actors[2, 1, 5], the environment is not factored in. While there are several successful commercial products, they too have certain limitations to either low speed and low clutter regimes (e.g. DJI Mavic ) or shorter planning horizons (e.g. Skydio R1 ). Even our previous work , despite handling environmental occlusions and collisions, assumes a prior elevation map and uses GPS to localize the actor. Such simplifications impose restrictions on the diversity of scenarios that the system can handle.
We address these challenges by building upon previous work that formulates the problem as an efficient real-time trajectory optimization . We make a key observation: we don’t need prior ground-truth information about the scene; our onboard sensors suffice to attain good performance. However, sensor data is noisy and needs to be processed in real-time; we develop robust and efficient algorithms to do so. To tackle actor localization, we use a visual tracking system. To map the environment, we use a long-range LiDAR and process it incrementally to build a signed distance field of the environment. Using these two, we plan over long horizons in unknown environments to film fast dynamic actors according to artistic guidelines.
In summary, our main contributions in this paper are threefold:
We develop an incremental signed distance transform algorithm for large-scale real-time environment mapping (Section IV-B);
We develop a complete system for autonomous cinematography that includes visual actor localization, online mapping, and efficient trajectory optimization that can deal with noisy measurements (Section IV);
We offer extensive quantitative and qualitative performance evaluations of our system both in simulation and field tests, while also comparing performance changes with scenarios with full map and actor knowledge (Section V).
Ii Problem Formulation
The overall task is to control a UAV to film an actor who is moving through an unknown environment. We formulate this as a trajectory optimization problem where the cost function measures shot quality, environmental occlusion of the actor, jerkiness of motion and safety. This cost function depends on the environment and the actor, both of which must be sensed on-the-fly. The changing nature of environment and actor trajectory also demands re-planning at a high frequency.
Let be the trajectory of the UAV, i.e., . Let be the trajectory of the actor, . The state of the actor, as sensed by onboard cameras, is fed into a prediction module that computes (Section IV-A).
be a voxel occupancy grid that maps every point in space to a probability of occupancy. Letbe the signed distance values of a point to the nearest obstacle. Positive sign is for points in free space, and negative sign is for points either in occupied or unknown space, which we assume to be potentially inside an obstacle. The UAV senses the environment with the onboard LiDAR, updates grid , and then updates (Section IV-B).
We briefly touch upon the four components of the cost function (refer to Section IV-C for mathematical expressions).
Smoothness : Penalizes jerky motions that may lead to camera blur and unstable flight;
Shot quality : Penalizes poor viewpoint angles and scales that deviate from the artistic guidelines
Safety : Penalizes proximity to obstacles that are unsafe for the UAV.
Occlusion : Penalizes occlusion of the actor by obstacles in the environment.
The objective is to minimize subject to initial boundary constraints .
The solution is then tracked by the UAV.
Iii Related Work
Camera control in virtual cinematography has been extensively examined by the computer graphics community as reviewed by . These methods typically reason about the utility of a viewpoint in isolation, following artistic principles and composition rules [10, 11] and employ either optimization-based approaches to find good viewpoints or reactive approaches to track the virtual actor. The focus is typically on through-the-lens control where a virtual camera is manipulated while maintaining focus on certain image features [12, 13, 14, 15]. However, virtual cinematography is free of several real-world limitations such as robot physics constraints and assumes full knowledge of the map.
Autonomous aerial cinematography
Several contributions on aerial cinematography focus on keyframe navigation. [16, 17, 18, 19, 20] provide user interface tools to re-time and connect static aerial viewpoints to provide smooth and dynamically feasible trajectories, as well as a visually pleasing images.  use key-frames defined on the image itself instead of world coordinates.
Other works focus on tracking dynamic targets, and employ a diverse set of techniques for actor localization and navigation. For example, [5, 22] detect the skeleton of targets from visual input, while others approaches rely on off-board actor localization methods from either motion-capture systems or GPS sensors [1, 3, 2, 4, 8]. These approaches have a varying level of complexity: [8, 4] can avoid obstacles and occlusions with the environment and with actors, while other approaches only handle collisions and occlusions caused by actors. We note differences in onboard versus off-board computing, and we observe distinct trajectory generation methods randing from trajectory optimization to search-based planners. Different contributions are summarized in Table I. It is important to notice that none of the previous approaches provide a solution for online environment mapping.
Online environment mapping
Dealing with imperfect representations of the world becomes a bottleneck for viewpoint optimization in physical environments. As the world is sensed online, it is usually incrementally mapped using voxel occupancy maps . To evaluate a viewpoint, methods typically raycast on such maps, which can be very expensive [24, 25]. Recent advances in mapping have led to better representations that can incrementally compute the truncated signed distance field (TSDF) [26, 27], i.e. return the distance and gradient to nearest object surface for a query. TSDFs are a suitable abstraction layer for planning approaches and have already been used to efficiently compute collision-free trajectories for UAVs [28, 29].
Visual target state estimation
Accurate object state estimation with monocular cameras is critical to many robot applications, including autonomous aerial filming. Deep networks have shown success in detecting objects [30, 31] and estimating 3D heading [32, 33] with several efficient architectures developed specifically for mobile applications [34, 35]
. However, many models do not generalize well to other tasks (e.g., aerial filming) due to data mismatch in terms of angles and scales. Our recent work in semi-supervised learning shows promise in increasing model generalizability with little labeled data by leveraging temporal continuity in training videos.
Our work exploits synergies at the confluence of several domains of research to develop an aerial cinematography platform that can follow dynamic targets in unknown and unstructured environments, as detailed next in our approach.
We now detail our approach for each sub-system of the aerial cinematography platform. At a high-level, three main sub-systems operate together: (A) Vision, required for localizing the target’s position and orientation and for recording the final UAV’s footage; (B) Mapping, required for creating an environment representation; and (C) Planning, which combines the actor’s pose and the environment to calculate the UAV trajectory. Fig. 2 shows a system diagram.
Iv-a Vision sub-system
We use only monocular images from the UAV’s gimbal and the UAV’s own state estimation to forecast the actor’s trajectory in the world frame. The vision sub-system counts with four main steps: actor bounding box detection and tracking, heading angle estimation; global ray-casting, and finally a filtering stage. Figure 3 summarizes the pipeline.
a) Detection and tracking: Our detection module is based on the MobileNet network architecture, due to its low memory usage and fast inference speed, which are well-suited for real-time applications on an onboard computer. We use the same network structure as detailed in our previous work . Our model is further trained with COCO  and fine-tuned on a custom aerial filming dataset. We limit the detection categories to person, car, bicycle, and motorcycle, which commonly appear in aerial filming. After a successful detection we use Kernelized Correlation Filters  to track the template over the next incoming frames. We actively position the independent camera gimbal with a PD controller to frame the actor on the desired screen position, following the commanded artistic principles (Fig. 5).
b) Heading Estimation: Accurate heading angle estimation is vital for the UAV to follow the correct shot type (front, back, left, right). As discussed in , human 2D pose estimation has been widely studied [39, 40], but 3D heading direction cannot be trivially recovered directly from 2D points because depth remains undefined. Therefore, we use the model architecture from  (Fig. 4), which takes as input a bounding box image, and outputs the cosine and sine values of the heading angle. This network uses a double loss during training, summing both errors in heading direction and temporal continuity. The latter loss is particularly useful to train the regressor on small datasets, following a semi-supervised approach.
c) Ray-casting: Given the actor’s bounding box and heading estimate on image space, we project the center-bottom of the bounding box onto the world’s ground plane and transform the actor’s heading using the camera’s state estimation to obtain the actor’s position and heading in world coordinates.
d) Motion Forecasting: The current actor pose updates a Kalman Filter (KF) to forecast the actor’s trajectory . We use separate KF models for people and vehicle dynamics.
Iv-B Mapping sub-system
As explained in Section II, the motion planner uses signed distance values in the optimization cost functions. The role of the mapping is to register LiDAR points from the onboard sensor, update the occupancy grid , and incrementally update the signed distance :
a) LiDAR registration: We receive approximately points per second from the laser mounted at the bottom of the aircraft. We register the points in the world coordinate system using a rigid body transform between sensor and the aircraft plus the UAV’s state estimation, which fuses GPS, barometer, internal IMUs and accelerometers. For each point we also store the corresponding sensor coordinate, used for the occupancy grid update.
b) Occupancy grid update: Points can represent either a hit (successful return) or a miss (return too far or too close from the sensor). We filter all expected misses caused by the aircraft’s structure, and then probabilistically update all voxels from between the sensor and LiDAR point. In practice, we use a grid size of , with square voxels that store an integer value between as the occupancy probability. Algorithm 1 covers the update process.
c) Incremental distance transform update: An important assumption we make in this work is that all unknown voxels are considered occupied (negative value) in the sign calculation of . Despite being pessimistic, this is a necessary assumption given the aerial cinematography occlusion cost function detailed in Section IV-C. Without this assumption, there is no way one can distinguish free regions from voxels which are inside of large obstacles and that would never be updated with sensor readings, making it impossible to evaluate the severity of occlusions.
As seen in Algorithm 1, voxel state changes are stored in and during the update. We use these changed cells as input to an incremental truncated signed distance transform (iTSDT) algorithm, modified from . This algorithm stores a separate data structure from , where all cells are considered initially free, and as voxel changes arrive, it incrementally updates using an efficient wavefront expansion technique.
The original iTSDT algorithm  only calculated the distance between each free cell and its closest occupied neighbor. Since our problem requires a signed version of DT, we introduced two major modifications. First, we only represent the borders of obstacles as occupied in the data structure. As seen in lines of Algorithm 1, border voxels are defined as all the unknown neighbors from a voxel that just became free, or when a cell becomes occupied (lines ). Second, when evaluating , we query the value of to attribute the sign, differentiating free (positive) from unknown or occupied (negative) cells.
Iv-C Planning sub-system
We want trajectories which are smooth, capture high quality viewpoints, avoid occlusion and are safe. Given the real-time nature of our problem, we desire fast convergence to locally optimal solutions rather than globally optimality taking a long time to obtain a solution. A popular approach is to cast the problem as an unconstrained optimization and apply covariant gradient descent . This is a quasi-Newton approach where some of the objectives have analytic Hessians that are easy to invert and are well-conditioned. Hence such methods exhibit fast convergence while being stable and computationally inexpensive.
For this implementation, we use a waypoint parameterization of trajectories, i.e., . The heading dimension is set to always point the drone from towards . We design a set of differentiable cost functions as follows:
We measure smoothness as the cumulative derivatives of the trajectory. Let be a discrete difference operator. The smoothness cost is:
where is a weight for different orders, and is the number of orders. We set . Note that smoothness is a quadratic objective where the Hessian is analytic.
Written in a quadratic form, shot quality measures the average squared distance between and an ideal trajectory that only considers positioning via cinematography parameters. can be computed analytically: for each point in the actor motion prediction, the ideal drone position lies on a sphere centered at with radius defined by the shot scale, relative yaw angle and relative tilt angle (Fig. 5):
We can then define the safety cost function (similar to )
Even though the concept of occlusion is binary, i.e, we either have or don’t have visibility of the actor, a major contribution of our past work  was defining a differentiable cost that expresses a viewpoint’s occlusion intensity among arbitrary obstacle shapes. Mathematically, we define occlusion as the integral of the TSDT cost over a 2D manifold connecting both trajectories and . The manifold is built by connecting each drone-actor position pair in time using the path .
Our objective is to minimize the total cost function (Eq. 1). We do so by covariant gradient descent, using the gradient of the cost function , and a analytic approximation of the Hessian :
This step is repeated till convergence. We follow conventional stopping criteria for descent algorithms, and limit the maximum number of iterations. Note that we only perform the matrix inversion once, outside of the main optimization loop, rendering good convergence rates . We use the current trajectory as initialization for the next planning problem.
V-a Experimental setup
Our platform is a DJI M210 drone, shown in Figure 6. All processing is done with a NVIDIA Jetson TX2, with 8GB of RAM and 6 CPU cores. An independently controlled gimbal DJI Zenmuse X4S records high-resolution images. Our laser sensor is a Velodyne Puck VLP-16 Lite, with vertical field of view and 100m max range.
V-B Field test results
Visual actor localization
We validate the precision of our pose and heading estimation modules in two experiments where the drone hovers and visually tracks the actor. First, the actor walks between two points over a straight line, and we compare the estimated and ground truth path lengths. Second, the actor walks on a circle at the center of a football field, and we verify the errors in estimated positioning and heading direction. Fig. 8 summarizes our findings.
Integrated field experiments
We test the real-time performance of our integrated system in several field test experiments. We use our algorithms in unknown and unstructured environments outdoors, following different types of shots and tracking different types of actors (people and bicycles) at both low and high speeds in unscripted scenes. Fig. 9 summarizes the most representative shots, and the supplementary video shows the final footages collected, along with visualizations of point clouds and the online map.
|Thread (%)||2.0cmRAM (MB)||2.0cmRuntime (ms)|
V-C Performance comparison with full information knowledge
An important hypothesis behind our system we can operate with a noisy estimate of the actor’s location and a partially known map with insignificant loss in performance. We compare our system with three assumptions present in previous work:
Ground-truth obstacles vs. online map
We compare average planning costs between results from a real-life test where the planner operated while mapping the environment in real time with planning results with the same actor trajectory but with full knowledge of the map beforehand. Results are averaged over 140 s of flight and approximately 700 planning problems. Table III shows a small increase in average planning costs, and Fig 10a shows that qualitatively both trajectories differ minimally. The planning time, however, doubles in the online mapping case due to extra load on CPU.
|time(ms)||1.5cmAvg. cost||1.5cmMedian cost|
|Noise in actor||30.2||0.1276||0.0953|
Ground-truth actor versus noisy estimate
We compare the performance between simulated flights where the planner has full knowledge of the actor’s position versus artificially noisy estimates with 1m amplitude. Results are also averaged over 140 s of flight and approximately 700 planning problems, and are displayed on Table III. The large cost difference is due to the shot quality cost, which relies on the actor’s position forecast and is severely penalized by the noise. However, if compared with the actor’s ground-truth trajectory, the difference in cost would be significantly smaller, as seen by the proximity of both final trajectories in Fig 10b. These results offer insight on the importance of our smoothness cost function when handling the noisy visual actor localization.
Height map assumption vs. 3D map
We present a complete system for autonomous aerial cinematography that can localize and track actors in unknown and unstructured environments with onboard computing in real time. Our platform uses a monocular visual input to localize the actor’s position, and a custom-trained network to estimate the actor’s heading direction. Additionally, it maps the world using a LiDAR and incrementally updates a signed distance map. Both of these are used by a camera trajectory planner that produces smooth and artistic trajectories while avoiding obstacles and occlusions. We evaluate the system extensively in different real-world tasks with multiple shot types and varying terrains.
We are actively working on a number of directions based on lessons learned from field trials. Our current approach assumes a static environment. Even though the our mapping can tolerate motion, a principled approach would track moving objects and forecast their motion. The TSDT is expensive to maintain because every-time unknown space is cleared, a large update is computed. We are looking into a just-in-time update that processes only the subset of the map queried by the planner, which is often quite small.
Currently, we do not close the loop between the image captured and the model used by the planner. Identifying model errors, such as actor forecasting or camera calibration, in an online fashion is a challenging next step. The system may also lose the actor due to tracking failures or sudden course changes. An exploration behavior to reacquire the actor is essential for robustness.
We thank Mirko Gschwindt, Xiangwei Wang, Greg Armstrong for the assistance in field experiments and robot construction.
-  N. Joubert, D. B. Goldman, F. Berthouzoz, M. Roberts, J. A. Landay, P. Hanrahan et al., “Towards a drone cinematographer: Guiding quadrotor cameras using visual composition principles,” arXiv preprint arXiv:1610.01691, 2016.
-  T. Nägeli, L. Meier, A. Domahidi, J. Alonso-Mora, and O. Hilliges, “Real-time planning for automated multi-view drone cinematography,” ACM Transactions on Graphics (TOG), vol. 36, no. 4, p. 132, 2017.
-  Q. Galvane, J. Fleureau, F.-L. Tariolle, and P. Guillotel, “Automated cinematography with unmanned aerial vehicles,” arXiv preprint arXiv:1712.04353, 2017.
-  Q. Galvane, C. Lino, M. Christie, J. Fleureau, F. Servant, F. Tariolle, P. Guillotel et al., “Directing cinematographic drones,” ACM Transactions on Graphics (TOG), vol. 37, no. 3, p. 34, 2018.
-  C. Huang, F. Gao, J. Pan, Z. Yang, W. Qiu, P. Chen, X. Yang, S. Shen, and K.-T. T. Cheng, “Act: An autonomous drone cinematography system for action scenes,” in 2018 IEEE International Conference on Robotics and Automation (ICRA). IEEE, 2018, pp. 7039–7046.
-  “Dji mavic,” https://www.dji.com/mavic, accessed: 2019-02-28.
-  Skydio. (2018) Skydio R1 self-flying camera. [Online]. Available: https://www.skydio.com/technology/
-  R. Bonatti, Y. Zhang, S. Choudhury, W. Wang, and S. Scherer, “Autonomous drone cinematographer: Using artistic principles to create smooth, safe, occlusion-free trajectories for aerial filming,” International Symposium on Experimental Robotics, 2018.
-  M. Christie, P. Olivier, and J.-M. Normand, “Camera control in computer graphics,” in Computer Graphics Forum, vol. 27, no. 8. Wiley Online Library, 2008, pp. 2197–2218.
-  D. Arijon, “Grammar of the film language,” 1976.
-  C. J. Bowen and R. Thompson, Grammar of the Shot. Taylor & Francis, 2013.
-  M. Gleicher and A. Witkin, “Through-the-lens camera control,” in ACM SIGGRAPH Computer Graphics, vol. 26, no. 2. ACM, 1992, pp. 331–340.
-  S. M. Drucker and D. Zeltzer, “Intelligent camera control in a virtual environment,” in Graphics Interface. Citeseer, 1994, pp. 190–190.
-  C. Lino, M. Christie, R. Ranon, and W. Bares, “The director’s lens: an intelligent assistant for virtual cinematography,” in Proceedings of the 19th ACM international conference on Multimedia. ACM, 2011, pp. 323–332.
-  C. Lino and M. Christie, “Intuitive and efficient camera control with the toric space,” ACM Transactions on Graphics (TOG), vol. 34, no. 4, p. 82, 2015.
-  M. Roberts and P. Hanrahan, “Generating dynamically feasible trajectories for quadrotor cameras,” ACM Transactions on Graphics (TOG), vol. 35, no. 4, p. 61, 2016.
-  N. Joubert, M. Roberts, A. Truong, F. Berthouzoz, and P. Hanrahan, “An interactive tool for designing quadrotor camera shots,” ACM Transactions on Graphics (TOG), vol. 34, no. 6, p. 238, 2015.
-  C. Gebhardt, S. Stevsic, and O. Hilliges, “Optimizing for aesthetically pleasing quadrotor camera motion,” 2018.
-  C. Gebhardt, B. Hepp, T. Nägeli, S. Stevšić, and O. Hilliges, “Airways: Optimization-based planning of quadrotor trajectories according to high-level user goals,” in Proceedings of the 2016 CHI Conference on Human Factors in Computing Systems. ACM, 2016, pp. 2508–2519.
-  K. Xie, H. Yang, S. Huang, D. Lischinski, M. Christie, K. Xu, M. Gong, D. Cohen-Or, and H. Huang, “Creating and chaining camera moves for quadrotor videography,” ACM Transactions on Graphics, vol. 37, p. 14, 2018.
-  Z. Lan, M. Shridhar, D. Hsu, and S. Zhao, “Xpose: Reinventing user interaction with flying cameras,” in Robotics: Science and Systems, 2017.
-  C. Huang, Z. Yang, Y. Kong, P. Chen, X. Yang, and K.-T. T. Cheng, “Through-the-lens drone filming,” in 2018 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). IEEE, 2018, pp. 4692–4699.
-  S. Thrun, W. Burgard, and D. Fox, Probabilistic robotics. MIT press, 2005.
-  S. Isler, R. Sabzevari, J. Delmerico, and D. Scaramuzza, “An information gain formulation for active volumetric 3d reconstruction,” in ICRA, 2016.
-  B. Charrow, G. Kahn, S. Patil, S. Liu, K. Goldberg, P. Abbeel, N. Michael, and V. Kumar, “Information-theoretic planning with trajectory optimization for dense 3d mapping,” in RSS, 2015.
-  R. A. Newcombe, S. Izadi, O. Hilliges, D. Molyneaux, D. Kim, A. J. Davison, P. Kohi, J. Shotton, S. Hodges, and A. Fitzgibbon, “Kinectfusion: Real-time dense surface mapping and tracking,” in Mixed and augmented reality (ISMAR), 2011 10th IEEE international symposium on. IEEE, 2011, pp. 127–136.
-  M. Klingensmith, I. Dryanovski, S. Srinivasa, and J. Xiao, “Chisel: Real time large scale 3d reconstruction onboard a mobile device using spatially hashed signed distance fields.” in Robotics: Science and Systems, vol. 4, 2015.
-  H. Oleynikova, Z. Taylor, M. Fehr, J. Nieto, and R. Siegwart, “Voxblox: Incremental 3d euclidean signed distance fields for on-board mav planning,” arXiv preprint arXiv:1611.03631, 2016.
-  H. Cover, S. Choudhury, S. Scherer, and S. Singh, “Sparse tangential network (spartan): Motion planning for micro aerial vehicles,” in 2013 IEEE International Conference on Robotics and Automation. IEEE, 2013, pp. 2820–2825.
-  J. Redmon and A. Farhadi, “Yolov3: An incremental improvement,” CoRR, vol. abs/1804.02767, 2018. [Online]. Available: http://arxiv.org/abs/1804.02767
-  S. Ren, K. He, R. B. Girshick, and J. Sun, “Faster R-CNN: towards real-time object detection with region proposal networks,” CoRR, vol. abs/1506.01497, 2015. [Online]. Available: http://arxiv.org/abs/1506.01497
M. Raza, Z. Chen, S.-U. Rehman, P. Wang, and P. Bao, “Appearance based pedestrians’ head pose and body orientation estimation using deep learning,”Neurocomputing, vol. 272, pp. 647–659, 2018.
S. Li and A. Chan, “3d human pose estimation from monocular images with deep convolutional neural network,” vol. 9004, 11 2014, pp. 332–347.
-  A. G. Howard, M. Zhu, B. Chen, D. Kalenichenko, W. Wang, T. Weyand, M. Andreetto, and H. Adam, “Mobilenets: Efficient convolutional neural networks for mobile vision applications,” arXiv preprint arXiv:1704.04861, 2017.
-  X. Zhang, X. Zhou, M. Lin, and J. Sun, “Shufflenet: An extremely efficient convolutional neural network for mobile devices,” in , June 2018.
-  W. Wang, A. Ahuja, Y. Zhang, R. Bonatti, and S. Scherer, “Improved generalization of heading direction estimation for aerial filming using semi-supervised regression,” Robotics and Automation (ICRA), 2019 IEEE International Conference on, 2019.
-  T.-Y. Lin, M. Maire, S. Belongie, J. Hays, P. Perona, D. Ramanan, P. Dollár, and C. L. Zitnick, “Microsoft coco: Common objects in context,” in European conference on computer vision. Springer, 2014, pp. 740–755.
-  J. F. Henriques, R. Caseiro, P. Martins, and J. Batista, “High-speed tracking with kernelized correlation filters,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 37, no. 3, pp. 583–596, 2015.
-  A. Toshev and C. Szegedy, “Deeppose: Human pose estimation via deep neural networks,” in Proceedings of the IEEE conference on computer vision and pattern recognition, 2014, pp. 1653–1660.
-  Z. Cao, T. Simon, S.-E. Wei, and Y. Sheikh, “Realtime multi-person 2d pose estimation using part affinity fields,” in CVPR, 2017.
-  M. Zucker, N. Ratliff, A. D. Dragan, M. Pivtoraiko, M. Klingensmith, C. M. Dellin, J. A. Bagnell, and S. S. Srinivasa, “CHOMP: Covariant hamiltonian optimization for motion planning,” The International Journal of Robotics Research, vol. 32, no. 9-10, pp. 1164–1193, 2013.