Fundamental Understanding of Tea Growth and Modeling of Precise Tea Shoot Picking Based on 3-D Coordinate Instrument

Tea is a popular beverage worldwide and also has great medical value. A fundamental understanding of tea shoot growth and a precision picking model should be established to realize mechanized picking of tea shoots with a small product loss. Accordingly, the terminal bud length (Lbud), tea stem length (Lstem), terminal bud angle (αbud), tea stem angle (αstem), and growth time (t) were considered as the key growth parameters; the sum of the vertical lengths of the terminal bud and stem (ξ), the picking radius (r), and the vertical length of the stem (Zstem) were considered as the picking indexes of the tea shoots. The variations in growth parameters with time were investigated using a 3-D coordinate instrument, and the relationships between the growth parameters and the picking indexes were established using an artificial neural network (ANN). The results indicated that the tea growth cycles for periods P1, P2, P3, P4, P5, and P6 were 14, 7, 6, 4, 4, and 6 d, respectively. A growth cycle diagram of the tea growth was established. Moreover, a 5-2-12-3 ANN model was developed. The best prediction of ξ, r, and Zstem was found with 16 training epochs. The MSE value was 0.0923 × 10−4, and the R values for the training, test, and validation data were 0.99976, 0.99871, and 0.99857, respectively, indicating that the established ANN model demonstrates excellent performance in predicting the picking indexes of tea shoots.


