arxiv: v1 [physics.soc-ph] 19 May 2018

Transcription

1 Classification algorithm for formations in football games Takuma Narizuka 1 and Yoshihiro Yamazaki 2 1 Department of Physics, Faculty of Science and Engineering, Chuo University, Bunkyo, Tokyo , Japan 2 Department of Physics, School of Advanced Science and Engineering, Waseda University, Shinjuku, Tokyo , Japan arxiv: v1 [physics.soc-ph] 19 May 2018 Abstract This paper develops a classification algorithm for formations in team sports, with a focus on football games. Our method first classifies formations into several average formations: 424, 4141, 433, 541, and 343. Then, each average formation is further divided into more specific patterns in which the configurations of players are slightly different. The latter step is based on hierarchical clustering and the Delaunay method, which defines a formation of a team as an adjacency matrix of Delaunay triangulation. A formation classified using our method is expressed in a form such as 424-C1. 1. Introduction In competitive team sports, such as football and basketball, each player coordinates with team members and interacts with opposing players. Throughout such interactions, players maintain a certain formation at the team level. Such a formation structure reflects a team s strategies for achieving effective attacks and defense in order to win [1 4]. A traditional method of characterizing formations employs notation such as 4-4-2, which indicates four defenders, four midfielders, and two forwards. Although this is a convenient means of roughly grasping formation structures, it is too simple to analyze real games. In fact, the following more quantitative methods have been introduced. The first example is based on a Voronoi region defined for each player, which is the set of field locations whose distances from the player are less than from any other [5]. Intuitively, this corresponds to the territory of the player on the field. The basic properties of the Voronoi region have been investigated for football and hockey games [6, 7], and modified version considering the velocity and acceleration of a player have also been proposed [8 12]. Bialkowski et al. developed another approach to formations, called role representation [13, 14]. Here, the role represents the relative position of each player in the formation, such as center forward or left wing. The key idea behind their method is that players are not distinguished by their uniform numbers, but rather by the role numbers assigned to them at each frame of the game data. Specifically, a heat map for a team is divided into 10 heat maps (roles) corresponding to roles, in order to yield a minimal overlap. Then, the set of rules is regarded as the formation of the team. It has been reported that this method can detect positional exchanges of players, and characteristic formations such as are distinguished more efficiently. address: pararel@gmail.com (T. Narizuka). 1

2 Along with these studies, we have proposed using the Delaunay method, which identifies the adjacency relationships of players Voronoi regions, i.e., the Delaunay network, with the formation of a team [15]. Because the formation at time t is quantified by an adjacency matrix in this method, dissimilarity measures between two different formations can be defined. On the basis of the Delaunay method, we have also proposed a classification algorithm for formations in a single game. This algorithm classifies Delaunay networks into clusters by means of hierarchical clustering. We have demonstrated that our method can characterize the differences and dynamics of football formations at different time scales within a game. The Delaunay method is useful for the time-series analysis and quantitative comparison of formations. However, the above classification algorithm for a single game cannot be straightforwardly extended to the case of multiple games. In this paper, we propose an extended algorithm that can classify formations over multiple games. 2. Method In the following analysis, we employ datasets comprising 45 football games provided by DataStadium Inc., Japan. These are the games from the 13th to 17th sections of the J1 League second stage in Because the J1 League comprises 18 teams, there are five games per team. Each dataset contains all player positions every 0.04 seconds. We focus on the 10 players (N = 10) other than the goalkeeper for each team. For simplicity, we only analyze the data of the first halves of games where player substitutions did not occur. Thus, seven teams were excluded, and our data consists of uniform numbers of players who are not changed during a game. The position of the j-th player of a team at time t is denoted as r j (t) = [x j (t), y j (t)]. The centroid position and gyration radius of a team respectively defined as follows: r c (t) = 1 N N r j (t), (1) j=1 σ(t) = 1 N N r c (t) r j (t) 2. (2) j=1 Using r c (t) and σ(t), the normalized coordinates R j (t) for the j-th player are calculated as 3. Classification algorithm for a single game R j (t) = r j(t) r c (t). (3) σ(t) Here, we present a brief review of the Delaunay method and the previous classification algorithm we introduced in [15]. As shown in Fig. 1(a), the Delaunay method quantifies a formation of a team using the adjacency matrix A(t) of the Delaunay network, whose components A ij (t) are given by 1 if the Voronoi regions of players i and j are adjacent with each other at t, A ij (t) = 0 otherwise. 2

