# Parabolic Detection Algorithm of Tennis Serve Based on Video Image Analysis Technology

Authors: Hongxin Tang
Journal: Security and Communication Networks
 (2021)

---

## Abstract

At present, the existing algorithm for detecting the parabola of tennis serves neglects the pre-estimation of the global motion information of tennis balls, which leads to great error and low recognition rate. Therefore, a new algorithm for detecting the parabola of tennis service based on video image analysis is proposed. The global motion information is estimated in advance, and the motion feature of the target is extracted. A tennis appearance model is established by sparse representation, and the data of high-resolution tennis flight appearance model are processed by data fusion technology to track the parabolic trajectory. Based on the analysis of the characteristics of the serve mechanics, according to the nonlinear transformation of the parabolic trajectory state vector, the parabolic trajectory starting point is determined, the parabolic trajectory is obtained, and the detection algorithm of the parabolic service is designed. Experimental results show that compared with the other two algorithms, the algorithm designed in this paper can recognize the trajectory of the parabola at different stages, and the detection accuracy of the parabola is higher in the three-dimensional space of the tennis service.


---

## Body

## 1. Introduction

With the development of the economy, tennis has become more and more popular in China. Every year, several large-scale tennis matches (professional and amateur) are held in each province, especially amateur matches. Serving is one of the most important techniques in modern tennis [1]. Serving is not only the beginning of a game, but also the only process of hitting that is controlled by oneself and not controlled by the opponent. Excellent serve can not only score directly, but also exert greater personal technical characteristics to control the opponent and maximize tactical intent. From the point of view of system theory, tennis service not only refers to the technical layer but also includes the tactical layer, psychological layer, and physical layer. Technical, tactical, psychological, and physical dimensions are interrelated, interdependent, and indivisible as a whole [2]. Parabolic trajectory is also the direct factor that affects the result of tennis match. Therefore, it is necessary to track and detect the parabolic trajectory of tennis service.Fu and Guan [3] proposed a particle filter tennis motion tracking method based on Kalman filter prediction. Based on the excellent localization characteristics of multi-scale wavelet transform in time domain and space domain, the adjacent video images are differentiated to extract the target feature information reflecting the foreground motion, so as to overcome the adverse factors of illumination change and tennis motion scale changing with time. At the same time, based on the structural characteristics of tennis court, the influence of adverse interference factors outside the court is excluded. On this basis, the Kalman filter is used to predict and modify the particles, integrate the current observation information into the particle filter process, estimate the mean and covariance of the predicted particle state, make the dynamic particles closer to the posterior probability distribution, and improve the tracking accuracy of tennis moving target. Hua et al. [4] proposed an improved multi-target detection framework LSTM-SSD, which is dedicated to traffic scene video multi-target detection. The single image detection framework is combined with recurrent neural network and LSTM network to form an interleaved cyclic convolution structure. By using a bottleneck LSTM layer to extract the feature mapping between propagation frames, the temporal correlation of network frame level information is realized, which greatly reduces the network computing cost. Combining the time perception information with the improved dynamic Kalman filter algorithm, the tracking and recognition of targets affected by strong interference such as illumination change and large-area occlusion in video is realized.However, the above two traditional methods ignore the pre-estimation of the global motion information of tennis, resulting in low accuracy of serve parabolic detection results and unsatisfactory trajectory recognition rate in different directions. Therefore, a tennis serve parabolic detection algorithm based on video image analysis technology is proposed.Our contribution is threefold:(1)
At present, the existing algorithm for detecting the parabola of tennis serves neglects the pre-estimation of the global motion information of tennis balls, which leads to great error and low recognition rate. Therefore, a new algorithm for detecting the parabola of tennis service based on video image analysis is proposed.(2)
A tennis appearance model is established by sparse representation, and the data of high-resolution tennis flight appearance model are processed by data fusion technology to track the parabolic trajectory.(3)
Based on the analysis of the characteristics of the serve mechanics, according to the nonlinear transformation of the parabolic trajectory state vector, the parabolic trajectory starting point is determined, the parabolic trajectory is obtained, and the detection algorithm of the parabolic service is designed.The remainder of this paper is organized as follows. Section2 introduces the construction of parabolic trajectory tracking model for tennis serve. Section 3 discusses the parabolic detection algorithm for tennis serve. Section 4 discusses experiment and analysis. Section 5 presents the conclusions of the study.


## 2. Construction of Parabolic Trajectory Tracking Model for Tennis Serve

### 2.1. Video Target Motion Feature Extraction

