Next Article in Journal
Precise Measurement of Stem Diameter by Simulating the Path of Diameter Tape from Terrestrial Laser Scanning Data
Next Article in Special Issue
Long-Term Post-Disturbance Forest Recovery in the Greater Yellowstone Ecosystem Analyzed Using Landsat Time Series Stack
Previous Article in Journal
Hurricane Wind Speed Estimation Using WindSat 6 and 10 GHz Brightness Temperatures
Previous Article in Special Issue
Landsat Imagery Spectral Trajectories—Important Variables for Spatially Predicting the Risks of Bark Beetle Disturbance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping Forest Health Using Spectral and Textural Information Extracted from SPOT-5 Satellite Images

1
Key Laboratory for Silviculture and Conservation of Ministry of Education, Beijing Forestry University, Beijing 100083, China
2
Institute of Forest Resource and Information Techniques, Chinese Academy of Forestry, Beijing 100091, China
3
Survey & Planning Institute of State Forestry Administration, Beijing 100714, China
4
School of Natural Resources, West Virginia University, Morgantown, WV 26506, USA
*
Authors to whom correspondence should be addressed.
Submission received: 22 May 2016 / Revised: 18 August 2016 / Accepted: 18 August 2016 / Published: 31 August 2016
(This article belongs to the Special Issue Remote Sensing of Forest Health)

Abstract

:
Forest health is an important variable that we need to monitor for forest management decision making. However, forest health is difficult to assess and monitor based merely on forest field surveys. In the present study, we first derived a comprehensive forest health indicator using 15 forest stand attributes extracted from forest inventory plots. Second, Pearson’s correlation analysis was performed to investigate the relationship between the forest health indicator and the spectral and textural measures extracted from SPOT-5 images. Third, all-subsets regression was performed to build the predictive model by including the statistically significant image-derived measures as independent variables. Finally, the developed model was evaluated using the coefficient of determination (R2) and the root mean square error (RMSE). Additionally, the produced model was further validated for its performance using the leave-one-out cross-validation approach. The results indicated that our produced model could provide reliable, fast and economic means to assess and monitor forest health. A thematic map of forest health was finally produced to support forest health management.

Graphical Abstract

1. Introduction

Forests are the largest terrestrial ecosystems on the earth and play a significant role in providing ecological, economic and social benefits [1,2]. However, the total service that they could provide significantly depends on forest stand conditions. Pfilf et al. [3] stressed that how forest disturbance events play out depend on stand condition, and ignoring it “is tantamount to ignoring the health of the forest community.” Forests of rich species composition and complex structure (good condition) are documented to be capable of providing much more ecological services as well as timber production compared with forests that have a simple structure (poor condition) [4,5,6]. This might be the reason why, currently, irregular forest management towards complex structure and high diversity is being widely adopted and prevails [7,8,9]. Forest health, a more formal and scientific term, is normally used in forestry to describe the forest stand condition. While this term first appeared in the forestry literature in the 1980s [10,11], there was no widely accepted definition for almost 10 years. The lack of a universal definition hindered the assessment of forest health as well as the monitoring of its dynamics. In this context, O’Laughlin et al. [12] integrated the definitions of forest, ecosystem and health and finally defined forest health as a condition of forest ecosystems that sustained their complexity while providing for human needs. This definition made a great effort to combine the social, ecological and economical perspectives [13] and was adopted by the US forest service and is frequently used in the forestry literature [13,14,15].
In addition to the precise definition, how to measure forest health or how forest health could be quantitatively represented arose when we wanted to assess the current forest health conditions and monitor their dynamics. A large number of studies have been dedicated to derive quantitative measures, i.e., forest health indicators. In Europe and North America, the tree crown condition, e.g., defoliation and crown dieback, is the most widely used indicator for forest health evaluation [16,17,18]. For instance, the USDA Forest Health Monitoring (FHM) Program has systematically produced a series of crown condition indicators to quantify the forest health condition, and the related materials can be found at [19]. To assess the forest health in East China, Wang et al. [20] first calculated the five crown condition indicators live crown, crown density, crown diameter, crown dieback, and foliar transparency, and then assigned these indicators to different health categories according to the USFS (United States Forest Service) Crown Condition Classification Guide (CCCG) standard. Similarly, Wang et al. [21] reported considerable forest damage of Masson pine trees in southern China by monitoring the crown condition. The importance of the tree crown condition for representing forest health can be attributed to its significant role in regulating forest ecosystem as well as its high sensitivity to natural or anthropogenic disturbance [22,23]. Zarnoch et al. [22] documented the significance of tree crowns in regulating solar energy, nutrient recycling, precipitation distribution and moisture retention in forests. Additionally, other indicators, e.g., soil properties [24], lichen communities [25], mycorrhizal mushroom diversity [26] and faunal taxa [27] were also documented to assess forest health. All these precedents indeed provided us a promising and quantitative means to understand and monitor forest health. These indicators, however, focused on only one aspect to describe the health condition. Actually, the mechanism that determines forest health is a complex process, which might include the complicated interaction between biotic and abiotic elements. As a result, forest health might not be fully represented if we examined only one aspect. Therefore, more aspects such as the site quality, forest species composition and forest structure, which might influence the forest ecosystem process as well as forest health, should be taken into consideration for better exploration of forest health.
Regardless of the indictors we employed for forest health evaluation, the data collection process for producing these indicators through a field survey was rather difficult, expensive, time-consuming and not spatially exhaustive. Fortunately, remote sensing provided us a rapid and economic approach to obtaining such data. In fact, remote sensing has been extensively employed in estimation of forest stand attributes such as forest biomass and carbon [26,28,29,30], forest diversity [31,32,33,34,35], leaf area index [35,36,37], forest age [38,39,40] and tree height [41,42,43,44]. In addition, it is also commonly used to produce forest maps, which serve as a fundamental basis for forest management and forest inventory [45,46,47,48]. For example, using waveform Lidar, Hyde [49] produced forest structure maps, which were prerequisites for wildlife habitat analysis. Chuvieco and Congalton [50] mapped forest fire hazards using remote sensing and geographic information systems.
The most commonly used remote sensing data in forestry are airborne Lidar data and optical multispectral satellite data. Although Lidar demonstrates a very promising capability to predict forest stand attributes, especially tree height estimation [51,52,53], its application is largely restricted due to its high cost [4,54]. In contrast, the optical remote sensing methods are more cost-effective and repetitive, and their performance is also reported to be as promising as airborne Lidar for most forest attributes. For instance, Wallis et al. [55] compared the performance of Lidar and optical satellite data in modeling the spatial distribution of bird biodiversity across a complex tropical mountain forest ecosystem and found that except for phylodiversity, the optical satellite data showed almost the same efficiency as the Lidar data in predicting the Shannon diversity and a measure of community structure. These authors therefore concluded that the optical multispectral satellite data could replace costly Lidar data for modeling certain aspects of biodiversity. Similarly, Maack et al. [56] found that Stereo-VHR images showed a great potential for canopy height model generation and could be an adequate alternative to Lidar and InSAR techniques. Using high spatial resolution aerial photos (1.0 m spatial resolution), Meng et al. [57] demonstrated the capability of Fourier-based textural ordination (FOTO) indices to obtain a higher prediction quality of forest aboveground biomass compared with Lidar.
In the present study, we investigated the potential of mapping forest health conditions using SPOT-5 satellite images. The objectives of the present study were to: (1) derive a comprehensive forest health indicator using forest stand attributes from field survey data; (2) investigate the correlations between the forest health indicator and the imagery-derived spectral and textural measures; (3) build a model predicting the forest health indicator using imagery-derived measures as potential independent variables; and (4) produce a forest health map based on the developed model.

2. Data and Methods

2.1. Data Source

2.1.1. Field Data

The field survey data were obtained from the eighth Chinese National Forest Inventory (CNFI) of the Guangxi Zhuang Autonomous Region in 2010. The sampling design of the CNFI was systematic sampling with a square grid of 4 × 6 km, i.e., one plot per 4 × 6 km grid cell (Figure 1). Square plots were employed with a size of 0.067 ha. In each plot, the DBH of all trees was measured, and each tree was identified to the species level. The coordinate information for each tree was recorded for spatial analysis. In addition to the sampled trees, information concerning the sample plots, e.g., elevation, slope, aspect, soil type, soil depth, and canopy closure, was also documented. The CNFI was carried out every 5 years to monitor the forest dynamics. In the present study, 233 sampling plots were included in the satellite images, of which only 48 plots were identified as forest stands. Unfortunately, 9 plots were found to lack plot information. Therefore, 39 plots were employed to derive forest stand variables. These forests are mostly degraded secondary forests caused by historic disturbance, which are now under strict protection for self-restoration. The main management objective of these forests is to reconstruct the mixed-species irregular forest structures.

2.1.2. Remote Sensing Data and Processing

