Classification of Tropical Forest Tree Species Using Meter-Scale Image Data

Accurate classification of tropical tree species is critical for understanding forest habitat, biodiversity, forest composition, biomass, and the role of trees in climate variability through carbon uptake. The aim of this study is to establish an accurate classification procedure for tropical tree species, specifically testing the feasibility of WorldView-3 (WV-3) multispectral imagery for this task. The specific study site is a defined arboretum within a well-known tropical forest research location in Costa Rica (La Selva Biological Station). An object-based classification is the basis for the analysis to classify six selected tree species. A combination of pre-processed WV-3 bands were inputs to the classification, and an edge segmentation process defined multi-pixel-scale tree canopies. WorldView-3 bands in the Green, Red, Red Edge, and Near-Infrared 2, particularly when incorporated in two specialized vegetation indices, provide high discrimination among the selected species. Classification results yield an accuracy of 85.37%, with minimal errors of commission (7.89%) and omission (14.63%). Shadowing in the satellite imagery had a significant effect on segmentation accuracy (identifying single-species canopy tops) and on classification. The methodology presented provides a path to better characterization of tropical forest species distribution and overall composition for improving biomass studies in a tropical environment.


Introduction
Accurate identification of tropical forest species would support a more accurate measure of several important species-dependent environmental variables, such as above-ground biomass and carbon uptake [1]. Currently, large uncertainties exist in estimates of biomass and carbon uptake for tropical forests [2]. An improved species inventory leading to better biomass estimates could improve current carbon budget measurements, leading to more accurate carbon offset programs [3]. The use of remote sensing imagery for forest analysis has a long history, from the use of the legacy Landsat systems [4], to contemporary imaging systems such as Landsat 8 [5], Sentinel 2, and high-resolution commercial imaging systems [1,6]. In addition, many studies have attempted to use remotely sensed imagery for tree type identification in complex tropical forest assemblages [4,7].
Many pixel-based classification studies have used several different classification processes (e.g., maximum likelihood, spectral angle mapper, support vector machine, random forest, etc.) to determine with a Wilks' Lambda test to analyze WorldView-3 bands and 14 distinct Spectral Vegetation Indexes (SVIs) in their discriminatory power to differentiate tropical forest tree types in Costa Rica.
Any process of classification of tropical vegetation would need to account for a variety of variables that control intra-species and inter-species variations [11]. Seasonality is an important consideration in any study of the tropical rainforest. Hesketh et al. [23] showed that all data used for classification purposes should be constrained to one of the seasons (wet or dry) as inter-seasonal variability is typically low (with a dry season preference due to showing more differences between features consistently through the season) but variability between seasons is high. Castro et al. [24] also reported this effect, where classifying using data from various seasons or sites can reduce classification accuracy substantially. In addition, lianas can significantly skew the tree crown spectral signature for a given species at the leaf level, and possibly at, the crown level [11,23].
In this study, we build on our previous research [20,25], which determines effective data types for species differentiation and validates WorldView-3 as a viable data source for differentiating tree species respectively, to move the current research forward by utilizing very high resolution image-derived data products in a rule-set object-based classification for accurate identification of tropical forest species. This process contrasts to the complex approaches stated above, where affordable multispectral imagery and a simple, straightforward approach is utilized to differentiate tree species. This approach employed the information from several WorldView-3 image bands and two WorldView-3 image spectral vegetation indices derived in Cross et al. [20] to identify six different tropical forest species in Costa Rica. All data and imagery collected was constrained to the dry season minimizing any variance [23]. Illumination and view angle corrections, and the application of an atmospheric compensation procedure, assisted in creating an accurate forest canopy reflectivity image data set. Guidance on the makeup of the segmentation settings [9,18,26,27] allowed an accurate classification of tree species within the study site.