In order to obtain the complete semantic information of the video moving target, the local motion features of targets in the video image are extracted to ensure the effective behavior recognition of subsequent video targets [5, 6]. In order to extract the local motion information of video object effectively, the global motion information must be estimated in advance, and the global motion information must be eliminated when analyzing the local target motion [7]. In this paper, a fast diamond search and matching method is used to perform the global motion extraction, and then the local dense motion vector field is deduced by optical flow analysis.The basic concept of optical flow field is to take the moving image functionfx,y,z as a continuous function of variables x, y, and z. In the case of period t+Δt, the point moves to a new position, the orientation in the image becomes x+Δx,y+Δy, and the gray value is fx+Δx,y+Δy,z+Δz. The optical flow convergence formula is constructed according to the image gray conservation theorem, and the following results are obtained:(1)fx,y,z=fx+Δx,y+Δy,z+Δz.Extend the right side of the equation using Taylor’s formula:(2)fx,y,z=fx,y,z+∂fx,y,z∂xΔx+∂fx,y,z∂yΔy+∂fx,y,z∂t+Δz+e,(3)ux,y,z=ΔxΔt,vx,y,z=ΔyΔt.Assuming that equation (3) is true, the relationship shown in equation (4) can be obtained.(4)∂fx,y,z∂x+∂fx,y,z∂y+∂fx,y,z∂z=0.Since∇f is the gradient of pixel x,y, equation (4) can be reduced to(5)∂f∂t=−∇f⋅u,v,where u,v represents the motion vector of pixel x,y in period t. Equation (5) is the optical flow convergence formula.The optical flow field caused by the same moving object should have the characteristics of continuity and smoothness, that is, if the rate of adjacent points in the same object is the same, the optical flow change projected onto the image should also be smooth [8, 9]. According to this characteristic, an additional convergence condition added to the optical flow field is designed in this paper:(6)minG=∂u∂x2+∂u∂y2+∂v∂x2+∂v∂y2.It can be seen from the above formula that the global smooth convergence can transform the optical flow field operation problem into an optimization problem. Replace the optical flow deviation and rate field gradient deviation [10] in equation (6) and record the two deviations as(7)e′=∂f∂xu+∂f∂yv+∂f∂t.


### 2.2. Appearance Model Based on Sparse Representation

Based on the sparse representation model, the appearance model of flying tennis was constructed by using the target feature extraction algorithm [11]. When extracting the features of tennis flight image, the local preserving mapping algorithm is used to reduce the dimension of tennis flight image data while retaining the original data features.In the feature extraction of tennis flight image, firstly, the K-nearest neighbor method is used to construct the interclass adjacency graph and intraclass adjacency graph [12, 13].Then, determine the edge weight. The algorithm used is as follows:(8)Si,j=1xiandxjareneighborstoeachother0otherwise,where Si,j represents the weight on the edge and xi and xj represent the edges of tennis flight image features, respectively.Then, build a flying tennis appearance model based on the sparse representation model, that is, build a flying tennis appearance model through a sparse reconstruction algorithm. The specific build process for the appearance model is as follows:(1)
Input the features of tennis flight images and train the dictionary pairs with higher resolution and lower resolution, which are represented byH1,H2 and U1,U2k=1N.(2)
For each image blockYiC (size is b×b) in the tennis flight image Y, each image block is processed sequentially from the image block in the upper left corner.(3)
Calculate theYi-means.(4)
Extract the feature block ofYi.(5)
Calculate dictionariesvecZi and H as follows:(9)vecZi=vecX˜ivecF1Yi⋮vecFkYi⋮vecFNYi,(10)H=H2h⊗H1hU21⊗U11⋮U2k⊗U1k⋮U2N⊗U1N,whereX˜i represents the reconstructed high-resolution tennis flight image of the reconstructed area.(6)
Build the image block of high-resolution tennis flight appearance model, as shown in the following formula:(11)Xi=BiTH1h⋅U2hT,whereXi represents the constructed high-resolution tennis flight appearance model image block and T represents the sparsity balance parameter.The repeated operation of partial weighted average is implemented on the image block of high-resolution tennis flight appearance model to construct high-resolution tennis flight appearance modelX0.(7)
Output the high-resolution tennis flight appearance modelX0.


### 2.3. Establishing Tracking Algorithm Model

