An Approach for High-Resolution Mapping of Hawaiian Metrosideros Forest Mortality Using Laser-Guided Imaging Spectroscopy

Rapid ‘Ōhi‘a Death (ROD) is a disease aggressively killing large numbers of Metrosideros polymorpha (‘ōhi‘a), a native keystone tree species on Hawaii Island. This loss threatens to deeply alter the biological make-up of this unique island ecosystem. Spatially explicit information about the present and past advancement of the disease is essential for its containment; yet, currently such data are severely lacking. To this end, we used the Carnegie Airborne Observatory to collect Laser-Guided Imaging Spectroscopy data and high-resolution digital imagery across >500,000 ha of Hawaii Island in June–July 2017. We then developed a method to map individual tree crowns matching the symptoms of both active (brown; desiccated ‘ōhi‘a crowns) and past (leafless tree crowns) ROD infection using an ensemble of two distinct machine learning approaches. Employing a very conservative classification scheme for minimizing false-positives, model sensitivity rates were 86.9 and 82.5, and precision rates were 97.4 and 95.3 for browning and leafless crowns, respectively. Across the island of Hawaii, we found 43,134 individual crowns suspected of exhibiting the active (browning) stage of ROD infection. Hotspots of potential ROD infection are apparent in the maps. The peninsula on the eastern side of Hawaii known as the Puna district, where the ROD outbreak likely originated, contained a particularly high density of brown crown detections. In comparison, leafless crown detections were much more numerous (547,666 detected leafless crowns in total) and more dispersed across the island. Mapped hotspots of likely ROD incidence across the island will enable scientists, administrators, and land managers to better understand both where and how ROD spreads and how to apply limited resources to limiting this spread.


Airborne Data Collection
We collected data with the CAO between 7 June and 5 July 2017 over the island of Hawaii with three co-aligned instruments: a visible-to-shortwave imaging spectrometer (VSWIR), a dual-channel airborne light-detection and ranging (LiDAR) scanner, and a high-resolution 3-band (red, green, and blue) 60-megapixel digital camera. For optimal fight coverage, the island was partitioned into 15 compartments of fairly uniform forest type ranging in size from 12,000 to 75,000 ha. During fight, aircraft speed was maintained within 10% of the 60 m s −1 nominal air speed and within 200 m of the nominal planned elevation of 2000 m above ground level. The LiDAR was operated at an effective pulse frequency of 200 kHz, a scan frequency of 34 Hz, 30 percent overlap, and total feld-of-view of 38 degrees to match the 34-degree FOV of the VSWIR after clipping 2 degrees from the scan edges. These settings yielded a mean density of 3.3 pulses m −2 and a VSWIR resolution of 2 m. VSWIR data were collected in 427 channels between 350 and 2485 nm at 5 nm increments (full-width at half max). The high-resolution digital camera had an effective resolution of 24 cm.

Data Processing and Orthorectifcation
Estimated location error in LiDAR point vertical positioning from planar surface analysis was less than 10 cm (root mean squared error). LiDAR point cloud data were processed to 1-m resolution digital surface maps of ground elevation, canopy surface elevation, and canopy surface height above ground using the LAStools software suite (Rapidlasso GmbH; Gilching, Germany).
The VSWIR radiance data were frst orthorectifed using LiDAR surface maps, positon, and orientation data from the on-board GPS-IMU navigation system, and known boresight and position offsets between the two instruments. To improve the signal-to-noise ratio, the radiance data were reduced by averaging each pair neighboring bands. This yielded a dataset of 214 bands at increments of 10 nm. Using computed observation angles and elevation, atmospheric correction was performed using ACORN-6LX software (Imspec LLC; Palmdale, CA, USA) to transform the VSWIR radiance data into surface-level refection spectra. For each fight line, we used an iterative procedure to fnd the optimal setting for the ACORN-6LX visibility parameter that is used to specify the amount of aerosols present at the time of collection. Beginning with a setting of 100 km, the visibility parameter was decreased until the resulting refectance at 420 nm was <1%. The individual VSWIR refectance data were mosaicked into one image for each of the 15 fight coverage compartments using the criteria of minimum phase angle (angle from sun to target to instrument) for overlap regions. These coverages were mosaicked into a full island map using the criteria of maximum normalized difference vegetation index (NDVI).