3 Owing to this quantification, a dissimilarity measure between two formations at different times can be introduced as D tt = N N [A ij (t) A ij (t )] 2. (4) i=1 j=1 D tt depends only on the adjacency relationship of players, and is independent of the deviation σ(t) of a team. Based on this dissimilarity measure, we introduced a classification algorithm for formations appearing in a single game through the following four steps (i)-(iv). (i) The Delaunay networks of all frames in a single game are computed. (ii) Hierarchical clustering is performed using Ward s method [16], by employing the dissimilarity measure defined by eq. (4). (iii) N c clusters are extracted by cutting the dendrogram at a certain height h c. (iv) Coarse-grained formations are visualized from each cluster as follows. For each Delaunay network in a cluster, the positional coordinates of each player are converted into normalized coordinates using eq. (3). Then, the average positions of players are visualized using ellipses. For example, we show the coarsegrained formations with N c = 3 obtained from the 13th section of Sendai in Fig. 1(b). Each cluster is distinguished by a cluster number, and the difference between them is that several pairs of players exchange their positions, for example, players 2 and 3, and players 5 and Extension of the previous classification algorithm to multiple games Let us consider the problem of classifying Delaunay networks over multiple games into clusters. In order to quantify a formation using an adjacency matrix A(t), uniform numbers U = [a, b,..., j] of players need to be assigned to the indexes I = [1, 2,..., 10] of A(t). If we classify Delaunay networks of a single game in which no player substitutions occur, then an arbitrary correspondence between U and I can be employed. However, classification over multiple games requires the assignment of multiple uniform numbers U 1, U 2,... for different games to one set of indices I. Because such an assignment is not uniquely determined, a criterion must be defined. Here, we adopt the framework of role representation introduced in [13, 14]. That is, we assume that players with similar average positions in a formation play the same roles in a team. Then, we label each player by a role number and identify them with the indexes I of A(t). In the following, we propose an extended classification algorithm based on this idea, consisting of two parts I and II. Part I: classification into average formation In part I, we assign the same index i of A(t) to players whose average positions in a formation are approximately the same. To estimate the average position of each player in a game, we compute the heat map of each game for each team in the normalized coordinates (3). We present the heat maps obtained for all teams and games in Fig. 2. In this figure, the average positions of players are expressed by ellipses. The direction and magnitude of an ellipse is determined by the eigenvector and eigenvalue of the covariance matrix for the corresponding player s position. These heat maps appear to be classified into several patterns. In fact, we find from our data that they belong to one of the following five patterns: 424, 4141, 433, 541, and 343 (these are referred to as average formations hereafter). A schematic representation of each 3