Introduction
Tea (Camellia sinensis O. Kuntz) belongs to the Camellia genus in the Dicotyledon class. It originated in southwest China and has become a popular beverage worldwide [1]. Tea can also be used in medicine, food, and industry. Tea polyphenols have antioxidant and antimutagenic effects and can effectively prevent cancer [2,3], resulting in their widespread use in medicine. Moreover, tea catechins can stimulate liver lipid catabolism and have a positive effect on diet-induced obesity (Murase et al., 2002); tea catechin products are popular in the food industry. In recent years, with an increase in the demand for tea, tea planting has also expanded. The quality of tea picking has become a key factor affecting its economic value [4].
Tea with one shoot with one leaf is often regarded as first-class tea; tea with one shoot with two leaves is regarded as second-class tea. The third and fourth classes are shown in Figure 1. Japan was the first country to conduct research on mechanized tea picking; tea picking and tea tree pruning using large scissors had been widely used [5]. Although mechanized picking is suitable for large-scale tea picking, the method is relatively rough and does not locate the tea picking area; a considerable number of leaves are broken, In recent decades, extensive research on the characteristics of tea shoots has been conducted to obtain a corresponding growth model. Takashi (1962) conducted morphological studies on the development of tea shoots as early as 1960 [6]; he found that there were many lateral buds on a long branch, and most young leaves were in the upper part of the branch. Sanui (1965) recorded lateral shoot differentiation data and the development of tea shoots in summer and further studied the morphology of tea shoot development [7]. To better harvest tea with scissors, Nakayama (1967) studied the vertical distribution of tea shoots and found that lateral buds had unique growth characteristics [8]; they grow at low positions and are inhibited by terminal buds. Based on this analysis, the morphological differences of tea shoots in different seasons are significant. Harada (1960) planted potted, young tea trees in greenhouses at temperatures of 5-35 °C to determine the factors influencing tea shoot development and further studied the development morphology and growth model of the tea shoots [9]. The results showed that the tea shoots germinated at 12.5 °C, and the growth rate reached the maximum value at 27.5 °C. The temperature had a smaller influence on the shoot length between 27.5 and 35 °C [10]. However, many problems remain if the temperature is the only variable used for tea shoot prediction.
Many researchers have extensively considered other variables affecting tea shoot prediction accuracy. Koichi (1988) and Takezawa (1994) improved the prediction accuracy of tea shoot growth at different temperatures using water as a variable [11,12]. Although there are many factors affecting the development of tea shoots, and the same batch of tea shoots differs in morphology, growth is generally similar in the same environment, which can provide a reference for follow-up research. To eliminate the influence of temperature and water on tea shoot development, Takayuki (1999) and  introduced growth time as a variable to investigate the variation in tea shoot growth over time in a natural environment [13,14]. Jayasinghe established models for predicting the leaf area and fresh weight of tea shoots at different development stages and locations. The results demonstrated that the models were in good agreement with the measured values at the validation stage (R 2 > 0.92 for the leaf area model and R 2 > 0.98 for the fresh weight model), In recent decades, extensive research on the characteristics of tea shoots has been conducted to obtain a corresponding growth model. Takashi (1962) conducted morphological studies on the development of tea shoots as early as 1960 [6]; he found that there were many lateral buds on a long branch, and most young leaves were in the upper part of the branch. Sanui (1965) recorded lateral shoot differentiation data and the development of tea shoots in summer and further studied the morphology of tea shoot development [7]. To better harvest tea with scissors, Nakayama (1967) studied the vertical distribution of tea shoots and found that lateral buds had unique growth characteristics [8]; they grow at low positions and are inhibited by terminal buds. Based on this analysis, the morphological differences of tea shoots in different seasons are significant. Harada (1960) planted potted, young tea trees in greenhouses at temperatures of 5-35 • C to determine the factors influencing tea shoot development and further studied the development morphology and growth model of the tea shoots [9]. The results showed that the tea shoots germinated at 12.5 • C, and the growth rate reached the maximum value at 27.5 • C. The temperature had a smaller influence on the shoot length between 27.5 and 35 • C [10]. However, many problems remain if the temperature is the only variable used for tea shoot prediction.
Many researchers have extensively considered other variables affecting tea shoot prediction accuracy. Koichi (1988) and Takezawa (1994) improved the prediction accuracy of tea shoot growth at different temperatures using water as a variable [11,12]. Although there are many factors affecting the development of tea shoots, and the same batch of tea shoots differs in morphology, growth is generally similar in the same environment, which can provide a reference for follow-up research. To eliminate the influence of temperature and water on tea shoot development, Takayuki (1999) and Jayasinghe (2018) introduced growth time as a variable to investigate the variation in tea shoot growth over time in a natural environment [13,14]. Jayasinghe established models for predicting the leaf area and fresh weight of tea shoots at different development stages and locations. The results demonstrated that the models were in good agreement with the measured values at the validation stage (R 2 > 0.92 for the leaf area model and R 2 > 0.98 for the fresh weight model), indicating that predicting the growth of tea shoots in the same environment considering growth time was relatively accurate. Similar findings have been reported for other plant growth, including rice yield [15], maize and soybean yield [16], and crop harvest date [17]. These studies have produced a relatively accurate tea shoot growth model. With technological innovation, an increasing number of researchers have developed more accurate plant growth models using advanced equipment. Bernotas et al. [18] studied the growth structure of Arabidopsis thaliana under different conditions using a low-cost portable three-dimensional (3-D) plant phenotypic platform (PS-Plant) based on photometric stereo imaging technology. The results demonstrated that PS-Plant is a powerful phenotypic tool that can accurately predict the spatial structure of plants during different growth periods. Bengocheaguevara et al. [19] developed a crop detection system using RGB-D sensors that automatically created accurate 3-D models (color point clouds) of woody crops. Golbach et al. [20] constructed a 3-D imaging space using ten cameras to produce 3-D imaging of a single tomato plant; the leaf length, width, area, and stem length were also obtained. Li et al. [21] introduced time analysis into inexpensive portable 3-D acquisition equipment to analyze plant growth from 4-D point cloud data. The research included detection and quantification of specific growth events (such as sprouting and branching) and further tracking their evolution over time. In another study, a laser scanner system was used to obtain the 3-D structural characteristics of vines [22]. These modeling methods must scan the entire plant, producing a large amount of original data. Partial data must be extracted from the original data when specific points are required for analysis, which increases the data processing difficulty. Whitfield et al. [23] quantitatively described the architectural characteristics (length, width, and angle) of individual tobacco plants using a spatial coordinate instrument; accurate leaf orientation distributions were ascertained by combining an ellipsoid model and the Fisher distribution. Sinoquet et al. [24] used ultrasonic and electromagnetic devices to conduct 3-D digitization of the canopy; architectural characteristics (length, width, and angle) of the canopy were obtained at each leaf position. Although these methods require point-by-point measurement, field operation is simple. In addition, the 3-D space coordinates of the target point can be obtained directly without damage. A method combining 3-D space coordinate technology with growth time analysis may be an effective means of predicting the growth of tea shoots in different growth periods.
An artificial neural network (ANN) is an anthropomorphic training method that is simple to use, especially in processing a large amount of data. Wu et al. [25] used a probabilistic neural network (PNN) combined with image and data processing technology to automatically identify leaves for plant classification. Twelve leaf features were extracted and orthogonalized into five principal variables used as the input neurons of the PNN. The PNN was trained with 1800 leaves to classify 32 types of plants; the results demonstrated a classification accuracy greater than 90%. Considering the wheat yield, associated plantation area, rainfall, and temperature as the input neurons, Guo et al. [26] established an ANN model for predicting wheat yield; the results demonstrated that the ANN model was highly reliable. Kucukonder et al. [27] aimed to develop the best leaf area estimation model for tomato plants grown in a plastic greenhouse. The ANN methodology and regression analysis were used to estimate the leaf area using the leaf width and leaf length as the key variables. The statistical indices of the ANN model (R 2 = 0.96, RMSE = 3.30, MAE = 1.94, and MAPE = 0.05) were lower than those of the regression model (R 2 = 0.92, RMSE = 4.71, MAE = 3.31, and MAPE = 0.08), indicating a better prediction performance than that of the regression model. Hu et al. [28] presented a photosynthetic rate model of a heuristic neural network for tomatoes based on a genetic algorithm. The test results showed that the training effects and accuracy of the genetic neural network prediction of the photosynthetic rate were good; the correlation coefficient between the predicted data and the measured data was 0.987, and the absolute error of the photosynthetic rate was less than ±0.5 µmol/(m 2 ·s). Green et al. [29] presented a spatial analysis neural network (SANN) algorithm combined with topographic attributes as explanatory variables for joint spatial interpolation and yield prediction. The results indicated the utility of SANN with topographic attributes containing implicit soil and water information for estimating the spatial patterns of dryland crop yield. From the analysis, the ANN methodology is an accurate and fast artificial intelligence algorithm for crop morphology and yield prediction and can be used for tea shoot growth prediction.
To realize precise picking mechanization for tea shoots, a fundamental understanding of tea shoot growth and a corresponding prediction model are necessary. This study focuses on the growth rule of tea shoots, applying a 3-D coordinate instrument to monitor key indices including terminal bud length (L bud ), stem length (L stem ), terminal bud angle (α bud ), and stem angle (α stem ) in different growth periods. The ANN methodology was applied to establish the growth rule prediction model considering L bud , L stem , α bud , α stem , and growth time (t) as the input neurons, and the sum of the vertical length of the shoot and stem (ξ), the picking radius (r), and the vertical length of the stem (Z stem ) as the output neurons. This study intends to establish a fundamental basis for guiding practical tea picking and the design of a precise picking machine.