Firstly, the data of high-resolution tennis flight appearance model are processed by data fusion technology [14, 15]. The data fusion algorithm used is D-S evidence theory, which mainly divides the evidence set into multiple irrelevant parts, independently judges the identification framework through the divided parts, and then recombines these parts through Dempster rules. The specific combination rules are as follows:(12)Z=m1Aim2Bjk⋅A≠0,mϕ=0Ai∩Bj=A,where Ai and Bj, respectively, represent two independent evidence sources; A stands for proposition; mϕ represents the universe set of trust functions; m represents trust function; and k stands for the number of propositions.Then, the parabolic trajectory tracking model of tennis serve is constructed by particle filter algorithm [16, 17]. In the model construction, the measurement covariance has an important impact on the final filter output. In order to avoid the influence of measurement covariance, dynamic correction is introduced, as shown in the following formula:(13)Rk=VI,where Rk represents the measurement covariance; V⋅ represents the dynamic correction function; and I represents the distance between the binocular camera and the tennis ball at a certain time in the three-dimensional space, that is, the parabolic trajectory tracking algorithm model of tennis ball serving, as shown in the following formula:(14)dk=xk−xc2+yk−yc2+zk−zc2,where xk, yk, and zk represents the coordinate measurement value of the tennis ball flying in the three-dimensional world coordinate system; xc, yc, and z represent the midpoint position when the two optical centers are connected with the left and right cameras.


## 2.1. Video Target Motion Feature Extraction

In order to obtain the complete semantic information of the video moving target, the local motion features of targets in the video image are extracted to ensure the effective behavior recognition of subsequent video targets [5, 6]. In order to extract the local motion information of video object effectively, the global motion information must be estimated in advance, and the global motion information must be eliminated when analyzing the local target motion [7]. In this paper, a fast diamond search and matching method is used to perform the global motion extraction, and then the local dense motion vector field is deduced by optical flow analysis.The basic concept of optical flow field is to take the moving image functionfx,y,z as a continuous function of variables x, y, and z. In the case of period t+Δt, the point moves to a new position, the orientation in the image becomes x+Δx,y+Δy, and the gray value is fx+Δx,y+Δy,z+Δz. The optical flow convergence formula is constructed according to the image gray conservation theorem, and the following results are obtained:(1)fx,y,z=fx+Δx,y+Δy,z+Δz.Extend the right side of the equation using Taylor’s formula:(2)fx,y,z=fx,y,z+∂fx,y,z∂xΔx+∂fx,y,z∂yΔy+∂fx,y,z∂t+Δz+e,(3)ux,y,z=ΔxΔt,vx,y,z=ΔyΔt.Assuming that equation (3) is true, the relationship shown in equation (4) can be obtained.(4)∂fx,y,z∂x+∂fx,y,z∂y+∂fx,y,z∂z=0.Since∇f is the gradient of pixel x,y, equation (4) can be reduced to(5)∂f∂t=−∇f⋅u,v,where u,v represents the motion vector of pixel x,y in period t. Equation (5) is the optical flow convergence formula.The optical flow field caused by the same moving object should have the characteristics of continuity and smoothness, that is, if the rate of adjacent points in the same object is the same, the optical flow change projected onto the image should also be smooth [8, 9]. According to this characteristic, an additional convergence condition added to the optical flow field is designed in this paper:(6)minG=∂u∂x2+∂u∂y2+∂v∂x2+∂v∂y2.It can be seen from the above formula that the global smooth convergence can transform the optical flow field operation problem into an optimization problem. Replace the optical flow deviation and rate field gradient deviation [10] in equation (6) and record the two deviations as(7)e′=∂f∂xu+∂f∂yv+∂f∂t.


## 2.2. Appearance Model Based on Sparse Representation

Based on the sparse representation model, the appearance model of flying tennis was constructed by using the target feature extraction algorithm [11]. When extracting the features of tennis flight image, the local preserving mapping algorithm is used to reduce the dimension of tennis flight image data while retaining the original data features.In the feature extraction of tennis flight image, firstly, the K-nearest neighbor method is used to construct the interclass adjacency graph and intraclass adjacency graph [12, 13].Then, determine the edge weight. The algorithm used is as follows:(8)Si,j=1xiandxjareneighborstoeachother0otherwise,where Si,j represents the weight on the edge and xi and xj represent the edges of tennis flight image features, respectively.Then, build a flying tennis appearance model based on the sparse representation model, that is, build a flying tennis appearance model through a sparse reconstruction algorithm. The specific build process for the appearance model is as follows:(1)
Input the features of tennis flight images and train the dictionary pairs with higher resolution and lower resolution, which are represented byH1,H2 and U1,U2k=1N.(2)
For each image blockYiC (size is b×b) in the tennis flight image Y, each image block is processed sequentially from the image block in the upper left corner.(3)
Calculate theYi-means.(4)
Extract the feature block ofYi.(5)
Calculate dictionariesvecZi and H as follows:(9)vecZi=vecX˜ivecF1Yi⋮vecFkYi⋮vecFNYi,(10)H=H2h⊗H1hU21⊗U11⋮U2k⊗U1k⋮U2N⊗U1N,whereX˜i represents the reconstructed high-resolution tennis flight image of the reconstructed area.(6)
Build the image block of high-resolution tennis flight appearance model, as shown in the following formula:(11)Xi=BiTH1h⋅U2hT,whereXi represents the constructed high-resolution tennis flight appearance model image block and T represents the sparsity balance parameter.The repeated operation of partial weighted average is implemented on the image block of high-resolution tennis flight appearance model to construct high-resolution tennis flight appearance modelX0.(7)
Output the high-resolution tennis flight appearance modelX0.