We used 3 scenes of SPOT-5 images to derive the imagery variables (independent variables). These images were captured on 21 September 2010, with the K-J numbers 275/300, 275/301, and 275/302. The multi-spectral images have 4 bands, i.e., near-infrared (0.78–0.89 mm), red (0.61–0.68 mm) and green (0.50–0.59 mm) bands at a resolution of 10 m and a shortwave infrared band (1.58–1.75 mm) at a resolution of 20 m. The panchromatic image was recorded at a resolution of 2.5 m.
Prior to analysis of the images, geometrial and atmospheric corrections were conducted by the Survey & Planning Institute of State Forestry Administration, China. Terrestrial control points were used for geometrical correction with a differential GPS. Atmospheric correction was conducted using the improved Dark Object Subtraction method proposed by Castillo-santiago et al. [58].

2.2. Forest Health Indicator

2.2.1. Forest Stand Attributes for Calculating the Forest Health Indicator

The candidate forest stand attributes used to derive the forest health indicator included both traditional stand attributes and forest structural diversity indices. The traditional forest variables, i.e., quadratic mean diameter (QMD), basal area (BA), number of trees (NT) and stand volume (SV), were calculated for each plot. Since forest structural diversity provided more details on forest structure and had an important underlying implication for formulating a forest management strategy [4,59,60], we also calculated structural diversity indices to derive the forest health indicator. Structural diversity can be subdivided into 3 categories: tree species diversity, tree dimension diversity and tree position diversity [61,62]. In the present study, we used the Shannon–Wiener index (SHI), Simpson index (SII) and Pielou’s evenness index (PI) to characterize species diversity. Tree size diversity was measured by the Gini coefficient (GC) and the standard deviation of the DBHs (SDDBH). Tree position diversity was represented by the uniform angle index (UAI), tree species intermingling index (TSII), DBH dominance index (DBHDI) and diameter differentiation index (DDI). A detailed description of these indices can be found in Meng et al. [63]. Additionally, we also included humus depth (HD) and canopy closure (CC) for the calculation of the forest health indicator. HD can be representative of site quality, whereas CC might represent the competition status of the forest stands.

2.2.2. Forest Health Indicator Derivation

We performed factor analysis to produce the forest health indicator using all 15 stand attributes mentioned above. Prior to factor analysis, Bartlett’s test of sphericity and the Kaiser–Meyer–Olkin index were conducted and calculated to investigate whether there was a correlation between these variables, which is a prerequisite for factor analysis [64,65].
Factor score and factor weight were calculated using the following formulas:
F j m = i = 1 n D m i U i j
W m = C m m = 1 n C m
where Fjm is the score for the m-th factor in plot j, Wm is the weight of the m-th factor, Dmi represents the score coefficient of the i-th variable for the m-th factor, Cm accounts for the variance explained by the m-th factor, and Uij represents the i-th standardized forest stand attributes in plot j.
The standardized forest stand attributes have a mean value of 0 and a variance of 1. The standardization is conducted using the following formula:
U i j = X i j X ¯ i S i
where Xij is the observed value for the i-th variable in plot j, and X ¯ i and Si represent the mean value and standard deviation of the i-th variable, respectively.
The forest health indicator for plot j (FHI), is calculated as follows:
FHI = F j m × W m
Because the produced FHI has both negative and positive values, it is not very convenient for comparison and classification. We therefore transformed this variable using the following formula to generate the transformed final forest health indicator (FFHI):
FFHI = H i H m i n H m a x H m i n × 10
where Hmin and Hmax are the minimum and maximum value of FHI, respectively.
The FFHI ranges from 0 to 10 and is therefore more convenient for use in assessing and monitoring forest health conditions.

2.3. Imagery-Derived Measures

2.3.1. Spectral Measures

The average surface reflectance of each plot was first extracted using the multiple spectral bands (mean_green, mean_red, mean_swir and mean_nir) as well as the panchromatic band (mean_pan). Additionally, we also calculated the 10 vegetation indices that are documented to be widely used in forestry researches. The formulas of these indices and their relevant applications in forestry are summarized in Table 1 [54,63,66,67,68,69].

2.3.2. Textural Measures

First- and second-order textural measures were extracted for each plot. We used the standard deviation of gray levels (SDGL) as the first-order textural statistics:
SDGL = 1 k i = 1 k ( x i μ ) 2
where k is the number of reflectance values in the window, and μ is the mean reflectance value. We calculated the SDGL for all multispectral reflectance bands, i.e., near-infrared, red, green and the shortwave infrared band.
In comparison to the first-order textural measures, the second-order textural measures were calculated using only the panchromatic band, which was documented to be well suited to textural analysis due to its relatively high spatial resolution [63,70]. In the present study, we employed the grey level co-occurrence matrix (GLCM) as the second-order textural measure to represent the textural features. The 8 GLCMs and their relevant studies in forestry are summarized in Table 2 [58,61,71,72].
The textural measures are multi-scale phenomena and therefore, the window size was a significant component when calculating these textural statistics [20,73]. A smaller window might not contain sufficient information about the area, whereas a larger window might result in edge effects or boarder problems [20,58]. The most common approach for determining the optimum window size is to compare the correlation coefficient between the textural statistics and the dependent variables (forest stand attributes) [63] or classification accuracy [20] at different window sizes. The optimum window size is the one that should represent a trade-off between a desirable high correlation coefficient or classification accuracy and a desirable minimum window size [58]. A window size of 9 × 9 pixels was reported to be optimal by Shaban and Dikshit [74], Castillo-Santiago et al. [58] and Meng et al. [63], who also employed SPOT-5 satellite images in their studies. Following them, we used a 9 × 9 pixel window as the optimum to derive the textural measures.

2.4. Statistical Methods

We first conducted pairwise correlation analysis to identify the image-derived measures that were significantly correlated with the FFHI. These significant image-derived variables were then employed as independent variables to build the predictive models using all-subsets regression. As multicollinearity normally occurs between remotely sensed variables [38,39,75], we used a cut-off value for the variance inflation factor (VIF) of less than 4 and also restricted the number of independent variables to 4 to avoid multicollinearity. In regression, it was assumed that there was no homoscedasticity and the residuals did not deviate significantly from normality [76]. We then produced residual plots. Additionally, the Shapiro-Wilk test and Breush-Pagan test were also performed to respectively investigate the normality and homoscedasticity of the residuals. The produced models were evaluated for precision using the coefficient of determination (R2) and the root mean square error (RMSE). The predictive models were further validated for their performance or robustness using the leave-one-out cross-validation procedure by calculating the corresponding cross-validated coefficient of determination (R2cv) and the root mean square error (RMSEcv).

2.5. Experimental Procedure

The workflow for the derivation of the forest health indicator, the extraction of imagery measures, the development of the predictive model and forest health mapping is shown in Figure 2. We first developed the forest health indicator with 15 stand variables using factor analysis. Second, spectral and textural measures were extracted from the SPOT-5 images with the optimum window size. Third, the predictive model predicating the forest health was built using the imagery-derived measures as independent variables by performing all-subsets regression. Finally, we produced the thematic map for forest health using the developed model.

3. Results

3.1. Forest Health Indicator Derivation

Bartlett’s test of sphericity (chisq = 717.23, df = 105, Sig. = 0.00) and the Kaiser–Meyer–Olkin index (MSA = 0.66) indicated that there was correlation between these 15 independent variables, which suggested factor analysis could be conducted to reduce the dimensions of the independent variables and form several comprehensive components (factors) to represent the forest health condition. The total variance explained by the factors is shown in Table 3. The cumulative variance explained by the first five components accounted for 83.9% of the total variance, which indicated that these five factors represented most of the original information (Table 3).
The factor loadings, which are the correlation coefficients between the factor and original variables (15 stand attributes), are listed in Table 4 for each factor. The variables with an absolute value of the coefficient above 0.5 are normally said to be dominating factors and allow for meaningful interpretation of the factor [77,78]. The first factor (F1) accounted for 29.4% of the total variance and was highly correlated with TSII, SHI, PI and SII. Therefore we regarded F1 as the indicator for species diversity. The second component (F2) represented 20.4% of the total variance and showed high correlations with UAI, BA, NT, SV and CC. F2 was considered to represent the competition status of the forest stand. The third component (F3) explained 19.9% of the total variance and was highly correlated with DDI, QMD and SDDBH. The fourth component (F4) accounted for 7.2% of the total variance and was only highly correlated with DBHDI. We might consider F3 and F4 as indicators of tree size diversity. Only 7.1% of the total variance was explained by the fifth component (F5), and it was only highly correlated with HD, which suggests F5 is a good indicator of site quality.
Table 5 depicts the score coefficients for each stand attribute, which can be used to calculate the factor score for each factor. Using Equations (1)–(5), the FFHI was derived for each plot, which is listed in Table 6.

3.2. Correlation Analyses