Sample Preparation
The experiment was conducted on 102 Yinghong #9 potted tea trees planted in the tea garden of the North building of the Engineering College of South China Agricultural University (longitude 113.3456 • , latitude 23.1587 • ). The height of the tea trees ranged from 0.9 to 1 m; the five-year-old trees were planted in full light. To reduce interference caused by uncontrollable factors, the tea trees were uniformly managed in terms of watering and fertilization during the data acquisition period. Tea shoots with one shoot and three intact leaves were used in the experiment.

Testing Apparatus
According to Tang et al. [30], the quality of picking is an important factor affecting the value of the tea; the location of the picking point is the decisive factor in the picking quality [31]. Thus, this study used a 3-D coordinate instrument to measure the spatial coordinates of tea shoots. The device was purchased from Polhemus and consisted of a system electronic unit (SEU), a power supply, a sensor, a transmitter, and a computer, as shown in Figure 2. The testing accuracy within a circular area with a radius of 2 m was 0.1 mm. A handheld anemometer (Xima AS836) was used to measure the wind speed.
accurate and fast artificial intelligence algorithm for crop morphology and yield tion and can be used for tea shoot growth prediction.
To realize precise picking mechanization for tea shoots, a fundamental und ing of tea shoot growth and a corresponding prediction model are necessary. Th focuses on the growth rule of tea shoots, applying a 3-D coordinate instrument to key indices including terminal bud length (Lbud), stem length (Lstem), terminal bu (αbud), and stem angle (αstem) in different growth periods. The ANN methodology plied to establish the growth rule prediction model considering Lbud, Lstem, αbud, α growth time (t) as the input neurons, and the sum of the vertical length of the sh stem (ξ), the picking radius (r), and the vertical length of the stem (Zstem) as th neurons. This study intends to establish a fundamental basis for guiding practical ing and the design of a precise picking machine.

Sample Preparation
The experiment was conducted on 102 Yinghong #9 potted tea trees plante tea garden of the North building of the Engineering College of South China Agr University (longitude 113.3456°, latitude 23.1587°). The height of the tea trees rang 0.9 to 1 m; the five-year-old trees were planted in full light. To reduce interferenc by uncontrollable factors, the tea trees were uniformly managed in terms of wate fertilization during the data acquisition period. Tea shoots with one shoot and thr leaves were used in the experiment.

Testing Apparatus
According to Tang et al. [30], the quality of picking is an important factor the value of the tea; the location of the picking point is the decisive factor in the quality [31]. Thus, this study used a 3-D coordinate instrument to measure the sp ordinates of tea shoots. The device was purchased from Polhemus and consisted tem electronic unit (SEU), a power supply, a sensor, a transmitter, and a comp shown in Figure 2. The testing accuracy within a circular area with a radius of 2 m mm. A handheld anemometer (Xima AS836) was used to measure the wind spee

Data Acquisition
Data was collected from 1 March to 23 October 2020 in the open roof tea garden with sufficient light. A total of 46 sets of data were collected during the experiment, between 7:00 and 11:00 a.m. Each set of data was measured three times, and the average value was used for calculation. To ensure the validity of the data and reduce the disturbance caused by wind, the data were recorded when the wind speed was at level 0 (0.0-0.2 m/s).

Biometric Analysis
The biological characteristics of the abstract structure of Yinghong #9 tea shoots are shown in Figure 3. The physical characteristics of the tea shoots are unusual; a reasonable statistical method should be used. The main indexes for tea picking are L bud , L stem , α bud , and α stem , which are important biological characteristics of tea shoots. These indexes are important in designing the shape and response time of the end-effector.

Data Acquisition
Data was collected from 1 March to 23 October 2020 in the open roof tea garden with sufficient light. A total of 46 sets of data were collected during the experiment, between 7:00 and 11:00 a.m. Each set of data was measured three times, and the average value was used for calculation. To ensure the validity of the data and reduce the disturbance caused by wind, the data were recorded when the wind speed was at level 0 (0.0-0.2 m/s).

