neon data data download eample

getNEONdata.R

library(tidyverse)
library(neonUtilities)
library(geoNEON)

sites = c('JORN','TALL','BART','BARR', 'HARV','YELL')

for(site_id in sites){
woody_structure=tryCatch(neonUtilities::loadByProduct('DP1.10098.001', site=site_id,
                                                      startdate = '2018-01', enddate = '2018-12',
                                                      check.size = F))

heights = woody_structure$vst_apparentindividual %>%
  group_by(individualID, plotID, eventID) %>%
  summarise(avg_height = mean(height),
            n_stems = n()) %>%
  ungroup()

write_csv(heights,  paste0('./',site_id,'.csv'))
}