The correlation analyses between the FFHI and spectral measures are summarized in Table 7. The average surface reflectance of all bands was significantly negatively correlated with the FFHI, with correlation coefficients ranging from −0.548 to −0.606. In contrast, only half of the 10 vegetation indices, i.e., brightness, NDVI, SR, VI and SAVI, indicated a significant correlation with the FFHI. Brightness and SAVI showed a negative correlation with the FFHI, whereas a positive pattern was observed for NDVI, SR and VI. Brightness had the largest absolute correlation coefficient with respect to the FFHI (0.656) followed by VI (0.633) SAVI (0.547) and SR (0.333).
Among the textural measures, the second-order statistics demonstrated much higher correlation with the FFHI than the first-order statistics. Five of the eight second-order statistics, i.e., Glcm_contrast, Glcm_correlation, Glcm_dissimilarity, Glcm_mean and Glcm_Variance, were significantly correlated with the FFHI, among which only Glcm_correlation had a positive correlation. The absolute correlation coefficient values were 0.607, 0.548, 0.543, 0.540 and 0.320 for Glcm_mean, Glcm_dissimilarity, Glcm_Variance, Glcm_contrast and Glcm_correlation, respectively. In comparison, the first-order statistics did not show any significant correlation with the FFHI expect for SDGL_red. In summary, there were a total of 12 imagery derived measures that were significantly correlated with the FHEI, i.e. Mean_green, Mean_swir, Mean_nir, Mean_red, Mean_pan, Brightness, Glcm_contrast, Glcm_correlation, Glcm_dissimilarity, Glcm_mean, Glcm_Variance and SDGL_red. All of these measures were considered as potential independent variables for producing the predictive model.

3.3. Model Establishment and Forest Health Mapping

The predictive model and its general statistics are listed in Table 8. While 12 imagery derived measures were significantly correlated with the FFHI (Table 7), only mean_swir and mean_pan were included as independent variables due to multicollinearity. We considered this predictive model to be appropriate for estimating FHEI since it was statistically significant (p < 0.01) and its correlation coefficient (R2) was 0.47, which was very close to 0.5, the cut-off value. The performance of this linear model was further substantiated using the cross-validation scores calculated from the leave-one-out cross-validation approach (R2 = 0.43, RMSEcv = 1.804).
The residual plot of the model is presented in Figure 3, and we did not observe any particular patterns or trends. Furthermore, the Shapiro–Wilk (W = 0.97506, p-value = 0.5282) and Breush–Pagan (BP = 0.0032012, df = 1, p-value = 0.9549) tests statistically demonstrated the normal distribution and homoscedasticity of the residuals. Based on this model, the thematic map of forest health was produced (Figure 4). In this thematic map, forest health was grouped into four categories with the same FFHI interval lengths, i.e., 0–2.5, 2.5–5.0, 5.0–7.5 and 7.5–10.0.

4. Discussion

4.1. Forest Health Indicator Derivation

Forest indicators are crucial for forest health evaluation and might directly influence the final evaluation results [12]. Consequently, many studies have been conducted to investigate the selection and formulation of these indicators. For instance, O’Laughlin et al. [79] subdivided the indicators into seven categories, i.e., soil, water, vegetation, animals, ecosystem cycling, landscape patterns, and non-native plants and animals. Additionally, they also identified individual indicators for some categories. For example, microbial activity, litter dynamics and the soil productivity index were suggested to represent the soil status. Visual symptom of foliar damage, tree growth efficiency and understory vegetation were considered to be good indicators of the vegetation status. O’Laughlin and Cook [12] systematically grouped the forest health indicators into two classes, i.e., indicators for measurement of stand structure (e.g., stand density, species composition, mortality/growth rate and growth-to-removal ratio) and other ecosystem structures and processes (e.g., soil and watershed processes, tree physiology, insect populations, tree resistance to insects and micro-environmental variables as well as nutrient cycles, energy flows, and ecological processes facilitating recovery from damage). Woodall et al. [80] reported that the forest indicators for the US Forest Inventory and Analysis (FIA) included tree crown condition, lichen communities, forest soils, vegetation diversity, downed woody material, and ozone injury, which represented a compromise between budgetary constraints and field sampling efficiency. In our present study, we selected 15 stand attributes to derive a comprehensive forest health indicator, and these stand attributes covered many aspects of the forest ecosystem. For instance, QMD, BA, NT and SV accounted for the tree growth efficiency; SHI, SII, PI, GC, SDDBH, UAI, TSII, DBHDI and DDI were representative of forest structural diversity; and HD represented the soil status or site productivity. In comparison to indicators only focusing on one aspect, e.g., tree crown condition [22], soil property [24], and mycorrhizal mushroom diversity and productivity [26], the forest health indicator generated here by factor analysis using these 15 stand attributes might contain slightly more information and might better represent the forest health condition. In factor analysis, it is generally understood that the use of large samples tends to provide results such that sample factor loadings are more precise estimates of population loadings and also more stable, or less variable, across repeated sampling [81,82]. For instance, Comrey [83] and Gorsuch [84] reported that 50 was a minimum sample size for factor analysis. In the present study, there were only 39 independent plot observations, and only 15 variables were used in the factor analysis. This small sample size might make our forest health indicator less precise and stable. However, many authors reported different opinions. For example, de Winter et al. [82] found that when data are well conditioned, i.e., high loadings, low number of factors, and high number of variables, factor analysis can yield reliable results for a sample size well below 50, even in the presence of small distortions. In fact, there is considerable divergence of opinion and evidence related to the minimum sample size for generating reliable results using factor analysis, and the recommendations and findings are diverse and even contradictory [81].
While our indicator can provide a general description of forest health status, it might not be suitable for health assessment for a particular purpose because the health condition of the same forest might differ significantly with different management objectives. For example, if we used our indicator to evaluate the health status of a eucalypt planation, which is mainly grown for timber production, an unhealthy status would be returned since eucalypt plantations have extremely low species diversity. However, it should be judged healthy from a timber production perspective. In fact, Tuominen et al. [12] already pointed out that based on different perspectives, the definition of a healthy forest could even appear contradictory. Therefore, O’Laughlin et al. [79] suggested that forest scientists and managers should work together with their customer to design forest health indicators for desired conditions and particular purposes.
Although our health indicator derived from factor analysis contains much information on forest structure and diversity and might well represent the forest status, we consider that our current definition of forest health, i.e., a more complex structure characterizes a healthy forest, might be rather broad and simplistic. Rather than a forest health indicator, this indicator might be more appropriate as an indicator of forest structure. However, since forest functions are determined by forest structure [63], a more complex structure might mean multiple functions and higher stability and hence represent a healthier forest status. Additionally, because the forests we studied are mainly degraded secondary forests, which are under strict protection for self-restoration, we therefore assumed that the more complex the forest structure, the healthier the forest, though it is not applicable or even wrong in most conditions. Regardless, our indicator allowed for a general assessment and comparison of forest status through a spatial and temporal horizon, which served as significant information for formulating forest management strategies. However, we would like to stress again that this health indicator might not be applicable under certain conditions, i.e., plantations grown mainly for timber production, and its validation should be judged by the management objectives. In these conditions, this so-called forest health indicator developed in the present study only represent indicator of stand structure and should be never used for forest health assessment. Therefore, we strongly recommend that both the management objectives and the forest current status should be carefully assessed before using the methods suggested in the present study for forest health assessment. Additionally, we simply classified our forests into four groups with the same FFHI interval rather than assigning them to different health categories. This was because it was extremely difficult to determine the category thresholds. However, the criteria to define a reasonable ecosystem threshold have been documented in detail by many authors [85,86,87]. Following these criteria, further research on the determination of a reasonable forest health category is therefore strongly encouraged.

4.2. Predictive Model Development

