# Design of Precise Estimation Algorithm for Sports Pose Based on Multifeature Fusion

Authors: Gao Songai; Zhang Ningning; Zheng Xiang
Journal: Mobile Information Systems
 (2022)

---

## Abstract

The estimation time and estimation precision of motion pose samples are problematic for the pose estimation algorithm of sports movements. This paper proposes a multifeature fusion-based algorithm for accurate posture estimation. The human rod model is constructed after analyzing the human pose estimation technology. Using the Kalman filter method, the degree of freedom and range of motion of the major joints of the human body were determined. The eight-star model was used to extract the sports posture features, and the weighted average method was used to process the grayscale images of sports. Using the multifeature fusion method, the extracted multisource feature vector information is thoroughly analyzed and processed, and a new group of fusion feature vectors is created. Using a mixture Gaussian distribution model, the posture estimation of an athlete’s body is accomplished. Experimental results indicate that when the amount of sports pose sample data is 900 GB, the accurate estimation time of the proposed method is 5.3 s, and its accuracy is 100 percent. Improve the estimation accuracy of samples of sports posture.


---

## Body

## 1. Introduction

With the advancement of society and the improvement of living standards, people are beginning to pay more attention to physical health and other aspects and are constantly improving their physical quality through sports. In the process of movement, standard movement posture can not only determine the effect of movement to a certain extent but also provide the greatest possible protection against injury [1]. In badminton, for instance, the use of standard badminton movements for daily exercise can not only effectively exercise the body but also improve a player's level of competition [2]. However, the definition of standard motor posture during movement is primarily based on images or verbal instruction, resulting in the absence of quantitative evaluation criteria for standard motor posture. Nevertheless, the estimation of human pose and recognition of human motion from the collected image or video sequence provide a theoretical foundation for measuring standard motion pose. Early human action recognition requires external equipment to perceive the change in human posture in order to recognize human action. With the development of machine learning and deep learning, there are numerous research directions in the academic field, including image processing, SVM classifier, and deep neural network [3]. These advancements have made it possible for computers to detect human motion using only devices like cameras, thereby drastically reducing the number of external perceptrons.Human pose estimation is the most critical technology in the process of sports pose recognition. In order to guide people’s posture during the movement process, a series of features derived from human pose estimation are compared with some standard actions. Consequently, the human motion assessment system has evolved significantly. Simultaneously, people’s focus is gradually shifting from the professional type to the portable type: a portable human motion assessment system that is unaffected by location and environment can enhance the quality of people’s daily exercise. Human pose recognition has garnered the interest of numerous domestic and international research institutes and university-level laboratories [4] due to its high theoretical value and vast application potential. A number of algorithms for human pose recognition and related specific problems have been proposed domestically and internationally. In addition, advances have been made in the areas of human object recognition, keyframe extraction, human pose estimation, and other specialized issues. The traditional human gesture recognition technology has low recognition rate, high difficulty, high cost, slow speed, and large storage space for human behavior recognition, which is not the requirement for the application of human-computer interaction technology in the future. Combining the keyframe selection technique with the human pose estimation technique drastically reduces computation and saves storage space. In recent years, the research status of human pose estimation at home and abroad is summarized as follows [5].The authors of [6] proposed a 3D human pose estimation method based on semisupervised learning convolutional neural network, collected video data containing motion pose, built a semisupervised learning convolutional neural network model, and estimated the corresponding 3D pose skeleton changing pose. This method can accurately estimate the 3D pose of the corresponding frame in a video, but sports pose estimation efficiency is low. The method proposed by the authors of [7] for acquiring specific motion frames in motion videos is based on human pose estimation and clustering. Select the HRNet pose estimation model as the basis, establish the model lightweight process and combine with the DARK data encoding, construct the Small-HRNet network model, extract from the video human body joint points, and use the human skeleton features of each frame of video as sample points. The standard skeleton feature of the motion frame is the cluster center, which clusters the entire video. The video’s specific motion frame is obtained through clustering, and the experiment is conducted on the motion dataset. This method can reduce the number of body pose estimation model parameters while maintaining accuracy, but the classification effect of the proposed method is inadequate. A human pose estimation optimization algorithm for deformable models is proposed with reference to [8]. The model parameters are extracted from each frame of a color image by a neural network, and the parameters are optimized and solved using human key points and contour as constraints. According to the interframe coherence of the video sequence, the error of the pose estimation results of all video frames is then corrected, making the motion sequence more fluid and smooth. The point cloud obtained from the depth map and the corresponding color map model were used as the joint input, followed by the use of the distance constraint between the point cloud and the corresponding point of the model to optimize and solve the problem. Eventually, a result resembling the actual human pose was achieved. Using point cloud data optimization, the algorithm can effectively correct the single-frame pose estimation results and greatly improve its accuracy. Nonetheless, the accuracy of posture estimation in sports is poor.This paper proposes an accurate estimation algorithm for sports pose based on the fusion of multiple features.


## 2. Acquisition of Sports Posture Features

### 2.1. Human Pose Estimation Techniques

According to the Oxford Dictionary, a posture is a particular posture of the body, and the way a person maintains his or her body. The human pose estimation method is to extract, classify, distinguish, and describe the pose features of the human body, which has been widely paid attention to in recent years. The application of human physiology, digital image processing, pattern recognition, and other disciplines is an interdisciplinary research field [9]. From the perspective of the collection method of human pose information, human pose estimation technology can be divided into two categories: contact recognition and noncontact recognition [10]. According to different classification criteria, it can be divided into a variety of human pose recognition algorithms. From the perspective of implementation methods, it is usually divided into three types: 3D model reconstruction method, which extracts 3D features from effective samples to construct 3D models, human appearance model method, by obtaining the shape features of the human body to establish a two-dimensional model and use the model matching method to complete the recognition, motion model method, according to the motion characteristics of the classification. From the perspective of pattern recognition, attitude recognition is a problem of classifying time-varying features, that is, matching the test sequence with the preset sequence according to the obtained feature information.


### 2.2. Construction of Human Rod Model

