Using upper quantile (default = 0.925) of multi-year MODIS data, determine the top NDVI for each id.

filter_top(DT, probs = 0.925, id = "id")

Arguments

DT

data.table of NDVI time series

probs

quantile probability to determine top. default is 0.925.

id

id column. default is 'id'. See details.

Value

filtered data.table with appended 'top' column of each id's top (quantile) NDVI value.

Details

The id argument is used to split between sampling units. This may be a point id, polygon id, pixel id, etc. depending on your analysis.

See also

Examples

# Load data.table library(data.table) # Read example data ndvi <- fread(system.file("extdata", "ndvi.csv", package = "irg")) filter_qa(ndvi, qa = 'SummaryQA', good = c(0, 1)) filter_winter(ndvi, probs = 0.025, limits = c(60L, 300L), doy = 'DayOfYear', id = 'id')
#> id yr DayOfYear NDVI SummaryQA filtered winter #> 1: 7652 2002 3 -1367 3 4099 4099 #> 2: 7035 2002 14 -304 3 5382 5382 #> 3: 4368 2002 1 374 2 3702 3702 #> 4: 5196 2002 15 635 3 5180 5180 #> 5: 12377 2002 9 685 2 4621 4621 #> --- #> 2526: 13322 2012 356 1257 2 4708 4708 #> 2527: 9537 2012 356 144 2 4536 4536 #> 2528: 18592 2012 362 1031 2 4011 4011 #> 2529: 2163 2012 356 2172 2 5022 5022 #> 2530: 2941 2012 356 -197 2 2839 2839
filter_roll(ndvi, window = 3L, id = 'id')
#> id yr DayOfYear NDVI SummaryQA filtered winter rolled #> 1: 7652 2002 3 -1367 3 4099 4099 4099 #> 2: 7035 2002 14 -304 3 5382 5382 5382 #> 3: 4368 2002 1 374 2 3702 3702 3702 #> 4: 5196 2002 15 635 3 5180 5180 5180 #> 5: 12377 2002 9 685 2 4621 4621 4621 #> --- #> 2526: 13322 2012 356 1257 2 4708 4708 4708 #> 2527: 9537 2012 356 144 2 4536 4536 4536 #> 2528: 18592 2012 362 1031 2 4011 4011 4011 #> 2529: 2163 2012 356 2172 2 5022 5022 5022 #> 2530: 2941 2012 356 -197 2 2839 2839 2839
filter_top(ndvi, probs = 0.925, id = 'id')
#> id yr DayOfYear NDVI SummaryQA filtered winter rolled top #> 1: 7652 2002 3 -1367 3 4099 4099 4099 7443.400 #> 2: 7035 2002 14 -304 3 5382 5382 5382 7462.900 #> 3: 4368 2002 1 374 2 3702 3702 3702 6709.800 #> 4: 5196 2002 15 635 3 5180 5180 5180 7288.100 #> 5: 12377 2002 9 685 2 4621 4621 4621 7645.100 #> --- #> 2526: 13322 2012 356 1257 2 4708 4708 4708 7776.250 #> 2527: 9537 2012 356 144 2 4536 4536 4536 7155.475 #> 2528: 18592 2012 362 1031 2 4011 4011 4011 7184.300 #> 2529: 2163 2012 356 2172 2 5022 5022 5022 7566.600 #> 2530: 2941 2012 356 -197 2 2839 2839 2839 5491.750