# Sports Auxiliary Training Based on Computer Digital 3D Video Image Processing

Authors: Saisai Xu
Journal: Computational Intelligence and Neuroscience
 (2022)

---

## Abstract

With the continuous development of social economy, sports have received more and more attention. How to improve the quality of sports has become the focus of research. The computer digital 3D video image processing is introduced in this paper, taking shooting as the starting point, in which computer digitization technology is used to collect images of sequence targets through combining the operation flow of shooting, monitor the results and data of shooting and process 3D video images, conduct the analyze and mine according to the corresponding statistical processing results, and evaluate the corresponding training. The simulation experiment proves that the computerized digital 3D video image processing is effective and can scientifically support sports-assisted training.


---

## Body

## 1. Introduction

With the development of social economy, sports have received more and more attention as an important way for people to exercise. In competitive competitions, how to improve the quality and effectiveness of sports players is extremely important [1, 2]. Traditional training methods are usually dominated by coaches, and they are judged and corrected through experience to achieve the integration of training rhythms and methods [3, 4]. But, most sports have strong requirements for athletes’ balance, attention, coordination, and sense of time. Therefore, how to quantitatively improve athletes’ training performance is the focus of further investigation [5, 6]. The development of computer technology has led to the development of many assistive technologies, and the computer-related technologies have been introduced into sports training within the industry to further improve the summary of exercise rules and exercise effectiveness, so as to achieve scientific and effective sports training [7, 8].During the process of sports training, computer technology can replace a variety of roles and functions, such as simulating sports and capturing corresponding 3D simulation sports (boxing, table tennis, etc.); 3D simulations and emulations of high jump actions are performed through 3D movies to realize training analysis and analyze the authenticity of the athlete’s movement [9–11]. Therefore, in the actual sports training process, it is necessary to integrate image capture and three-dimensional simulation. Through image processing, cleaning, and analysis, the three-dimensional sports simulation of athletes is realized. At the same time, the corresponding dynamic data and equations are used to simulate athletes’ movements, realize synchronous perspective and synchronous training according to training requirements and needs, and provide reference for sports training [12–14].For shooting sports, it pays more attention to attentiveness, critical state, etc. Traditional training usually uses visual artificial judgment. In actual training, there are often the problems of inaccurate judgments, long time-consuming, and inability to analyze relevant data. Therefore, in response to this need, the computer digital 3D video image processing is introduced in this paper, computer images are used to identify the shooting target ring through the analysis of the shooting movement process, and the shooting results parameter changes are calculated real time to achieve the evaluation of the training results, aiming at providing an auxiliary reference for sports training, so as to improve the quality and effect of training.


## 2. Computer Digital 3D Video Image Processing

For computer digital 3D video image processing, the specific principle is shown in Figure1. First, the simulated data need to be visualized. Secondly, according to the needs, the typical characteristic area shall be selected as the relevant sample data for classification. These samples are stored in the corresponding training network to save the results, and corresponding data recognition is performed through feature acquisition [15, 16]. During the actual processing process, feature samples can be selected according to corresponding needs and demands, and training iterations can be performed to realize feature visualization.Figure 1
Principle of multiresolution visualization.


### 2.1. Feature Detection and Recognition