Biometric Analysis
The biological characteristics of the abstract structure of Yinghong #9 tea shoots are shown in Figure 3. The physical characteristics of the tea shoots are unusual; a reasonable statistical method should be used. The main indexes for tea picking are Lbud, Lstem, αbud, and αstem, which are important biological characteristics of tea shoots. These indexes are important in designing the shape and response time of the end-effector. As shown in Figure 3, the variation trends of the length and inclination angle for the stem internode and terminal bud are different in different growth periods. As reported by [30], one shoot with one leaf is one type of tea grading. From the morphological description of tea shoots, the terminal bud of the previous growth period differentiates into the terminal bud and leaf of the next growth period after a period of incubation and varies with the seasons. Normally, without human intervention, the branches circulate until the tea shoots become folium dichotomous, which marks the end of the branching process [32]. According to the branching process, the tea picking period is between Figure 4a,c. To obtain sufficient yield in actual production, experienced workers often pick in the late stage of one shoot with three leaves (late period c).  As shown in Figure 3, the variation trends of the length and inclination angle for the stem internode and terminal bud are different in different growth periods. As reported by [30], one shoot with one leaf is one type of tea grading. From the morphological description of tea shoots, the terminal bud of the previous growth period differentiates into the terminal bud and leaf of the next growth period after a period of incubation and varies with the seasons. Normally, without human intervention, the branches circulate until the tea shoots become folium dichotomous, which marks the end of the branching process [32]. According to the branching process, the tea picking period is between Figure 4a,c. To obtain sufficient yield in actual production, experienced workers often pick in the late stage of one shoot with three leaves (late period c).

Data Acquisition
Data was collected from 1 March to 23 October 2020 in the open roof tea garden with sufficient light. A total of 46 sets of data were collected during the experiment, between 7:00 and 11:00 a.m. Each set of data was measured three times, and the average value was used for calculation. To ensure the validity of the data and reduce the disturbance caused by wind, the data were recorded when the wind speed was at level 0 (0.0-0.2 m/s).

Biometric Analysis
The biological characteristics of the abstract structure of Yinghong #9 tea shoots are shown in Figure 3. The physical characteristics of the tea shoots are unusual; a reasonable statistical method should be used. The main indexes for tea picking are Lbud, Lstem, αbud, and αstem, which are important biological characteristics of tea shoots. These indexes are important in designing the shape and response time of the end-effector. As shown in Figure 3, the variation trends of the length and inclination angle for the stem internode and terminal bud are different in different growth periods. As reported by [30], one shoot with one leaf is one type of tea grading. From the morphological description of tea shoots, the terminal bud of the previous growth period differentiates into the terminal bud and leaf of the next growth period after a period of incubation and varies with the seasons. Normally, without human intervention, the branches circulate until the tea shoots become folium dichotomous, which marks the end of the branching process [32]. According to the branching process, the tea picking period is between Figure 4a,c. To obtain sufficient yield in actual production, experienced workers often pick in the late stage of one shoot with three leaves (late period c).

Experimental Plan
The development of tea shoot branching is not uniform owing to different planting environments and pruning methods. Accordingly, in this study, light pruning was performed after each picking to reduce the differences between individual tea shoots. Light pruning can also stimulate the germination of tea shoots, relieve the suppression of the side shoots by the top shoots, keep the crown surface flat, and adjust the number and thickness of production branches. In this study, the Yinghong #9 growing data in spring, summer, and autumn were collected. Generally, spring, summer, and autumn teas can be picked twice, three times, and once, respectively, owing to the hot weather in Guangdong. According to the actual tea picking conditions in 2020, the coordinate data corresponding to the period of one shoot with three leaves (Figure 4c) 2 March-15 March (P 1 ), 2 May-9 May (P 2 ), 19 June-25 June (P 3 ), 15 July-19 July (P 4 ), 18 August-22 August (P 5 ), and 16 September-23 October (P 6 ) were collected. To obtain the growth characteristics of the entire process for autumn tea shoots (from dormant to the picking period of one shoot with three leaves), the data were recorded immediately after pruning. The detailed experimental plan is presented in Table 1.

Key Parameters for Mathematical Model of Precision Tea Picking
To establish a mathematical model of precision tea picking, the terminal buds of the tea shoots shown in Figure 5 were abstracted into spatial line AB; the spatial coordinates of points A and B were determined to be (X a , Y a , Z a ) and (X b , Y b , Z b ), as shown in Figure 5a. Point B was translated to axis Y, and the spatial coordinates of point B became (0, Y b , Z b ). The sets for different terminal buds are shown in Figure 5b. Point B was translated to the intersection of the X and Z axes, as shown in Figure 5c. Different terminal bud points (B i ) were translated following the same method; the sets of abstract lines are summarized in Figure 5d. From a spatial perspective, all of the terminal buds are clustered. Based on the analysis, a plane π was constructed to minimize the sum of the deviations from point set A to the π plane.