Materials and Methods
We have developed a processing and validation methodology for this study and our previous work [20,25] that applies corrected WorldView-3 multispectral imagery to extract a species canopy map. Tree species identification is achieved through an object-based classification. This is summarized in a flowchart ( Figure 1) and discussed in detail below.
Remote Sens. 2019, 11, x FOR PEER REVIEW 3 of 18 consideration in any study of the tropical rainforest. Hesketh et al. [23] showed that all data used for classification purposes should be constrained to one of the seasons (wet or dry) as inter-seasonal variability is typically low (with a dry season preference due to showing more differences between features consistently through the season) but variability between seasons is high. Castro et al. [24] also reported this effect, where classifying using data from various seasons or sites can reduce classification accuracy substantially. In addition, lianas can significantly skew the tree crown spectral signature for a given species at the leaf level, and possibly at, the crown level [11,23].
In this study, we build on our previous research [20,25], which determines effective data types for species differentiation and validates WorldView-3 as a viable data source for differentiating tree species respectively, to move the current research forward by utilizing very high resolution imagederived data products in a rule-set object-based classification for accurate identification of tropical forest species. This process contrasts to the complex approaches stated above, where affordable multispectral imagery and a simple, straightforward approach is utilized to differentiate tree species. This approach employed the information from several WorldView-3 image bands and two WorldView-3 image spectral vegetation indices derived in Cross et al. [20] to identify six different tropical forest species in Costa Rica. All data and imagery collected was constrained to the dry season minimizing any variance [23]. Illumination and view angle corrections, and the application of an atmospheric compensation procedure, assisted in creating an accurate forest canopy reflectivity image data set. Guidance on the makeup of the segmentation settings [9,18,26,27] allowed an accurate classification of tree species within the study site.

Materials and Methods
We have developed a processing and validation methodology for this study and our previous work [20,25] that applies corrected WorldView-3 multispectral imagery to extract a species canopy map. Tree species identification is achieved through an object-based classification. This is summarized in a flowchart ( Figure 1) and discussed in detail below.

Study Site
La Selva Biological Station in Costa Rica (hereafter, 'La Selva') was chosen as the study area for this research due to its accessible location, extensive trail and tower infrastructure, and excellent support staff. La Selva also contains the 3.5 ha Holdridge Arboretum (hereafter, 'Arboretum'), a managed research area within La Selva established in 1968 (Figure 2), which is the study site for this research effort. The Arboretum contains approximately 929 cataloged plants and trees, with the most recent census occurring in August 2016-March 2017. Within the managed area, 727 trees represent 185 native species [28]. The Arboretum is the focus of many research efforts at La Selva, and provides a baseline for taxonomic studies of tropical tree species within Central America.

Study Site
La Selva Biological Station in Costa Rica (hereafter, 'La Selva') was chosen as the study area for this research due to its accessible location, extensive trail and tower infrastructure, and excellent support staff. La Selva also contains the 3.5 ha Holdridge Arboretum (hereafter, 'Arboretum'), a managed research area within La Selva established in 1968 (Figure 2), which is the study site for this research effort. The Arboretum contains approximately 929 cataloged plants and trees, with the most recent census occurring in August 2016-March 2017. Within the managed area, 727 trees represent 185 native species [28]. The Arboretum is the focus of many research efforts at La Selva, and provides a baseline for taxonomic studies of tropical tree species within Central America.
A logical data choice for achieving the stated goals in this study would be either a hyperspectral system [29,30] and/or a LiDAR system [31,32], as these imaging systems have proven to be excellent choices for species discrimination/identification, especially in complex forest areas [11,33]. The purpose for using WorldView-3 (a multi-spectral sensor imaging system) for this study is to determine if a more cost-effective image data set can be effective in a tropical forest setting. The advantages of the WorldView-3 sensor are a large area collection, good revisit times over broad areas of interest, and overall cost advantages over specialized aircraft-based sensor systems. While the WorldView-3 sensor does not match the spatial resolution of LiDAR or spectral resolution of a hyperspectral sensor, there is sufficient spatial and radiometric resolution to collect many intra-crown pixels within a tree crown, and a sufficient number of image bands for spectral characterization of a particular species [25].  Table 1 for imagery specifications.

WorldView-3 Imagery Acquisition and Preparation
DigitalGlobe provided the WorldView-3 imagery for this study, and the acquisition from November 11, 2014 was a nearly cloud-free scene from the end of the dry season (Table 1), maximizing our potential of tree species identification [23,24]. We used DigitalGlobe's proprietary Atmospheric Compensation process (AComp), a physics-based compensation schema that uses observed in-scene pixel spectra for its correction procedure [34] that accounts for water vapor in the atmospheric column, to atmospherically correct the image. The derived pixel-based Aerosol Optical Depth (AOD) information was applied into a radiative transfer schema for the imagery bands collected and the output  Table 1 for imagery specifications.