## 2.3. Establishing Tracking Algorithm Model

Firstly, the data of high-resolution tennis flight appearance model are processed by data fusion technology [14, 15]. The data fusion algorithm used is D-S evidence theory, which mainly divides the evidence set into multiple irrelevant parts, independently judges the identification framework through the divided parts, and then recombines these parts through Dempster rules. The specific combination rules are as follows:(12)Z=m1Aim2Bjk⋅A≠0,mϕ=0Ai∩Bj=A,where Ai and Bj, respectively, represent two independent evidence sources; A stands for proposition; mϕ represents the universe set of trust functions; m represents trust function; and k stands for the number of propositions.Then, the parabolic trajectory tracking model of tennis serve is constructed by particle filter algorithm [16, 17]. In the model construction, the measurement covariance has an important impact on the final filter output. In order to avoid the influence of measurement covariance, dynamic correction is introduced, as shown in the following formula:(13)Rk=VI,where Rk represents the measurement covariance; V⋅ represents the dynamic correction function; and I represents the distance between the binocular camera and the tennis ball at a certain time in the three-dimensional space, that is, the parabolic trajectory tracking algorithm model of tennis ball serving, as shown in the following formula:(14)dk=xk−xc2+yk−yc2+zk−zc2,where xk, yk, and zk represents the coordinate measurement value of the tennis ball flying in the three-dimensional world coordinate system; xc, yc, and z represent the midpoint position when the two optical centers are connected with the left and right cameras.


## 3. Parabolic Detection Algorithm for Tennis Serve

### 3.1. Stress Analysis of Tennis

When analyzing the force on a tennis ball, firstly, the rotating axis of the tennis ball is defined, and the object of the study is a certain state of the tennis ball stagnating in the parabolic trajectory of the actual service. Combined with the flow velocity of air during the course of the tennis ball rotation, the coordinate system is established as follows. TheZ axis of the rotating tennis ball is the coordinate system, and the force coordinate system shown in Figure 1 is obtained.Figure 1
Tennis coordinate system established.From the tennis coordinate system established in Figure1, the velocity of the tennis is set to the X axis, and the pressure on the tennis is set to the Y axis. According to the direction of the coordinate system, it can be analyzed that the tennis ball is mainly affected by gravity, air resistance, and Magnus force in the air. Define time t of the tennis ball shown in Figure 1 above. The speed of the tennis ball is(15)vt=vxt,vyt,vztT,where T represents the flight time of the tennis ball and vxt,vyt,vzt represent the movement speeds in different axes, respectively. Then, the rotation angular velocity of the tennis ball is calculated, and the calculation formula is(16)W=ωx,ωy,ωzT,where ωx,ωy,ωz, respectively, represent the angular velocities of tennis balls in different axes. Based on the above values of angular velocity and rotation speed, the gravity, air resistance, and MAG effort of the tennis ball at that time are calculated:(17)Fg=−m00gTFa=12CdρAvtFm=12CmρrAW×vt,where m is the weight of the tennis ball, g is the gravitational acceleration, Cd is the air resistance coefficient, Cm is the Magnus force coefficient, ρ is the air density under the standard air pressure, A is the cross-sectional area of the tennis ball, r is the radius of the tennis ball, and the meaning of other parameters remains unchanged.From the above formula, we can see that the air resistance is proportional to the square of the speed of the tennis. When the speed of one side of the tennis ball is not parallel to the direction of the spin speed, under the action of air on one side of the tennis ball, the superposition of the flying speed and the spin speed becomes larger, and the speed of the other side of the tennis ball is offset by the action of air. The parabolic trajectory of the serve is to track the tennis movement process obtained by the force analysis above and, finally, obtain the parabolic trajectory and mechanical analysis of different tennis balls.


