Example usage of iNaturalist API to create leaflet map for a taxon within radius of a point

create_inat_leaflet_map.R

get_inat_obs_radius <- function(taxon_id, lat, lng, radius) {
  api_response <- httr2::request(
    "https://api.inaturalist.org/v1/observations/"
  ) %>%
    httr2::req_url_query(
      taxon_id = taxon_id,
      lat = lat,
      lng = lng,
      radius = radius
    ) %>%
    httr2::req_retry() %>%
    httr2::req_perform() %>%
    httr2::resp_body_json(simplifyVector = TRUE)

  tibble::as_tibble(api_response$results)
}

get_inat_obs_radius(
  taxon_id=157924,
  lat=50.86190,
  lng=4.36038,
  radius=50
)

create_inat_leaflet_map <- function(taxon_id, lat, lng, radius) {
  
  # 1. Fetch observations
  observations <- get_inat_obs_radius(
    taxon_id = taxon_id,
    lat      = lat,
    lng      = lng,
    radius   = radius
  )
  
  # 2. Extract lat/lon from the geojson column
  observations_map_data <- 
    observations %>%
    dplyr::mutate(
      lon = purrr::map_dbl(.data$geojson$coordinates, ~ purrr::chuck(.x, 1)),
      lat = purrr::map_dbl(.data$geojson$coordinates, ~ purrr::chuck(.x, 2))
    )

  # 3. Make a map
  leaflet::leaflet(observations_map_data) %>%
    leaflet::addTiles() %>%
    leaflet::addMarkers(
      lng   = ~lon,
      lat   = ~lat
    )
}

create_inat_leaflet_map(
  taxon_id=157924,
  lat=50.86190,
  lng=4.36038,
  radius=50
)