WorldView-3 Imagery and Ground Truth Data Preparation
A logical data choice for achieving the stated goals in this study would be either a hyperspectral system [29,30] and/or a LiDAR system [31,32], as these imaging systems have proven to be excellent choices for species discrimination/identification, especially in complex forest areas [11,33]. The purpose for using WorldView-3 (a multi-spectral sensor imaging system) for this study is to determine if a more cost-effective image data set can be effective in a tropical forest setting. The advantages of the WorldView-3 sensor are a large area collection, good revisit times over broad areas of interest, and overall cost advantages over specialized aircraft-based sensor systems. While the WorldView-3 sensor does not match the spatial resolution of LiDAR or spectral resolution of a hyperspectral sensor, there is sufficient spatial and radiometric resolution to collect many intra-crown pixels within a tree crown, and a sufficient number of image bands for spectral characterization of a particular species [25].

WorldView-3 Imagery Acquisition and Preparation
DigitalGlobe provided the WorldView-3 imagery for this study, and the acquisition from November 11, 2014 was a nearly cloud-free scene from the end of the dry season (Table 1), maximizing our potential of tree species identification [23,24]. We used DigitalGlobe's proprietary Atmospheric Compensation process (AComp), a physics-based compensation schema that uses observed in-scene pixel spectra for its correction procedure [34] that accounts for water vapor in the atmospheric column, to atmospherically correct the image. The derived pixel-based Aerosol Optical Depth (AOD) information was applied into a radiative transfer schema for the imagery bands collected and the output was an atmospherically compensated image of surface reflectance values per pixel [35]. This was an essential step in acquiring a true surface reflectance measure from the WorldView-3 imagery.
Due to the non-Lambertian characteristics of tree canopies, and with the imagery collected at a non-nadir imaging view angle (Table 1), a Bidirectional Reflectance Factor (BRF) correction was required to convert data to an on-nadir view [36][37][38]. Reflectance anisotropy can vary by up to 30% within a closed forest canopy [39] and high spatial resolution imagery with variable view angles (WorldView-3) can be especially affected [40]. Average BRF correction values ( Table 2) for different wavelength regions derived from Breunig et al. [41] for a typical mixed tropical tree canopy, and supporting information from [42], provided the correction parameters needed for the imagery. Because of the off-nadir image collection and the inherent displacement error due to overall tree canopy height, a georectification process was performed to improve the locational accuracy of trees within the study area. A Trimble GeoExplorer 2008 GPS with a Trimble Zephyr 2 external backpack antenna, utilizing differential corrections through the Satellite Based Augmentation System and employing multipath correction, provided high positional accuracy ground control points (most points were better than 1-m positional accuracy). An image-to-map nearest-neighbor rectification procedure within ENVI provided appropriate geometrically corrected imagery.

