Next Article in Journal
Research on Data Poisoning Attack against Smart Grid Cyber–Physical System Based on Edge Computing
Previous Article in Journal
FA-Net: A Fused Feature for Multi-Head Attention Recoding Network for Pear Leaf Nutritional Deficiency Diagnosis with Visual RGB-Image Depth and Shallow Features
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Shooting Events in Soccer Videos Using Complete Bipartite Graphs and Players’ Spatial-Temporal Relations

1
Graduate School of Information Science and Technology, Hokkaido University, N-14, W-9, Kita-ku, Sapporo 060-0814, Hokkaido, Japan
2
Faculty of Information Science and Technology, Hokkaido University, N-14, W-9, Kita-ku, Sapporo 060-0814, Hokkaido, Japan
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(9), 4506; https://doi.org/10.3390/s23094506
Submission received: 3 April 2023 / Revised: 1 May 2023 / Accepted: 3 May 2023 / Published: 5 May 2023
(This article belongs to the Section Sensing and Imaging)

Abstract

:
In soccer, quantitatively evaluating the performance of players and teams is essential to improve tactical coaching and players’ decision-making abilities. To achieve this, some methods use predicted probabilities of shoot event occurrences to quantify player performances, but conventional shoot prediction models have not performed well and have failed to consider the reliability of the event probability. This paper proposes a novel method that effectively utilizes players’ spatio-temporal relations and prediction uncertainty to predict shoot event occurrences with greater accuracy and robustness. Specifically, we represent players’ relations as a complete bipartite graph, which effectively incorporates soccer domain knowledge, and capture latent features by applying a graph convolutional recurrent neural network (GCRNN) to the constructed graph. Our model utilizes a Bayesian neural network to predict the probability of shoot event occurrence, considering spatio-temporal relations between players and prediction uncertainty. In our experiments, we confirmed that the proposed method outperformed several other methods in terms of prediction performance, and we found that considering players’ distances significantly affects the prediction accuracy.

1. Introduction