Training Data
Using orthorectifed color images from the on-board digital camera as a reference, we identifed 688 VSWIR pixels within crowns exhibiting the brown state of ROD infection, labeled "Brown", and 903 VSWIR pixels within leafess crowns, labeled "Leafess", across the entirety of 'ōhi'a stand distribution on Hawaii Island (Figure 1). Because signal from background vegetation can have a large impact on observed spectra, especially in leafess crowns, care was taken to select training pixels from central regions of crowns and that tree branching was dense enough to ensure that the resulting spectra would be dominated by the class of interest. We also identifed 13,946 pixels representing other entities we were not interested in mapping; these were primarily living canopies of other plant species, along with a small number of non-vegetation targets such as bare ground, manmade structures, and other surfaces, all of which provided the model with a third target class, labeled "Live". The non-vegetated pixels make up <1% of the target set; therefore, despite their spectral characteristics being different from each other and living vegetation, we chose to include them in the living category, because we were not interested in separating samples classifed into this category. A few shaded pixels were included for each of the three classes to reduce any artifact caused by decreased signal-to-noise ratios in darker pixels. The full VSWIR refectance spectrum was extracted for each of these training pixels and labeled by class ( Figure 2). Refectance values of wavelengths known to be absorbed or scattered by atmospheric constituents were removed, leaving 159 bands representing four regions of the measured spectrum (Table 1). Cleaned individual spectra were subsequently brightness-normalized (BN), such that the magnitude (Euclidean norm) of each BN spectrum was 1.0. This was done for each spectrum (pixel) using the formula in which v r is the original refectance spectrum as a vector, v n is the brightness-normalized refectance spectrum.

Modelling
Both false detections and false omissions have the potential to waste resources when fnal maps are used in the feld; we sought to maximize the model results under dual criteria of high sensitivity and high precision. To meet this goal, we combined two model forms into an ensemble model for differentiating pixels within the Brown and Leafess classes from those in the Live class. The frst model form was a Support Vector Machine (SVM), which fnds optimal divisions between the data points in n-dimensional vector space [14]. This machine learning tool has been demonstrated to be an effective classifcation tool with hyperspectral data [15,16] and has performed well with CAO data in particular [17][18][19]. SVMs naturally incorporate information from all wavelengths of each input spectrum and offer the beneft of being computationally simple to train and apply. The second model form used was the Gradient Boosting Machine (GBM) [20,21]. GBM models are created by generating an ensemble of fairly weak (shallow) decision trees. Because each successive tree is ft to the remaining residuals from the previous combination of trees, a process known as 'boosting', GBMs incorporate both bias and variance reduction, distinguishing them from other decision-tree based approaches such as a Random Forest (RF), which are focused on variance reduction. GBM models are able to incorporate nonlinear interactions between spectral features while also being robust to spectral collinearity. We found GBM models outperformed RF models during exploratory testing in this study. We found optimal model parameters for both the SVM and GBM fts using a multi-stage grid search technique. The frst stage grid covered a broadest region of the available parameter space with the intent of surrounding the most-optimal parameters. The search grid for each successive stage used fner granularity, covering a small region of the grid that proved most optimal in the last stage. A 10-fold cross-validation was performed for each set of parameters in the grid to get a reasonable estimate of how well the model would perform on new data. All model training was done using the package scikit-learn [22] for Python. As the same tree canopy can simultaneously contain both brown foliage and leafess branches, leaving the canopy state rather ambiguous, we focused less on misclassifcations between these two states and more on misclassifed living trees. In addition, errors of commission regarding living trees would be much more problematic to users of the fnal maps. Thus, the criteria used during model parameter selection for both models was the number of correctly classifed Brown and Leafess pixels (combined) while holding the false positive rate of each of the Brown and Leafess classes to less than 1% (independently). Error rates were assessed on a per-pixel, rather than a per-crown, basis because of the relatively large number of small, single-pixel instances in the training data. Once optimal parameter values for each model were obtained, we trained fnal models with the full set of training pixels using these parameters to build models for subsequent full-island map application.
We next devised a way to combine the two models in order to improve the applied model predictions. Both models, when applied at the whole-island scale, left patches of false-positives for either the Brown or Leafess classes due to heavily site-specifc vegetation or observation conditions. We sought to eliminate these artifacts, while preserving model accuracy. Initial testing suggested that the SVM model was better at detecting Brown and Leafess pixels yet more prone to false positives. We found that the best performance with the least amount of artifact occurred when the GBM model was able to 'veto' perceived spurious predictions of the SVM. For each pixel and for each class, each of the two models produced a likelihood value between 0 and 1, which effectively portrays the model's confdence that a given pixel belongs to a given class. We found that the vast majority of false positives in maps produced by the SVM model had only moderate likelihood values, yet corresponding likelihood values from the GBM model were extremely low. Using this feature to our advantage, any SVM likelihood value for Brown or Leafess classes of less than 0.975 was replaced by the associated GBM likelihood, if the GBM likelihood value was less than 0.01. For a given class i, this model combination rule can be written as in which L i is the combined model likelihood value for class i ∈ {Brown, Leafess}, L i,svm is the SVM model likelihood value for class i, and L i,gbm is the GBM model likelihood value for class i. For the likelihood values produced by the above combination approach, still in the range 0 to 1, we used the following decision rule to assign a class to each pixel: in which L V , L B , and L D are likelihood values for classes Live, Brown, and Leafess, respectively, T is a specifed likelihood threshold value, and C is the predicted classifcation of a given pixel.