Ground Truth Data Collection and Processing
We established ground truth of the tree species studied by finding well-exposed examples of the selected tree species crowns in the field outside of the Arboretum study area. Six select tree species  (Table 3). All tree samples for ground truth were independent from the Arboretum study site (but still within the confines of La Selva). Tree crown samples were of sufficient size to be at the top of the forest canopy structure and visible in the WorldView-3 image. The availability of each species within the Arboretum and the clear identification of the tree crown within the imagery determined the tree species selection for this study. The ground truth data collected outside of the Arboretum (Table 3), determined which canopy trees would be a part of the classification analysis within the Arboretum. Figure 3 illustrates two typical crown data extractions from the WorldView-3 imagery for ground truth purposes.
Remote Sens. 2019, 11, x FOR PEER REVIEW 6 of 18 The availability of each species within the Arboretum and the clear identification of the tree crown within the imagery determined the tree species selection for this study. The ground truth data collected outside of the Arboretum (Table 3), determined which canopy trees would be a part of the classification analysis within the Arboretum.    GPS positions of the selected tree observations, converted to an ESRI Shapefile, provided accurate geolocation for matching the field-identified tree to the proper canopy visible in the WorldView-3 imagery. Additional ESRI Shapefiles provided by the staff at La Selva (trail locations, trail signs, streams and rivers, etc.) assisted in locating selected tree species within the image. An extraction of the full-crown pixel clusters yielded mean image-derived spectra for each crown. Figure 4 illustrates the average reflectance per tree crown for each WorldView-3 band of the ground truth of the six tree species studied, displaying the variance in the reflectivity values between species. These average reflectivity data values and their specific characteristics are typical of trees within the dry season [25].
Remote Sens. 2019, 11, x FOR PEER REVIEW 7 of 18 imagery. Additional ESRI Shapefiles provided by the staff at La Selva (trail locations, trail signs, streams and rivers, etc.) assisted in locating selected tree species within the image. An extraction of the fullcrown pixel clusters yielded mean image-derived spectra for each crown. Figure 4 illustrates the average reflectance per tree crown for each WorldView-3 band of the ground truth of the six tree species studied, displaying the variance in the reflectivity values between species. These average reflectivity data values and their specific characteristics are typical of trees within the dry season [25]. In addition to the spectral reflectivity for each species, an extracted Canopy Average of the Arboretum was also included in Figure 4. This represents the collective response of canopy vegetation over the study area, and is an important data input to a vegetation index used in this study. A sampling of most of the non-shadow area of the Arboretum yielded a collection of 9750 pixels. Avoiding shadow areas was a priority during the extraction of the Canopy Average to ensure that the value calculated included mostly tree canopies.
An additional data set used for ground truth was an Arboretum inventory data file, which includes precise locations of canopy tree species [28], and was used for determining the accuracy of the classification in this study. A detailed catalog for the Arboretum describes all trees according to species, size measured in Diameter at Breast Height (DBH), and precise geolocation referenced from a permanent 25 m x 25 m grid (azimuth and distance from each grid post). The staff at La Selva produced In addition to the spectral reflectivity for each species, an extracted Canopy Average of the Arboretum was also included in Figure 4. This represents the collective response of canopy vegetation over the study area, and is an important data input to a vegetation index used in this study. A sampling of most of the non-shadow area of the Arboretum yielded a collection of 9750 pixels. Avoiding shadow areas was a priority during the extraction of the Canopy Average to ensure that the value calculated included mostly tree canopies.
An additional data set used for ground truth was an Arboretum inventory data file, which includes precise locations of canopy tree species [28], and was used for determining the accuracy of the classification in this study. A detailed catalog for the Arboretum describes all trees according to species, size measured in Diameter at Breast Height (DBH), and precise geolocation referenced from a permanent 25 m × 25 m grid (azimuth and distance from each grid post). The staff at La Selva produced an ESRI Shapefile (UTM Zone16N grid, meters) of the Arboretum catalog for easy integration with the imagery data [28].

Analysis
A series of preliminary statistical evaluation procedures helped determine which image data would be most effective for classification. Next, an imagery segmentation process defined tree crowns within the Arboretum, and a rule-set object-based classification approach with the appropriate data input guided the identification of individual tree species. Finally, an error matrix provided a determination of the accuracy of the classification process.