According to the knowledge of human anatomy, there are more than 200 bones in the human body, and the bones are connected by joints, which makes the whole system extremely complex, and the actual data acquisition is very complicated. Due to the particularity of the acquisition unit used in this system, we could not collect all the data for each bone [11].At present, most of the studies on the human body are simplified modeling and analysis with the rod model; that is, the main joint of the human body is abstracted into a point, and the limbs between the joints are abstracted into a link, mainly considering elbow, shoulder, hip, hip, knee, ankle, and other joints. In this way, the whole human body is simplified into a rod-like model, which is also the premise for the analysis of human posture, as shown in Figure1.Figure 1
Schematic diagram of rod mannequin.After the simplified rod human body model has been established, the range of motion of each joint must be constrained in order to bring the model closer to the normal human posture. In rehabilitation medicine, the initial posture of human movement is defined as follows: the body is upright and facing forward, the eyes are flat, the feet are standing together, the toes are pointing forward, the upper limbs are hanging at the side of the body, and the palm is attached to the outer thigh of the side of the body. The joint has three components: the articular surface, the articular capsule, and the articular cavity [12]. The shape and structure of the articular surface determine the axis of joint activity, and the freedom of joint movement is closely associated with the movement around the axis. There are several degrees of freedom if there are multiple directions of joint activity. All joints with more than two degrees of freedom are capable of producing circumferential motion, and the limb typically rotates about the joint axis. The degrees of freedom and motion range of each major human joint are as follows:(1)
Shoulder joint: there are three degrees of freedom, and the upper arm sags in a neutral position. Its range of motion is forward flexion: 70∼90°, posterior extension: 40∼45°, forward flexion: 150∼170°, upward lift: 160∼180°, abduction: 80∼90°, adduction: 20∼40°, internal rotation: 70∼90°, external rotation: 40∼50°.(2)
Elbow joint: there is one degree of freedom. The neutral position of the elbow joint is forearm extension. Flexion: 135∼150°, overextension: 10°, pronation: 80∼90°, supination: 80∼90°.(3)
Wrist joint: there are two degrees of freedom. The neutral position is the hand in line with the forearm, and the palm is downward. Dorsal extension: 80∼90°, palmar flexion: 50∼60°, radial tilt: 25∼30°, ulnar tilt: 30∼40°.(4)
Hip joint: there are three degrees of freedom. The neutral position is hip joint extension and patella forward. Flexion can reach 130∼140°, posterior extension can reach 10∼15°, abduction can reach 30∼45°, adduction can reach 20∼30°, knee flexion 90°, normal external rotation 30∼40°, internal rotation 40∼50°.(5)
Knee joint: there are two degrees of freedom. The neutral position is knee joint extension, flexion: 120∼150°, hyperextension: 5∼10°, internal rotation about 10°, and external rotation: 20° when bending the knee.(6)
Ankle joint: there are two degrees of freedom. The neutral position is a 90 degree Angle between the foot and the calf, without varus or valgus. Dorsiflexion is about 20∼30°, plantar flexion is about 40∼50°, varus is 30°, and valgus is 30∼35°.In this paper, attitude solution refers to the use of existing nine-axis sensor data and relevant theories and algorithms to get attitude Angle information from these data. Attitude solving algorithm belongs to the category of data fusion, that is, attitude solving by fusing multiple sensor data [13]. The commonly used algorithms for attitude resolution include Kalman filter, extended Kalman filter, complementary filter, and gradient descent algorithm.The model-based linear minimum variance estimation of Kalman filter is a very widely used attitude calculation algorithm. As an important optimal estimation theory, it is widely used in various fields such as automatic control and aerospace [14]. Since discrete Kalman filter algorithm can be directly implemented on computer to estimate multidimensional stochastic processes, it will be introduced in the following.Let the system noise sequenceWk−1 act on the estimated state Xk at time tk and can be described as follows:(1)Xk=Φk,k−1Xk−1+Γk−1Wk−1.The measurement equation is a linear function, as shown in the following equation:(2)Zk=HkXk+Vk,where Φk,k−1is the one-step state transition matrix from time tk−1 to time tk, Γk−1 is the system noise control matrix, Hk is the observation matrix, Vk is the measurement noise sequence, and wk is the disturbance noise vector of the system [15].Assuming that the system noise variance matrixQk is nonnegative definite and the observation noise variance matrix Rk is positive definite, the estimator can be solved in the following way:Single-step state prediction equation is as follows:(3)xk.k−1′=Φk,k−1Xk−1′.State estimation equation is as follows:(4)xk′=xk.k−1′+KkZk−HkXk,k−1′.Mean square error estimation matrix is as follows:(5)Pk=I−KkHkPk,k−1,where the initial values P0 and x0′ need to be given in advance, so that the state estimate xk′ at time k can be obtained from the measurement equation Zk at time k recursively.


### 2.3. Attitude Feature Extraction Based on Eight-Star Model

The image of the human object is a two-dimensional array. The extracted features of the human object are transformed into a one-dimensional feature vector, which greatly reduces the feature vector’s dimension and computational complexity. This article focuses on walking, running, and jumping in three postures for feature extraction and recognition, through the analysis of three kinds of attitude, and the legs move in arm swinging amplitude, arising in the process of adopting the star model to the human body modeling, namely through the calculation of center of mass of human movement target image points and local contour pole and its mutual relations, in order to obtain the attitude. In this paper, an eight-star model is utilized to extract a set of feature vectors containing seventeen feature values to characterize the human pose in each frame [16]. Compared to the six-star model, it can describe the human body’s contours, including the head and limbs, more precisely, and it is not redundant.Eight-star model is an improvement of the traditional star model, which is a feature description model formed according to the eight local contour poles and centroid points of the human object. The human pose model is established as shown in the following equation:(6)Fet=D,A,e,where D is the distance between the eight local contour poles of the eight-star model and the centroid of the human target. A is the minimum Angle formed by the straight line and horizontal line from the eight local contour poles to the centroid point of the eight-star model; e is the eccentricity of the human target. The specific feature extraction process is as follows:Calculate the centroid of human objectxc,yc:(7)xc=1Nb∑i=1Nbxi,yc=1Nb∑i=1Nbyi,where Nb is the number of target pixels.The target is divided into four parts by horizontal lines and vertical lines passing through the center of mass. Based on the overall contour of the target, eight local contour pole coordinatesxi,yi of the rightmost, uppermost, leftmost, and bottommost contour in each part are extracted, respectively.The Euclidean distances between the eight local contour polesxi,yi and the centroid point xc,yc are calculated, which are, respectively, expressed as di, and the calculation formula is shown in the following equation:(8)di=xi−xc2+yi−yc2,i=1,2,⋯,8.Then,D=d1,d2,⋯,d8. The extraction process is shown in Figure 2.Figure 2
Feature extraction process of pole distance of local contour of human eight-star model.After feature extraction and fusion, the feature vector of the human pose in the video at framet is shown in the following equation:(9)Fett=d1,d2,⋯,d8,a8,e.The feature vector is a 17-dimensional feature vector constructed from the local contour pole distance feature, local contour pole Angle feature, and eccentricity feature of the human eight-star model, wheredi is the pole distance feature of the eight-star local contour, ai is the Angle feature of the human local contour pole, and e is the eccentricity of the human body. The computation complexity and memory of the gesture description operator are relatively low, so the system can achieve a good real-time recognition effect.In practice, the video image will contain some noise and affect the effect of detection and recognition. Therefore, the primary task of video image processing is image preprocessing, and its effect directly affects the subsequent detection and recognition.