Model Application
We used a combination of two metrics to evaluate potential values of the likelihood threshold value, T, for application of the combined model to full maps of Hawaii Island. For a given value of T, we computed classifcation sensitivity (producer accuracy) and classifcation precision (user accuracy) for the Brown and Leafess classes using the SVM, GBM, and combined model predictions from the 10-fold cross-validation training under the optimal parameters for each model. We repeated this calculation for each of 19 values, T ∈ {0.05, 0.10, . . . , 0.95}, to produce response curves of both metrics for both Brown and Leafess classes.
Ultimately, we selected threshold of 0.65, with which we applied the combined model decision rule C(T, L V , L B , L D ) to the model prediction likelihood values of the full CAO VSWIR coverage of the island. Because VSWIR pixels are 4 m 2 in area, light refecting from brown vegetation on the ground can blend with light refecting from the edges of tree vegetation, causing the spectra for a number of edge pixels to be identifed as Brown or Leafess. We fltered out false detections due to such edges by using the canopy height data built from the LiDAR data. If the LiDAR-derived height of the top of the vegetation surface was <1.5 m in any of the four 1 m × 1 m quadrant subregions within a VSWIR pixel, then the pixel was re-classifed as Live when building output maps.

Additional Spectral Filtering
Field verifcation of early versions of the fnal map indicated that excessive false-positives remained for the Brown class. There were three scenarios in which a crown in the map was not observed to be a leafess '¯Each of these scenarios below ohi'a in the feld later in the same year. involved vegetation not represented in the training data sets.

1.
A small number of instances of crowns in the far north region, Kohala, were classifed as Brown but observed to be green three months after the CAO fight coverage. The most likely cause of this was that these were '¯ohi'a fowers are typically ohi'a exhibiting unusually heavy fowering (' b right red) along with relatively low foliage. This situation may have been exacerbated by an extremely low sun angle at the time of overfight. 2.
In the region of Hawaii Volcanoes National Parks, some instances of tall or climbing brown Dicranopteris linearis fern (uluhe) were found to be misclassifed as Brown by the model.