Data Selection
Typically, the most important imagery bands for vegetation are in the mid-to long-wavelength visible range showing the variation in chlorophyll activity, and near infrared, showing the variation in structure and water content. [44]. WorldView-3 has additional image bands in the Blue, Yellow, Red Edge, and two bands in the near-infrared (NIR-1 and NIR-2), providing a greater potential for differentiating vegetation [45]. Chlorophyll absorption is pronounced in the Blue and Red bands [46,47], and the Yellow band holds promise in providing additional information for species identification [48]. The Red Edge band characterizes plant health and is sensitive to Leaf Area Index [47,49,50]. An additional near-infrared band allows more information about vegetation water content and overall structure.
In addition to the sensor band information, we analyzed two Spectral Vegetation Indices (SVIs) to determine their ability to classify vegetation. The use of SVIs has been well documented for determining a variety of vegetation parameters, including chlorophyll production, gross primary productivity, leaf area index, vegetation type, and biomass estimates [51][52][53][54][55]. The creation of these particular SVIs were driven by better tree differentiation within a dense tropical forest canopy [20].
The WorldView Red Edge Slope Weighted Index (WV-RESWI) provides a measure of Red Edge band intensity and a measure of the reflectivity slope between the Red and the Near-IR1 bands of Worldview-3 (Equation (1)). Variation in red absorption by chlorophyll and infrared reflectivity by plant structure are important reflective characteristics that define specific tree species [46,51]: The second specialized SVI used in this study is the WorldView Average Canopy Reference Index (WV-ACRI). It is a measure of the differentiation of a specific tree species from the overall reflective response from a complex forest canopy (Equation (2)). The overall infrared response is calculated using the average of the near-infrared bands in WorldView-3 (Near-IR1, Near-IR2), and the visible response is a combination of the Green, Yellow, and Red Edge band reflective measurements as they compare to the Canopy Average for each of those bands. In the equations below, ACnir1 refers to average canopy Near-IR1 band, ACgrn refers to average canopy Green band, ACre refers to the average canopy Red Edge Band, etc. Canopy Average values for each WorldView-3 band provided the data input necessary for the calculation of WV-ACRI for the Arboretum. A discriminant analysis (DA) was performed to determine which of the eight multispectral WorldView-3 imagery bands (independent values) and the two SVIs above, had the most discriminatory Remote Sens. 2019, 11, 1411 9 of 18 power for the tree types (specific classes) studied [56]. As a part of the DA, a Wilks' Lambda test evaluated the discriminatory power of the independent variables [12,29], providing objective clarity as to which WorldView-3 bands and Vegetation Indices are most important for characterizing and differentiating tropical tree species. This test can indicate an initial assessment of group membership from the independent samples and insight into the importance of each independent variable. When the value of Wilks' Lambda for a value is small, a higher discriminatory ability is realized and corresponds to a statistically greater separability between the classified groups [57]. This is a common use of the Wilks' Lambda procedure in image analysis [12,29]. Table 4 shows the results of the Wilks' Lambda test for the six tree species studied on all eight bands of WorldView-3 imagery and the two SVIs. The Wilks' Lambda scores demonstrated significant discriminatory power in the traditional remote sensing bands typically used for vegetation analysis [44]. All bands and SVIs were significant and were below the threshold of α = 0.05. Both of the specialized SVIs constructed specifically for the WorldView-3 sensor performed as well or better than individual imagery bands in their discriminatory power. A statistical correlation (Table 5) provided additional information of the importance of each band/index studied. The optimum independent data grouping for a classification would be data with highest discriminatory values and smallest correlation between the variables [12,57].  Figures 5 and 6 show the WorldView-3 band reflectivity and SVI values for the six tree species. Figure 5 shows the visible portion of the spectrum, from the Coastal Band to the Red Band, and Figure 6 illustrates the range of values from the Red Edge Band through the two SVIs analyzed. The graphs visualize the reflective variance, as the visible portion of the spectrum ( Figure 5) is approximately one order of magnitude smaller than the other reflectivity values and SVI values ( Figure 6).  Both SVIs provide an additional important set of information for species differentiation. The necessity of the additional SVIs is evident in their ability to separate tropical forest species that are  Both SVIs provide an additional important set of information for species differentiation. The necessity of the additional SVIs is evident in their ability to separate tropical forest species that are Both SVIs provide an additional important set of information for species differentiation. The necessity of the additional SVIs is evident in their ability to separate tropical forest species that are extremely close in spectral response through the original WorldView-3 bands, especially in the visible portion of the spectrum and overlap of responses in the near-infrared bands.

Object-Based Classification
A full crown measure of reflectivity for each species (Table 3) is the basis for the classification within the Arboretum. Each of the tree crowns in the study area were organized into specific object segments comprising similar digital pixel values (Figure 7), having defined edges that separate each group from other distinct groupings [17].
Remote Sens. 2019, 11, x FOR PEER REVIEW 11 of 18 extremely close in spectral response through the original WorldView-3 bands, especially in the visible portion of the spectrum and overlap of responses in the near-infrared bands.

