eBird Status and Trends divides space and time into variably sized "stixels" within which individual base models are fit. The process of stixelization is performed many times and the prediction at any given point is the median of the predictions from all the stixels that that point falls in. This function loads summary statistics for each stixel, for example, the size of the stixels, the number of observations within each stixel, and a suite of predictive performance metrics (PPMs) for the model fit within that stixel.

load_stixels(path, ext, return_sf = FALSE)



character; directory that the Status and Trends data for a given species was downloaded to. This path is returned by ebirdst_download() or get_species_path().


ebirdst_extent object; the spatiotemporal extent to filter the data to. The spatial component of the extent object must be provided in unprojected, latitude-longitude coordinates.


logical; whether to return an sf object of spatial points rather then the default data frame.


Data frame, or sf object if return_sf = TRUE, containing stixel summary data. Data are organized with one stixel per row and each stixel identified by a unique stixel_id, the centroid of each stixel in space and time is specified by lat, lon, and date.


if (FALSE) {
# download example data
path <- ebirdst_download("example_data", tifs_only = FALSE)
# or get the path if you already have the data downloaded
path <- get_species_path("example_data")

# load stixel summary information
stixels <- load_stixels(path)