Multivariate Least Squares Linear Regression Analysis to Analyze Parameters
It is assumed that there are n points of ( ) )was determined to minimize the sum of the deviations from the point sets of Pj to π. Accordingly, the problem can be transformed into the following optimization problem: Similarly, a mathematical model for tea shoot picking can be obtained by analyzing and combining the BC segment of the stem, as shown in Figure 5e. The sum of the vertical lengths of terminal buds and stems is an important index for picking. The vertical length of the stem is regarded as a potential picking area and is directly related to whether the picked shoots have one shoot with one leaf. In theory, a longer Z stem indicates a greater probability of one shoot with one leaf being picked. Based on the statistics, it was found that the length and deflection angle of the terminal bud was always greater than those of the stems. Thus, the projection distance of the terminal bud was greater than that of the stem in the horizontal plane and is regarded as the picking radius. To protect the main body of the shoots during the picking process, r is also an important parameter in the subsequent design of the end-actuator. The key parameters in the precision tea picking model were determined to be ξ, Z stem , and r.

Multivariate Least Squares Linear Regression Analysis to Analyze Parameters
It is assumed that there are n points of P j (x j , y j , z j ), j = 1, 2, 3, . . . n (position of node B). A plane π(Z − ξ = Ax + By) was determined to minimize the sum of the deviations from the point sets of P j to π. Accordingly, the problem can be transformed into the following optimization problem: After finding the stagnation point, the system can be solved. Thus, Equation (10) is the regression plane, which can be solved by substituting the values of (x j , y j ) to obtain When the plane is parallel to the ground, the problem can be simplified as Node A produces the value closest to the observed value. Plane π may not be parallel to the ground, which may be a key issue in practical production. To simplify, plane π was assumed to be parallel to the ground, and the picking actuator was perpendicular to the ground. The algorithm represents the derivation process for the entire tea bud picking plane. Based on the algorithm, the picking length of tea shoots in each period can be calculated, and several important parameters can be obtained. The following equations can be obtained by combining the picking model shown in Figure 5.
α bud = arg cos Z bud L bud (21) α stem = arg cos Z stem L stem (22) where L bud , α bud , and Z bud are the actual length, inclination angle, and vertical length of the terminal bud, respectively; L stem , α stem , and Z stem are the actual length, inclination angle, and vertical length of the stem, respectively; r is the projection distance of the terminal bud on the horizontal plane. Using this method, we processed the original data in Table 1 to obtain the tea shoot parameters. Partial experimental data are presented in Table 2.

ANN Modeling for Precise Tea Picking Prediction
To investigate the relationships between the key growth parameters (L bud , L stem , α stem , α bud , t) of the tea, and the picking indexes (ξ, r, and Z stem ), a feed-forward backprop artificial neural network model was established; the key parameters of the model can be determined through the following steps.

Input and Output Layer Neurons
Based on the analysis, L bud , L stem , α bud , α stem , and t were considered as the input neurons of the ANN model; ξ, r, and Z stem were considered as the output neurons. Accordingly, there were four input and three output layers.

Hidden Layer Neurons
The number of hidden layers directly affects the training accuracy of the ANN model. According to Dai et al. [33], the number of hidden layers of an ANN model can be determined by the following equation: where m is the number of hidden layer neurons; a is the number of input layer neurons, which was five, as shown in Figure 6a; b is the number of output layer neurons, which was three; k is a constant between 1 and 10. Accordingly, m can be calculated to be in the range of 3-13.

Transfer and Train Function
The Tansig function provides a good tradeoff for neural networks, where speed is important, and the exact shape of the transfer function is not [34]. The Tansig function was used as the transfer function. The Levenberg-Marquardt algorithm was applied as the training function for the ANN model because it has a fast gradient descent and meets the error requirement with fewer epochs [35].

Physical Structure of Model
Based on the analysis, the established ANN and the corresponding training process are summarized in Figure 6b. esses 2021, 9, x FOR PEER REVIEW 10 of where m is the number of hidden layer neurons; a is the number of input layer neuro which was five, as shown in Figure 6a; b is the number of output layer neurons, whi was three; k is a constant between 1 and 10. Accordingly, m can be calculated to be in t range of 3-13.
(a) (b) Figure 6. Established ANN (a) and the corresponding training process (b).

Transfer and Train Function
The Tansig function provides a good tradeoff for neural networks, where speed important, and the exact shape of the transfer function is not [34]. The Tansig function w used as the transfer function. The Levenberg-Marquardt algorithm was applied as t training function for the ANN model because it has a fast gradient descent and meets t error requirement with fewer epochs [35].

Physical Structure of Model
Based on the analysis, the established ANN and the corresponding training proce are summarized in Figure 6b.