Object-Based Classification
A full crown measure of reflectivity for each species (Table 3) is the basis for the classification within the Arboretum. Each of the tree crowns in the study area were organized into specific object segments comprising similar digital pixel values (Figure 7), having defined edges that separate each group from other distinct groupings [17].  (Table 4), maximizing the differentiation in the imagery between features.
The ENVI edge segmentation procedure utilized region-based algorithms which specify spectral similarity (scale setting) and spatial similarity (merge setting). This segmentation procedure is typically less sensitive to slight variations in texture, which can be a significant advantage when using highresolution imagery [58]. In addition, the edge segmentation is overall a more accurate process for characterizing tree crowns [59,60]. Settings for scale level and merge level must be defined appropriately based on the image resolution and object complexity in the landscape [13]. User-defined thresholds for the scale level and merge level define the segments (with a standard values for both settings from 0 to 100), and the success of the segmentation is dependent on the complexity of the landscape features [58]. Incorrect settings of the scale level and merge level will lead to segments that are too complex, which subdivide known features, or assimilation of multiple features into large  (Table 4), maximizing the differentiation in the imagery between features.
The ENVI edge segmentation procedure utilized region-based algorithms which specify spectral similarity (scale setting) and spatial similarity (merge setting). This segmentation procedure is typically less sensitive to slight variations in texture, which can be a significant advantage when using high-resolution imagery [58]. In addition, the edge segmentation is overall a more accurate process for characterizing tree crowns [59,60]. Settings for scale level and merge level must be defined appropriately based on the image resolution and object complexity in the landscape [13]. User-defined thresholds for the scale level and merge level define the segments (with a standard values for both settings from 0 to 100), and the success of the segmentation is dependent on the complexity of the landscape features [58]. Incorrect settings of the scale level and merge level will lead to segments that are too complex, which subdivide known features, or assimilation of multiple features into large polygons [13]. The settings for scale and merge level are defined at the discretion of the investigator, as there is no perfectly objective approach or process [14]. Several tropical forest studies have defined a range between 10 and 40 for the scale level (creating more segments) and between 70 and 90 for the merge level (merging several segments) to ensure a good characterization of features within an image and optimize tree crown definition [9,18,26,27,61]. Through testing of several different combinations of scale and merge values, a setting of 11 for the scale level and 86 for the merge level ensured a good object differentiation in the Arboretum study area (Figure 7). The two SVIs described in this study were the data inputs used for the segmentation process, as each of the SVIs maximize differentiation between species in the study area ( Figure 6).
A rule-set object-based classification identified the tree species based only on their distinctive mean crown response from the imagery bands and SVIs from the segmentation [61]. This is in contrast to traditional pixel-based classifications, which classify on individual pixel values and provide a class for each pixel in the imagery. The rule-set process allows the user to define specific data inputs to each class, and each data set can be constrained to a specific range of values that represent a certain class [61]. In order to determine the ability of the imagery to classify tree species, no other values (surface texture, segment metrics, etc.) were included in the classification procedure. We assigned weights to each band or SVI based on their Wilks' Lambda score within the rule-set process, ensuring the application of each band and SVI was at a level commensurate with their ability to differentiate tree species.
For the rule-set classification process, a value range of +/-5% from each mean value for each WorldView-3 band and SVI defined a particular tree species. This value range chosen encompasses the variability that can exist within a particular species allowing for any inter-species variation while still maintaining a species-based consistency in mean spectral response within the defined object [25].
An accuracy analysis utilizing an error matrix compared the classification results to known point locations of the same canopy tree species within the Arboretum, determining if the classifier was identifying the existence and location of the tree type studied. The locations of individual tree species within the study area were acquired from the Arboretum inventory data file [28], described earlier in the Materials and Methods Section 2.2.2. The error matrix was applied to Arboretum trees with a DBH >50 cm to ensure that the tree was large enough to be a part of the visible Arboretum canopy [62]. Trees that fell in this category were checked in situ to verify they were a part of the top canopy of the Arboretum.