For feature detection and recognition, real-time interactive processing is required. The feature detection and recognition algorithm is carried out on the CPU, and in order to improve the speed, the method of taking the area of the critical point as the candidate unit is adopted. In order to realize the real-time interaction of feature detection and recognition, this paper designs a GPU-based feature detection and recognition algorithm. The basic idea is to convert the flow field into texture fragment blocks and use the high parallel characteristics and programmability of GPU to convert BP neural network feature recognition into processing texture fragments. The basic flow is shown in Figure2. The algorithm mainly includes the following steps:Figure 2
GPU feature recognition algorithm flow chart.Step 1.
Texture conversion: it is responsible for converting the flow field data into a color texture that is easy to process by the GPU.Step 2.
GPU processing: it is responsible for feature recognition of the area where the current fragment is located.Step 3.
Saving the results: it is responsible for reading the recognition result from the GPU and saving it to the corresponding data structure.
Because the GPU feature recognition process is a parallel process, this paper does not adopt the critical point region candidate unit method, but uses the sequential traversal method. The reason is analyzed as follows: When using the critical point candidate unit method or the traversal method on the GPU, assuming that the texture conversion process time isT¯1 and T1, the GPU processing process is T¯2 and T2, and the result saving process is T¯3 and T3, the entire pipeline The processing time is T¯=T¯1+T¯2+T¯3 and T=T1+T2+T3, respectively. Since Step 1 and Step 3 require the same time for the two methods, namely T¯1=T1, T¯3=T3, the length of the pipeline processing time depends on T¯2 and XT2. The critical point area candidate unit method needs to discriminate the fragment type in the fragment shader. If the current fragment corresponds to the critical point, then the identification judgment is performed; otherwise, the process is skipped. Let's suppose that the fragment shader Fi corresponds to the critical point fragment type, and the processing time is T¯Fi, and the fragment shader Fj corresponds to the non-critical point type, and the processing time is T¯Fj. Although T¯Fi<T¯Fj, because the GPU is a parallel processing process, the shader with the longest calculation time among the fragment shaders constitutes the bottleneck of the recognition algorithm, thus T¯2=maxT¯Fi,T¯Fj=T¯Fj. Similarly, there is T2=maxTFi,0≤i≤N for the traversal method, where TFi is the processing time of each fragment. Because the segment type judgment process sentence is added in the critical point area candidate unit method, there must be T¯2=maxT¯Fi>maxTFi=T2, so T¯>T. That is, in the GPU processing process, the speed of the traversal method is better than that of the critical point candidate unit method. (Algorithm 1).
First, obtain the node data of the characteristic area according to the current texture coordinates, then perform the characteristic recognition calculation, and set the current fragment color according to the recognition result. The implementation code of the feature recognition process is basically the same as that on the CPU, but the texture-speed reverse calculation is first required; in addition, to ensure that the texture conversion process does not lose data accuracy, the texture format uses a 32 bit floating point type.
This paper also uses the pressure field data to correct the characteristics of the BP neural network. In the experiment, the characteristics of cyclones and anticyclones were mainly extracted, and cyclones and anticyclones correspond to the low and high pressure centers in the pressure field, respectively. There is a pressureP≤Pthreshhold for the center of the cyclone and P≥Pthreadhold′ for the center of the anticyclone. For wind field data, assume that the position obtained after feature detection using BP neural network is Pi, and the position obtained after detecting the corresponding pressure field data using the pressure amplitude method is Pi′. If DPi,Pi′≤dthreshhold, the detection result is considered correct and output pi; otherwise, the detection result is considered to be wrong, where dthreshhold is the Euler distance error threshold.Algorithm1: Feature detection and recognition algorithm.
Input: Image input data, network weight, etc.Output: Feature position array F.Step 1: Texture conversion. For the wind field dataV, the three direction components ri of the XYZ axis are respectively corresponded, and the vector field velocity value is converted into a texture. The mapping function θi is defined as shown in the formula: C=c1,c2,…,c10.Among them,α ∈ (0,255) is the segmentation parameter.Step 2: Fragment processing.(1) Calculate the speed value of the texture fragment in the sample texture in reverse, as the BP neural network input layer data; calculate the hidden layer and output layer value according toC=c1,c2,…,c10, where rc=1/10∑i=110riθc=1/10∑i=110θi is the output value of the layer above the i-th layer.(2) Calculate the error value between the characteristic texture output value and each specified class according torc−ri; among them, rc represents the i-th ideal output of the k-th standard class, and rc represents the actual i-th output.(3) Choose a smallestC=C1,C2,…,Cn. If C=ci,1,ci,2,…,ci,10 is less than the specified error threshold, the fragment is considered to be the specified v-th flow field feature, and the fragment color is set to the specified color Cv; otherwise, the fragment color is set to background color B.Step 3: Save the result.After the vector field is converted into a GPU easy-to-process texture using texture processing, the core code of the relevant feature recognition fragment program is as follows:uniform vec2 tc_off set[243]; //Adjacent grid point texture coordinates.uniform sampler 3D velTex; //vector field texturevoid main(void){vec3 sample[243];for(inti = 0; i＜9 ∗ 9 ∗ 3; i++)sample[i] = texture 3D(velTex, gl_MultiTexCoord0.stp + tc_off set[i]);int res = BP(sample); //BP calculationvec3 color = TF(res); //Calculate the colorgl_FragColor = ve4(color,1.0); }


### 2.2. Multiresolution Rendering