4 average formation is shown in Fig. 3(a). It should be noted that we manually classified the heat maps into average formations. We also present the average formations by team in Fig. 3(b). Although most teams used the same average formation throughout the five games, several teams changed. For a certain average formation, the ellipses (average positions of players in a game) are distinguished by serial numbers, as shown in Fig. 3(b). It is considered that players belonging to the same average formation with the same serial number play the same role in the team (e.g., player 1 in 4141 is interpreted as a center forward ). Therefore, we identify these serial numbers with the indexes I of A(t), and a one-to-one correspondence between U 1, U 2,... and I is obtained for each average formation. Part II: classification into clusters As shown in Fig. 2, the ellipses of some players in a heat map overlap, indicating that these players exchange their positions or move close to each other in the game. In addition, the configurations of players are slightly different even within the same average formation. In order to distinguish such patternsin part II we classify all the Delaunay networks belonging to the same average formation into clusters using the previous classification algorithm (see Sec. 3 for details). Figures 4 and 5 present typical examples of clustering results, with N c = 4 and N c = 8, respectively, for the five games of Sendai. Because Sendai adopted 424 in all five games (see Fig.3(b)), the coarse-grained formation obtained using this method is expressed as 424-C1, where the former number denotes the average formation and the latter is the cluster number. Furthermore, each ellipse in a cluster in Fig. 5(b) consists of all the positions of players with the same index in the five teams. We find that each cluster exhibits a more specific pattern compared with the corresponding average formations. It should be noted that C4 in N c = 4 and C6 in N c = 8 only include irregular patterns, which could be associated with transitional situations such as competition in front of goal or counter attacks. Here, the value of N c depends on the cutting height h c of the dendrogram in Fig 4(b), where the height represents the distance between two merged clusters in the clustering process. By varying N c, we can control the degree of coarse graining of formations: more fine (coarse) patterns are obtained by increasing (decreasing) N c (see Fig 5 for a typical example). It should be noted that N c = 1 is the most coarse pattern, corresponding to the superposition of all average formations of the five games. 5. Discussion We have proposed an extended classification algorithm based on role representation (part I) and hierarchical clustering (part II). Here, we compare our classification algorithm with the method introduced by Bialkowski et al. [13, 14]. In that method, a 2D probability distribution H( r) (heat map) for a team is divided into 10 heat maps, H( r) = 10 r=1 H r( r), and the set F = {H r ( r); r = 1,, 10} is regarded as the formation. Each H r ( r) is computed to achieve a minimal overlap with others, under the condition that each player belongs to a different r at each frame. Because each player is labeled by a role number r instead of uniform number u at each frame, this method is called role representation. In the role representation approach, H r ( r) consists of various players at different frames, and patterns in which two players exchange 4

5 their positions are regarded as the same. In contrast, our algorithm describes an entire heat map H( r) as the sum of players heat maps, H( r) = 10 u=1 H u( r), where u denotes the uniform number. The set F = {H u ( r); u = 1,, 10} is called a average formation. This decomposition does not achieve the minimal overlap, namely, players with different u values can exchange their positions during a game. Instead, our method distinguishes such position-exchanged patterns as different formations, based on the Delaunay method and hierarchical clustering. In this sense, our method realizes more detailed characterization of formations compared with that in [13, 14]. While our decomposition of the entire heat map H( r) does not achieve the minimal overlap, the average positions of players, expressed by ellipses, are still clearly separated (see Fig. 2). That is, each player carries out an individual role in a football game. This feature of football games allows us to label players not only by uniform numbers U, but also by role numbers (role representation). Furthermore, it provides a criterion for the correspondence between multiple uniform numbers U 1, U 2,... and the indexes I of A(t), and allows hierarchical clustering to be realized over multiple games. We note that our method can be applied to specific sports in which players average positions are almost fixed, because it relies on the one-to-one correspondence between U and I. The variation in average formations and switches among them are a reflection of teams strategies [1, 2]. It has been reported that football teams adopt a so-called win-stay lose-shift strategy for formation changes between games [2]: they tend to adopt the same (a different) formation after a win (loss). Our method has the potential to provide a more detailed characterization of strategies or game flow by focusing on formation changes within a game. As a real-world example, we present the transition diagram between the clusters of Fig. 5(b) in Fig. 6. Here, a transition from one cluster to another represents a change in the configuration of the players in the formation, e.g., C1 C2 indicates that players 2 and 3 exchange their positions. We find that transition diagrams of five games exhibit different transition patterns, even though they belong to the same average formation of 424. We expect the analysis of the cluster transitions for different N c values to provide insights into the characterization of team styles. Regarding this type of analysis, a further extension of the Delaunay network could also be considered. In fact, the present Delaunay network lacks information on opposing teams. This means that the edges do not always represent pass routes, because opposing players may be present on these edges. We can address this problem by introducing a Delaunay triangulation method including an opposing team. Further dynamical analyses of formation structures incorporating ball passes or interactions with opposing players by employing this extended Delaunay network will be a topic of future research. In conclusion, we have developed a framework for clustering formations in team sports. A formation classified using our method is expressed in a form such as 424-C1, indicating that a football formation can be understood as the combination of a average formation and the deviation from it. We expect that our method can be applied to various team sports and provide a common tool for formation analysis. 5