Key Training Parameters of the ANN Model
Scientifically, an ANN model should ascertain the number of neurons in the inp and outlet layers and the number of hidden layers but should also determine the traini parameters of the model, including the number of hidden neurons, momentum coef cients, and training epochs. According to Reed, R. D [36], for a simple dataset, one or tw hidden layers are usually enough. Accordingly, the present work considers the numb of hidden layers is in the range of 1-3. The momentum coefficient and training epochs a usually no more than 0.4 and 1500, respectively, as recommended by Chokphoemphu S. et al. [37] and Zeng, Z. et al. [38], accordingly, the ranges of momentum coefficient a training epochs are ascertained to be in the ranges of 0.1-0.4 and 300-1500, respective Based on the analysis in Sections 3.3.1 and 3.3.2, the configuration of the ANN model c be summarized, as shown in Table 3. Table 3. Configurations of ANN model.

Key Training Parameters of the ANN Model
Scientifically, an ANN model should ascertain the number of neurons in the input and outlet layers and the number of hidden layers but should also determine the training parameters of the model, including the number of hidden neurons, momentum coefficients, and training epochs. According to Reed, R. D [36], for a simple dataset, one or two hidden layers are usually enough. Accordingly, the present work considers the number of hidden layers is in the range of 1-3. The momentum coefficient and training epochs are usually no more than 0.4 and 1500, respectively, as recommended by Chokphoemphun, S. et al. [37] and Zeng, Z. et al. [38], accordingly, the ranges of momentum coefficient and training epochs are ascertained to be in the ranges of 0.1-0.4 and 300-1500, respectively. Based on the analysis in Section 3.3.1 and Section 3.3.2, the configuration of the ANN model can be summarized, as shown in Table 3. Table 3. Configurations of ANN model. In order to ascertain the architecture of the ANN model, the parameters in Table 3 are tuned by training the 134 sets of experimental data. The mean squared error (MSE) is adopted to evaluate the training results.

Statistical Analysis
The predictive performance of the established ANN model was estimated using statistical indexes including the coefficient of determination (R 2 ), the mean squared error (MSE), and the mean absolute error (MAE), which can be calculated using Equations (24)- (26), respectively [39].
where n is the number of data points; α i is the i-th experimental result; β i is the i-th result predicted by the ANN model.

Variation Rules for L bud , L stem , and ξ with Growth Time for Different Periods
To ascertain the variation rules of the picking indexes with growth time, the interval time between two adjacent terminal bud branches was considered as the growth cycle of the tea shoots. The variations in L bud , L stem , and ξ with t for different periods are shown in Figure 7; the growth cycles of the tea shoots are different for different periods. The growth cycles of P 1 , P 2 , P 3 , P 4 , P 5 , and P 6 were 13, 7, 6, 4, 4, and 6 d, respectively. Although the tea shoot growth cycle was significantly affected by seasonal differences, the change rate of L stem for all periods was essentially the same. The growth rates of L bud , L stem , and ξ for P 2 , P 3 , and P 6 were significantly reduced near the end of the growth cycle, as shown in Figure 7b,c,f. Liu et al. [35] reported that the terminal bud in the next growth cycle must be inoculated in advance for the nutrients in the terminal bud to be consumed, which may explain the growth cycle changes. It was found that L bud is larger than L stem in each period; the tea shoots sorted by length in descending order for different periods are represented by: P 1 , P 2 , P 3 , P 6 , P 4 , and P 5 .

Variation Rules for αbud and αstem with Growth Time for Different Periods
In Figure 5e, αbud and αstem are the key factors affecting r, and r is a key parameter Figure 7. Variation rules for L bud , L stem , and ξ with growth time for P 1 , P 2 , P 3 , P 4 , P 5 and P 6 (a-f).

Variation Rules for α bud and α stem with Growth Time for Different Periods
In Figure 5e, α bud and α stem are the key factors affecting r, and r is a key parameter for tea picking. Thus, the variation rules for α bud and α stem with respect to t were investigated. Figure 8 shows the variation of α bud and α stem with growth time. The terminal bud in the previous cycle grows into a new terminal bud and a stem in the next cycle; consequently, α bud is equal to α stem at the initial point of the growth cycle for all periods. It was observed that α stem decreased with an increase in growth time and tended to change gradually in the final stage of the growth cycle. However, α bud decreased rapidly during the initial stage of the growth cycle and further increased at a relatively rapid rate, perhaps because when the terminal bud reaches a critical point in the vertical state, it continues to incline in the same direction. The minimum values (all close to 0 • ) of α bud for P 1 , P 2 , P 3 , P 4 , P 5 , and P 6 were at 4, 2, 2, 1, 1, and 2 d of the growth cycle, respectively. In the final stage, α bud tended to change gradually, similar to the L bud trends, as shown in Figure 7. It is also observed in Figure 8 that α bud in the final stage is essentially equal to α bud at the initial point of the growth cycle for each period, possibly due to different biological characteristics.