Due to the unevenness of the image, an octree is needed to divide the corresponding space. The specific principle is shown in Figure3.Figure 3
The feature area is located in the octree partition subspace.Voronoi diagram is a space division structure generated based on the principle of nearest neighbors, and its definition is as follows: suppose S is a two-dimensional planeq which is any geometric point on S, and XO=O1,O2,…,On,n≥3 is a set of discrete points on the Euclidean plane. The area VOi is a set that satisfies the following conditions: VOi=p|p∈S,anddp,Oi≤dp,Oj,j≠i,j=1,2,…,n, then VOi is called the Voronoi area associated with the object Oi, and Oi is called the growth object of this area. Let V=VO1,VO2,…,VOn; call V the Voronoi diagram on S generated by O. The Voronoi diagram technology is spatially divided to form a polygon set, where each polygon area corresponds to a point target. The distance from each point of the polygon to the corresponding point target is smaller than other point targets. Voronoi diagram can be divided into vector method and grid method according to its generation method, considering the experimental data as a regular grid structure.A feature-based Voronoi diagram data organization method is proposed based on the grid method in this paper. The steps are as follows:Step 4.
According to the chessboard distance in Figure4, extract the corresponding point distance definition.Figure 4
Checkerboard distance.Step 5.
Local distance propagation calculation is performed for each point target in turn, as shown in the following formula:(1)Di,j=minDi−1,j−1+1,Di−1,j+1,Di−1,j+1+1,Di,j−1+1,Di,j,Di,j+1+1,Di+1,j−1+1,Di+1,j+1,Di+1,j+1+1.
The distance from the surrounding nodes to the point target is calculated. In Algorithm1, D(i, j) represents the distance from the node with the serial number (i, j) to a certain point of target.Step 6.
According to Figure5, organize each feature area and adjacent area nodes into a tree structure, which is called a feature tree. If the superscript is used to indicate the layer number of the node in the feature tree, the subscript indicates the sequence number in the layer, such as Nim represents the i-th node in the m-th layer of the feature tree, then the specific construction process of the feature tree is as follows:Figure 5
Feature tree structure based on Voronoi diagram.Step 7.
Initialize the root node and set the child nodes to be blank.Step 8.
Iterate 2.1–2.3 until all image processing is completed.Step 9.
Create a new nodeNi1. Set the corresponding attributes according to the feature category, and set the node as the i-th child node of R.Step 10.
Create a new nodeposture1t1=p01t1,q01t1,…,qn1t1 and join the feature tree as the first child node of the node posture2t2=p02t2,q02t2,…,qn2t2; obtain the nodes with a distance of D ≤ 1 in the graph as child nodes and join the node q01t1 in turn, as shown in Figure 5.Step 11.
Obtain the neighboring area nodes in the distance graph 2≤D≤3, form groups in turn according to the screening scale factor a, and filter out the parent nodes q02t2 to join the nodes distq01t1,q01t2=logq01t1−1×q02t2 according to the screening rules.
For the node screening rule for neighboring area, the method of taking the minimum dimension node is adopted for the experiment in this paper, that is, if the neighboring area node groupM contains nodes qi1t and qi2t, and ri1t, then ri2t is selected [17, 18]. After adopting the feature tree method, any feature area corresponds to a unique feature tree subnode, which effectively solves the problem of low efficiency in drawing the feature area when the data field is represented by an octree. The generation and extinction of time-varying field features only correspond to the addition and deletion of a single node in the feature tree. The data structure is easier to maintain. After the feature tree and the global octree are generated, the fisheye view technology is used for multiresolution rendering. The fisheye view technology was first proposed by Furnas. Its basic idea is to display fine-grained information to the user’s attention area and coarse-grained information to the background area. After obtaining the feature tree and the global octree, the multiresolution rendering process mainly includes two steps: (1) drawing the nodes in the global octree according to the background field detail control parameter β; (2) drawing the corresponding nodes in the feature tree. In order to ensure the authenticity of the visualization of the data field, the original image data are restored appropriately by keeping the data visualization graphs of the focus area and the background area in the same size.


## 2.1. Feature Detection and Recognition

For feature detection and recognition, real-time interactive processing is required. The feature detection and recognition algorithm is carried out on the CPU, and in order to improve the speed, the method of taking the area of the critical point as the candidate unit is adopted. In order to realize the real-time interaction of feature detection and recognition, this paper designs a GPU-based feature detection and recognition algorithm. The basic idea is to convert the flow field into texture fragment blocks and use the high parallel characteristics and programmability of GPU to convert BP neural network feature recognition into processing texture fragments. The basic flow is shown in Figure2. The algorithm mainly includes the following steps:Figure 2
GPU feature recognition algorithm flow chart.Step 1.
Texture conversion: it is responsible for converting the flow field data into a color texture that is easy to process by the GPU.Step 2.
GPU processing: it is responsible for feature recognition of the area where the current fragment is located.Step 3.
Saving the results: it is responsible for reading the recognition result from the GPU and saving it to the corresponding data structure.
Because the GPU feature recognition process is a parallel process, this paper does not adopt the critical point region candidate unit method, but uses the sequential traversal method. The reason is analyzed as follows: When using the critical point candidate unit method or the traversal method on the GPU, assuming that the texture conversion process time isT¯1 and T1, the GPU processing process is T¯2 and T2, and the result saving process is T¯3 and T3, the entire pipeline The processing time is T¯=T¯1+T¯2+T¯3 and T=T1+T2+T3, respectively. Since Step 1 and Step 3 require the same time for the two methods, namely T¯1=T1, T¯3=T3, the length of the pipeline processing time depends on T¯2 and XT2. The critical point area candidate unit method needs to discriminate the fragment type in the fragment shader. If the current fragment corresponds to the critical point, then the identification judgment is performed; otherwise, the process is skipped. Let's suppose that the fragment shader Fi corresponds to the critical point fragment type, and the processing time is T¯Fi, and the fragment shader Fj corresponds to the non-critical point type, and the processing time is T¯Fj. Although T¯Fi<T¯Fj, because the GPU is a parallel processing process, the shader with the longest calculation time among the fragment shaders constitutes the bottleneck of the recognition algorithm, thus T¯2=maxT¯Fi,T¯Fj=T¯Fj. Similarly, there is T2=maxTFi,0≤i≤N for the traversal method, where TFi is the processing time of each fragment. Because the segment type judgment process sentence is added in the critical point area candidate unit method, there must be T¯2=maxT¯Fi>maxTFi=T2, so T¯>T. That is, in the GPU processing process, the speed of the traversal method is better than that of the critical point candidate unit method. (Algorithm 1).
First, obtain the node data of the characteristic area according to the current texture coordinates, then perform the characteristic recognition calculation, and set the current fragment color according to the recognition result. The implementation code of the feature recognition process is basically the same as that on the CPU, but the texture-speed reverse calculation is first required; in addition, to ensure that the texture conversion process does not lose data accuracy, the texture format uses a 32 bit floating point type.
This paper also uses the pressure field data to correct the characteristics of the BP neural network. In the experiment, the characteristics of cyclones and anticyclones were mainly extracted, and cyclones and anticyclones correspond to the low and high pressure centers in the pressure field, respectively. There is a pressureP≤Pthreshhold for the center of the cyclone and P≥Pthreadhold′ for the center of the anticyclone. For wind field data, assume that the position obtained after feature detection using BP neural network is Pi, and the position obtained after detecting the corresponding pressure field data using the pressure amplitude method is Pi′. If DPi,Pi′≤dthreshhold, the detection result is considered correct and output pi; otherwise, the detection result is considered to be wrong, where dthreshhold is the Euler distance error threshold.Algorithm1: Feature detection and recognition algorithm.
Input: Image input data, network weight, etc.Output: Feature position array F.Step 1: Texture conversion. For the wind field dataV, the three direction components ri of the XYZ axis are respectively corresponded, and the vector field velocity value is converted into a texture. The mapping function θi is defined as shown in the formula: C=c1,c2,…,c10.Among them,α ∈ (0,255) is the segmentation parameter.Step 2: Fragment processing.(1) Calculate the speed value of the texture fragment in the sample texture in reverse, as the BP neural network input layer data; calculate the hidden layer and output layer value according toC=c1,c2,…,c10, where rc=1/10∑i=110riθc=1/10∑i=110θi is the output value of the layer above the i-th layer.(2) Calculate the error value between the characteristic texture output value and each specified class according torc−ri; among them, rc represents the i-th ideal output of the k-th standard class, and rc represents the actual i-th output.(3) Choose a smallestC=C1,C2,…,Cn. If C=ci,1,ci,2,…,ci,10 is less than the specified error threshold, the fragment is considered to be the specified v-th flow field feature, and the fragment color is set to the specified color Cv; otherwise, the fragment color is set to background color B.Step 3: Save the result.After the vector field is converted into a GPU easy-to-process texture using texture processing, the core code of the relevant feature recognition fragment program is as follows:uniform vec2 tc_off set[243]; //Adjacent grid point texture coordinates.uniform sampler 3D velTex; //vector field texturevoid main(void){vec3 sample[243];for(inti = 0; i＜9 ∗ 9 ∗ 3; i++)sample[i] = texture 3D(velTex, gl_MultiTexCoord0.stp + tc_off set[i]);int res = BP(sample); //BP calculationvec3 color = TF(res); //Calculate the colorgl_FragColor = ve4(color,1.0); }


