Detection of Yunnan Pine Shoot Beetle Stress Using UAV-Based Thermal Imagery and LiDAR

: Infestations of Tomicus spp. have caused the deaths of millions of Yunnan pine forests in Southwest China; consequently, accurate monitoring methods are required to assess the damage caused by these pest insects at an early stage. Considering the limited sensitivity of optical reﬂectance on the early stage of beetle stress, the potential of thermal infrared (TIR) can be exploited for monitoring forest health on the basis of the change of canopy surface temperature (CST). However, few studies have investigated the impact of the leaf area index (LAI) on the accuracy of TIR data-based SDR assessments. Therefore, the current study used unmanned airborne vehicle (UAV)-based TIR and light detection and ranging (LiDAR) data to assess the capacity of determining the potential for using TIR data for determining SDR under different LAI conditions. The feasibility of using TIR for monitoring SDRs at the tree level and plot scales were analyzed using the relationship between SDR and canopy temperature. Results revealed that: (1) prediction accuracy of SDR from CST is promising at high LAI values and decreases quickly with LAI, and is higher at the single tree scale (R 2 = 0.7890) than at the plot scale (R 2 = 0.5532); (2) at either single tree or plot scale, a signiﬁcant negative correlation can be found between CST and LAI ( − 0.9121 at tree scale and − 0.5902 at plot scale); (3) LAI affects the transmission paths of sunlight and sensor, which mainly disturbs the relationship between CST and SDR. This article evaluated the high possibility of using TIR data to monitor SDRs at both tree and plot levels and assessed the negative impact of a low LAI (<1) on the relationship between temperature and SDR. Accordingly, when measuring forest health using TIR data, additional data sources are required to eliminate the negative impact of low LAIs and to improve the monitoring accuracy.


Introduction
Forest pests, such as Tomicus spp., have seriously affected the economic benefits and ecological environment of coniferous forests [1,2].Moreover, forest pests also cause the issues of forest degradation and disturbed ecosystem services and have an impact on forest carbon sequestration and sustainable forest management [3,4].The appearance and spread of pests and diseases can be reflected quickly and objectively using data acquired from high-precision satellite technology [5].Therefore, it is becoming increasingly necessary to use spatial data at various scales to rapidly locate and assess affected areas and the degrees of pest stress in forests [6].Optical remote sensing, which uses spectral characteristics to Appl.Sci.2022, 12, 4372 2 of 14 monitor forest pests, is limited to red and gray attacks (late stage) [7][8][9].Accordingly, some researchers use the water content sensitivity of thermal infrared (TIR) remote sensing to monitor forest pests during the early stage of an attack [10,11].
The canopy surface temperature (CST) reflects the transfer and exchange of energy within a soil-plant-atmosphere continuum (SPAC).The changes in CST monitored by TIR imaging exhibit responses to environment stress, such as drought [12,13], disease [14], and pest infestation [15].With improvements in the method and accuracy of the technique, the TIR satellite images can provide surface temperatures to determine CSTs at the regional scale [1,16].However, CSTs obtained from satellite images are unsuitable for the precise analysis of small areas because of their coarse resolution compared with optical images.
The resolution of thermal data was limited to 1-10 m for manned airborne sensors or to centimeter-scale resolutions for infrared thermal images via a settled flux tower.The advent of unmanned aerial vehicles (UAVs) supports the application of high-resolution thermal images with low-weight thermal cameras at centimeter-scale resolutions [17][18][19].These images can describe in detail the transpiration from canopies according to the basic relationship between canopy temperatures and evaporation in plant ecophysiological research [20,21].High-resolution thermal images can help to quantify the link between CST and canopy health and clarify the response mechanism of canopy physiological parameters to pest stress [14].
Forest pests can disrupt water equilibrium in plants and decrease the water content of the canopy.Whether by defoliating needles or boring into shoots, bark, or wood, pests may cause temporary or permanent damage to physiological and biochemical parameters of forests [22].The main species of pests responsible for destroying 1.5 million hectares of Yunnan pine forest in Southwest China over the past 20 years is Tomicus spp.[23,24].The beetles of these species attack the shoots of the Yunnan pine from May and commence wood boring in November.Over several weeks, the affected needle shoots gradually change from green to yellow and, typically, the needles turn red within several months (red attack).During December, the host trees bored by beetles lose their needles (gray attack) (this occurs within 3 years of the initial attack) [25].The spectral changes that occur during the late stage are easily identified by optical remote sensing [26].However, it is difficult to recognize the spectral characteristics of a green attack (early stage) [27].Despite the proposal of various remote sensing methods for evaluating forest damage, the forecasting of beetle attacks remains a challenge, particularly during the green attack stage [22,28].
Thermal infrared imaging techniques can be used to monitor abnormal canopy temperatures due to the disruption to the water cycle of canopy during the green attack stage [29].However, the use of thermal remote sensing for monitoring forest pests is uncommon.One reason is the difficulty of interpreting the CST, which is heavily influenced by environmental factors and physiological vegetation parameters [13].In agricultural research, the underlying surface of crops is homogeneous, and the temperature of the crop canopy receives little disturbance, making it easy to obtain.
The complex structure and irregular shape of forestry canopies result in the heterogeneity of the canopy [30], which increases the uncertainty of canopy temperatures extracted from UAV-based thermal images.Light detection and ranging (LiDAR) data provide a method for obtaining detailed three-dimensional (3D) structural parameters of tree crowns with dense cloud points, which are used to obtain the spatial distribution of foliage [31][32][33].The intensity of laser signal can produce a reflection characteristic in the near-infrared region of foliage within a canopy.Therefore, LiDAR has significant potential for monitoring forest health by detecting the severity of damage from forest pest infestations [34].However, obtaining fine-grained measurements that reflect the health of an individual canopy is limited when using TIR or LiDAR only.The mixed radiation from trees and the background presents difficulties for accuracy accessing the canopy structures using thermal remote sensing techniques.The data obtained by LiDAR serves as an additional source of information for capturing the fine structures of tree crowns, such as the leaf area index (LAI), fractional cover, and canopy density [35,36].
Few studies have used TIR for monitoring forest infestations at the canopy scale to assess the damage caused by pine-shoot beetles in Yunnan pine forests.Accordingly, this article investigated the capacity of thermal imagery for monitoring the shoot damage ratio (SDR) of the crowns of Yunnan pine trees.The degree of damage to Yunnan pine crowns from attack by Tomicus spp. was assessed using the SDR [23,37].The study also presents and discusses the impact of the LAI on the CST-SDR relationship using TIR and LiDAR datasets at canopy and plot scales in addition to corresponding field measurements.