Changes in Tea Growth Cycle with Time
To fully understand the growth rule for the Yinghong #9 tea shoot, the growth characteristics from the tea branch pruning period to the tea picking period were monitored; the growth cycle diagram is shown in Figure 9. It can be intuitively observed in Figure 9 that the transition period of the tea is nearly 20 d, nearly half of the entire process. This is because the transition period requires preparation of the nutrients for subsequent growth of the tea shoots. The next stage is the budding period, during which the bud gradually changes from dark brown to light green, and the bud length changes continuously. In the final stage of the budding period, the bud begins to branch and grow rapidly. Tea shoots then grow into the first stage of the germination period (a stage), where many top spring teas of greater value originate; the tea production yield is small, and the picking operation is difficult in this stage. Stages b and c are the stages of one shoot with two leaves and one shoot with three leaves, respectively. From a morphological perspective, the terminal buds in these two stages are very similar; the change trends of αbud and αstem in the c stage are essentially the same. In stage c, Lstem and Lbud increased with increasing growth time. In practical production, to increase production yield, workers often pick tea shoots in stage c. It is concluded from the figure that the c stage is 5-6 d. The appropriate picking time is near the end of stage c, which is approximately 1-2 d. Most of the shoots stopped branching when they had three leaves; the tea leaves became darker, and the stems gradually changed from green to brown, finally growing into folium dichotomous.
The growth of tea shoots is affected by different factors; the most influential is "top Variation rules for α bud and α stem with growth time for P 1 , P 2 , P 3 , P 4 , P 5 and P 6 (a-f).

Changes in Tea Growth Cycle with Time
To fully understand the growth rule for the Yinghong #9 tea shoot, the growth characteristics from the tea branch pruning period to the tea picking period were monitored; the growth cycle diagram is shown in Figure 9. It can be intuitively observed in Figure 9 that the transition period of the tea is nearly 20 d, nearly half of the entire process. This is because the transition period requires preparation of the nutrients for subsequent growth of the tea shoots. The next stage is the budding period, during which the bud gradually changes from dark brown to light green, and the bud length changes continuously. In the final stage of the budding period, the bud begins to branch and grow rapidly. Tea shoots then grow into the first stage of the germination period (a stage), where many top spring teas of greater value originate; the tea production yield is small, and the picking operation is difficult in this stage. Stages b and c are the stages of one shoot with two leaves and one shoot with three leaves, respectively. From a morphological perspective, the terminal buds in these two stages are very similar; the change trends of α bud and α stem in the c stage are essentially the same. In stage c, L stem and L bud increased with increasing growth time. In practical production, to increase production yield, workers often pick tea shoots in stage c.
It is concluded from the figure that the c stage is 5-6 d. The appropriate picking time is near the end of stage c, which is approximately 1-2 d. Most of the shoots stopped branching when they had three leaves; the tea leaves became darker, and the stems gradually changed from green to brown, finally growing into folium dichotomous.

Number of Hidden Layers and Momentum Coefficient
The mean square error (MSE) for different momentum coefficients with different numbers of hidden layers and neurons was used as the index to determine the optimal number of hidden layers and the momentum coefficient. As shown in Table 4, the lowest MSE value (0.033 × 10 −4 ) was found with a momentum coefficient of 0.4, when the number of hidden layers was two, and the number of neurons was seven; the ANN model achieved the best training performance with these parameters. The training test was conducted with 1000 training epochs and a learning rate of 0.1. To determine the optimal number of training epochs and neurons in the hidden layer, the MSE for different numbers of training epochs (300-1500) with different numbers of hidden layer neurons were investigated; the results are shown in Figure 10. In the figure, the optimal training performance (MSE = 0.079 × 10 −4 ) was found with 1200 training epochs and 11 neurons. Thus, the optimal configuration of the ANN model included two hidden layers, 11 neurons, a momentum coefficient of 0.4, and 1200 training epochs; the optimal architecture of the established model was 5-2-11-3. A similar methodology was reported by Fang et al. [43] for ANN modeling of the corn drying process in an industrial drying system and by Chokphoemphun et al. [39] for modeling of the moisture content prediction of paddy drying in a fluidized-bed drier with a vortex flow generator.  The growth of tea shoots is affected by different factors; the most influential is "top advantage" [40]. Many terminal buds branch and grow continuously without human intervention, greatly impacting the growth and development of side branches [41]. In actual production, to produce a more uniform growth of tea shoots and to facilitate precision picking, the "top advantage" of tea trees should be removed [42]. The growth cycle diagram presented in this study may only be aimed at the growth of Yinghong #9 tea in autumn, but to some extent, the basic understanding of this tea has been realized for the first time. The growth cycle diagram may provide a theoretical foundation and a technical reference for guiding picking and managing tea shoots in tea gardens and is a highlight of this study.

