A multi-sensor benchmark dataset for detecting individual trees in airborne RGB, Hyperspectral and LIDAR point clouds
Maintainer: Ben Weinstein - University of Florida.
This benchmark dataset is the first dataset to have consistent annotation approach across a variety of ecosystems. This repo is the R package for reproducible evaluation against the dataset. The benchmark dataset contains over 6,000 image-annotated crowns, 424 field-annotated crowns, and 3,777 overstory stem points from a wide range of forest types.
library(devtools)
install_github("Weecology/NeonTreeEvaluation_package")
To download evaluation data from the Zenodo archive (1GB), use the download() function to place the data in the correct package location. Download the much larger training data, set training=TRUE.
library(NeonTreeEvaluation)
download()
The package contains two vignettes. The ‘Data’ vignette describes each datatype and how to interact with it in R. The ‘Evaluation’ vignette shows how to submit predictions to the benchmark.
The format of the submission is as follows
- A csv file
- 5 columns: Plot Name, xmin, ymin, xmax, ymax
Each row contains information for one predicted bounding box.
The plot column should be named the same as the files in the dataset without extension (e.g. SJER_021 not SJER_021.tif) and not the full path to the file on disk. Not all evaluation data are available for all plots. Functions like evaluate_field_crowns and evaluate_image_crowns will look for matching plot name and ignore other plots.Depending on the speed of the algorithm, the simplest thing to do is predict all images in the RGB folder (see list_rgb()) and the package will handle matching images with the correct data to the correct evaluation procedure.
The package contains a sample submission file.
library(raster)
library(dplyr)
library(NeonTreeEvaluation)
head(submission)
#> xmin ymin xmax ymax score label plot_name
#> 1 217.24730 265.0254 290.4331 337.5056 0.7930149 Tree ABBY_020
#> 2 284.37094 339.3117 338.5593 388.7977 0.7155186 Tree ABBY_020
#> 3 197.60432 350.3243 249.3276 400.0000 0.7128302 Tree ABBY_020
#> 4 25.21722 186.6763 123.3826 283.8552 0.6598154 Tree ABBY_020
#> 5 332.43198 293.3734 377.3869 343.7137 0.6279798 Tree ABBY_020
#> 6 198.98254 181.8541 266.0995 249.9332 0.5793932 Tree ABBY_020
Author | Precision | Recall | Cite/Code |
---|---|---|---|
Weinstein et al. 2020 | 0.66 | 0.79 | https://deepforest.readthedocs.io/ |
Silva et al. 2016 | 0.34 | 0.47 | lidR package |
The main data source are image-annotated crowns, in which a single observer annotated visible trees in 200 40m x 40m images from across the United States. This submission has bounding boxes in image coordinates. To get the benchmark score image-annotated ground truth data.
#Get a three sample plots to run quickly, ignore to run the entire dataset
df<-submission %>% filter(plot_name %in% c("SJER_052"))
#Compute total recall and precision for the overlap data
results<-evaluate_image_crowns(submission = df,project = T, show=F, summarize = T)
#> [1] SJER_052
#> 1292 Levels: 2018_SJER_3_252000_4104000_image_628 ...
results[1:3]
#> $overall
#> # A tibble: 1 x 2
#> precision recall
#> <dbl> <dbl>
#> 1 1 0.778
#>
#> $by_site
#> # A tibble: 1 x 3
#> # Groups: Site [1]
#> Site recall precision
#> <chr> <dbl> <dbl>
#> 1 SJER 0.778 1
#>
#> $plot_level
#> # A tibble: 1 x 3
#> # Groups: plot_name [1]
#> plot_name submission ground_truth
#> <fct> <int> <int>
#> 1 SJER_052 7 9
For a list of NEON site abbreviations: https://www.neonscience.org/field-sites/field-sites-map
Author | Recall | Cite/Code |
---|---|---|
Weinstein et al. 2020 | 0.61 | https://deepforest.readthedocs.io/ |
The second data source is a small number of field-deliniated crowns from three geographic sites. These crowns were drawn on a tablet while physically standing in the field, thereby reducing the uncertainty in crown segmentation.
df <- submission %>% filter(plot_name=="OSBS_95_competition")
results<-evaluate_field_crowns(submission = df,project = T)
#> [1] OSBS_95_competition
#> 1292 Levels: 2018_SJER_3_252000_4104000_image_628 ...
results[1:3]
#> $overall
#> # A tibble: 1 x 2
#> precision recall
#> <dbl> <dbl>
#> 1 0.029 1
#>
#> $by_site
#> # A tibble: 1 x 3
#> # Groups: Site [1]
#> Site recall precision
#> <chr> <dbl> <dbl>
#> 1 OSBS_95 1 0.029
#>
#> $plot_level
#> # A tibble: 1 x 3
#> # Groups: plot_name [1]
#> plot_name submission ground_truth
#> <fct> <int> <int>
#> 1 OSBS_95_competition 34 1
Author | Recall | Cite/Code |
---|---|---|
Weinstein et al. 2020 | 0.74 | https://deepforest.readthedocs.io/ |
The third data source is the NEON Woody Vegetation Structure Dataset. Each tree stem is represented by a single point. This data has been filtered to represent overstory trees visible in the remote sensing imagery.
df <- submission %>% filter(plot_name=="JERC_049")
results<-evaluate_field_stems(submission = df,project = F, show=T, summarize = T)
#> [1] "JERC_049"
results
#> $overall
#> recall
#> 1 0.5555556
#>
#> $by_site
#> # A tibble: 1 x 2
#> Site recall
#> <fct> <dbl>
#> 1 JERC 0.556
#>
#> $plot_level
#> siteID plot_name recall n
#> 1 JERC JERC_049 0.5555556 9
If you would prefer not to clone this repo, a static version of the benchmark is here: https://zenodo.org/record/3723357#.XqT_HlNKjOQ
library(raster)
library(NeonTreeEvaluation)
#Read RGB image as projected raster
rgb_path<-get_data(plot_name = "SJER_021",type="rgb")
rgb<-stack(rgb_path)
#Find path and parse
xmls<-get_data("SJER_021",type="annotations")
annotations<-xml_parse(xmls)
#View one plot's annotations as polygons, project into UTM
#copy project utm zone (epsg), xml has no native projection metadata
xml_polygons <- boxes_to_spatial_polygons(annotations,rgb)
plotRGB(rgb)
plot(xml_polygons,add=T)
To access the draped lidar hand annotations, use the “label” column. Each tree has a unique integer.
library(lidR)
path<-get_data("TEAK_052",type="lidar")
r<-readLAS(path)
trees<-lasfilter(r,!label==0)
plot(trees,color="label")
We elected to keep all points, regardless of whether they correspond to tree annotation. Non-tree points have value 0. We highly recommend removing these points before predicting the point cloud. Since the annotations were made in the RGB and then draped on to the point cloud, there will naturally be some erroneous points at the borders of trees.
Hyperspectral surface reflectance (NEON ID: DP1.30006.001) is a 426 band raster covering visible and near infrared spectrum.
path<-get_data("MLBS_071",type="hyperspectral")
g<-stack(path)
nlayers(g)
#> [1] 426
#Grab a three band combination to view as false color
f<-g[[c(52,88,117)]]
plotRGB(f,stretch="lin")
To add score to this benchmark, please submit a pull request to this README with the scores and the submission csv for confirmation.
This benchmark is currently in review. Either cite this repo, or the original article using these data: 1 Weinstein, Ben G., et al. “Individual tree-crown detection in RGB imagery using semi-supervised deep learning neural networks.” Remote Sensing 11.11 (2019): 1309. https://www.mdpi.com/2072-4292/11/11/1309