Integrating Spectral and Textural Information for Monitoring the Growth of Pear Trees Using Optical Images from the UAV Platform

: With the recent developments of unmanned aerial vehicle (UAV) remote sensing, it is possible to monitor the growth condition of trees with the high temporal and spatial resolutions of data. In this study, the daily high-throughput RGB images of pear trees were captured from a UAV platform. A new index was generated by integrating the spectral and textural information using the improved adaptive feature weighting method (IAFWM). The inter-relationships of the air climatic variables and the soil’s physical properties (temperature, humidity and conductivity) were ﬁrstly assessed using principal component analysis (PCA). The climatic variables were selected to independently build a linear regression model with the new index when the cumulative variance explained reached 99.53%. The coefﬁcient of determination (R 2 ) of humidity (R 2 = 0.120, p = 0.205) using linear regression analysis was the dominating inﬂuencing factor for the growth of the pear trees, among the air climatic variables tested. The humidity (%) in 40 cm depth of soil (R 2 = 0.642, p < 0.001) using a linear regression coefﬁcient was the largest among climatic variables in the soil. The impact of climatic variables on the soil was commonly greater than those in the air, and the R 2 grew larger with the increasing depth of soil. The effects of the ﬂuctuation of the soil-climatic variables on the pear trees’ growth could be detected using the sliding window method (SWM), and the maximum absolute value of coefﬁcients with the corresponding day of year (DOY) of air temperature, soil temperature, soil humidity, and soil conductivity were conﬁrmed as 221, 227, 228, and 226 (DOY), respectively. Thus, the impact of the ﬂuctuation of climatic variables on the growth of pear trees can last 14, 8, 7, and 9 days, respectively. Therefore, it is highly recommended that the adoption of the integrated new index to explore the long-time impact of climate on pears growth be undertaken.


Introduction
Climate change has induced variations in temperature, precipitation, and solar radiation and aggregated variables, such as chilling and degree-day units with high impact on the vegetation growth and phenology shifts, and changes of growing degree days [1][2][3][4]. Climate change may significantly affect the phenology of vegetation, although there remains great uncertainty in spatial patterns and different species. Vegetation is very sensitive to these climatic changes and even an increase of 1 • C may lead to the advance and the delay of several days for leaf on and leaf off, respectively [5][6][7]. Thus, investigating the impact of climatic variables on the growth of vegetation is not only to investigate to what extent has the climate affected the growth of vegetation, but also to understand how the vegetation responds to these changes and to make effective adaptive measures. China contains the most plentiful Pyrus germplasm resources in the world, being the world's leading country in pear (Pyrus) production, accounting for 65.7% of world production and 72.3% of the world's growing area. Pears rank third among fruit industries in China, after oranges and apples. Therefore, timely monitoring of the growth condition of pear trees is of vital importance [8]. Since good growth of the pear trees is very important to guarantee agricultural production, thus, timely monitoring of the growth condition of the pear trees is of vital importance.
Currently, the field-scale destructive sampling method (FDSM) is a very precise approach for investigating the growth condition of pear trees. However, it is also very labor-intensive and time-consuming [9]. The Earth Observation Satellite (EOS) providing multi-spectral images can be applied to monitor the growth status of vegetations. However, EOS is limited by clouds, the specific revisiting periods, and the coarse spatial and spectral resolutions. In recent years, there have been important developments in remote sensing technology and sensors in terms of mass reduction, cost and operability (e.g., preprogramed flights) with simplicity of operation, which make use of unmanned aerial vehicles (UAV) or drones viable for large-scale agricultural applications [10][11][12]. Nowadays, UAVs present an affordable means to obtain images with very high spatial (centimeter-level) resolution, as well as a high temporal resolution acquisition adequate to monitor the various phases of the crop development throughout the growing season [13,14]. This programable, high spatial and temporal data are valuable data for high throughput plant phenotyping, which is the description of the observable properties or behavior of plants [15]. Therefore, drones equipped with common RGB cameras provide a great opportunity for advanced crop characterization since these optical sensors are comparatively inexpensive [16,17]. However, it is clear that there is still a considerable gap between the data and information provided by RGB cameras and the great potential agricultural-phenotyping offered by UAVs data is still yet to be fully realized [18]. The high-throughput images from UAV contain plenty of spectral and textural information that can be applied for image classification as well as for the detection of growth change and damage, among others, which is paramount to support efficient and sustainable agronomic practices in line with precision agriculture. In remote sensing domains, PCA is commonly used to transfer the multilayers of remote sensing data into fewer layers with almost the same useful information for knowledge extraction [19][20][21]. PCA can also be applied to determine the weight of multiple factors and to rank the importance of each influencing factor [22].
The spectral information can be extracted using the spectral vegetation indices (VIs) and the textural properties (TPs) can be calculated from the gray level co-occurrence matrix (GLCM). The VIs were commonly adopted for their convenience of usage while the TPs have only raised limited interest for vegetation monitoring. Zhang et al. calculated some commonly applied VIs, based on the UAV RGB images and the vegetation information was extracted with the iterative method based on optimal threshold [23]. Li et al. assessed the crop lodging using the object-based classification, and the feature selection of spectral information was extracted using various VIs calculated from RGB images [24]. Bendig et al. integrated the selected VIs and plant height information by developing the allometric model to estimate biomass in a summer barley experiment [25]. Lu et al. combined the coverage of the canopy, the plant height, and VIs extracted from RGB images captured from the UAV platform to predict the nitrogen concentration in the leaves of summer maize [26]. The TPs are an important implementation for the VIs, and the TPs contain relevant information for the improvement of the VIs. Then TPs are the dominating indicator for image classification and object detection especially when the VIs are saturated [27][28][29]. Therefore, there is a need to generate improved indices for integrating the spectral and textural information in remote sensing domains.
In this study, the daily high-throughput RGB images from the UAV platform were captured concerning the pear trees in an orchard. The daily sequence of RGB images were processed and analyzed with the climatic variables. This work focuses on addressing the following questions: (1) how to combine spectral and textural information to monitor the growth of pear trees; (2) how to monitor the growth of pear trees using the new index; and (3) how to explore the effects of fluctuation in climatic variables on the pear trees' growth.

