1. Introduction
Pine wilt disease (PWD), caused by the pine wood nematode (PWN), is a highly destructive disease that threatens the forests of East Asia [
1,
2,
3]. In China, since its first discovery at the Sun Yat-sen Mausoleum of NanJing in 1982, PWD has rapidly spread and caused irreparable losses [
4,
5]. The pathogenic pine wood nematode causes pine wilt rapidly and its pathogenesis is complex, which has not been fully understood [
6].
Traditional field survey methods to detect PWD are time-consuming and labor-intensive [
7]. Therefore, a faster and more effective PWD detection method is needed. Remote sensing are non-contact detection technologies that use spectral information from sensors to indirectly obtain characteristic of target objects [
8]. Remote sensing has been widely used in many fields, including flood monitoring [
9], land use and land cover survey [
10], and plant disease identification [
11,
12,
13].
After being infected with PWD, changes in physiological variables and leaf cell structure are reflected in the spectral characteristics, which can be used to infer PWD infection status [
8,
14]. On this basis, there have been numerous studies using multispectral remote sensing to identify trees with late stage PWD infections [
15,
16,
17]. However, in the early infection stage, multispectral remote sensing cannot distinguish the subtle differences in spectral responses between infected and non-infected trees. Compared with multispectral remote sensing, hyperspectral remote sensing has a much higher spectral resolution, which gives it a significant advantage in distinguishing among subtle spectral changes. Kim et al. used hyperspectral data to determine the optimal vegetation indices (VIs) for detecting PWD in the early infection stage [
7]. The Green-Red Spectral Area Index (GRSAI) they established was able to identify PWD much earlier than other VIs. Zhang et al. established a genetic algorithm-partial least squares regression (GA-PLSR) model based on the hyperspectral data to predict PWD [
16]. Yu et al. integrated unmanned aerial vehicle (UAV) based hyperspectral images and ground-based data to determine the accuracy of PWD detection models established by VIs, red edge parameters, moisture indices, and their combinations [
18]. Their results showed that the model combining all parameters had the highest accuracy.
Continuous monitoring of trees infected with PWD using remote sensing, especially in the early infection stage, can capture dynamic changes in spectral information and help detect PWD before it spreads to other pine trees. Traditional satellite remote sensing is not suitable for the early detection of PWD due to its fixed return visit period and lack of flexibility [
18]. In the past decade, UAV-based remote sensing technology has been widely used in forest disease detection due to its high flexibility, efficiency, and spatial resolution [
19,
20,
21]. Zhou et al. used UAV-based RGB images to segment and detect individual trees infected by PWD. By using adaptive local threshold selection methods, infected trees in grayscale images could be automatically segmented according to the vegetation index (VEG) with an accuracy of 90% [
22]. Deng et al. set up a deep learning framework using faster region convolutional neural networks to detect PWD and the model accuracy reached 90% [
23]. The emergence of UAV enables continuous monitoring which facilitates the early detection of PWD, and also enables researchers to obtain images with different spatial resolutions.
Although many studies have used UAV-based data to detect forest diseases, few have focused on the impact of spatio-temporal scales on forest disease monitoring, especially in PWD detection. Guo et al. used UAV-based hyperspectral images to detect wheat yellow rust disease at the field scale [
24]. They resampled 1.2 cm spatial resolution images at 3 cm, 5 cm, 7 cm, 10 cm, 15 cm, and 20 cm, and determined that the 10 cm spatial resolution was optimal for detecting wheat yellow rust disease. Jonathan et al. used UAV images, including original images and resampled images, to monitor disease outbreak in mature Pinus radiata D. simulated using herbicides [
25]. They found that 1 m was the optimal spatial resolution for detecting simulated forest disease in both small and large tree clusters. Zeng et al. used the RGB images of unusual dead pines (dead not due to PWD) collected by UAV at 430 m and 700 m altitudes to monitor PWD [
26]. Their results showed that images collected at 430 m and 700 m could both be used to accurately identify unusual dead trees. However, the above studies have only used natural color or multi-spectral images for forest disease detection. In addition, they used resampling to obtain images with different spatial resolutions instead of collecting images at different flight altitudes. However, resampling just simulates the real value through the algorithm, which is still different from the real value and may produce erroneous results.
In this study, UAV-based hyperspectral images collected at different spatio-temporal scales (different spatial resolutions and infection stages) were used to assess the development of PWD. The objectives are: (1) To analyze the impact of spatio-temporal scales on PWD detection; (2) determine the optimal spectral bands, vegetation indices (VIs), and spatial resolution for PWD detection in early, mid-, and late infection stages.
4. Discussion
As the PWD develops, the J-M distance of wavelengths to detect PWD increased. Moreover, this study found that the 401–430 nm, 614–624 nm, 655–679 nm, and 755–766 nm wavelengths had better separability than other wavelength regions in PWD detection. These wavelength regions are located in the blue, red, and red edge regions. Similarly, Iordache et al. also identified the 400 nm and 670 nm wavelengths as suitable for PWD detection, and their proposed 750 nm was close to the 755–766 nm observed in this study [
8]. Because of the small change in the near infrared region, there were no selected wavelengths in that region. In addition, any data redundancy in the hyperspectral imagery will increase the costs of data processing during practical applications. This can be addressed by programming the camera to scan only the characteristic bands with high J-M distances identified in this study. A wavelength specific camera would greatly simplify data processing and reduce costs.
We selected eight existing VIs and two modified VIs to create a random forest-based PWD detection model which showed high detection accuracy in the early (72.05–79.48%), mid- (83.71–89.59%), and late infection stages (96.81–99.28%). However, the random forest is not a trustworthy machine learning algorithm for variable importance. Compared with manually extracted shallow features, deep learning can automatically explore higher-dimensional information and features of hyperspectral images, which will be useful for future research seeking to improve the early detection of PWD. Indeed, there have been many studies using deep learning for PWD detection [
17,
48,
49].
Outbreaks of PWD in Nanjing, China usually occur from 6 June to 17 October each year. To obtain hyperspectral data of different infection stages, we manually inoculated Masson pine with PWN to simulate the natural infection by PWN and the development of PWD. Considering the strong spread ability of PWD, we only selected four sample trees for PWN inoculation to prevent the uninfected local coniferous forests from acquiring the PWD disease inoculated by this experiment. In addition, more attention should be paid to the influence of other factors such as tree species, terrain, and background on detection.
The optimal spatial resolution for PWD detection increased as PWD developed. In this study, the 10 cm, 8 cm, and 4 cm spatial resolutions were selected as the optimal spatial resolution in early, mid-, and late infection stages, respectively. In the early infection stage, the PWD detection accuracy was low and susceptible to influence from similar objects. Increases in spatial resolution at this stage did not effectively improve the PWD detection accuracy. However, in the mid- and late infection stages, the differences between infected and non-infected samples became more obvious, and higher spatial resolutions made it easier to capture these differences. While the highest spatial resolution in our research was 2 cm, it did not have the highest accuracy for any PWD infection stage, which may have been related to the canopy size of infected trees. Jonathan et al. and Guo et al. also produced similar results [
24,
25]. Limited by the permitted flying height in the study area, we were only able to collect hyperspectral images below 300 m. Future work should obtain a larger variety of hyperspectral images containing different canopy sizes to explore the impact of canopy size on the selection of optimal spatial resolution in PWD detection, especially in the early infection stage before symptoms are visible.
5. Conclusions
The UAV-based hyperspectral images collected at different infection stages and at different spatial resolutions were used to assess the impact of spatio-temporal scales on PWD detection. The following results were obtained. (1) The separability in the visible region was higher than that of near infrared and red edge regions at different spatio-temporal scales. There was a valley in the J-M distance in the green region and two J-M distance peaks in blue and red regions. The J-M distances at the lower spatial resolutions (8 cm, 10 cm, and 12 cm) were smaller than those of higher spatial resolutions (2 cm, 4 cm, and 6 cm) in the early infection stage, but greater in the mid-infection stage. In the late infection stage, the J-M distances of higher spatial resolutions were greater than those of lower spatial resolutions in near infrared region, but smaller in the visible region. (2) The PWD detection accuracies were 72.05–79.48%, 83.71–89.59%, and 96.81–99.28% with peak accuracies at 10 cm, 8 cm, and 4 cm spatial resolutions during the early, mid-, and late infection stages, respectively. (3) The GNDVI and REP were confirmed as the optimal VIs in early and mid-infection stages, respectively. However, in late infection stage, most VIs had high feature importance, so there was no optimal VI. This research focused on the impact of spatial-temporal scale on PWD detection, and has provided new ideas for future studies. This will play an important role in improving the effectiveness of early PWD detection. Identifying the optimal VIs for PWD detection will help in reducing the loss of forest resources caused by PWD. To apply these results at a larger scale, future research should pay more attention to the impact of other tree species, terrain, and background on forest disease detection.