The coefficient of determination (R2) is normally used to evaluate the modeling efficiency [88]. Ozdemir and Karnieli [61] stated that 0.5 is the threshold for selecting a reliable predictive model and argued that models with R2 values lower than 0.5 are incapable of providing reliable predictions. However, many studies reported different opinions. For instance, Murfitt et al. [89] developed seven models predicting a forest health score using seven individual vegetation indices as independent variables, and the R2 values ranged from 0.23 to 0.38. They used a model with an R2 of 0.38 to map ash health for an entire area. The individual basal area growth model for Pinus halepensis produced by Condés and Sterba [90] had an R2 of 0.36, which was increased to 0.47 by including additional random effects. Therefore, we concluded that we succeeded in producing a model capable of predicting a forest health indicator, although the R2 was only 0.47. Furthermore, the predictive ability of our model was also supported by the leave-one-out cross-validation.
While 16 imagery-derived measures (10 spectral measures and six textural measures), which were significantly correlated with the health indicator, were employed as potential independent variables to build the model, the final predictive model contained only two variables, i.e., mean_swir and mean_pan, which made the model a simple measure of the amount of forest cover. This might be due to the problem of multicollinearity. Moreover, it is also noteworthy that these two variables were spectral measures, and no textural variables entered into the model. In contrast, many researches have reported the promising performance of textural measures in forest classification as well as in predicting forest stand attributes. For instance, Kayitakire et al. [91] produced models predicting top height, circumference, stand density and age variables using the second-order texture derived from IKONOS-2 imagery, and the correlation coefficients ranged from 0.76 to 0.82. Pu and Chen [92] demonstrated that texture-based features had higher capability than spectrum-based features for estimating and mapping forest LAI using WorldView-2 images. Johansen et al. [93] reported that the classification accuracy of vegetation structure was increased by 2%–19% with the inclusion of textural measures derived from QuickBird images. These studies had at least one aspect in common, i.e., the satellite images used were of very high spatial resolution. The reason for the absence of textural measures in our model might be attributed to the relatively lower spatial resolution of the SPOT-5 images, which could not detect the spatial variability of the forest health indicator. A similar result was reported by Castillo-Santiago et al. [58], who also employed both textural and spectral measures extracted from SPOT-5 images to predict stand variables; these authors reported that textural measures were not included in the predictive models. An even worse correlation between the textural measures extracted from TM images (much coarser spatial resolution) and forest stand variables was documented by Cohen and Spies [94]. In fact, considerable literature has already demonstrated that the significance of textural measures increased with the spatial resolution of the image [95,96]. This is because the higher the spatial resolution, the more pixels represent the ground objects and therefore, the textural information becomes increasingly significant [97].
In this study, we used a regression technique to build a predictive model. Although this statistical approach has been extensively used in remote sensing studies of the environment and forestry [63,67,98,99,100], many authors do not recommend it due to its rigid assumptions of the data. In regression analysis, there are four principal assumptions, i.e., linearity and additivity of the relationship between dependent and independent variables, statistical independence of the errors, homoscedasticity of the errors and normality of the error distribution [76]. However, under most conditions, the assumptions cannot be met since ecological and remotely sensed data are highly complex and in most cases, nonlinear relationships are often observed between these two types of datasets [39]. Furthermore, in order to correct the multicorrelation of the independent variables, certain variables, which significantly account for the total variance, might be excluded from the predictive model. For instance, the textural measures in the present study, which did not enter into our final predictive model, might be excluded due to multicollinearity concerns, though some of them showed significant correlations with the derived FFHI. Therefore, using simple linear regression, we might run the risk of losing information richness. Additionally, in the case of possible non-linearity, we might only use a simple linear regression approach. In this present study, only 39 plots were used to derive the predictive model. However, many studies used a similar or even lower number of plots to produce models predicting forest attributes using variables extracted from satellite images, e.g., Means et al. [101] used 19 plots, Ozdemir and Karnieli [61] used 29 plots, and Cohen and Spies used 41 plots [93]. Consequently, we consider that 39 plots were sufficient to produce a model of moderate precision and generality.
Instead of simple linear regression, more robust statistical methods, which do not make any assumptions about the data, have been extensively documented to explore the relationship between forest stand attributes and remotely sensed data. The most frequently used approaches in forestry include regression and decision trees [102,103], artificial neutral networks [70,104] random forests [105,106,107], and support vectors [108,109]. The size of the training data set for machine learning greatly influences the stability and accuracy of the trained model [110]. Koprinska [111] recommended a ratio of at least 10 times more training instances than features. In the present study, we had 15 features, and the minimum training data set should contain at least 150 features according to the recommendation of Koprinska [111]. Unfortunately, we only had 39 plots, which prevented us from using a machine learning algorithm to produce the predictive model. If the training data set is sufficient, these robust statistical approaches are encouraged, and the minimum training data sets could be defined to create a learning curve that generates (average) model performance as a function of the training sample size [112,113].

4.3. Implication for Forest Management

Forest health management is widely applied to avoid catastrophic forest disturbance or to maintain forest stability, productivity and vitality. For instance, Oliver et al. [114] documented that thinning an overly dense stand by removing excess stems would reduce the susceptibility of the remaining trees to insects, diseases, and fires, thus reducing the potential for catastrophic fires. Using stand density management diagrams, López-Sánchez and Rodríguez-Soalleiro [115] and Castedo-Dorado et al. [116] designed optimal quantitative thinning schedules for Pseudotsuga menziesii and Pinus radiata plantations, respectively, to avoid or reduce the potential of windthrow or forest fire. A detailed review article on forest health management was proposed by Oliver et al. [114], who demonstrated the necessity of forest health management, explored the way in which forests should be managed to restore and maintain health, described the available tools and barriers associated with forest health management, and discussed the potential way to reduce the management cost and to overcome the barriers.
In the present study, we succeeded in producing the quantitative tools, i.e., the predictive model and forest health map, for forest health management. These tools provide forest managers a detailed and quantitative picture of the distribution of forest health status and therefore, provide significant insights into formulating forest management strategies. Similarly, Haywood and Stone [117] have already stated that the identification of forest health status could assist forest managers in prioritizing and formulating management strategies. For instance, forest fire monitoring should be intensified in poor health areas since these forests are more susceptible to forest fire. The thematic map could also be used to determine the area with top priority when conducting restoration activities. With the chronological thematic maps, the change in forest health can be detected, and the corresponding forest management strategies and policy can be formulated. Similarly, many studies have been conducted where remote sensing techniques have been used to assist in forest health management. For example, Xiao [118] used multispectral remote sensing data and GIS techniques to determine tree health at the University of California, Davis. Using WorldView-2 (WV2) imagery, Murfitt et al. [89] produced a remote sensing-based method for mapping ash trees undergoing various infestation stages. Both of these efforts could directly indicate or estimate the exact health status of the trees and therefore, provide solid and straightforward information for forest managers.
While our predictive model is a promising tool for prescribing forest management strategies, it should be used with great care. For example, incorrect results might be generated if we used the model to predict forest health in an area beyond which the model was developed. Further, the prediction might not be reliable if we used the model in a forest with extremely poor health, for which the model has demonstrated poor predictability. Budget-permitting, the model is encouraged to be further improved by including additional ground plots and by using very high spatial solution imagery as well as more robust statistical approaches.

5. Conclusions

Forest health assessments are of great importance for the formulation of forest management strategies as well as forest policy. However, the lack of rapid, economic and quantitative measures of forest health hinders our understanding of forest health status. In the present study, we derived a quantitative forest health indicator using 15 forest stand attributes. Additionally, we also succeeded in building a model to predict the forest health indicator using imagery-derived variables. This predictive model allowed for a rapid, economic and quantitative assessment of the forest health status, which could facilitate the decision making process of forest health management planning. Budget-permitting, the predictive model is strongly encouraged to be improved by increasing the number of survey plots and employing Lidar and other satellite imagery with very high spatial resolution.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (31370635), the project National Science and Technology Major Projects of China (21-Y30B05-9001-13/15) and the National Natural Science Foundation of China (31300532). We thank the Survey & Planning Institute of the State Forestry Administration, China, who provided data support during our research.

Author Contributions