Results and Discussion
The first classification analysis performed (Classification 1) assessed all bands and SVIs available for the process, excluding the Coastal band and Blue band, as both of these bands are unreliable for use in a humid tropical environment as they are susceptible to severe atmospheric attenuation [46]. The data input for Classification 1 included only the top three image bands (Red Edge, Near-IR1 and Near-IR2) and the top SVI (WV-RESWI), based on the Wilk's Lambda scores and Correlation values. Table 6 shows the results of the error matrix of Classification 1. An "Unknown" classification result was specified to represent trees or ground cover not included in the error matrix analysis.
Classification 2 included additional bands and SVIs based on information from Tables 4 and 5, and the visual separation of values as shown in Figures 5 and 6. For each tree species, rule weights varied for each data input in the classification rule-set. These were based on the combination of bands or SVIs that best defined a particular species, improving the classification accuracy. The error matrix results of Classification 2 are shown in Table 7.
As is evident from the error matrix results, the WorldView-3 imagery captured the spectral variations between the tree species studied. The low errors of commission (7.89%) and omission (14.63%) in the error matrix for Classification 2 (Table 7) suggest that the segmentation parameters and rule-sets chosen maximized the unique species reflectivity response in the imagery. Figure 8 illustrates the distribution of species in and around the Arboretum as defined by Classification 2.  Miss-classifications could be the result of shadowing within the Arboretum image due to the off-nadir WorldView-3 image acquisition at 26.2 • zenith view angle (Table 1), as foreground tree crowns likely obscured or shadowed other possible canopy crowns. There is some evidence of this by the low Canopy Average values derived from the Arboretum (Figure 4). Other variables, such as lianas suspended within the canopies studied, could have also affected classification accuracy [23,24], potentially adding to miss-classifications and increasing the unknown fraction in the error matrix.
Both SVIs significantly improved the segmentation process and overall classification accuracy outperforming all other bands available ( Figure 6, Table 4, Table 5). The WV-ACRI showed promise, but its overall performance was likely affected in this study due to suppressed Canopy Average values ( Figure 4). The WV-RESWI performed well in showing variability between the chlorophyll production and plant structure because of its focus on the Red Edge variability between species [46].
Over-classifying individual crowns to a specific species was likely, as the rule-set process attempts to fit a tree type to a defined segment that represents all or part of a known tree crown within the study. The high classification accuracy from Classification 2 (85.37%) likely reflects this potential error, as the tree assemblage within the study area is extremely complex, as tree crowns can easily overlap creating confusion in the classification procedure. Figure 8, which shows the geographical distribution of Classification 2 illustrates this effect, where C. alliodora and C. odorata are shown as widespread as P. macroloba in the Arboretum based on the classification, but in truth P. macroloba is the dominant species in this area [28].
attempts to fit a tree type to a defined segment that represents all or part of a known tree crown within the study. The high classification accuracy from Classification 2 (85.37%) likely reflects this potential error, as the tree assemblage within the study area is extremely complex, as tree crowns can easily overlap creating confusion in the classification procedure. Figure 8, which shows the geographical distribution of Classification 2 illustrates this effect, where C. alliodora and C. odorata are shown as widespread as P. macroloba in the Arboretum based on the classification, but in truth P. macroloba is the dominant species in this area [28]. Figure 8. The output of Classification 2 using custom rule-sets for each species. Legend: C. elasticagreen, C. odorata-blue, C. alliodora-red, P. macroloba-yellow, P. sp. A-teal, S. microstachyummagenta.

Conclusions
We sought to establish a process of differentiating selected crown-canopy tree species within the tropical forest regime through the implementation of an object-based rule set classification schema. The complexity of a tropical forest assemblage, both in species diversity and inter-species variability, poses a great challenge to identifying individual tree species [2]. Successfully achieving a deeper understanding of the species assemblage within the tropical forest could improve our overall understanding of the role of forests in climate [1]. The results of this study show that a simple objectbased rule-set classification, using readily available multispectral data, can yield accurate results in a complex tropical rainforest.
Characterizing the intra-crown reflectivity response (using the collective pixels in a segment that defines the tree crown) is the key to identification of species in a complex forest [63], as all of the Figure 8. The output of Classification 2 using custom rule-sets for each species. Legend: C. elastica-green, C. odorata-blue, C. alliodora-red, P. macroloba-yellow, P. sp. A-teal, S. microstachyum-magenta.

Conclusions
We sought to establish a process of differentiating selected crown-canopy tree species within the tropical forest regime through the implementation of an object-based rule set classification schema. The complexity of a tropical forest assemblage, both in species diversity and inter-species variability, poses a great challenge to identifying individual tree species [2]. Successfully achieving a deeper understanding of the species assemblage within the tropical forest could improve our overall understanding of the role of forests in climate [1]. The results of this study show that a simple object-based rule-set classification, using readily available multispectral data, can yield accurate results in a complex tropical rainforest.
Characterizing the intra-crown reflectivity response (using the collective pixels in a segment that defines the tree crown) is the key to identification of species in a complex forest [63], as all of the information (and variability) within the tree crown is necessary information to identify tree-species within the forest assemblage. This provides a significant advantage as compared to pixel-based approaches, as the natural variability in the object (the tree canopy) is included in the characterization of the object [12]. This study identified that the necessary components needed to achieve that goal are high spatial and radiometric resolution imagery (e.g., WorldView-3 or similar) with appropriate bands and/or data products [20], careful corrections for illumination and attenuation-absorption effects [25], and a segmentation process that properly assigns single-tree canopies to one (or a few) segmented pixel cluster(s). Continuing to improve the inputs to the object-based rule-set classification parameters (more precise segmentation, better data inputs through field measurements of more tree species) will be important steps in future analyses. This will ultimately contribute to accurately defining more tree species with their unique spectral signatures in other diverse tropical rainforest locations.
It is envisioned that this simple, straightforward process can be expanded to more regional scale tropical environments.