### 2.4. Sports Color Image Grayscale

The color value of the true color image pixel is composed of three components: red R, green G, and blue B. The value of each component has 256 levels of brightness, which is represented by numbers from 0 to 255. Therefore, there are 224 colors, which is relatively complicated to process directly. The gray image is a special color image, whose R, G, and B components are all equal, so there are only 28 colors left, which is convenient for image processing [17]. Therefore, in order to realize the real-time performance of the algorithm, the original image is converted into gray image before the video image processing, so as to reduce the computation amount of the algorithm. In general, the gray-level methods used include component method, maximum method, average method, and weighted average method.In the component method, one of the three components of RGB is used as the gray value of the pixel, as shown in the following equation:(10)I=RorI=GorI=B,where I represents the pixel brightness of gray image, and R, G, and B, respectively, represent the brightness values of RGB components. A gray image can be selected according to the application needs. Taking the 7th frame of LENA_WALK1 video, the 10th frame of PETS2006 video set, and the 21st frame of test video set in Weizmann database as examples, the grayscale is realized by component method.The above method is simple to implement but lacks consideration of the importance of different components. According to the importance of components and other needs of images, the weighted average of RGB three components can be taken as the gray value of pixels [18]. From the perspective of human physiology, the human eye is sensitive to green and insensitive to blue. Therefore, according to the visual model of human eyes, the weight of G component is set to be large, while the weight of B component is small, and its commonly used formula is shown in the following equation:(11)I=0.3R+0.5G+0.11B.The gray weights adopted by OpenCV open source computer vision library also follow the visual model of human eyes, and the weighted average formula is shown in the following equation:(12)I=0.212671R+0.715160G+0.72169B.From the perspective of theory and practice, the weighted average method is more reasonable and more in line with the needs of practical application, so the weighted average method is the most commonly used method of image graying. In this paper, the weighted average method is used to gray-scale the image [19], and the cvCvtColor function in the OpenCV computer vision open source library is used to complete the gray-scale of the image by setting its third parameter as CV_BGR2GRAY. Or using cvLoadImage function, by setting its second parameter to 0 can realize the color image gray.


## 2.1. Human Pose Estimation Techniques

According to the Oxford Dictionary, a posture is a particular posture of the body, and the way a person maintains his or her body. The human pose estimation method is to extract, classify, distinguish, and describe the pose features of the human body, which has been widely paid attention to in recent years. The application of human physiology, digital image processing, pattern recognition, and other disciplines is an interdisciplinary research field [9]. From the perspective of the collection method of human pose information, human pose estimation technology can be divided into two categories: contact recognition and noncontact recognition [10]. According to different classification criteria, it can be divided into a variety of human pose recognition algorithms. From the perspective of implementation methods, it is usually divided into three types: 3D model reconstruction method, which extracts 3D features from effective samples to construct 3D models, human appearance model method, by obtaining the shape features of the human body to establish a two-dimensional model and use the model matching method to complete the recognition, motion model method, according to the motion characteristics of the classification. From the perspective of pattern recognition, attitude recognition is a problem of classifying time-varying features, that is, matching the test sequence with the preset sequence according to the obtained feature information.


## 2.2. Construction of Human Rod Model