Study Sites
The study sites (25 • 14 -25 • 29 N, 100 • 48 -101 • 3 E) were located on Tianfeng Mountain in Yunnan Province, Southwestern China (Figure 1) at an elevation of between 1720 and 2570 m.In this area, continuous drought and infestations of Tomicus spp.have led to a high number of Yunnan pine deaths.The annual average temperature in this area is 14.  the background presents difficulties for accuracy accessing the canopy structures using thermal remote sensing techniques.The data obtained by LiDAR serves as an additional source of information for capturing the fine structures of tree crowns, such as the leaf area index (LAI), fractional cover, and canopy density [35,36].Few studies have used TIR for monitoring forest infestations at the canopy scale to assess the damage caused by pine-shoot beetles in Yunnan pine forests.Accordingly, this article investigated the capacity of thermal imagery for monitoring the shoot damage ratio (SDR) of the crowns of Yunnan pine trees.The degree of damage to Yunnan pine crowns from attack by Tomicus spp. was assessed using the SDR [23,37].The study also presents and discusses the impact of the LAI on the CST-SDR relationship using TIR and LiDAR datasets at canopy and plot scales in addition to corresponding field measurements.

Study Sites
The study sites (25°14′-25°29′ N, 100°48′-101°3′ E) were located on Tianfeng Mountain in Yunnan Province, Southwestern China (Figure 1) at an elevation of between 1720 and 2570 m.In this area, continuous drought and infestations of Tomicus spp.have led to a high number of Yunnan pine deaths.The annual average temperature in this area is 14.2 °C, and the annual precipitation is 783.7 mm.Field measurements of tree parameters were obtained in September 2018.As shown in Figure 1, the study areas contained plots of 50 × 50 m.The coordinates of each plot and trees were measured using a real-time kinematic (RTK) device (HI-TARGET A8 GNSS) with an accuracy of ±2.5 mm.Variables, such as tree height (H), diameter at breast height (DBH), crown diameter, and SDR, were measured on each tree (DBH > 4 cm).As an indicator of the degree of Yunnan pine infestation damage, the SDR was defined as the ratio of the number of damaged shoots to total shoots for each crown (a value between 0% and 100%).The plot SDR was the average SDR of the single trees within a plot.Field measurements of tree parameters were obtained in September 2018.As shown in Figure 1, the study areas contained plots of 50 × 50 m.The coordinates of each plot and trees were measured using a real-time kinematic (RTK) device (HI-TARGET A8 GNSS) with an accuracy of ±2.5 mm.Variables, such as tree height (H), diameter at breast height (DBH), crown diameter, and SDR, were measured on each tree (DBH > 4 cm).As an indicator of the degree of Yunnan pine infestation damage, the SDR was defined as the ratio of the number of damaged shoots to total shoots for each crown (a value between 0% and 100%).The plot SDR was the average SDR of the single trees within a plot.

/ 100%
Each plot was divided into 25 smaller plots of 10 m × 10 m.The LAI of each small plot was obtained using plant canopy analyzer with a 270-degree view restrictor (LAI-2200, Li-Cor, Lincoln, NE, USA), which was located 1.2 m above the ground (shrubs were excluded).To prevent interference from sunlight, the LAI measurement of each small plot was performed in the afternoon or under cloudy conditions.The summary statistics for each tree variable of two plots are shown in Table 1.

