This routine combines a forward and backward run of transition_dates function to calculate the phenophases in both rising and falling parts of a PhenoCam time series.
phenophases(data, mat, internal = TRUE, out_dir = tempdir(), ...)Arguments
- data
a PhenoCam data file (or data frame)
- mat
mean annual temperature
- internal
return PhenoCam data file or data frame
- out_dir
output directory
- ...
pass parameters to the transition_dates() function
Value
Estimates of transition dates for both rising and falling parts of a PhenoCam time series. All time series are evaluated (gcc_90, gcc_75, etc). The function returns a nested list with UNIX time based values including uncertainties on these estimates and their associated thresholds. When written to disk UNIX dates are converted to YYYY-MM-DD. The nested list has named locations rising and falling, or location 1 and 2 in the list respectivelly.
Examples
if (FALSE) {
# downloads a time series
download_phenocam(site = "harvard$",
veg_type = "DB",
roi_id = "1000",
frequency = "3")
# read in data as data frame and calculate phenophases
df <- read_phenocam(file.path(tempdir(),"harvard_DB_1000_3day.csv"))
my_dates <- phenophases(df, internal = TRUE)
# print results
print(my_dates)
}