## 2.2. Multiresolution Rendering

Due to the unevenness of the image, an octree is needed to divide the corresponding space. The specific principle is shown in Figure3.Figure 3
The feature area is located in the octree partition subspace.Voronoi diagram is a space division structure generated based on the principle of nearest neighbors, and its definition is as follows: suppose S is a two-dimensional planeq which is any geometric point on S, and XO=O1,O2,…,On,n≥3 is a set of discrete points on the Euclidean plane. The area VOi is a set that satisfies the following conditions: VOi=p|p∈S,anddp,Oi≤dp,Oj,j≠i,j=1,2,…,n, then VOi is called the Voronoi area associated with the object Oi, and Oi is called the growth object of this area. Let V=VO1,VO2,…,VOn; call V the Voronoi diagram on S generated by O. The Voronoi diagram technology is spatially divided to form a polygon set, where each polygon area corresponds to a point target. The distance from each point of the polygon to the corresponding point target is smaller than other point targets. Voronoi diagram can be divided into vector method and grid method according to its generation method, considering the experimental data as a regular grid structure.A feature-based Voronoi diagram data organization method is proposed based on the grid method in this paper. The steps are as follows:Step 4.
According to the chessboard distance in Figure4, extract the corresponding point distance definition.Figure 4
Checkerboard distance.Step 5.
Local distance propagation calculation is performed for each point target in turn, as shown in the following formula:(1)Di,j=minDi−1,j−1+1,Di−1,j+1,Di−1,j+1+1,Di,j−1+1,Di,j,Di,j+1+1,Di+1,j−1+1,Di+1,j+1,Di+1,j+1+1.
The distance from the surrounding nodes to the point target is calculated. In Algorithm1, D(i, j) represents the distance from the node with the serial number (i, j) to a certain point of target.Step 6.
According to Figure5, organize each feature area and adjacent area nodes into a tree structure, which is called a feature tree. If the superscript is used to indicate the layer number of the node in the feature tree, the subscript indicates the sequence number in the layer, such as Nim represents the i-th node in the m-th layer of the feature tree, then the specific construction process of the feature tree is as follows:Figure 5
Feature tree structure based on Voronoi diagram.Step 7.
Initialize the root node and set the child nodes to be blank.Step 8.
Iterate 2.1–2.3 until all image processing is completed.Step 9.
Create a new nodeNi1. Set the corresponding attributes according to the feature category, and set the node as the i-th child node of R.Step 10.
Create a new nodeposture1t1=p01t1,q01t1,…,qn1t1 and join the feature tree as the first child node of the node posture2t2=p02t2,q02t2,…,qn2t2; obtain the nodes with a distance of D ≤ 1 in the graph as child nodes and join the node q01t1 in turn, as shown in Figure 5.Step 11.
Obtain the neighboring area nodes in the distance graph 2≤D≤3, form groups in turn according to the screening scale factor a, and filter out the parent nodes q02t2 to join the nodes distq01t1,q01t2=logq01t1−1×q02t2 according to the screening rules.
For the node screening rule for neighboring area, the method of taking the minimum dimension node is adopted for the experiment in this paper, that is, if the neighboring area node groupM contains nodes qi1t and qi2t, and ri1t, then ri2t is selected [17, 18]. After adopting the feature tree method, any feature area corresponds to a unique feature tree subnode, which effectively solves the problem of low efficiency in drawing the feature area when the data field is represented by an octree. The generation and extinction of time-varying field features only correspond to the addition and deletion of a single node in the feature tree. The data structure is easier to maintain. After the feature tree and the global octree are generated, the fisheye view technology is used for multiresolution rendering. The fisheye view technology was first proposed by Furnas. Its basic idea is to display fine-grained information to the user’s attention area and coarse-grained information to the background area. After obtaining the feature tree and the global octree, the multiresolution rendering process mainly includes two steps: (1) drawing the nodes in the global octree according to the background field detail control parameter β; (2) drawing the corresponding nodes in the feature tree. In order to ensure the authenticity of the visualization of the data field, the original image data are restored appropriately by keeping the data visualization graphs of the focus area and the background area in the same size.


