phenoPhase {phenex}

Phenological Phase Extraction

Description

Extracts phenological metrics based on modelled NDVI values.

Usage

phenoPhase(x, phase, method, threshold, n)

Arguments

x

An object of class ‘NDVI’ containing modelled NDVI values.

phase

Determines which phase will be extracted:

“max”: Day of the year with highest NDVI value is returned.

“maxval”: The highest modelled NDVI value and its standard deviation is returned.

“min”: Day of the year (before day with maximum NDVI value) with lowest NDVI value is returned.

“minval”: The lowest modelled NDVI value and its standard deviation is returned.

“greenup”: Day of the year at which greenup takes place and its standard deviation is returned. This day refers to the point where the function of modelled NDVI values exceeds a certain threshold.

“senescence”: Day of the year at which senescence takes place and its standard deviation is returned. This day refers to the point where the function of modelled NDVI values exceeds a certain threshold.

method

Determines whether a global or local threshold is used for greenup and senescence extraction.

“global” threshold: The day of the year is returned, where NDVI values are first equal or higher as the value of ‘threshold’. If the threshold is higher than the values of the timeseries, ‘-1’ will be returned.

“local” threshold: The day of the year is returned, for which NDVI values first reach the value of ‘threshold’ (interpreted as percentage) between lowest and highest NDVI value of timeserie. The lowest NDVI value is chosen depending on phase selected. For “greenup”, the lowest value before day of maximum NDVI value is used. For “senescence”, the lowest value after day of maximum NDVI value is used.

threshold

Threshold for local or global greenup/senescence-extraction method. Should have numerical value between ‘0’ and ‘1’. Global thresholds refer to a fix NDVI value independent of actual NDVI magnitudes. Local thresholds are based on (max - min) * e.g. 0.55. Hence, the respective NDVI value will vary dependent on e.g. land-cover.

n

The number ‘n’ of normal distributed values around the threshold for estimation of ‘sd’. The normal distribution uses the threshold as mean and a standaed deviation consisting of satellite error and standard deviation of fit residuals. The satellite error is considered as 0.02+0.02*value.

Value

A list containing the julian day of the year at which the phenologial phase occurs as list entry ‘mean’ and its standard deviation as list entry ‘sd’. The list contains vectors in ‘mean’ and ‘sd’ if multiple seasons are available in ‘NDVI’ object.

Author(s)

Daniel Doktor, Maximilian Lange

See Also

NDVI, detectSeasons, seasons

Examples

	# load data
	data(avhrr)

	# create NDVI object, correct and model NDVI values
	ndvi <- modelNDVI(ndvi.values=avhrr.ndvi/10000, year.int=1995, 
				correction="bise", method="LinIP", MARGIN=2, 
				doParallel=FALSE, slidingperiod=40)[[1]]

	# extract greenup DOY
	greenup <- phenoPhase(ndvi, phase="greenup", method="local", threshold=0.55, n=1000)