6 Acknowledgements The authors are very grateful to DataStadium Inc., Japan for providing the player tracking data. The authors thank H. Kuninaka for a fruitful discussion. This work was partially supported by the Data Centric Science Research Commons Project of the Research Organization of Information and Systems, Japan, a Grant-in-Aid for Young Scientists 18K18013 from the Japan Society for the Promotion of Science (JSPS), and Hayao Nakayama Foundation for Science, Technology and Culture. References [1] Hirotsu, N., Ito, M., Miyaji, C., Hamano, K., Taguchi, A.: Modeling tactical changes of formation in association football as a non-zero-sum game. Journal of Quantitative Analysis in Sports 5 (2009) [2] Tamura, K., Masuda, N.: Win-stay lose-shift strategy in formation changes in football. EPJ Data Science 4, 9 (2015) [3] Memmert, D., Lemmink, K.A., Sampaio, J.: Current approaches to tactical performance analyses in soccer using position data. Sports Medicine 47(1), 1 10 (2017) [4] Sumpter, D.: Soccermatics: Mathematical Adventures in the Beautiful Game. Bloomsbury Sigma, London (2017) [5] Okabe, A., Boots, B., Sugihara, K., Nok-Chiu, S.: Spatial Tessellations: Concepts and Applications of Voronoi Diagrams. John Wiley & Sons, New York (2000) [6] Kim, S.: Voronoi analysis of a soccer game. Nonlinear Analysis: Modelling and Control 9, (2004) [7] Fonseca, S., Milho, J., Travassos, B., Araújo, D.: Spatial dynamics of team sports exposed by Voronoi diagrams. Human Movement Science 31, (2012) [8] Taki, T., Hasegawa, J., Fukumura, T.: Development of motion analysis system for quantitative evaluation of teamwork in soccer games. Proceedings of the 3rd IEEE International Conference on Image Processing 3, (1996) [9] Taki, T., Hasegawa, J.: Visualization of dominant region in team games and its application to teamwork analysis. Proceedings of the Computer Graphics International 2000, (2000) [10] Fujimura, A., Sugihara, K.: Geometric analysis and quantitative evaluation of sport teamwork. Systems and Computers in Japan 36, (2005) [11] Gudmundsson, J., Wolle, T.: Football analysis using spatio-temporal tools. Computers, Environment and Urban Systems 47, (2014) [12] Gudmundsson, J., Horton, M.: Spatio-temporal analysis of team sports. ACM Computing Surveys (CSUR) 50, 22 (2017) 6

7 [13] Bialkowski, A., Lucey, P., Carr, P., Yue, Y., Sridharan, S., Matthews, I.: Large-scale analysis of soccer matches using spatiotemporal tracking data. Proceedings of the 2014 IEEE International Conference on Data Mining, (2014) [14] Bialkowski, A., Lucey, P., Carr, P., Matthews, I., Sridharan, S., Fookes, C.: Discovering team structures in soccer from spatiotemporal data. IEEE Transactions on Knowledge and Data Engineering 28, (2016) [15] (in Japanese) Narizuka, T., Yamazaki, Y.: Characterization of the formation structure in team sports. Proceedings of the Institute of Statistical Mathematics, Special Topic: New Challenges to Statistical Science in Sports 65, (2017). [English version: arxiv: ] [16] Pang-Ning, T., Steinbach, M., Kumar, V.: Introduction to Data Mining. Addison Wesley, Boston (2005) 7