3.
Also in the region of Hawaii Volcanoes National Park, a large number Morella faya (fretree) trees are known to be dying in the area. These trees also appear dry and brown-similar, but not identical, to ' ōhi'a in the brown-state of ROD infection-and the model classifed map pixels, covering them as Brown with high likelihood values.
We created a flter to reduce detected crowns under these false-positive scenarios using an analysis of the spectra of the affected pixels. The pixels from brown D. linearis fern and fowering 'ōhi'a in scenarios 1 and 2 all had highly-negative slopes in the red-green region of the visible spectrum (600-680 nm), whereas brown ' ōhi'a and M. faya both consistently had positive slopes in the same spectral region. We therefore used the following criteria to identify and remove such pixels from the map.
(r 650 + r 660 ) − (r 600 + r 610 ) < −0.004 (2) in which r w is the refectance at the indicated wavelength w in nm. M. faya pixels classifed as Brown (scenario 3 above) were separable from verifed Brown 'ōhi'a pixels using a combination of two features: (1) the peak height difference between the infrared and frst shortwave regions of the spectrum was greater for 'ōhi'a pixels and (2) the slope of the peak in the frst shortwave region (1640-1680 nm). We employed the following criteria to identify and remove brown M. faya pixels from the map.
While we do not believe any examples of brown D. linearis fern, brown M. faya, or fowering 'ōhi'a were represented in the training data, some of the Brown training pixels were caught by one or more of the two spectral flter criteria. For the purposes of evaluating map accuracy, it was unclear how to represent these ambiguous pixels in a confusion matrix. Applying the flter criteria to just the predicted classes would show degraded mapping performance over the pre-flter results, while applying the flter criteria to both observed and predicted classes would show improvement. We therefore simply kept track of the number of pixels affected, and performed an analysis of the potential effects of this fltering on model sensitivity of Brown pixels.
We made fnal maps of the locations of suspected ROD-infected brown 'ōhi'a tree crowns (as clusters of predicted Brown pixels) and leafess tree crowns (as clusters of predicted Leafess pixels) from the remaining pixel-level detections. Centroid points of all contiguous, adjacent (no diagonal connections) clusters of Brown or Leafess pixels greater than or equal to two pixels (8 m 2 ) in size were used as spatial coordinates for suspected brown-state 'ōhi'a or leafess tree crowns. Although only a small number of points were found in areas known to not contain 'ōhi'a, these points were considered false detections and were therefore masked out in the fnal maps.

