• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Innovative deep learning artificial intelligence applications for predicting relationships between individual tree height and diameter at breast height

    2020-07-16 07:19:14lkerErcanl
    Forest Ecosystems 2020年2期

    ?lker Ercanl?

    Abstract

    Keywords: Artificial intelligence,Prediction, Deep learning algorithms, Individual tree height

    Introduction

    The significant components of forest inventory, which is the first phase of forest planning, are the measurement of the individual tree heights (ITH) and the diameter at breast height (DBH). These individual tree attributes are used to predict total and merchant volume and biomass,forest site index, especially for uneven-aged stand, and also these attributes have the roles of significant input and independent variable in yield and growth models(Vanclay 1994; Kv and Hui, 1999). The measurements of the individual tree heights are more difficult and time consuming than those of DBH (Huang et al. 1992;Martin and Flewelling 1998) and so the ITH of all trees in sampling units cannot be measured in forest managements (Loetsch et al. 1973; Van Laar and Ak?a 2007).The ITH, whose could not be measured in forest inventory applications, can be predicted by the stand height curves which show the statistical relationships between the ITH and DBH (Avery and Burkhart 1983; Van Laar and Ak?a 2007).

    In forest biometric studies, the empirical relationships between the ITH and DBH are represented by the statistical equations and these relations are modelled by using Nonlinear Regression Models (NLRM) owing to the sigmoid or “S” shaped trend to be evident to these ITH and DBH relations (Wykoff et al. 1982; Huang et al. 1992;Robinson and Wykoff 2004). In forest areas, various stand growing conditions with stand age, site quality and stand stocking have significant effects on the relationships between the ITH and DBH. Thus, the models with only an independent variable such as DBH can remain incapable in successful and effective predicting these relations. In this regard, different statistical prediction techniques have been proposed and used in modelling these relationships between the ITH and DBH that were sampled from different stand growing structures. Ferguson and Leech (1978), Krumland and Wensel (1988),Larsen and Hann (1987) and Parresol (1992) proposed an approach which comprises the prediction of parameter values of these regression models separately for different stand structures at first phase and subsequently developed linear regression models for the relationships between the parameters of this regression model and some stand attributes such as stand age, site index and stocking index at second phase. As a more common approach, the multivariate nonlinear regression models which comprise various stand attributes such as stand basal area, site index, stand age or stocking index in addition to the DBH were developed by various studies such as Huang et al. (2000), Sharma and Zhang (2004),Temesgen and Gadow (2004),Dorado et al.(2005),Trincado et al. (2007), Adame et al. (2008), Paulo et al.(2011). These multivariate ITH models with supplemental stand attributes are also called as “generalized heightdiameter models”. Nanos et al. (2004) analyzed the spatial pattern of the height models and offered the“geostatistical” modelling.

    Another statistical modelling technique that has been used widely to predict the ITH in forestry literature is the Nonlinear Mixed Effect Regression Modelling Approach.This regression modeling technique has been frequently used in modeling empirical relationships between the ITH and DBH,because hierarchically correlated data with clustered and hierarchically sample plots that have been measured to develop the ITH-DBH models may cause serious fitting problem in modeling these relations (Dorado et al.2006; Sharma and Parton 2007). These hierarchical data structures can be evident in the sample plots measured from the stands with different growing structures owing to different stand site quality, stocking and stand age(Calama and Montero 2004; Budhathoki et al. 2008).These highly correlated data violated the assumption of independence of data which is one of the basic assumptions in developing the regression models. The violation of this assumption is called as “autocorrelation” or “serial correlation” (Littell et al. 1996; Lappi 1997). The usage of approaches of nonlinear regression models, especially for the hierarchical data structures, causes biased predictions of the confidence intervals of model parameters in regression models (Searle et al. 1992; Grégoire et al. 1995).This situation negatively affects the reliability of the results of the regression models and as a result incorrect results can be obtained in the height predictions. Especially, in forest studies including the development of the ITH models,the Nonlinear Mixed Effect (NLME) Regression Models have been commonly proposed and used as a solution to deal with this“autocorrelation”problem(Calama and Montero 2004; Meht?talo 2004; Lynch et al. 2005; Dorado et al.2006; Sharma and Parton 2007; Trincado et al. 2007;Adame et al. 2008; Budhathoki et al. 2008; Crecente-Campo et al.2010).

    Beside these statistical modeling techniques with the NLRM and NLME, Artificial Neural Networks (ANN),which is a part of Artificial Intelligence (AI), have become popular as another modelling methodology for predicting the individual tree and stand yield and growth. Especially, significant studies related to the ANN have been conducted since the beginning of 2000s.Numerous prediction models based on AI, especially ANNs, have been developed for modeling various individual tree and stand attributes such as tree volume(Diamantopoulou 2005a, 2006; ?z?elik et al. 2008; Diamantopoulou and Milios 2010; ?z?elik et al. 2010;Soares et al. 2011; Miguel et al. 2016), tree taper (Diamantopoulou 2005b; Leite et al. 2011; Nunes and G?rgens 2016), tree height (Diamantopoulou and ?z?elik 2012; ?z?elik et al. 2013), tree mortality (Hasenauer et al. 2001), survival model (Guan and Gertner 1991),regeneration establishment and height growth (Hasenauer and Kindermann 2002), bark volume (Diamantopoulou 2005a), biomass prediction (?z?el?k et al. 2017),basal area and volume increment growth model (Ashraf et al. 2013). In addition to many ANN studies, Deep Learning Algorithms (DLA) stand out as another prominent AI technique. Although there are a number of significant ANN studies predicting the yield and growth of tree and stand in forestry literature, the DLA models seem to be an innovative technique in front of forest biometrics since 2010. Especially, DLA can be successfully used in analyzing the data clouds (structures which consist of millions or billions of data) and in data mining. The DLA models are basically multi-layer ANN models with at least 3 hidden layers, and this artificial intelligence technique tries to approach the learning and decision-making capacity of the human brain to a certain extent with its complex structure that can contain 5-10 or tens of layers and hundreds and thousands of neurons. Although there is a certain number of studies consisting of ANN for forest yield and growth predictions in today’s literature, the modelling studies with DLA models are in the beginning phase. With the development of computer systems which consist of highly effective graphic processing units, DLA models become more applicable and accessible in today’s world. Its application examples such as the diagnosis of plant illnesses and plant specification have been conducted in agriculture areas (Lee et al. 2015; Mohanty et al. 2016;Sladojevic et al. 2016; Carranza-Rojas et al. 2017; Sun et al. 2017; Ferentinos 2018; Ubbens et al. 2018). Furthermore, there is a need to address the evaluation of new AI techniques for investigating the capability and obtainability in predicting tree and forest attributes that have been important in forest management applications.According to the knowledge of the forest biometric studies including growth and yield models, no studies have been achieved to develop the DLA models to predict individual tree attributes, especially tree height and so the issue of the capability of DLA in predicting tree attributes has been uncertain and needs to be clarified.By widespread of AI techniques such as the DLA models, these scientific evaluations based on the comparative methods have been received remarkable interest and require further modelling studies in forest literature.In this study, it is aimed to evaluate the capability of the usability of the DLA models in predicting empirical relationships between the ITH and DBH as a leading and innovative application. To that end, (1) the DLA models in order to predict relationships between the ITH and DBH measured from stands with different growing structures were trained, (2) the success status of these predictions obtained from DLA models was compared with those of nonlinear regression (NLRM) models,nonlinear mixed effect regression (NLME) and artificial neural network (ANN) models, and (3) the ideal and optimal DLA model structure in the prediction of the ITH was decided by comparing the DLA network structures with various numbers of layers and neurons alternatives.Thus,this study presents scientific the clarification about the issue of whether DLA models can be evaluated as an alternative technique for statistical methods in predicting the individual tree height.

    Materials and methods

    Materials

    In this study, the research material was the data obtained from temporary 150 sample plots as a result of measuring of stands of even aged and pure Anatolian Crimean Pine(Pinus nigra J.F.Arnold ssp.pallasiana(Lamb.)Holmboe)in Konya Forest Enterprise. The studied Pure Anatolian Crimean Pine Stands covered Ak?ehir,Ilg?n and A?a??cigil Forest district areas (Fig.1). This tree species is the most common and dominant species in this region, and so this species of Anatolian Crimean Pine was selected to model relationships between the ITH and DBH as the particular tree species. The characteristic of these studied stands is even aged and pure forest stands with the dominant species of Anatolian Crimean Pine. The altitudes of studied area varied from 250 to 1050 m and the slope ranged between 5% and 60%. The areas studied were characterized geomorphologically as being high mountainous land with moderate and steep slopes.The mean annual temperature is between -5.8°C and 24.8°C, respectively. The climatic regime is a typical a semi-arid continental climate characterized by hot, dry summers and cold, snowy winters.Most of the region usually has low precipitation throughout the year. The mean annual rainfall varies from 400 to 850 mm with a relatively homogeneous precipitation.

    These sample plots were selected by random sampling in terms of different stand age, site quality and density.The sample plots were in the shape of a circle and their size varied from 400 to 800 m2depending on the structure of the stand. At each sample plot, DBH was measured to 0.1 cm precision using calipers at every living tree with a DBH >8 cm. Individual Tree height (ITH)was measured in a subset of trees, selecting two-three trees for each of 4 cm diameter class using Blume-Less Altimeter (0.1 m precision). In addition, the ITH and DBH measurements were obtained from dominant and co-dominant trees, which were selected based on the 100 dominant and co-dominant highest trees per unit area (e.g. Four highest trees in a 0.04-ha plot).

    Totally, 2024 pairs of height-diameter measurements were obtained by the measurement,which realized in these sample plots.These data were divided into two groups randomly in order to use in the training of DLA and ANN models and in developing NLRM and NLME models (1st group data set) and in validating ITH predictions obtained by these methods (2nd group data set). There are approximately 85% (1720 sample trees) of the total data in the 1st group and approximately 15% (304 sample trees) of the total data in the 2nd group. Various statistical information related to the data is provided in Table 1.

    Methods

    Nonlinear regression models (NLRM)

    In order to model the empirical relationships between the ITH and DBH obtained from the stands in different growing structure, various regression models including various stands attributes, further to DBH, have been proposed and used (Huang et al. 1992; Fang and Bailey 1998; Peng 1999;Temesgen and Gadow 2004). Peng et al. (2001) expressed that some model attributes such as the number of parameters, the biological explanation and the validity of model prediction of these models. The model which was chosen to model the relationships between the ITH and DBH possesses some mathematical characteristics such as(i)monotonic increment (ii) inflection point and (iii) horizontal asymptote (Peng et al. 2001). Therefore, seven commonly used functions (M1, M2, M3, M4, M5, M6 and M7) were selected to model the relationships between the ITH and DBH and develop generalized height-diameter model(Table 2). These tested ITH-DBH functions which were proposed by Meyer (1940) modified by Ca?adas et al.(1999) (M1), Loetsch et al. (1973) modified by Ca?adas et al.(1999)(M2),Prodan(1965)modified by Tomé(1989)(M3), Hui and Kv (1993) (M4), Soares and Tomé (2002)(M5), Richards (1959) modified by Sharma and Parton(2007) (M6) and Schnute (1981) modified by Dorado et al.(2006) (M7), have desirable characteristics such as asymptotic with inflection point models.They are biologically reasonable and can provide biological growth curves. These ITH-DBH functions were chosen owing to their desired properties and to commonly preferred in numerous studies modeling relationships between ITH and DBH.

    The nonlinear functions were fitted using the 1st group data set (1720 trees). Based on the Nonlinear Least Squares (NLS), which uses the Levenberg-Marquardt algorithm, the parameters of ITH-DBH functions waspredicted by using NLS package available in the R statistical environment(R Development Core Team 2018).

    Table 1 Summary statistics of the sample trees originated for 1st group and 2nd group data set

    Nonlinear mixed effect (NLME) regression models

    To deal with this “autocorrelation” problem originating from the hierarchical data structures, a Nonlinear Mixed Effect (NLME) modeling procedure was applied to the best predictive height-diameter model by simultaneously predicting both fixed and random parameters.Different from the NLRM, the model parameters of the NLME are divided into two groups as fixed effects and random effect parameters in its model structure. While the fixed effect parameter reveals ITH trend which is common to overall stands, the random effect parameter represents the variance between the stands and defines the variability in the relationships between the ITH and DBH along various stands (Lappi 1997; Calama and Montero 2004; Meht?talo 2004; Dorado et al. 2006;Crecente-Campo et al. 2010).

    The procedure of NLME package available in the R statistical environment, which is based on the Maximum Likelihood method, was used to obtain the parameter predictions of NLME that presents the best predictive height-diameter model. To decide the best predictive random-fixed parameter alternative for model structure based on NLME, the ITH-DBH models including one,two or three random parameters were fitted and compared based on some statistical comparison criteria. The adaptive Gaussian quadrature was used in the computation of the integral over the random effects as described by Pinheiro and Bates (2000). Furthermore, this NLME procedure was performed assuming the homogenous within-tree variance and uncorrelated residuals.

    Artificial neural network models

    Being as Artificial Intelligence (AI) prediction technique,Artificial Neural Network (ANN) based on the Feed Forward Backprop (FFB) and Cascade Correlation (CC)training algorithms with training function of Levenberg-Marquardt were used to model the relationship between ITH and DBH. These training algorithms including the Feed Forward Backprop (FFB) and Cascade Correlation(CC) have commonly been used to predict tree and forest attributes in forest literature. The reason for choosing these training functions from different training algorithms is its intensive use in forestry. When training the ANN models with FFB and CC algorithms, the individual tree height values, ITH were predicted as target variable. In these ANN models, DBH and the best predictor variables selected from preliminary analyses including a trial and error procedure using different combinations of these stand attributes, such as basal area, number of trees of sample plots, quadratic mean diameter, the dominant DBH and ITH of the sample plot, were used as input variables. The standard ANN models can include three layers such as input layer, hidden layer and output layer. Especially, the activation functions including Hyperbolic tangent sigmoid (tansig), logistic sigmoid function (log-sig) and linear function (Pure-lin) connect these network layers in ANN structures. These activation function alternatives have significant effects on fitting performance of neural network. In this study, alternatives including some activation functions in the connection between input, hidden and output layers were compared to decide the best predictive one: (A1) tan-sig function between input layer and hidden layer and tan-sig function between hidden layer and output layer, (A2) tan-sig function between input layer and hidden layer and log-sig function between hidden layer and output layer, (A3) tan-sig function between input layer and hidden layer and pure-lin function between hidden layer and output layer, (A4) log-sig function between input layer and hidden layer and logsig function between hidden layer and output layer, (A5)log-sig function between input layer and hidden layer and tan-sig function between hidden layer and output layer, (A6) log-sig function between input layer and hidden layer and pure-lin function between hidden layer and output layer, (A7) pure-lin function between input layer and hidden layer and pure-lin function between hidden layer and output layer, (A8) pure-lin function between input layer and hidden layer and log-sig function between hidden layer and output layer and (A9) pure-lin function between input layer and hidden layer and tansig function between hidden layer and output layer.Other important parameter of the network structure is the number of neurons in hidden layer. Thus, some alternatives for the number of neurons which ranged from 1 to 100; 1, 2, 3, …, 20, 30, 50, 70, 90 and 100 number of neurons were compared to select the best predictive neuron alternative in this study. As a result, a total of 900 network alternatives including 100 number neurons and 9 transfer function alternatives (100×9=900 alternatives) based on the Feed Forward Backprop (FFB)and Cascade Correlation (CC) training algorithms, totally 1800 alternatives for FFB and CC-ANN models,were trained and used to obtain the ITH predictions.Being as other significant parameters for ANN structure,the value of 3000 for epochs, the value of 1×10-10for performance goal, the value of 1×10-10for Minimum performance gradient and 1×10-8for epsilon gave the best predictive results to train these FFB and CC-ANN models in the preliminary of this study and so, these parameters were used to obtain the ITH predictions and to compare with those by other predictions methods such as NLRM, NLME and DLA models. These network trainings for 1800 network alternatives for FFB and CCANN models were carried out using newff syntax for the feed-forward backpropagation network and newcf syntax for the cascade-forward backpropagation network codded in MATLAB software (MATLAB 2014).

    Table 2 The nonlinear ITH-DBH functions tested by this study

    Deep learning algorithms

    Deep Learning Algorithm (DLA) models are an artificial intelligence technique which has remained on the agenda since 2010. The DLA has shown quite successful results in various applications such as image classification, video analysis, speech recognition, natural language learning process in recent time. The Artificial Neural Network (ANN) models which is another Artificial Intelligence (AI) type have been usually developed to the input layer, hidden layer (two hidden layers in some cases), output layers in its model structure. However, the DLA models have a quite complex structure comprising many (5, 10 or even tens of) hidden layers compared with ANN models. Especially, the use of Graphics Processing Units (GPU) of computers in the training of DLA models provides this DLA model to be more accessible and usable with effective and successful results in various applications, especially visual and speech recognition by modern day computer technologies, that have not been seen in the history of humanity. In addition to all these successful and efficient use of DLA models in computer systems, the use of DLA models in forestry applications, specifically some tree and forest attributes prediction practices,has been quite limited.

    As the calculations in applications of DLAs are quite complex and intense, obtaining the predictions for tree and forest attributes by DLA model requires intensive use of computer software. Despite some DLA applications and platforms were developed in various languages,the H2O package (R Development Core Team 2018),which has been coded in R software language, becomes prominent in terms of its characteristics such as userfriend, the ability of finding a solution and comparison of different network alternatives. Basically, the H2O package, which operates on the R software platform, is an open source coded artificial intelligence library and comprises different artificial intelligence applications such as “Generalized Linear Models”,“Gradient Boosting Machines”, “Random Forests”, “Deep Neural Networks(Deep Learning)”, “Stacked Ensembles”, “Naive Bayes”,“Cox Proportional Hazards”, “K-Means”, “PCA” and“Word2Vec” (R Development Core Team 2018).

    In this study, the H2O package was used to train the network models based on the DLA models which predicted the individual tree height values, ITH (target variable). In order to determine the input variables in DLA model structure, the trial and error method were used by comparing some alternatives including various independent variables such as the DBH and stand attributes similar to variable determination method in ANN models. Also, the network parameters such as number of layers, number of neuron and type of transfer function are important attributes that affect the success of prediction results in obtaining the predictions with the DLAs in these training DLAs. From various transfer functions, the “Rectifier” function was selected as a transfer function in DLAs’ structure owing to its successful fitting results in our preliminary analyses. The H2O package uses the adaptive learning rate algorithm(ADADELTA) in the trainings of DLA (Zeiler 2012).The rho describes the rate of ADADELTA and epsilon expresses learning rate for DLA models. In the present study, the value of 0.999 for rho and 1×10-8for epsilon were used to train DLA models. Also, the value of 1000 for the epochs, the number of iterations to be accepted in training networks, was used in the training of DLA models, since the best predictive results have been obtained with 1000 in various neural network studies. As a training algorithm, the Gradient Descent Function with the Gaussian distribution model based on the Mean Squared Error function type was used.

    In addition to these parameters of DLAs, the number of hidden layers and the number of neurons in these hidden layers are network parameters that need special attention in training DLA models. In training DLA models, 8 numbers of hidden layers (with 3, 4, 5, 6, 7, 8,9 and 10 layers) starting from 3 layers (that is the minimum number of layers of DLA) to 10 layers and 10 different neuron alternatives ranging from 10 to 100 by increasing 10 at each step (10, 20, 30, 40, 50, 60, 70, 80,90 and 100 neurons) were considered as important network parameters. Thus, 80 different DLA models, 8 different numbers of hidden layers and 10 different numbers of neuron alternatives, were trained to obtain the predictions of ITH.

    The K-Fold Cross Validation method was used in trainings of DLA models,because this method“cross validation up to k number”may reduce“overfitting errors”in obtaining the predictions by the network models. In this study,the value of“cross validation up to k number”was applied on the basis of as(nfolds=10)with the“nfolds”parameter of the H2O.ai Team package (R Development Core Team 2018).

    Comparison criteria

    In this study, various statistical fitting criterion values were used to compare and evaluate the predictions of ITH that were obtained by the NLRM, NLME, FFBANN and CC-ANN and DLA models. These fitting criteria are (1) average absolute error (AAE), (2) the maximum absolute error (max. AE), (3) the root mean squared error (RMSE), (4) % root mean squared error(RMSE%), (5) the average Bias (Bias), (6) % average Bias(Bias%), (7) the fit index (FI), (8) Akaike’s information criterion (AIC) and (9) Bayesian information criterion(BIC).These criteria are calculated as follows:

    This study has carried out two-level comparisons in evaluating many prediction models including NLME models including one (five alternatives), two (ten alternatives) and three (ten alternatives) random parameters(non-convergence fitting results for four and five random parameter alternatives) for the best predictive function from seven ITH-DBH functions tested, 900 # FFBANN models and 900 # CC-ANN models including 100 number neurons and 9 transfer function alternatives and 80 # DLA models including 8 different numbers of hidden layers and 10 different numbers of neuron alternatives. This two-stage evaluation process was carried out to determine the best predictive one from different prediction methods: (1) Firstly, the performance of ITH predictions obtained by the NLME (different random and fixed effect parameter alternatives), DLA (80 different models), FFB-ANN (900 different models) and CCANN (900 different models) methods were compared based on the Relative Rank Values proposed by Poudel and Cao (2013) within each prediction methods, (2) in the second stage, the best predictive model alternative at each AI model level such as DLA, FFB-ANN and CCANN, and NLME including various random and fixed effect parameters were compared with those by NLR.Thus, it has been possible to evaluate about 1900 various model alternatives obtained by various modeling techniques such as NLRM, NLME, DLA, FFB-ANN and CCANN and to determine the best predictive model.

    The validation of prediction methods

    Results

    In first level comparisons, the best predictive models from NLME,FFB-ANN,CC-ANN and DLA model alternatives were selected based on the Relative Rank Methods proposed by Poudel and Cao (2013) and these best predictive models of different prediction methods with NLME, FFB-ANN, CC-ANN and DLA were evaluated with NLRM models in second level comparison. As input variables in these AI models, the stand attributes including diameter at breast height (cm), the dominant height (h0, m), and (cm) (Dg) gave best predictive fitting results in FFB-ANN, CC-ANN models. Nevertheless,DLA models including diameter at breast height (d, cm),the dominant height (h0, m) and dominant diameter (d0,cm) as predictor variables resulted in best predictive ITH from various input variables.

    As a result of the second level comparison, the fitting criteria of AAE, max. AE, RMSE, RMSE%, Bias, Bias%,FI, AIC and BIC for various prediction methods with NLRM, NLME, FFB-ANN, CC-ANN and DLA models are given in Table 3. The relative rank values (Poudel and Cao 2013) related to these goodness-of-fit criteria values and total relative rank values were shown in Table 4. In these fitting criterion values, RMSE ranged from 0.5575 to 0.8306, RMSE% ranged from 4.9504% to 7.3750%, AIC ranged from -998.9540 to -313.3060,BIC ranged from 884.6591 to 1570.3072, FI ranged from 0.8749 to 0.9436, AAE ranged from 0.4077 to 0.6170,max. AE ranged from 2.3696 to 4.4859,Bias ranged from-0.0006 to -0.2695 and Bias% ranged from -0.0050%to -2.3927%. From various ITH-DBH functions tested,the function of Soares and Tomé (2002), M5, gave the best predictive fitting results a RMSE value of 0.7621,RMSE% value of 6.7672%, the AIC value of -461.2447,BIC value of 1422.3685, FI of 0.8946 values, AAE value of 0.6132, max. AE value of 3.8927, Bias value of -0.005

    and Bias% value of -0.0047%. This best predictive ITHDBH function, M5, was selected to apply NLME procedure including some random-fixed evaluations based on the relative rank values. From these different randomfixed effect parameter alternatives, the NLME model of M5 with one random parameters, f parameter, and four fixed effect parameters, resulted in the best predictive fitting statistics with a RMSE value of 0.7073, RMSE%value of 6.2807%, the AIC value of -589.567, BIC value of 1294.0461, FI of 0.9092 values, AAE value of 0.5769,max. AE value of 3.4091, Bias value of -0.0006 and Bias% value of -0.005%. According to these results obtained with NLRM and NLME prediction models, FFBANN and CC-ANN models have partially improved the prediction success. In ANN models, CC-ANN based on A3 activation function alternative and 85 # neuron and FFB-ANN based on A3 activation function alternative and 73 # neuron gave the best predictive fitting results with a RMSE values of 0.7110 and 7160, RMSE% values of 6.3132% and 6.3576%, the AIC values of -580.6896 and-568.6347, BIC values of 1302.9236 and 1314.9787,FI values of 0.9083 and 0.9070, AAE values of 0.5638 and 0.5775, max. AE values of 2.9023 and 2.5863, Bias values of -0.0144 and-0.0175 and Bias% values of -0.1282% and-0.1557%, respectively. Nonetheless, the DLA models showed better predictive performance in explaining the variation in ITH and resulted lower total ranks (ranging from 11.073 to 57.538) than those by NLRM, NLME, FFB-ANN, CC-ANN (ranging from 83.736 to 176.923). In Tables 3 and 4, the results of DLA models with the best predictive number of neuron alternative according to each hidden layer choices from 80 various DLA models was presented. On the basis of the total relative rank values for these prediction methods, the DLA model structure with 9 hidden layers and 100 neurons showed the better predictive results in the prediction of the ITH than those by the other prediction models (Table 4). This DLA model structure alternative has a significant predictive ability with a RMSE value of 0.5575, RMSE% value of 4.9504%, the AIC value of-998.9540,BIC value of 884.6591,FI of 0.9436 values,AAE value of 0.4077,max.AE value of 2.5106,Bias value of 0.0057 and Bias% value of 0.0502% compared to other prediction models.

    Table 3 The goodness-of-fit statistics r, AAE, max.AE, RMSE,RMSE%,Bias, Bias%, FI,AIC and BIC for the best predictive DLA models with best predictive number of neuron alternative according to each hidden layer choices,the ITH-DBH functions based on NLRM,M5 based on NLME with f random,FFB-ANN and CC-ANN

    Table 4 The relative rank values of r, AAE, max.AE, RMSE,RMSE%,Bias, Bias%, FI,AIC and BIC for the best predictive DLA models with best predictive number of neuron alternative according to each hidden layer choices,the ITH-DBH functions based on NLRM,M5 based on NLME with f random,FFB-ANN and CC-ANN

    Figure 1 showed the relationships obtained between observed and predicted height values by network models including (a) the M5 based on NLRM, (b) the M5 based with f random on NLME,(c)FFB-ANN based on A3 activation function alternative and 85#neuron,(d) CC-ANN based on A3 activation function alternative and 73 #neuron,(e)DLA with 100#neurons in nine hidden layers.When these graphs were examined (Fig.2), it is seen that the best predictive DLA network model (DLA with 100 #neurons in nine hidden layers) evidenced more correlated relationships between predicted and measured values around the 1:1-line than those for other prediction models with NLRM, NLME, FFB-ANN and CC-ANN. Thus, ITH predictions which were obtained by this best predictive DLA network model more precise than those of other prediction methods including NLRM, NLME, FFB-ANN and CC-ANN. These graphical results about predictive ability of this best predictive DLA network model were propped with the relationships between these residual and prediction values which were presented in Fig.2. This graph(Fig.3) presented random trend of residual around zero and no important relations, suggesting that there is no serious failure of homoscedasticity, violations of the assumption of constant variance, for those by this best predictive DLA model. For a further analysis of residuals of the best predictive DLA, NLRM, NLME, FFB-ANN and CC-ANN models, Fig.4 presented the plot of residuals against lagged residuals by (a) the M5 based on NLRM,(b) the M5 based on NLME with f random, (c) the FFBANN based on A3 activation function alternative and 85#neuron,(d)the CC-ANN based on A3 activation function alternative and 73 # neuron, (e) DLA with 100 # neurons in nine hidden layers. This plot shows a significant autocorrelation in residuals from the ITH prediction by NLRM of M5 function. A moderate improvement was obtained in predictions with the NLME of M5 function including f random parameters. This improvement about autocorrelation quite clearly obtained by this best predictive DLA model, give no trends in the lag-residuals, suggesting that no-autocorrelation problem was the case for the height predictions by this network model(Fig.3e).

    In this study, it was pointed out the effect of alternatives for different numbers of hidden layers and neurons on the fitting ability of the ITH predictions and so judged the ideal and optimal DLA model structure in these predictions. The results related to this evaluation are presented as the average fitting criteria of RMSE,RMSE%, AIC, BIC, FI and AAE according to the alternatives for different numbers of layers and neurons in Tables 5 and 6. When the changes of these fitting criterion values according to the alternatives for the number of layers and neurons, it was seen that there was a progress in the criteria values from the 3rd layer to 8th layer generally; however, there was a worsening in 7th, 9th and 10th layers in these prediction success values. On the other hand, it is seen that the increase in the number of neurons causes a general improvement in these fitting criteria,which except the number of 50 and 90 neurons.

    The present study validated the NLRM, NLME, FFBANN, CC-ANN and DLA models to the independent data set by using “Equivalence” test and the results related to this test were shown in Table 7. Consistent with these analysis results, the h0hypothesis which pronounce that “the constant is different from 10.8421 cm(the average observed ITH values) and the slope coefficient (b1) (except the DLA model with 8 layers) is different from 1” has been rejected. Thus, it can be concluded that the aforementioned DLA models (except the DLA with 8 layers) can be accepted and used statistically 95%in the ITH predictions of the stands in the study areas.Also, the fitting criteria values related to the prediction obtained with different DLA models from these 304 trees are shown in Table 8.

    Discussion

    This study is the first attempt to model individual tree height-diameter relationships by using Deep Learning Algorithms (DLA) that have been another application of Artificial Intelligence Techniques. The main topic of this research is the question whether the DLA model, as an alternative, will offer predictive results compared as the classical regression models, which have been in use for many years in modelling the growth of trees, and ANN models, another type of AI technique. In addition, various network alternatives were evaluated to determine the optimal network structure based on the statistical criteria and, for this purpose, 80 different DLA models were trained by using the data collected from different forest stands. When considering the evaluation results based on the Relative Rank Methods (Poudel and Cao 2013) seen in Tables 3 and 4, these DLA models offer better statistical performance than those by the NLRM,NLME, FFB-ANN, CC-ANN and DLA models in the predictions of tree heights. Especially, the DLA network model with 9 layers and 100 neurons resulted in the best predictive tree heights in this study. This DLA network model gave significant improvement in the values of RMSE, AIC, BIC, FI, AAE, max. AE with the rates of 26.85%, 116.58%, 37.80%, 5.48%, 33.52%, 35.51%, respectively, compared as those of NLRM.

    Table 5 The average of fitting criteria of RMSE,RMSE%,AIC, BIC, FI and AAE according to the number of hidden layers

    Table 6 The average of fitting criteria of RMSE,RMSE%,AIC, BIC, FI and AAE according to the number of Neurons

    Considering the predictive capability of ITH obtained by these DLA models,it can be observed that the DLA model with 9 layers and 100 neurons produced higher prediction precisions than those by the NLRM, NLME and FFBANN and CC-ANN (Fig.2), which this DLA model gave the tree height predictions that were very close to the observed ones.Also,the graphical analysis of the scatter plot of the residuals against to predicted heights(Fig.3)shows a uniform distribution around zero with approximately constant variance, indicating that the homoscedastic model provides a good representation of the data. Moreover, this uniform and random distribution in the errors obtained by the DLA model with 9 layers and 100 neurons is more distinct (Fig.2e). When considering the residuals against lagged residuals obtained by the DLA model with 9 layers and 100 neurons and others(Fig.4),it is seen that this DLA network model provides no trends in the lagresiduals (Fig.4e) and more desirable qualities for autocorrelation problem than those by the NLR model. Based on all these results with fitting performance criteria, it is concluded that the DLA network models, especially the network model with 9 layers and 100 neurons, have been considered as an alternative prediction method to traditional regression models such as NLRM or NLME and other AI technique including FFB-ANN and CC-ANN to model individual tree height-diameter relationships. In this research area about modeling height-diameter relations, Brandao (2007), de A Silva et al. (2008), Diamantopoulou and ?z?elik (2012) and ?z?elik et al. (2013)compared Artificial Neural Network models with NLR for predicting tree heights, and these studies found that the ANN is superior to NLRM in terms of many statistical criteria. Similarly, Lee et al. (2015), Mohanty et al. (2016),Sladojevic et al. (2016), Carranza-Rojas et al. (2017), Sun et al. (2017), Ferentinos (2018) and Ubbens et al. (2018)successfully used the DLA to determine plant disease diagnosis in agriculture applications. Beyond all these studies including the ANN models in forestry and the DLA models in the agriculture area, this study presents a first DLA model for predicting the relationships between individual tree height and diameter at breast height that have been an important individual tree measurement in forest inventory. When evaluated the results obtained by the present study, it is seen that the DLA models which are a leading and innovative artificial intelligence technique can be used as an alternative method for regression models whose applications has started in the 1940s such as Metzler (1940), Samuelson (1942), Tintner (1944) and which have problems in providing various statistical assumptions mentioned in many studies nowadays. Although the regression models have provided a certain extent successful fitting results for predicting the relationships between ITH and DBH, the DLA models stand out with some important and attractive features: (1) its strong nonlinear modeling capability without predetermined any statistical functions and(2)no assumptions needed for independence, normal distribution, and homoscedasticity of residuals; and multicollinearity among variables, andspatial and longitudinal autocorrelations in data. In this respect, as an alternative to traditional regression models,the use of DLA models for predicting these ITH-DBH relations and other possible tree and forest attributes can be highlighted.

    Table 7 The results of equivalence tests for the best predictive DLA of the number of neuron alternatives regarding the numbers of hidden layer and M5 function based on NLRM, M5 based on NLME with f random, FFB-ANN,CC-ANN

    Table 8 The goodness-of-fit statistics of DLA in validation data set

    In addition to the satisfactory findings by the DLA to training data, another issue that should be considered is the analysis of the fitting ability in the simulation data group, especially later uses of the trained model, which were not used in the training process.In the simulation of Artificial Intelligence (AI) models to other forest areas or new measured data, the predictive performance may substantially decrease and the “overfitting” problem may occur in the AI applications.In this regard,the analysis of the success status in independent data is an issue which should be given particular importance in the evaluation of the applicability of AI models. In this study, the DLA models were evaluated in terms of “overfitting” problem by using“Equivalence”test in the independent data.When the “Equivalence” test in Table 7 and fitting criteria in Table 8 are evaluated,it is seen that the DLA models provided acceptable results for these independent data and produced the fitting criteria similar to those of training data. These better predictive results obtained for these independent data compared as those for the training data set suggest that the DLA models may not have a problem of“overfitting”.These predictive results of the DLA,especially for independent data, with no “overfitting” problem can be explained by the fact that the DLA models were trained with the appropriate number of iterations to represent successfully the relationships in the data,which detailed information were provided by Ruder (2017). In this regard, the determination of ideal and optimum DLA has a significant effect on not only increasing the predictive ability of DLA models, but also overcoming the “overfitting” problem during the simulation of independent data in the trainings of DLA models.

    In this study, various alternatives with the number of layers and neurons included in network structure were compared and evaluated to decide the optimal network structure for DLA models, because another issue that should be considered in studies about DLA models was the determination of the optimal network structure.While a significant improvement in fitting criteria can be seen in the average of these criteria from 3 to 8 (3, 4, 5, 6, 7 and 8)number of layers,thus upgrading could not be observed in the average success criteria in 7, 9 and 10 layers(Table 5).With respect to increase in number of neurons,from 10 to 100 numbers of neurons, consistent progress in the average success criteria was obtained in general(Table 6). This worsening in the success criteria depending on 7,9,and 10 numbers of layers can be explained by the failure of a DLA model structure to represent and model the height-diameter relations, owing to unsatisfactory solution of parameter values related to a DLA model structure which is complicated by excessive increase of the number of layers.On the other hand,when the change in the success criteria due to the increase in the number of neurons is evaluated, it can be explained that the complex model structure that was formed with the increase in the number of neurons in DLA model structure,even with 100 neurons, do not cause a data representation failure in the parameter values. However, another issue that should be considered is the interaction of layer and neuron number changes in the DLA model structure. When these interactive changes of the number of layers and neurons are evaluated, the best predictive results were obtained with 8 layers (Table 5), nevertheless; the best results predictions were obtained with the DLA model which has 9 layers and 100 neurons due to the mutual interaction of numbers of layers and neurons (Table 3). As it is seen in this study, the number of both layers and neurons,if possible, the other parameters of the DLA model structure,should be evaluated together to decide the ideal and optimal DLA model structures and these evaluations can be carried out by comprising the mutual interaction of these factors. These preliminary findings about the number of layers and neurons for a DLA model structure which were firstly obtained by present study are important results that will make significant contributions to future DLA studies.

    Besides the predictive ability of the DLA models in predicting individual tree height-diameter relationships,some features restricting the applicability of these models should also be taken into consider while evaluating the applicability of DLA models. In general, the regression models where the equation structures and parameter values can be given together are preferred in modelling studies.Also,the DLA models,which are consisted of tens of layers and neurons, can have the model structures which comprise hundreds even thousands of weight values. In this regard, it will not be possible to give the equation structure of the DLA that has many weight values and to use the applications such as excel,etc. Thus, the applications of DLA models are only possible with the support of various computer software and programs, which it comes insight clearly that it will not be very difficult given that we live in the computer era.Especially, the R software platform, which becomes prominent with its applications and usage nowadays, will allow the forest planners initially and other various applicators to use the DLA and various AI models.The applications of DLA models, which were trained by various researchers and applicators, should be prepared in R platform, which is free and open for all, and shared with various stakeholders and other users in forest management.

    The study provided the R syntax file of the best predictive DLA network model with 9 layers and 100 neurons as the supplementary file and the downloadable link from Google Drive Link (https://drive.google.com/open?id=1ewzoB0-0G89rZLkKHVqdkFSLhMjnR9JP) so that other forest practitioner can use this best predictive DLA model,which similar applications were applied for validation data of 304 trees in this study. This DLA model can be downloaded and can be used by future forest practitioner to obtain the ITH predictions for other tree species in other parts of the world.In the use of these best predictive DLA models in other species and areas, it is an important requirement that the tree species and area for future use are similar to the study area in which the species and data included in this study. As this present study has shown by training the DLA models and providing R syntax codes of the best predictive DLA models, artificial intelligence studies should provide more innovate network tools for different users, as well as including comparisons with other classical methods. This study provides a presentation of R syntax code file for artificial intelligence models to give the opportunity to other forest practitioners to use artificial intelligence model developed in this study.

    The data in this study were limited in the sample size of this study is 2024, of which 1720 were used for training, and so the effectiveness and success of artificial intelligence models in modeling big data may not have been obtained sufficiently, or a limited number of data may have negative effects on iteration success. However,while data pools in the forest growth and yield modeling studies such as this study remain limited the sample size,data analysis which may consist of millions or even millions of data, also called as big data, may be involved in applications such as forestry image processing such as Hamdi et al. (2019), Fricker et al. (2019) and Sylvain et al. (2019). In the analysis of forestry image processing data based on big data, the effectiveness of deep learning techniques will be even more apparent.

    This study has introduced innovative Deep Learning Algorithms (DLA), being as another application of Artificial Intelligence Techniques, which were resulted in superior fitting statistics compared as conventional regression models. The weakness of this study is that the fitting results are obtained by modeling only one species form pure stands. However, the future applications of DLA models need to be realized for mixed stands or uneven forest stands. Thus, the acceptability of the results for the DLA models will become even more apparent and the availability of other models can be achieved.However,more scientific studies are needed to compare DLA models with other convenient models.As an artificial intelligence technique, the present study is a preliminary step and contribution to the evaluation process regarding the future usability of deep learning technique and its scientific acceptability.

    Conclusion

    We have been experiencing the fourth Industrial Revolution with the proliferation of the use of artificial intelligence nowadays and the evaluation of the Deep Learning Algorithms, one of the Artificial Intelligence Techniques that has come up since 2010, stands out as an important requirement in forest yield and growth modelling studies.This paper presents the DLA models, as innovative prediction technique,to predict the relationships between individual tree heights and diameter at breast height,which are an important growth parameter of trees and so, the usability and capability of the DLA were evaluated based on some fitting criteria in both training and simulation datasets.The fitting results obtained by the DLA models underlined that the DLA models can be assessed as an alternative prediction method for traditional regression models to obtain individual tree heights in forest inventory. This paper introduces the abilities of the DLA models that have been a novel neural network model in the field of Artificial Intelligence to predict the individual tree heights from the diameter at breast height measured in the sample plots.Besides predictive applications of the DLA models in modelling tree height-diameter relations in this study, the fitting ability and usability of the DLA models should be evaluated in predicting the other individual tree attributes such as tree volume,taper and growth and so stand attributes such as stand volume, basal area, biomass and carbon. It is confronted as an important need that the realization of different studies related to the evaluation of the DLA models being as novel Artificial Intelligence Application, which found a place newly in the forestry literature,as an alternative for conventional statistical methods in predicting various stands and individual tree attributes.

    Supplementary information

    Supplementary information accompanies this paper at https://doi.org/10.1186/s40663-020-00226-3.

    Abbreviations

    AAE: the Average Absolute Error; ADADELTA:the Adaptive Learning Rate Algorithm; AI:the Artificial Intelligence;AIC: Akaike’s Information Criterion;ANN: the Artificial Neural Networks; Bias: the average Bias; Bias%:percent of average Bias; BIC: Bayesian Information Criterion;DBH: the diameter at breast height; DLA: Deep Learning Algorithms;FI: the Fit Index; ITH: the individual total height; max. AE: the Maximum Absolute Error; NLME: the Nonlinear Mixed Effect; NLRM: the Nonlinear Regression Models; NLS:the Nonlinear Least Squares; RMSE: the Root Mean Squared Error; RMSE%: percent of Root Mean Squared Error; TOST: Two One-Sided Test Strategy

    Acknowledgements

    Not Applicable.

    Authors’contributions

    ?lker ERCANLI contributed the manuscript as supervising the work,conception, administrative, design, material support, critical revision and statistical analysis. The author(s) read and approved the final manuscript.

    Authors’information

    ?lker ERCANLI is an associated professor at ?ank?r? Karatekin University, Forest Faculty, Department of Forest Yield Studies. He is a forest biometrician,focused on forest growth modeling.

    Funding

    Not Applicable.

    Availability of data and materials

    Available on request.

    Ethics approval and consent to participate

    Not Applicable.

    Consent for publication

    Not Applicable.

    Competing interests

    The author declare that they have no competing interests.

    Received: 3 August 2019 Accepted: 2 March 2020

    后天国语完整版免费观看| 黄网站色视频无遮挡免费观看| 无限看片的www在线观看| 手机成人av网站| 免费在线观看影片大全网站| 国产一卡二卡三卡精品| 欧美精品一区二区大全| 首页视频小说图片口味搜索| 日日夜夜操网爽| 99国产精品一区二区三区| 精品国产超薄肉色丝袜足j| 99re在线观看精品视频| 在线看a的网站| 国产欧美日韩综合在线一区二区| 国产在线精品亚洲第一网站| 亚洲欧美激情在线| 国产精品久久久人人做人人爽| 黑丝袜美女国产一区| 超色免费av| 可以免费在线观看a视频的电影网站| 国产精品一区二区免费欧美| 国产精品98久久久久久宅男小说| 巨乳人妻的诱惑在线观看| 国产成人精品在线电影| 三级毛片av免费| av天堂久久9| 美女国产高潮福利片在线看| 丁香六月天网| 中文字幕另类日韩欧美亚洲嫩草| 少妇的丰满在线观看| 日日摸夜夜添夜夜添小说| 91精品三级在线观看| 欧美在线一区亚洲| 2018国产大陆天天弄谢| 午夜激情av网站| av线在线观看网站| 免费高清在线观看日韩| 成人黄色视频免费在线看| 男女高潮啪啪啪动态图| 麻豆成人av在线观看| 国产在线观看jvid| 午夜福利在线免费观看网站| 久久这里只有精品19| 日本av手机在线免费观看| 亚洲欧美激情在线| 欧美老熟妇乱子伦牲交| 国产欧美日韩一区二区精品| 国产av国产精品国产| 国产日韩欧美视频二区| 亚洲,欧美精品.| 午夜福利在线观看吧| 婷婷丁香在线五月| 曰老女人黄片| 日韩精品免费视频一区二区三区| 亚洲七黄色美女视频| aaaaa片日本免费| 日日爽夜夜爽网站| 韩国精品一区二区三区| 丝袜美腿诱惑在线| 两人在一起打扑克的视频| 中文字幕av电影在线播放| 激情视频va一区二区三区| 久久天堂一区二区三区四区| 国产黄频视频在线观看| 黑人欧美特级aaaaaa片| 婷婷成人精品国产| 大香蕉久久网| 三级毛片av免费| 99国产精品99久久久久| 91字幕亚洲| 久久久久精品人妻al黑| 欧美日韩成人在线一区二区| 我的亚洲天堂| 999久久久国产精品视频| 中文字幕人妻丝袜一区二区| 免费人妻精品一区二区三区视频| 日本av手机在线免费观看| 亚洲精品在线美女| 大香蕉久久网| 亚洲天堂av无毛| 熟女少妇亚洲综合色aaa.| 露出奶头的视频| 美女高潮到喷水免费观看| av超薄肉色丝袜交足视频| 国产一区有黄有色的免费视频| 日本精品一区二区三区蜜桃| 99精国产麻豆久久婷婷| 国产欧美日韩一区二区三| 国产在线免费精品| 老鸭窝网址在线观看| 无人区码免费观看不卡 | 黄色丝袜av网址大全| 大片免费播放器 马上看| 国产成+人综合+亚洲专区| 黄色视频在线播放观看不卡| 午夜免费成人在线视频| 一区二区日韩欧美中文字幕| 五月开心婷婷网| 亚洲欧美一区二区三区黑人| 久久精品亚洲av国产电影网| 麻豆国产av国片精品| 女人精品久久久久毛片| 欧美另类亚洲清纯唯美| 久久久久久久久免费视频了| 天天躁日日躁夜夜躁夜夜| 激情在线观看视频在线高清 | 欧美 日韩 精品 国产| 午夜精品国产一区二区电影| 丁香六月天网| 国产一区二区激情短视频| 高清毛片免费观看视频网站 | 啦啦啦免费观看视频1| 最近最新中文字幕大全免费视频| 亚洲精品乱久久久久久| 一级黄色大片毛片| 精品乱码久久久久久99久播| 黄色片一级片一级黄色片| 中文字幕高清在线视频| 国产精品久久久av美女十八| 色精品久久人妻99蜜桃| 黄色怎么调成土黄色| 国产欧美日韩一区二区三区在线| 日韩视频一区二区在线观看| 亚洲欧洲日产国产| 蜜桃国产av成人99| 久久国产精品大桥未久av| 波多野结衣一区麻豆| 国产精品欧美亚洲77777| 中文字幕精品免费在线观看视频| 另类精品久久| 国产日韩欧美亚洲二区| 久久久国产欧美日韩av| 男女无遮挡免费网站观看| 国产欧美日韩一区二区三| 午夜福利免费观看在线| 人成视频在线观看免费观看| 别揉我奶头~嗯~啊~动态视频| 国产成人欧美在线观看 | 色尼玛亚洲综合影院| 啦啦啦中文免费视频观看日本| 亚洲中文日韩欧美视频| 精品少妇一区二区三区视频日本电影| 99国产极品粉嫩在线观看| 一进一出好大好爽视频| cao死你这个sao货| 女人精品久久久久毛片| 黄色怎么调成土黄色| 精品少妇一区二区三区视频日本电影| 欧美日本中文国产一区发布| 可以免费在线观看a视频的电影网站| 中文字幕精品免费在线观看视频| 午夜福利欧美成人| 国产黄频视频在线观看| 国产1区2区3区精品| 精品人妻1区二区| 久久国产精品大桥未久av| 国产单亲对白刺激| 一本大道久久a久久精品| 久久国产精品男人的天堂亚洲| 大型av网站在线播放| 久久久精品94久久精品| 亚洲专区中文字幕在线| 精品福利永久在线观看| 日日夜夜操网爽| 女人高潮潮喷娇喘18禁视频| 久久精品成人免费网站| 精品卡一卡二卡四卡免费| 91成人精品电影| 国产成人精品久久二区二区91| 成人18禁在线播放| 久久久久久久国产电影| a级毛片在线看网站| 最近最新中文字幕大全电影3 | 少妇粗大呻吟视频| 国产伦人伦偷精品视频| 怎么达到女性高潮| 美女午夜性视频免费| 1024视频免费在线观看| 色综合婷婷激情| 搡老岳熟女国产| 最新的欧美精品一区二区| 脱女人内裤的视频| 国产淫语在线视频| 亚洲五月婷婷丁香| 欧美性长视频在线观看| 十八禁网站免费在线| 亚洲国产av新网站| 又紧又爽又黄一区二区| 国产一区二区三区在线臀色熟女 | 欧美日韩一级在线毛片| 成年版毛片免费区| 在线天堂中文资源库| 日本vs欧美在线观看视频| 老司机亚洲免费影院| 久久毛片免费看一区二区三区| 香蕉国产在线看| 大片免费播放器 马上看| 精品一区二区三卡| 一进一出抽搐动态| 欧美 亚洲 国产 日韩一| 精品一品国产午夜福利视频| 国产黄频视频在线观看| 欧美激情久久久久久爽电影 | 日日夜夜操网爽| 少妇精品久久久久久久| 男女边摸边吃奶| 国精品久久久久久国模美| 免费黄频网站在线观看国产| 制服诱惑二区| 午夜激情av网站| 91精品三级在线观看| 性少妇av在线| 久久婷婷成人综合色麻豆| 肉色欧美久久久久久久蜜桃| 一级毛片电影观看| 国产野战对白在线观看| 国产黄色免费在线视频| 亚洲国产欧美网| 一区二区日韩欧美中文字幕| 最新在线观看一区二区三区| 国产精品98久久久久久宅男小说| 国产午夜精品久久久久久| 操美女的视频在线观看| 成人特级黄色片久久久久久久 | 老司机福利观看| 国产亚洲一区二区精品| 免费观看人在逋| 最近最新中文字幕大全电影3 | 天天躁夜夜躁狠狠躁躁| 久久亚洲真实| 亚洲成人手机| 久久天堂一区二区三区四区| 两人在一起打扑克的视频| 在线观看免费视频日本深夜| 狠狠婷婷综合久久久久久88av| 成年女人毛片免费观看观看9 | 99国产精品一区二区三区| 免费看a级黄色片| 狠狠狠狠99中文字幕| 国产区一区二久久| 建设人人有责人人尽责人人享有的| 国产男女超爽视频在线观看| 人人澡人人妻人| 午夜福利视频在线观看免费| 美女主播在线视频| 国产日韩欧美视频二区| 一夜夜www| 精品国产乱子伦一区二区三区| 无遮挡黄片免费观看| 十八禁人妻一区二区| 色综合婷婷激情| 亚洲视频免费观看视频| 老司机福利观看| 三上悠亚av全集在线观看| 国产成人精品无人区| 久久影院123| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩视频精品一区| 免费日韩欧美在线观看| 美女高潮到喷水免费观看| 免费av中文字幕在线| 天天躁夜夜躁狠狠躁躁| 大香蕉久久成人网| 亚洲一区中文字幕在线| 人人妻人人澡人人爽人人夜夜| 欧美一级毛片孕妇| 大型黄色视频在线免费观看| 最新的欧美精品一区二区| 窝窝影院91人妻| av欧美777| 欧美日韩视频精品一区| 国产精品一区二区精品视频观看| 国内毛片毛片毛片毛片毛片| 午夜福利视频精品| 国产精品久久久人人做人人爽| 两人在一起打扑克的视频| 水蜜桃什么品种好| 国产欧美亚洲国产| 最近最新中文字幕大全电影3 | 无遮挡黄片免费观看| 成人黄色视频免费在线看| 欧美日韩av久久| 天天添夜夜摸| 亚洲欧美日韩高清在线视频 | 日本av手机在线免费观看| 精品视频人人做人人爽| 操美女的视频在线观看| 狠狠精品人妻久久久久久综合| 久久香蕉激情| 热re99久久精品国产66热6| 日韩欧美免费精品| av国产精品久久久久影院| 精品一区二区三卡| 国产精品亚洲av一区麻豆| 老司机靠b影院| 国产精品 国内视频| 国产亚洲精品一区二区www | 亚洲一卡2卡3卡4卡5卡精品中文| 成人特级黄色片久久久久久久 | 丝袜美足系列| 99re6热这里在线精品视频| av网站在线播放免费| 男人舔女人的私密视频| 美女主播在线视频| 蜜桃国产av成人99| 中文欧美无线码| 亚洲一区中文字幕在线| 国产伦人伦偷精品视频| 9191精品国产免费久久| 久久精品人人爽人人爽视色| 国产精品一区二区在线不卡| 免费少妇av软件| 精品国产乱码久久久久久小说| 在线观看免费视频日本深夜| 亚洲精品自拍成人| 精品福利观看| 51午夜福利影视在线观看| 在线观看免费高清a一片| 国产不卡一卡二| 丝袜人妻中文字幕| 亚洲专区字幕在线| 国产一区二区三区在线臀色熟女 | 国产黄频视频在线观看| 欧美乱码精品一区二区三区| 一个人免费在线观看的高清视频| 欧美日韩亚洲国产一区二区在线观看 | 成在线人永久免费视频| 亚洲成人手机| 久久久国产精品麻豆| 黄色a级毛片大全视频| 日韩欧美一区视频在线观看| 成年女人毛片免费观看观看9 | 日本欧美视频一区| 最新在线观看一区二区三区| 国产日韩欧美亚洲二区| 日本vs欧美在线观看视频| 久久久久久久大尺度免费视频| 一边摸一边抽搐一进一小说 | 啦啦啦免费观看视频1| 日本撒尿小便嘘嘘汇集6| 美女视频免费永久观看网站| tocl精华| 热re99久久精品国产66热6| 午夜福利在线观看吧| 岛国在线观看网站| 男女之事视频高清在线观看| 中文字幕精品免费在线观看视频| 欧美激情高清一区二区三区| 午夜91福利影院| 女性生殖器流出的白浆| 中文字幕av电影在线播放| 欧美黑人欧美精品刺激| svipshipincom国产片| 天天影视国产精品| 少妇的丰满在线观看| 亚洲少妇的诱惑av| 夜夜夜夜夜久久久久| 伊人久久大香线蕉亚洲五| tocl精华| 12—13女人毛片做爰片一| 国产成人精品无人区| 亚洲 欧美一区二区三区| 自线自在国产av| 人成视频在线观看免费观看| 国产精品久久久久久精品电影小说| 日韩欧美免费精品| 国产精品久久久久成人av| 免费高清在线观看日韩| 欧美乱码精品一区二区三区| 99国产综合亚洲精品| 免费一级毛片在线播放高清视频 | 精品卡一卡二卡四卡免费| 国产精品九九99| 亚洲欧美日韩高清在线视频 | 亚洲av片天天在线观看| 天天添夜夜摸| 高清欧美精品videossex| 露出奶头的视频| 国产免费视频播放在线视频| 五月开心婷婷网| 国精品久久久久久国模美| 国产欧美日韩精品亚洲av| 久久影院123| 国产不卡av网站在线观看| 黄色丝袜av网址大全| av视频免费观看在线观看| 老司机福利观看| 国产xxxxx性猛交| 国产成人欧美在线观看 | 国产精品av久久久久免费| 亚洲av成人一区二区三| 99精国产麻豆久久婷婷| 我要看黄色一级片免费的| 91av网站免费观看| 欧美人与性动交α欧美软件| 国产伦人伦偷精品视频| 一边摸一边抽搐一进一出视频| 国产成人精品久久二区二区91| 757午夜福利合集在线观看| 肉色欧美久久久久久久蜜桃| 国产老妇伦熟女老妇高清| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲精品国产色婷小说| 欧美激情高清一区二区三区| 在线 av 中文字幕| 一级,二级,三级黄色视频| 一二三四在线观看免费中文在| 欧美在线一区亚洲| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久人妻精品电影 | 美女高潮喷水抽搐中文字幕| 色视频在线一区二区三区| 国产精品二区激情视频| 精品少妇久久久久久888优播| 国产欧美日韩精品亚洲av| 99精品在免费线老司机午夜| 国产成人精品在线电影| 久久精品亚洲精品国产色婷小说| 高清黄色对白视频在线免费看| 老熟妇乱子伦视频在线观看| 色精品久久人妻99蜜桃| 亚洲精品在线美女| 国产极品粉嫩免费观看在线| 国产日韩欧美亚洲二区| 午夜激情久久久久久久| 一区二区日韩欧美中文字幕| 最近最新中文字幕大全电影3 | 久久久水蜜桃国产精品网| xxxhd国产人妻xxx| 成人亚洲精品一区在线观看| 高清欧美精品videossex| 亚洲精品中文字幕一二三四区 | 大香蕉久久成人网| 欧美日韩亚洲高清精品| 无人区码免费观看不卡 | www.精华液| 久久毛片免费看一区二区三区| 免费av中文字幕在线| 精品一品国产午夜福利视频| 十八禁高潮呻吟视频| 久久久精品国产亚洲av高清涩受| 无遮挡黄片免费观看| 欧美日韩福利视频一区二区| av线在线观看网站| 久久久久国产一级毛片高清牌| av在线播放免费不卡| 久久久国产一区二区| 80岁老熟妇乱子伦牲交| 黄色a级毛片大全视频| 欧美激情极品国产一区二区三区| 亚洲欧美一区二区三区久久| 69av精品久久久久久 | 国产黄频视频在线观看| 午夜精品国产一区二区电影| 久久久久久亚洲精品国产蜜桃av| 日韩精品免费视频一区二区三区| 国产日韩一区二区三区精品不卡| 欧美激情极品国产一区二区三区| 亚洲成人手机| 黄频高清免费视频| 啦啦啦中文免费视频观看日本| 久久久欧美国产精品| 中文字幕最新亚洲高清| 热99久久久久精品小说推荐| 我的亚洲天堂| 国产精品自产拍在线观看55亚洲 | 99久久人妻综合| 日韩大片免费观看网站| 精品人妻在线不人妻| 黄片播放在线免费| 精品久久久久久久毛片微露脸| 亚洲va日本ⅴa欧美va伊人久久| 精品国产乱码久久久久久男人| 欧美精品啪啪一区二区三区| 天天操日日干夜夜撸| 欧美激情久久久久久爽电影 | 视频区图区小说| 成人特级黄色片久久久久久久 | 高清在线国产一区| 日本欧美视频一区| 亚洲免费av在线视频| 亚洲全国av大片| 成人黄色视频免费在线看| 成人特级黄色片久久久久久久 | 久久ye,这里只有精品| 国产一区有黄有色的免费视频| 精品久久久久久电影网| 另类精品久久| 中国美女看黄片| 水蜜桃什么品种好| 亚洲一码二码三码区别大吗| 一区二区三区乱码不卡18| 国产精品一区二区在线观看99| 日日夜夜操网爽| 欧美日韩视频精品一区| 日韩中文字幕欧美一区二区| 精品国产一区二区久久| 国产日韩欧美在线精品| 国产亚洲午夜精品一区二区久久| 99re6热这里在线精品视频| 99久久99久久久精品蜜桃| 纵有疾风起免费观看全集完整版| 麻豆乱淫一区二区| 精品国产乱码久久久久久小说| tube8黄色片| 欧美精品av麻豆av| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品中文字幕在线视频| e午夜精品久久久久久久| 国产主播在线观看一区二区| 精品国内亚洲2022精品成人 | av超薄肉色丝袜交足视频| 美女视频免费永久观看网站| 一本色道久久久久久精品综合| 人妻 亚洲 视频| 日本av手机在线免费观看| 极品教师在线免费播放| 一级黄色大片毛片| 三级毛片av免费| avwww免费| 精品国内亚洲2022精品成人 | 亚洲精品自拍成人| 男人舔女人的私密视频| 欧美另类亚洲清纯唯美| 日本黄色日本黄色录像| 我要看黄色一级片免费的| 欧美亚洲 丝袜 人妻 在线| 丝瓜视频免费看黄片| 欧美一级毛片孕妇| 亚洲久久久国产精品| 日韩免费高清中文字幕av| 18在线观看网站| 天天躁夜夜躁狠狠躁躁| 国产精品久久电影中文字幕 | 五月开心婷婷网| 欧美激情久久久久久爽电影 | 亚洲国产成人一精品久久久| 大型av网站在线播放| 成人永久免费在线观看视频 | 国产免费视频播放在线视频| 日韩熟女老妇一区二区性免费视频| 国产精品一区二区免费欧美| 一二三四社区在线视频社区8| 淫妇啪啪啪对白视频| 美女高潮喷水抽搐中文字幕| 久久人人爽av亚洲精品天堂| 亚洲久久久国产精品| av网站在线播放免费| 国产欧美日韩一区二区三| 高潮久久久久久久久久久不卡| 久久久国产精品麻豆| 久久精品亚洲精品国产色婷小说| 91成年电影在线观看| 精品卡一卡二卡四卡免费| 免费观看av网站的网址| 色老头精品视频在线观看| 中亚洲国语对白在线视频| 建设人人有责人人尽责人人享有的| 久久九九热精品免费| 男男h啪啪无遮挡| 国产伦理片在线播放av一区| 我的亚洲天堂| 老司机影院毛片| 在线 av 中文字幕| 亚洲精品自拍成人| 老司机靠b影院| 菩萨蛮人人尽说江南好唐韦庄| 波多野结衣av一区二区av| 婷婷丁香在线五月| 国产高清激情床上av| 国产aⅴ精品一区二区三区波| 精品少妇黑人巨大在线播放| 精品一区二区三区av网在线观看 | 久久久久久亚洲精品国产蜜桃av| 精品一区二区三卡| 欧美另类亚洲清纯唯美| 亚洲第一av免费看| 自拍欧美九色日韩亚洲蝌蚪91| 日韩 欧美 亚洲 中文字幕| 欧美国产精品一级二级三级| 国产精品免费一区二区三区在线 | 午夜激情久久久久久久| 亚洲成人国产一区在线观看| 91精品国产国语对白视频| 国产男女超爽视频在线观看| 成人永久免费在线观看视频 | 我要看黄色一级片免费的| 少妇裸体淫交视频免费看高清 | 欧美国产精品va在线观看不卡| 一夜夜www| 97人妻天天添夜夜摸| 午夜福利在线观看吧| 国产一区二区三区视频了| 热99re8久久精品国产| tube8黄色片| 多毛熟女@视频| 国产精品二区激情视频| 久久国产亚洲av麻豆专区| 91av网站免费观看| 国产精品一区二区精品视频观看| 亚洲五月婷婷丁香| 午夜两性在线视频| 两性夫妻黄色片| 久久精品国产a三级三级三级| 久久久久视频综合| 欧美变态另类bdsm刘玥| av电影中文网址| 69精品国产乱码久久久| 80岁老熟妇乱子伦牲交| 美国免费a级毛片| 日本av免费视频播放| 久久亚洲精品不卡| 午夜日韩欧美国产| 久久九九热精品免费| 精品一品国产午夜福利视频|