Number of Hidden Layers and Momentum Coefficient
The mean square error (MSE) for different momentum coefficients with different numbers of hidden layers and neurons was used as the index to determine the optimal number of hidden layers and the momentum coefficient. As shown in Table 4, the lowest MSE value (0.033 × 10 −4 ) was found with a momentum coefficient of 0.4, when the number of hidden layers was two, and the number of neurons was seven; the ANN model achieved the best training performance with these parameters. The training test was conducted with 1000 training epochs and a learning rate of 0.1. To determine the optimal number of training epochs and neurons in the hidden layer, the MSE for different numbers of training epochs (300-1500) with different numbers of hidden layer neurons were investigated; the results are shown in Figure 10. In the figure, the optimal training performance (MSE = 0.079 × 10 −4 ) was found with 1200 training epochs and 11 neurons. Thus, the optimal configuration of the ANN model included two hidden layers, 11 neurons, a momentum coefficient of 0.4, and 1200 training epochs; the optimal architecture of the established model was 5-2-11-3. A similar methodology was reported by Fang et al. [43] for ANN modeling of the corn drying process in an industrial drying system and by Chokphoemphun et al. [39] for modeling of the moisture content prediction of paddy drying in a fluidized-bed drier with a vortex flow generator.

Training Results of the ANN Model
To investigate the fitting performance and the predictive capability of the established 5-2-12-3 ANN model, 134 sets of experimental data ( Table 2) were used to train the ANN; 70% were used for the training, 15% were used for the test, and 15% were used for the validation. The training results are shown in Figure 11. The best training performance was found with 16 training epochs; the MSE was 0.0923 × 10 −4 , and the R values for the training data, test data, and validation data were 0.99976, 0.99871, and 0.99857, respectively. The training results of the established 5-2-12-3 ANN model demonstrated a perfect fitting performance and can be used to predict the picking indexes ξ, r, and Zstem. The weights of the

Training Results of the ANN Model
To investigate the fitting performance and the predictive capability of the established 5-2-12-3 ANN model, 134 sets of experimental data ( Table 2) were used to train the ANN; 70% were used for the training, 15% were used for the test, and 15% were used for the validation. The training results are shown in Figure 11. The best training performance was found with 16 training epochs; the MSE was 0.0923 × 10 −4 , and the R values for the training data, test data, and validation data were 0.99976, 0.99871, and 0.99857, respectively. The training results of the established 5-2-12-3 ANN model demonstrated a perfect fitting performance and can be used to predict the picking indexes ξ, r, and Z stem . The weights of the ANN model are presented in Table 5.

Validation of the Established ANN Model
To validate the prediction capacity of the established ANN model, 46 sets of experimental data collected during the tea picking period were used for validation. The statistical parameters introduced in Section 3.2 were used to estimate the prediction capacity of the model. The validation results for ξ, r, and Zstem are shown in Figure 12a-c, respectively. The statistical results of the 5-2-12-3 ANN model are presented in Table 6. In Figure 12, the predicted values of the picking indexes (ξ, r, and Zstem) are in good agreement with the experimental data, indicating that the established ANN model has a good prediction capacity for ξ, r, and Zstem. According to the statistical results in Table 6, the determination coefficients (R 2 ) for the indexes are all greater than 0.99, and the mean square errors (MSE) for ξ, r, and Zstem are 1.35183, 0.12726, and 0.27117, respectively; the mean absolute errors (MAE) for ξ, r, and Zstem, are 0.41474, 0.03799, and 0.13437, respectively. The statistical results indicate that the established 5-2-12-3 ANN model has a good prediction capacity and can be used to predict tea picking indexes.

Validation of the Established ANN Model
To validate the prediction capacity of the established ANN model, 46 sets of experimental data collected during the tea picking period were used for validation. The statistical parameters introduced in Section 3.2 were used to estimate the prediction capacity of the model. The validation results for ξ, r, and Z stem are shown in Figure 12a-c, respectively. The statistical results of the 5-2-12-3 ANN model are presented in Table 6. In Figure 12, the predicted values of the picking indexes (ξ, r, and Z stem ) are in good agreement with the experimental data, indicating that the established ANN model has a good prediction capacity for ξ, r, and Z stem . According to the statistical results in Table 6, the determination coefficients (R 2 ) for the indexes are all greater than 0.99, and the mean square errors (MSE) for ξ, r, and Z stem are 1.35183, 0.12726, and 0.27117, respectively; the mean absolute errors (MAE) for ξ, r, and Z stem , are 0.41474, 0.03799, and 0.13437, respectively. The statistical results indicate that the established 5-2-12-3 ANN model has a good prediction capacity and can be used to predict tea picking indexes.  Based on the established 5-2-12-3 ANN model, an application scheme is proposed, as shown in Figure 13. During the picking process, the input neurons including Lbud, Lstem, αbud, αstem, and t are measured by the 3-D coordinate instrument, and the measured datasets are used to train the ANN model, the predicted indexes (output neurons) including ξ, r, and Zstem are transferred to an intelligent optimized controller, and further inputted into the tea picking system after compensation. The output signals are executed by an actuator and achieve the goal of precision tea picking.  Based on the established 5-2-12-3 ANN model, an application scheme is proposed, as shown in Figure 13. During the picking process, the input neurons including L bud , L stem , α bud , α stem , and t are measured by the 3-D coordinate instrument, and the measured datasets are used to train the ANN model, the predicted indexes (output neurons) including ξ, r, and Z stem are transferred to an intelligent optimized controller, and further inputted into the tea picking system after compensation. The output signals are executed by an actuator and achieve the goal of precision tea picking.