Study Area
The study area is located in Nanpi County, Cangzhou city, Hebei province, North China Plain (NCP), China (38.00 • N, 116.40 • E) ( Figure 1). The NCP is one of China's major food production regions with irrigated agriculture [30,31]. Meanwhile, fruit and vegetables also accounted for a large proportion of the agricultural production in this region. This study was carried out in a pear orchard of 28.81 ha that included three pear cultivars (Xinliqi, Yuluxiang, and Huangjinyoupan) ( Figure 1). The trees were planted during 2018 and were well managed by the Agricultural and Rural Bureau (ARB) of Nanpi County. The study area is a typical orchard located on the NCP, and the mean annual temperature of this region is 12 • C with the average annual precipitation being 500 mm (http://he.cma.gov.cn/cz/, accessed on 5 March 2021). The soil of this pear orchard is representative of the NCP, containing three layers: A-AB-BK. The A, AB, and BK layers represent the depth from 0 to 20, 21 to 48, and 49 to 100 cm of the soil and these data were obtained from the National Earth System Science Data Center (NESSDC), (http://soil.geodata.cn/, accessed on 5 March 2021).

Data Collection and Preprocessing
The DJI Phantom 4 Pro V 2.0 was selected as the data collection platform for acquiring the RGB images to monitor the growth condition of pear trees in the orchard. The resolution of the image was 5472 × 3648 pixels with the focal lens being 8.8 mm. The flight altitude was set to 75 meters and the flight speed was strictly controlled as 7.5 m/s within the commercial software Altizure (V4.7.0.196, https://www.altizure.com, accessed on 5 March 2021, Shenzhen, China). Furthermore, the forward and side overlap of flights with each being 80% and 70% to ensure the high spatial resolution. The daily RGB images of UAV were captured and collected from 23rd August (235, day of year, DOY) to 6th September (249, DOY) including 15 different dates in the year 2020. The single images were processed within Pix4D Mapper (Lausanne, Switzerland) under a standard procedure, and the spatial resolution of the images was 2.1 cm. The tagged image format files (TIF) were generated and saved for each date and the image was clipped using the same regions of interest (ROI) within ENVI 5.3 (https://www.l3harrisgeospatial.com, accessed on 5 March 2021, Broomfield, CO, USA). Thus, the 15 mosaiced images of the same region were obtained and applied for further analysis.

Soil-Climatic Data
The daily temperature and humidity in the air and the soil were recorded and obtained from the national weather station of Nanpi (code 54719) that was managed by the China Meteorological Data Service Center (CMDSC, http://data.cma.cn/en, accessed on 5 March 2021). The hourly climatic data containing the temperature ( • C) in the air (average temperature, maximum temperature, and minimum temperature), and soil data at 20, 40, 60 cm depth of temperature ( • C), humidity (%), conductivity (us/cm) were obtained from this weather station. These climatic data were processed under a standard procedure and the daily climatic variables were obtained for further analysis.

Vegetation Index and Grey Level Co-Occurrence Matrix
The spectral and textural information were extracted using the VI and the properties of GLCM through the RGB images, respectively. To make the results comparable in quantitative remote sensing, the original RGB bands were each normalized from 0 to 1 by dividing the sum of original red, green and blue bands. The images of the studied orchard represent an anisotropic cover (or discontinuous) with complex mix reflectance from both the row and inter-row of pear trees (grass and background of soil) ( Figure 1). Thus, to mitigate the strong reflectance differences, the pixels in the images were classified into green pixels (rows of pear trees) and non-green pixels (inter-row grass and background of soil) using the difference of vegetation index as referred to in [32]. To be more specific, the difference of excess green (EGVI = 2G − R − B) and excess red (ERVI = 1.4 × R − B) were calculated for the image of each date [33]. The pixels of positive value were defined as green pixels (pear trees) and those of negative value were defined as nongreen pixels (soil), respectively. Thus, the binary images were made based on the normalized RGB image using this method for each date. Then, the commonly applied VIs were calculated based only on the green pixels of images for each date using the formulations presented in Table 1.

Vegetation Indices Formulation Reference
Newly built R, G, and B represent the normalized red (R), green (G), and blue (B) bands, respectively.
The texture is one of the most important characteristics for knowledge extraction in digital image processing domains [38][39][40]. The gray level co-occurrence matrix (GLCM) functions, which characterize the texture of an image by calculating how often pairs of pixels with specific values and in a specified spatial relationship occur in an image, was used to create a GLCM, to extract the contrast, correlation, energy, and homogeneity of images [24]. The normalized RGB image of each subsample image was converted into the grey-scale levels of the image. The GLCM of each date was calculated using these grey-scale levels of images. The textural information of an image can be derived using the properties of contrast, correlation, energy, and homogeneity extracted from GLCM, using the formulations presented in Table 2. Table 2. Formulations of the commonly used textural features based on the gray level co-occurrence matrix (GLCM).

Textural Properties Formulation
where, i and j are the number of rows and columns of the image. The p(i,j) is the normalized co-occurrence matrix.
The σi and σj are the standard deviation of row i and column j, respectively.
To make sure the reliable results of the analysis, the VI and GLCM were extracted and shown using all three cultivars together, and independently, a single cultivar. The temporal dynamic changes of VI were assessed and compared with each other. If there was little difference in the VI, and then a further analysis of the impacts of the climatic variables on the growth of pear trees was conducted, based on VI and GLCM calculated from the whole image containing all three cultivars of pear trees. If there was much difference in VI and GLCM, further analysis was independently conducted for the three cultivars of pear trees.

The New Index Based on Spectral and Textural Data
The improved adaptive feature weighting method (IAFWM) automatically measures the distance and adjusts the weights of spectral and textural features [41][42][43]. The weighting factors of spectral and textural information were confirmed using the IAFWM, and thus, the new index was generated, based on the updated weighting factors. The whole process can be divided into four sections: (1) Connecting the variables of six VIs (spectral information) and four features (textural information) of GLCM in series.
(2) Calculating the mean and variance values of each variable, and measuring the standardized distance of every two variables using the following equation: where, d norm is the standardized distance, ϕ i and ϕ j are the mean values, and ω i and ω j are the standard deviation of each variable, respectively. (3) Calculating the sum of the standardized distance between every two features and the weighting factors using the total of spectral and texture features in proportion. The detailed procedure is as follows: Assuming the feature expressions of all variables x (2)n , and m and n are the numbers of spectral and textural information, respectively. If w 1 and w 2 are the weighting factors of spectral information and textural information, then the following formulations can be obtained: where, the d ijk is the standardized distance of the kth feature between i and j. The new index by integrating features in series can be expressed using updated Equation (5).
The new index was generated and normalized (0 to 1), and the sequence was adopted for further analysis.

Impact of Climatic Variables on the Pear Trees' Growth
Before assessing the impacts of climatic variables on the growth condition of pear trees, the PCA was applied to identify the main influencing climatic factors, both in soil and in the air. The detailed procedures included the following steps: (1) removing the average value separately for each sample; (2) calculating the covariance matrix; (3) calculating the eigenvalues and eigenvectors of the covariance matrix; (4) sorting the eigenvectors in descending order with respect to their eigenvalues to form the transformation matrix; (5) obtaining the principal components through the transformation matrix. In this study, the PCA assessed the inter-relationships of the climatic variables and the dimension whose cumulative variance was greater than 99.53% was selected as the feature after PCA dimensionality reduction.
The confirmed main influencing climatic variables using PCA were independently applied to build a linear relationship with the new index generated using IAFWM. The coefficients of linear regression analysis were separately obtained and compared with each other. The maximum value of the coefficient of determination (R 2 ) with the corresponding soil-climatic variable was confirmed as the main (dominating) influencing factor. The dominating soil-climatic variables that influence the growth condition of pear trees can be confirmed using this method. Therefore, the R 2 between the climatic variables in the air, climatic variables in the soil, and the new index could be obtained and compared with each other.
For assessing the effects of fluctuation of climatic variables on the growth condition of pear trees, the first date of soil-climatic variable observation varies from 214 to 235 DOY with the length each being 15 days (235 to 249 DOY). Thus, the data of the soil-climatic variable was a sequence that was the same length as the new index. Then, the linear regression analysis was built between the independent changing soil-climatic variables and the new index. Through this process, the sliding window (SW) was applied, and a one-day interval was set as the sliding step. The temporal changes of coefficients using linear regression analysis between the soil-climatic variables and the new index could thus be obtained. There were, in total, 22 kinds of coefficients as the first date of climatic variable varied from 214 to 235 DOY with each step being one day to confirm the most optimal DOY that was most correlated with the growth condition of pear trees (new index). The maximum absolute value of the linear regression coefficient was confirmed for each soil-climatic variable and the corresponding DOY was confirmed. Therefore, the effects of fluctuation of soil-climatic variables on the growth of pear trees can be detected using the maximum absolute value and the corresponding DOY. Since there were four types of soil-climatic data (air temperature, soil temperature, soil humidity, and soil conductivity), thus, the averages of the coefficients of the different types were calculated and compared with each other.

Changes in Spectral and Textural Information during the Monitoring Period
Single images acquired from the UAV platform concerning the pear trees from 23rd August (235, DOY) to 6th September (249, DOY) were obtained and processed using the method introduced in Section 2.   The original bands were first normalized within 0 to 1 using the method introduced in Section 2.3 Method. The spectral information was extracted using the equations introduced in Table 1, based on these normalized images. The VIs calculated, based on the whole image including three cultivars of pear trees, were shown in a different color for each date during the period from 235 to 249 DOY (Figure 4). It can be observed that the trend of all VIs using the whole image including three cultivars was similar, of which the values of the red green blue vegetation index (RGBVI) and the gray were the highest and lowest, respectively. In order to eliminate the uncertaintities that arise from the different cultivars of pear trees, the VI was calculated, based on the independent three cultivars of pear trees and is shown in Appendix A Figure A1. Compared with the dynamic changes of VI in Figures 3 and A1, it can be noted that there was little difference in the VI calculated, based on the three cultivars of pear trees. Thus, further analysis was based on the VI calculated from the whole image with all three cultivars of pear trees included. The textural information was extracted based on the GLCM calculated, based on the image of each date using the Equations presented in Table 2. The four daily properties containing contrast, correlation, energy, and homogeneity from the whole image were extracted using the GLCM for the period from 235 to 249 DOY ( Figure 5). It can be seen that the correlation and the energy dramatically varied during the observation period, while the contrast and the homogeneity changed little. The four properties containing contrast, correlation, energy, and homogeneity were also independently calculated for the three cultivars of pear trees.

The Linear Regression between Soil-Climatic Variables and the New Index
Using the improved adaptive feature weighting method (IAFWM), the weighting factors using the total of spectral and texture features in proportion were obtained. The ratio between spectral and texture features was set as 5.616. Based on the ratio, the new index was generated by multiplying all spectral information by 5.616 and adding the original textural information. Then the new index was normalized by dividing the largest value in the new index ( Figure 6). The climatic variables, including the temperature ( • C) in the air and temperature ( • C), humidity (%), and conductivity (us/cm) at 20, 40, 60 cm depth in the soil, were analyzed using the PCA method introduced in Section 2.3.3. The variance explained with the corresponding numbers of climatic variables using PCA are shown in Figure 7. There were 4 variables in the air and 9 variables in the soil, respectively. It can be noted that when the number of variables reached 7, the variance was explained by 99.21% of the total components. The ranking was as follows: soil humidity in 60 cm, soil conductivity in 40 cm, soil humidity in 40 cm, soil conductivity in 20 cm, soil conductivity in 60 cm, soil temperature in 60 cm, and soil temperature in 40 cm, since the climatic variables of the different depths in the soil all contained three properties: humidity, conductivity, and temperature. From the variance explained, when the number of variables reached 11 (9 variables in the soil and minimum temperature and humidity in the air), the variance explained 99.53% compared with the original dataset. Since the climatic variables in the air were meant to be investigated, thus only the two variables of minimum temperature and humidity in the air were added. Therefore, further analysis was based on the minimum temperature and humidity in the air, and all the variables in the soil. The linear regression analysis was separately applied between the selected climatic variables in the air and the new index. The coefficients of the climatic variables of minimum temperature and humidity in the air were 0.046 and 0.120, respectively. The linear regression analysis was also separately applied between the climatic variable in the soil and the new index. The R 2 between the climatic variables of temperature and soil humidity, and conductivity in the soil (20, 40, and 60 cm) and the new index were calculated and shown ( Figure 8). It can be observed that the coefficients of temperature and humidity became larger with the increase in the depth of soil. However, the R 2 of conductivity increased from 20 to 40 cm depth of soil, and then decreased to 60 cm depth of soil. It also can be concluded that humidity and conductivity were the dominating influencing factors for the growth of pear trees. Besides, the coefficients of climatic variables in soil were commonly larger than those of climatic variables in the air. Thus, the climatic variables in the soil may be more important than those in the air.

The Effects of Fluctuation in Climatic Variables on the Growth of Pear Trees
The fluctuation of climatic variables on the growth of pear trees was evaluated using the linear regression method and the SW method with a step of one day, as introduced in Section 2.3.2. The temporal changes of coefficients using linear regression analysis were calculated for climatic variables in the air and in the soil. To better show the different changes in impact, the climatic variables were divided into four classes containing air temperature, soil temperature, soil humidity, and soil conductivity. The average of the coefficients for the four classes were calculated and are shown in Figure 9. The maximum absolute value of coefficients for air temperature, soil temperature, soil humidity, and soil conductivity occurred on 221, 227, 228, 226 (DOY), respectively. Since the first date of the UAV collection of images happened on 235 DOY and thus, the new index was most correlated with the air temperature, soil temperature, soil humidity, and soil conductivity 14, 8, 7, and 9 days before the data collection of the UAV image. The impact of fluctuations in climatic variables on the growth of pear trees was assessed, and it can be concluded that the impacts of air temperature, soil temperature, soil humidity, and soil conductivity can last 14, 8, 7, and 9 days, respectively.

Discussion
A UAV platform mounted with a high spatial resolution optical cameras can be used to collect images concerning fruit trees at high frequency. The high-throughput images contain plenty of information on the spectrum and texture [11,44]. The spectral and textural information can be applied in ecological and agricultural applications such as the identification of the phenological events in agriculture and forest, monitoring the growth condition of crops and trees, and predicting agricultural production [45][46][47][48][49][50][51][52]. Commonly, the spectral information can be extracted using various VIs composed by linear and nonlinear equations, and here the texture was extracted using the GLCM. In this study, the temporal changes of spectral information using various VIs were applied to identify the dynamic changes in the growth of the pear trees. The trends of the spectral information were similar and they have shown a strong inter-correlation. It was reported that the textural information was a very important measure for the VI when the spectral information was being saturated [53][54][55]. However, textural information was commonly ignored in agricultural and ecological applications. The temporal changes of textural information including the contrast, correlation, energy, and homogeneity were investigated. Correlation and energy changed dramatically during the observations while contrast and the homogeneity changed little. This may be because the growth condition of the pear trees was closely correlated with the types and components of the pigments in the leaves [56,57]. The changes in correlation and energy of the textural information in the images were reflections of the variations of pigments in the vegetation, including the pear trees [58][59][60]. Similar results were obtained from previous studies. The vegetation indices and textures were combined for estimating the leaf area index (LAI) of rice based on a UAV RGB image leaf area index (LAI) [61]. The imagery was obtained from UAV fixed-wing concerning the whole growth stages of wheat and the texture, color, and vegetation indices were integrated to monitor the growth conditions [62]. Similar research has proven that the integrated vegetation indices and the textures can better capture the dynamic changes of vegetation such as trees and crops. Thus, textural information should be included in the monitoring of fruit trees and related agricultural and ecological applications. Besides, the new index may have a potential usage for monitoring the growth of maize, wheat, and rice et al.
The R 2 between the new index and the soil temperature, soil humidity was relatively high and the R 2 grew with the depth of the soil (from 20, 40, to 60 cm). The humidity in the soil was most closely correlated with the new index, which meant that more humidity promoted the growth of the pear trees and higher greenness would be achieved. Meanwhile, proper nurturing of pear trees would lead to less bare soil and the avoidance of sun exposure, and thus, the soil might contain more humidity. Transpiration far exceeds evaporation when the plant roots are extensive in the soil. Conversely, less humidity will lead to the negative growth of pear trees, and thus, less greenness on the ground will increase sun exposure (and soil evaporation), and thus, the soil would have lower humidity. The data collection in this study was based on an RGB camera mounted on the UAV platform and it is commonly believed that advanced sensors such as multispectral and hyper-spectral sensors would have greater potential in assessing the vegetation-related applications. The normalized difference vegetation index (NDVI) can be calculated from multispectral images and solar-induced chlorophyll fluorescence (SIF) can be extracted from hyperspectral images [63,64]. The NDVI and SIF have greater capabilities to detect changes in vegetation than that of the knowledge extracted from only RGB images [65,66]. Therefore, it is recommended that UAV platforms mounted with advanced multispec-tral and hyperspectral sensors should be applied to collect the images concerning fruit trees [67]. The NDVI calculated from the red and near-infrared bands from multispectral and hyperspectral images might be more correlated with, and more sensitive to, the response of vegetation [68,69]. Besides, the data from light detection and ranging (LiDAR) can be applied to extract more information on texture, such as the length, width, and height [70,71]. Therefore, it is suggested that more advanced sensors such as multispectral cameras, high-spectral cameras, and LiDAR be applied to explore and investigate the growth conditions of fruit trees, but the cost−benefit ratio must also be assessed [72][73][74]. Besides, the fusion method of the multisource image should be explored, and thus, the different sources of data can be integrated and hidden knowledge may be exposed when the data fusion has been conducted. More detailed information concerning the vegetation may be derived from the combined data. Moreover, more advanced approaches such as machine learning and deep learning containing multilayers for data analysis should be explored as to their excellent performance in classification and regression models [75][76][77]. Machine learning and deep learning can be used to build multilayers using the original images and more precise information such as the growth conditions and agricultural yields which can be obtained from this [78,79].
With the development of new electronics, mini-computers and light-sensors, the use of UAVs mounted with light sensors such as UAV-based multispectral cameras, hyper-spectral cameras, and (LiDAR) has emerged as a promising alternative for ecological and agricultural applications. There have been plenty of related applications, such as monitoring the growth condition of forests, estimation of forest structural parameters, measurement of within-season tree height, change detection of grass and trees, and the agricultural production predictions [80][81][82][83]. UAVs can fly at relatively low altitudes acquiring high resolution images at high frequency for the users. Commonly, the flights were conducted using commercial software such as the Dji pro software (https://www.dji.com/cn/downloads/djiapp/dji-go-3, accessed on 5 March 2021, Shen Zhen, China) and these flights were strictly controlled and seldomly encountered yawing. However, the UAV is dangerous as there may be mechanical failure, signal interruption, and even collision with flying birds. The first UAV regulation was proposed in 1944, and the United Kingdom and Australia first published their UAV regulations in 2012. These days, with the development of 5 G and the higher precision of the Global Positioning System (GPS), new regulations must be developed to keep pace with the rapid emergence of UAVs [84].

Conclusions
The new index was generated using the improved adaptive feature weighting method by integrating the spectral and textural information derived from a sequence of highthroughput images. The monitoring of the growth condition of pear trees using the new index indicated that it can precisely extract the dynamic temporal changes of the pear trees. The new index may also have potential abilities in monitoring crops such as maize, wheat, and rice. The results of PCA showed that climatic variables in the soil were the main components, compared with the climatic variables in the air. The further effects assessed using linear regression analysis between the selected main components and the new index also showed that the impact of temperature and humidity in the soil was greater than that in the air. The impact from humidity and conductivity in the soil were the dominating influencing factors for the growth of pear trees. The effects of fluctuation in climatic variables on the growth of pear trees were assessed and evaluated, and it can be concluded that the impact of air temperature, soil temperature, soil humidity, and soil conductivity can last 14, 8, 7, and 9 days, respectively. It is strongly recommended that the new index with more advanced approaches such as machine learning and deep learning be applied for monitoring the growth condition of pear trees.

Conflicts of Interest:
The authors declare no conflict of interest.