user: text: Write a concise one-sentence TLDR summary of the scientific paper content below.

Title: Stochastic Prediction of Multi-Agent Interactions from Partial Observations

Paper content:
We present a method which learns to integrate temporal information, from a learned dynamics model, with ambiguous visual information, from a learned vision model, in the context of interacting agents.
Our method is based on a graph-structured variational recurrent neural network, which is trained end-to-end to infer the current state of the (partially observed) world, as well as to forecast future states.
We show that our method outperforms various baselines on two sports datasets, one based on real basketball trajectories, and one generated by a soccer game engine.
Imaging watching a soccer game on television.
At any given time, you can only see a subset of the players, and you may or may not be able to see the ball, yet you probably have some reasonable idea about where all the players currently are, even if they are not in the field of view.
(For example, the goal keeper is probably close to the goal.)
Similarly, you cannot see the future, but you may still be able to predict where the "agents" (players and ball) will be, at least approximately.
Crucially, these problems are intertwined: we are able to predict future states by using a state dynamics model, but we can also use the same dynamics model to infer the current state of the world by extrapolating from the last time we saw each agent.
In this paper, we present a unified approach to state estimation and future forecasting for problems of this kind.
More precisely, we assume the observed data consists of a sequence of video frames, v 1:T , obtained from a stationary or moving camera.
The desired output is a (distribution over a) structured representation of the scene at each time step, p(s t |v 1:t ), as well as a forecast into the future, p(s t+∆ |v 1:t ), where s k t encodes the state (e.g., location) of the k'th agent and s t = {s DISPLAYFORM0 The classical approach to this problem (see, e.g., BID3 ) is to use state-space models, such as Kalman filters, for tracking and forecasting, combined with heuristics, such as nearest neighbor, to perform data association (i.e., inferring the mapping from observations to latent objects).
Such generative approaches require a dynamical model for the states, p(s t |s t−1 ), and a likelihood model for the pixels, p(v t |s t ).
These are then combined using Bayes' rule.
However, it is hard to learn good generative model of pixels, and inverting such models is even harder.
By contrast, our approach is discriminative, and learns an inference network to compute the posterior belief state p(s t |v 1:t ) directly.
In particular, our model combines ideas from graph networks, variational autoencoders, and RNNs in a novel way, to create what we call a graph-structured variational recurrent neural network (Graph-VRNN).We have tested our approach on two datasets: real basketball trajectories, rendered as a series of (partially observed) bird's eye views of the court; and a simple simulated soccer game, rendered using a 3d graphics engine, and viewed from a simulated moving camera.
We show that our approach can infer the current state more accurately than other methods, and can also make more accurate future forecasts.
We also show that our method can vary its beliefs in a qualitatively sensible way.
For DISPLAYFORM1 Figure 1: Illustration of visual VRNN with a single agent.
Dotted edges are not used.
Dashed edges are non-standard edges that we add.example, it "knows" the location of the goalie even if it is not visible, since the goalie does not move much (in our simulation).
Thus it learns to "see beyond the pixels".In summary, our main contribution is a unified way to do state estimation and future forecasting at the level of objects and relations directly from pixels using Graph-VRNN.
We believe our technique will have a variety of other applications beyond analysing sports videos, such as self-driving cars (where inferring and predicting the motion of pedestrians and vehicles is crucial), and human-robot interaction.
In this section, we briefly review some related work.
Graph-structured and stochastic RNNs.
There are several papers that combine standard RNNs (recurrent neural networks) with graph-structured representations, such as BID9 BID4 .
There are also several papers that extend standard RNNs by adding stochastic latent variables, notably the variational RNN approach of BID10 , as well as recent extensions, such BID27 BID15 ; BID22 BID17 ; BID8 .
However, as far as we know, combining VRNNs with graphs is novel.
Predicting future pixels from past pixels.
There are many papers that try to predict the future at the pixel level (see e.g., BID25 for a review).
Some use a single static stochastic variable, as in a conditional VAE BID23 , which is then "decoded" into a sequence using an RNN, either using a VAE-style loss BID43 BID47 or a GAN-style loss BID42 BID32 .
More recent work, based on VRNNs, uses temporal stochastic variables, e.g., the SV2P model of BID2 and the SVGLP model of BID11 .
There are also various GAN-based approaches, such as the SAVP approach of BID28 and the MoCoGAN approach of BID38 .
The recent sequential AIR (attend, infer, repeat) model of BID26 uses a variable-sized, object-oriented latent state space rather than a fixed-dimensional unstructured latent vector.
This is trained using a VAE-style loss, where p(v t |s t ) is an RNN image generator.
Forecasting future states from past states.
There are several papers that try to predict future states given past states.
In many cases the interaction structure between the agents is assumed to be known (e.g., using fully connected graphs or using a sparse graph derived from spatial proximity), as in the social LSTM BID0 and the social GAN BID18 , or methods based on inverse optimal control BID25 BID29 .
In the "neural relational inference" method of BID24 , they infer the interaction graph from trajectories, treating it as a static latent variable in a VAE framework.
BID13 proposed to predict a Gaussian distribution for the uncertain future states of a sphere.
It is also possible to use graph attention networks BID40 for this task (see e.g., BID21 ).
Forecasting future states from past pixels.
Our main interest is predicting the hidden state of the world given noisy visual evidence.
BID41 ) train a conditional CNN to predict future feature vectors, and several papers (e.g., BID30 BID33 BID45 BID39 BID16 BID44 ) predict future object states.
Our work differs by allowing the vision and dynamics model to share and update a common belief space; also, our model is stochastic, not deterministic.
Our method is related to the "deep tracking" approach of BID12 .
However, they assume the input is a partially observed binary occupancy grid , and the output is a fully observed binary occupancy grid, whereas we predict the state of individual objects using a graph-structured model.
Our method is also related to the "backprop Kalman filtering" approach of BID19 .
However, they assume the structure and dynamics of the latent state is known, and they use the Kalman filter equations to integrate information from the dynamical prior and (a nonlinear function of) the observations, whereas we learn the dynamics, and also learn how to integrate the two sources of information using attention.
We start by reviewing the VRNN approach of BID10 .
This model has three types of variable: the observed output x t , the stochastic VAE state z t , and the deterministic RNN hidden state h t , which summarizes x ≤t and z ≤t .
We define the following update equations: DISPLAYFORM0 DISPLAYFORM1 DISPLAYFORM2 where ϕ are neural networks (see section 3.4 for the details).VRNNs are trained by maximizing the evidence lower bound (ELBO), which can be decomposed as follows: DISPLAYFORM3 We use Gaussians for the prior and posterior, so we can leverage the reparameterization trick to optimize this objective using SGD, as explained in BID23 .
In practice, we scale the KL divergence term by β, which we anneal from 0 to 1, as in BID7 .
Also, we start by using the ground truth values of x t−1 (teacher forcing), and then we gradually switch to using samples from the model; this technique is known as scheduled sampling , and helps the method converge.
We can create a conditional generative model of the form p(s 1:t |v 1:t ) by using a VRNN (where x t = s t is the generated output) by making all the equations in section 3.1 depend on v t as input.
This is similar to a standard sequence-to-sequence model, but augmented to the graph setting, and with additional stochastic noise injected into the latent dynamics.
We choose to only condition the outputs s t on the visual input v t , rather than making the hidden state, h t , depend on v t .
The reason is that we want to be able to perform future prediction in the latent state space without having to condition on images.
In addition, we want the latent noise z t to reflect variation of the dynamics, not the appearance of the world.
See fig. 1 for a simplified illustration of our model (ignoring the graph network component).
All parameters are shared across all agents except for the visual encoder; they have a shared feature extraction "backbone", but then use agent-specific features to identify the specific agent.
(Agents have the same appearance across videos.)
Thus the model learns to perform data association.
The key issue is how to define the output decoder, p(s t |z t , h t−1 , v t ), in such a way that we properly combine the information from the current (partially observed) visual input, v t , with our past beliefs about the object states, h t−1 , as well as any stochastic noise z t (used to capture any residual uncertainty).
In a generative model, this step would be performed with Bayes' rule, combining the dynamical prior with the visual likelihood.
In our discriminative setting, we can learn how to weight the information sources using attention.
In particular, we define DISPLAYFORM0 where ϕ DV is the visible decoder, ϕ DH is the hidden decoder, and the α i terms are attention weights, computed as follows: DISPLAYFORM1 where i ∈ {V, H}.This is similar to the gating mechanism in an LSTM, where we either pass through the current observations, or we pass through the prior hidden state, or some combination of the two.
If the visual input is uninformative (e.g., for future frames, we set v t+∆ = 0), the model will rely entirely on its dynamics model.
To capture uncertainty more explicitly, we also pass in s heat t into the ϕ Si functions that compute attention, where s heat t = softmax(ϕ dec (z t , h t−1 , v t )) is the set of "heatmaps" over object locations.
(This is illustrated by the s t−1 → s t edges in fig. 1.)
We also replace the sample z t with its sufficient statistics, (µ t , Σ t ), computed using the state dependent prior, ϕ prior (h t−1 ).To encourage the model to learn to forecast future states, in addition to predicting the current state, we modify the above loss function to maximize a lower bound † on log p(s 1:T +∆ |v 1:T ), computed as follows: DISPLAYFORM2 where λ t = max(t/T, 1) is a weighting factor that trades off state estimation and forecasting losses, DISPLAYFORM3 is the variational posterior, the expectation is over z 1:T , and where we set v τ = 0 if τ > t.
For the image encoder, we use convolutional neural networks with random initialization.
In the case when there are multiple frames per step, we combine the image network with an S3D BID46 inspired network, which adds additional temporal convolutions on spatial feature maps.
For the recurrent model, we tried both vanilla RNNs and GRUs, and found that results were similar.
We report results using the GRU model.
For the graph network, we use relation networks BID37 , i.e. a fully connected topology with equal edge weights.
During training, both state estimation and future prediction losses are important, and need to be weighted properly.
When the future prediction loss weight is too small, the dynamics model is not trained well (reflected by log-likelihood).
We find that scaling the losses using an exponential discount factor for future predictions is a good heuristic.
We normalize the total weights for state estimation losses and future predictions losses to be the same.
In this section, we show results on two datasets, comparing our full model to various baselines.
Since we are interested in inferring and forecasting the state of a system composed of multiple interacting objects, based on visual evidence, analysing sports videos is a natural choice.
We focus on basketball and soccer, since these are popular games for which we can easily get data.
The states here are the identities of the objects (players and ball), and their (x, y) locations on the ground plane.
Although specialized solutions to state estimation for basketball and soccer already exist (e.g., based on wearable sensors BID36 , multiple calibrated cameras BID31 , or complex monocular vision pipelines BID34 ), we are interested in seeing how far we can get with a pure learning based approach, since such a method will be more generally applicable.
Basketball.
We use the basketball data from BID48 .
Each example includes the trajectories of 11 agents (5 offensive players, 5 defensive players and 1 ball) for 50 steps.
We follow the standard setup from BID48 and just model the 5 offensive players, ignoring the defense.
However, we also consider the ball, since it has very different dynamics than the players.
Overall, there are 107,146 training and 13,845 test examples.
We generate bird's-eye view images based on the trajectories (since we do not have access to the corresponding video footage of the game), where each agent is represented as a circle color-coded by its identity.
To simulate partial observation, we randomly remove one agent from the rendered image every 10 steps.
Some example rendered images can be found in FIG3 .Soccer.
To evaluate performance in a more visually challenging environment, we consider soccer, where it is much harder to get a bird's-eye view of the field, due to its size.
Instead most soccer videos are recorded from a moving camera which shows a partial profile view of the field.
Since we do not have access to ground truth soccer trajectories, we decided to make our own soccer simulator, which we call Soccer World, using the Unity game engine.
The location of each player is determined by a hand-designed "AI", which is a probabilistic decision tree based on the player's current state, the states of other players and the ball.
The players can take actions such as kicking the ball and tackling.
The ball is driven by the players' actions.
Each game lasts 5 minutes.
For each game, a player with the same id is assigned to a random position to play, except for the two goal keepers.
We render the players using off-the-shelf human models from Unity.
The identities of the players are color-coded on their shirts.
The camera tracks the ball, and shows a partial view of the field at each frame.
See FIG3 for some visualizations.
We create 700 videos for training, and 300 videos for test.
We apply a random sliding window of 10 seconds to sample the training data.
The test videos are uniformly segmented into 10-second clips, resulting in 9000 total test examples.
We plan to release the videos along with the game engine after publication of the paper.
‡
The first task we evaluate is inferring the current state (which we define to be the 2d location) of the all the objects (players and ball).
To do this, we replace the discrete predictions with the corresponding real-valued coordinates by using weighted averaging.
We then compute normalized 2 distance between the ground truth locations and the predicted locations of all objects at each step.
The second task is predicting future states.
Since there are many possible futures, using 2 loss does not make sense, either for training or evaluation (c.f., BID32 ).
Instead we evaluate the negative log-likelihood on the discretized predictions, which can better capture the multi-modal nature of the problem. (This is similar to the perplexity measure used to evaluate language models, except we apply it to each object separately and sum the results.)For each task, we consider the following models: Visual only: standalone visual encoder without any recurrent neural network as backbone.
For future prediction, we follow BID14 and directly predict future locations from last observed visual features.
RNN: standard RNN, the hidden states are globally shared by all agents.
VRNN: standard VRNN.
Indep-RNN: one RNN per agent, each RNN runs independently without any interaction.
Social-RNN: one RNN per agent, agents interact via the pooling mechanism in Social GAN BID18 .
Graph-RNN: one RNN per agent, with a graph interaction network.
Graph-VRNN: the full model, which adds stochasticity to Graph-RNN.
All models have the same architectures for visual encoders and state decoders.
Our visual encoder is based on ResNet-18 BID20 , we use the first two blocks of ResNet to maintain spatial resolution, and then aggregate the feature map with max pooling.
The encoder is pre-trained on visible players, and then fine-tuned for each baseline.
We find this step important to stabilize training.
For the soccer data, we down-sample the video to 4 FPS, and treat 4 frames (1 second) as one step.
We consider 10 steps in total, 6 observed, 4 unobserved.
We set the size of GRU hidden states to 128 for all baselines.
The state decoder is a 2-layer MLP.
For basketball data, we set every 5 frames as one step, and consider 10 steps as well.
The size of GRU hidden states is set to 128.
The parameters of the VRNN backbone and state decoders are shared across all agents, while each agent has its own visual encoder as "identifier".
For all experiments, we use the standard momentum optimizer.
The models are trained on 6 V100 GPUs with synchronous training with batch size of 8 per GPU, we train the model for 80K steps on soccer and 40K steps on basketball.
We use a linear learning rate warmup schedule for the first 1K steps, followed by a cosine learning rate schedule.
‡ Video samples can be found at bit.ly/2E3qg6FMethod t = 1 t = 2 t = 3 t = 4 t = 5 t = 1 t = 2 t = 3 t = 4 t = 5 Table 1 : Normalized 2 distances to ground truth locations at different steps on (left) soccer data and (right) basketball data.
The hyper parameters, such as the base learning rate and the KL divergence weight β, are tuned on a hold-out validation set.
Table 3 : Future prediction performance as measured by log-likelihood ratio over random guess.
Basketball.
The right half of Table 1 shows the average normalized 2 distance between the true and predicted location of all the agents for the basketball data.
(Error bars are not shown, but variation across trials is quite small.)
We see that the Graph-VRNN error generally goes down over time, as the system integrates evidence from multiple frames. (Recall that one of the agents becomes "invisible" every 10 frames or 2 steps, so just looking at the current frame is insufficient.)The visual-only baseline has more or less constant error, which is expected since it does not utilize information from earlier observations.
All other methods outperform the visual-only baseline.
We can see that stochasticity helps (VRNN better than RNN, Graph-VRNN better than Graph-RNN), and Graph-RNN is better than vanilla RNN.
As most of the players are visible in the basketball videos, we also report performance for visible and hidden (occluded) agents in TAB1 .
As expected, the 2 distances for visible agents are very low, since the localization task is trivial for this data.
When the agents are hidden, Graph-VRNN significantly outperforms other baselines.
Again, we observe that stochasticity helps.
We also find that graph network is a better interaction model than social pooling, both outperform Indep-RNN.Soccer.
The left half of Table 1 shows the average normalized 2 distance between the true and predicted location of all the agents for the soccer data as a function of time.
The results are qualitatively similar to the basketball case, although in this setting the 2 distances are higher since the vision problem is much harder.
However, the gains from adding stochasticity to the Graph-RNN are smaller in this setting.
We believe the reason for this is that the dynamics of the agents in soccer world is much more predictable than in the basketball case, because our simulator is not very realistic. (This will become more apparent in Section 4.4.)Forecasting.
In this section, we assess the ability of the models to predict the future.
In particular, the input is the last 6 steps, and we forecast the next 4 steps.
Since the future is multimodal, we do not use 2 error, but instead we compute the log-likelihood of the discretized ground truth locations, normalized by the performance of random guessing § , which is a uniform prior over all locations.
The prior baseline trains the decoder based on constant zero input vectors, thus reflecting the prior marginal distribution over locations.
The visual only baseline predicts all 4 future steps based on last observed visual feature, as is in BID14 .
The results are shown in Table 3 .
Not surprisingly, we see that modeling the interaction between the agents helps, and adding stochasticity to the latent dynamics also helps.
This section shows some qualitative visualizations of the belief state of the system (i.e., the marginal distributions p(s Soccer.
FIG7 (left) shows the belief states for the soccer domain for three different kinds of agents: a regular player (top row), the goal keeper (middle row), and the ball (bottom row).
For the player, we see that the initial belief state reflects the 4-4-2 pattern over possible initial locations of the players.
In frame 3, enough visual evidence has been accumulated to localize the player to one of two possible locations.
For the future prediction, we draw a single stochastic sample (to break symmetry), and visualize the resulting belief states.
We see that the model predicts the player will start moving horizontally (as depicted by the heatmap diffusing), since he is already at the top edge of the field.
The goal keeper beliefs are always deterministic, since in the simulated game, the movement of the goalie is below the quantization threshold.
The ball is tracked reliably (since the camera is programmed to track the ball), but future forecasts start to diffuse.
FIG7 (right) shows similar results for the basketball domain.
We see that the dynamics of the players and ball are much more complex than in our soccer simulator.
We have presented a method that learns to integrate temporal information with partially observed visual evidence, based on graph-structured VRNNs, and shown that it outperforms various baselines on two simple datasets.
In the future, we would like to consider more challenging datasets, such as real sports videos.
We would also like to reduce the dependence on labeled data, perhaps by using some form of self-supervised learning.
A.1 SAMPLED SOCCER TRAJECTORIES Figure A1 : Sampled trajectories for Soccer World for 11 "home" players.
Figure A1 shows the sampled trajectories for Soccer World.
Unlike basketball, Soccer World has a moving camera with limited field of view.
We observe that the trajectories for the first few observed steps are quite shaky since only a few players have been observed.
We thus show the trajectories from t = 5.
From fig. A1 , we can see that the trajectories generated by RNN are very shaky.
Graph-VRNN generates much better trajectories, but some of the players are assigned incorrect identities, or incorrect locations.
We conjecture that this issue can be mitigated by providing longer visual inputs to the model, such that most of the players could be observed at some point of the videos.
A.2 SAMPLED BASKETBALL TRAJECTORIES FIG1 and Figure

assistant: text: The authors propose Graph VRNN which models the interaction of multiple agents by deploying a VRNN for each agent