According to the knowledge of human anatomy, there are more than 200 bones in the human body, and the bones are connected by joints, which makes the whole system extremely complex, and the actual data acquisition is very complicated. Due to the particularity of the acquisition unit used in this system, we could not collect all the data for each bone [11].At present, most of the studies on the human body are simplified modeling and analysis with the rod model; that is, the main joint of the human body is abstracted into a point, and the limbs between the joints are abstracted into a link, mainly considering elbow, shoulder, hip, hip, knee, ankle, and other joints. In this way, the whole human body is simplified into a rod-like model, which is also the premise for the analysis of human posture, as shown in Figure1.Figure 1
Schematic diagram of rod mannequin.After the simplified rod human body model has been established, the range of motion of each joint must be constrained in order to bring the model closer to the normal human posture. In rehabilitation medicine, the initial posture of human movement is defined as follows: the body is upright and facing forward, the eyes are flat, the feet are standing together, the toes are pointing forward, the upper limbs are hanging at the side of the body, and the palm is attached to the outer thigh of the side of the body. The joint has three components: the articular surface, the articular capsule, and the articular cavity [12]. The shape and structure of the articular surface determine the axis of joint activity, and the freedom of joint movement is closely associated with the movement around the axis. There are several degrees of freedom if there are multiple directions of joint activity. All joints with more than two degrees of freedom are capable of producing circumferential motion, and the limb typically rotates about the joint axis. The degrees of freedom and motion range of each major human joint are as follows:(1)
Shoulder joint: there are three degrees of freedom, and the upper arm sags in a neutral position. Its range of motion is forward flexion: 70∼90°, posterior extension: 40∼45°, forward flexion: 150∼170°, upward lift: 160∼180°, abduction: 80∼90°, adduction: 20∼40°, internal rotation: 70∼90°, external rotation: 40∼50°.(2)
Elbow joint: there is one degree of freedom. The neutral position of the elbow joint is forearm extension. Flexion: 135∼150°, overextension: 10°, pronation: 80∼90°, supination: 80∼90°.(3)
Wrist joint: there are two degrees of freedom. The neutral position is the hand in line with the forearm, and the palm is downward. Dorsal extension: 80∼90°, palmar flexion: 50∼60°, radial tilt: 25∼30°, ulnar tilt: 30∼40°.(4)
Hip joint: there are three degrees of freedom. The neutral position is hip joint extension and patella forward. Flexion can reach 130∼140°, posterior extension can reach 10∼15°, abduction can reach 30∼45°, adduction can reach 20∼30°, knee flexion 90°, normal external rotation 30∼40°, internal rotation 40∼50°.(5)
Knee joint: there are two degrees of freedom. The neutral position is knee joint extension, flexion: 120∼150°, hyperextension: 5∼10°, internal rotation about 10°, and external rotation: 20° when bending the knee.(6)
Ankle joint: there are two degrees of freedom. The neutral position is a 90 degree Angle between the foot and the calf, without varus or valgus. Dorsiflexion is about 20∼30°, plantar flexion is about 40∼50°, varus is 30°, and valgus is 30∼35°.In this paper, attitude solution refers to the use of existing nine-axis sensor data and relevant theories and algorithms to get attitude Angle information from these data. Attitude solving algorithm belongs to the category of data fusion, that is, attitude solving by fusing multiple sensor data [13]. The commonly used algorithms for attitude resolution include Kalman filter, extended Kalman filter, complementary filter, and gradient descent algorithm.The model-based linear minimum variance estimation of Kalman filter is a very widely used attitude calculation algorithm. As an important optimal estimation theory, it is widely used in various fields such as automatic control and aerospace [14]. Since discrete Kalman filter algorithm can be directly implemented on computer to estimate multidimensional stochastic processes, it will be introduced in the following.Let the system noise sequenceWk−1 act on the estimated state Xk at time tk and can be described as follows:(1)Xk=Φk,k−1Xk−1+Γk−1Wk−1.The measurement equation is a linear function, as shown in the following equation:(2)Zk=HkXk+Vk,where Φk,k−1is the one-step state transition matrix from time tk−1 to time tk, Γk−1 is the system noise control matrix, Hk is the observation matrix, Vk is the measurement noise sequence, and wk is the disturbance noise vector of the system [15].Assuming that the system noise variance matrixQk is nonnegative definite and the observation noise variance matrix Rk is positive definite, the estimator can be solved in the following way:Single-step state prediction equation is as follows:(3)xk.k−1′=Φk,k−1Xk−1′.State estimation equation is as follows:(4)xk′=xk.k−1′+KkZk−HkXk,k−1′.Mean square error estimation matrix is as follows:(5)Pk=I−KkHkPk,k−1,where the initial values P0 and x0′ need to be given in advance, so that the state estimate xk′ at time k can be obtained from the measurement equation Zk at time k recursively.


## 2.3. Attitude Feature Extraction Based on Eight-Star Model

The image of the human object is a two-dimensional array. The extracted features of the human object are transformed into a one-dimensional feature vector, which greatly reduces the feature vector’s dimension and computational complexity. This article focuses on walking, running, and jumping in three postures for feature extraction and recognition, through the analysis of three kinds of attitude, and the legs move in arm swinging amplitude, arising in the process of adopting the star model to the human body modeling, namely through the calculation of center of mass of human movement target image points and local contour pole and its mutual relations, in order to obtain the attitude. In this paper, an eight-star model is utilized to extract a set of feature vectors containing seventeen feature values to characterize the human pose in each frame [16]. Compared to the six-star model, it can describe the human body’s contours, including the head and limbs, more precisely, and it is not redundant.Eight-star model is an improvement of the traditional star model, which is a feature description model formed according to the eight local contour poles and centroid points of the human object. The human pose model is established as shown in the following equation:(6)Fet=D,A,e,where D is the distance between the eight local contour poles of the eight-star model and the centroid of the human target. A is the minimum Angle formed by the straight line and horizontal line from the eight local contour poles to the centroid point of the eight-star model; e is the eccentricity of the human target. The specific feature extraction process is as follows:Calculate the centroid of human objectxc,yc:(7)xc=1Nb∑i=1Nbxi,yc=1Nb∑i=1Nbyi,where Nb is the number of target pixels.The target is divided into four parts by horizontal lines and vertical lines passing through the center of mass. Based on the overall contour of the target, eight local contour pole coordinatesxi,yi of the rightmost, uppermost, leftmost, and bottommost contour in each part are extracted, respectively.The Euclidean distances between the eight local contour polesxi,yi and the centroid point xc,yc are calculated, which are, respectively, expressed as di, and the calculation formula is shown in the following equation:(8)di=xi−xc2+yi−yc2,i=1,2,⋯,8.Then,D=d1,d2,⋯,d8. The extraction process is shown in Figure 2.Figure 2
Feature extraction process of pole distance of local contour of human eight-star model.After feature extraction and fusion, the feature vector of the human pose in the video at framet is shown in the following equation:(9)Fett=d1,d2,⋯,d8,a8,e.The feature vector is a 17-dimensional feature vector constructed from the local contour pole distance feature, local contour pole Angle feature, and eccentricity feature of the human eight-star model, wheredi is the pole distance feature of the eight-star local contour, ai is the Angle feature of the human local contour pole, and e is the eccentricity of the human body. The computation complexity and memory of the gesture description operator are relatively low, so the system can achieve a good real-time recognition effect.In practice, the video image will contain some noise and affect the effect of detection and recognition. Therefore, the primary task of video image processing is image preprocessing, and its effect directly affects the subsequent detection and recognition.


## 2.4. Sports Color Image Grayscale

