tMoveRes: tMoveRes

Description

Tool to support the selection of an adequate satellite temporal resolution. It evaluates how the change in temporal resolution changes the amount of samples and sample regions based on a set of coordinate pairs and their observation dates.

Usage

tMoveRes(xy, obs.date, time.res, pixel.res)

Arguments

Object of class SpatialPoints or SpatialPointsDataFrame.

obs.date

Object of class Date with xy observation dates.

time.res

Vector of temporal resolutions (expressed in days).

pixel.res

Spatial resolution (unit depends on spatial projection).

Value

A list object reporting on the amount and distribution of unique pixels and connected pixel regions per temporal resolution.

Details

Given a base spatial resolution (pixel.res and a vector of temporal resolutions (time.res), the function determines the number of unique pixels and unique pixel regions after their temporal aggregation. For each temporal resolution, the function starts by converting xy to unique pixel coordinates and labels them based on their spatial aggregation. Then, the function counts the number of samples and sample regions. The output of the function consists of:

stats - Summary statistics reporting on the number of temporal widows, unique samples and unique sample regions per temporal resolution.
plot - Plot representing the change in number of samples and sample regions per temporal resolution.
indices - Indices for each sample in xy based on their spatial aggregation within each temporal resolution.

Examples

Run this code

# NOT RUN {
{

 require(raster)

 # reference data
 data(longMove)

 # test function for intervals of 1, 8 and 16 days (e.g. of MODIS data)
 obs.date <- as.Date(longMove@data$timestamp)
 a.res <- tMoveRes(longMove, obs.date, c(1,8,16), 0.01)

}
# }