## 3. Visual Aid Training System

For shooting sports training, it can rely on the computer to interpret and read the shooting indication results and meanwhile realize the shooting process backtracking, accounting shooting distribution law, shooting deviation error, and historical shooting data analysis, through the overall analysis of these data; to achieve the assistance of sports training, the specific system block diagram is shown in Figure6.Figure 6
Computer vision aided training system.By setting up multiple 3D camera instruments, the target image can be collected, and the multiple processing results and data of the shooting location can be monitored. After the computerized collection of multiple 3D images, the unified processing of video data is realized, and finally 3D image preprocessing is realized, such as deformation correction, image segmentation, image calculation, target recognition, and orientation determination, which are processed according to the corresponding results before the shooting data are obtained.For the processing of the target image, the recognition, statistics, and analysis of the shooting target can be realized. Meanwhile, the deviation calculation is carried out according to the existing design results, and the corresponding shooting correction is given. Ultimately, the quality and effectiveness of shooting training are improved.


## 4. Shooting Data Processing

For shooting data processing, it mainly analyzes and recognizes the digitalized 3D video image, calculates the corresponding shooting ring number, and realizes the unification and display of the results. 3D video image data processing can be divided into image filtering, geometric correction, image segmentation, calculation processing, data storage, and other steps. The specific processing is shown in Figure7.Figure 7
Principles of shooting data processing.The preprocessing of the image is mainly realized by performing grayscale change and filtering processing on the 3D video image. The purpose of the image grayscale transformation is to reduce the dimension of the image data and improve the processing speed. The role of image filtering is to eliminate noise signals in the image and improve the reliability of subsequent processing. Currently, commonly used image processing algorithms are generally applicable to RGB images. Because the target is simple and the shape is basically fixed in the image processing of this system, the processing is performed on grayscale images.


### 4.1. Image Pretreatment

The 3D video imageI0 (RGB) is transformed into a grayscale image I through preprocessing, specifically as shown in the following formula:(2)Ix,y=0.299I0r,x,y+0.587I0g,x,y+0.114I0b,x,y,where I (x, y) is the gray value of the gray image at the coordinates (x, y), I0r,x,y is the gray value of the red component of the RGB image at the coordinates (x, y), and I0g,x,y is the gray value of the RGB image at the coordinates (x, y). The gray value of the green component at the coordinates (x, y) and I0b,x,y is the gray value of the blue component of the RGB image at the coordinates (x, y).Equation (2) is generally used in grayscale conversion occasions, but for images with different tones and brightness, the 3D images obtained are different. Under certain circumstances, the 3D image features obtained are not necessarily the most prominent. Therefore, for the specific environment generally adopts the weighted summation method to calculate, the specific is shown in the following formula:(3)Ix,y=w1I0t,x,y+w2I0g,x,y+w3I0b,x,yw1+w2+w3=1,where w1, w2, and w3 are the weights of the RGB components of the in-color image, which can be obtained through experiments.Suppose the image weighted filtering template is shown in the following formula:(4)hl=w11⋯w1l⋮⋮wl1⋯wll∑wij=1.The image weighted filtering algorithm is shown in the following formula:(5)I′x,y=∑∑wx+i,y+iIx+i,y+j,i∈−l−12,l−12,i∈Z,j∈−l−12,l−12,j∈Z,where wx+i,y+j is the element value of hl at (x + i, y + j).Generally,wij=wi′j′=1/l2 is taken.


### 4.2. Image Segmentation and Model Correction