Results
We found both the SVM and GBM models to be capable of detecting a high number of Brown and Leafess training pixels with minimal false detection. Use of a radial basis function (RBF) for the SVM model yielded the best results. Under this kernel, the optimal parameter values were a gamma of 0.0075 and cost of 200. For the GBM approach, optimal parameter values were a learning rate of 0.225, a maximum tree depth of 6, and the assembly of 1750 trees.
Cross validation results indicated that the SVM was more sensitive to spectral features that defne both Brown and Leafess classes ( Table 2). Sensitivity values were 90.1% and 88.0% for Brown and Leafess classes, respectively. The sensitivity rate for both these classes combined, compared to the Live class, was 90.5%. GBM was considerably more conservative during training, with many more false-negatives for the Brown class, though it seemed to be slightly better at distinguishing Brown from Leafess pixels. The GBM Brown sensitivity was 70.8%, Leafess sensitivity was 65.1%, and combined Brown and Leafess sensitivity was 69.5%. Upon application, the SVM was much more prone to erroneous detections of both Brown and Leafess pixels due to lighting or forest canopy conditions not represented in the training dataset. When analyzing the effect of threshold value on the different models, we found a generally slow decline in true detection rates for all models as the likelihood threshold T increased from 0 to approximately 0.8 ( Figure 3). Thus, selection of any value within this range will still result in a high number of true positives. Additionally, the rate of decline was signifcantly lower for the GBM model compared to the SVM model. This analysis showed that the combined model retained nearly all of the superior sensitivity of the SVM model, indicating that for most of the truly Brown or truly Leafess training pixels the GBM was in at least partial agreement. Similarly, precision rapidly increased with T to >90% for all models in the range 0 to ca. 0.4 ( Figure 4). However, compared to the values for sensitivity, the difference between models was not as large. The rate of precision for Brown pixels was lower for the GBM than for the SVM. However, with respect to Leafess pixels, the two models had nearly the same values for precision. We found the combined model not only kept the properties of the best of its component models, it improved upon them. Here, the combined model Brown pixel precision values most rapidly improve with increasing threshold, reaching a semi-plateau by a threshold value of 0.2. To a lesser degree the same was true for Leafess pixels. Given the above results, and upon inspection of several sample map areas, we selected a likelihood threshold of 0.65. When applied to the combined model, this threshold maintained desirable detection rates (i.e., sensitivity values of 86.9% and 82.5% for Brown and Leafess pixels, respectively), while greatly reducing the false positives (i.e., precision values >95% for both classes) to acceptable levels ( Table 3). The kappa statistic for this classifcation was 0.895. Only 42 of the 13,946 Live training samples were misclassifed as class Brown or Leafess, indicating the model is more likely to classify ambiguous trees as Live. The spectral fltering to remove D. linearis fern, M. faya tree, and fowering 'ōhi'a affected a noticeable number of Brown training pixels. We found that 48 of the 688 Brown training pixels met the criteria for one or both of the two spectral flters (Equations (1) and (2)). Thus, assuming our training pixels truly contain no examples of D. linearis fern, M. faya tree, or fowering 'ōhi'a, the model may lose up to 7.0% of Brown pixel sensitivity by applying these flters. However, 19 of these 48 pixels were not even identifed as Brown by the combined model before fltering, indicating the actual loss in sensitivity might be closer to 4.2%. Of the 8316 Live pixels that met one or more of the flter criteria, none were mistakenly classifed as Brown. This indicates that, outside the region of Hawaii Volcanoes National Park, this particular type of misclassifcation is probably rare.
The LiDAR-derived canopy height and the spectral fltering both reduced the number of isolated Brown pixel detections, resulting in a statistically robust map with very little visual artifact present ( Figure 5). Across the island of Hawaii, we found 43,134 individual crowns suspected of exhibiting the brown state of ROD infection. Hotspots of potential ROD infection are apparent in the maps. The peninsula on the eastern side of Hawaii known as the Puna district, where the ROD outbreak likely originated, contained a very high density of brown-state crown detections. In comparison, leafess crown detections were much more numerous (547,666 detected leafess crowns in total) and more dispersed across the island (Figure 6). It was apparent in the map ( Figure 6) that detection of leafess crowns was more sensitive to cross-fightline observation angle effects, leading to some linear artifact in the maps. However, we saw no evidence of an increase in false-positive rates of leafess crowns due to these viewing angle effect, but rather we observed that a decrease in detection on one side of each fightline was the cause of the artifacts. As viewing angle and solar angles approach nadir, green undergrowth beneath a thin, leafess crown has a more dominant effect on refectance, leading to this decreased detection. Since we were unable to distinguish the cause of death for so many trees, we cannot ascertain what percentage of leafess crowns were the result of ROD. However, dense areas of detected leafess crowns can be found in regions containing higher counts of brown-state crown detections, and ROD is expected to be the dominant cause of death in these locations.