Thermal Imagery
The TIR data were acquired using a high-sensitivity thermal sensor (TAU2, FLIR, Wilsonville, OR, USA) on a fixed-wing UAV between 13:00 and 13:30 on 17 September 2018.The TIR sensor was a focal plane array with an uncooled vanadium oxide microbolometer and heat sensitivity below 50 mK.The flight took place at a height of 280 m, and the flight path was designed with overlaps of 60% (side) and 80% (forward).The pixel value of the TIR images obtained from the thermal sensor was the digital number (DN) of the land surface thermal radiation, which required correction by temperature.The produced resolution of the TIR data was 0.23 m.

The Correction of Thermal Imagery
The temperatures of four ground features (tile, polyvinyl chloride board, wood, and asphalt road) with known emissivity (Table 2), were measured by an infrared thermal imager (T-420, FLIR Systems, Inc., Wilsonville, OR, USA).The coordinates of the features were measured using an RTK global positioning satellite device.Then, the linear regression relationship between the temperatures and DN values of the four features in the TIR images was established; consequently, the gray TIR image was converted into brightness temperature (Figure 2).

The Temperature of Ground Features within Plots
According to the actual sample investigation, the ground features were divided into five types, namely, bare land (BL), falling dead trees (FDT), dead trees with red crowns (RDT), living trees (LT), and shrubbery (SHR).The shapes and locations of the five types were extracted from visible UAV-based images with 0.02 m resolution that were shot on the same LiDAR data flight.The average temperatures of features at the plots were extracted from TIR images.

The Average Temperatures of the Small Plots
The 50 small plots of 10 m × 10 m were subdivided from two larger plots to analyze the effect of the LAI on the relationship between temperature and the SDR.The plot temperature (Tp) was the production of mixed of thermal radiation from all the ground features in the plot.When calculating the average temperature of a sample plot, if many shrubs are present, the contribution of shrub temperature should be considered because of the lower temperature of shrubs compared with pine-tree crowns.Therefore, LiDAR data were used to reverse LAI of a plot with a large number of shrubs.

Individual Tree Segmentation from LiDAR
The UAV-based LiDAR data were obtained from the core area of the affected Yunnan pine forests (including the two plots) using a UAV-mounted system (LiAir 200, GreenVlley Inc., Beijing, China) with a 40-channel Pandar40 laser sensor.The average flight altitude was 70 m above ground level, with a pulse rate of 10 HZ and a maximum viewing angle of 10° (from channels 2 to 30).The georeferenced coordinate of the LiDAR data was generated in the UTM 50 N zone using Word Geodetic System data.The LiDAR cloud points were classified into ground and above-ground points using the improved progressive TIN densification algorithm in the LiDAR360 software package.
The canopy height model with a resolution of 0.5 m was derived from the classified point cloud using a point cloud segmentation (PCS) algorithm, and the accuracy assessment of canopu segmentation was evaluated using three parameters, which included precision (p), recall (r), and F-score (F) [30].The structural parameters of the tree canopy (e.g.,

The Temperature of Ground Features within Plots
According to the actual sample investigation, the ground features were divided into five types, namely, bare land (BL), falling dead trees (FDT), dead trees with red crowns (RDT), living trees (LT), and shrubbery (SHR).The shapes and locations of the five types were extracted from visible UAV-based images with 0.02 m resolution that were shot on the same LiDAR data flight.The average temperatures of features at the plots were extracted from TIR images.

The Average Temperatures of the Small Plots
The 50 small plots of 10 m × 10 m were subdivided from two larger plots to analyze the effect of the LAI on the relationship between temperature and the SDR.The plot temperature (T p ) was the production of mixed of thermal radiation from all the ground features in the plot.When calculating the average temperature of a sample plot, if many shrubs are present, the contribution of shrub temperature should be considered because of the lower temperature of shrubs compared with pine-tree crowns.Therefore, LiDAR data were used to reverse LAI of a plot with a large number of shrubs.

LiDAR Data 2.3.1. Individual Tree Segmentation from LiDAR
The UAV-based LiDAR data were obtained from the core area of the affected Yunnan pine forests (including the two plots) using a UAV-mounted system (LiAir 200, GreenVlley Inc., Beijing, China) with a 40-channel Pandar40 laser sensor.The average flight altitude was 70 m above ground level, with a pulse rate of 10 HZ and a maximum viewing angle of 10 • (from channels 2 to 30).The georeferenced coordinate of the LiDAR data was generated in the UTM 50 N zone using Word Geodetic System data.The LiDAR cloud points were classified into ground and above-ground points using the improved progressive TIN densification algorithm in the LiDAR360 software package.
The canopy height model with a resolution of 0.5 m was derived from the classified point cloud using a point cloud segmentation (PCS) algorithm, and the accuracy assessment of canopu segmentation was evaluated using three parameters, which included precision (p), recall (r), and F-score (F) [30].The structural parameters of the tree canopy (e.g., locations, heights, crown areas, and volumes) were extracted from the point clouds of individual trees (Figure 3).The TIR images and the LiDAR canopy areas were registered using visible images via ground control points, which were measured by the RTK device.p = TP/(TP + FP), (3) where TP is the number of detected trees in a plot; FP is the number of trees that are falsely detected; FN is the number of trees that are omitted; and p, r, and F refer to the detection accuracy, the detection rate, and the overall accuracy, respectively.
Appl.Sci.2022, 12, x FOR PEER REVIEW 6 of 14 locations, heights, crown areas, and volumes) were extracted from the point clouds of individual trees (Figure 3).The TIR images and the LiDAR canopy areas were registered using visible images via ground control points, which were measured by the RTK device.
where TP is the number of detected trees in a plot; FP is the number of trees that are falsely detected; FN is the number of trees that are omitted; and p, r, and F refer to the detection accuracy, the detection rate, and the overall accuracy, respectively.