8 <latexit sha1_base64="x0vvtaelfyxek+dgehkmhm6lpnw=">aaaeznichvhbbtnaej0aa8ufmsificuktfv6idyicysevmsfy9qstljdrfzmk67q2ja9sqiwxrvx4qiie0gcej/bhr/gup4acsxslxwyb6yantistd7zmfdm38w4gssisenjgnzop3/h4uily+nylavltzvru5e/dblvm9/1w33hjrgrpn6wqrp8pwi5pxbcvuccp8riezepi+f7t+qo4icdu++jnmc2rfdnrdmenuqgeuamy+u92taisrzef15ih6e9shm2wsngckfknzi4awkxfzalwoamofkf/hc4izgfnpw5+eojxnzblfksuqyu75zfykte+exnuf0fnqvavx2zwwl1ellpcq+bz9awqte/khudwp02qvpk2jbzow75avm172bbf3xgmiqbcpbaou2cdrf+b2achqb9h5cgl0rlqdihfazkdhhfewgj9xj/fbwd5f4p71noslezvulidpfjyi1+o5/okf1kp9mf9hdlrktlylsm8htgxb50ll/e2dn7l2uap4sjv6y5mix04j7sklb7odxzfwzmj5+9pt25v72wrnmp9cfqf09p6beswit/sy9bfpvdhd0oahl3riv3hp4xngxfydpsxz4n6msq9ogpsh3vr+6khlwfwsusky7zla912ti93trp09y6u9+k+cax4sbcggbmuqub8bha0aamce2v9kz7q7f0we/152ootpbzrknh6s/+afek+gm=</latexit> <latexit sha1_base64="x0vvtaelfyxek+dgehkmhm6lpnw=">aaaeznichvhbbtnaej0aa8ufmsificuktfv6idyicysevmsfy9qstljdrfzmk67q2ja9sqiwxrvx4qiie0gcej/bhr/gup4acsxslxwyb6yantistd7zmfdm38w4gssisenjgnzop3/h4uily+nylavltzvru5e/dblvm9/1w33hjrgrpn6wqrp8pwi5pxbcvuccp8riezepi+f7t+qo4icdu++jnmc2rfdnrdmenuqgeuamy+u92taisrzef15ih6e9shm2wsngckfknzi4awkxfzalwoamofkf/hc4izgfnpw5+eojxnzblfksuqyu75zfykte+exnuf0fnqvavx2zwwl1ellpcq+bz9awqte/khudwp02qvpk2jbzow75avm172bbf3xgmiqbcpbaou2cdrf+b2achqb9h5cgl0rlqdihfazkdhhfewgj9xj/fbwd5f4p71noslezvulidpfjyi1+o5/okf1kp9mf9hdlrktlylsm8htgxb50ll/e2dn7l2uap4sjv6y5mix04j7sklb7odxzfwzmj5+9pt25v72wrnmp9cfqf09p6beswit/sy9bfpvdhd0oahl3riv3hp4xngxfydpsxz4n6msq9ogpsh3vr+6khlwfwsusky7zla912ti93trp09y6u9+k+cax4sbcggbmuqub8bha0aamce2v9kz7q7f0we/152ootpbzrknh6s/+afek+gm=</latexit> <latexit sha1_base64="x0vvtaelfyxek+dgehkmhm6lpnw=">aaaeznichvhbbtnaej0aa8ufmsificuktfv6idyicysevmsfy9qstljdrfzmk67q2ja9sqiwxrvx4qiie0gcej/bhr/gup4acsxslxwyb6yantistd7zmfdm38w4gssisenjgnzop3/h4uily+nylavltzvru5e/dblvm9/1w33hjrgrpn6wqrp8pwi5pxbcvuccp8riezepi+f7t+qo4icdu++jnmc2rfdnrdmenuqgeuamy+u92taisrzef15ih6e9shm2wsngckfknzi4awkxfzalwoamofkf/hc4izgfnpw5+eojxnzblfksuqyu75zfykte+exnuf0fnqvavx2zwwl1ellpcq+bz9awqte/khudwp02qvpk2jbzow75avm172bbf3xgmiqbcpbaou2cdrf+b2achqb9h5cgl0rlqdihfazkdhhfewgj9xj/fbwd5f4p71noslezvulidpfjyi1+o5/okf1kp9mf9hdlrktlylsm8htgxb50ll/e2dn7l2uap4sjv6y5mix04j7sklb7odxzfwzmj5+9pt25v72wrnmp9cfqf09p6beswit/sy9bfpvdhd0oahl3riv3hp4xngxfydpsxz4n6msq9ogpsh3vr+6khlwfwsusky7zla912ti93trp09y6u9+k+cax4sbcggbmuqub8bha0aamce2v9kz7q7f0we/152ootpbzrknh6s/+afek+gm=</latexit> <latexit sha1_base64="x0vvtaelfyxek+dgehkmhm6lpnw=">aaaeznichvhbbtnaej0aa8ufmsificuktfv6idyicysevmsfy9qstljdrfzmk67q2ja9sqiwxrvx4qiie0gcej/bhr/gup4acsxslxwyb6yantistd7zmfdm38w4gssisenjgnzop3/h4uily+nylavltzvru5e/dblvm9/1w33hjrgrpn6wqrp8pwi5pxbcvuccp8riezepi+f7t+qo4icdu++jnmc2rfdnrdmenuqgeuamy+u92taisrzef15ih6e9shm2wsngckfknzi4awkxfzalwoamofkf/hc4izgfnpw5+eojxnzblfksuqyu75zfykte+exnuf0fnqvavx2zwwl1ellpcq+bz9awqte/khudwp02qvpk2jbzow75avm172bbf3xgmiqbcpbaou2cdrf+b2achqb9h5cgl0rlqdihfazkdhhfewgj9xj/fbwd5f4p71noslezvulidpfjyi1+o5/okf1kp9mf9hdlrktlylsm8htgxb50ll/e2dn7l2uap4sjv6y5mix04j7sklb7odxzfwzmj5+9pt25v72wrnmp9cfqf09p6beswit/sy9bfpvdhd0oahl3riv3hp4xngxfydpsxz4n6msq9ogpsh3vr+6khlwfwsusky7zla912ti93trp09y6u9+k+cax4sbcggbmuqub8bha0aamce2v9kz7q7f0we/152ootpbzrknh6s/+afek+gm=</latexit> (a) A(t) = B C A (b) C1 C2 C3 Fig. 1. A typical example of (a) a Delaunay network and (b) coarse-grained formations obtained from the 13th section of Sendai. 8