On the basis of target determination, the single-threshold segmentation method is used for calculation, and the specifics are shown in the following formula:(6)Bx,y=1,I′x,y≥Gth,0,I′x,y<Gth,where B is the segmented binary image, and Gth is the segmentation threshold.In order to further quantitatively calculate and simulate and realize the effectiveness of recognition and calculation, the model replaces the 3D video image and integrates the actual target, which can not only reduce the workload but also reduce the corresponding interference, to achieve accurate simulation. Therefore, the image needs to be segmented first to realize the effective identification of the position of the ring.By dividing the target, the center and the ring are distinguished, and the side view of the target is realized, and the distortion of the 3D video image is further corrected to obtain a more accurate target model. The specific imaging diagram is shown in Figure8.Figure 8
Schematic diagram of imaging.Due to the deviation between the position of the equipment and the target, there will be deformation when acquiring images. As shown in Figure8, the larger the included angle θ, the greater the deformation may be.


## 4.1. Image Pretreatment

The 3D video imageI0 (RGB) is transformed into a grayscale image I through preprocessing, specifically as shown in the following formula:(2)Ix,y=0.299I0r,x,y+0.587I0g,x,y+0.114I0b,x,y,where I (x, y) is the gray value of the gray image at the coordinates (x, y), I0r,x,y is the gray value of the red component of the RGB image at the coordinates (x, y), and I0g,x,y is the gray value of the RGB image at the coordinates (x, y). The gray value of the green component at the coordinates (x, y) and I0b,x,y is the gray value of the blue component of the RGB image at the coordinates (x, y).Equation (2) is generally used in grayscale conversion occasions, but for images with different tones and brightness, the 3D images obtained are different. Under certain circumstances, the 3D image features obtained are not necessarily the most prominent. Therefore, for the specific environment generally adopts the weighted summation method to calculate, the specific is shown in the following formula:(3)Ix,y=w1I0t,x,y+w2I0g,x,y+w3I0b,x,yw1+w2+w3=1,where w1, w2, and w3 are the weights of the RGB components of the in-color image, which can be obtained through experiments.Suppose the image weighted filtering template is shown in the following formula:(4)hl=w11⋯w1l⋮⋮wl1⋯wll∑wij=1.The image weighted filtering algorithm is shown in the following formula:(5)I′x,y=∑∑wx+i,y+iIx+i,y+j,i∈−l−12,l−12,i∈Z,j∈−l−12,l−12,j∈Z,where wx+i,y+j is the element value of hl at (x + i, y + j).Generally,wij=wi′j′=1/l2 is taken.


## 4.2. Image Segmentation and Model Correction

On the basis of target determination, the single-threshold segmentation method is used for calculation, and the specifics are shown in the following formula:(6)Bx,y=1,I′x,y≥Gth,0,I′x,y<Gth,where B is the segmented binary image, and Gth is the segmentation threshold.In order to further quantitatively calculate and simulate and realize the effectiveness of recognition and calculation, the model replaces the 3D video image and integrates the actual target, which can not only reduce the workload but also reduce the corresponding interference, to achieve accurate simulation. Therefore, the image needs to be segmented first to realize the effective identification of the position of the ring.By dividing the target, the center and the ring are distinguished, and the side view of the target is realized, and the distortion of the 3D video image is further corrected to obtain a more accurate target model. The specific imaging diagram is shown in Figure8.Figure 8
Schematic diagram of imaging.Due to the deviation between the position of the equipment and the target, there will be deformation when acquiring images. As shown in Figure8, the larger the included angle θ, the greater the deformation may be.


## 5. Computer-Aided Training Simulation Experiment

After obtaining the data of a single shot, the system can perform microanalysis according to the characteristics of the detected data and after a complete shot, according to the data distribution characteristics and changes. Meanwhile, it can also analyze the changes in athlete performance based on the changes in historical data, evaluate the training effects of athletes and coaches, and give reference training programs based on the characteristics of the data.


### 5.1. Single Data-Aided Training

The calculation of the correction variable is to give a reference correction value based on the current deviation and the last deviation. Under normal circumstances, it is considered that this shooting is the result of the correction based on the previous shooting data, and the correction deviation is calculated according to the theoretical correction situation, and then the appropriate correction plan is estimated based on the current deviation. If the current shot is the first shot, the data at the center position is the last shot data, and the deviation caused after correction is not made. After the callback, there is a system deviation, so the last excellent target is the adjustment point for adjustment.


### 5.2. Complete Process Assisted Training

Suppose that 10 shots have been completed, and the ordered dataset isC=c1,c2,…,c10; according to the dataset, the ring count change, the position change, the bullet point distribution, the data validity change, and the bullet point system deviation are macro analyzed, and the bullet point dispersion assessment is performed.


#### 5.2.1. Analysis of Data Changes

Through the bullet point data curve, you can observe the changes in the athlete’s performance during the entire shooting process, analyze the best and worst points of the state, and provide a reference for the adjustment of the state during the shooting process.


#### 5.2.2. Analysis of Data Statistical Characteristics

The analysis of the statistical characteristics of the data includes the deviation data of the average center point, the statistical circle coordinates and radius of the bullet point set, the dispersion of the data, and the credibility of each data. The center point deviation data is calculated using the mean value of the number of points.(7)rc=110∑i=110ri,θc=110∑i=110θi.The credibility of the data can be calculated by using the distance between the shooting point and the statistical circle. The closer to the statistical circle center, the higher the credibility of the shot. The farther away from the statistical circle, the lower the credibility of the shot.