The Calculated Leaf Area Indices of Individual Trees
To analyze the impact of the background, the LAI of individual trees was selected as the parameter for evaluating the relationship between SDR and CST, which was calculated using the point clouds from the LiDAR data.The point clouds of individual trees were extracted from the segmentation results using a circumcircle of the crown and divided into canopy and ground parts.Then, the porosity of the canopy was calculated on the basis of the principle that a specific ground point was the reflection point at which the laser penetrated the canopy and hit the ground.Finally, the LAI of an individual tree was calculated according to the exponential relationship between the probability of finding canopy P and LAI in the random Beer-Lambert law of discontinuous canopy proposed by Shabanov [38]: * ln(1 ) where cloudtree is the point clouds of canopy, cloudtotal is the total point cloud of the circumcircle, and k is the empirical coefficient that describes the relevance between the LAI and the porosity of the tree crown (typically, k = 0.5 for a conifer forest).

The Calculated Leaf Area Indices of Individual Trees
To analyze the impact of the background, the LAI of individual trees was selected as the parameter for evaluating the relationship between SDR and CST, which was calculated using the point clouds from the LiDAR data.The point clouds of individual trees were extracted from the segmentation results using a circumcircle of the crown and divided into canopy and ground parts.Then, the porosity of the canopy was calculated on the basis of the principle that a specific ground point was the reflection point at which the laser penetrated the canopy and hit the ground.Finally, the LAI of an individual tree was calculated according to the exponential relationship between the probability of finding canopy P and LAI in the random Beer-Lambert law of discontinuous canopy proposed by Shabanov [38]: where cloud tree is the point clouds of canopy, cloud total is the total point cloud of the circumcircle, and k is the empirical coefficient that describes the relevance between the LAI and the porosity of the tree crown (typically, k = 0.5 for a conifer forest).

Correlation Analysis Method
To detect the shoot beetle stress using TIR and LiDAR data, the correlations between SDR and temperature at tree and plot scales were analyzed using Origin software (Origin-Lab, Northampton, MA, USA).At the individual tree scale, the data of trees (15 samples) with different LAIs were selected from the plots to calculate their LAI from the LiDAR point cloud dataset and evaluate the impact of LAI on the correlation between SDR and temperature.After weakening the impact of LAI, a linear regression method with a 95% level of confidence was applied to build the correlation between SDR and temperature.
To further analyze the correlation between SDR and temperature on a large scale, a total data of 50 small plots were analyzed with the same method, which was the same as the tree scale.The coefficient of determination (R 2 ) between SDR and temperature was used to assess the accuracy of the regression model.The Pearson correlation coefficient (r) was used to indicate the relationship between LAI and temperature.

CST Characteristics of Features within Plots
The temperatures of the five ground features were obtained using TIR imagery (Figure 4).BL had the highest temperature, and the temperature of SHR with high transpiration was lower than LT.Despite a large number of dead needles in RDTs, their RDT temperature was similar to that of FDTs due to the cessation of transpiration by the dead needles in the canopy.Considering the reduced contribution of SHR temperature to plot temperature (T p ), SHR features should be eliminated when the plot average temperature is calculated to analyze the correlation with plot SDR.
To detect the shoot beetle stress using TIR and LiDAR data, the correlations between SDR and temperature at tree and plot scales were analyzed using Origin software (OriginLab, Northampton, MA, USA).At the individual tree scale, the data of trees (15 samples) with different LAIs were selected from the plots to calculate their LAI from the LiDAR point cloud dataset and evaluate the impact of LAI on the correlation between SDR and temperature.After weakening the impact of LAI, a linear regression method with a 95% level of confidence was applied to build the correlation between SDR and temperature.
To further analyze the correlation between SDR and temperature on a large scale, a total data of 50 small plots were analyzed with the same method, which was the same as the tree scale.The coefficient of determination (R 2 ) between SDR and temperature was used to assess the accuracy of the regression model.The Pearson correlation coefficient (r) was used to indicate the relationship between LAI and temperature.