Discussion
Using CAO laser-guided imaging spectroscopy, we obtained very high classifcation accuracy for detecting brown-state, potential ROD-infected, and leafess trees with both the GBM and SVM models and their combination. This suggests that the spectra of these two classes were highly distinguishable from surrounding vegetation. While it was possible to visually identify particular brown and leafess crowns with certainty in the three-band (red, green, and blue) 24 cm aerial imagery due to the location and pattern of occurrence, this was not the case for all images captured across the island. Due to differences in lighting conditions and aircraft movement, our RGB imagery was insuffcient to determine whether some trees were truly exhibiting the brown state of ROD infection. This is why hyperspectral data are often sought for diffcult classifcation problems [7,23,24]. As their canopies turn brown, ROD-infected 'ōhi'a trees also experience a characteristic change in leaf chemical composition and shape, both of which highly affect the refection of light across the full spectrum measurable by imaging spectroscopy [25] (Figure 2). The ensemble approach of mixing two model forms greatly improved whole-island brown-state and leafess detection maps. SVM models have been used extensively in remote sensing classifcation applications [16], and GBM models have more recently shown promise in the feld [18,26,27]. Yet, individually, each of the two model forms was prone to misclassifcation when applied to regions of the island not represented in the training data. An increase in training data would likely help address this issue, but generating a training set that completely spans the extremely diverse climate and biogeographical conditions represented on the island of Hawaii [28,29] presents a major challenge. Fortunately, each model form was overactive in different parts of the island and under different conditions. We were able to take advantage of this effect by combining the models into an ensemble. Ensemble models in general have been shown to improve performance over individual component models [30][31][32]. In fact, a GBM itself is a form of ensemble modeling, as classifcation trees are sequentially added during boosting to incrementally improve ft [33]. Additionally, more diversity in the component models of an ensemble has been shown to improve ensemble model performance and further reduce overftting [34]. Here, we mixed two distinct model forms and were able to reduce overftting (manifested most-critically as artifacts in the fnal maps) while retaining high target detection capabilities.
The range of spatial scale of this project was extreme in the sense of applied geospatial classifcation models. VSWIR data collected by the CAO in 2017 covered an area of approximately 525,000 ha, about 320,000 ha of which is considered forested. This suggests that our maps cover nearly 800 million pixels potentially able to contain infected 'ōhi'a. There are few examples of classifcation maps covering such a large area at high spatial resolution. This is not surprising given the challenges encountered in such a project. In addition to the challenging problem of minimizing fight coverage costs subject to unpredictable tropical weather conditions, merging data from numerous fights occurring under diverse lighting conditions is quite diffcult. The boundary between two different collections can easily result in visible artifacts in maps [35,36]. Neighboring areas may be covered by fight lines collected under extremely different conditions. Additionally, effects of the bi-directional refectance distribution function (BRDF) may signifcantly alter the shape of recorded spectra [37] and affect model predictions, although the "pushbroom" style of collection of the VSWIR instrument essentially keeps this a cross-track problem. In this project we chose to handle these challenges with a custom fight line mosaic scheme, manual cloud masking, and most importantly, training data identifcation across regions of diverse lighting and view angle conditions. While not entirely free of viewing angle artifact (see Figures 5 and 6), the resulting maps of the crown classes of interest, with surprisingly little artifact, across nearly a full month of data collection.
We chose a pixel-based approach for mapping brown and leafess trees, primarily because our target of interest was partial crowns. ROD symptoms typically spread from whole branch to whole branch [8]; as such, we believed a focus on pixel-sized was essential to keep model sensitivity high and consistent across the landscape. Working at the crown level was an alternative approach to this problem, and there are several tools available for the automatic delineation of individual tree crowns from LiDAR data [38][39][40], but had we incorporated this methodology, there would have been several setbacks. First, these tools are most-often used for small regions of fairly simple, consistent vegetation growthform (e.g., boreal forest). In contrast, the variation in growthform among living 'ōhi'a across Hawaii Island is incredibly high. Adding more variation from dead or dying trees would only increase this complexity. In addition, there are many instances of multiple tree crowns blending together into a single uniform canopy, which severely limits the effcacy of individual crown segmentation [38,39].
Our resulting maps detailing the extent and distribution of both brown-state crowns and leafess crowns across broad regions of Hawaii Island reveal a startling extent of mortality. The extent and distribution of brown-state canopies presented here agree well with Hawaii Island-wide helicopter surveys conducted in January 2017 that estimated that ca. 30,000 ha of ' ōhi'a forest stands had recently experienced levels of Ceratocystis-induced mortality >10% (unpublished data, Hawaii Department of Land and Natural Resources). Although approximately half a million leafess crowns documented across Hawaii Island cannot and should not be uniformly attributed to Ceratocystis-induced Rapid 'Ōhi'a Death (ROD), it is likely that the majority of the 43,134 trees exhibiting brown canopies in our 2017 imagery were experiencing the lethal effects of such Ceratocystis infection. Furthermore, these brown-state canopies constitute relatively recent mortality, certainly within the last year, and likely within 6 months prior to image collection. We employed numerous measures to reduce false-positive detection rates, such as a high likelihood threshold of 0.65 for class inclusion and minimum cluster size of 8 m 2 . These measures likely affect small trees more than large trees. However, this may not be a large problem, because prior results suggest that medium-sized trees are the most affected by ROD [8]. In addition, instances of isolated small brown-state trees should be minimal, because the disease spreads so quickly [9]. If small individual crowns were not detected, we feel it is very likely that many infected crowns in the surrounding area were detected, and fnal maps should accurately refect the relative density of ROD infection.
A prior, well-acknowledged, and well-studied process of 'ohi'a ¯mortality characterized as "'ohi'a ¯dieback" was responsible for high levels of mortality across an estimated 50,000 ha of windward portions of Hawaii Island during the thirty years between the mid-1950's to the mid-1980's [41]. By comparison, Ceratocystis-induced ROD has affected more than half that amount of area (i.e., 30,000 ha) in approximately one-tenth the time. Our brown-state tree crown maps clearly indicate that mortality from the disease has spread well beyond the area of infection that in 2014 was restricted to the easternmost Puna District of Hawaii Island; currently, 'ohi'a ¯stands experiencing Ceratocystis-induced mortality are found on each of the fve main volcanoes of the Island. Preliminary results from forest inventory plot data indicate that, on average, 'ohi'a stands ¯in which infected trees are present exhibit annual mortality rates of 10%, with some stands experiencing rates as high as 42% [42].
The fnal brown-state and leafess crown maps will support efforts to combat ROD throughout Hawaii Island, and on other islands if the disease is transported westward. Final maps show several hotspots where either the extent or the existence of infection was previously unknown. Recent developments have allowed the cause of death of these crowns to be rapidly confrmed as ROD using feld kits that can genetically screen for Ceratocystis fungus [43]. The location of confrmed hotspots of ROD infections can help pathologists, ecologists, and land managers to better understand mechanisms of how ROD spreads, what regions are most vulnerable, and where and how to best allocate resources to reduce damage [8]. Additionally, research is underway to determine the feasibility of restoring 'ōhi'a to infected areas using genetic resistance techniques [44]. These maps will help to identify examples of ' ōhi'a trees that are that may be resistant to the Ceratocystis species causing the disease. As a single snapshot in time, these maps integrate the past progression of ROD across the ' ōhi'a forests of Hawaii Island. However, remapping the island using laser-guided imaging spectroscopy will enable a more-refned and up-to-date understanding of the rate of ROD proliferation. While the general methodology presented here can be readily transferred to such a remapping or any new applications, several specifc facets of this methodology should be transferred with care. Just as highly-optimized models should not be blindly applied to new data-whether they come from different instruments, different time periods, or collected over different locations-the transferability of specifc values in the classifcation rule or in the different post-processing flters will need to be reevaluated during each application as well.

Conclusions
Our results demonstrate that it is feasible to produce detailed, high-resolution maps of the likely extent of ROD outbreak at the level of individual tree crowns throughout native Metrosideros polymorpha ('ōhi'a) forests of Hawaii Island using airborne laser-guided imaging spectroscopy. By identifying [10] and compiling LGIS refectance spectra of individual ROD-infected brown-state tree crowns and leafess crowns identifed in the high-resolution color imagery, we found that these two tree conditions exhibit high spectral separability both from living tree canopies, as well as from one another. Via an ensemble of two different machine learning algorithms, we accurately identifed locations of potential ROD infection, both current (brown-state) and past (leafess), with minimal false detection. Such maps open a door to new research to aid in the effort to reduce the impact of Rapid ' Ō hi'a Death and to potentially restore 'ōhi'a forests already lost to this deadly pathogen.