### 3.2. Parabolic Trajectory Detection of Tennis Serve

Under the tennis motion expression obtained from the above analysis, the tennis stop rotation process is defined to represent the end of the iterative process, and the noise generated by the iterative process is calculated. The calculation formula is(18)K=12π⋅ρm⋅C⋅D3,where C represents the lift coefficient, D represents the diameter of the tennis ball, and the meaning of other parameters remains unchanged. Under the influence of the above noise value, predict the state vector of the tennis ball and take the above calculated noise value as a priori estimation value, and the state vector can be expressed as [18, 19](19)Xk|k−1=FXk−1,Uk−1Pk|k−1=Ak−1Pk−1+Qk,where Xk|k−1 represents a priori estimated value of tennis state variable, Pk|k−1 represents a priori estimated value of covariance of state variable, Ak represents partial derivative of motion process function F to state variable, and Qk represents noise vector. According to the above state vector values, the number of state data of different tennis balls is divided, as shown in Table 1.Table 1
Data vectors in different rotation states.
Serial numberRotation attributeOutput data vector1Internal rotation042External rotation013Irrotational07


### 3.3. Parabolic Trajectory Acquisition of Tennis Serve

According to the nonlinear transformed state vector of the serving parabolic trajectory, the starting point of the serving parabolic trajectory is determined. In order to simplify the calculation process in the tracking process of the serving parabolic trajectory, the distribution distance is regarded as a constant, and the force on the tennis ball can be expressed as(20)Ftotal=Ft,κk2,where κk represents the level parameter of class in the process of continuous sampling motion. According to the tracking results of the above collection points, the coordinates of the tennis ball at the sampling point are determined, and the coordinate values of the collection points are summarized as shown in Table 2.Table 2
Collection point coordinates of parabolic trajectory of tennis serve.
Collection point nameX-axis value (m)Y-axis value (m)Z-axis value (m)17.1351.1880.39327.3221.0934.20433.7161.9180.77344.2370.3992.72174.1849.7011.03832.3498.0333.79374.0127.1937.31884.4053.3019.78797.3797.4213.728100.2747.0721.029113.1982.0282.707127.1929.3227.823138.9157.0792.002147.3160.7730.724177.7483.1331.373133.4832.8733.819177.1973.7037.729187.1398.3417.818197.2799.3084.878207.7077.8777.897The coordinates of the tennis ball in three directions shown in Table2 are taken as the parabolic trajectory dataset, the coordinate values in the above table are processed by nonlinear filtering, and the parabolic trajectory of the tennis ball is continuously updated. The update expression is as follows:(21)Kn=FSxSySzcosθ,where Sx, Sy, and Sz represent the values of three axes, respectively, and θ represents the included angle between the parabolic trajectory of the actual tennis serve and the three axes. According to the calculation formula of sampling points formed by the above integration, the parabolic trajectory of tennis serve is simulated in a three-dimensional coordinate, as shown in Figure 2.Figure 2
Projectile trajectory of a tennis serve.As shown in Figure2, the projectile trajectories of different tennis balls in tennis sports are finally studied on the premise of mechanical analysis.


## 3.1. Stress Analysis of Tennis

When analyzing the force on a tennis ball, firstly, the rotating axis of the tennis ball is defined, and the object of the study is a certain state of the tennis ball stagnating in the parabolic trajectory of the actual service. Combined with the flow velocity of air during the course of the tennis ball rotation, the coordinate system is established as follows. TheZ axis of the rotating tennis ball is the coordinate system, and the force coordinate system shown in Figure 1 is obtained.Figure 1
Tennis coordinate system established.From the tennis coordinate system established in Figure1, the velocity of the tennis is set to the X axis, and the pressure on the tennis is set to the Y axis. According to the direction of the coordinate system, it can be analyzed that the tennis ball is mainly affected by gravity, air resistance, and Magnus force in the air. Define time t of the tennis ball shown in Figure 1 above. The speed of the tennis ball is(15)vt=vxt,vyt,vztT,where T represents the flight time of the tennis ball and vxt,vyt,vzt represent the movement speeds in different axes, respectively. Then, the rotation angular velocity of the tennis ball is calculated, and the calculation formula is(16)W=ωx,ωy,ωzT,where ωx,ωy,ωz, respectively, represent the angular velocities of tennis balls in different axes. Based on the above values of angular velocity and rotation speed, the gravity, air resistance, and MAG effort of the tennis ball at that time are calculated:(17)Fg=−m00gTFa=12CdρAvtFm=12CmρrAW×vt,where m is the weight of the tennis ball, g is the gravitational acceleration, Cd is the air resistance coefficient, Cm is the Magnus force coefficient, ρ is the air density under the standard air pressure, A is the cross-sectional area of the tennis ball, r is the radius of the tennis ball, and the meaning of other parameters remains unchanged.From the above formula, we can see that the air resistance is proportional to the square of the speed of the tennis. When the speed of one side of the tennis ball is not parallel to the direction of the spin speed, under the action of air on one side of the tennis ball, the superposition of the flying speed and the spin speed becomes larger, and the speed of the other side of the tennis ball is offset by the action of air. The parabolic trajectory of the serve is to track the tennis movement process obtained by the force analysis above and, finally, obtain the parabolic trajectory and mechanical analysis of different tennis balls.