CST Characteristics of Features within Plots
The temperatures of the five ground features were obtained using TIR imagery (Figure 4).BL had the highest temperature, and the temperature of SHR with high transpiration was lower than LT.Despite a large number of dead needles in RDTs, their RDT temperature was similar to that of FDTs due to the cessation of transpiration by the dead needles in the canopy.Considering the reduced contribution of SHR temperature to plot temperature (Tp), SHR features should be eliminated when the plot average temperature is calculated to analyze the correlation with plot SDR.

Crown Segmentation from LiDAR Data
The 0.23 m resolution of the TIR images made it difficult to confirm the edges of the canopy and to locate individual affected trees, particularly when canopies containing several trees were adjacent.Therefore, this study used LiDAR point clouds to distinguish nearby trees.

Crown Segmentation from LiDAR Data
The 0.23 m resolution of the TIR images made it difficult to confirm the edges of the canopy and to locate individual affected trees, particularly when canopies containing several trees were adjacent.Therefore, this study used LiDAR point clouds to distinguish nearby trees.
Segmenting the crown of a single plant is the first step in extracting CST to obtain the temperature characteristics.Unfortunately, due to the coarser resolution of the TIR band compared with the optical bands for the UAV-based platform, the degeneration of the canopy boundary made it difficult to accurately extract the CST (Figure 5a).Because of the mixed radiation from the canopy and surrounding objects, the degeneration of pixels at the canopy boundary occurred in the TIR image.To segment the canopy using TIR data, it was necessary to determine the temperature threshold to calculate the canopy boundary.The canopy boundary of individual trees can be extracted from high-resolution UAV-based optical images.However, the accurate separation of overlapping crowns is impossible from optical images without using the vertical information.
Segmenting the crown of a single plant is the first step in extracting CST to obtain the temperature characteristics.Unfortunately, due to the coarser resolution of the TIR band compared with the optical bands for the UAV-based platform, the degeneration of the canopy boundary made it difficult to accurately extract the CST (Figure 5a).Because of the mixed radiation from the canopy and surrounding objects, the degeneration of pixels at the canopy boundary occurred in the TIR image.To segment the canopy using TIR data, it was necessary to determine the temperature threshold to calculate the canopy boundary.The canopy boundary of individual trees can be extracted from high-resolution UAVbased optical images.However, the accurate separation of overlapping crowns is impossible from optical images without using the vertical information.The multiple point clouds in LiDAR data provide additional 3D information on a canopy structure, allowing its edge to be demarcated in detail.To extract the canopy boundary, the canopy was segmented by LiDAR using the PCS algorithm, on the basis of the 3D information of the canopy (Figure 5b).The canopy segmentation of an individual tree had a good accuracy according to LiDAR data for the two plots (Table 3).The multiple point clouds in LiDAR data provide additional 3D information on a canopy structure, allowing its edge to be demarcated in detail.To extract the canopy boundary, the canopy was segmented by LiDAR using the PCS algorithm, on the basis of the 3D information of the canopy (Figure 5b).The canopy segmentation of an individual tree had a good accuracy according to LiDAR data for the two plots (Table 3).

The Impact of LAI on the Relationship between Shoot Damage Ratio and CST in Individual Trees
Figure 6a shows the significant negative relationship between LAI and CST (r = −0.9121).Under a forest canopy, the CST was easily affected by background soil, as was canopy spectrum.As LAI decreased, the soil covered by the crown was exposed to the thermal sensor.The thermal radiation from the canopy received by the sensor was the mixed result of soil, background, and canopy radiation.As the radiation from the canopy was lower than from the soil, in the TIR images, the mixed result caused the CST to increase with a rising LAI.

The Impact of LAI on the Relationship between Shoot Damage Ratio and CST in Individual Trees
Figure 6a shows the significant negative relationship between LAI and CST (r = −0.9121).Under a forest canopy, the CST was easily affected by background soil, as was canopy spectrum.As LAI decreased, the soil covered by the crown was exposed to the thermal sensor.The thermal radiation from the canopy received by the sensor was the mixed result of soil, background, and canopy radiation.As the radiation from the canopy was lower than from the soil, in the TIR images, the mixed result caused the CST to increase with a rising LAI.An insignificant correlation was revealed for the same trees (the red line in Figure 6b) because of the impact of the abnormal data with a low LAI (0.81, 0.61, and 0.62), which is shown in Figure 6b (the red points).After the elimination of the abnormal CST with a low LAI, the correlation between other CSTs and SDRs was significantly improved, and the determination coefficient (R 2 ) increased from 0.0907 to 0.6936.This confirmed that when the canopy had a low LAI, the impact of the LAI weakened the relationship between the CST and SDR by enhancing the ground radiation from the soil, resulting in an inaccurate CST.Therefore, to analyze the correlation between CST and SDR, field measurement data on Yunnan pine with high LAIs were used to ensure that the influence of LAI was weakened.