#### 5.2.3. Auxiliary Training

The auxiliary training content mainly includes the following:(i)
Correction of system deviation(ii)
Psychological adjustment reference during shooting(iii)
Reference for posture adjustment during shooting(iv)
Reference for breathing adjustment during shooting(v)
Suggestions for further improving performance


### 5.3. Tracking and Evaluation of Training Process

The evaluation of the training process includes two levels: athletes and coaches.


#### 5.3.1. Evaluation of the Athlete’s Training Process

If the long-term training does not significantly reduce this deviation and shows regular changes, the training method or the athlete’s suitability for the sport should be reassessed.


#### 5.3.2. Evaluation of Coach Training Process

In addition to the abovementioned auxiliary training, the auxiliary training system can also statistically analyze the impact of shooting performance and other factors.(i)
The correlation between shooting performance and shooting environment temperature(ii)
The correlation between shooting performance and sunny and rainy(iii)
Changes in shooting performance with the four seasons(iv)
The correlation between shooting performance and time period, etc.The specific analysis of the abovementioned situation provides a reference for enhancing the strengths and avoiding weaknesses and strengthening the training purposefully.Assuming that there are two action segments mi(t) and m2(t), they can be connected into a new action sequence using motion mirroring and motion transition techniques. The last posture of m1(t) and the first posture of m2(t) are set to be(8)posture1t1=p01t1,q01t1,…,qn1t1,posture2t2=p02t2,q02t2,…,qn2t2.According to the differencedistq01t1,q01t2=logq01t1−1×q02t2 between the direction Z2−Z1 and Z3−Z1 to determine whether the action m2(t) needs to be mirrored, the result may still be recorded as m2(t).Assume that the long side of the virtual trampoline is consistent with theX direction, and the wide side is consistent with the Z direction, and its coordinate system OXY B is defined according to the right-handed system, and its initial position is assumed to coincide with the global coordinate system. Select the three vertices p1u1,v1, p2u2,v2, and p3u3,v3 of the trampoline from the training video, and assume that p2 is the common vertices of the long side and the wide side, then the camera orthogonal projection model is shown in the following formula:(9)uv=s100010XYZ.The pointsp1, p2, and p3 can be mapped to the points P1 (X1, Y1, Z1), P2 (X2, Y2, Z2), and P3 (X3, Y3, Z3) in the three-dimensional space, respectively, and the relative depth Z2−Z1 and Z3−Z1 can be solved, so as to determine the location of the virtual shooting. The simulation experiment results show that the computer digitized 3D video image is effective.


## 5.1. Single Data-Aided Training

The calculation of the correction variable is to give a reference correction value based on the current deviation and the last deviation. Under normal circumstances, it is considered that this shooting is the result of the correction based on the previous shooting data, and the correction deviation is calculated according to the theoretical correction situation, and then the appropriate correction plan is estimated based on the current deviation. If the current shot is the first shot, the data at the center position is the last shot data, and the deviation caused after correction is not made. After the callback, there is a system deviation, so the last excellent target is the adjustment point for adjustment.


## 5.2. Complete Process Assisted Training

Suppose that 10 shots have been completed, and the ordered dataset isC=c1,c2,…,c10; according to the dataset, the ring count change, the position change, the bullet point distribution, the data validity change, and the bullet point system deviation are macro analyzed, and the bullet point dispersion assessment is performed.


### 5.2.1. Analysis of Data Changes

Through the bullet point data curve, you can observe the changes in the athlete’s performance during the entire shooting process, analyze the best and worst points of the state, and provide a reference for the adjustment of the state during the shooting process.


### 5.2.2. Analysis of Data Statistical Characteristics

The analysis of the statistical characteristics of the data includes the deviation data of the average center point, the statistical circle coordinates and radius of the bullet point set, the dispersion of the data, and the credibility of each data. The center point deviation data is calculated using the mean value of the number of points.(7)rc=110∑i=110ri,θc=110∑i=110θi.The credibility of the data can be calculated by using the distance between the shooting point and the statistical circle. The closer to the statistical circle center, the higher the credibility of the shot. The farther away from the statistical circle, the lower the credibility of the shot.


### 5.2.3. Auxiliary Training

The auxiliary training content mainly includes the following:(i)
Correction of system deviation(ii)
Psychological adjustment reference during shooting(iii)
Reference for posture adjustment during shooting(iv)
Reference for breathing adjustment during shooting(v)
Suggestions for further improving performance


## 5.2.1. Analysis of Data Changes

Through the bullet point data curve, you can observe the changes in the athlete’s performance during the entire shooting process, analyze the best and worst points of the state, and provide a reference for the adjustment of the state during the shooting process.


## 5.2.2. Analysis of Data Statistical Characteristics

The analysis of the statistical characteristics of the data includes the deviation data of the average center point, the statistical circle coordinates and radius of the bullet point set, the dispersion of the data, and the credibility of each data. The center point deviation data is calculated using the mean value of the number of points.(7)rc=110∑i=110ri,θc=110∑i=110θi.The credibility of the data can be calculated by using the distance between the shooting point and the statistical circle. The closer to the statistical circle center, the higher the credibility of the shot. The farther away from the statistical circle, the lower the credibility of the shot.