## 3.2. Parabolic Trajectory Detection of Tennis Serve

Under the tennis motion expression obtained from the above analysis, the tennis stop rotation process is defined to represent the end of the iterative process, and the noise generated by the iterative process is calculated. The calculation formula is(18)K=12π⋅ρm⋅C⋅D3,where C represents the lift coefficient, D represents the diameter of the tennis ball, and the meaning of other parameters remains unchanged. Under the influence of the above noise value, predict the state vector of the tennis ball and take the above calculated noise value as a priori estimation value, and the state vector can be expressed as [18, 19](19)Xk|k−1=FXk−1,Uk−1Pk|k−1=Ak−1Pk−1+Qk,where Xk|k−1 represents a priori estimated value of tennis state variable, Pk|k−1 represents a priori estimated value of covariance of state variable, Ak represents partial derivative of motion process function F to state variable, and Qk represents noise vector. According to the above state vector values, the number of state data of different tennis balls is divided, as shown in Table 1.Table 1
Data vectors in different rotation states.
Serial numberRotation attributeOutput data vector1Internal rotation042External rotation013Irrotational07


## 3.3. Parabolic Trajectory Acquisition of Tennis Serve

According to the nonlinear transformed state vector of the serving parabolic trajectory, the starting point of the serving parabolic trajectory is determined. In order to simplify the calculation process in the tracking process of the serving parabolic trajectory, the distribution distance is regarded as a constant, and the force on the tennis ball can be expressed as(20)Ftotal=Ft,κk2,where κk represents the level parameter of class in the process of continuous sampling motion. According to the tracking results of the above collection points, the coordinates of the tennis ball at the sampling point are determined, and the coordinate values of the collection points are summarized as shown in Table 2.Table 2
Collection point coordinates of parabolic trajectory of tennis serve.
Collection point nameX-axis value (m)Y-axis value (m)Z-axis value (m)17.1351.1880.39327.3221.0934.20433.7161.9180.77344.2370.3992.72174.1849.7011.03832.3498.0333.79374.0127.1937.31884.4053.3019.78797.3797.4213.728100.2747.0721.029113.1982.0282.707127.1929.3227.823138.9157.0792.002147.3160.7730.724177.7483.1331.373133.4832.8733.819177.1973.7037.729187.1398.3417.818197.2799.3084.878207.7077.8777.897The coordinates of the tennis ball in three directions shown in Table2 are taken as the parabolic trajectory dataset, the coordinate values in the above table are processed by nonlinear filtering, and the parabolic trajectory of the tennis ball is continuously updated. The update expression is as follows:(21)Kn=FSxSySzcosθ,where Sx, Sy, and Sz represent the values of three axes, respectively, and θ represents the included angle between the parabolic trajectory of the actual tennis serve and the three axes. According to the calculation formula of sampling points formed by the above integration, the parabolic trajectory of tennis serve is simulated in a three-dimensional coordinate, as shown in Figure 2.Figure 2
Projectile trajectory of a tennis serve.As shown in Figure2, the projectile trajectories of different tennis balls in tennis sports are finally studied on the premise of mechanical analysis.


## 4. Analysis of Experimental Results

### 4.1. Tennis Parabolic Motion Recognition Rate