The color value of the true color image pixel is composed of three components: red R, green G, and blue B. The value of each component has 256 levels of brightness, which is represented by numbers from 0 to 255. Therefore, there are 224 colors, which is relatively complicated to process directly. The gray image is a special color image, whose R, G, and B components are all equal, so there are only 28 colors left, which is convenient for image processing [17]. Therefore, in order to realize the real-time performance of the algorithm, the original image is converted into gray image before the video image processing, so as to reduce the computation amount of the algorithm. In general, the gray-level methods used include component method, maximum method, average method, and weighted average method.In the component method, one of the three components of RGB is used as the gray value of the pixel, as shown in the following equation:(10)I=RorI=GorI=B,where I represents the pixel brightness of gray image, and R, G, and B, respectively, represent the brightness values of RGB components. A gray image can be selected according to the application needs. Taking the 7th frame of LENA_WALK1 video, the 10th frame of PETS2006 video set, and the 21st frame of test video set in Weizmann database as examples, the grayscale is realized by component method.The above method is simple to implement but lacks consideration of the importance of different components. According to the importance of components and other needs of images, the weighted average of RGB three components can be taken as the gray value of pixels [18]. From the perspective of human physiology, the human eye is sensitive to green and insensitive to blue. Therefore, according to the visual model of human eyes, the weight of G component is set to be large, while the weight of B component is small, and its commonly used formula is shown in the following equation:(11)I=0.3R+0.5G+0.11B.The gray weights adopted by OpenCV open source computer vision library also follow the visual model of human eyes, and the weighted average formula is shown in the following equation:(12)I=0.212671R+0.715160G+0.72169B.From the perspective of theory and practice, the weighted average method is more reasonable and more in line with the needs of practical application, so the weighted average method is the most commonly used method of image graying. In this paper, the weighted average method is used to gray-scale the image [19], and the cvCvtColor function in the OpenCV computer vision open source library is used to complete the gray-scale of the image by setting its third parameter as CV_BGR2GRAY. Or using cvLoadImage function, by setting its second parameter to 0 can realize the color image gray.


## 3. An Accurate Estimation Algorithm for Sports Posture

### 3.1. Median Filtering Method for Posture Image Processing

In this method, the arithmetic mean or median of the pixels of the same coordinate in the continuous image frame is taken as the pixel value of the point in the background model, so as to realize the background modeling. The background modeling process is shown in (1314), respectively.(13)Btx,y=1k∑i=1kIt=ix,y.(14)Btx,y=medIi(x,yi=1,2,⋯,k).Among them,med means taking the median operation, Itx,y means the video frame at time t, and Btx,y means the updated background model at time t. However, it can be seen from formula (14) that the memory requirement increases because k frames of images need to be cached for updating the background model. Therefore, in order to reduce memory and computation, the background model parameters are usually updated according to the following equation:(15)Bt+1x,y=λItx,y+1−λBtx,y,where λ is called the update rate, which controls the update rate of the background model. Through experimental verification, in order to reduce the influence of foreground change on the model, the value of λ is usually small, generally set at 0.05.


### 3.2. Multifeature Fusion of Human Motion Posture

Multifeature fusion refers to the comprehensive analysis and processing of multisource feature vector information extracted from different feature extraction methods through feature fusion algorithm to form a new fusion feature vector group. In general, it can be divided into two categories: direct feature combination and feature selection combination. Direct feature combination method is to directly synthesize new feature vector groups by all feature vectors according to some simple rules, for example, serial and parallel fusion method. The method of feature selection and combination is to put all the original feature vectors together, and according to some selection rules, select and retain some feature vectors from the original feature vectors as a new combination of feature vectors.Serial feature fusion method: the feature vectors in the sample space are directly merged into new feature vectors in turn. Assume thatA and B are standardized features in sample space Ω, any sample ξ∈Ω, feature vector is expressed as α∈Ω, and β∈Ω fused feature is expressed as(16)γ=αβ.Parallel feature fusion method [20]: two groups of feature vectors in the sample space are merged into a new feature vector by a complex vector. Assuming that A and B are standardized features in sample space Ω, and any sample ξ∈Ω, feature vectors are expressed as α∈Ω and β∈Ω, and then, the fused features are expressed as(17)γ=α+iβ.In (17), i is an imaginary unit. If the dimensions of features α and β are not equal, feature fusion will be carried out after zero-complementing for features of low dimension.The serial fusion method is the simplest and most effective method for fusing features. The dimension of the fused feature is increased, all feature information is retained, and there are numerous redundant feature vectors. The parallel fusion method guarantees the invariance of feature dimension, but the use of complex space vectors increases the operation’s complexity, which has an effect on the recognition’s real-time performance. In order to improve accuracy and real-time performance, the feature selection method generates new features that can ensure the effectiveness within and between feature classes without generating an excessive amount of redundancy. Feature level fusion method is roughly divided into five categories: probability statistical method (Bayesian estimation, Kalman filter, etc.), logic reasoning method (D-S evidence theory, fuzzy logic, etc.), neural network methods, fusion method based on feature extraction (PCA, LDA, etc.), and fusion method based on search (genetic algorithm, particle swarm optimization (pso), etc.).


### 3.3. Human Pose Estimation in Sports Based on Mixed Gaussian Distribution Model