## 5.2.3. Auxiliary Training

The auxiliary training content mainly includes the following:(i)
Correction of system deviation(ii)
Psychological adjustment reference during shooting(iii)
Reference for posture adjustment during shooting(iv)
Reference for breathing adjustment during shooting(v)
Suggestions for further improving performance


## 5.3. Tracking and Evaluation of Training Process

The evaluation of the training process includes two levels: athletes and coaches.


### 5.3.1. Evaluation of the Athlete’s Training Process

If the long-term training does not significantly reduce this deviation and shows regular changes, the training method or the athlete’s suitability for the sport should be reassessed.


### 5.3.2. Evaluation of Coach Training Process

In addition to the abovementioned auxiliary training, the auxiliary training system can also statistically analyze the impact of shooting performance and other factors.(i)
The correlation between shooting performance and shooting environment temperature(ii)
The correlation between shooting performance and sunny and rainy(iii)
Changes in shooting performance with the four seasons(iv)
The correlation between shooting performance and time period, etc.The specific analysis of the abovementioned situation provides a reference for enhancing the strengths and avoiding weaknesses and strengthening the training purposefully.Assuming that there are two action segments mi(t) and m2(t), they can be connected into a new action sequence using motion mirroring and motion transition techniques. The last posture of m1(t) and the first posture of m2(t) are set to be(8)posture1t1=p01t1,q01t1,…,qn1t1,posture2t2=p02t2,q02t2,…,qn2t2.According to the differencedistq01t1,q01t2=logq01t1−1×q02t2 between the direction Z2−Z1 and Z3−Z1 to determine whether the action m2(t) needs to be mirrored, the result may still be recorded as m2(t).Assume that the long side of the virtual trampoline is consistent with theX direction, and the wide side is consistent with the Z direction, and its coordinate system OXY B is defined according to the right-handed system, and its initial position is assumed to coincide with the global coordinate system. Select the three vertices p1u1,v1, p2u2,v2, and p3u3,v3 of the trampoline from the training video, and assume that p2 is the common vertices of the long side and the wide side, then the camera orthogonal projection model is shown in the following formula:(9)uv=s100010XYZ.The pointsp1, p2, and p3 can be mapped to the points P1 (X1, Y1, Z1), P2 (X2, Y2, Z2), and P3 (X3, Y3, Z3) in the three-dimensional space, respectively, and the relative depth Z2−Z1 and Z3−Z1 can be solved, so as to determine the location of the virtual shooting. The simulation experiment results show that the computer digitized 3D video image is effective.


## 5.3.1. Evaluation of the Athlete’s Training Process

If the long-term training does not significantly reduce this deviation and shows regular changes, the training method or the athlete’s suitability for the sport should be reassessed.


## 5.3.2. Evaluation of Coach Training Process

In addition to the abovementioned auxiliary training, the auxiliary training system can also statistically analyze the impact of shooting performance and other factors.(i)
The correlation between shooting performance and shooting environment temperature(ii)
The correlation between shooting performance and sunny and rainy(iii)
Changes in shooting performance with the four seasons(iv)
The correlation between shooting performance and time period, etc.The specific analysis of the abovementioned situation provides a reference for enhancing the strengths and avoiding weaknesses and strengthening the training purposefully.Assuming that there are two action segments mi(t) and m2(t), they can be connected into a new action sequence using motion mirroring and motion transition techniques. The last posture of m1(t) and the first posture of m2(t) are set to be(8)posture1t1=p01t1,q01t1,…,qn1t1,posture2t2=p02t2,q02t2,…,qn2t2.According to the differencedistq01t1,q01t2=logq01t1−1×q02t2 between the direction Z2−Z1 and Z3−Z1 to determine whether the action m2(t) needs to be mirrored, the result may still be recorded as m2(t).Assume that the long side of the virtual trampoline is consistent with theX direction, and the wide side is consistent with the Z direction, and its coordinate system OXY B is defined according to the right-handed system, and its initial position is assumed to coincide with the global coordinate system. Select the three vertices p1u1,v1, p2u2,v2, and p3u3,v3 of the trampoline from the training video, and assume that p2 is the common vertices of the long side and the wide side, then the camera orthogonal projection model is shown in the following formula:(9)uv=s100010XYZ.The pointsp1, p2, and p3 can be mapped to the points P1 (X1, Y1, Z1), P2 (X2, Y2, Z2), and P3 (X3, Y3, Z3) in the three-dimensional space, respectively, and the relative depth Z2−Z1 and Z3−Z1 can be solved, so as to determine the location of the virtual shooting. The simulation experiment results show that the computer digitized 3D video image is effective.


## 6. Conclusions

Physical training is an important way and means to improve sports performance. Therefore, reasonable and effective physical training is extremely important. Relying on computer digital 3D video image processing, the assistance of the design training system is achieved through combing the shooting process. Through the analysis of 3D image processing, data processing and analysis under different series of shooting levels are realized to achieve data statistics and mining and finally provide support for sports training. The results of the simulation experiment show that the computer digitized 3D video images are effective and can support sports-assisted training.



---

*Source: 2105790-2022-01-10.xml*