In the experiment, the flying image of the robot tennis ball is collected, and the image is processed in a series of ways, which is used to detect the trajectory of the tennis ball in the experiment. Experimental images are shown in Figure3.Figure 3
Tennis projectile motion simulation diagram.In order to verify the effectiveness of the algorithm in this paper, the algorithm in this paper, the tennis moving target tracking algorithm based on Kalman predictive particle filter proposed in literature [3], and the identification results of different stages of tennis projectile motion by video multi-target detection technology based on recursive neural network proposed in literature [4] are tested successively, as shown in Figure 4.Figure 4
Comparison of recognition rates of different algorithms in different stages of tennis projectile motion.Figure4 shows that the tennis parabolic motion recognition result generated by this algorithm is better than the other two documents. This is because the algorithm fully considers the use of camera calibration in the process of generating the neutral parabolic trajectory and the target dynamic analysis of the neutrality of the target motion state parabolic trajectory to improve the accuracy of video action recognition.


### 4.2. Parabolic Trajectory Tracking Results of Tennis Serve inX, Y, and Z Directions

#### 4.2.1.X-Axis Test Results

Using the designed tennis serve parabolic trajectory detection algorithm, the serve parabolic trajectory of tennis thrown by the experimental tennis robot is tracked. In the tracking, the parabolic trajectory detection of tennis serve is tested from three angles ofX axis, Y axis, and Z axis. The test results of parabolic trajectory detection of tennis serve on X axis are shown in Figure 5.Figure 5
Test results of parabolic trajectory detection of tennis serve onX axis.From theX-axis tennis serve parabolic trajectory detection experimental results in Figure 5, it can be seen that there is almost no significant difference between the tennis serve parabolic trajectory detection output and the actual flight serve parabolic trajectory of tennis. With the passage of time, the tracking result of the design algorithm is still relatively stable and the output noise is small, which proves that the detection effect of X-axis tennis serve parabolic trajectory of the design algorithm is better.


#### 4.2.2.Y-Axis Test Results

Then, the parabolic trajectory of theY-axis tennis service is tested, and the test results are shown in Figure 6.Figure 6
Parabolic trajectory detection ofY-axis tennis serve.Based on the detection of the parabolic trajectory of theY-axis tennis service in Figure 6, it is found that the algorithm based on data fusion and sparse representation model can track the Y axis accurately in 3D world coordinate system. It is proved that the designed algorithm can detect the parabolic trajectory of the Y-axis service with high precision because the result is very close to the reality.


#### 4.2.3.Z-Axis Test Results

Finally, theZ-axis tennis serve parabolic trajectory detection of the designed algorithm is tested, and the specific test results are shown in Figure 7.Figure 7
Parabolic trajectory detection ofZ-axis tennis serve.From the parabolic trajectory detection of theZ-axis tennis ball in Figure 7, we can see that the accuracy of the parabolic trajectory tracking of the Z-axis tennis ball is lower than that of the X-axis and Y-axis tennis ball, and the whole parabolic trajectory tracking has some errors, but the accuracy is still high.Combining the detection results of the three direction axes, it can be found that the output of the tennis serve parabolic detection algorithm is very close to the reality, that is, the designed tennis serve parabolic trajectory tracking algorithm based on data fusion and sparse representation model can detect the tennis serve parabolic trajectory more accurately.


## 4.1. Tennis Parabolic Motion Recognition Rate

In the experiment, the flying image of the robot tennis ball is collected, and the image is processed in a series of ways, which is used to detect the trajectory of the tennis ball in the experiment. Experimental images are shown in Figure3.Figure 3
Tennis projectile motion simulation diagram.In order to verify the effectiveness of the algorithm in this paper, the algorithm in this paper, the tennis moving target tracking algorithm based on Kalman predictive particle filter proposed in literature [3], and the identification results of different stages of tennis projectile motion by video multi-target detection technology based on recursive neural network proposed in literature [4] are tested successively, as shown in Figure 4.Figure 4
Comparison of recognition rates of different algorithms in different stages of tennis projectile motion.Figure4 shows that the tennis parabolic motion recognition result generated by this algorithm is better than the other two documents. This is because the algorithm fully considers the use of camera calibration in the process of generating the neutral parabolic trajectory and the target dynamic analysis of the neutrality of the target motion state parabolic trajectory to improve the accuracy of video action recognition.


## 4.2. Parabolic Trajectory Tracking Results of Tennis Serve inX, Y, and Z Directions

### 4.2.1.X-Axis Test Results

Using the designed tennis serve parabolic trajectory detection algorithm, the serve parabolic trajectory of tennis thrown by the experimental tennis robot is tracked. In the tracking, the parabolic trajectory detection of tennis serve is tested from three angles ofX axis, Y axis, and Z axis. The test results of parabolic trajectory detection of tennis serve on X axis are shown in Figure 5.Figure 5
Test results of parabolic trajectory detection of tennis serve onX axis.From theX-axis tennis serve parabolic trajectory detection experimental results in Figure 5, it can be seen that there is almost no significant difference between the tennis serve parabolic trajectory detection output and the actual flight serve parabolic trajectory of tennis. With the passage of time, the tracking result of the design algorithm is still relatively stable and the output noise is small, which proves that the detection effect of X-axis tennis serve parabolic trajectory of the design algorithm is better.