The Correlation Analysis of SDR and CST for an Individual Tree Canopy
Due to the influence of LAI, the SDR data of 30 trees with high LAIs (>1) were used to determine the relationship of SDR with CST.The canopy temperatures were extracted from TIR images using canopy areas generated from LiDAR data.The actual field data revealed the crown SDRs of the affected Yunnan pines was mainly between 10% and 20%, and there was fewer affected pines when the SDR was between 50% and 80%.This was because the survival time of an affected tree was shorter when the crown contained many dead needles as this reduced the volume of water obtained from the roots as the draw of An insignificant correlation was revealed for the same trees (the red line in Figure 6b) because of the impact of the abnormal data with a low LAI (0.81, 0.61, and 0.62), which is shown in Figure 6b (the red points).After the elimination of the abnormal CST with a low LAI, the correlation between other CSTs and SDRs was significantly improved, and the determination coefficient (R 2 ) increased from 0.0907 to 0.6936.This confirmed that when the canopy had a low LAI, the impact of the LAI weakened the relationship between the CST and SDR by enhancing the ground radiation from the soil, resulting in an inaccurate CST.Therefore, to analyze the correlation between CST and SDR, field measurement data on Yunnan pine with high LAIs were used to ensure that the influence of LAI was weakened.

The Correlation Analysis of SDR and CST for an Individual Tree Canopy
Due to the influence of LAI, the SDR data of 30 trees with high LAIs (>1) were used to determine the relationship of SDR with CST.The canopy temperatures were extracted from TIR images using canopy areas generated from LiDAR data.The actual field data revealed the crown SDRs of the affected Yunnan pines was mainly between 10% and 20%, and there was fewer affected pines when the SDR was between 50% and 80%.This was because the survival time of an affected tree was shorter when the crown contained many dead needles as this reduced the volume of water obtained from the roots as the draw of transpiration decreased.Therefore, trees with a SDR exceeding 50% can die quickly, leading to a low number of samples in this stage.The CST increased with a rising SDR, as shown by the significant correlation in Figure 7.This strong correlation (R 2 = 0.7890) indicated that UAV-based thermal remote sensing could effectively monitor the SDR of the canopy in individual trees.transpiration decreased.Therefore, trees with a SDR exceeding 50% can die quickly, leading to a low number of samples in this stage.The CST increased with a rising SDR, as shown by the significant correlation in Figure 7.This strong correlation (R 2 = 0.7890) indicated that UAV-based thermal remote sensing could effectively monitor the SDR of the canopy in individual trees.

The Effect of LAI on the Relationship between Temperature and SDR on Plot Scale
The correlation analysis of LAI and Tp revealed that LAI was also negatively correlated with temperature (r = −0.5902),which was the same as that at the canopy scale (Figure 8a).Stronger transpiration with a high number of leaves in a plot caused the Tp to decrease with a rising LAI.Additionally, the plot SDRs were calculated using the average SDR of the single trees within a plot.Figure 8b shows the non-significant correlation of SDR and Tp, on the basis of the data of 50 plots, including different LAIs from 0 to 4 (R 2 = 0.0105).
To further understand the relationship between SDR and Tp on a large scale and to explore the role of the LAI, the 20 plot data were filtered using the expression LAI > 1 for

The Effect of LAI on the Relationship between Temperature and SDR on Plot Scale
The correlation analysis of LAI and T p revealed that LAI was also negatively correlated with temperature (r = −0.5902),which was the same as that at the canopy scale (Figure 8a).Stronger transpiration with a high number of leaves in a plot caused the T p to decrease with a rising LAI.transpiration decreased.Therefore, trees with a SDR exceeding 50% can die quickly, ing to a low number of samples in this stage.The CST increased with a rising SD shown by the significant correlation in Figure 7.This strong correlation (R 2 = 0.7890) cated that UAV-based thermal remote sensing could effectively monitor the SDR o canopy in individual trees.

