annotate | R Documentation |
Annotate point data with rasters based on matching dates associated with points to dates associated with rasters. Specifically, we're thinking of the points as species observations and the rasters as remotely sensed environmental layers, but they can represent any points and rasters with dates.
annotate(datedOccs, env, envDates, dateScale)
datedOccs |
a 'SpatialPointsDataFrame' of occurrence localities (generally longitude and latitude in decimal degrees) paired with dates. One column must be labeled 'date' and have class 'POSIXct', e.g., as obtained from using 'lubridate::parse_date_time' |
env |
a raster stack |
envDates |
a vector of dates the same length as 'env'. The vector should have class 'POSIXct', e.g., as obtained from using 'lubridate::parse_date_time' |
dateScale |
string: 'year', 'month', or 'day' |
See Examples.
a SpatialPointsDataFrame
Cory Merow <cory.merow@gmail.com>,
r1 <- raster::raster(nrows=50, ncols=50, xmn=-50, xmx=50) raster::values(r1)<- runif(n = (50*50)) r2 <- raster::raster(nrows=50, ncols=50, xmn=-50, xmx=50) raster::values(r2)<- runif(n = (50*50)) env <- raster::stack(r1,r2) names(env) <- c("1995","1996") datedOccs <- data.frame(cbind(c(0,10), c(-10,15))) colnames(datedOccs) <- c("long", "lat") datedOccs$date <- c("1995", "1996") datedOccs$date <- lubridate::parse_date_time(datedOccs$date, orders = c("Y", "Ym")) sp::coordinates(datedOccs) <- c("long", "lat") raster::projection(datedOccs) <- raster::projection(env) dateScale = "year" envDates <- c("1995","1996") annotate(datedOccs = datedOccs, env = env, envDates = envDates, dateScale = dateScale)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.