### 4.2.2.Y-Axis Test Results

Then, the parabolic trajectory of theY-axis tennis service is tested, and the test results are shown in Figure 6.Figure 6
Parabolic trajectory detection ofY-axis tennis serve.Based on the detection of the parabolic trajectory of theY-axis tennis service in Figure 6, it is found that the algorithm based on data fusion and sparse representation model can track the Y axis accurately in 3D world coordinate system. It is proved that the designed algorithm can detect the parabolic trajectory of the Y-axis service with high precision because the result is very close to the reality.


### 4.2.3.Z-Axis Test Results

Finally, theZ-axis tennis serve parabolic trajectory detection of the designed algorithm is tested, and the specific test results are shown in Figure 7.Figure 7
Parabolic trajectory detection ofZ-axis tennis serve.From the parabolic trajectory detection of theZ-axis tennis ball in Figure 7, we can see that the accuracy of the parabolic trajectory tracking of the Z-axis tennis ball is lower than that of the X-axis and Y-axis tennis ball, and the whole parabolic trajectory tracking has some errors, but the accuracy is still high.Combining the detection results of the three direction axes, it can be found that the output of the tennis serve parabolic detection algorithm is very close to the reality, that is, the designed tennis serve parabolic trajectory tracking algorithm based on data fusion and sparse representation model can detect the tennis serve parabolic trajectory more accurately.


## 4.2.1.X-Axis Test Results

Using the designed tennis serve parabolic trajectory detection algorithm, the serve parabolic trajectory of tennis thrown by the experimental tennis robot is tracked. In the tracking, the parabolic trajectory detection of tennis serve is tested from three angles ofX axis, Y axis, and Z axis. The test results of parabolic trajectory detection of tennis serve on X axis are shown in Figure 5.Figure 5
Test results of parabolic trajectory detection of tennis serve onX axis.From theX-axis tennis serve parabolic trajectory detection experimental results in Figure 5, it can be seen that there is almost no significant difference between the tennis serve parabolic trajectory detection output and the actual flight serve parabolic trajectory of tennis. With the passage of time, the tracking result of the design algorithm is still relatively stable and the output noise is small, which proves that the detection effect of X-axis tennis serve parabolic trajectory of the design algorithm is better.


## 4.2.2.Y-Axis Test Results

Then, the parabolic trajectory of theY-axis tennis service is tested, and the test results are shown in Figure 6.Figure 6
Parabolic trajectory detection ofY-axis tennis serve.Based on the detection of the parabolic trajectory of theY-axis tennis service in Figure 6, it is found that the algorithm based on data fusion and sparse representation model can track the Y axis accurately in 3D world coordinate system. It is proved that the designed algorithm can detect the parabolic trajectory of the Y-axis service with high precision because the result is very close to the reality.


## 4.2.3.Z-Axis Test Results

Finally, theZ-axis tennis serve parabolic trajectory detection of the designed algorithm is tested, and the specific test results are shown in Figure 7.Figure 7
Parabolic trajectory detection ofZ-axis tennis serve.From the parabolic trajectory detection of theZ-axis tennis ball in Figure 7, we can see that the accuracy of the parabolic trajectory tracking of the Z-axis tennis ball is lower than that of the X-axis and Y-axis tennis ball, and the whole parabolic trajectory tracking has some errors, but the accuracy is still high.Combining the detection results of the three direction axes, it can be found that the output of the tennis serve parabolic detection algorithm is very close to the reality, that is, the designed tennis serve parabolic trajectory tracking algorithm based on data fusion and sparse representation model can detect the tennis serve parabolic trajectory more accurately.


## 5. Conclusion

In order to solve the problem that the traditional algorithm for detecting the parabola of tennis serve cannot track the track of the ball in different directions, a new algorithm based on video image analysis is proposed. Through sparse representation and data fusion technology, the high-resolution tennis flying appearance model data are processed and the parabolic trajectory of tennis serve is tracked. We determine the starting point of the parabolic trajectory and complete the parabolic testing of the tennis service. Experimental results show that the proposed algorithm can identify the trajectory of tennis service parabola in different stages with high accuracy, which provides a reliable theoretical support for further research.



---

*Source: 7901677-2021-11-29.xml*