The Effect of LAI on the Relationship between Temperature and SDR on Plot Scale
The correlation analysis of LAI and Tp revealed that LAI was also negatively c lated with temperature (r = −0.5902),which was the same as that at the canopy scale ure 8a).Stronger transpiration with a high number of leaves in a plot caused the decrease with a rising LAI.Additionally, the plot SDRs were calculated using the average SDR of the single within a plot.Figure 8b shows the non-significant correlation of SDR and Tp, on the of the data of 50 plots, including different LAIs from 0 to 4 (R 2 = 0.0105).
To further understand the relationship between SDR and Tp on a large scale a explore the role of the LAI, the 20 plot data were filtered using the expression LAI > Additionally, the plot SDRs were calculated using the average SDR of the single trees within a plot.Figure 8b shows the non-significant correlation of SDR and T p , on the basis of the data of 50 plots, including different LAIs from 0 to 4 (R 2 = 0.0105).
To further understand the relationship between SDR and T p on a large scale and to explore the role of the LAI, the 20 plot data were filtered using the expression LAI > 1 for analysis (Figure 9).The results revealed a significant correlation in the relationship between SDR and T p after the influence of a low LAI (R 2 = 0.5532) was weakened.Therefore, the influence of the LAI on plot temperature should be considered when using TIR data to monitor forest infestation with the SDR at the pixel scale.In TIR images, pixels with a low LAI may include an excessive surface temperature due to the background radiation of soil; this will affect the correlation between the SDR and temperature, resulting in errors when monitoring forest infestations.In particular, in sparse coniferous forests where the crown needles have a small surface area and a non-aggregated canopy structure, background thermal radiation generated below the canopy can easily be captured by sensors.
to monitor forest infestation with the SDR at the pixel scale.In TIR images, pixels with a low LAI may include an excessive surface temperature due to the background radiation of soil; this will affect the correlation between the SDR and temperature, resulting in errors when monitoring forest infestations.In particular, in sparse coniferous forests where the crown needles have a small surface area and a non-aggregated canopy structure, background thermal radiation generated below the canopy can easily be captured by sensors.

The Need for Thermal Infrared Imaging Combined with High-Resolution Optical Data and LiDAR Data
The technical limitations of TIR instruments mean that the resolution of TIR images is coarser than optical and LiDAR data generated from UAV-based sensors [39].Only temperature is contained in the canopy information of TIR data, and additional data sources are needed to provide the required data regrading the tree crown.In TIR data, temperature information is available for reverse water content, stomatal conductance, evapotranspiration, and other ground-feature parameters [40][41][42].However, optical data contain more information on the physiological and biochemical parameters of plants, which could cause spectral changes [43,44].In contrast, LiDAR data contain rich spatial structural characteristics for distinguishing the geometric information of surface features, which is used mainly for extracting the structural parameters of plants at the canopy level [33,45].
Therefore, using a combination of different remote sensing data can exploit the advantages of each source and the different plant information, providing a more accurate evaluation of the health status of vegetation.The improved accuracy for monitoring forest infestations via the comprehensive use of multi-source, high-resolution remote sensing data may be conducive to quickly estimating the outbreak of forest pests in a large area and taking appropriate actions to limit the spread of the pests [30,46].The technical limitations of TIR instruments mean that the resolution of TIR images is coarser than optical and LiDAR data generated from UAV-based sensors [39].Only temperature is contained in the canopy information of TIR data, and additional data sources are needed to provide the required data regrading the tree crown.In TIR data, temperature information is available for reverse water content, stomatal conductance, evapotranspiration, and other ground-feature parameters [40][41][42].However, optical data contain more information on the physiological and biochemical parameters of plants, which could cause spectral changes [43,44].In contrast, LiDAR data contain rich spatial structural characteristics for distinguishing the geometric information of surface features, which is used mainly for extracting the structural parameters of plants at the canopy level [33,45].
Therefore, using a combination of different remote sensing data can exploit the advantages of each source and the different plant information, providing a more accurate evaluation of the health status of vegetation.The improved accuracy for monitoring forest infestations via the comprehensive use of multi-source, high-resolution remote sensing data may be conducive to quickly estimating the outbreak of forest pests in a large area and taking appropriate actions to limit the spread of the pests [30,46].

The Adverse Effects of Low LAI on the Correlation between SDR and Temperature
The analysis of the data from 15 trees revealed that the LAI was significantly negatively correlated with CST (r = −0.9121) in individual trees, and the LAI had a major impact on CST.Under a low LAI, the CST of a tree crown was affected by the radiation from the soil below the canopy; this led to mixed of crown/soil radiation that was not representative of the real temperature of the canopy.After the abnormal data was removed from three trees with low LAIs (<1), the relationship between the SDR and CST increased from 0.0907 to 0.6936.Hence, 30 trees with different SDRs were selected from the two plots to analyze the correlation between the extracted canopy temperature and the actual SDR.The great coefficient of determination (R 2 = 0.7890) indicated that TIR was able to monitor the effect of bark beetles, as was also the case at the leaf scale [29].Moreover, it performed better using TIR data than the spectral index approach, including normalized differential vegetation index (NDVI) and leaf pigment activity indices to estimate SDR (R 2 = 0.67) [22].
The negative correlation between the LAI and temperature was also expressed at the plot scale (r = −0.5902).At this scale, the relatively weak correlation was influenced by shrubs within the plot.The lower temperature of the shrubs contributed to the overall temperature of the plot; therefore, if a plot included more shrubs, the average temperature decreased.As it is difficult to measure the LAI of shrubs, LiDAR is required for correction, although the results may not be entirely accurate.As such, the effects of non-target objects on monitoring results still exist on a large scale.It is necessary to eliminate this influence using other additional data sources such as LiDAR and hyperspectral images.The combination of structure and spectral information generated from LiDAR and hyperspectral data can take advantage of different sources of information to retrieve the parameters, such as LAI and vegetation indices, for monitoring forest health for improving the accuracy of SDR detection [47,48].