All authors made significant contributions to the manuscript: Jinghui Meng and Shiming Li conceived, designed and performed the experiments; Shiming Li, Wei Wang, Qingwang Liu and Shiqin Xie processed the remote sensing data and derived the remote sensing parameters corresponding to the field plots; Jinghui Meng, Qingwang Liu, and Shiqin Xie analyzed the data and results; Wei Wang contributed reagents/materials/analysis tools; Jinghui Meng and Shiming Li are the main authors who developed and revised the manuscript; and Wu Ma provided important suggestions for the framework of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Davis, L.S.; Norman Johnson, K.; Bettinger, P.S.; Howard, T.E.; Alván Encinas, L.; Salazar, M.; Gretzinger, S.; Lange, G.; Schmithusen, F.; Hyde, W.; et al. Forest Management: To Sustain Ecological, Economic, and Social Values; Universidad Nacional Agraria La Molina: Lima, Perú, 2001. [Google Scholar]
  2. Buschbacher, R.J. Natural forest management in the humid tropics: Ecological, social, and economic considerations. Ambio (Sweden) 1990, 19, 253–258. [Google Scholar]
  3. Pfilf, R.J.; Marker, J.; Averill, R.D. Forest Health and Fire: An Overview and Evaluation; National Association of Forest Service Retirees: Fort Collins, CO, USA, 2002. [Google Scholar]
  4. Meng, J.; Lu, Y.; Zeng, J. Transformation of a degraded pinus massoniana plantation into a mixed-species irregular forest: Impacts on stand structure and growth in southern China. Forests 2014, 5, 3199–3221. [Google Scholar] [CrossRef]
  5. Knoke, T.; Ammer, C.; Stimm, B.; Mosandl, R. Admixing broadleaved to coniferous tree species: A review on yield, ecological stability and economics. Eur. J. For. Res. 2008, 127, 89–101. [Google Scholar] [CrossRef]
  6. Redondo-Brenes, A.; Montagnini, F. Growth, productivity, aboveground biomass, and carbon sequestration of pure and mixed native tree plantations in the caribbean lowlands of Costa Rica. For. Ecol. Manag. 2006, 232, 168–178. [Google Scholar] [CrossRef]
  7. O’Hara, K.L. The silviculture of transformation—A commentary. For. Ecol. Manag. 2001, 151, 81–86. [Google Scholar] [CrossRef]
  8. Puettmann, K.J.; Wilson, S.M.; Baker, S.C.; Donoso, P.J.; Drössler, L.; Amente, G.; Harvey, B.D.; Knoke, T.; Lu, Y.; Nocentini, S.; et al. Silvicultural alternatives to conventional even-aged forest management—What limits global adoption? For. Ecosyst. 2015, 2, 1–16. [Google Scholar] [CrossRef]
  9. Filotas, E.; Parrott, L.; Burton, P.J.; Chazdon, R.L.; Coates, K.D.; Coll, L.; Haeussler, S.; Martin, K.; Nocentini, S.; Puettmann, K.J.; et al. Viewing forests through the lens of complex systems science. Ecosphere 2014, 5, 1–23. [Google Scholar] [CrossRef]
  10. Waring, R.H. Forest, fresh perspectives from ecosystem analysis. In Vital Signs of Forest Ecosystems; Oregon State University Press: Corvallis, OR, USA, 1980; pp. 131–136. [Google Scholar]
  11. Smith, W.H. Health of north american forests: Stress and risk assessment. J. For. USA 1990, 88. [Google Scholar]
  12. O’Laughlin, J.; Cook, P.S. Inventory-based forest health indicators: Implications for national forest management. J. For. 2003, 101, 11–17. [Google Scholar]
  13. Tuominen, J.; Haapanen, R.; Lipping, T.; Kuosmanen, V. Remote Sensing of Forest Health; INTECH Open Access Publisher: Rijeka, Croatia, 2009. [Google Scholar]
  14. Lim, S.S. Development of Forest Aesthetic Indicators in Sustainable Forest Management Standards; University of British Columbia: Vancouver, BC, Canada, 2012. [Google Scholar]
  15. Turnblom, K.W. Private Forests in The Wildland-Urban Interface: Using Geographic Information Systems (GIS) to Identify Management Challenges in Eastern Washington, United States. Master’s Thesis, Utah State University, Logan, UT, USA, 2015. [Google Scholar]
  16. Lorenz, M. International co-operative programme on assessment and monitoring of air pollution effects on forests-ICP forests. Water Air Soil Pollut. 1995, 85, 1221–1226. [Google Scholar] [CrossRef]
  17. Alexander, S.A.; Palmer, C.J. Forest health monitoring in the united states: First four years. Environ. Monit. Assess. 1999, 55, 267–277. [Google Scholar] [CrossRef]
  18. D’Eon, S.; Magasi, L.; Lachance, D.; DesRochers, P. Canada’s National Forest Health Monitoring Plot Network Manual on Plot Establishment and Monitoring (Revised); Petawawa National Forestry Institute: Ontario, CA, USA, 1994. [Google Scholar]
  19. Crown Indicator Homepage. Available online: http://srsfia2.fs.fed.us/crowns/ (accessed on 11 April 2016).
  20. Wang, H.; Zhao, Y.; Pu, R.; Zhang, Z. Mapping Robinia pseudoacacia forest health conditions by using combined spectral, spatial, and textural information extracted from IKONOS imagery and random forest classifier. Remote Sens. 2015, 7, 9020–9044. [Google Scholar] [CrossRef]
  21. Wang, Y.; Solberg, S.; Yu, P.; Myking, T.; Vogt, R.D.; Du, S. Assessments of tree crown condition of two masson pine forests in the acid rain region in South China. For. Ecol. Manag. 2007, 242, 530–540. [Google Scholar] [CrossRef]
  22. Zarnoch, S.J.; Bechtold, W.A.; Stolte, K. Using crown condition variables as indicators of forest health. Can. J. For. Res. 2004, 34, 1057–1070. [Google Scholar] [CrossRef]
  23. Schomaker, M.E.; Zarnoch, S.J.; Bechtold, W.A.; Latelle, D.J.; Burkman, W.G.; Cox, S.M. Crown-Condition Classification: A Guide to Data Collection and Analysis; US Department of Agriculture, Forest Service, Southern Research Station: Asheville, NC, USA, 2007.
  24. O’Neill, K.P.; Amacher, M.C.; Perry, C.H. Soils as An Indicator of Forest Health; United States Forest Service, North Central Research Station: St. Paul, MN, USA, 2005.
  25. McCune, B. Lichen communities as indicators of forest health. Bryologist 2000, 103, 353–356. [Google Scholar] [CrossRef]
  26. Egli, S. Mycorrhizal mushroom diversity and productivity—An indicator of forest health? Ann. For. Sci. 2011, 68, 81–88. [Google Scholar] [CrossRef]
  27. Hilty, J.; Merenlender, A. Faunal indicator taxa selection for monitoring ecosystem health. Biol. Conserv. 2000, 92, 185–197. [Google Scholar] [CrossRef]
  28. Zhang, J.; Huang, S.; Hogg, E.; Lieffers, V.; Qin, Y.; He, F. Estimating spatial variation in alberta forest biomass from a combination of forest inventory and remote sensing data. Biogeosciences 2014, 11, 2793–2808. [Google Scholar] [CrossRef] [Green Version]
  29. Andersen, H.-E.; Strunk, J.; Temesgen, H.; Atwood, D.; Winterberger, K. Using multilevel remote sensing and ground data to estimate forest biomass resources in remote regions: A case study in the boreal forests of Interior Alaska. Can. J. Remote Sens. 2012, 37, 596–611. [Google Scholar] [CrossRef]
  30. Du, L.; Zhou, T.; Zou, Z.; Zhao, X.; Huang, K.; Wu, H. Mapping forest biomass using remote sensing and national forest inventory in China. Forests 2014, 5, 1267–1283. [Google Scholar] [CrossRef]
  31. Powers, R.P.; Coops, N.C.; Morgan, J.L.; Wulder, M.A.; Nelson, T.A.; Drever, C.R.; Cumming, S.G. A remote sensing approach to biodiversity assessment and regionalization of the Canadian boreal forest. Progr. Phys. Geogr. 2013, 37, 36–62. [Google Scholar] [CrossRef]
  32. Gairola, S.; Procheş, Ş.; Rocchini, D. High-resolution satellite remote sensing: A new frontier for biodiversity exploration in Indian Himalayan forests. Int. J. Remote Sens. 2013, 34, 2006–2022. [Google Scholar] [CrossRef]
  33. Laurin, G.V.; Chan, J.C.-W.; Chen, Q.; Lindsell, J.A.; Coomes, D.A.; Guerriero, L.; Del Frate, F.; Miglietta, F.; Valentini, R. Biodiversity mapping in a tropical West African forest with airborne hyperspectral data. PLoS ONE 2014, 9, e97910. [Google Scholar]
  34. Getzin, S.; Wiegand, K.; Schöning, I. Assessing biodiversity in forests using very high-resolution images and unmanned aerial vehicles. Method. Ecol. Evol. 2012, 3, 397–404. [Google Scholar] [CrossRef]
  35. Peduzzi, A.; Wynne, R.H.; Fox, T.R.; Nelson, R.F.; Thomas, V.A. Estimating leaf area index in intensively managed pine plantations using airborne laser scanner data. For.Ecol. Manag. 2012, 270, 54–65. [Google Scholar] [CrossRef]
  36. Liu, Q.; Liang, S.; Xiao, Z.; Fang, H. Retrieval of leaf area index using temporal, spectral, and angular information from multiple satellite data. Remote Sens. Environ. 2014, 145, 25–37. [Google Scholar] [CrossRef]
  37. Gao, F.; Anderson, M.C.; Kustas, W.P.; Houborg, R. Retrieving leaf area index from landsat using MODIS LAI products and field measurements. IEEE Geosci. Remote Sens. Lett. 2014, 11, 773–777. [Google Scholar]
  38. Jensen, J.; Qiu, F.; Ji, M. Predictive modelling of coniferous forest age using statistical and artificial neural network approaches applied to remote sensor data. Inter. J. Remote Sens. 1999, 20, 2805–2822. [Google Scholar]
  39. Dye, M.; Mutanga, O.; Ismail, R. Combining spectral and textural remote sensing variables using random forests: Predicting the age of Pinus patula forests in Kwazulu-Natal, South Africa. J. Spat. Sci. 2012, 57, 193–211. [Google Scholar] [CrossRef]
  40. Tooke, T.R.; Coops, N.C.; Webster, J. Predicting building ages from lidar data with random forests for building energy modeling. Energy Build. 2014, 68, 603–610. [Google Scholar] [CrossRef]
  41. Donoghue, D.; Watt, P. Using Lidar to compare forest height estimates from IKONOS and Landsat ETM+ data in Sitka spruce plantation forests. Int. J. Remote Sens. 2006, 27, 2161–2175. [Google Scholar] [CrossRef]
  42. Naesset, E. Determination of mean tree height of forest stands using airborne laser scanner data. ISPRS J. Photogramm. Remote Sens. 1997, 52, 49–56. [Google Scholar] [CrossRef]
  43. Næsset, E.; Økland, T. Estimating tree height and tree crown properties using airborne scanning laser in a boreal nature reserve. Remote Sens. Environ. 2002, 79, 105–115. [Google Scholar] [CrossRef]
  44. McCombs, J.W.; Roberts, S.D.; Evans, D.L. Influence of fusing lidar and multispectral imagery on remotely sensed estimates of stand density and mean tree height in a managed loblolly pine plantation. For. Sci. 2003, 49, 457–466. [Google Scholar]
  45. Wilkinson, G.; Folving, S.; Kanellopoulos, I.; McCormick, N.; Fullerton, K.; Megier, J. Forest mapping from multi-source satellite data using neural network classifiers—An experiment in Portugal. Remote Sens. Rev. 1995, 12, 83–106. [Google Scholar] [CrossRef]
  46. Van Coillie, F.M.; Verbeke, L.P.; De Wulf, R.R. Feature selection by genetic algorithms in object-based classification of IKONOS imagery for forest mapping in Flanders, Belgium. Remote Sens. Environ. 2007, 110, 476–487. [Google Scholar] [CrossRef]
  47. Pax-Lenney, M.; Woodcock, C.E.; Macomber, S.A.; Gopal, S.; Song, C. Forest mapping with a generalized classifier and Landsat TM data. Remote Sens. Environ. 2001, 77, 241–250. [Google Scholar] [CrossRef]
  48. Gjertsen, A.K. Accuracy of forest mapping based on Landsat TM data and a KNN-based method. Remote Sens. Environ. 2007, 110, 420–430. [Google Scholar] [CrossRef]
  49. Hyde, P.; Dubayah, R.; Peterson, B.; Blair, J.; Hofton, M.; Hunsaker, C.; Knox, R.; Walker, W. Mapping forest structure for wildlife habitat analysis using waveform lidar: Validation of montane ecosystems. Remote Sens. Environ. 2005, 96, 427–437. [Google Scholar] [CrossRef]
  50. Chuvieco, E.; Congalton, R.G. Application of remote sensing and geographic information systems to forest fire hazard mapping. Remote Sens. Environ. 1989, 29, 147–159. [Google Scholar] [CrossRef]
  51. Dubayah, R.; Sheldon, S.; Clark, D.; Hofton, M.; Blair, J.; Hurtt, G.; Chazdon, R. Estimation of tropical forest height and biomass dynamics using lidar remote sensing at La Selva, Costa Rica. J. Geophys. Res. Biogeosci. 2010, 115. [Google Scholar] [CrossRef]
  52. Nilsson, M. Estimation of tree heights and stand volume using an airborne lidar system. Remote Sens. Environ. 1996, 56, 1–7. [Google Scholar] [CrossRef]
  53. Simard, M.; Pinto, N.; Fisher, J.B.; Baccini, A. Mapping forest canopy height globally with spaceborne lidar. J. Geophys. Res. Biogeosci. 2011, 116. [Google Scholar] [CrossRef]
  54. Wolter, P.T.; Townsend, P.A.; Sturtevant, B.R. Estimation of forest structural parameters using 5 and 10 meter SPOT-5 satellite data. Remote Sens. Environ. 2009, 113, 2019–2036. [Google Scholar] [CrossRef]
  55. Wallis, C.I.; Paulsch, D.; Zeilinger, J.; Silva, B.; Fernández, G.F.C.; Brandl, R.; Farwig, N.; Bendix, J. Contrasting performance of lidar and optical texture models in predicting avian diversity in a tropical mountain forest. Remote Sens. Environ. 2016, 174, 223–232. [Google Scholar] [CrossRef]
  56. Maack, J.; Kattenborn, T.; Ewald Fassnacht, F.; Enssle, F.; Hernández Palma, J.; Corvalán Vera, P.; Koch, B. Modeling forest biomass using Very-High-Resolution data—Combining textural, spectral and photogrammetric predictors derived from spaceborne stereo images. Eur. J. Remote Sens. 2015, 48, 245–261. [Google Scholar] [CrossRef]
  57. Meng, S.; Pang, Y.; Zhang, Z.; Jia, W.; Li, Z. Mapping aboveground biomass using texture indices from aerial photos in a temperate forest of northeastern China. Remote Sens. 2016, 8, 230. [Google Scholar] [CrossRef]
  58. Castillo-Santiago, M.A.; Ricker, M.; de Jong, B.H. Estimation of tropical forest structure from SPOT-5 satellite images. Int. J. Remote Sens. 2010, 31, 2767–2782. [Google Scholar] [CrossRef]
  59. Lexerød, N.L.; Eid, T. An evaluation of different diameter diversity indices based on criteria related to forest management planning. For. Ecol. Manag. 2006, 222, 17–28. [Google Scholar] [CrossRef]
  60. Bettinger, P.; Tang, M. Tree-level harvest optimization for structure-based forest management based on the species mingling index. Forests 2015, 6, 1121–1144. [Google Scholar] [CrossRef]
  61. Ozdemir, I.; Karnieli, A. Predicting forest structural parameters using the image texture derived from Worldview-2 multispectral imagery in a dryland forest, Israel. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 701–710. [Google Scholar] [CrossRef]
  62. Pommerening, A. Approaches to quantifying forest structures. Forestry 2002, 75, 305–324. [Google Scholar] [CrossRef]
  63. Meng, J.; Li, S.; Wang, W.; Liu, Q.; Xie, S.; Ma, W. Estimation of forest structural diversity using the spectral and textural information derived from SPOT-5 satellite images. Remote Sens. 2016, 8, 125. [Google Scholar] [CrossRef]
  64. Bartlett, M.S. Tests of significance in factor analysis. Br. J. Stat. Psychol. 1950, 3, 77–85. [Google Scholar] [CrossRef]
  65. Kaiser, H.F. An index of factorial simplicity. Psychometrika 1974, 39, 31–36. [Google Scholar] [CrossRef]
  66. Germany, T. Object features: Layer values. In Ecognition Developer 7 Reference Book; Trimble Germany: Munich, Germany, 2011; pp. 262–272. [Google Scholar]
  67. Wallner, A.; Elatawneh, A.; Schneider, T.; Knoke, T. Estimation of forest structural information using rapideye satellite data. Forestry 2014, 88, 96–107. [Google Scholar] [CrossRef]
  68. Huete, A.R. A soil-adjusted vegetation index (SAVI). Remote Sens. Environ. 1988, 25, 295–309. [Google Scholar] [CrossRef]
  69. Verstraete, M.M.; Pinty, B. Designing optimal spectral indexes for remote sensing applications. IEEE Trans. Geosci. Remote Sens. 1996, 34, 1254–1265. [Google Scholar] [CrossRef]
  70. Gebreslasie, M.; Ahmed, F.; Van Aardt, J.A. Extracting structural attributes from IKONOS imagery for Eucalyptus plantation forests in Kwazulu-Natal, South Africa, using image texture analysis and artificial neural networks. Int. J. Remote Sens. 2011, 32, 7677–7701. [Google Scholar] [CrossRef]
  71. Gallardo-Cruz, J.A.; Meave, J.A.; González, E.J.; Lebrija-Trejos, E.E.; Romero-Romero, M.A.; Pérez-García, E.A.; Gallardo-Cruz, R.; Hernández-Stefanoni, J.L.; Martorell, C. Predicting tropical dry forest successional attributes from space: Is the key hidden in image texture? PLoS ONE 2012, 7, e30506. [Google Scholar] [CrossRef] [PubMed]
  72. St-Louis, V.; Pidgeon, A.M.; Radeloff, V.C.; Hawbaker, T.J.; Clayton, M.K. High-resolution image texture as a predictor of bird species richness. Remote Sens. Environ. 2006, 105, 299–312. [Google Scholar] [CrossRef]
  73. Moskal, L.M.; Franklin, S.E. Classifying multilayer forest structure and composition using high resolution, compact airborne spectrographic imager image texture. In Proceedings of the American Society of Remote Sensing and Photogrammetry Annual Conference, St. Louis, MI, USA, 23–27 April 2001.
  74. Shaban, M.; Dikshit, O. Improvement of classification in urban areas by the use of textural features: The case study of Lucknow City, Uttar Pradesh. Int. J. Remote Sens. 2001, 22, 565–593. [Google Scholar] [CrossRef]
  75. De’ath, G.; Fabricius, K.E. Classification and regression trees: A powerful yet simple technique for ecological data analysis. Ecology 2000, 81, 3178–3192. [Google Scholar] [CrossRef]
  76. Sheather, S. A Modern Approach to Regression with R; Springer Science & Business Media: Berlin, Germany, 2009. [Google Scholar]
  77. Lise, W. Factors influencing people’s participation in forest management in India. Ecol. Econ. 2000, 34, 379–392. [Google Scholar] [CrossRef]
  78. Everitt, B.; Hothorn, T. An Introduction to Applied Multivariate Analysis with R; Springer Science & Business Media: Berlin, Germany, 2011. [Google Scholar]
  79. O’Laughlin, J.; Livingston, R.L.; Thier, R.; Thornton, J.P.; Toweill, D.E.; Morelan, L. Defining and measuring forest health. J. Sustain. For. 1994, 2, 65–85. [Google Scholar] [CrossRef]
  80. Woodall, C.W.; Amacher, M.C.; Bechtold, W.A.; Coulston, J.W.; Jovan, S.; Perry, C.H.; Randolph, K.C.; Schulz, B.K.; Smith, G.C.; Tkacz, B.; et al. Status and future of the forest health indicators program of the USA. Environ. Monit. Assess. 2011, 177, 419–436. [Google Scholar] [CrossRef] [PubMed]
  81. MacCallum, R.C.; Widaman, K.F.; Zhang, S.; Hong, S. Sample size in factor analysis. Psychol. Method. 1999, 4, 84–99. [Google Scholar] [CrossRef]
  82. De Winter, J.C.F.; Dodou, D.; Wieringa, P.A. Exploratory factor analysis with small sample sizes. Multivar. Behav. Res. 2009, 44, 147–181. [Google Scholar] [CrossRef] [PubMed]
  83. Comrey, A.L. A First Course in Factor Analysis; Psychology Press: Hove, UK, 1973. [Google Scholar]
  84. Gorsuch, R.L. Factor Analysis; W.B. Saunders Company: Philadelphia, PA, USA, 1974. [Google Scholar]
  85. Schaeffer, D.J.; Cox, D.K. Establishing Ecosystem Threshold Criteria; Ecosystem Health-New Goals for Environmental Management: Washington, DC, WA, USA, 1992. [Google Scholar]
  86. Costanza, R.; Norton, B.G.; Haskell, B.D. Ecosystem Health: New Goals for Environmental Management; Island Press: Washington, DC, WA, USA, 1992. [Google Scholar]
  87. Srebotnjak, T.; Polzin, C.; Giljum, S.; Herbert, S.; Lutter, S. Establishing Environmental Sustainability Thresholds and Indicators Final Report; Ecologic Institute and SERI: Brussels, Belgium, 2010. [Google Scholar]
  88. Burkhart, H.E.; Tomé, M. Modeling Forest Trees and Stands; Springer Science & Business Media: Berlin, Germany, 2012. [Google Scholar]
  89. Murfitt, J.; He, Y.; Yang, J.; Mui, A.; de Mille, K. Ash decline assessment in emerald ash borer infested natural forests using high spatial resolution images. Remote Sens. 2016, 8, 256. [Google Scholar] [CrossRef]
  90. Condés, S.; Sterba, H. Comparing an individual tree growth model for Pinus halepensis Mill. in the spanish region of Murcia with yield tables gained from the same area. Eur. J. For. Res. 2008, 127, 253–261. [Google Scholar] [CrossRef]
  91. Kayitakire, F.; Hamel, C.; Defourny, P. Retrieving forest structure variables based on image texture analysis and IKONOS-2 imagery. Remote Sens. Environ. 2006, 102, 390–401. [Google Scholar] [CrossRef]
  92. Pu, R.; Cheng, J. Mapping forest leaf area index using reflectance and textural information derived from Worldview-2 imagery in a mixed natural forest area in Florida, US. Int. J. Appl. Earth Obs. Geoinf. 2015, 42, 11–23. [Google Scholar] [CrossRef]
  93. Johansen, K.; Coops, N.C.; Gergel, S.E.; Stange, Y. Application of high spatial resolution satellite imagery for riparian and forest ecosystem classification. Remote Sens. Environ. 2007, 110, 29–44. [Google Scholar] [CrossRef]
  94. Cohen, W.B.; Spies, T.A. Estimating structural attributes of Douglas-fir/western hemlock forest stands from Landsat and SPOT imagery. Remote Sens. Environ. 1992, 41, 1–17. [Google Scholar] [CrossRef]
  95. Lu, D.; Weng, Q. A survey of image classification methods and techniques for improving classification performance. Int. J. Remote sens. 2007, 28, 823–870. [Google Scholar] [CrossRef]
  96. Franklin, S.; Wulder, M.; Gerylo, G. Texture analysis of IKONOS panchromatic data for Douglas-fir forest age class separability in British Columbia. Int. J. Remote sens. 2001, 22, 2627–2632. [Google Scholar] [CrossRef]
  97. St-Onge, B.A.; Cavayas, F. Automated forest structure mapping from high resolution imagery based on directional semivariogram estimates. Remote Sens. Environ. 1997, 61, 82–95. [Google Scholar] [CrossRef]
  98. Lefsky, M.A.; Harding, D.; Cohen, W.; Parker, G.; Shugart, H. Surface lidar remote sensing of basal area and biomass in deciduous forests of eastern Maryland, USA. Remote Sens. Environ. 1999, 67, 83–98. [Google Scholar] [CrossRef]
  99. Gong, C.; Yu, S.; Joesting, H.; Chen, J. Determining socioeconomic drivers of urban forest fragmentation with historical remote sensing images. Landsc. Urban Plan. 2013, 117, 57–65. [Google Scholar] [CrossRef]
  100. Tian, X.; Su, Z.; Chen, E.; Li, Z.; van der Tol, C.; Guo, J.; He, Q. Reprint of: Estimation of forest above-ground biomass using multi-parameter remote sensing data over a cold and arid area. Int. J. Appl. Earth Obs. Geoinf. 2012, 17, 102–110. [Google Scholar] [CrossRef]
  101. Means, J.E.; Acker, S.A.; Fitt, B.J.; Renslow, M.; Emerson, L.; Hendrix, C.J. Predicting forest stand characteristics with airborne scanning lidar. Photogramm. Eng. Remote Sens. 2000, 66, 1367–1372. [Google Scholar]
  102. Donmez, C.; Berberoglu, S.; Erdogan, M.A.; Tanriover, A.A.; Cilek, A. Response of the regression tree model to high resolution remote sensing data for predicting percent tree cover in a Mediterranean ecosystem. Environ. Monit. Assess. 2015, 187, 1–12. [Google Scholar] [CrossRef] [PubMed]
  103. Gómez, C.; Wulder, M.A.; Montes, F.; Delgado, J.A. Modeling forest structural parameters in the Mediterranean pines of central Spain using QuickBird-2 imagery and classification and regression tree analysis (CART). Remote Sens. 2012, 4, 135–159. [Google Scholar] [CrossRef]
  104. Gopal, S.; Woodcock, C. Remote sensing of forest change using artificial neural networks. IEEE Trans. Geosci. Remote Sens. 1996, 34, 398–404. [Google Scholar] [CrossRef]
  105. Mutanga, O.; Adam, E.; Cho, M.A. High density biomass estimation for wetland vegetation using Worldview-2 imagery and random forest regression algorithm. Int. J. Appl. Earth Obs. Geoinf. 2012, 18, 399–406. [Google Scholar] [CrossRef]
  106. Baccini, A.; Friedl, M.; Woodcock, C.; Warbington, R. Forest biomass estimation over regional scales using multisource data. Geophys. Res. Lett. 2004, 31. [Google Scholar] [CrossRef]
  107. Guo, L.; Chehata, N.; Mallet, C.; Boukir, S. Relevance of airborne lidar and multispectral image data for urban scene classification using random forests. ISPRS J. Photogramm. Remote Sens. 2011, 66, 56–66. [Google Scholar] [CrossRef]
  108. Singh, S.K.; Srivastava, P.K.; Gupta, M.; Thakur, J.K.; Mukherjee, S. Appraisal of land use/land cover of mangrove forest ecosystem using support vector machine. Environ. Earth Sci. 2014, 71, 2245–2255. [Google Scholar] [CrossRef]
  109. Ueyama, M.; Ichii, K.; Iwata, H.; Euskirchen, E.S.; Zona, D.; Rocha, A.V.; Harazono, Y.; Iwama, C.; Nakai, T.; Oechel, W.C. Upscaling terrestrial carbon dioxide fluxes in alaska with satellite remote sensing and support vector regression. J. Geophys. Res. Biogeosci. 2013, 118, 1266–1281. [Google Scholar] [CrossRef]
  110. Wieland, M.; Pittore, M. Performance evaluation of machine learning algorithms for urban pattern recognition from multi-spectral satellite images. Remote Sens. 2014, 6, 2912–2939. [Google Scholar] [CrossRef]
  111. Koprinska, I. Feature selection for brain-computer interfaces. In New Frontiers in Applied Data Mining; Springer: Berlin, Germany, 2009; pp. 106–117. [Google Scholar]
  112. Figueroa, R.L.; Zeng-Treitler, Q.; Kandula, S.; Ngo, L.H. Predicting sample size required for classification performance. BMC Med. Inf. Decis. Mak. 2012, 12, 8. [Google Scholar] [CrossRef] [PubMed]
  113. Beleites, C.; Neugebauer, U.; Bocklitz, T.; Krafft, C.; Popp, J. Sample size planning for classification models. Anal. Chim. Acta 2013, 760, 25–33. [Google Scholar] [CrossRef] [PubMed]
  114. Oliver, C.D.; Ferguson, D.E.; Harvey, A.E.; Malany, H.S.; Mandzak, J.M.; Mutch, R.W. Managing ecosystems for forest health: An approach and the effects on uses and values. J. Sustain. For. 1994, 2, 113–133. [Google Scholar] [CrossRef]
  115. López-Sánchez, C.; Rodríguez-Soalleiro, R. A density management diagram including stand stability and crown fire risk for Pseudotsuga menziesii (Mirb.) Franco in Spain. Mt. Res. Dev. 2009, 29, 169–176. [Google Scholar] [CrossRef]
  116. Castedo-Dorado, F.; Crecente-Campo, F.; Álvarez-Álvarez, P.; Barrio-Anta, M. Development of a stand density management diagram for radiata pine stands including assessment of stand stability. Forestry 2009, 82, 1–16. [Google Scholar] [CrossRef]
  117. Haywood, A.; Stone, C. Mapping eucalypt forest susceptible to dieback associated with bell miners (Manorina melanophys) using laser scanning, SPOT 5 and ancillary topographical data. Ecol. Model. 2011, 222, 1174–1184. [Google Scholar] [CrossRef]
  118. Xiao, Q.; McPherson, E.G. Tree health mapping with multispectral remote sensing data at UC Davis, California. Urban Ecosyst. 2005, 8, 349–361. [Google Scholar] [CrossRef]