The method of mixed Gaussian background modeling is to update the parameters dynamically, so that the background model established by it has better adaptability. Its basic principle is to assume that every pixel in the image is independent and can be mixed by K Gaussian distributions, whereK is generally 3∼5.At timet, the pixel value at position i,j is xt, and the probability of each pixel value can be represented by the weighted sum of K Gaussian probability density functions, as shown in the following equation:(18)pxt=∑i=1Kωi,t×ηxt,μi,t,∑i,t,where ωi,t is the weighting coefficient of the ith Gaussian distribution at time t, the mean and covariance matrix of the ith Gaussian distribution at time t are μi,t, ∑i,t, and ∑i=1Kωi,t=1, and the probability density function of its Gaussian distribution is ηxt,μi,t,∑i,t. The specific calculation is shown in the following equation:(19)ηxt,μi,t,∑i,t=12πn/2∑i,t1/2e−1/2xt−μi,tT∑t−1xt−μi,t,when a new video frame is input, according to (19), the pixel value xt of each newly input image is successively compared with the K models already established at the current pixel point of this position. Where D is the confidence parameter, usually set at 2.5.(20)xt−μi,t−1≤D×σi,t−1.If the new input pixel valuext satisfies (20), then xt is judged as the background point after successful matching; otherwise, it is the foreground point. After the new video frame input pixel matching is completed, the model parameters of this point are updated and Mk,t=1 is taken, and the parameter update method is shown in equations (18) to (20), where α is the learning rate is a fixed value, and the value is 0<α<1, and ρ can be updated through the following equation:(21)ωk,t=1−αωk,t−1+α×Mk,t,(22)μt=1−ρ×μt−1+ρ×xt,σt2=1−ρ×σt−12+ρ×xt−μtTxt−μt,(23)ρ=α×ηxt|μk,σk.If no matching is achieved,Mk,t=0 is taken, which is only newer (23), and the rest is unchanged. If none of the models is matched successfully, then the model with the lowest weight in the mixed Gaussian model is replaced, and then, the new parameter mean is xt, the standard deviation is the initial maximum value σ0, and the weight is the small value ω0.In order to obtain the model with good correlation between background and Gaussian model, that is, the model with heavy weight and small standard deviation, all Gaussian models are arranged from large to small according toω/α2 value to obtain the sports pose estimation function B, which is expressed as(24)B=argminb∑k=1bωk>T.In the formula, the parameterT represents the proportion of background. If xt in the input video frame successfully matches any of the best description model, xt is judged to be the background point. Otherwise, it is judged as the foreground point. At this time, the obtained sports pose estimation function value is the result of sports pose estimation, so as to achieve accurate sports pose estimation.


## 3.1. Median Filtering Method for Posture Image Processing

In this method, the arithmetic mean or median of the pixels of the same coordinate in the continuous image frame is taken as the pixel value of the point in the background model, so as to realize the background modeling. The background modeling process is shown in (1314), respectively.(13)Btx,y=1k∑i=1kIt=ix,y.(14)Btx,y=medIi(x,yi=1,2,⋯,k).Among them,med means taking the median operation, Itx,y means the video frame at time t, and Btx,y means the updated background model at time t. However, it can be seen from formula (14) that the memory requirement increases because k frames of images need to be cached for updating the background model. Therefore, in order to reduce memory and computation, the background model parameters are usually updated according to the following equation:(15)Bt+1x,y=λItx,y+1−λBtx,y,where λ is called the update rate, which controls the update rate of the background model. Through experimental verification, in order to reduce the influence of foreground change on the model, the value of λ is usually small, generally set at 0.05.


## 3.2. Multifeature Fusion of Human Motion Posture

Multifeature fusion refers to the comprehensive analysis and processing of multisource feature vector information extracted from different feature extraction methods through feature fusion algorithm to form a new fusion feature vector group. In general, it can be divided into two categories: direct feature combination and feature selection combination. Direct feature combination method is to directly synthesize new feature vector groups by all feature vectors according to some simple rules, for example, serial and parallel fusion method. The method of feature selection and combination is to put all the original feature vectors together, and according to some selection rules, select and retain some feature vectors from the original feature vectors as a new combination of feature vectors.Serial feature fusion method: the feature vectors in the sample space are directly merged into new feature vectors in turn. Assume thatA and B are standardized features in sample space Ω, any sample ξ∈Ω, feature vector is expressed as α∈Ω, and β∈Ω fused feature is expressed as(16)γ=αβ.Parallel feature fusion method [20]: two groups of feature vectors in the sample space are merged into a new feature vector by a complex vector. Assuming that A and B are standardized features in sample space Ω, and any sample ξ∈Ω, feature vectors are expressed as α∈Ω and β∈Ω, and then, the fused features are expressed as(17)γ=α+iβ.In (17), i is an imaginary unit. If the dimensions of features α and β are not equal, feature fusion will be carried out after zero-complementing for features of low dimension.The serial fusion method is the simplest and most effective method for fusing features. The dimension of the fused feature is increased, all feature information is retained, and there are numerous redundant feature vectors. The parallel fusion method guarantees the invariance of feature dimension, but the use of complex space vectors increases the operation’s complexity, which has an effect on the recognition’s real-time performance. In order to improve accuracy and real-time performance, the feature selection method generates new features that can ensure the effectiveness within and between feature classes without generating an excessive amount of redundancy. Feature level fusion method is roughly divided into five categories: probability statistical method (Bayesian estimation, Kalman filter, etc.), logic reasoning method (D-S evidence theory, fuzzy logic, etc.), neural network methods, fusion method based on feature extraction (PCA, LDA, etc.), and fusion method based on search (genetic algorithm, particle swarm optimization (pso), etc.).


## 3.3. Human Pose Estimation in Sports Based on Mixed Gaussian Distribution Model

The method of mixed Gaussian background modeling is to update the parameters dynamically, so that the background model established by it has better adaptability. Its basic principle is to assume that every pixel in the image is independent and can be mixed by K Gaussian distributions, whereK is generally 3∼5.At timet, the pixel value at position i,j is xt, and the probability of each pixel value can be represented by the weighted sum of K Gaussian probability density functions, as shown in the following equation:(18)pxt=∑i=1Kωi,t×ηxt,μi,t,∑i,t,where ωi,t is the weighting coefficient of the ith Gaussian distribution at time t, the mean and covariance matrix of the ith Gaussian distribution at time t are μi,t, ∑i,t, and ∑i=1Kωi,t=1, and the probability density function of its Gaussian distribution is ηxt,μi,t,∑i,t. The specific calculation is shown in the following equation:(19)ηxt,μi,t,∑i,t=12πn/2∑i,t1/2e−1/2xt−μi,tT∑t−1xt−μi,t,when a new video frame is input, according to (19), the pixel value xt of each newly input image is successively compared with the K models already established at the current pixel point of this position. Where D is the confidence parameter, usually set at 2.5.(20)xt−μi,t−1≤D×σi,t−1.If the new input pixel valuext satisfies (20), then xt is judged as the background point after successful matching; otherwise, it is the foreground point. After the new video frame input pixel matching is completed, the model parameters of this point are updated and Mk,t=1 is taken, and the parameter update method is shown in equations (18) to (20), where α is the learning rate is a fixed value, and the value is 0<α<1, and ρ can be updated through the following equation:(21)ωk,t=1−αωk,t−1+α×Mk,t,(22)μt=1−ρ×μt−1+ρ×xt,σt2=1−ρ×σt−12+ρ×xt−μtTxt−μt,(23)ρ=α×ηxt|μk,σk.If no matching is achieved,Mk,t=0 is taken, which is only newer (23), and the rest is unchanged. If none of the models is matched successfully, then the model with the lowest weight in the mixed Gaussian model is replaced, and then, the new parameter mean is xt, the standard deviation is the initial maximum value σ0, and the weight is the small value ω0.In order to obtain the model with good correlation between background and Gaussian model, that is, the model with heavy weight and small standard deviation, all Gaussian models are arranged from large to small according toω/α2 value to obtain the sports pose estimation function B, which is expressed as(24)B=argminb∑k=1bωk>T.In the formula, the parameterT represents the proportion of background. If xt in the input video frame successfully matches any of the best description model, xt is judged to be the background point. Otherwise, it is judged as the foreground point. At this time, the obtained sports pose estimation function value is the result of sports pose estimation, so as to achieve accurate sports pose estimation.