Conclusions
This study confirmed the feasibility of using TIR imagery for monitoring infestations of bark beetles in Yunnan pines using combined UAV-based LiDAR data at the individual tree and plot scales.The correlation between SDR and temperature was greatly affected by low LAIs.After the data of low LAI were removed, an obvious correlation between SDR and temperature was found to exist at the single plant and plot scales.Determination coefficients of 0.7890 and 0.5532 were calculated, respectively, which confirmed the feasibility of using TIR to estimate SDRs.In this study, the detection accuracy of SDR in infected forest still had room for improvement if a three-dimensional radiative transfer (RT) model was to be used, such as RAPID, DART, and LESS.These RT models can provide more detailed data for reducing the influence of environment parameters and furthering the application of multi-source remote sensing.In further research, the scale will be expanded to pixel scale using a combination of thermal infrared and optical satellite images to improve the accuracy of monitoring large-scale regional forest infestations.

Figure 1 .
Figure 1.The location of the study sites and the flight area of the unmanned airborne vehicle (UAV).The red rectangles represent the location of two field plots, and the blue rectangle represents the UAV-based thermal and light detection and ranging (LiDAR) fight area.

Figure 1 .
Figure 1.The location of the study sites and the flight area of the unmanned airborne vehicle (UAV).The red rectangles represent the location of two field plots, and the blue rectangle represents the UAV-based thermal and light detection and ranging (LiDAR) fight area.

Figure 2 .
Figure 2. The correction of the thermal infrared (TIR) data: (a) the original TIR gray image with a digital number value, and (b) the TIR image in a plot after calibration and emissivity correction.

Figure 2 .
Figure 2. The correction of the thermal infrared (TIR) data: (a) the original TIR gray image with a digital number value, and (b) the TIR image in a plot after calibration and emissivity correction.

Figure 3 .
Figure 3.The segmentation of individual trees in the plot from LiDAR point cloud with the specific viewing angles (45°, 45°).Gray color means background, white color means no points, and other colors present segmented trees.

Figure 3 .
Figure 3.The segmentation of individual trees in the plot from LiDAR point cloud with the specific viewing angles (45 • , 45 • ).Gray color means background, white color means no points, and other colors present segmented trees.

Figure 4 .
Figure 4.The temperatures of different ground features in two plots (leaf region).The visible images of each feature (right region) are cropped from the high-resolution visible image obtained by the unmanned airborne vehicle.

Figure 4 .
Figure 4.The temperatures of different ground features in two plots (leaf region).The visible images of each feature (right region) are cropped from the high-resolution visible image obtained by the unmanned airborne vehicle.

Figure 5 .
Figure 5. (a) Thermal infrared data of the plot and temperature profiles of the canopies.(b) LiDAR data of the same plot and the extracted canopy boundaries.

Figure 5 .
Figure 5. (a) Thermal infrared data of the plot and temperature profiles of the canopies.(b) LiDAR data of the same plot and the extracted canopy boundaries.

Figure 6 .
Figure 6.The correlation between leaf area index (LAI) and canopy surface temperature (CST) in 15 trees in the plots (a).The regression analysis of shoot damage ratio and CST with the same sample of trees (b).The red points represent the abnormal CST points of three trees with low LAIs.

Figure 6 .
Figure 6.The correlation between leaf area index (LAI) and canopy surface temperature (CST) in 15 trees in the plots (a).The regression analysis of shoot damage ratio and CST with the same sample of trees (b).The red points represent the abnormal CST points of three trees with low LAIs.

Figure 7 .
Figure 7.The correlation between SDR and CST of 30 affected pines with high LAI in plots.

Figure 8 .
Figure 8.The correlation analysis between plot temperature and LAI (a) and SDR (b) in 50 plots, which were demarcated from two large plots.

Figure 7 .
Figure 7.The correlation between SDR and CST of 30 affected pines with high LAI in plots.

Figure 7 .
Figure 7.The correlation between SDR and CST of 30 affected pines with high LAI in plots.

Figure 8 .
Figure 8.The correlation analysis between plot temperature and LAI (a) and SDR (b) in 50 which were demarcated from two large plots.

Figure 8 .
Figure 8.The correlation analysis between plot temperature and LAI (a) and SDR (b) in 50 plots, which were demarcated from two large plots.

Figure 9 .
Figure 9.The correlation analysis of plot SDR and Tp with a leaf area index > 1.

Figure 9 .
Figure 9.The correlation analysis of plot SDR and T p with a leaf area index > 1.
. The Need for Thermal Infrared Imaging Combined with High-Resolution Optical Data and LiDAR Data • C, and the annual precipitation is 783.7 mm.

Table 1 .
The statistics for each tree variables of two plots (n = 409).

Table 2 .
The characteristic parameters and temperatures of four ground features.

Table 3 .
Accuracy assessment of tree crown segmentation using LiDAR data in two plots.