9 fukuoka hiroshima section 13 (541) section 14 (343) section 15 (343) section 16 (541) section 17 (424) section 16 (343) section 17 (343) iwata kashima kashiwa kawasaki section 13 (343) section 14 (343) section 15 (others) section 16 (others) section 17 (424) section 16 (others) section 17 (424) kobe kofu section 13 (541) section 14 (541) section 15 (others) section 16 (541) section 17 (541) 9

10 nagoya section 13 (4141) section 14 (4141) section 15 (4141) niigata section 13 (343) section 14 (424) section 15 (541) section 16 (4141) section 17 (4141) section 16 (424) section 17 (541) omiya section 13 (424) section 14 (424) section 15 (others) osaka sendai shonan section 13 (343) section 14 (343) section 15 (343) section 16 (343) section 17 (424) tokyo section 13 (424) section 14 (others) section 15 (others) tosu section 13 (433) section 14 (433) section 15 (433) section 16 (433) section 17 (433) 10

11 urawa section 13 (343) section 14 (343) section 15 (343) yokohama section 13 (4141) section 14 (4141) section 15 (424) section 16 (343) section 17 (343) section 16 (424) section 17 (4141) Fig. 2. Heat maps (average formations) of five games for 18 teams. The direction of offense is upward in each panel. Each heat map belongs to one of the following five average formations: 424, 4141, 433, 541, and 343. The label (others) means that player substitutions occurred in the first half of the game, or the average formation could not be identified. (a) Section 13 Sendai Section 13 Nagoya Section 13 Tosu Section 13 Kofu Section 13 Shonan (b) Fig. 3. (a) Schematic representation of each average formation (upper panels) and typical examples thereof (lower panels). (b) average formations of each team throughout five games. The label others means that player substitutions occurred in the first half or the average formation could not be identified. 11

12 (a) (b) Fig. 4. (a) Dendrogram and (b) the relation between the number of clusters and height obtained from the hierarchical clustering for five games of Sendai. The height of the vertical axes corresponds to the distance between two merged clusters in the clustering process. (a) N c = 4 C1 C2 C3 C4 (b) N c = 8 C1 C2 C3 C4 C5 C6 C7 C8 Fig. 5. Visualization of coarse-grained formations for five games of Sendai where (a) N c = 4 and (b) N c = 8. Each cluster is distinguished by a cluster number. 12

13 Section 13 Section 14 Section 15 Section 16 Section 17 C3 C3 C3 C3 C3 C4 C2 C4 C2 C4 C2 C4 C2 C4 C2 C5 C1 C5 C1 C5 C1 C5 C1 C5 C1 C6 C8 C6 C8 C6 C8 C6 C8 C6 C8 C7 C7 C7 C7 C7 Fig. 6. Transition diagram between clusters for all games of Sendai where N c = 8. 13