## 4. Experiment

### 4.1. Experimental Design

The experiment was conducted on two benchmark datasets, JPL First Person Interaction Dataset (tJPL) and DogCentric Activity Dataset (Dogcentric), containing 84 and 209 videos, respectively. The former is the human-computer interaction data set, and the latter is the animal-human interaction data set from the first perspective. The video in both data sets has a resolution of 320 240 and a frame rate of 30 FPS.In this paper, multiple evaluation indicators are used to evaluate the model, including sample accuracy, sample level category accuracy, sample level average category accuracy, and video accuracy. Sample accuracy, also known as overall accuracy, refers to the proportion of correctly classified samples in the total number of samples and is the simplest and intuitive evaluation index in classification problems. Specifically, assuming thatN sample segments are obtained from the video according to the sampling strategy, and the number of segments correctly classified by the model is Ncor, and then, the formula for calculating the sample accuracy is as follows:(25)Acc=NcorN.However, when the proportion of categories is seriously unbalanced, it is not reliable to rely only on the sample accuracy as the evaluation standard. Suppose an extreme case: when the proportion of negative samples is 99% of the total samples, the classifier can also obtain a high accuracy by predicting all samples as negative samples. When the data are skewed, the categories with a large proportion become the main factors affecting the assessment criteria. Therefore, we also use the sample level category accuracy and sample level average category accuracy. Among them, the sample-level category accuracy is to group all sample segments according to categories and obtain the probability of correct prediction in each category. The average category accuracy at sample level is averaged according to the weighted summation of each category accuracy. Specifically, assuming that the number of classes isnum_classes (num_classes>0), each group is grouped according to the category, and the number of samples of num_classesnum_classes>0 in each group is Nncor0≤Nncor≤Nn, and then,(26)classAcc=NncorNn,MclassAcc=∑n=0num_class−1Nncor/Nnnum-classes.In order to obtain the video level prediction from the prediction metrics of the sample segments, the method described in Section3.2 is used for voting fusion, that is, the fusion of the maximum number of category predictions, the maximum category scores, and the maximum weighted scores, and the estimation accuracy of sports posture is calculated according to the fusion results.


### 4.2. Experimental Result

#### 4.2.1. Accurate Estimation of Athletic Posture Takes Time

In order to verify the precision estimation efficiency of the proposed method, the method in literature [6], the method in literature [7], and the method in this paper are used to estimate the precise estimation time of sports posture, and the results are shown in Figure 3.Figure 3
Time consumption for accurate estimation of sports posture.Analyzing Figure3 reveals that when the sample data amount of sports posture is 100 GB, the estimation time of the method in reference [6] is 43.5 s, the estimation time of the method in reference [7] is 35.2 s, and the estimation time of this method is 1.5 s. When the sample size of sports posture data is 300 GB, the estimation time of the method in reference [6] is 42.1 s, reference [7] method is 31.6 s, and the proposed method takes 2.3 s. When the sample data amount of sports posture is 900 GB, the accurate estimation time of the method in reference [6] is 46.5 s, the method in reference [7] is 58.3 s, and this method’s accurate estimation time is 5.3 s. The posture accuracy of the proposed method is significantly higher than that of other methods, indicating that it can effectively improve the posture accuracy of sports.


#### 4.2.2. Sample Accuracy of Posture Estimation in Sports

In order to verify the sample validity of the proposed motion pose estimation method, literature [6], literature [7] and the methods proposed in the paper are used for comparison, and the results are shown in Figure 4.Figure 4
Estimation accuracy of sports posture samples.According to Figure4, when the amount of sports pose sample data is 200 GB, the estimation accuracy of the method in reference [6] is 72.5 percent, the method in reference [7] is 48.6 percent, and the proposed method is 95 percent. When the amount of data for the sports pose sample is 600 GB, the estimation accuracy of the method in reference [6] is 69.9%, and that of the method in reference [7] is 71.5%. The proposed method has a 99.9% estimation precision. The estimation accuracy of the sports pose sample is 69.2% for the method in reference [6] and 89.6% for the method in reference [7] when the sample data size is 900 GB. The estimation precision of the proposed method is one hundred percent. This indicates that the proposed method can effectively improve the estimation accuracy of sports pose samples.


## 4.1. Experimental Design