Data mining and machine learning technologies have made data analytics increasingly crucial in the field of sports field [1]. Quantitative evaluation of players’ and teams’ performance is a major task for enhancing competitiveness, and several methods have been proposed for various sports such as baseball [2,3] and basketball [4,5]. In soccer, such quantitative evaluation is expected to improve tactical coaching [6,7] and help players in making decisions [8,9,10]. In particular, the importance of data analytics in soccer has increased further since the International Football Association Board (IFAB) (https://www.theifab.com/, accessed on 1 February 2023) allowed the use of electronic communication devices by the coaching staff for tactical coaching during matches in 2018 [11], making it even more crucial. However, as soccer is a low-scoring and dynamic team sport, it is challenging to determine whether individual player actions contribute to match results or not [12]. Consequently, there is a need for quantitative evaluation methods that consider the complexity of soccer analytics. Soccer analytics researchers have conventionally proposed several methods to quantify the actions of each player and the team strategy by utilizing the occurrence probability of specific events such as shots and passes [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. Since the performance metrics based on event probability are robust in the face of soccer-specific chance, they represent the team’s performance more reliably than traditional performance metrics such as ball possession, the number of shots on goal and pass rates [13]. In particular, the probability of a shoot event occurrence is often used as a substitute for the quantitative evaluation of both offensive and defensive team play since shoot events have the potential to directly impact the scoreline [12]. The conventional methods have mainly achieved the prediction of the event occurrence probability using tracking data, which consists of time-series data showing the two-dimensional positions of players on the field. Studies have shown that prediction methods based on both video and tracking data for sports such as water polo and basketball are more accurate than those relying on only one type of data [28]. Despite these findings, traditional prediction methods in soccer have not integrated visual information obtained from video data. Thus, incorporating both video and tracking data in a new prediction method is likely to lead to more accurate predictions in soccer. Moreover, while conventional deep learning methods have achieved remarkable success in various fields, some of these methods do not consider the uncertainty in their predictions [24,25,26,27]. Consequently, the generalization ability of these methods can decrease due to overfitting, leading to decreased reliability in the prediction results and overconfidence in their accuracy. Therefore, it becomes difficult to make correct decisions when introducing unreliable prediction results into tactical coaching in soccer. Since uncertainty has several aspects, such as the lack of knowledge about the predicted event, it is challenging to estimate prediction uncertainty through statistical inferences based on stochastic models. Thus, there is a need for a method that enables the evaluation of prediction uncertainty. Based on the aforementioned issues, the problems that need to be addressed are as follows.
(1)
Conventional methods solely concentrate on the positional relations among players, without incorporating any visual information from video data.
(2)
Conventional methods fail to account for the reliability of event occurrence probabilities derived from deep learning-based methods.
In this paper, we introduce a novel approach that predicts shoot events by considering players’ spatial-temporal relations in soccer videos through a complete bipartite graph. The occurrence of shoot events in soccer matches is closely linked to the position and movement of the players, and modeling spatio-temporal relations can aid in predicting them.
Therefore, our method utilizes a graph to represent players’ relations in each video frame, effectively incorporating the soccer domain knowledge into the treatment of visual information obtained from video data. We employ a graph convolutional recurrent neural network (GCRNN) [29,30], which which has fused a graph convolutional network (GCN) [31,32] and a recurrent neural network (RNN) [33,34], to learn the spatial and temporal relations between players in the graph. The proposed method uses visual information obtained from the video data as nodes and constructs a complete bipartite graph based on players’ distances and team information. It is worth noting that broadcast videos, which emphasize the ball carrier and include dissolves and replay scenes, may not capture the complete spatio-temporal information of players. Therefore, we focused on unedited videos that provide a birds-eye view of the entire match, which are commonly referred to as scouting videos. Overall, our method can accurately predict shooting events by effectively incorporating the spatio-temporal relations of almost all players in soccer videos. Moreover, to address the issue of prediction uncertainty, we incorporate a Bayesian neural network (BNN) [35,36] into our method. This approach allows for the quantitative evaluation of both model- and data-dependent uncertainties by treating the neural network’s weight parameters as random variables. The proposed method enables the prediction of the probability of shoot event occurrence and also allows for the calculation of the degree of uncertainty associated with making the prediction. Note that this study is an extension of our previous work [37], in which we proposed a more accurate prediction model by incorporating a complete bipartite graph with soccer-specific domain knowledge, including players’ team information. The main contributions of this study are as follows.
(1)
Our novel method for predicting shoot events in soccer videos incorporates player-specific features and spatio-temporal relations using GCRNN.
(2)
We propose a method that aims to enhance the robustness of predictions by providing a measure of prediction uncertainty through the use of BNN.

2. Pre-Processing

This section describes the pre-processing methods used to construct complete bipartite graphs based on the detection of players in a soccer video. The pre-processing flow is illustrated in Figure 1. Initially, we employ a camera calibration technique to localize players on the soccer field, as explained in (Section 2.1). Furthermore, we present a technique for classifying team classification of detected players in the video (in Section 2.2).

2.1. Camera Calibration

To obtain the two-dimensional coordinates of the players detected in the soccer video on the field, it is necessary to estimate the camera parameters for the homography transformation. Then we use stroke width transform (SWT) [38] as an edge detector to detect the field lines on the soccer field. SWT is effective not only in preventing sparse detections but also in being robust to difficult viewpoints and noise effects such as shadows, and field stripes (both bright and dark green stripes). Therefore, it has already been confirmed to be effective in extracting the white line region on the soccer field [39]. After applying the SWT to a full-frame image (a), the resultant image (b) contains numerous edges that are irrelevant to the field lines that require removal. Consequently, we perform certain processes to eliminate the edges that are not associated with the field lines, as follows:
  • Since field lines are believed to maintain a small and consistent stroke width, we apply a filter to include only strokes that are longer than 10 pixels.
  • To segment strokes based on a the field color, we extract only those strokes within the 36 < h < 90 range in the HSV color space, as the field lines are typically found on the grass.
  • We use a field mask (b) to mask off-field areas such as spectator seats and benches.
  • To eliminate player regions, we use a Mask R-CNN detector [40].
The field mask (c) is obtained using the segmentation image output of the generative network from the two-conditional GAN model for detecting soccer field lines [41] with the full-frame (a) as input. The generative network of the two-conditional GAN model is based on the pix2pix network architecture [42], and it uses modules of the form convolution BatchNorm-Relu and follows U-Net [43]. Although this model is capable of generating the edge image of the field line from broadcast images, it did not work well with the scouting video frame used in this study. Therefore, we only used the segmentation image as the field mask. Finally, we can obtain an edge image (e) of the soccer field by conducting a dilation operation to remove noise.
As a solution for capturing camera parameters necessary for homography transformation, we propose an approach based on HOG features [44] to search for nearest neighbors in a synthetic edge image dictionary. The proposed architecture, involves finding the most similar edge image from the dictionary, giving an input edge image (e), and outputting the corresponding camera parameters. For this study, we use a perspective projection camera model (refer to Figure 2) to generate the projection-transformed image of the field template. This approach considers the virtual image plane in the pinhole camera model as the true image plane and generates images that are used as synthetic field edge images.
For the field template, we have used the standard field (refer to Figure 3), which is commonly used in international matches like the World Cup and the Olympic Games, and is defined by the Federation Internationale de Football Association (FIFA) (https://www.fifa.com/, accessed on 1 February 2023). The equation for the projection transformation is presented below:
Z n x n y n 1 = f l x 0 u 0 f l y v 0 0 1 x c y c z c = f l x 0 u 0 f l y v 0 0 1 R ( θ ) x w c x w y w c y w z w c z w ,
where x w , y w , z w and c x w , c y w , c z w represent the field edge coordinates and camera coordinates in the world coordinate system, respectively. R ( θ ) θ = θ x , θ y , θ z denotes the Rodrigues rotation matrix with the camera rotation angle θ . The camera parameters obtained above are used to calculate the field edge coordinates x c , y c , z c in the camera coordinate system. The value of f l x , f l y and u , v are determined based on the camera scaling and the size of the projection-transformed image, respectively. As a result of the aforementioned projection transformation, we obtain the field edge coordinates x n , y n in the image coordinate system, which are then used to retrieve the field edge image. These camera parameters are constrained within specific ranges based on the environmental conditions to be captured. In this study, our focus is on scouting soccer videos during a home game of a team in the Japan J1 League. Therefore, we experimentally determined the camera position to be c x w , c y w , c z w = 7 , 105 , 36 . For the camera rotation angles, we limited the ranges to θ x [ 75 , 65 ] , θ y [ 20 , 20 ] and θ z [ 20 , 20 ] , respectively. f l x , f l y were restricted to the range f l x = f l y 1200 , 2000 to ensure that at least 1/3 of the field fits within the template projection image. u , v were set to 640 and 360, respectively. In addition, each camera parameter θ x , θ y , θ z , f l x , f l y is varied in increments of δ ( θ x ) = δ ( θ y ) = δ ( θ z ) = 1 and δ ( f l x ) = δ ( f l y ) = 50 , respectively, within the limited range ( θ x [ 75 , 65 ] , θ y [ 20 , 20 ] , θ z [ 20 , 20 ] , and f l x = f l y [ 1200 , 2000 ] ). Here, δ ( · ) represents the stride of each camera parameter. After computing the HOG features of the template projection images and the edge image (d), we obtain the camera parameters of the most similar template projection image. That is the optimal ones, by performing nearest neighbor search using the Euclidean distance of those features.

2.2. Team Classification

This subsection describes the method used to classify players into teams, which is necessary for constructing a complete bipartite graph. To detect players in each frame, we use the well-known object detection method Mask R-CNN [40]. Mask R-CNN is a single-task framework that extends Faster R-CNN [45] and is capable of performing instance segmentation and object detection and instance simultaneously. Instance segmentation is a more challenging task than semantic segmentation as it involves classifying each pixel in an image according to both the object class and the specific instance it belongs to, whereas semantic segmentation only classifies each pixel according to its object class. Thus, we can obtain the region of each player by applying Mask R-CNN, and from the pixels of those players’ region, our method identifies the team to which each player belongs using a color-based classification method [46] that is based on the HSV uniform color. To minimize the effects of unrelated people such as referees, bench players, coaches, and photographers, we only include the field players (i.e., excluding goalkeepers) of the home and away teams, respectively, in the construction of the complete bipartite graph. These players are identified from the detected persons in the video frame. Here, it is worth nothing that while referees and goalkeepers can be distinguished based on their uniform color, it may not be possible to identify substitute players and coaches using this approach. To address this issue, our method determines whether a detected person is a field player or not by utilizing the coordinates of the detected person on the field, which are obtained through camera calibration.
Thus, the field players’ coordinates on the field can be obtained by using the obtained camera parameters obtained from the previous steps and mapping the field player’s coordinates in the video frame to the field (as shown in the lower right of Figure 1).

3. Method for Prediction of Shoot Events

In this section, we provide a detailed explanation of the proposed method for predicting shoot events, as illustrated in Figure 4. Specifically, our approach involves constructing a complete bipartite graph using visual features of player regions and a full-frame at each time step, as described in Section 3.1. Next, we use a graph convolutional recurrent neural network (GCRNN) [29,30] to compute latent features by incorporating spatio-temporal relations between players, as described in Section 3.2. Finally, we make use of a Bayesian neural network (BNN) [35,36] to predict the probability of a shooting event occurring while considering the prediction uncertainty. This is further elaborated in Section 3.3.

3.1. Construction of Complete Bipartite Graphs

This subsection explains the construction process of a complete bipartite graph. To effectively incorporate visual information obtained from soccer video data along with domain knowledge, our proposed method represents the relationships between players in each video frame as a complete bipartite graph G t . Initially, we detect players in soccer videos using the widely-used object detection and segmentation method, Mask R-CNN [40]. After detection, we classify players into teams based on their uniform colors, as described in See Section 2.2). Next, to extract the visual features of each player and the full frame at time step t, we feed the rectangular region of each player detected by Mask R-CNN and the full-frame at time step t into ResNet-101 [47], which is a pre-trained feature extractor with ImageNet [48]. The resulting features are denoted as P t and F t , respectively. To enhance the feature representation capacity of the graph data, we obtain low-dimensional features P ^ t and F ^ t of each player and full-frame using a fully connected layer. Then, we concatenated these features and introduce the resulting feature vectors as a graph node X t according to the following formula:
P ^ t = FC ( P t ) , F ^ t = FC ( F t ) ,
X t ( i ) = P ^ t ( i ) , F ^ t ,
where FC ( · ) denotes the fully connected layer, and P ^ t ( i ) are i -th player’s low-dimensional visual features at time step t.
It has been reported that in soccer matches, other players in the target player’s action range are especially crucial to make quick and accurate decisions [49]. We represent the graph G t with the adjacency matrix A t at time step t using edge weights that consider players’ distances and team information, which are crucial factors in soccer analysis. The proposed method specifically introduces a complete bipartite graph G t where edge weights are only between players from different teams. The adjacency matrix A t at time step t is then calculated as follows:
A t ( i j ) = 0 ( team ( i ) = team ( j ) ) exp { d ( r t i , r t j ) / w } i j exp { d ( r t i , r t j ) / w } ( team ( i ) team ( j ) ) ,
where r t denotes the position on the field of each player detected by Mask R-CNN (See Section 2.1) and d ( r t i , r t j ) represents the Euclidean distance between i-th and j-th players. The value w is a constant used to balance the edge weights, and the function team ( · ) is used to classify the teams of the detected players (refer to Section 2.2). The proposed method incorporates the distances and team information of both players into the graph, allowing for the utilization of player relations that accurately reflect the intricate match situation.

3.2. Calculation of Spatio-Temporal Relational Features

This subsection outlines the calculation method of spatio-temporal features between players using GCRNN, which combines a graph convolutional network (GCN) [31,32] and a recurrent neural network (RNN) [33,34]. Figure 5 illustrates the architecture of GCRNN. To learn the spatial relations between players, the proposed method utilizes stacked GCN layers. Additionally, the hidden state h t 1 of the RNN is introduced to enhance the second GCN layer, as shown bellow:
Z t = GCN ( [ GCN ( X t , A t ) , h t 1 ] , A t ) ,
where X t and A t represent the graph node and adjacency matrix, respectively, of the complete bipartite graph constructed in Section 3.1. The proposed method incorporates h t 1 into the final GCN layer, enabling the calculation of the latent feature Z t while considering both temporal and spatial relations. Furthermore, the proposed method inputs the latent feature Z t , which is fused with the corresponding graph node features X t , into the RNN. The hidden state h t is then updated for the next time step. As the number of players in a frame of the soccer video may vary across different time steps, it is essential to consider the dynamic changes of node features in the graph. To capture temporal relations, the proposed method utilizes a Graph Convolutional Recurrent Network (GCRN) [34], which is one of the RNN models. The hidden state h t is updated recurrently according to the following formula:
h t = GCRN ( [ Z t , X t ] , h t 1 ) .
Through the input of X t and Z t into the GCRN, the proposed method can obtain the hidden state h t while taking into account both player-specific and relational features. This iterative process allows our model to capture latent features Z t while considering spatio-temporal relations at each time step.

3.3. Event Prediction with Its Uncertainty

This subsection demonstrates the prediction of shoot event occurrences with their uncertainty based on BNN. BNN allows for the estimation of both model- and data-dependent uncertainty through the inference result by treating the weights of the neural network as random variables. As a result, whereas a simple neural network model is susceptible to incorrect predictions for extrapolation due to the overfitting problem, BNN can mitigate this issue. In particular, since soccer matches may not always follow the same pattern as the training data, overfitting the training data could lead to a decrease in the model’s ability to make accurate predictions. Therefore, by incorporating BNN, the proposed method enhances the accuracy of shoot event predictions in soccer matches. With the latent features Z ^ t of unknown data as input, the probability of a shoot event occurrence a ^ t can be calculated by the predictive distribution p ( a ^ t | Z ^ t ) using the following equation:
p ( a ^ t | Z ^ t ) = p ( a ^ t | Z ^ t , w ) p ( w | D ) d w ,
where p ( a ^ t | Z ^ t , w ) and p ( w | D ) represent the predictive distribution for a given value of network parameters w and posterior distribution of network parameters w , respectively. Network parameters w are sampled from a Gaussian distribution p ( w ) N ( μ , σ 2 ) , where the mean μ and variance σ 2 are learned using the training dataset D = ( Z t , a t ) .
To compute the true posterior distribution p ( w | D ) of network parameters using Bayes’ theorem, it is necessary to calculate the marginal distribution p ( D | w ) p ( w ) d w , in addition to the likelihood p ( D | w ) and the prior distribution p ( w ) of w . Because using a posterior distribution p ( w | D ) to obtain predicted probabilities is equivalent to employing an ensemble of numerous neural networks, this approach is not feasible for a neural network of a practical size. To address this issue, the proposed method utilizes Bayes-by-Backprop [36], a variational inference method that approximates the posterior distribution p ( w | D ) with the following equation:
α * = arg min α k = 1 K [ log q ( w k | α ) log p ( w k ) log p ( D | w k ) ] ,
where w k represents the k -th Monte Carlo sample drawn from the variational posterior q ( w k | α ) . The variational posterior parameters α = ( μ , σ 2 ) are updated using the re-parameterization trick and gradient descent. This technique is a well-known variance reduction method in Bayesian statistics and uncertainty quantification [50]. It is important to note that each term in this approximation cost is dependent on a specific weight obtained from the variational posterior. The proposed method consider the first and second terms as loss functions, denoted as L VPOS and L PRI , respectively. The third term refers to the negative log-likelihood, and minimizing this term results is equivalent to minimizing the mean squared error. Therefore, the proposed method uses the exponential binary cross entropy as the loss function L EXP , which is expressed as:
L E X P = t = 1 T e max ( 0 , T o t fps ) log a t ( p ) + t = 1 T log ( 1 a t ( n ) ) ,
where fps represents the frame rate of the input soccer videos, and T o represents the time of shoot events. Here, a t = ( a t ( n ) , a t ( p ) ) represents the negative and positive probability of the training dataset at time step t, respectively, measured in units of probability.
The proposed method involves multiple forwards passes at each time step to obtain multiple prediction distribution results, which are then used to calculate a prediction distribution. The mean and variance of this distribution are considered as the expected value and the prediction uncertainty, respectively. Furthermore, the predictive uncertainty, that is, the variance, can be decomposed into the aleatoric uncertainty U t a l t and the epistemic uncertainty U t e p t [51], allowing for a better understanding of the source of the uncertainty.
U t a l t = 1 M m = 1 M [ diag ( a ^ m ) a ^ m a ^ m ] ,
U t e p t = 1 M m = 1 M ( a ^ m a ¯ ) ( a ^ m a ¯ ) ,
where M denotes the number of forward passes. a ^ m = ( a ^ t ( n ) , a ^ t ( p ) ) m represents the output value obtained from the m-th layer in BNN, consisting of the predicted negative probability a ^ t ( n ) and positive probability a ^ t ( p ) of the unknown data at time step t. a ¯ is the average of output a ^ m from each layer of the BNN. The aleatoric uncertainty U t a l t is a measure of the inherent variability in the latent features calculated by the GCRNN, which serves as input to the BNN, and is independent of the specific model used. In contrast, the epistemic uncertainty U t e p t is a measure of the model’s uncertainty about its own predictions, and decreases as the model observes more training data, and confidence in the model is expected to increase. In the proposed method, we utilize the uncertainty-guided ranking loss L RANK [30] that leverages the epistemic uncertainty U t e p t to estimate the predictive capability of the model for the input data.
L R A N K = max ( 0 , trace ( U t e p t U t 1 e p t ) ) .
Finally, the prediction model is trained to minimize the coupling loss function L given by the following equation:
L = L EXP + λ 1 · ( L VPOS L PRI ) + λ 2 · L RANK ,
where the constants λ 1 and λ 2 are used to balance the values of the loss functions. By utilizing BNN, the proposed method can make predictions of shoot events while accounting for uncertainty. The proposed method combines latent features that consider the detailed relations between players with uncertainty modeling to accurately predict shoot event occurrences.

4. Experimental Results

This section presents experimental results that verify the effectiveness of the proposed method fin predicting shoot events. Subsequently, we will describe the experimental settings (in Section 4.1), followed by a performance evaluation and discussion (in Section 4.2).

4.1. Experimental Settings

For this experiment, we created a dataset consisting of 400 video clips extracted from scouting soccer videos of 31 matches played during the 2019 and 2020 seasons of the Japan J1 League. Each video clip contained an attacking scene captured at 10 fps for 6 s, with each frame having a size of 1280 × 720 pixels. Positive samples were defined as shooting scenes, while negative samples were defined as non-shooting attacking scenes. We randomly selected 80 video clips (40 positive and 40 negative samples) for the test dataset and used the remaining 320 video clips (160 positive and 160 negative samples) for training. The time step T o in which the shoot event occurs in positive samples is determined manually through our observations. To determine the probability a t = ( a t ( n ) , a t ( p ) ) of shoot events in the training data, we set a t = ( 1 , 0 ) at all time steps for negative samples. For positive samples, we set the probability as follows:
a t ( n ) = 1 a t ( p ) , a t ( p ) = 0 ( t < T o ) 1 ( t T o ) .
To balance the edge weights and loss scale values, we experimentally set w to 1050 in Equation (4) and λ 1 and λ 2 to 0.001 and 10, respectively, in Equation (13). The mixing ratio, which is a hyperparameter of the prior distribution, was set to 0.5. The number of forward passes, M, was determined by Bayes-by-Backprop [36] to be 2 during the training phase and 10 during the testing phase. In the training phase, we initialized the learning rate to 5 × 10 4 and used the ReduceLROnPlateau scheduler in conjunction with the Adam optimizer [52] to train the model with a batch size of 16 for 60 epochs. Moreover, we used the neural network with one layer as the FC layer in the graph construction for calculating P ^ t and F ^ t that had the 256 dimensions, that is, the dimension of X t became 512. In addition, we adopted the two GCN with one layer. The first GCN transformed graph node X t and adjacency matrix A t into the 256-dimensional feature, and the second GCN transformed the outputs of the first GCN with the hidden state h t 1 and adjacency matrix A t into the latent feature Z t . We set the dimensions of h t and Z t as 256 and 128, respectively. The GCRN in the GCRNN also had one layer that transformed the graph node X t , latent feature Z t and hidden state h t 1 into the next hidden state h t . Finally, we input the latent feature Z t and calculated the prediction results by using BNN with two layers that had 64 and 2 (probabilities of positive or negative shoot event occurrence) nodes, respectively. To evaluate the efficacy of our proposed method (Ours) in predicting shoot events, we implemented several comparative methods, namely ASs1-4 and CMs1-2.
AS1: 
This refers to our previous approach [37], which is similar to the proposed method but uses a complete graph without any team information.
AS2: 
This method builds upon our previous method but utilizes a complete bipartite graph that does not consider the distances between players.
AS3: 
This method is similar to our proposed method but employs a complete graph without incorporating edge weights in the graph representation.
AS4: 
This method uses the full-frame visual features, which applies a field mask obtained in Section 2.1 to the full-frame, as node features in the proposed method.
To our knowledge, our previous study [37] is the only method that predicts shoot event occurrence using visual information. The following additional methods were adopted as comparison methods, in addition to the comparative methods mentioned earlier.
CM1: 
This is a method for predictiong traffic accidents in dashcam videos using Dynamic-Spatial-Attention Recurrent Neural Network (DSA-RNN), as presented in the work of [53]. The DSA-RNN facilitates prediction by dynamically allocating soft-attention [54] to candidate objects in each frame, collecting cues, and learning the temporal relationships of all cues.
CM2: 
This is a state-of-the-art method for video classification that uses the Video Vision Transformer (ViViT) [55]. ViViT tokenizes the video by dividing it into either 2-D spatial or 3-D spatio-temporal grids. The video can be classified by feeding the tokens extracted from it into the transformer encoder.
CM1 is a prediction method for traffic accidents that considers the spatio-temporal relationships between objects in the video, and thus we included it as a benchmark method for comparison in this study. Furthermore, as the classification performance of the model is indicative of its ability to accurately predict shoot events, we also used CM2 as another benchmark method for comparison.
In our experiment, we assessed the prediction performance of shoot event occurrence using Average Precision (AP) and F1-score as indicators of prediction accuracy, and mean Time to Event (mTTE) as an indicator of prediction earliness for positive samples. The AP and F1-score are calculated as follows:
AP = 0 1 Precision ( s ) d Recall ( s ) ,
F 1 score ( s ) = 2 × Recall ( s ) × Precision ( s ) Precision ( s ) + Recall ( s ) ,
where s denotes the threshold used to determine whether the model predicted a future shoot event. Based on the prediction results, Recall ( s ) and Precision ( s ) were calculated as follows:
Precision ( s ) = TP ( s ) TP ( s ) + FP ( s ) ,
Recall ( s ) = TP ( s ) TP ( s ) + FN ( s ) ,
where TP ( s ) , FP ( s ) , and FN ( s ) represent the numbers of videos predicted as true-positive, false-positive and false-negative at threshold s, respectively. Details of the mTTE calculation are provided below.
mTTE = E [ TTE ( s ) ] ,
TTE ( s ) = max { T o t | r e s ( t ) s } ,
r e s ( t ) = exp ( a ^ t ( p ) ) exp ( a ^ t ( n ) ) + exp ( a ^ t ( p ) ) .

4.2. Performance Evaluation

4.2.1. Quantitative Results

The experimental results with varying video clip lengths are presented in Table 1. According to Table 1, the proposed method outperforms all other methods in terms of prediction performance for many video clip lengths. Comparing the proposed method with ASs1-3, it can be concluded that the use of player relations reflecting more detailed match conditions, and the inclusion of players’ distances and team information in the graph representation of the proposed architecture, is effective in predicting shoot event occurrence. Comparing the proposed method with AS4, it can be concluded that the inclusion of off-field visual information is also effective in predicting shoot event occurrence. The comparison of the proposed method with CM1 indicates that the proposed method is more effective in predicting shoot event occurrence than CM1, a video-based event prediction method that achieved high performance in a different task. Here, CM1 achieves high prediction accuracy, although not as high as the proposed method. CM1 adopts DSA-RNN, which applies dynamic soft attention to the features by fused frame- and object-level features of each video frame and obtaining the weighted attention level for each object. Therefore, it is suggested that CM1 was able to accurately predict shoot events by appropriately considering the relationships between players in the video. The comparison between the proposed method and CM2 indicates that the proposed method achieves higher classification accuracy than the state-of-the-art video classification method. Here, the results of CM2 were not as good as those of all other methods, and there may be two possible reasons for this problem. The first reason may be the insufficient scale of the dataset. ViViT used in CM2 is a transformer-based method, and it was reported that its generalization performance was low in the case of a small dataset due to weak inductive bias, and a large dataset was required to achieve a high-performance model [56,57]. In this experiment, the training dataset consisted of only 320 video clips, which may have contributed to the suboptimal classification results. The other reason is the potential difficulty of predicting the occurrence of the shooting event based solely on full-frame visual information. One study has reported that video-based methods alone are not as effective as location-based methods in prediction sports match events [28]. Moreover, recent research on event classification tasks in soccer matches, such as shoot, pass, and corner kicks, have demonstrated high performance by jointly leveraging not only full-frame visual information but also audio and players’ visual information [58]. Thus, it can be concluded that CM2 did not achieve good classification performance.
In the experiment, since Ours, ASs1-4, and CM1 predict the probability of a shoot event occurrence at each time step, we also evaluated the prediction accuracy at the time step T o when the shooting event occurs and the prediction earliness for positive samples. These experimental results are shown in Table 2. As shown in Table 2, although the proposed method exhibits the best performance in terms of AP and F1-score, it has been observed that ASs2-3 predicted positive samples correctly at an earlier stage compared to the proposed method. It should be noted that if the model predicts that most video clips are positive at an early stage, the resulting AP may not be accurately reflect the model’s performance due to overfitting on positive samples, despite obtaining a high mTTE. Predicting shoot events at an early stage may support players in making timely decisions based on the outcome. However, obtaining a high TTA without ensuring high accuracy does not hold much significance. In our model, we use the loss function L EXP in Equation (9), where the exponential increase in loss with time for positive samples is intended to strike a balance between the earliness and accuracy of event prediction during training. Table 1 shows that the proposed method and AS1 exhibit stable AP as the video length increases, whereas ASs2-3 shows a decreasing trend in AP. Therefore, our model confirms that methods that consider the distance of players (that is, Ours and AS1) can predict the occurrence of shoot events more meticulously and precisely compared to ASs2-3. From our experiments, we indicate the effectiveness of the proposed method in the perspectives of the ablation studies and the comparison with CMs1-2. Besides, although there are some other event prediction methods, e.g., Graph and Spatio-temporal Continuity based framework (GSC) [59] that is one of the state-of-the art methods, the two contributions mentioned in Section 1 have been well demonstrated by comparing the proposed method with CMs1-2. Thus, we consider the comparisons for the other event prediction methods in the future work.

4.2.2. Qualitative Results

In this subsection, we present the visualization of the shoot event prediction results obtained by the proposed method (refer to Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10) and provide a qualitative assessment of the prediction performance. Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 depict the prediction examples of the proposed method. Specifically, Figure 6 and Figure 7 show successful predictions for positive samples, while Figure 8 is an example of successful negative sample predictions. Figure 9 and Figure 10 are false positive and false negative examples, respectively. To depict each outcome, a threshold value of 0.5 was established to determine whether the model predicted the occurrence of a shooting event in the future. Each example displays several frames of the video clip at the top of the figure. The bottom section of the figure portrays the probability of the shooting prediction (blue curve) and the corresponding aleatoric and epistemic uncertainties (orange and yellow regions, respectively) in the time series. Figure 6 is a brief counter-attack scene that may result in shooting events, and the probability of shooting prediction has reached the threshold value of 3.54 s before the time of the shooting event. In the case of Figure 7, a long forward pass was made, and the predicted probability rapidly increases around the successful pass, resulting in a TTA of 0.89 s. In the example of Figure 6, the model provides predictions with uncertainties during the increasing predicted probability phase. However, in the later phase, the model is relatively certain about its predictions, resulting in low uncertainties. In example of Figure 7, the uncertainty region demonstrates that the model provides predictions with uncertainties during the probability increase phase (when it is not yet known whether the long pass will be successful or not), in comparison to both the initial and later phases. In the case of Figure 8, which is an instance of a true negative depicting a centering scene, the probability does not attain the threshold value. Nonetheless, it suggests that there are significant uncertainties in the predicted outcomes. In the case of Figure 9, which is an instance of false positives and also a centering scene that does not culminate in a shooting event, the probability surpasses the threshold value with considerable uncertainty. In the case of Figure 9, an example of false positive, which is also a centering scene that does not lead to a shoot, the probability exceeds the threshold value with a large uncertainty. In soccer matches, it has been reported that the probability of a centering attack leading to a shoot is lower than that of a middle attack [60]. Therefore, in cases of uncertainty and false prediction results, this fact can be taken into consideration. In soccer matches, the probability that a centering attack leads to a shoot is lower than a middle attack [60], and thus it can be considered that it is confirmed in the case of uncertainty and false prediction results. In the case of Figure 10, which is a false negative example, several defensive players were positioned in front of the attacking players around the goal. It should be noted that the limitation of this study is that the ball carrier is unknown, thus the probability of shoot event occurrence may not have reached the threshold value, even though it occurred.
The proposed method demonstrates effectiveness in predicting the occurrence of shoot events, as evidenced by both qualitative and quantitative evaluations.

5. Conclusions

This paper proposes a shoot event prediction method that considers the spatial-temporal relationships between players in soccer videos using complete bipartite graphs. Specifically, the proposed method represents player relations as a complete bipartite graph using various visual features at each time step. The method then applies GCRNN to calculate latent features that consider the spatio-temporal relations between players. In the proposed method, the use of complete bipartite graphs enables the consideration of both players’ distance and team information, resulting in a more detailed representation of the match situation through the incorporation of player relations. The latent features obtained through GCRNN are then fed into the BNN, which uses them to make accurate predictions of the probability of a shooting event occurring. The BNN takes into account prediction uncertainty to produce reliable estimates. The experimental results using real soccer video clips demonstrate the effectiveness of the proposed method. Our shoot event prediction model is constructed for soccer videos. In soccer, players’ positions and their temporal changes are critical factors for scoring. Besides, to introduce shoot event prediction techniques into tactical coaching in soccer, the reliability of the prediction results should be uncovered. From these standpoints, the contributions of this paper are summarized as follows:
(1)
We developed a new method for predicting shoot events in soccer videos, incorporating player-specific features and spatio-temporal relations with GCRNN.
(2)
The proposed method provides a measure of prediction uncertainty by using BNN, which is intended to enhance the robustness of the predictions.
On the other hand, our shoot event prediction model still has the following limitations:
(1)
Our model cannot consider the skill levels of each player that could change the development of the game. Annotations for each player are needed for introducing such skill levels, but manual annotation is not a realistic solution since it takes a lot of effort and time.
(2)
Reflection of the ball position, which is the direct method for grasping the situations of the game, should be performed in our model. However, when using the scouting soccer videos for recognizing more players, the ball region is too small to be recognized by the object detector. The use of GPS Information is one possible solution, but it is difficult to implement in minor games.
As the applicability, we consider that our event prediction model can easily be applied into other sports. Here, soccer is the specific sport with the property that there is a large number of players and has a high degree of freedom for each player, and the event prediction task in soccer is more difficult than that of other sequentially sports, e.g., basketball, water polo, and rugby. Hence, our prediction model is likely to provide accurate event predictions for such sports, as it has been successful in soccer by capturing the more complex relation between players. Therefore, our model has the high applicability and great potential to provide well benefit for other sports.

Author Contributions

Conceptualization, R.G., Y.M., K.M., T.O. and M.H.; methodology, R.G., Y.M., K.M. and T.O.; software, R.G.; validation, R.G.; data curation, R.G.; writing—original draft preparation, R.G.; writing—review and editing, Y.M., K.M., T.O. and M.H.; visualization, R.G.; funding acquisition, K.M., T.O. and M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by JSPS KAKENHI Grant Numbers JP20K19856, JP21H03456.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Experimental data cannot be disclosed.

Acknowledgments

In this research, we used the data provided by Consadole Sapporo.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Assunção, R.; Pelechrinis, K. Sports analytics in the era of big data: Moving toward the next frontier. Big Data 2018, 6, 237–238. [Google Scholar] [CrossRef]
  2. Mercier, M.A.; Tremblay, M.; Daneau, C.; Descarreaux, M. Individual factors associated with baseball pitching performance: Scoping review. BMJ Open Sport Exerc. Med. 2020, 6, e000704. [Google Scholar] [CrossRef]
  3. Mizels, J.; Erickson, B.; Chalmers, P. Current state of data and analytics research in baseball. Curr. Rev. Musculoskelet. Med. 2022, 15, 283–290. [Google Scholar] [CrossRef] [PubMed]
  4. Sarlis, V.; Tjortjis, C. Sports analytics—Evaluation of basketball players and team performance. Inf. Syst. 2020, 93, 101562. [Google Scholar] [CrossRef]
  5. Cervone, D.; D’Amour, A.; Bornn, L.; Goldsberry, K. A multiresolution stochastic process model for predicting basketball possession outcomes. J. Am. Stat. Assoc. 2016, 111, 585–599. [Google Scholar] [CrossRef]
  6. Kubayi, A.; Larkin, P. Analysis of teams’ corner kicks defensive strategies at the FIFA World Cup 2018. Int. J. Perform. Anal. Sport 2019, 19, 809–819. [Google Scholar] [CrossRef]
  7. Beavan, A.; Fransen, J. The influence of coaching instructions on decision-making in soccer. J. Sport Exerc. Sci. 2021, 5, 3–12. [Google Scholar]
  8. Siddharth, S.; Saurav, S.; Kan, J.; Bimlesh, W.; Song, D.J. Model Driven Inputs to aid Athlete’s Decision Making. In Proceedings of the Asia-Pacific Software Engineering Conference, Singapore, 1–4 December 2020; pp. 485–489. [Google Scholar]
  9. Van Roy, M.; Robberechts, P.; Yang, W.C.; De Raedt, L.; Davis, J. Leaving goals on the pitch: Evaluating decision making in soccer. In Proceedings of the MIT Sloan Sports Analytics Conference, Virtual, 8–9 April 2021; pp. 1–25. [Google Scholar]
  10. Baca, A.; Dabnichki, P.; Hu, C.W.; Kornfeind, P.; Exel, J. Ubiquitous Computing in Sports and Physical Activity—Recent Trends and Developments. Sensors 2022, 22, 8370. [Google Scholar] [CrossRef]
  11. Laws of the Game 2018/2019. Available online: https://digitalhub.fifa.com/m/50518593a0941079/original/khhloe2xoigyna8juxw3-pdf (accessed on 11 October 2022).
  12. Anzer, G.; Bauer, P. A goal scoring probability model for shots based on synchronized positional and event data in football (soccer). Front. Sports Act. Living 2021, 3, 624475. [Google Scholar] [CrossRef] [PubMed]
  13. Link, D.; Lang, S.; Seidenschwarz, P. Real time quantification of dangerousity in football using spatiotemporal tracking data. PLoS ONE 2016, 11, e0168768. [Google Scholar] [CrossRef]
  14. Decroos, T.; Van Haaren, J.; Dzyuba, V.; Davis, J. STARSS: A spatio-temporal action rating system for soccer. In Proceedings of the ECML/PKDD Workshop on Machine Learning and Data Mining for Sports Analytics, Skopje, North Macedonia, 18 September 2017; Volume 1971, pp. 11–20. [Google Scholar]
  15. Spearman, W.; Basye, A.; Dick, G.; Hotovy, R.; Pop, P. Physics-based modeling of pass probabilities in soccer. In Proceedings of the MIT Sloan Sports Analytics Conference, Boston, MA, USA, 3–4 March 2017; pp. 1–14. [Google Scholar]
  16. Spearman, W. Beyond expected goals. In Proceedings of the MIT Sloan Sports Analytics Conference, Boston, MA, USA, 23–24 February 2018; pp. 1–17. [Google Scholar]
  17. Goldner, K. A Markov model of football: Using stochastic processes to model a football drive. J. Quant. Anal. Sports 2012, 8, 1–18. [Google Scholar] [CrossRef]
  18. Lucey, P.; Bialkowski, A.; Monfort, M.; Carr, P.; Matthews, I. Quality vs quantity: Improved shot prediction in soccer using strategic features from spatiotemporal data. In Proceedings of the Annual MIT Sloan Sports Analytics Conference, Boston, MA, USA, 28 February–1 March 2014; pp. 1–9. [Google Scholar]
  19. Dick, U.; Brefeld, U. Learning to rate player positioning in soccer. Big Data 2019, 7, 71–82. [Google Scholar] [CrossRef] [PubMed]
  20. Liu, G.; Luo, Y.; Schulte, O.; Kharrat, T. Deep soccer analytics: Learning an action-value function for evaluating soccer players. Data Min. Knowl. Discov. 2020, 34, 1531–1559. [Google Scholar] [CrossRef]
  21. Szczepański, Ł.; McHale, I. Beyond completion rate: Evaluating the passing ability of footballers. J. R. Stat. Soc. Ser. Stat. Soc. 2016, 179, 513–533. [Google Scholar] [CrossRef]
  22. Power, P.; Ruiz, H.; Wei, X.; Lucey, P. Not all passes are created equal: Objectively measuring the risk and reward of passes in soccer from tracking data. In Proceedings of the ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, Halifax, NS, Canada, 13–17 August 2017; pp. 1605–1613. [Google Scholar]
  23. Decroos, T.; Davis, J. Interpretable prediction of goals in soccer. In Proceedings of the AAAI Workshop on Artificial Intelligence in Team Sports, Honolulu, HI, USA, 27–28 January 2019. [Google Scholar]
  24. Fernández, J.; Bornn, L.; Cervone, D. Decomposing the immeasurable sport: A deep learning expected possession value framework for soccer. In Proceedings of the MIT Sloan Sports Analytics Conference, Boston, MA, USA, 1–2 March 2019. [Google Scholar]
  25. Fernández, J.; Bornn, L. PassNet: Learning Pass Probability Surfaces from Single-Location Labels. An Architecture for Visually-Interpretable Soccer Analytics. In Proceedings of the International Conference on Learning Representations, New Orleans, LA, USA, 6–9 May 2019; Available online: https://openreview.net/forum?id=r1xxKJBKvr (accessed on 2 April 2023).
  26. Fernández, J.; Bornn, L. Soccermap: A deep learning architecture for visually-interpretable analysis in soccer. In Proceedings of the Joint European Conference on Machine Learning and Knowledge Discovery in Databases, Ghent, Belgium, 14–18 September 2020; pp. 491–506. [Google Scholar]
  27. Simpson, I.; Beal, R.J.; Locke, D.; Norman, T.J. Seq2Event: Learning the Language of Soccer using Transformer-based Match Event Prediction. In Proceedings of the ACM SIGKDD Conference on Knowledge Discovery and Data Mining, Washington, DC, USA, 14–18 August 2022; pp. 3898–3908. [Google Scholar]
  28. Felsen, P.; Agrawal, P.; Malik, J. What will happen next? Forecasting player moves in sports videos. In Proceedings of the IEEE International Conference on Computer Vision, Venice, Italy, 22–29 October 2017; pp. 3342–3351. [Google Scholar]
  29. Ruiz, L.; Gama, F.; Ribeiro, A. Gated graph convolutional recurrent neural networks. In Proceedings of the European Signal Processing Conference, A Coruna, Spain, 2–6 September 2019; pp. 1–5. [Google Scholar]
  30. Bao, W.; Yu, Q.; Kong, Y. Uncertainty-based traffic accident anticipation with spatio-temporal relational learning. In Proceedings of the ACM International Conference on Multimedia, Seattle, WA, USA, 12–16 October 2020; pp. 2682–2690. [Google Scholar]
  31. Defferrard, M.; Bresson, X.; Vandergheynst, P. Convolutional neural networks on graphs with fast localized spectral filtering. Adv. Neural Inf. Process. Syst. 2016, 29, 3844–3852. [Google Scholar]
  32. Kipf, T.N.; Welling, M. Semi-supervised classification with graph convolutional networks. arXiv 2016, arXiv:1609.02907. [Google Scholar]
  33. Lipton, Z.C.; Berkowitz, J.; Elkan, C. A critical review of recurrent neural networks for sequence learning. arXiv 2015, arXiv:1506.00019. [Google Scholar]
  34. Seo, Y.; Defferrard, M.; Vandergheynst, P.; Bresson, X. Structured sequence modeling with graph convolutional recurrent networks. In Proceedings of the International Conference on Neural Information Processing, Siem Reap, Cambodia, 13–16 December 2018; pp. 362–373. [Google Scholar]
  35. Neal, R.M. Bayesian Learning for Neural Networks; Springer Science & Business Media: Cham, Switzerland, 2012. [Google Scholar]
  36. Blundell, C.; Cornebise, J.; Kavukcuoglu, K.; Wierstra, D. Weight uncertainty in neural network. In Proceedings of the International Conference on Machine Learning, Lille, France, 6–11 July 2015; pp. 1613–1622. [Google Scholar]
  37. Goka, R.; Moroto, Y.; Maeda, K.; Ogawa, T.; Haseyama, M. Shoot event prediction from soccer videos by considering players’ spatio-temporal relations. In Proceedings of the IEEE Global Conference on Consumer Electronics, Osaka, Japan, 18–21 October 2022; pp. 193–194. [Google Scholar]
  38. Epshtein, B.; Ofek, E.; Wexler, Y. Detecting text in natural scenes with stroke width transform. In Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, San Francisco, CA, USA, 13–18 June 2010; pp. 2963–2970. [Google Scholar]
  39. Sharma, R.A.; Bhat, B.; Gandhi, V.; Jawahar, C. Automated top view registration of broadcast football videos. In Proceedings of the IEEE Winter Conference on Applications of Computer Vision, Lake Tahoe, NV, USA, 12–15 March 2018; pp. 305–313. [Google Scholar]
  40. He, K.; Gkioxari, G.; Dollár, P.; Girshick, R. Mask r-cnn. In Proceedings of the IEEE International Conference on Computer Vision, Venice, Italy, 22–29 October 2017; pp. 2961–2969. [Google Scholar]
  41. Chen, J.; Little, J.J. Sports camera calibration via synthetic data. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition Workshops, Long Beach, CA, USA, 16–17 June 2019; pp. 2497–2504. [Google Scholar]
  42. Isola, P.; Zhu, J.Y.; Zhou, T.; Efros, A.A. Image-to-image translation with conditional adversarial networks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, HI, USA, 21–26 July 2017; pp. 1125–1134. [Google Scholar]
  43. Ronneberger, O.; Fischer, P.; Brox, T. U-net: Convolutional networks for biomedical image segmentation. In Proceedings of the International Conference on Medical Image Computing and Computer-Assisted Intervention, Munich, Germany, 5–9 October 2015; pp. 234–241. [Google Scholar]
  44. Dalal, N.; Triggs, B. Histograms of oriented gradients for human detection. In Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, San Diego, CA, USA, 20–26 June 2005; Volume 1, pp. 886–893. [Google Scholar]
  45. Ren, S.; He, K.; Girshick, R.; Sun, J. Faster r-cnn: Towards real-time object detection with region proposal networks. Adv. Neural Inf. Process. Syst. 2015, 28, 91–99. [Google Scholar] [CrossRef] [PubMed]
  46. Sun, S.B.; Cui, R.y. Player classification algorithm based on digraph in soccer video. In Proceedings of the IEEE Joint International Information Technology and Artificial Intelligence Conference, Chongqing, China, 20–21 December 2014; pp. 459–463. [Google Scholar]
  47. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep residual learning for image recognition. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 27–30 June 2016; pp. 770–778. [Google Scholar]
  48. Deng, J.; Dong, W.; Socher, R.; Li, L.J.; Li, K.; Fei-Fei, L. Imagenet: A large-scale hierarchical image database. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Miami, FL, USA, 20–25 June 2009; pp. 248–255. [Google Scholar]
  49. González-Víllora, S.; García-López, L.M.; Gutiérrez-Díaz, D.; Pastor-Vicedo, J.C. Tactical awareness, decision making and skill in youth soccer players (under-14 years). J. Hum. Sport Exerc. 2013, 8, 412–426. [Google Scholar] [CrossRef]
  50. Owen, A.B. Monte Carlo theory, methods and examples. In Handbooks in Operations Research and Management Science; Elsevier: Amsterdam, The Netherlands, 2013. [Google Scholar]
  51. Kendall, A.; Gal, Y. What uncertainties do we need in bayesian deep learning for computer vision? Adv. Neural Inf. Process. Syst. 2017, 30, 5574–5584. [Google Scholar]
  52. Diederik, P.K.; Ba, J.L. Adam: A method for Stochastic Optimization. In Proceedings of the International Conference on Learning Representations, San Diego, CA, USA, 7–9 May 2015; pp. 1–15. [Google Scholar]
  53. Chan, F.H.; Chen, Y.T.; Xiang, Y.; Sun, M. Anticipating accidents in dashcam videos. In Proceedings of the Asian Conference on Computer Vision, Taipei, Taiwan, 20–24 November 2016; pp. 136–153. [Google Scholar]
  54. Xu, K.; Ba, J.; Kiros, R.; Cho, K.; Courville, A.; Salakhudinov, R.; Zemel, R.; Bengio, Y. Show, attend and tell: Neural image caption generation with visual attention. In Proceedings of the International Conference on Machine Learning, Lille, France, 6–11 July 2015; pp. 2048–2057. [Google Scholar]
  55. Arnab, A.; Dehghani, M.; Heigold, G.; Sun, C.; Lučić, M.; Schmid, C. Vivit: A video vision transformer. In Proceedings of the IEEE/CVF International Conference on Computer Vision, Montreal, BC, Canada, 11–17 October 2021; pp. 6836–6846. [Google Scholar]
  56. Dosovitskiy, A.; Beyer, L.; Kolesnikov, A.; Weissenborn, D.; Zhai, X.; Unterthiner, T.; Dehghani, M.; Minderer, M.; Heigold, G.; Gelly, S.; et al. An image is worth 16x16 words: Transformers for image recognition at scale. arXiv 2020, arXiv:2010.11929. [Google Scholar]
  57. Raghu, M.; Unterthiner, T.; Kornblith, S.; Zhang, C.; Dosovitskiy, A. Do vision transformers see like convolutional neural networks? Adv. Neural Inf. Process. Syst. 2021, 34, 12116–12128. [Google Scholar]
  58. Haruyama, T.; Takahashi, S.; Ogawa, T.; Haseyama, M. User-selectable event summarization in unedited raw soccer video via multimodal bidirectional LSTM. ITE Trans. Media Technol. Appl. 2021, 9, 42–53. [Google Scholar] [CrossRef]
  59. Wang, T.; Chen, K.; Chen, G.; Li, B.; Li, Z.; Liu, Z.; Jiang, C. GSC: A Graph and Spatio-temporal Continuity Based Framework for Accident Anticipation. IEEE Trans. Intell. Veh. 2023; online early access. [Google Scholar] [CrossRef]
  60. Izquierdo, J.M.; Redondo, J.C. Offensive difference styles and technical situational variables between european and south american elite football leagues. MHSalud 2022, 19, 25–37. [Google Scholar] [CrossRef]
Figure 1. This section provides an overview of our pre-processing approach. Firstly, we use a full-frame image (a) as input and obtain a Stroke Width Transform (SWT) image (b) and a field mask (c) through SWT and a pix2pix-cGan segmentation generator network, respectively. This leads to an edge image (d) of the field. Next, we capture the corresponding camera parameters by searching for the nearest neighbor field edge image (e) based on HOG features. In another step, we classify the detected players into their respective teams based on HSV color. Finally, using both pieces of information, we obtain the players’ positions on the field as well as their teams.
Figure 1. This section provides an overview of our pre-processing approach. Firstly, we use a full-frame image (a) as input and obtain a Stroke Width Transform (SWT) image (b) and a field mask (c) through SWT and a pix2pix-cGan segmentation generator network, respectively. This leads to an edge image (d) of the field. Next, we capture the corresponding camera parameters by searching for the nearest neighbor field edge image (e) based on HOG features. In another step, we classify the detected players into their respective teams based on HSV color. Finally, using both pieces of information, we obtain the players’ positions on the field as well as their teams.
Sensors 23 04506 g001
Figure 2. Perspective projection camera model. First, the coordinates of the field edge ( x w , y w , z w ) in the world coordinate system are transformed to their coordinates ( x c , y c , z c ) in the camera coordinate system, which originates from the camera center. Then, we obtain the projected coordinates ( x n , y n ) in the normalized image plane of size u × v .
Figure 2. Perspective projection camera model. First, the coordinates of the field edge ( x w , y w , z w ) in the world coordinate system are transformed to their coordinates ( x c , y c , z c ) in the camera coordinate system, which originates from the camera center. Then, we obtain the projected coordinates ( x n , y n ) in the normalized image plane of size u × v .
Sensors 23 04506 g002
Figure 3. Standard soccer field dimensions. These are expressed in meters for each measurement.
Figure 3. Standard soccer field dimensions. These are expressed in meters for each measurement.
Sensors 23 04506 g003
Figure 4. Overview of the proposed method. Our method starts by extracting visual features of each detected player and full-frame at each time step form a given soccer video input. A complete bipartite graph is constructed with these visual features as node features. We then use GCRNN to calculate latent features. Next, our method predicts the probability of shoot event occurrence based on BNN, and we obtain the resulting predictions along with their corresponding prediction uncertainties. The overall process is illustrated in Figure.
Figure 4. Overview of the proposed method. Our method starts by extracting visual features of each detected player and full-frame at each time step form a given soccer video input. A complete bipartite graph is constructed with these visual features as node features. We then use GCRNN to calculate latent features. Next, our method predicts the probability of shoot event occurrence based on BNN, and we obtain the resulting predictions along with their corresponding prediction uncertainties. The overall process is illustrated in Figure.
Sensors 23 04506 g004
Figure 5. The GCRNN architecture comprises stacked GCN layers that initially learn the spatial relations. By inputting the RNN’s hidden state h t 1 into the second GCN layer, the method can calculate the latent feature Z t , considering both temporal and spatial relations. This process in repeated recursively to capture latent features at each time step, as illustrated in Figure.
Figure 5. The GCRNN architecture comprises stacked GCN layers that initially learn the spatial relations. By inputting the RNN’s hidden state h t 1 into the second GCN layer, the method can calculate the latent feature Z t , considering both temporal and spatial relations. This process in repeated recursively to capture latent features at each time step, as illustrated in Figure.
Sensors 23 04506 g005
Figure 6. The qualitative prediction results of the proposed method for shoot event occurrence at a threshold of s = 0.5 . It displays an example of a positive sample with successful predictions. The proposed method predicts the probability of shoot event occurrence at each time step (indicated by the blue curve) and estimates both aleatoric uncertainty (shown in the orange region) and epistemic uncertainty (shown in the yellow region) for each time step.
Figure 6. The qualitative prediction results of the proposed method for shoot event occurrence at a threshold of s = 0.5 . It displays an example of a positive sample with successful predictions. The proposed method predicts the probability of shoot event occurrence at each time step (indicated by the blue curve) and estimates both aleatoric uncertainty (shown in the orange region) and epistemic uncertainty (shown in the yellow region) for each time step.
Sensors 23 04506 g006
Figure 7. This figure illustrates an instance of accurate predictions for the positive sample. The details depicted are identical to those in Figure 6.
Figure 7. This figure illustrates an instance of accurate predictions for the positive sample. The details depicted are identical to those in Figure 6.
Sensors 23 04506 g007
Figure 8. This figure presents an example of precise predictions for the negative sample. The specifics depicted are identical to those in Figure 6.
Figure 8. This figure presents an example of precise predictions for the negative sample. The specifics depicted are identical to those in Figure 6.
Sensors 23 04506 g008
Figure 9. This figure demonstrates an instance of inaccurate predictions for the negative sample. The particulars shown are identical to those in Figure 6.
Figure 9. This figure demonstrates an instance of inaccurate predictions for the negative sample. The particulars shown are identical to those in Figure 6.
Sensors 23 04506 g009
Figure 10. This figure displays an instance of erroneous predictions for the positive sample. The specifics depicted are identical to those in Figure 6.
Figure 10. This figure displays an instance of erroneous predictions for the positive sample. The specifics depicted are identical to those in Figure 6.
Sensors 23 04506 g010
Table 1. Comparison of AP and F1-score among different video clip lengths for the proposed and comparison methods. Each video clip length is counted from the start of the video clip. The best performance at each video clip length is shown in bold.
Table 1. Comparison of AP and F1-score among different video clip lengths for the proposed and comparison methods. Each video clip length is counted from the start of the video clip. The best performance at each video clip length is shown in bold.
   t = 10      t = 20      t = 30      t = 40      t = 50      t = 60   
APF1-ScoreAPF1-ScoreAPF1-ScoreAPF1-ScoreAPF1-ScoreAPF1-Score
Ours0.9600.8800.9660.9090.9640.9230.9640.9230.9700.9210.9630.914
AS1 [37]0.9550.8890.9610.8970.9620.8970.9590.9090.9580.9000.9540.880
AS20.9540.8920.9610.8810.9320.8860.9290.8920.9240.8920.9190.892
AS30.9560.8650.9600.8890.9340.8830.9140.8950.9090.8950.9040.895
AS40.9510.8810.9640.8890.9610.9140.9620.9040.9630.8970.9560.897
CM1 [53]0.9310.8780.9580.9020.9660.9020.9590.9000.9610.8920.9580.881
CM2 [55]0.7520.7140.7760.7290.7780.7250.7820.7470.7840.7590.7900.722
Table 2. AP, F1-score and mTTE at the time step T o . Please, note that CM2 cannot handle varying video lengths, so it is denotes as ‘‘-’’. The best performance in terms of prediction accuracy is denoted in bold.
Table 2. AP, F1-score and mTTE at the time step T o . Please, note that CM2 cannot handle varying video lengths, so it is denotes as ‘‘-’’. The best performance in terms of prediction accuracy is denoted in bold.
APF1-ScoremTTE(s)
Ours0.9670.9143.88
AS1 [37]0.9410.8803.60
AS20.9170.8924.32
AS30.9030.8954.66
AS40.9500.8863.67
CM1 [53]0.9310.8533.62
CM2 [55]---
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Goka, R.; Moroto, Y.; Maeda, K.; Ogawa, T.; Haseyama, M. Prediction of Shooting Events in Soccer Videos Using Complete Bipartite Graphs and Players’ Spatial-Temporal Relations. Sensors 2023, 23, 4506. https://doi.org/10.3390/s23094506

AMA Style

Goka R, Moroto Y, Maeda K, Ogawa T, Haseyama M. Prediction of Shooting Events in Soccer Videos Using Complete Bipartite Graphs and Players’ Spatial-Temporal Relations. Sensors. 2023; 23(9):4506. https://doi.org/10.3390/s23094506

Chicago/Turabian Style

Goka, Ryota, Yuya Moroto, Keisuke Maeda, Takahiro Ogawa, and Miki Haseyama. 2023. "Prediction of Shooting Events in Soccer Videos Using Complete Bipartite Graphs and Players’ Spatial-Temporal Relations" Sensors 23, no. 9: 4506. https://doi.org/10.3390/s23094506

APA Style

Goka, R., Moroto, Y., Maeda, K., Ogawa, T., & Haseyama, M. (2023). Prediction of Shooting Events in Soccer Videos Using Complete Bipartite Graphs and Players’ Spatial-Temporal Relations. Sensors, 23(9), 4506. https://doi.org/10.3390/s23094506

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop