ebird-occurrences.R

library(arrow)
library(dplyr)
library(terra)

path <- arrow::s3_bucket("ebird/observations", endpoint_override = "minio.carlboettiger.info")
obs <- arrow::open_dataset(path) |> to_duckdb()
df <- obs |> mutate(latitude = round(latitude, digits=1),
                   longitude = round(longitude, digits=1)) |> 
  count(longitude, latitude) |> collect() |> 
  mutate(n = log(n))

r <- rast(df, crs="epsg:4326")
y <- r; y[] <- 0; r <- cover(r,y) # NA fill
plot(r, col= viridisLite::viridis(100) )