The experiment was conducted on two benchmark datasets, JPL First Person Interaction Dataset (tJPL) and DogCentric Activity Dataset (Dogcentric), containing 84 and 209 videos, respectively. The former is the human-computer interaction data set, and the latter is the animal-human interaction data set from the first perspective. The video in both data sets has a resolution of 320 240 and a frame rate of 30 FPS.In this paper, multiple evaluation indicators are used to evaluate the model, including sample accuracy, sample level category accuracy, sample level average category accuracy, and video accuracy. Sample accuracy, also known as overall accuracy, refers to the proportion of correctly classified samples in the total number of samples and is the simplest and intuitive evaluation index in classification problems. Specifically, assuming thatN sample segments are obtained from the video according to the sampling strategy, and the number of segments correctly classified by the model is Ncor, and then, the formula for calculating the sample accuracy is as follows:(25)Acc=NcorN.However, when the proportion of categories is seriously unbalanced, it is not reliable to rely only on the sample accuracy as the evaluation standard. Suppose an extreme case: when the proportion of negative samples is 99% of the total samples, the classifier can also obtain a high accuracy by predicting all samples as negative samples. When the data are skewed, the categories with a large proportion become the main factors affecting the assessment criteria. Therefore, we also use the sample level category accuracy and sample level average category accuracy. Among them, the sample-level category accuracy is to group all sample segments according to categories and obtain the probability of correct prediction in each category. The average category accuracy at sample level is averaged according to the weighted summation of each category accuracy. Specifically, assuming that the number of classes isnum_classes (num_classes>0), each group is grouped according to the category, and the number of samples of num_classesnum_classes>0 in each group is Nncor0≤Nncor≤Nn, and then,(26)classAcc=NncorNn,MclassAcc=∑n=0num_class−1Nncor/Nnnum-classes.In order to obtain the video level prediction from the prediction metrics of the sample segments, the method described in Section3.2 is used for voting fusion, that is, the fusion of the maximum number of category predictions, the maximum category scores, and the maximum weighted scores, and the estimation accuracy of sports posture is calculated according to the fusion results.


## 4.2. Experimental Result

### 4.2.1. Accurate Estimation of Athletic Posture Takes Time

In order to verify the precision estimation efficiency of the proposed method, the method in literature [6], the method in literature [7], and the method in this paper are used to estimate the precise estimation time of sports posture, and the results are shown in Figure 3.Figure 3
Time consumption for accurate estimation of sports posture.Analyzing Figure3 reveals that when the sample data amount of sports posture is 100 GB, the estimation time of the method in reference [6] is 43.5 s, the estimation time of the method in reference [7] is 35.2 s, and the estimation time of this method is 1.5 s. When the sample size of sports posture data is 300 GB, the estimation time of the method in reference [6] is 42.1 s, reference [7] method is 31.6 s, and the proposed method takes 2.3 s. When the sample data amount of sports posture is 900 GB, the accurate estimation time of the method in reference [6] is 46.5 s, the method in reference [7] is 58.3 s, and this method’s accurate estimation time is 5.3 s. The posture accuracy of the proposed method is significantly higher than that of other methods, indicating that it can effectively improve the posture accuracy of sports.


### 4.2.2. Sample Accuracy of Posture Estimation in Sports

In order to verify the sample validity of the proposed motion pose estimation method, literature [6], literature [7] and the methods proposed in the paper are used for comparison, and the results are shown in Figure 4.Figure 4
Estimation accuracy of sports posture samples.According to Figure4, when the amount of sports pose sample data is 200 GB, the estimation accuracy of the method in reference [6] is 72.5 percent, the method in reference [7] is 48.6 percent, and the proposed method is 95 percent. When the amount of data for the sports pose sample is 600 GB, the estimation accuracy of the method in reference [6] is 69.9%, and that of the method in reference [7] is 71.5%. The proposed method has a 99.9% estimation precision. The estimation accuracy of the sports pose sample is 69.2% for the method in reference [6] and 89.6% for the method in reference [7] when the sample data size is 900 GB. The estimation precision of the proposed method is one hundred percent. This indicates that the proposed method can effectively improve the estimation accuracy of sports pose samples.


## 4.2.1. Accurate Estimation of Athletic Posture Takes Time

In order to verify the precision estimation efficiency of the proposed method, the method in literature [6], the method in literature [7], and the method in this paper are used to estimate the precise estimation time of sports posture, and the results are shown in Figure 3.Figure 3
Time consumption for accurate estimation of sports posture.Analyzing Figure3 reveals that when the sample data amount of sports posture is 100 GB, the estimation time of the method in reference [6] is 43.5 s, the estimation time of the method in reference [7] is 35.2 s, and the estimation time of this method is 1.5 s. When the sample size of sports posture data is 300 GB, the estimation time of the method in reference [6] is 42.1 s, reference [7] method is 31.6 s, and the proposed method takes 2.3 s. When the sample data amount of sports posture is 900 GB, the accurate estimation time of the method in reference [6] is 46.5 s, the method in reference [7] is 58.3 s, and this method’s accurate estimation time is 5.3 s. The posture accuracy of the proposed method is significantly higher than that of other methods, indicating that it can effectively improve the posture accuracy of sports.


## 4.2.2. Sample Accuracy of Posture Estimation in Sports

In order to verify the sample validity of the proposed motion pose estimation method, literature [6], literature [7] and the methods proposed in the paper are used for comparison, and the results are shown in Figure 4.Figure 4
Estimation accuracy of sports posture samples.According to Figure4, when the amount of sports pose sample data is 200 GB, the estimation accuracy of the method in reference [6] is 72.5 percent, the method in reference [7] is 48.6 percent, and the proposed method is 95 percent. When the amount of data for the sports pose sample is 600 GB, the estimation accuracy of the method in reference [6] is 69.9%, and that of the method in reference [7] is 71.5%. The proposed method has a 99.9% estimation precision. The estimation accuracy of the sports pose sample is 69.2% for the method in reference [6] and 89.6% for the method in reference [7] when the sample data size is 900 GB. The estimation precision of the proposed method is one hundred percent. This indicates that the proposed method can effectively improve the estimation accuracy of sports pose samples.


## 5. Conclusion

In this paper, a human rod model is constructed using a precise estimation algorithm for sports posture based on the fusion of multiple features. Based on the Kalman filter method to estimate the motion posture, determine the degree of freedom and range of motion of the main body joints, use the eight-star model to extract the motion posture features, use the multifeature fusion method to extract the multisource vector information features, and use the mixed Gaussian distribution model to realize the motion body posture. Experimental results indicate that when the amount of sports pose sample data is 900 GB, the accurate estimation time of the proposed method is 5.3 s, and its accuracy is 100 percent. Improve the estimation accuracy of samples of sports posture.In the future, we will introduce a multihead attention mechanism and more advanced 6D recognition technology to further improve the accuracy of motion gesture recognition.



---

*Source: 1224939-2022-10-11.xml*