Figure 1. Geographical location of the study area and the images and sampling plots used in the study.
Figure 1. Geographical location of the study area and the images and sampling plots used in the study.
Remotesensing 08 00719 g001
Figure 2. A flow chart of the forest health indicator derivation, imagery measures extraction, prediction model development and forest health mapping.
Figure 2. A flow chart of the forest health indicator derivation, imagery measures extraction, prediction model development and forest health mapping.
Remotesensing 08 00719 g002
Figure 3. Plot of predicted FFHI against observed FFHI and the residual plot.
Figure 3. Plot of predicted FFHI against observed FFHI and the residual plot.
Remotesensing 08 00719 g003
Figure 4. Thematic map of the final forest health index (FFHI) for a county in the Guangxi Autonomous Region.
Figure 4. Thematic map of the final forest health index (FFHI) for a county in the Guangxi Autonomous Region.
Remotesensing 08 00719 g004
Table 1. Spectral vegetation indices derived from the SPOT-5 images.
Table 1. Spectral vegetation indices derived from the SPOT-5 images.
Vegetation IndicesFormulaReference
Brightness Brightness = RED + GREEN + NIR + SWIR 4 [66,67]
Maximum Difference Max. Diff. = m a x | R i R j | B r i g h t n e s s [63,66]
Normalized Difference Vegetation Index NDVI = N I R R E D N I R + R E D [63,67]
Simple Ratio SR = N I R R E D [63,67]
Ratio of NIR to GREEN GR = N I R G R E E N [63,67]
Ratio of GREEN to RED VI = G R E E N R E D [63,67]
Soil Adjusted Vegetation Index SAVI = N I R R E D N I R + R E D + L [63,68]
Moisture Stress Index MSI = S W I R N I R [54,63]
Standardized Vegetation Index SVR = S W I R ( R E D + G R E E N 2 ) [54,63]
Global Environment Monitoring Index GEMI = η ( 1 0.25 η ) R E D 0.125 1 R E D
η = 2 ( N I R 2 R E D 2 ) + 1.5 N I R + 0.5 R E D N I R + R E D + 0.5
[63,69]
RED, GREEN, NIR and SWIR are the surface reflectance of the red, green, near-infrared and shortwave infrared bands, respectively, and L was set to 0.5 for SAVI.
Table 2. The grey level co-occurrence matrix (GLCM) extracted from SPOT-5 images.
Table 2. The grey level co-occurrence matrix (GLCM) extracted from SPOT-5 images.
Textural MeasuresFormulaReference
Mean MEAN = i , j = 0 N 1 i P i , j [70]
Variance VAR = i , j = 0 N 1 P i , j ( i M E A N ) 2 [60,70]
Correlation COR = i , j = 0 N 1 P i , j [ ( i M E A N ) ( j M E A N ) V A R ] [70,71]
Contrast CONT = i , j = 0 N 1 P i , j ( i j ) 2 [57,60,70,71]
Dissimilarity DISS = i , j = 0 N 1 P i , j | i j | [57,70]
Homogeneity HOM = i , j = 0 N 1 P i , j 1 + ( i j ) 2 [57,70]
Angular second moment ASM = i , j = 0 N 1 P i , j 2 [57,70,71]
Entropy ENT = i , j = 0 N 1 P i , j l n P i , j [57,60,70]
Table 3. Total variance explained by the components generated using factor analysis.
Table 3. Total variance explained by the components generated using factor analysis.
ComponentTotal VariancePercentage VarianceCumulative Percentage
14.40329.429.4
23.05920.449.7
32.98519.969.6
41.0737.276.8
51.0687.183.9
Table 4. Rotated component matrix generated by factor analysis.
Table 4. Rotated component matrix generated by factor analysis.
Stand AttributesComponent
12345
TSII0.950−0.0400.197−0.061−0.008
DBHDI−0.1350.1460.2370.810−0.083
DDI0.5560.0160.690−0.013−0.069
UAI0.385−0.7090.192−0.1710.158
QMD0.1130.0870.9080.0160.061
BA0.0980.8300.4510.1190.138
NT0.0540.872−0.1120.199−0.014
SV0.0950.7260.5770.0650.180
SDDBH0.4160.0050.8770.0500.069
GC0.1250.137−0.4120.5600.157
SHI0.969−0.0110.1260.0260.002
PI0.946−0.1770.158−0.066−0.013
SII0.975−0.0020.1420.034−0.004
HD−0.0430.0290.0630.0070.973
CC−0.1290.7060.007−0.0490.005
Numbers in bold font denote a dominating indicator (factor loading ≥0.5 or ≤−0.5).
Table 5. Component score coefficient matrix generated by factor analysis.
Table 5. Component score coefficient matrix generated by factor analysis.
Stand AttributesComponents
12345
TSII0.2400.039−0.063−0.051−0.010
DBHDI−0.070−0.1510.1460.834−0.108
DDI0.049−0.0270.2200.005−0.113
UAI0.035−0.2490.083−0.0200.168
QMD−0.108−0.0630.3730.029−0.012
BA0.0120.2600.088−0.0420.069
NT0.0910.324−0.1460.015−0.037
SV−0.0160.2140.151−0.0720.104
SDDBH−0.021−0.0800.3190.0810.005
GC0.112−0.020−0.2110.5470.175
SHI0.2580.040−0.0990.0320.003
PI0.236−0.011−0.065−0.027−0.007
SII0.2580.041−0.0930.039−0.004
HD−0.010−0.035−0.041−0.0010.926
CC0.0090.286−0.056−0.207−0.022
Table 6. Final forest health indicator (FFHI) for each plot.
Table 6. Final forest health indicator (FFHI) for each plot.
PlotFFHIPlotFFHIPlotFFHIPlotFFHI
15.12114.06211.72313.62
25.94126.46222320
36.95130.49233.01332.75
46.15145.27245344.44
55.49156.56257.36352.52
610164.47265.46368.24
73.51173.84273.99377.41
84.36185.94283.38382.32
98.26192.58291.49391.12
107.81205.46304.45
Table 7. Pearson correlation coefficients between the image-derived measures and the FFHI.
Table 7. Pearson correlation coefficients between the image-derived measures and the FFHI.
Image-Derived MeasuresMeasuresCorrelation Coefficientp Value
Spectral measuresMean_green−0.599 **0
Mean_swir−0.567 **0
Mean_nir−0.548 **0
Mean_red−0.612 **0
Mean_pan−0.606 **0
Brightness−0.656 **0
Max_diff0.2450.133
NDVI0.327 *0.042
SR0.333 *0.038
GR0.1400.395
VI0.633 **0
SAVI−0.547 **0
MSI−0.2680.099
SVR−0.0200.905
GEMI−0.0220.892
Textural measuresSDGL_green−0.2970.067
SDGL_nir0.2680.010
SDGL_swir0.0620.707
SDGL_red−0.380 *0.017
SDGL_pan−0.1200.467
Glcm_contrast−0.540 **0
Glcm_correlation0.320 *0.047
Glcm_dissimilarity−0.548 **0
Glcm_entropy-0.2960.068
Glcm_homogeneity0.2950.068
Glcm_mean−0.607 **0
Glcm_ASM0.3030.061
Glcm_Variance−0.543 **0
* Correlation is significant at the 0.05 level (2-tailed); ** Correlation is significant at the 0.01 level (2-tailed).
Table 8. Regression model predicting the FFHI derived using all-subsets regression.
Table 8. Regression model predicting the FFHI derived using all-subsets regression.
Predictive ModelR2RMSEp
Model: FFHI = 16.840 − 0.053·mean_swir − 0.053·mean_pan0.471.6740

Share and Cite

MDPI and ACS Style

Meng, J.; Li, S.; Wang, W.; Liu, Q.; Xie, S.; Ma, W. Mapping Forest Health Using Spectral and Textural Information Extracted from SPOT-5 Satellite Images. Remote Sens. 2016, 8, 719. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8090719

AMA Style

Meng J, Li S, Wang W, Liu Q, Xie S, Ma W. Mapping Forest Health Using Spectral and Textural Information Extracted from SPOT-5 Satellite Images. Remote Sensing. 2016; 8(9):719. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8090719

Chicago/Turabian Style

Meng, Jinghui, Shiming Li, Wei Wang, Qingwang Liu, Shiqin Xie, and Wu Ma. 2016. "Mapping Forest Health Using Spectral and Textural Information Extracted from SPOT-5 Satellite Images" Remote Sensing 8, no. 9: 719. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8090719

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop