diff --git a/.Rbuildignore b/.Rbuildignore index a7d7e9a..6cf7893 100644 --- a/.Rbuildignore +++ b/.Rbuildignore @@ -24,3 +24,8 @@ examples_x64 ^LICENSE\.md$ ^dtwSat\.Rproj$ ^\.github$ +^doc$ +^Meta$ + +vignettes.awk +_pkgdown.yml \ No newline at end of file diff --git a/.github/workflows/pkgdown.yaml b/.github/workflows/pkgdown.yaml new file mode 100644 index 0000000..ed7650c --- /dev/null +++ b/.github/workflows/pkgdown.yaml @@ -0,0 +1,48 @@ +# Workflow derived from https://github.com/r-lib/actions/tree/v2/examples +# Need help debugging build failures? Start at https://github.com/r-lib/actions#where-to-find-help +on: + push: + branches: [main, master] + pull_request: + branches: [main, master] + release: + types: [published] + workflow_dispatch: + +name: pkgdown + +jobs: + pkgdown: + runs-on: ubuntu-latest + # Only restrict concurrency for non-PR jobs + concurrency: + group: pkgdown-${{ github.event_name != 'pull_request' || github.run_id }} + env: + GITHUB_PAT: ${{ secrets.GITHUB_TOKEN }} + permissions: + contents: write + steps: + - uses: actions/checkout@v3 + + - uses: r-lib/actions/setup-pandoc@v2 + + - uses: r-lib/actions/setup-r@v2 + with: + use-public-rspm: true + + - uses: r-lib/actions/setup-r-dependencies@v2 + with: + extra-packages: any::pkgdown, local::. + needs: website + + - name: Build site + run: pkgdown::build_site_github_pages(new_process = FALSE, install = FALSE) + shell: Rscript {0} + + - name: Deploy to GitHub pages 🚀 + if: github.event_name != 'pull_request' + uses: JamesIves/github-pages-deploy-action@v4.4.1 + with: + clean: false + branch: gh-pages + folder: docs diff --git a/.gitignore b/.gitignore index 527b631..29aaef1 100644 --- a/.gitignore +++ b/.gitignore @@ -37,3 +37,5 @@ CRAN-SUBMISSION src/symbols.rds *.o *.so +/doc/ +/Meta/ diff --git a/DESCRIPTION b/DESCRIPTION index 4e08634..f4ba85e 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -2,7 +2,7 @@ Package: dtwSat Type: Package Title: Time-Weighted Dynamic Time Warping for Satellite Image Time Series Analysis Version: 1.0.0 -Date: 2023-09-03 +Date: 2023-09-21 Authors@R: c(person(given = "Victor", family = "Maus", @@ -24,18 +24,19 @@ Authors@R: ) Description: Provides a robust approach to land use mapping using multi-dimensional (multi-band) satellite image time series. By leveraging the Time-Weighted Dynamic - Time Warping (TWDTW) distance metric in tandem with a 1-NN classifier, the package - provides functions to produce land use maps based on distinct seasonality patterns, - typically observed in the phenological cycles of vegetation. The TWDTW algorithm is - described in Maus et al. (2016) and - Maus et al. (2019) . A key strength of TWDTW is its ability - to recognize patterns with only a minimal training set, achieving notable accuracy. - The package features tools for generating temporal patterns for various land cover types, - conducting land use mapping, and visualizing the outcomes. + Time Warping (TWDTW) distance metric in tandem with a 1 Nearest-Neighbor (1-NN) Classifier, + this package offers functions to produce land use maps based on distinct seasonality patterns, + commonly observed in the phenological cycles of vegetation. The approach is described in + Maus et al. (2016) and Maus et al. (2019) . + A primary advantage of TWDTW is its capability to handle irregularly sampled and noisy time series, + while also requiring minimal training sets. The package includes tools for training the 1-NN-TWDTW model, + visualizing temporal patterns, producing land use maps, and visualizing the results. License: GPL (>= 3) URL: https://github.com/vwmaus/dtwSat/ BugReports: https://github.com/vwmaus/dtwSat/issues/ Maintainer: Victor Maus +VignetteBuilder: + knitr Encoding: UTF-8 Roxygen: list(markdown = TRUE) RoxygenNote: 7.2.3 @@ -47,11 +48,17 @@ Depends: Imports: mgcv, stats, - scales, - reshape2, - rlang + tidyr, + proxy Suggests: - rbenchmark, - stringr, + knitr, + rmarkdown, testthat (>= 3.0.0) Config/testthat/edition: 3 +Collate: + 'plot.R' + 'predict.R' + 'prepare_time_series.R' + 'shift_dates.R' + 'train.R' + 'zzz.R' diff --git a/NAMESPACE b/NAMESPACE index fb17b02..4c99aae 100644 --- a/NAMESPACE +++ b/NAMESPACE @@ -1,8 +1,10 @@ # Generated by roxygen2: do not edit by hand -export(create_patterns) -export(plot_patterns) +S3method(plot,twdtw_knn1) +S3method(predict,twdtw_knn1) +S3method(print,twdtw_knn1) export(shift_dates) +export(twdtw_knn1) import(ggplot2) import(sf) import(stars) @@ -10,11 +12,11 @@ import(twdtw) importFrom(mgcv,gam) importFrom(mgcv,predict.gam) importFrom(mgcv,s) -importFrom(reshape2,melt) -importFrom(rlang,.data) -importFrom(scales,date_format) -importFrom(scales,percent) -importFrom(scales,pretty_breaks) +importFrom(proxy,dist) importFrom(stats,as.formula) importFrom(stats,predict) importFrom(stats,setNames) +importFrom(tidyr,nest) +importFrom(tidyr,pivot_longer) +importFrom(tidyr,pivot_wider) +importFrom(tidyr,unnest) diff --git a/R/create_patterns.R b/R/create_patterns.R deleted file mode 100644 index 856d1dd..0000000 --- a/R/create_patterns.R +++ /dev/null @@ -1,146 +0,0 @@ -#' Create a Pattern Using GAM -#' -#' This function creates a pattern based on Generalized Additive Models (GAM). -#' It uses the specified formula to fit the model and predict values. -#' -#' @param x A three dimensions stars object (x, y, time) with the satellite image time series. -#' @param y An sf object with the coordinates of the training points. -#' @param formula A formula for the GAM. Default is \code{band ~ \link[mgcv]{s}(time)}. -#' @param start_column Name of the column in y that indicates the start date. Default is 'start_date'. -#' @param end_column Name of the column in y that indicates the end date. Default is 'end_date'. -#' @param label_colum Name of the column in y that contains land use labels. Default is 'label'. -#' @param sampling_freq The time frequency for sampling including unit, e.g '16 day'. If NULL, the function will infer it. -#' @param ... Additional arguments passed to the GAM function. -#' -#' @return A list containing the predicted values for each label. -#' -#' -#' -#' @export -create_patterns = function(x, y, formula = band ~ s(time), start_column = 'start_date', - end_column = 'end_date', label_colum = 'label', - sampling_freq = NULL, ...){ - - # Check if x is a stars object with a time dimension - if (!inherits(x, "stars") || dim(x)['time'] < 1) { - stop("x must be a stars object with a 'time' dimension") - } - - # Check if y is an sf object with point geometry - if (!inherits(y, "sf") || !all(st_is(y, "POINT"))) { - stop("y must be an sf object with point geometry") - } - - # Check for required columns in y - required_columns <- c(start_column, end_column, label_colum) - missing_columns <- setdiff(required_columns, names(y)) - if (length(missing_columns) > 0) { - stop(paste("Missing required columns in y:", paste(missing_columns, collapse = ", "))) - } - - # Check if formula has two - if(length(all.vars(formula)) != 2) { - stop("The formula should have only one predictor!") - } - - # Convert columns to date-time - y[ , start_column] <- to_date_time(y[[start_column]]) - y[ , end_column] <- to_date_time(y[[end_column]]) - - # Extract time series from stars - y_ts <- extract_time_series(x, y) - y_ts$geom <- NULL - - # Shift dates - unique_ids <- unique(y_ts$id) - for (id in seq_along(unique_ids)) { - idx <- y_ts$id == unique_ids[id] - y_ts[idx, 'time'] <- shift_dates(y_ts[idx, 'time']) - } - - # Split data frame by label and remove label column - y_ts <- lapply(split(y_ts, y_ts$label), function(df) { - df$label <- NULL - return(df) - }) - - # Determine sampling frequency - if (is.null(sampling_freq)) { - sampling_freq <- get_stars_time_freq(x) - cat("Sampling frequency inferred from the stars object:", as.numeric(sampling_freq), attr(sampling_freq, "units"), "\n") - } - - # Define GAM function - gam_fun <- function(y, t, formula, ...){ - df <- data.frame(y, t = as.numeric(t)) - df <- setNames(list(y, as.numeric(t)), all.vars(formula)) - fit <- mgcv::gam(data = df, formula = formula, ...) - pred_t <- setNames(list(as.numeric(seq(min(t), max(t), by = sampling_freq))), all.vars(formula)[2]) - predict(fit, newdata = pred_t) - } - - # Apply GAM function - patterns <- lapply(y_ts, function(ts){ - y_time <- ts$time - ts$time <- NULL - ts$id <- NULL - sapply(as.list(ts), function(y) { - gam_fun(y, y_time, formula, ...) - }) - }) - - return(patterns) - -} - - - -#' Extract Time Series from a Stars Object for Specified Points -#' -#' This function extracts a time series from a stars object for each specified point in the sf object. -#' Each extracted sample is then labeled with an ID and the label from the sf object. -#' -#' @param x A stars object containing the raster time series data. -#' @param y An sf object containing the point geometries and their associated labels. -#' -#' @return A data.frame with the extracted time series for each point in the sf object, -#' with additional columns for the ID and label of each sample. -#' -#' -extract_time_series <- function(x, y) { - ts_samples <- st_extract(x, y) - ts_samples$id <- 1:dim(ts_samples)["geom"] - ts_samples$label <- y$label - as.data.frame(ts_samples) -} - - - -#' Compute the Most Common Sampling Frequency in a Stars Object -#' -#' This function calculates the most common difference between consecutive time points in a stars object. -#' This can be useful for determining the sampling frequency of the time series data. -#' -#' @param x A stars object containing time series data. -#' -#' @return A difftime object representing the most common time difference between consecutive samples. -#' -#' -get_stars_time_freq <- function(x) { - - # Extract the time dimension - time_values <- st_get_dimension_values(x, "time") - - # Compute the differences between consecutive time points - time_diffs <- diff(time_values) - - # Convert differences to days (while retaining the difftime class) - time_diffs <- as.difftime(time_diffs, units = "days") - - # Identify the mode - mode_val_index <- which.max(tabulate(match(time_diffs, unique(time_diffs)))) - freq <- diff(time_values[mode_val_index:(mode_val_index+1)]) - - return(freq) -} - diff --git a/R/plot.R b/R/plot.R new file mode 100644 index 0000000..70a2fc9 --- /dev/null +++ b/R/plot.R @@ -0,0 +1,48 @@ +#' Plot Patterns from twdtw-knn1 model +#' +#' This function visualizes time series patterns from the \code{"twdtw_knn1"} model. +#' It produces a multi-faceted plot, where each facet represents a different time series +#' label from the model's data. Within each facet, different bands or indices (attributes) +#' are plotted as distinct lines, differentiated by color. +#' +#' @param x A model of class \code{"twdtw_knn1"}. +#' +#' @param bands A character vector specifying the bands or indices to plot. +#' If NULL (default), all available bands or indices in the data will be plotted. +#' +#' @param ... Additional arguments passed to \code{\link[ggplot2]{ggplot}}. Currently not used. +#' +#' @return A \code{\link[ggplot2]{ggplot}} object displaying the time series patterns. +#' +#' @seealso twdtw_knn1 +#' +#' @inherit twdtw_knn1 examples +#' +#' @export +plot.twdtw_knn1 <- function(x, bands = NULL, ...) { + + # Convert the list of time series data into a long-format data.frame + df <- x$data + df$id <- 1:nrow(df) + df <- unnest(df, cols = 'observations') + + # Select bands + if(!is.null(bands)){ + df <- df[c('id', 'time', 'label', bands)] + } + + # Pivote data into long format for ggplot2 + df <- pivot_longer(df, !c('id', 'label', 'time'), names_to = "band", values_to = "value") + + # Construct the ggplot + gp <- ggplot(df, aes(x = .data$time, y = .data$value, colour = .data$band, group = interaction(.data$id, .data$band))) + + geom_line() + + facet_wrap(~label) + + theme(legend.position = "bottom") + + guides(colour = guide_legend(title = "Bands")) + + ylab("Value") + + xlab("Time") + + return(gp) + +} diff --git a/R/plot_patterns.R b/R/plot_patterns.R deleted file mode 100644 index dd983b8..0000000 --- a/R/plot_patterns.R +++ /dev/null @@ -1,41 +0,0 @@ -#' Plot Patterns from Time Series Data -#' -#' This function takes a list of time series data and creates a multi-faceted plot -#' where each facet corresponds to a different time series from the list. -#' Within each facet, different attributes (columns of the time series) are -#' plotted as lines with different colors. -#' -#' @param x A list where each element is a data.frame representing a time series. -#' Each data.frame should have the same number of rows and columns, -#' with columns representing different attributes (e.g., bands or indices) -#' and rows representing time points. -#' The name of each element in the list will be used as the facet title. -#' -#' @param ... Not used. -#' -#' @return A ggplot object displaying the time series patterns. -#' -#' @export -plot_patterns = function(x, ...) { - - # Convert the list of time series data into a long-format data.frame - df.p = do.call("rbind", lapply(names(x), function(p) { - ts = x[[p]] - # Create a new data.frame with a 'Time' column and a 'Pattern' column - # representing the name of the current time series (facet name). - data.frame(Time = 1:nrow(ts), ts, Pattern = p) # Assuming the time series are evenly spaced - })) - - # Melt the data into long format suitable for ggplot2 - df.p = melt(df.p, id.vars = c("Time", "Pattern")) - - # Construct the ggplot - gp = ggplot(df.p, aes(x = .data$Time, y = .data$value, colour = .data$variable)) + - geom_line() + - facet_wrap(~Pattern) + - theme(legend.position = "bottom") + - guides(colour = guide_legend(title = "Bands")) + - ylab("Value") - - return(gp) -} \ No newline at end of file diff --git a/R/predict.R b/R/predict.R new file mode 100644 index 0000000..3e25220 --- /dev/null +++ b/R/predict.R @@ -0,0 +1,42 @@ +#' Predict using the twdtw_knn1 model +#' +#' This function predicts the classes of new data using the Time Warped Dynamic Time Warping (TWDTW) +#' method with a 1-nearest neighbor approach. The prediction is based on the minimum TWDTW distance +#' to the known patterns stored in the `twdtw_knn1` model. +#' +#' @param object A `twdtw_knn1` model object generated by the `twdtw_knn1` function. +#' @param newdata A data frame or similar object containing the new observations +#' (time series data) to be predicted. +#' @param ... Additional arguments passed to the \link[twdtw]{twdtw} function. +#' If provided, they will overwrite twdtw arguments previously passed to \link[dtwSat]{twdtw_knn1}. +#' +#' @return A vector of predicted classes for the `newdata`. +#' +#' @seealso twdtw_knn1 +#' +#' @inherit twdtw_knn1 examples +#' +#' @export +predict.twdtw_knn1 <- function(object, newdata, ...){ + + # Update twdtw_args with new arguments passed via ... + new_twdtw_args <- list(...) + matching_twdtw_args <- intersect(names(new_twdtw_args), names(object$twdtw_args)) + object$twdtw_args[matching_twdtw_args] <- new_twdtw_args[matching_twdtw_args] + + # Convert newdata to time series + newdata_ts <- prepare_time_series(newdata) + + # Compute TWDTW distances + distances <- sapply(object$data$observations, function(pattern){ + sapply(newdata_ts$observations, function(ts) { + do.call(proxy::dist, c(list(x = as.data.frame(ts), y = as.data.frame(pattern), method = 'twdtw'), object$twdtw_args)) + }) + }) + + # Find the nearest neighbor for each observation in newdata + nearest_neighbor <- apply(distances, 1, which.min) + + # Return the predicted label for each observation based on the nearest neighbor + return(factor(object$data$label[nearest_neighbor])) +} diff --git a/R/prepare_time_series.R b/R/prepare_time_series.R new file mode 100644 index 0000000..ed76a56 --- /dev/null +++ b/R/prepare_time_series.R @@ -0,0 +1,83 @@ +#' Prepare a Time Series Tibble from a 2D stars Object with Bands and Time Attributes +#' +#' This function reshapes a data frame, which has been converted from a stars object, into a nested wide tibble format. +#' The stars object conversion often results in columns named in formats like "band.YYYY.MM.DD", "XYYYY.MM.DD.band", or "YYYY.MM.DD.band". +#' +#' @param x A data frame derived from a stars object containing time series data in wide format. +#' The column names should adhere to one of the following formats: "band.YYYY.MM.DD", "XYYYY.MM.DD.band", or "YYYY.MM.DD.band". +#' +#' @return A nested tibble in wide format. Each row of the tibble corresponds to a unique 'ts_id' that maintains the order from the original stars object. +#' The nested structure contains observations (time series) for each 'ts_id', including the 'time' of each observation, and individual bands are presented as separate columns. +#' +#' +prepare_time_series <- function(x) { + + # Remove the 'geom' column if it exists + x$geom <- NULL + x$x <- NULL + x$y <- NULL + var_names <- names(x) + var_names <- var_names[!var_names %in% 'label'] + + # Extract date and band information from the column names + date_band <- do.call(rbind, lapply(var_names, function(name) { + + # Replace any hyphen with periods for consistent processing + name <- gsub("-", "\\.", name) + + # Extract date and band info based on different name patterns + if (grepl("^.+\\.\\d{4}\\.\\d{2}\\.\\d{2}$", name)) { + date_str <- gsub("^.*?(\\d{4}\\.\\d{2}\\.\\d{2})$", "\\1", name) + band_str <- gsub("^(.+)\\.\\d{4}\\.\\d{2}\\.\\d{2}$", "\\1", name) + } + else if (grepl("^X\\d{4}\\.\\d{2}\\.\\d{2}\\..+$", name)) { + date_str <- gsub("^X(\\d{4}\\.\\d{2}\\.\\d{2})\\..+$", "\\1", name) + band_str <- gsub("^X\\d{4}\\.\\d{2}\\.\\d{2}\\.(.+)$", "\\1", name) + } + else if (grepl("^\\d{4}\\.\\d{2}\\.\\d{2}\\..+$", name)) { + date_str <- gsub("^(\\d{4}\\.\\d{2}\\.\\d{2})\\..+$", "\\1", name) + band_str <- gsub("^\\d{4}\\.\\d{2}\\.\\d{2}\\.(.+)$", "\\1", name) + } + else { + stop(paste("Unrecognized format in:", name)) + } + + # Convert the date string to Date format + date <- to_date_time(gsub("\\.", "-", date_str)) + + return(data.frame(time = date, band = band_str)) + })) + + # Construct tiem sereis + ns <- nrow(x) + x$ts_id <- 1:ns + if (!'label' %in% names(x)) { + x$label <- NA + } + x <- pivot_longer(x, !c('ts_id', 'label'), names_to = 'band_date', values_to = 'value') + x$band <- rep(date_band$band, ns) + x$time <- rep(date_band$time, ns) + x$band_date <- NULL + result_df <- pivot_wider(x, id_cols = c('ts_id', 'label', 'time'), names_from = 'band', values_from = 'value') + result_df <- nest(result_df, .by = c('ts_id', 'label'), .key = 'observations') + + return(result_df) + +} + + +#### TO BE REMOVED: twdtw package will export this fucntion +to_date_time <- function(x){ + if (!inherits(x, c("Date", "POSIXt"))) { + # check if all strings in the vector include hours, minutes, and seconds + if (all(grepl("\\d{4}-\\d{2}-\\d{2} \\d{2}:\\d{2}:\\d{2}", x))) { + x <- try(as.POSIXct(x), silent = TRUE) + } else { + x <- try(as.Date(x), silent = TRUE) + } + if (inherits(x, "try-error")) { + stop("Some elements of x could not be converted to a date or datetime format") + } + } + return(x) +} diff --git a/R/utils.R b/R/shift_dates.R similarity index 76% rename from R/utils.R rename to R/shift_dates.R index b9bb556..a37e060 100644 --- a/R/utils.R +++ b/R/shift_dates.R @@ -43,22 +43,8 @@ shift_dates <- function(x, origin = "1970-01-01") { return(shifted_dates) } -#### TO BE REMOVED: twdtw package will export this fucntion -to_date_time <- function(x){ - if (!inherits(x, c("Date", "POSIXt"))) { - # check if all strings in the vector include hours, minutes, and seconds - if (all(grepl("\\d{4}-\\d{2}-\\d{2} \\d{2}:\\d{2}:\\d{2}", x))) { - x <- try(as.POSIXct(x), silent = TRUE) - } else { - x <- try(as.Date(x), silent = TRUE) - } - if (inherits(x, "try-error")) { - stop("Some elements of x could not be converted to a date or datetime format") - } - } + +shift_ts_dates <- function(x) { + x$time <- shift_dates(x$time) return(x) } - - - - diff --git a/R/train.R b/R/train.R new file mode 100644 index 0000000..e6a8454 --- /dev/null +++ b/R/train.R @@ -0,0 +1,274 @@ +#' +#' Train a KNN-1 TWDTW model with optional GAM resampling +#' +#' This function prepares a KNN-1 model with the Time Warp Dynamic Time Warping (TWDTW) algorithm. +#' If a formula is provided, the training samples are resampled using Generalized Additive Models (GAM). +#' +#' @param x A three-dimensional stars object (x, y, time) with bands as attributes. +#' @param y An sf object with the coordinates of the training points. +#' @param time_weight A numeric vector with length two (steepness and midpoint of logistic weight) or a function. +#' See details in \link[twdtw]{twdtw}. +#' @param cycle_length The length of the cycle, e.g. phenological cycles. Details in \link[twdtw]{twdtw}. +#' @param time_scale Specifies the time scale for the observations. Details in \link[twdtw]{twdtw}. +#' @param formula Either NULL or a formula to reduce samples of the same label using Generalized Additive Models (GAM). +#' Default is \code{band ~ s(time)}. See details. +#' @param start_column Name of the column in y that indicates the start date. Default is 'start_date'. +#' @param end_column Name of the column in y that indicates the end date. Default is 'end_date'. +#' @param label_colum Name of the column in y containing land use labels. Default is 'label'. +#' @param sampling_freq The time frequency for sampling, including the unit (e.g., '16 day'). +#' If NULL, the function will infer the frequency. This parameter is only used if a formula is provided. +#' @param ... Additional arguments passed to the \link[mgcv]{gam} function and to \link[twdtw]{twdtw} function. +#' +#' @details If \code{formula} is NULL, the KNN-1 model will retain all training samples. If a formula is passed (e.g., \code{band ~ \link[mgcv]{s}(time)}), +#' then samples of the same label (land cover class) will be resampled using GAM. +#' Resampling can significantly reduce prediction processing time. +#' +#' @return A 'twdtw_knn1' model containing the trained model information and the data used. +#' +#' @examples +#' \dontrun{ +#' +#' # Read training samples +#' samples <- +#' system.file("mato_grosso_brazil/samples.gpkg", package = "dtwSat") |> +#' st_read(quiet = TRUE) +#' +#' # Get satellite image time sereis files +#' tif_files <- system.file("mato_grosso_brazil", package = "dtwSat") |> +#' dir(pattern = "\\.tif$", full.names = TRUE) +#' +#' # Get acquisition dates +#' acquisition_date <- regmatches(tif_files, regexpr("[0-9]{8}", tif_files)) |> +#' as.Date(format = "%Y%m%d") +#' +#' # Create a 3D datacube +#' dc <- read_stars(tif_files, +#' proxy = FALSE, +#' along = list(time = acquisition_date), +#' RasterIO = list(bands = 1:6)) |> +#' st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")) |> +#' split(c("band")) +#' +#' # Create a knn1-twdtw model +#' m <- twdtw_knn1(x = dc, +#' y = samples, +#' cycle_length = 'year', +#' time_scale = 'day', +#' time_weight = c(steepness = 0.1, midpoint = 50), +#' formula = band ~ s(time)) +#' +#' print(m) +#' +#' # Visualize model patterns +#' plot(m) +#' +#' # Classify satellite images +#' system.time(lu <- predict(dc, model = m)) +#' +#' # Visualise land use classification +#' ggplot() + +#' geom_stars(data = lu) + +#' theme_minimal() +#' +#' } +#' @export +twdtw_knn1 <- function(x, y, time_weight, cycle_length, time_scale, + formula = NULL, start_column = 'start_date', + end_column = 'end_date', label_colum = 'label', + sampling_freq = NULL, ...){ + + # Check if x is a stars object with a time dimension + if (!inherits(x, "stars") || dim(x)['time'] < 1 || length(dim(x)) != 3) { + stop("x must be a three-dimensional stars object with a 'time' dimension") + } + + x <- split(x, c("time")) + + # Check if y is an sf object with point geometry + if (!inherits(y, "sf") || !all(st_is(y, "POINT"))) { + stop("y must be an sf object with point geometry") + } + + # check for minimum set of twdtw arguments + if (!(is.function(time_weight) || (is.numeric(time_weight) && length(time_weight) == 2))) stop("'time_weight' should be either a function or a numeric vector with length two") + if (is.null(cycle_length)) stop("The 'cycle_length' argument is missing.") + if (is.null(time_scale)) stop("The 'time_scale' argument is missing for 'cycle_length' type character.") + + # Check for required columns in y + required_columns <- c(start_column, end_column, label_colum) + missing_columns <- setdiff(required_columns, names(y)) + if (length(missing_columns) > 0) { + stop(paste("Missing required columns in y:", paste(missing_columns, collapse = ", "))) + } + + # Convert columns to date-time + y[, start_column] <- to_date_time(y[[start_column]]) + y[, end_column] <- to_date_time(y[[end_column]]) + + # Prepare time series samples from stars object + ts_data <- st_extract(x, y) + ts_data$label <- y[[label_colum]] + ts_data <- prepare_time_series(as.data.frame(ts_data)) + ts_data$ts_id <- NULL + + if(!is.null(formula)) { + + # Check if formula has two + if(length(all.vars(formula)) != 2) { + stop("The formula should have only one predictor!") + } + + # Determine sampling frequency + if (is.null(sampling_freq)) { + sampling_freq <- get_time_series_freq(ts_data) + } + + # Shift dates + ts_data$observations <- lapply(ts_data$observations, shift_ts_dates) + + # Split data frame by label + ts_data <- unnest(ts_data, cols = 'observations') + ts_data <- nest(ts_data, .by = 'label', .key = "observations") + + # Define GAM function + gam_fun <- function(band, t, pred_t, formula, ...){ + df <- setNames(list(band, as.numeric(t)), all.vars(formula)) + pred_t[[all.vars(formula)[2]]] <- as.numeric(pred_t[[all.vars(formula)[2]]]) + fit <- mgcv::gam(data = df, formula = formula, ...) + predict(fit, newdata = pred_t) + } + + # Apply GAM function + ts_data$observations <- lapply(ts_data$observations, function(ts){ + y_time <- ts$time + ts$time <- NULL + pred_time <- setNames(list(seq(min(y_time), max(y_time), by = sampling_freq)), all.vars(formula)[2]) + cbind(pred_time, as.data.frame(sapply(as.list(ts), function(band) { + gam_fun(band, y_time, pred_time, formula, ...) + }))) + }) + + } + + model <- list() + model$call <- match.call() + model$formula <- formula + model$data <- ts_data + # add twdtw arguments to model + model$twdtw_args <- list(time_weight = time_weight, + cycle_length = cycle_length, + time_scale = time_scale, + origin = NULL, + max_elapsed = Inf, + version = "f90") + new_twdtw_args <- list(...) + matching_twdtw_args <- intersect(names(model$twdtw_args), names(new_twdtw_args)) + model$twdtw_args[matching_twdtw_args] <- new_twdtw_args[matching_twdtw_args] + + class(model) <- "twdtw_knn1" + + return(model) + +} + +#' Print method for objects of class twdtw_knn1 +#' +#' This method provides a structured printout of the important components +#' of a `twdtw_knn1` object. +#' +#' @param x An object of class `twdtw_knn1`. +#' @param ... ignored +#' +#' @return Invisible `twdtw_knn1` object. +#' +#' @export +print.twdtw_knn1 <- function(x, ...) { + cat("\nModel of class 'twdtw_knn1'\n") + cat("-----------------------------\n") + + # Printing the call + cat("Call:\n") + print(x$call) + + # Printing the formula, if available + cat("\nFormula:\n") + print(x$formula) + + # Printing the data summary + cat("\nData:\n") + print(x$data) + + # Printing twdtw arguments + cat("\nTWDTW Arguments:\n") + pretty_arguments(x$twdtw_args) + + invisible(x) # Returns the object invisibly, so it doesn't print twice +} + +#' Print Pretty Arguments +#' +#' Display a list of arguments of a given function in a human-readable format. +#' +#' @param args A list of named arguments to display. +#' +#' @return Invisible NULL. The function is mainly used for its side effect of printing. +#' +#' @examples +#' \dontrun{ +#' pretty_arguments(formals(twdtw_knn1)) +#' } +#' +pretty_arguments <- function(args) { + + if (is.null(args)) { + cat("Arguments are missing.\n") + return(invisible(NULL)) + } + + for (name in names(args)) { + default_value <- args[[name]] + + if (is.symbol(default_value)) { + default_value <- as.character(default_value) + + } else if (is.null(default_value)) { + default_value <- "NULL" + + } else if (is.vector(default_value) && !is.null(names(default_value))) { + # Handle named vectors + values <- paste(names(default_value), default_value, sep = "=", collapse = ", ") + default_value <- paste0("c(", values, ")") + } + cat(paste0(" - ", name, ": ", default_value, "\n")) + } +} + + + +#' Compute the Most Common Sampling Frequency across all observations +#' +#' This function calculates the most common difference between consecutive time points. +#' This can be useful for determining the aproximate sampling frequency of the time series data. +#' +#' @param x A data frame including a column called `observations`` with the time series +#' +#' @return A difftime object representing the most common time difference between consecutive samples. +#' +get_time_series_freq <- function(x) { + + # Extract the time dimension + time_values <- unlist(lapply(x$observations, function(ts) ts$time)) + + # Compute the differences between consecutive time points + time_diffs <- unlist(lapply(x$observations, function(ts) diff(ts$time))) + + # Convert differences to days (while retaining the difftime class) + time_diffs <- as.difftime(time_diffs, units = "days") + + # Identify the mode + mode_val_index <- which.max(tabulate(match(time_diffs, unique(time_diffs)))) + freq <- diff(time_values[mode_val_index:(mode_val_index+1)]) + + return(freq) + +} diff --git a/R/zzz.R b/R/zzz.R index 499fd6c..90fa6ae 100644 --- a/R/zzz.R +++ b/R/zzz.R @@ -3,9 +3,8 @@ #' @import stars #' @import ggplot2 #' @importFrom stats as.formula predict setNames -#' @importFrom mgcv gam s predict.gam -#' @importFrom scales pretty_breaks date_format percent -#' @importFrom reshape2 melt -#' @importFrom rlang .data +#' @importFrom mgcv gam s predict.gam +#' @importFrom tidyr pivot_longer pivot_wider nest unnest +#' @importFrom proxy dist #' NULL diff --git a/README.md b/README.md index 4db524f..031f02c 100644 --- a/README.md +++ b/README.md @@ -7,33 +7,42 @@ [![CRAN](https://www.r-pkg.org/badges/version/dtwSat)](https://cran.r-project.org/package=dtwSat) [![Downloads](https://cranlogs.r-pkg.org/badges/dtwSat?color=brightgreen)](https://www.r-pkg.org/pkg/dtwSat) [![total](http://cranlogs.r-pkg.org/badges/grand-total/dtwSat)](http://www.r-pkg.org/pkg/dtwSat) +[![cran checks](https://badges.cranchecks.info/worst/dtwSat.svg)](https://cran.r-project.org/web/checks/check_results_dtwSat.html) +[![status](https://tinyverse.netlify.com/badge/dtwSat)](https://CRAN.R-project.org/package=dtwSat) + +Provides a robust approach to land use mapping using multi-dimensional +(multi-band) satellite image time series. By leveraging the Time-Weighted Dynamic +Time Warping (TWDTW) distance metric in tandem with a 1 Nearest-Neighbor (1-NN) Classifier, +this package offers functions to produce land use maps based on distinct seasonality patterns, +commonly observed in the phenological cycles of vegetation. The approach is described in +Maus et al. (2016) and Maus et al. (2019). +A primary advantage of TWDTW is its capability to handle irregularly sampled and noisy time series, +while also requiring minimal training sets. The package includes tools for training the 1-NN-TWDTW model, +visualizing temporal patterns, producing land use maps, and visualizing the results. -### Time-Weighted Dynamic Time Warping for satellite image time series analysis - -The package *dtwSat* provides an implementation of the Time-Weighted -Dynamic Time Warping (TWDTW) method for land cover mapping using -multi-band satellite image time series (Maus et al. 2016, 2019). -*dtwSat* provides full cycle of land cover classification using image -time series, ranging from selecting temporal patterns to visualizing, -and assessing the results. +## Getting Started -## Installing - -Install either from CRAN +You can install dtwSat from CRAN using the following command: ``` r install.packages("dtwSat") ``` -or install the development versions from GitHub +Alternatively, you can install the development version from GitHub: ``` r -library(devtools) devtools::install_github("vwmaus/dtwSat") ``` +After installation, you can read the vignette for a quick start guide: + +``` r +vignette("landuse-mapping", "dtwSat") +``` + + ## References
diff --git a/inst/REFERENCES.bib b/inst/REFERENCES.bib index c2dacc5..d1931b3 100644 --- a/inst/REFERENCES.bib +++ b/inst/REFERENCES.bib @@ -23,47 +23,6 @@ @article{Maus:2019 month={Jan} } -@article{Olofsson:2014, - title = {Good Practices for Estimating Area and Assessing Accuracy of Land Change}, - journal = {Remote Sensing of Environment}, - volume = {148}, - pages = {42--57}, - year = {2014}, - issn = {0034--4257}, - doi = {10.1016/j.rse.2014.02.015}, - author = {Pontus Olofsson and Giles M. Foody and Martin Herold and Stephen V. Stehman and Curtis E. Woodcock and Michael A. Wulder} -} - -@article{Olofsson:2013, - author = {Pontus Olofsson and Giles M. Foody and Stephen V. Stehman and Curtis E. Woodcock}, - title = {Making Better Use of Accuracy Aata in Land Change Studies: Estimating Accuracy and Area and Quantifying Uncertainty Using Stratified Estimation}, - journal = {Remote Sensing of Environment}, - volume = {129}, - pages = {122--131}, - year = {2013}, - doi = {10.1016/j.rse.2012.10.031} -} - -@article{Giorgino:2009, - author = {Toni Giorgino}, - title = {Computing and Visualizing Dynamic Time Warping Alignments in {R}: The dtw Package}, - journal = {Journal of Statistical Software}, - year = {2009}, - volume = {31}, - pages = {1--24}, - number = {7}, - doi = {10.18637/jss.v031.i07} -} - -@book{Muller:2007, - title = {Information Retrieval for Music and Motion}, - publisher = {Springer-Verlag}, - year = {2007}, - author = {Meinard Muller}, - pages = {64-84}, - address = {London} -} - @misc{Didan:2015, author = {Kamel Didan}, title = {{MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global diff --git a/inst/lucc_MT/assessment.RData b/inst/lucc_MT/assessment.RData deleted file mode 100644 index f32243f..0000000 Binary files a/inst/lucc_MT/assessment.RData and /dev/null differ diff --git a/inst/lucc_MT/cross_validation.RData b/inst/lucc_MT/cross_validation.RData deleted file mode 100644 index 13044dd..0000000 Binary files a/inst/lucc_MT/cross_validation.RData and /dev/null differ diff --git a/inst/lucc_MT/data/blue.tif b/inst/lucc_MT/data/blue.tif deleted file mode 100644 index a0ad167..0000000 Binary files a/inst/lucc_MT/data/blue.tif and /dev/null differ diff --git a/inst/lucc_MT/data/doy.tif b/inst/lucc_MT/data/doy.tif deleted file mode 100644 index 2ad04d5..0000000 Binary files a/inst/lucc_MT/data/doy.tif and /dev/null differ diff --git a/inst/lucc_MT/data/evi.tif b/inst/lucc_MT/data/evi.tif deleted file mode 100644 index ded4628..0000000 Binary files a/inst/lucc_MT/data/evi.tif and /dev/null differ diff --git a/inst/lucc_MT/data/mir.tif b/inst/lucc_MT/data/mir.tif deleted file mode 100644 index bc7e021..0000000 Binary files a/inst/lucc_MT/data/mir.tif and /dev/null differ diff --git a/inst/lucc_MT/data/ndvi.tif b/inst/lucc_MT/data/ndvi.tif deleted file mode 100644 index bf88dae..0000000 Binary files a/inst/lucc_MT/data/ndvi.tif and /dev/null differ diff --git a/inst/lucc_MT/data/nir.tif b/inst/lucc_MT/data/nir.tif deleted file mode 100644 index 4715ff5..0000000 Binary files a/inst/lucc_MT/data/nir.tif and /dev/null differ diff --git a/inst/lucc_MT/data/red.tif b/inst/lucc_MT/data/red.tif deleted file mode 100644 index 6ab2386..0000000 Binary files a/inst/lucc_MT/data/red.tif and /dev/null differ diff --git a/inst/lucc_MT/data/samples.csv b/inst/lucc_MT/data/samples.csv deleted file mode 100644 index 9e9706e..0000000 --- a/inst/lucc_MT/data/samples.csv +++ /dev/null @@ -1,604 +0,0 @@ -"longitude","latitude","from","to","label" --55.9881860661,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9911845738,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9860559013,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.985621884,-12.0343749989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9847540916,-12.0302083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9839257365,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.976707366,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9809276839,-12.0322916656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9928805432,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9886201806,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9813615874,-12.0343749989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9869241782,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9877131401,-12.0239583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9894886521,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9834917357,-12.0343749989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9420173552,-11.9885416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9817955717,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9907503619,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9933147716,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.989054376,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9838866222,-12.0260416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9864899993,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.984359818,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9843203165,-12.0281249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9954449694,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9958792951,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9855830744,-12.0239583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9916188664,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9860167045,-12.0260416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9788375803,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9821902178,-12.0281249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.944147141,-11.9885416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9800601191,-12.0281249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.98175654,-12.0260416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9864504153,-12.0281249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9784033869,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9937490807,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9809677947,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9732749126,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9647542535,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9711447478,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9668844182,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9656222,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9677523978,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9634920021,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9609278238,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9694485489,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9754050774,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9720127934,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9698825956,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9707108461,-12.0343749989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9728409944,-12.0343749989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9613618043,-12.0406249989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9737089115,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.969014583,-12.0364583323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9651881863,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.963058005,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.970316723,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9715787302,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9681865086,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9673183676,-12.0385416656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9617958656,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9660562943,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9596656513,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9745372886,-12.0322916656,"2011-09-01","2012-09-01","Cotton-fallow" --55.9749711426,-12.0343749989,"2011-09-01","2012-09-01","Cotton-fallow" --55.9766674204,-12.0322916656,"2011-09-01","2012-09-01","Cotton-fallow" --55.96392608,-12.0427083323,"2011-09-01","2012-09-01","Cotton-fallow" --55.9415804282,-12.0382604187,"2008-09-01","2009-09-01","Forest" --55.9254548352,-12.0411572318,"2007-09-01","2008-09-01","Forest" --55.9254548352,-12.0411572318,"2009-09-01","2010-09-01","Forest" --55.9254548352,-12.0411572318,"2010-09-01","2011-09-01","Forest" --55.9254548352,-12.0411572318,"2011-09-01","2012-09-01","Forest" --55.9192749672,-12.0337220781,"2007-09-01","2008-09-01","Forest" --55.9192749672,-12.0337220781,"2008-09-01","2009-09-01","Forest" --55.9192749672,-12.0337220781,"2009-09-01","2010-09-01","Forest" --55.9192749672,-12.0337220781,"2010-09-01","2011-09-01","Forest" --55.9192749672,-12.0337220781,"2011-09-01","2012-09-01","Forest" --55.9192749672,-12.0337220781,"2012-09-01","2013-09-01","Forest" --55.9169575167,-12.0332392759,"2007-09-01","2008-09-01","Forest" --55.9169575167,-12.0332392759,"2008-09-01","2009-09-01","Forest" --55.9169575167,-12.0332392759,"2009-09-01","2010-09-01","Forest" --55.9169575167,-12.0332392759,"2010-09-01","2011-09-01","Forest" --55.9169575167,-12.0332392759,"2012-09-01","2013-09-01","Forest" --55.9168609563,-12.035460166,"2007-09-01","2008-09-01","Forest" --55.9168609563,-12.035460166,"2008-09-01","2009-09-01","Forest" --55.9168609563,-12.035460166,"2009-09-01","2010-09-01","Forest" --55.9400354612,-12.0368120121,"2007-09-01","2008-09-01","Forest" --55.9400354612,-12.0368120121,"2008-09-01","2009-09-01","Forest" --55.9400354612,-12.0368120121,"2009-09-01","2010-09-01","Forest" --55.9400354612,-12.0368120121,"2010-09-01","2011-09-01","Forest" --55.9400354612,-12.0368120121,"2011-09-01","2012-09-01","Forest" --55.9169575167,-12.0332392759,"2011-09-01","2012-09-01","Forest" --55.9400354612,-12.0368120121,"2012-09-01","2013-09-01","Forest" --55.9254548352,-12.0411572318,"2008-09-01","2009-09-01","Forest" --55.9181162419,-12.0384535395,"2007-09-01","2008-09-01","Forest" --55.9181162419,-12.0384535395,"2009-09-01","2010-09-01","Forest" --55.9181162419,-12.0384535395,"2010-09-01","2011-09-01","Forest" --55.9181162419,-12.0384535395,"2011-09-01","2012-09-01","Forest" --55.9181162419,-12.0384535395,"2012-09-01","2013-09-01","Forest" --55.9393595381,-12.0387432209,"2007-09-01","2008-09-01","Forest" --55.9393595381,-12.0387432209,"2008-09-01","2009-09-01","Forest" --55.9393595381,-12.0387432209,"2009-09-01","2010-09-01","Forest" --55.9393595381,-12.0387432209,"2010-09-01","2011-09-01","Forest" --55.9393595381,-12.0387432209,"2011-09-01","2012-09-01","Forest" --55.9393595381,-12.0387432209,"2012-09-01","2013-09-01","Forest" --55.9403251425,-12.0421228362,"2007-09-01","2008-09-01","Forest" --55.9403251425,-12.0421228362,"2008-09-01","2009-09-01","Forest" --55.9403251425,-12.0421228362,"2009-09-01","2010-09-01","Forest" --55.9403251425,-12.0421228362,"2010-09-01","2011-09-01","Forest" --55.9403251425,-12.0421228362,"2011-09-01","2012-09-01","Forest" --55.9412907469,-12.0319839902,"2008-09-01","2009-09-01","Forest" --55.9435116369,-12.0408675505,"2007-09-01","2008-09-01","Forest" --55.9435116369,-12.0408675505,"2008-09-01","2009-09-01","Forest" --55.9435116369,-12.0408675505,"2009-09-01","2010-09-01","Forest" --55.9435116369,-12.0408675505,"2010-09-01","2011-09-01","Forest" --55.9435116369,-12.0408675505,"2011-09-01","2012-09-01","Forest" --55.9435116369,-12.0408675505,"2012-09-01","2013-09-01","Forest" --55.9422563513,-12.0409641109,"2007-09-01","2008-09-01","Forest" --55.9422563513,-12.0409641109,"2008-09-01","2009-09-01","Forest" --55.9422563513,-12.0409641109,"2009-09-01","2010-09-01","Forest" --55.9422563513,-12.0409641109,"2010-09-01","2011-09-01","Forest" --55.9199508902,-12.0360395286,"2007-09-01","2008-09-01","Forest" --55.9181162419,-12.0384535395,"2008-09-01","2009-09-01","Forest" --55.9415804282,-12.0382604187,"2007-09-01","2008-09-01","Forest" --55.9415804282,-12.0382604187,"2009-09-01","2010-09-01","Forest" --55.9415804282,-12.0382604187,"2010-09-01","2011-09-01","Forest" --55.9415804282,-12.0382604187,"2011-09-01","2012-09-01","Forest" --55.9415804282,-12.0382604187,"2012-09-01","2013-09-01","Forest" --55.9168609563,-12.035460166,"2010-09-01","2011-09-01","Forest" --55.9397457799,-12.0403847483,"2007-09-01","2008-09-01","Forest" --55.9397457799,-12.0403847483,"2008-09-01","2009-09-01","Forest" --55.9397457799,-12.0403847483,"2009-09-01","2010-09-01","Forest" --55.9397457799,-12.0403847483,"2010-09-01","2011-09-01","Forest" --55.9397457799,-12.0403847483,"2011-09-01","2012-09-01","Forest" --55.9397457799,-12.0403847483,"2012-09-01","2013-09-01","Forest" --55.9179231211,-12.0401916274,"2009-09-01","2010-09-01","Forest" --55.9220752199,-12.0413503527,"2007-09-01","2008-09-01","Forest" --55.9220752199,-12.0413503527,"2008-09-01","2009-09-01","Forest" --55.9220752199,-12.0413503527,"2009-09-01","2010-09-01","Forest" --55.9220752199,-12.0413503527,"2010-09-01","2011-09-01","Forest" --55.9220752199,-12.0413503527,"2011-09-01","2012-09-01","Forest" --55.9220752199,-12.0413503527,"2012-09-01","2013-09-01","Forest" --55.9412907469,-12.0319839902,"2007-09-01","2008-09-01","Forest" --55.9412907469,-12.0319839902,"2009-09-01","2010-09-01","Forest" --55.9254548352,-12.0411572318,"2012-09-01","2013-09-01","Forest" --55.9179231211,-12.0401916274,"2012-09-01","2013-09-01","Forest" --55.9168609563,-12.035460166,"2011-09-01","2012-09-01","Forest" --55.9168609563,-12.035460166,"2012-09-01","2013-09-01","Forest" --55.9179231211,-12.0401916274,"2007-09-01","2008-09-01","Forest" --55.9179231211,-12.0401916274,"2008-09-01","2009-09-01","Forest" --55.9403251425,-12.0421228362,"2012-09-01","2013-09-01","Forest" --55.9179231211,-12.0401916274,"2010-09-01","2011-09-01","Forest" --55.9179231211,-12.0401916274,"2011-09-01","2012-09-01","Forest" --55.9382973733,-12.0419297153,"2007-09-01","2008-09-01","Forest" --55.9382973733,-12.0419297153,"2008-09-01","2009-09-01","Forest" --55.9382973733,-12.0419297153,"2009-09-01","2010-09-01","Forest" --55.9382973733,-12.0419297153,"2010-09-01","2011-09-01","Forest" --55.9382973733,-12.0419297153,"2011-09-01","2012-09-01","Forest" --55.9382973733,-12.0419297153,"2012-09-01","2013-09-01","Forest" --55.9412907469,-12.0319839902,"2010-09-01","2011-09-01","Forest" --55.9412907469,-12.0319839902,"2011-09-01","2012-09-01","Forest" --55.9412907469,-12.0319839902,"2012-09-01","2013-09-01","Forest" --55.9381042525,-12.0403847483,"2007-09-01","2008-09-01","Forest" --55.9381042525,-12.0403847483,"2008-09-01","2009-09-01","Forest" --55.9381042525,-12.0403847483,"2009-09-01","2010-09-01","Forest" --55.9381042525,-12.0403847483,"2010-09-01","2011-09-01","Forest" --55.9381042525,-12.0403847483,"2011-09-01","2012-09-01","Forest" --55.9381042525,-12.0403847483,"2012-09-01","2013-09-01","Forest" --55.9423529117,-12.0352670451,"2008-09-01","2009-09-01","Forest" --55.9423529117,-12.0352670451,"2009-09-01","2010-09-01","Forest" --55.9409045051,-12.0352670451,"2007-09-01","2008-09-01","Forest" --55.9409045051,-12.0352670451,"2008-09-01","2009-09-01","Forest" --55.9409045051,-12.0352670451,"2009-09-01","2010-09-01","Forest" --55.9409045051,-12.0352670451,"2010-09-01","2011-09-01","Forest" --55.9409045051,-12.0352670451,"2011-09-01","2012-09-01","Forest" --55.9409045051,-12.0352670451,"2012-09-01","2013-09-01","Forest" --55.9199508902,-12.0360395286,"2010-09-01","2011-09-01","Forest" --55.9199508902,-12.0360395286,"2011-09-01","2012-09-01","Forest" --55.9199508902,-12.0360395286,"2012-09-01","2013-09-01","Forest" --55.9212061759,-12.0384535395,"2008-09-01","2009-09-01","Forest" --55.9423529117,-12.0352670451,"2007-09-01","2008-09-01","Forest" --55.9422563513,-12.0409641109,"2012-09-01","2013-09-01","Forest" --55.9212061759,-12.0384535395,"2007-09-01","2008-09-01","Forest" --55.9212061759,-12.0384535395,"2009-09-01","2010-09-01","Forest" --55.939552659,-12.0349773638,"2007-09-01","2008-09-01","Forest" --55.939552659,-12.0349773638,"2008-09-01","2009-09-01","Forest" --55.939552659,-12.0349773638,"2009-09-01","2010-09-01","Forest" --55.939552659,-12.0349773638,"2010-09-01","2011-09-01","Forest" --55.939552659,-12.0349773638,"2011-09-01","2012-09-01","Forest" --55.939552659,-12.0349773638,"2012-09-01","2013-09-01","Forest" --55.9365592855,-12.0406744296,"2007-09-01","2008-09-01","Forest" --55.9365592855,-12.0406744296,"2008-09-01","2009-09-01","Forest" --55.9365592855,-12.0406744296,"2010-09-01","2011-09-01","Forest" --55.9365592855,-12.0406744296,"2011-09-01","2012-09-01","Forest" --55.9365592855,-12.0406744296,"2012-09-01","2013-09-01","Forest" --55.9212061759,-12.0384535395,"2010-09-01","2011-09-01","Forest" --55.9212061759,-12.0384535395,"2012-09-01","2013-09-01","Forest" --55.9422563513,-12.0409641109,"2011-09-01","2012-09-01","Forest" --55.9423529117,-12.0352670451,"2010-09-01","2011-09-01","Forest" --55.9423529117,-12.0352670451,"2011-09-01","2012-09-01","Forest" --55.9423529117,-12.0352670451,"2012-09-01","2013-09-01","Forest" --55.9199508902,-12.0360395286,"2008-09-01","2009-09-01","Forest" --55.9212061759,-12.0384535395,"2011-09-01","2012-09-01","Forest" --55.9365592855,-12.0406744296,"2009-09-01","2010-09-01","Forest" --55.9199508902,-12.0360395286,"2009-09-01","2010-09-01","Forest" --55.9646793633,-12.0156249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.976629328,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9685425484,-12.0239583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9715397241,-12.0281249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9693725546,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9625493635,-12.0156249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9711061287,-12.0260416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9689760465,-12.0260416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9630193291,-12.0281249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.977459362,-12.0156249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9728026799,-12.0239583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9706726141,-12.0239583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9672795266,-12.0281249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9613228239,-12.0302083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.963380371,-12.0093749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9595533541,-12.0114583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9574233873,-12.0114583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9492708486,-11.9927083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9770628114,-12.0239583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.973236211,-12.0260416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9570625934,-12.0302083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.966809363,-12.0156249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.962947535,-12.0072916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9795893618,-12.0156249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9736325871,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9740658926,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9668459643,-12.0260416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.964715882,-12.0260416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9694096253,-12.0281249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9651494278,-12.0281249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9608892303,-12.0281249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9672425383,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9651125221,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9655457616,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9548955977,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9574588847,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9646446971,-12.0052083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9616833209,-12.0114583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.955726352,-12.0135416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9518993646,-12.0156249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9488387402,-11.9906249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9787593773,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9570256305,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9378553357,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9590871297,-11.9989583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9459423756,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9399853684,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9361583922,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9591927086,-12.0302083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9629825058,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9608524895,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9634157288,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9578922197,-12.0239583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9612504207,-12.0093749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9578563353,-12.0135416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9438123593,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9442454339,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.950600669,-12.0093749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9530983446,-11.9906249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.94714103,-11.9927083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9421154012,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9382884415,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9404184907,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9591556633,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9497693648,-12.0156249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9480723919,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9514006672,-11.9927083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9595196266,-12.0010416656,"2011-09-01","2012-09-01","Soybean-cotton" --55.959588934,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.94890352,-12.0114583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9587224732,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.956592457,-12.0177083323,"2011-09-01","2012-09-01","Soybean-cotton" --55.961285696,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9638490325,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9617189832,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9509685424,-11.9906249989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9676757943,-12.0197916656,"2011-09-01","2012-09-01","Soybean-cotton" --55.9553288355,-12.0218749989,"2011-09-01","2012-09-01","Soybean-cotton" --55.9600222854,-12.0239583323,"2011-09-01","2012-09-01","Soybean-cotton" --55.9786559709,-11.9906249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9766674204,-12.0322916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9745372886,-12.0322916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9707108461,-12.0343749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9748286723,-11.9927083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9710013886,-11.9947916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9718663378,-11.9989583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9795893618,-12.0156249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9736325871,-12.0177083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9740658926,-12.0197916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9787593773,-12.0218749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9864504153,-12.0281249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9843203165,-12.0281249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9800601191,-12.0281249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9928407115,-12.0281249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9915783427,-12.0322916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9754050774,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9765261687,-11.9906249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9855830744,-12.0239583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9860167045,-12.0260416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9838866222,-12.0260416656,"2010-09-01","2011-09-01","Soybean-maize" --55.98175654,-12.0260416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9821902178,-12.0281249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9847540916,-12.0302083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9839257365,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9928805432,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.989054376,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9803534987,-11.9885416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9497352753,-12.0052083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9522978657,-12.0072916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9807857731,-11.9906249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9790883096,-11.9927083323,"2010-09-01","2011-09-01","Soybean-maize" --55.976958491,-11.9927083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9726988537,-11.9927083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9752610588,-11.9947916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9735636744,-11.9968749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9714338228,-11.9968749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9782559418,-11.9989583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9735972044,-12.0072916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9719001724,-12.0093749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9877131401,-12.0239583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9907106128,-12.0281249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9932745526,-12.0302083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9911444374,-12.0302083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9894482109,-12.0322916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9809276839,-12.0322916656,"2010-09-01","2011-09-01","Soybean-maize" --55.985621884,-12.0343749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9834917357,-12.0343749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9881860661,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9860559013,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9907503619,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9886201806,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.984359818,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9954449694,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9911845738,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9869241782,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9916188664,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9786885868,-12.0010416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9774242015,-12.0052083323,"2010-09-01","2011-09-01","Soybean-maize" --55.977459362,-12.0156249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9693039713,-11.9968749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9697364698,-11.9989583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9676066018,-11.9989583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9701690489,-12.0010416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9680391644,-12.0010416656,"2010-09-01","2011-09-01","Soybean-maize" --55.96590928,-12.0010416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9820830806,-11.9968749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9748615106,-12.0031249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9757271383,-12.0072916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9958792951,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9937490807,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9749711426,-12.0343749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9728409944,-12.0343749989,"2010-09-01","2011-09-01","Soybean-maize" --55.969014583,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9715787302,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9694485489,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.963058005,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9677523978,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9656222,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.970316723,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9681865086,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9731312237,-11.9947916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9842129321,-11.9968749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9744631219,-12.0114583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9898060848,-12.0135416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9889761219,-12.0197916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9404184907,-12.0218749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9813615874,-12.0343749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9817955717,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9864899993,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9933147716,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9894886521,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9784033869,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9809677947,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9788375803,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.976707366,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9920123288,-12.0343749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9518993646,-12.0156249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9769914116,-12.0031249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9731643667,-12.0052083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9714672705,-12.0072916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9740301228,-12.0093749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9497693648,-12.0156249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9438123593,-12.0177083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9378553357,-12.0197916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9399853684,-12.0197916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9514326003,-12.0031249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9480379979,-12.0072916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9480723919,-12.0177083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9459423756,-12.0177083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9442454339,-12.0197916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9684718078,-12.0031249989,"2010-09-01","2011-09-01","Soybean-maize" --55.979953229,-11.9968749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9765587023,-12.0010416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9752942841,-12.0052083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9668844182,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9382884415,-12.0218749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9837803992,-11.9947916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9421154012,-12.0197916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9803858098,-11.9989583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9720127934,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.976629328,-12.0218749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9518651927,-12.0052083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9476053579,-12.0052083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9501679318,-12.0072916656,"2010-09-01","2011-09-01","Soybean-maize" --55.9711447478,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9673183676,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9732749126,-12.0364583323,"2010-09-01","2011-09-01","Soybean-maize" --55.9651881863,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.96392608,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9484707186,-12.0093749989,"2010-09-01","2011-09-01","Soybean-maize" --55.9698825956,-12.0406249989,"2010-09-01","2011-09-01","Soybean-maize" --55.9660562943,-12.0427083323,"2010-09-01","2011-09-01","Soybean-maize" --55.9737089115,-12.0385416656,"2010-09-01","2011-09-01","Soybean-maize" --55.9250412648,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9229113145,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9186514138,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9165214634,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9143915131,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9318637854,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9297338186,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9276038518,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9233439181,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9212139513,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9190839845,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9148240508,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.928036569,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9237766024,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9199493678,-12.0156249989,"2011-09-01","2012-09-01","Soybean-millet" --55.917819368,-12.0156249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9156893682,-12.0156249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9135593684,-12.0156249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9348947306,-12.0260416656,"2012-09-01","2013-09-01","Soybean-millet" --55.9285044838,-12.0260416656,"2012-09-01","2013-09-01","Soybean-millet" --55.9310678479,-12.0281249989,"2012-09-01","2013-09-01","Soybean-millet" --55.9289377491,-12.0281249989,"2012-09-01","2013-09-01","Soybean-millet" --55.9268076504,-12.0281249989,"2012-09-01","2013-09-01","Soybean-millet" --55.9336313256,-12.0302083323,"2012-09-01","2013-09-01","Soybean-millet" --55.9315012103,-12.0302083323,"2012-09-01","2013-09-01","Soybean-millet" --55.9293710951,-12.0302083323,"2012-09-01","2013-09-01","Soybean-millet" --55.9272409798,-12.0302083323,"2012-09-01","2013-09-01","Soybean-millet" --55.9340647852,-12.0322916656,"2012-09-01","2013-09-01","Soybean-millet" --55.9319346535,-12.0322916656,"2012-09-01","2013-09-01","Soybean-millet" --55.9298045217,-12.0322916656,"2012-09-01","2013-09-01","Soybean-millet" --55.9276743899,-12.0322916656,"2012-09-01","2013-09-01","Soybean-millet" --55.9255442582,-12.0322916656,"2012-09-01","2013-09-01","Soybean-millet" --55.9323681773,-12.0343749989,"2012-09-01","2013-09-01","Soybean-millet" --55.930238029,-12.0343749989,"2012-09-01","2013-09-01","Soybean-millet" --55.9281078807,-12.0343749989,"2012-09-01","2013-09-01","Soybean-millet" --55.9259777324,-12.0343749989,"2012-09-01","2013-09-01","Soybean-millet" --55.9285414522,-12.0364583323,"2012-09-01","2013-09-01","Soybean-millet" --55.9264112874,-12.0364583323,"2012-09-01","2013-09-01","Soybean-millet" --55.9212480474,-12.0218749989,"2012-09-01","2013-09-01","Soybean-millet" --55.9195510362,-12.0239583323,"2012-09-01","2013-09-01","Soybean-millet" --55.9174209705,-12.0239583323,"2012-09-01","2013-09-01","Soybean-millet" --55.9199841548,-12.0260416656,"2012-09-01","2013-09-01","Soybean-millet" --55.9178540726,-12.0260416656,"2012-09-01","2013-09-01","Soybean-millet" --55.9157239904,-12.0260416656,"2012-09-01","2013-09-01","Soybean-millet" --55.9182872554,-12.0281249989,"2012-09-01","2013-09-01","Soybean-millet" --55.9161571566,-12.0281249989,"2012-09-01","2013-09-01","Soybean-millet" --55.9306345661,-12.0260416656,"2012-09-01","2013-09-01","Soybean-millet" --55.9203488576,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.913959056,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9271712152,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9199493678,-12.0156249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9135593684,-12.0156249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9203821804,-12.0177083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9182521641,-12.0177083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9161221478,-12.0177083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9139921315,-12.0177083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9211817061,-12.0010416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9190518216,-12.0010416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9190518216,-12.0010416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9169219372,-12.0010416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9169219372,-12.0010416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9147920527,-12.0010416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9126621682,-12.0010416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9343933929,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9322634919,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9322634919,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.930133591,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9216139873,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9194840863,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9194840863,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9173541854,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9173541854,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9152242845,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9152242845,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9130943835,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9130943835,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9326959361,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9326959361,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9305660186,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9284361012,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9241762664,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.922046349,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9199164316,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9177865142,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9156565968,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9135266794,-12.0052083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9331284608,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9309985269,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9288685931,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9267386592,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9267386592,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9246087253,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9224787914,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9203488576,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9182189237,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9182189237,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9160889898,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9335610662,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9335610662,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9314311159,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9314311159,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9293011655,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9250412648,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9229113145,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9207813641,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9165214634,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9143915131,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9339937522,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9297338186,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9276038518,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.925473885,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9233439181,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9212139513,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9190839845,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9169540177,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9301665523,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9301665523,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.928036569,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9259065857,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9216466191,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9216466191,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9195166358,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9195166358,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9173866525,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9152566692,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9152566692,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9263393671,-12.0156249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9263393671,-12.0156249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9242093673,-12.0156249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9242093673,-12.0156249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9220793676,-12.0156249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9220793676,-12.0156249989,"2011-09-01","2012-09-01","Soybean-millet" --55.917819368,-12.0156249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9156893682,-12.0156249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9203821804,-12.0177083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9182521641,-12.0177083323,"2011-09-01","2012-09-01","Soybean-millet" --55.9161221478,-12.0177083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9139921315,-12.0177083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9165550081,-12.0197916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9165550081,-12.0197916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9144249753,-12.0197916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9323314307,-12.0239583323,"2012-09-01","2013-09-01","Soybean-millet" --55.9302013649,-12.0239583323,"2012-09-01","2013-09-01","Soybean-millet" --55.9327646483,-12.0260416656,"2012-09-01","2013-09-01","Soybean-millet" --55.9353280454,-12.0281249989,"2012-09-01","2013-09-01","Soybean-millet" --55.9331979466,-12.0281249989,"2012-09-01","2013-09-01","Soybean-millet" --55.9144249753,-12.0197916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9224787914,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.913959056,-12.0072916656,"2011-09-01","2012-09-01","Soybean-millet" --55.9271712152,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9160889898,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9207813641,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.925473885,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9259065857,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9237766024,-12.0135416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9173866525,-12.0135416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9246087253,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9211817061,-12.0010416656,"2011-09-01","2012-09-01","Soybean-millet" --55.9147920527,-12.0010416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9126621682,-12.0010416656,"2010-09-01","2011-09-01","Soybean-millet" --55.9309985269,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9339937522,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9169540177,-12.0114583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9148240508,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9237438882,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9135266794,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9331284608,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.9186514138,-12.0093749989,"2011-09-01","2012-09-01","Soybean-millet" --55.9343933929,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9284361012,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9199164316,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9293011655,-12.0093749989,"2010-09-01","2011-09-01","Soybean-millet" --55.9318637854,-12.0114583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9241762664,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9288685931,-12.0072916656,"2010-09-01","2011-09-01","Soybean-millet" --55.922046349,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9156565968,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.9177865142,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" --55.930133591,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9186196375,-11.9989583323,"2011-09-01","2012-09-01","Soybean-millet" --55.9237438882,-12.0031249989,"2010-09-01","2011-09-01","Soybean-millet" --55.9216139873,-12.0031249989,"2011-09-01","2012-09-01","Soybean-millet" --55.9186196375,-11.9989583323,"2010-09-01","2011-09-01","Soybean-millet" --55.9305660186,-12.0052083323,"2010-09-01","2011-09-01","Soybean-millet" diff --git a/inst/lucc_MT/data/samples_projection b/inst/lucc_MT/data/samples_projection deleted file mode 100644 index 4128610..0000000 --- a/inst/lucc_MT/data/samples_projection +++ /dev/null @@ -1 +0,0 @@ -"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" \ No newline at end of file diff --git a/inst/lucc_MT/data/timeline b/inst/lucc_MT/data/timeline deleted file mode 100644 index 43db978..0000000 --- a/inst/lucc_MT/data/timeline +++ /dev/null @@ -1,137 +0,0 @@ -2007-09-14 -2007-09-30 -2007-10-16 -2007-11-01 -2007-11-17 -2007-12-03 -2007-12-19 -2008-01-01 -2008-01-17 -2008-02-02 -2008-02-18 -2008-03-05 -2008-03-21 -2008-04-06 -2008-04-22 -2008-05-08 -2008-05-24 -2008-06-09 -2008-06-25 -2008-07-11 -2008-07-27 -2008-08-12 -2008-08-28 -2008-09-13 -2008-09-29 -2008-10-15 -2008-10-31 -2008-11-16 -2008-12-02 -2008-12-18 -2009-01-01 -2009-01-17 -2009-02-02 -2009-02-18 -2009-03-06 -2009-03-22 -2009-04-07 -2009-04-23 -2009-05-09 -2009-05-25 -2009-06-10 -2009-06-26 -2009-07-12 -2009-07-28 -2009-08-13 -2009-08-29 -2009-09-14 -2009-09-30 -2009-10-16 -2009-11-01 -2009-11-17 -2009-12-03 -2009-12-19 -2010-01-01 -2010-01-17 -2010-02-02 -2010-02-18 -2010-03-06 -2010-03-22 -2010-04-07 -2010-04-23 -2010-05-09 -2010-05-25 -2010-06-10 -2010-06-26 -2010-07-12 -2010-07-28 -2010-08-13 -2010-08-29 -2010-09-14 -2010-09-30 -2010-10-16 -2010-11-01 -2010-11-17 -2010-12-03 -2010-12-19 -2011-01-01 -2011-01-17 -2011-02-02 -2011-02-18 -2011-03-06 -2011-03-22 -2011-04-07 -2011-04-23 -2011-05-09 -2011-05-25 -2011-06-10 -2011-06-26 -2011-07-12 -2011-07-28 -2011-08-13 -2011-08-29 -2011-09-14 -2011-09-30 -2011-10-16 -2011-11-01 -2011-11-17 -2011-12-03 -2011-12-19 -2012-01-01 -2012-01-17 -2012-02-02 -2012-02-18 -2012-03-05 -2012-03-21 -2012-04-06 -2012-04-22 -2012-05-08 -2012-05-24 -2012-06-09 -2012-06-25 -2012-07-11 -2012-07-27 -2012-08-12 -2012-08-28 -2012-09-13 -2012-09-29 -2012-10-15 -2012-10-31 -2012-11-16 -2012-12-02 -2012-12-18 -2013-01-01 -2013-01-17 -2013-02-02 -2013-02-18 -2013-03-06 -2013-03-22 -2013-04-07 -2013-04-23 -2013-05-09 -2013-05-25 -2013-06-10 -2013-06-26 -2013-07-12 -2013-08-13 -2013-08-29 diff --git a/inst/lucc_MT/patterns/araguaia.csv b/inst/lucc_MT/patterns/araguaia.csv deleted file mode 100644 index 5963c44..0000000 --- a/inst/lucc_MT/patterns/araguaia.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2003-08-29,0.14739560800314064,0.03375004786473038,0.16061986869675754,0.06745353123598899,0.1744863791459053,0.3993501997794543 -2003-09-06,0.157352603470522,0.03283999361940011,0.15760991443727204,0.06761722512983084,0.16909019940948328,0.3908705414590904 -2003-09-14,0.166365304317564,0.03228892645263741,0.15600700391187805,0.06796599367343695,0.165777889375191,0.38518743214636003 -2003-09-22,0.1734582725724762,0.03242836176730144,0.15714924688969617,0.0686589554917648,0.16657126021120217,0.3850443219951707 -2003-09-30,0.17779525989035463,0.03335403459044552,0.16163057855506222,0.06967886328863515,0.1725760599478699,0.3920953908862747 -2003-10-08,0.17870440189667308,0.035026416772141515,0.1695738045221412,0.07091837574100046,0.18426725424949741,0.40719887006264455 -2003-10-16,0.17608685024863752,0.037109866694218584,0.17995824907191404,0.07210093281851691,0.20047225236465455,0.4290373382994531 -2003-10-24,0.17013006507800504,0.03918608169908758,0.19136827279395896,0.07292006549863311,0.21938179028890187,0.4554032562772271 -2003-11-01,0.16154606808416627,0.04087412505004342,0.20217471990296168,0.07317250127345427,0.23855967868985076,0.4830280792289939 -2003-11-09,0.1512338960461685,0.041873144066829654,0.2108630020902044,0.07274613984144263,0.25555856025201784,0.5085005720149143 -2003-11-17,0.14021953314069602,0.042204927574101075,0.21676337555121764,0.07180918689659588,0.2687784419502823,0.5292618793865862 -2003-11-25,0.12950402634263466,0.042033045186983066,0.21965148893933137,0.07063862685959069,0.27715768558976833,0.5433908207714471 -2003-12-03,0.11978880032162414,0.041723439864045066,0.22018250307003648,0.06961945565349638,0.2809486495633471,0.5507410945514427 -2003-12-11,0.1115861441706858,0.04166675366768011,0.2192806618282189,0.06912008870845215,0.28092144850123757,0.5519489713213668 -2003-12-19,0.10494851334437523,0.04207862834271925,0.2178721202425613,0.06927526452532279,0.27839458457848393,0.5487837440398329 -2003-12-27,0.09972902208799964,0.043049174214916615,0.2167413269286631,0.07007865292829155,0.2747628298136268,0.5433268428879617 -2004-01-04,0.09550298845254201,0.044328907077615186,0.2160014402972491,0.07119573052188678,0.27098913190402385,0.5374967692659045 -2004-01-12,0.0917642532782145,0.04548986091984521,0.21530166748080942,0.07213989571118518,0.26756406342199635,0.532740288074907 -2004-01-20,0.08801521799953277,0.04598525039890601,0.21387231158483536,0.07234351303425546,0.26442081066111617,0.5298077556922617 -2004-01-28,0.08391386876877482,0.04540769401237883,0.2110775568888287,0.07139973233629925,0.2613606440408012,0.5290190919525677 -2004-02-05,0.07937434015443924,0.04360095182683821,0.2065748695420763,0.06918180451013763,0.25806176882809717,0.5301001495579897 -2004-02-13,0.07450057344826219,0.04065092921378689,0.20043627467576508,0.06580936955217019,0.254378579397443,0.532648027491538 -2004-02-21,0.06966689375763725,0.03700022154188866,0.19336685892795216,0.06176098048712692,0.25047533797265775,0.5361452241054652 -2004-02-29,0.06531495503354413,0.033171318896810895,0.18625024908997712,0.057588566654946514,0.2466781292493536,0.5401556138668995 -2004-03-08,0.061951413708870334,0.029757206808600157,0.18015518473798128,0.053903459859318714,0.24352563705864286,0.5443950725206345 -2004-03-16,0.05998092421547459,0.027191389840628393,0.17587807624881563,0.0511507298730356,0.2414642349843929,0.5486610116184384 -2004-03-24,0.05962786580424056,0.025631940616367834,0.17373382655777153,0.04949060628896217,0.2407514474109368,0.5528600117231096 -2004-04-01,0.060934981320364,0.02502931407661026,0.17364000799499646,0.04887022143567699,0.2414207282822007,0.5569214944244473 -2004-04-09,0.0636749767468329,0.025038481719891915,0.17494100053886008,0.04893365113014275,0.24316838413993908,0.5607952372878434 -2004-04-17,0.06750530445765626,0.025283589122385597,0.17689735557291725,0.04928881980244361,0.24560086162138384,0.5644413565528826 -2004-04-25,0.07193885743704981,0.025388340767777292,0.1787289661726774,0.04953572252460384,0.2482334205190675,0.5678333279597554 -2004-05-03,0.07652401889659972,0.025123786950867012,0.17990969657064151,0.04942245547914148,0.25067122875850845,0.5709203051114895 -2004-05-11,0.08089285332226885,0.024503818559083825,0.1803383460858965,0.048942656930490394,0.25268186615664284,0.5736069511722656 -2004-05-19,0.0848243279340348,0.023669082108743834,0.18013614753072504,0.04823050089515993,0.25410666903781187,0.5756679601621144 -2004-05-27,0.08832644259444705,0.022927565149191385,0.17971794456725929,0.04760996684017513,0.25487986252944683,0.5766693645222478 -2004-06-04,0.09154106927142404,0.022548973710341917,0.17942378032528716,0.04738633769974419,0.2548682601400938,0.5760070091751255 -2004-06-12,0.09479248965629977,0.022711410984285214,0.17942076954594297,0.0478017787147809,0.2538124568906927,0.5728324226775636 -2004-06-20,0.0984198696396643,0.023462461310021228,0.17966074341709878,0.04897223769212221,0.2513813984720936,0.566315060103165 -2004-06-28,0.1027556609381923,0.024650478844696903,0.1797715878616486,0.05081646083978328,0.24714983138853244,0.5556991999607237 -2004-07-06,0.10802579349893218,0.026074293497807193,0.17935598215516363,0.05318057685840513,0.24081307365893567,0.5405610759148308 -2004-07-14,0.11428102096408499,0.02748518394086073,0.1780312768437244,0.055824490845978,0.23228331147694947,0.5210151521251677 -2004-07-22,0.12144745459042529,0.028710461348427935,0.1756200556110752,0.058551267752860674,0.22171866188868067,0.4975866164133348 -2004-07-30,0.12927541277095625,0.029710329749120672,0.17227961390656615,0.06124927716924343,0.20964499520251292,0.47138565636258273 -2004-08-07,0.13750509358429996,0.03051311873502884,0.16828041463943977,0.06387106940984308,0.19663102641092445,0.4435262266159805 diff --git a/inst/lucc_MT/patterns/campo_cerrado.csv b/inst/lucc_MT/patterns/campo_cerrado.csv deleted file mode 100644 index c7525fa..0000000 --- a/inst/lucc_MT/patterns/campo_cerrado.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2003-09-14,0.21413687134018405,0.048672053929892906,0.23336044543861054,0.1139625352737606,0.19341421281796733,0.3526594875417768 -2003-09-22,0.21362667186245604,0.04973961018982551,0.23980984777402398,0.11416749095994455,0.20381472473003215,0.36598947241844165 -2003-09-30,0.2125613235486052,0.050830033007930805,0.24634028525910143,0.11428666631525122,0.21461004873753484,0.3796537727110462 -2003-10-08,0.21039223926835268,0.05196241889045884,0.2530048478727791,0.11422943037340932,0.22615806302942137,0.3939635085287554 -2003-10-16,0.2067682713656191,0.05313141232071317,0.2597020779348156,0.11390835980778735,0.23852587643222545,0.4090211893469124 -2003-10-24,0.2014899247357263,0.05431796822255233,0.2662563789206229,0.11324866649057108,0.251607884563542,0.424792240154226 -2003-11-01,0.19478702532232584,0.0554738602697076,0.27236634053893816,0.11223053692047544,0.2649049100078465,0.4408972494582877 -2003-11-09,0.18707552357501578,0.056541963536654825,0.27771477745723155,0.11086720491950054,0.277789123924952,0.45682005715733787 -2003-11-17,0.17903023010487965,0.057465267520634286,0.28209174181575386,0.10924857335536704,0.2895936307710538,0.4718967720582977 -2003-11-25,0.17137898957685568,0.058192682621602126,0.2853732750293922,0.10749557016255464,0.29971440579030567,0.48545835998629744 -2003-12-03,0.16472711634690618,0.058700782496924786,0.2876856916943161,0.10576485920332387,0.30789931016585986,0.49702419413127974 -2003-12-11,0.1595593616033224,0.05897937017289364,0.2892534005433843,0.10421252966775993,0.3140824945531364,0.5062524174322326 -2003-12-19,0.15596900508610945,0.05904067494274494,0.2904020104026127,0.10294036395510205,0.3185685594619343,0.5131838609203984 -2003-12-27,0.15385942443839642,0.05890240101101236,0.2914497352522534,0.10200912466008208,0.32178434505232173,0.5180373353385004 -2004-01-04,0.15280211762151288,0.05857560334278675,0.29253213490403185,0.10136368153295455,0.3241978143747844,0.5213201617658131 -2004-01-12,0.15222856918029656,0.05805811541555839,0.2935892316005838,0.10086588412051482,0.32617741221594343,0.5236568234625288 -2004-01-20,0.15151003626426734,0.05732940985160185,0.2943320219893709,0.10030002603593136,0.32790586374665287,0.5257168044516656 -2004-01-28,0.15015715766522272,0.05636341968646091,0.2944173595151486,0.09945740416636709,0.32941985570672544,0.5280484220868996 -2004-02-05,0.14792173541505171,0.05512824054181848,0.29345261374737486,0.09815449523317178,0.3305468024361334,0.530991348857883 -2004-02-13,0.14477139574729558,0.053605575466021384,0.2911665781217506,0.09628074176968703,0.3310820276562491,0.5346877548971765 -2004-02-21,0.14098966703352853,0.051800770318743736,0.2874974442989778,0.09384220553213161,0.3308032151599817,0.5389879105818669 -2004-02-29,0.13696424332154863,0.049744580822072976,0.2825563356015002,0.09092390629388465,0.32956392795613276,0.5436029180649412 -2004-03-08,0.13320114927659948,0.0475058094052369,0.27670409636373805,0.08772162469178323,0.3273453755696549,0.548062679757431 -2004-03-16,0.1301289398907495,0.04517024199629483,0.2703461976314658,0.08445315331672937,0.3241889495583395,0.551880208134506 -2004-03-24,0.12802798219993405,0.04284527118901761,0.26391906778550667,0.08135386058757288,0.3202133589932805,0.5545709918669182 -2004-04-01,0.12702192710903568,0.04063468156362301,0.25774892570653357,0.07861925507050072,0.3155131130518545,0.5557235863154986 -2004-04-09,0.12700085424316485,0.03863259876031449,0.25197505216508626,0.07637606285935802,0.3101300873350482,0.555047052473798 -2004-04-17,0.12778001398507027,0.03691512940499616,0.24662067470795024,0.07470001707047849,0.30406105972172093,0.5523414996770774 -2004-04-25,0.12909090705385068,0.035530174661156655,0.2415497979702182,0.07359429115441674,0.2972382836831107,0.5475378109502248 -2004-05-03,0.13072814800989752,0.03450789212183736,0.2366423976223756,0.0730460990340305,0.28962266010588356,0.5406349707072269 -2004-05-11,0.13260463925981472,0.03385376074933691,0.23183199038061553,0.07302837543671904,0.2812340308560191,0.5317234711282044 -2004-05-19,0.13473856053119806,0.03356309980154111,0.22716177011713246,0.07352696627536312,0.27218574831653436,0.5209463261488404 -2004-05-27,0.13730125871889132,0.0336184263703944,0.22284615532060714,0.07455100572436149,0.2627354155588773,0.5085189237495332 -2004-06-04,0.14047795437483093,0.03399626021454609,0.21914420227868722,0.07611052565937816,0.25319330566843123,0.4947093450386692 -2004-06-12,0.14445015097051325,0.0346652192318844,0.21635855183617747,0.07821267593515466,0.24393153357447386,0.4798609883857846 -2004-06-20,0.1492993497833582,0.03558792986834457,0.21469439573139848,0.08083506786508972,0.23526489778673165,0.46435244429886124 -2004-06-28,0.1549424642057409,0.03671841981720148,0.21418824244663487,0.08390694443236946,0.22740258843044076,0.4486073283496693 -2004-07-06,0.16119377997026402,0.03800770968265471,0.21473817496418934,0.08732684992908879,0.22043947988233856,0.43302419408322324 -2004-07-14,0.16772658233295445,0.03940338407208868,0.21604275630271502,0.09095109640301441,0.2142994761861164,0.41795154648488675 -2004-07-22,0.1742272114558916,0.04085836946513127,0.21777417069815358,0.0946439420452857,0.20884274821303928,0.4036378235188046 -2004-07-30,0.18042803380413525,0.04233509005663477,0.21959528236540543,0.0982902019916366,0.20385306451610746,0.3901736805790519 -2004-08-07,0.18617760277820739,0.04380753278623462,0.22126961874909412,0.10181818265503631,0.19913652766424364,0.377530319422114 -2004-08-15,0.19150815580995495,0.045267194285248086,0.22272800597390977,0.10522444392198034,0.19454868124450575,0.3655084331883162 -2004-08-23,0.19652345751709577,0.04671297880911857,0.22398367861536944,0.10853332033349779,0.18999413127918283,0.35388173097518827 diff --git a/inst/lucc_MT/patterns/cerradao.csv b/inst/lucc_MT/patterns/cerradao.csv deleted file mode 100644 index 8e172e2..0000000 --- a/inst/lucc_MT/patterns/cerradao.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2014-09-14,0.10486548793474104,0.02457024257294535,0.26807279429696595,0.051281688200426126,0.3927871893841736,0.6721307221483572 -2014-09-22,0.10684493006770285,0.025910671316964035,0.27898527858272343,0.052437852960509276,0.40706022514662993,0.6790892506196199 -2014-09-30,0.1083414479858966,0.027348504408938115,0.2898025361620075,0.053580585930292315,0.4213713259407666,0.686318924363666 -2014-10-08,0.10891567868653959,0.028970691032589077,0.30038598097289004,0.054695249708194435,0.43568064589072425,0.6940187867482321 -2014-10-16,0.10827463976018523,0.0308292333335144,0.31045515077594166,0.05576340350368547,0.4496919005388834,0.7021489094178377 -2014-10-24,0.10629918596164378,0.032941383418332984,0.31968797561678874,0.05677170359459448,0.4629789475797468,0.7104930673707655 -2014-11-01,0.10333571082188467,0.03523528395415447,0.32772522785770236,0.05772612697397507,0.47487000147345704,0.7184257945871666 -2014-11-09,0.09987590007277256,0.03761230620507549,0.334257723440106,0.05864315876327113,0.4846992946348014,0.7252502961321585 -2014-11-17,0.09656265715832044,0.039948232577996005,0.33916638714223313,0.059560208482194225,0.4919925309561932,0.7303219311146 -2014-11-25,0.09397210103156939,0.042130906269543494,0.3424606742803833,0.06051810816771808,0.4964823086593536,0.7331734992532681 -2014-12-03,0.09235712265771621,0.04410528672630091,0.34443846913592396,0.06155429377713925,0.4984451455675524,0.7338819520540345 -2014-12-11,0.09178423886897982,0.04584503584851089,0.3454824939168496,0.06269436741522366,0.4983727740024509,0.7327706613969387 -2014-12-19,0.09193543893887912,0.04737781911563496,0.3460522197661012,0.0639274590651133,0.4970825809843745,0.730590486451149 -2014-12-27,0.09240328673965761,0.04873388603968766,0.3465458736662347,0.06521728774481261,0.4953878983702356,0.7281455353952043 -2015-01-04,0.09275972236745157,0.04991256244793424,0.3471005597125128,0.06646840185060705,0.49385252376802274,0.726109681025306 -2015-01-12,0.09271949320168124,0.05085432789424521,0.34758133250650863,0.06753709096901994,0.49268491860280844,0.7248937968834501 -2015-01-20,0.09217553469118167,0.05143677809984666,0.34759558429411513,0.06823715211081396,0.49173199194230494,0.7246280786973067 -2015-01-28,0.09114468466699603,0.05150593898856005,0.34669553786219875,0.06837361293041361,0.4907096679360583,0.7253579104751444 -2015-02-05,0.0897742229791679,0.05087718828375527,0.34439501319403737,0.06774612122019073,0.48921367743552224,0.7270557408783634 -2015-02-13,0.08815676840615733,0.049416605877531015,0.3404131914018668,0.06620490788459055,0.48703433454083966,0.7298548543186222 -2015-02-21,0.08630758341264651,0.04705628071713094,0.33472301927782966,0.0636640389405455,0.4842071727086752,0.7340843121240666 -2015-03-01,0.08414812844501042,0.043840639263868476,0.3275546539593694,0.06014547406623669,0.48097337718907335,0.7401058772773447 -2015-03-09,0.08149851372656607,0.03994448140204279,0.31942349861346314,0.0557973937162853,0.4777928443994194,0.7482903970391883 -2015-03-17,0.07826487513469702,0.035626716260357905,0.31088366456203126,0.050876286892618125,0.4749902538877112,0.7586067392085017 -2015-03-25,0.07447012657496685,0.03123499882184307,0.3025066495860907,0.045757671032954784,0.4727019328664949,0.770530606568934 -2015-04-02,0.07035491437204436,0.027099726598480294,0.294668099043091,0.04084008489006129,0.4707079635734086,0.7830092349568474 -2015-04-10,0.0664092423277786,0.02352069694656215,0.28750563937869894,0.03653672929773383,0.4683803665826403,0.7944053399993192 -2015-04-18,0.0631726124428744,0.02069332516466418,0.280980485040056,0.03315912439417543,0.46501202988812634,0.8031072053118702 -2015-04-26,0.06121439050963985,0.018688685569265846,0.274865974425115,0.030892191093367022,0.4598486880582345,0.8076045266410968 -2015-05-04,0.06083312405809231,0.01747285808584342,0.2689905543259659,0.029758740222167353,0.45256630202567033,0.8071584364002172 -2015-05-12,0.06199463184367691,0.01689653954938422,0.2632684074394348,0.029596904149479106,0.44336921183765815,0.8019595575372646 -2015-05-20,0.06434007069234274,0.0167826918674132,0.25780694769941614,0.03015534599121155,0.43297329688755015,0.7930014924140224 -2015-05-28,0.06715479971911041,0.016932061929615407,0.25293616527256585,0.031093473049914888,0.42265476539269814,0.782150426701038 -2015-06-05,0.06976808559181613,0.01719935273003564,0.24903170091199645,0.03213774233466209,0.41362155768628883,0.7711738087177733 -2015-06-13,0.07160381368407212,0.017507021685280733,0.24649842150504553,0.03310605123705532,0.4069334787331696,0.761620836868981 -2015-06-21,0.07254291062514376,0.017851501992067333,0.2455225330309055,0.033981563705372114,0.40292538956205515,0.7540048751163599 -2015-06-29,0.07301140270872195,0.01831223539038236,0.24602122114691177,0.03493881042023697,0.40106611961297195,0.7476039643669076 -2015-07-07,0.07382919794995524,0.01899574123354671,0.2476482589727907,0.036253436452359046,0.4001644926960537,0.7407953417022393 -2015-07-15,0.07623324192993793,0.020032909154659873,0.24976744566398235,0.03830316111097278,0.3983232012989928,0.7309989883676266 -2015-07-23,0.08129216583204041,0.021516311610222515,0.25174379194757396,0.04138882683194642,0.39382260845547157,0.7159473490986479 -2015-07-31,0.08982264742476251,0.023489274215923305,0.2529755007632126,0.0457062539410117,0.38523738469301494,0.6938594302839713 -2015-08-08,0.10200601774813328,0.025929171676227096,0.25314519296614413,0.05125702150961883,0.372053969491603,0.6642818668996597 -2015-08-16,0.11728040134336662,0.028738362610052116,0.25227599118918304,0.05781809429461421,0.3548321818109576,0.6283180915606503 -2015-08-24,0.1346931383533912,0.03179047033672381,0.250620219415458,0.0650618916875826,0.3347454617079458,0.5879125327441632 diff --git a/inst/lucc_MT/patterns/cerrado.csv b/inst/lucc_MT/patterns/cerrado.csv deleted file mode 100644 index 6d48938..0000000 --- a/inst/lucc_MT/patterns/cerrado.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2016-09-13,0.18298918491901764,0.03885305916797465,0.24325999767904288,0.09281133369139333,0.25072563802475056,0.4510244423558842 -2016-09-21,0.18138078380950254,0.04013064264594449,0.2514960257198116,0.09257489821019926,0.2654181562105211,0.46711007243052904 -2016-09-29,0.17928460122168352,0.04136783662793099,0.259691593260739,0.0921913261169291,0.2804580865096713,0.4836218673746996 -2016-10-07,0.17622628927693065,0.042522937812212294,0.2677792278193608,0.09151795533087594,0.29613309761016915,0.5009348684255963 -2016-10-15,0.17194096820846982,0.0435603027246404,0.27557789694192314,0.09047729199950669,0.31234809467462854,0.5190577215800363 -2016-10-23,0.1663121268846292,0.04445048677137976,0.2828608057295408,0.08903707685879486,0.3287989040030691,0.5377843488306101 -2016-10-31,0.15961957123583692,0.045184425653889984,0.2893495888697367,0.08728396865797006,0.34474931217605154,0.5564261853769862 -2016-11-08,0.15230186202681248,0.045762055520745984,0.2947762824287796,0.08535121144753283,0.3593455042847912,0.5741364329888559 -2016-11-16,0.14497482402673384,0.04619640034535784,0.2989965999103642,0.0834202950999115,0.37180245453446614,0.5900179604642742 -2016-11-24,0.13827053054468325,0.046505421554232304,0.30195160024514267,0.08167579670014263,0.38147329614152303,0.6032529410958035 -2016-12-02,0.1326429005817493,0.04671162523757717,0.3038063665105309,0.08024758547222223,0.3882625949186302,0.6134692090907994 -2016-12-10,0.12841998384695275,0.04684002356019983,0.30481811287059135,0.07922935661846232,0.39234754892518,0.6205405396309353 -2016-12-18,0.12559368665307422,0.04692096708034673,0.30535534407642106,0.07862327050883267,0.39439524350659666,0.6248574622884235 -2016-12-26,0.12401024228923921,0.046986336227380646,0.3057936186689042,0.0783920311272501,0.39521582652275056,0.6270024735271167 -2017-01-03,0.12331315512417564,0.04706724331308916,0.3063793536752066,0.07844389149617018,0.39558520932367564,0.6277180405169064 -2017-01-11,0.12308273451847876,0.04718263335688684,0.30718771189440813,0.0786619686216636,0.3960583741808504,0.6277150201859385 -2017-01-19,0.12289886479484101,0.04733131970423883,0.3080765041794923,0.0789154239361927,0.3968420933432869,0.6275700950774494 -2017-01-27,0.12242826997839217,0.04748199605428784,0.30879202487662505,0.07907327818566663,0.3979012782487554,0.6277073793591393 -2017-02-04,0.12148124816177552,0.04755684721901366,0.30893358710100494,0.07900656703850047,0.3988609509330849,0.628321942567538 -2017-02-12,0.11996764435386673,0.0474517112651253,0.30812404403176785,0.07859165322746839,0.3993100326287613,0.6295249536919673 -2017-02-20,0.11793798418219947,0.0470270690934767,0.30605066236785317,0.07771378716165964,0.3988354469455417,0.6313267032171151 -2017-02-28,0.11549676404063249,0.04615652203349541,0.3025371551368125,0.07628499321703687,0.39717464628150534,0.6337198103395375 -2017-03-08,0.11281461506670445,0.04474193863123022,0.2976221542696675,0.07426081803006766,0.39431323133665425,0.6366825768047623 -2017-03-16,0.11007675500533592,0.042746910340314144,0.2914733418454173,0.07165697044138954,0.3903599893897095,0.6401299870738106 -2017-03-24,0.10748176250031516,0.040232686572447135,0.28443538626375503,0.06858179544309453,0.3855843508854543,0.64387296671069 -2017-04-01,0.1052238176580875,0.03733097954735127,0.27687988533988994,0.06520937455437095,0.3802165783996265,0.647584872306862 -2017-04-09,0.10348733427052105,0.03427338346002464,0.26919313350729707,0.06180671692873349,0.37439532519332963,0.6507292811440112 -2017-04-17,0.10243010928967988,0.031310874486717974,0.2616991823458007,0.05865723220982574,0.3681467801423657,0.6526842233446282 -2017-04-25,0.10216697986660327,0.028709158424805783,0.2546169770627527,0.056055206981936175,0.36132815596967194,0.6527235124505237 -2017-05-03,0.10276046041439311,0.02668074330808759,0.24809826863076054,0.05423744325318395,0.3537862327249869,0.6502230810344863 -2017-05-11,0.10419352931244594,0.025347366018594764,0.24220161346375257,0.05334344080194322,0.345380207117472,0.6447471560193706 -2017-05-19,0.10640298259468453,0.024748001152889068,0.2369704287842807,0.05342324968772539,0.3360940789788963,0.6361139500795719 -2017-05-27,0.10926300101547709,0.0248001279948309,0.23243855008284425,0.05439545434752498,0.32611576730059955,0.6245304386187589 -2017-06-04,0.11265441414573896,0.02537754203618641,0.2286582844964608,0.056130475285925595,0.3157484556289821,0.6104146857873088 -2017-06-12,0.11648003304458363,0.026302480798052738,0.22571350101323312,0.05844146464282989,0.30545028782672917,0.5944677678932532 -2017-06-20,0.12070173937185405,0.027419962688206824,0.22368764474342376,0.06116511079769345,0.2956475897489557,0.5774033264422902 -2017-06-28,0.12537521527779438,0.028621329855802687,0.2226576261639141,0.06418658855376597,0.28668728648584874,0.559907734760478 -2017-07-06,0.13060158042373554,0.0298490486256732,0.22266146960361177,0.06744347572812125,0.27877181740705803,0.5425172201031618 -2017-07-14,0.13654050270625256,0.031122005493374703,0.22367380203394865,0.07095140290482818,0.27187848277595605,0.5255245666206467 -2017-07-22,0.14331821398725325,0.03247959389873803,0.22561999488165474,0.07474149184502549,0.26587938668974065,0.5090732811219957 -2017-07-30,0.1509952097743944,0.0339820147408855,0.228354658113779,0.07885703070572471,0.26051108799720224,0.49310113933321814 -2017-08-07,0.15953929942927073,0.03566419184285347,0.23170942473214767,0.08330671745806364,0.25552519314838756,0.477502050594214 -2017-08-15,0.1687765301420638,0.037515971389211505,0.23548994811585416,0.08804030650737087,0.250720611232391,0.46213140033741285 -2017-08-23,0.17848519939639781,0.039499319203539815,0.23951794831068252,0.09297778847970586,0.24596328239926663,0.44687152581561723 diff --git a/inst/lucc_MT/patterns/cerrado_rupestre.csv b/inst/lucc_MT/patterns/cerrado_rupestre.csv deleted file mode 100644 index efa3ab2..0000000 --- a/inst/lucc_MT/patterns/cerrado_rupestre.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2004-08-28,0.16803610264828206,0.042324438815742285,0.21400416970292646,0.1079630571507581,0.17134922968041633,0.32929107511306355 -2004-09-05,0.17555257460987575,0.043134901994650886,0.21690089010194144,0.11129184040139427,0.1692318722550313,0.3197568310723091 -2004-09-13,0.1821166877829931,0.04401805318516535,0.21966562204018658,0.11425557584101566,0.16776558900492416,0.31184040677526254 -2004-09-21,0.18676553830513645,0.04501734203373825,0.22215240085409824,0.11646491356017848,0.16761938658062658,0.30721938926898246 -2004-09-29,0.18877366539194163,0.046018227571200054,0.22418998815809538,0.11752719666073491,0.16934999186230976,0.30736317928850776 -2004-10-07,0.1876267564317335,0.04683705323778552,0.22561878652264128,0.11709475164182057,0.17339522384775996,0.3134641065886451 -2004-10-15,0.18346206145533694,0.04719329985341725,0.22637095858521444,0.11503687011034555,0.17979516434092402,0.32570982319085173 -2004-10-23,0.17671732430366038,0.04681580401363805,0.22644983458753604,0.11135944400831957,0.18839556952655018,0.34375685814614776 -2004-10-31,0.1683083556457451,0.04561761207182022,0.22606235896795693,0.10640022683892134,0.19869881333940687,0.36620419742271054 -2004-11-08,0.15928997723414412,0.04363612933871917,0.22550918985872115,0.10064265786953652,0.21008617574047017,0.39124549652446927 -2004-11-16,0.1506846106792662,0.041228914318594975,0.22523906619094028,0.09477888259946485,0.22186187649615283,0.4166858463861882 -2004-11-24,0.14340928399795821,0.03886735583494597,0.22572256020787995,0.08953737847422169,0.2333436423211057,0.44033125951134644 -2004-12-02,0.137931660212459,0.03709726418541179,0.22731922906303112,0.08551695596211864,0.24402313161846342,0.4605512049940992 -2004-12-10,0.13447910772639562,0.03641829496495025,0.23027891803877082,0.08318793622803608,0.25349789801889305,0.47612038868502 -2004-12-18,0.1327935042301796,0.03700534491765605,0.23449402629306712,0.08260167174505864,0.2616139887935616,0.4869120257070938 -2004-12-26,0.13242948596935772,0.038846123691761915,0.23968096775366565,0.08360249999849194,0.2683246857584337,0.4933016555937485 -2005-01-03,0.132758424913687,0.0415135464163046,0.24524135129993133,0.08566187909536663,0.2737315093950955,0.49645158848885784 -2005-01-11,0.13315642036133218,0.044382477108008786,0.2504584064244347,0.08809594701100994,0.2779930640634103,0.49780289109359444 -2005-01-19,0.1331020284894823,0.046715190133039416,0.25458281113109466,0.09015673933518878,0.28128650266986754,0.4988516961959615 -2005-01-27,0.13224973300428097,0.04795002668619585,0.25704241644906617,0.09126474953699248,0.2837477404239786,0.5007140754346411 -2005-02-04,0.13052776777051908,0.0478359509868845,0.25756419390866586,0.09113028173434651,0.28544160422642273,0.5038681836241675 -2005-02-12,0.12802002972337195,0.04639947489382567,0.25612298643365594,0.08970885752487481,0.28638412308640593,0.5083200706008232 -2005-02-20,0.12502168320428325,0.04407221758594865,0.25305254760657725,0.08731189697920673,0.2865228030358324,0.5134042741399504 -2005-02-28,0.12186268385546725,0.04137115164265633,0.24878746764179266,0.08434517840317837,0.2857943634839996,0.5183229078402672 -2005-03-08,0.11889492134075186,0.03888467315011958,0.24386802354980908,0.08130980988756822,0.2841293954479244,0.5221734315317936 -2005-03-16,0.11639776632117553,0.03702617604088738,0.2387382089302926,0.07859814483389371,0.2814898407730972,0.5243012533692305 -2005-03-24,0.11452706689983405,0.03590815492157088,0.23366553854115313,0.07640708770334145,0.27789065217979486,0.5244782825944435 -2005-04-01,0.11334601310071829,0.03543410721791615,0.22877831732496554,0.07478151667768092,0.27339286031855076,0.5227609752216664 -2005-04-09,0.11278552783462997,0.03521456631908164,0.2240056557051444,0.0735491045071807,0.2681248384816574,0.519602074520454 -2005-04-17,0.11274024998849085,0.03486379428182098,0.2192642446838401,0.07252012197105982,0.2622437084099787,0.5154549006371255 -2005-04-25,0.11306718838371904,0.03405258575704504,0.21448431861901668,0.07151325009838276,0.25594276206270267,0.5107061503078357 -2005-05-03,0.1136550155011623,0.03264720557192684,0.20970377726002157,0.07046230566833264,0.24941121561190333,0.5055093922713665 -2005-05-11,0.11445104805042956,0.030826446944371637,0.20513529831526162,0.06948923401753113,0.24282532804747056,0.4996458757856064 -2005-05-19,0.11545681088992249,0.02890572914475424,0.20106036387275564,0.06880246765632864,0.23633114960690593,0.4927579195449367 -2005-05-27,0.11675339008460309,0.027371755445549466,0.19784324416213322,0.06872036329629296,0.23002679943116744,0.4843235972332026 -2005-06-04,0.11844470331047324,0.026618397018290274,0.19576834881433355,0.06950141740664792,0.22397467394477374,0.47397563308466556 -2005-06-12,0.12065876480157772,0.026856690100934767,0.19497351255861972,0.07128616785552942,0.21818603956611807,0.4616307731461985 -2005-06-20,0.1234964173303709,0.028088069238573772,0.19544103509494787,0.07407345967182051,0.21265016345062215,0.4474863912497309 -2005-06-28,0.12700933879301315,0.02999966364841647,0.19692146657664086,0.07764853457141635,0.20732843763805942,0.4321539805689269 -2005-07-06,0.13119733168096703,0.03220333627271604,0.1990994495332455,0.08173490369211492,0.20217892762962802,0.41632018256028924 -2005-07-14,0.13598270892992617,0.034240898137866994,0.20159054944900473,0.08599090408940188,0.19715876599533907,0.40072390072733055 -2005-07-22,0.14125472066471118,0.03579158561995119,0.20409180089688897,0.09015370882773219,0.19223036539526062,0.3858795967359355 -2005-07-30,0.1468611769500946,0.036772999150998,0.20644625752473272,0.09410151238695116,0.18736492902458235,0.3719223931040041 -2005-08-07,0.15266119177467535,0.037257014670789394,0.20860661842183886,0.097817179552961,0.1825367148648554,0.35874343298078726 diff --git a/inst/lucc_MT/patterns/ciliary_forest.csv b/inst/lucc_MT/patterns/ciliary_forest.csv deleted file mode 100644 index 4d10167..0000000 --- a/inst/lucc_MT/patterns/ciliary_forest.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2015-09-14,0.1730252110514431,0.03774239172130268,0.2606665155515582,0.07739066617825621,0.3210740421633871,0.5425203930902326 -2015-09-22,0.17195406413794967,0.03915698311181836,0.2685809700726744,0.07835991000413925,0.3319952559856766,0.5514635696525839 -2015-09-30,0.17061103334965388,0.04062626096203241,0.27646159896724004,0.07930809455369134,0.3429948257480173,0.5605645556242084 -2015-10-08,0.16872655573042794,0.04219324981092722,0.28425087549527106,0.08020312732293186,0.35413130251594643,0.569975160888747 -2015-10-16,0.16612361568010214,0.04383024575501085,0.29179117364432733,0.08096790972174275,0.36534750209101785,0.5797717580489126 -2015-10-24,0.16269456514084377,0.045471258355967645,0.29888489557180153,0.08150704369818866,0.37652791975364025,0.5899822654560063 -2015-11-01,0.1585247367659139,0.04697205448509318,0.30529071744509606,0.08170388809250337,0.3874495379008145,0.6005154385310666 -2015-11-09,0.15377979098130123,0.04816740727439125,0.310776716721543,0.08144562749648815,0.39786563184088064,0.6112351570493695 -2015-11-17,0.14872605211490594,0.048906031746056636,0.3152185294872718,0.0806682852305077,0.4075726694870524,0.6219652518094847 -2015-11-25,0.143649609378167,0.049065449326820595,0.3185633836460529,0.07934510607606662,0.4164082306496505,0.6325198850024554 -2015-12-03,0.1387851740075578,0.04863890108674353,0.32094082768650267,0.0775545557348781,0.42434085550931216,0.6427245357141473 -2015-12-11,0.13432740576572194,0.047681208638630276,0.3225509506536031,0.0754251738299602,0.4313860229645438,0.6524005411047172 -2015-12-19,0.13035659650515002,0.04637069892340309,0.32366340847144387,0.07316996723410156,0.43759126875204446,0.6613328739170579 -2015-12-27,0.12690464953914293,0.04492736186168716,0.3245413302640499,0.07102436121998414,0.44298787583326904,0.6692825616169509 -2016-01-04,0.12394002804924042,0.04358886241166668,0.3253151939667332,0.06920679848181863,0.4474835743137828,0.6759510221670952 -2016-01-12,0.12141010130522241,0.0425493058679681,0.32596547839097667,0.06787286682710636,0.4508837915545919,0.6810472097653935 -2016-01-20,0.11925956764278334,0.04192103064296524,0.3262927034479784,0.06708217667496108,0.45288942698098605,0.684325455727692 -2016-01-28,0.11744240696256082,0.04172530418511746,0.3260317850806916,0.06680425373117749,0.45323550906620447,0.6856944431813315 -2016-02-05,0.11592794614989685,0.041846574399823526,0.32484045051299876,0.06687577838455772,0.4517326349399708,0.6853118429998343 -2016-02-13,0.1146788535141675,0.0421088867021722,0.322440049109112,0.06707990883294768,0.4483510702440134,0.6835443655999941 -2016-02-21,0.11364410790124618,0.04225985831347469,0.3186707766146458,0.06713410482655087,0.44331256090657295,0.6810616050755124 -2016-02-29,0.11275534496345459,0.04206901076082469,0.3135136282669886,0.06678518555715561,0.43698482034994013,0.6785947961045631 -2016-03-08,0.11191821089381046,0.04135541881305995,0.3071623118097339,0.06584014585093277,0.42993339483110415,0.676923633390467 -2016-03-16,0.11104439516684839,0.040024470268298404,0.29990474543590273,0.06420017046873548,0.4227037821052507,0.6766255813690896 -2016-03-24,0.11006048243560124,0.03811610343528612,0.2921484048628566,0.06191235176577515,0.41577383515455013,0.677933666520054 -2016-04-01,0.10893680791956595,0.03575071910665331,0.2842801171470629,0.05911306028071059,0.40944427607082956,0.6807433709045329 -2016-04-09,0.10771378691796242,0.03315894681952403,0.27663235102194944,0.05606162694286817,0.40373435479904834,0.6844502406546452 -2016-04-17,0.10648386015737006,0.030585435472649605,0.2694485072627337,0.05303822420113981,0.39850043655149336,0.6882587046730497 -2016-04-25,0.10541356877662242,0.028279196174783586,0.26283418337992914,0.050336516737043244,0.3933779392085337,0.691142998206772 -2016-05-03,0.1046859935928899,0.026435637903060732,0.25683476295878566,0.0481971212056622,0.3880074641690315,0.6921958335979996 -2016-05-11,0.10449883608315141,0.025160697159723404,0.25141950855721407,0.046770514715414074,0.38207608243901003,0.6907457982722294 -2016-05-19,0.10501860998651152,0.02449152426127434,0.24656919204646674,0.04613202218006321,0.3754208169310879,0.6864304141687178 -2016-05-27,0.1063576406665929,0.024366659833595566,0.24229619257028323,0.046247894391370906,0.3681184803135134,0.6793581776961487 -2016-06-04,0.10857750283002593,0.024695019585322976,0.23865348042863443,0.047045400939115,0.3603769521321101,0.6698617219856993 -2016-06-12,0.11166487949071749,0.025353177797291702,0.23575523607784998,0.04840570636685616,0.3525796381773046,0.6585643810745311 -2016-06-20,0.11557127482753647,0.026236595636069775,0.23371969715485602,0.05022095869692364,0.3450997692101112,0.6460757495122167 -2016-06-28,0.12020245616078379,0.0272772640650532,0.23266101028788994,0.052409015760074056,0.3382638662914396,0.6329339706501879 -2016-07-06,0.12545438844432574,0.02843906311782664,0.2326450488246496,0.05491457106917877,0.3322824528405641,0.6195216304930667 -2016-07-14,0.1312143105650785,0.02973219392460348,0.23365837470507184,0.05772204344834635,0.3271824279847293,0.6059771270713047 -2016-07-22,0.13737302774741306,0.0311734373567488,0.23562676194999765,0.060819714861953394,0.322895369293491,0.5923396907564713 -2016-07-30,0.1438286649562781,0.03278254978899412,0.23838887174143547,0.0641954537572264,0.31922641044975747,0.5785261351647004 -2016-08-07,0.1504876887380759,0.0345593267847064,0.2417574506874674,0.06781758950859211,0.31597560641750433,0.5644796704934313 -2016-08-15,0.1572693521988186,0.03647033095967021,0.24551734337504683,0.07162255458586206,0.3129556681789383,0.5502067327957032 -2016-08-23,0.16411010963225026,0.03847177035892044,0.24948415532830887,0.07554214711388682,0.3100284143400543,0.5357560459584966 diff --git a/inst/lucc_MT/patterns/dunas.csv b/inst/lucc_MT/patterns/dunas.csv deleted file mode 100644 index 2f5a409..0000000 --- a/inst/lucc_MT/patterns/dunas.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2005-08-29,0.4629688219030038,0.2226423352501511,0.5079676914688456,0.42703830190039616,0.08432125917021596,0.08527651037921133 -2005-09-06,0.4738510788694299,0.23123868967610883,0.5172463967699534,0.43512246006719735,0.08599218503204956,0.08526462513362515 -2005-09-14,0.4845790275319588,0.23985506369975843,0.5264365969313871,0.4431433141310613,0.0876390608444868,0.08520390852463905 -2005-09-22,0.49495961365762264,0.24849871346653057,0.5354313653721181,0.4510211677393378,0.08923749343419643,0.08505084286450143 -2005-09-30,0.5046733842420393,0.2571347211762355,0.5440637252344954,0.45862278565475467,0.09076211222602179,0.08477970700719226 -2005-10-08,0.5133630218852937,0.2656932461050963,0.552147490266262,0.4657945235286248,0.09219062569314293,0.0843842265639477 -2005-10-16,0.5206347458768701,0.2740193414745895,0.5594750324610754,0.4723529244449296,0.09351137927908046,0.08390922701950934 -2005-10-24,0.5261024570433886,0.28191790106485726,0.5658346771226026,0.47810713993702414,0.09471330969795602,0.08341437527710664 -2005-11-01,0.5294305612943101,0.28912144270802587,0.5710163347500548,0.48286281907390544,0.09577981099911145,0.08297228104717042 -2005-11-09,0.5303321050907215,0.29534676702183843,0.5748190182983692,0.4864356352651386,0.09668539046493507,0.0826436623529683 -2005-11-17,0.5286446093481283,0.30030959582975103,0.5770705898967505,0.48867666639496565,0.09737931458839044,0.08244130529081324 -2005-11-25,0.5242835993627087,0.3037541801085056,0.5776265301788971,0.48946909090251756,0.09779984596183941,0.08234838610271517 -2005-12-03,0.5173339591455839,0.3055191378296938,0.5764129780559005,0.48877683029080843,0.0978668288494316,0.0822845148342235 -2005-12-11,0.5079656919798462,0.30551097367253105,0.5734036263572866,0.4866130542739263,0.09750076044373647,0.08214963486428943 -2005-12-19,0.49651023005381845,0.3037931508140577,0.5686790433805681,0.4830962983508774,0.0966371411367733,0.0818234922327119 -2005-12-27,0.4833615969543999,0.30050533719478417,0.5623734214514161,0.4783921960684581,0.09522842249383365,0.08119224711582473 -2006-01-04,0.46900721491911745,0.2959293405921672,0.5547237181628464,0.4727480793935433,0.09327324327537655,0.08017824481886257 -2006-01-12,0.45396353545050633,0.29042042666995116,0.5460258551664434,0.46644502113885977,0.0908119121677319,0.0787472125440149 -2006-01-20,0.43875634446199263,0.2843836195908201,0.5366191662443904,0.4597821485587568,0.08792313525433805,0.07690819616655922 -2006-01-28,0.4238686910598898,0.278186620776189,0.526830473775949,0.45302264257856956,0.08470706417247391,0.07470299571037328 -2006-02-05,0.4097303725433442,0.27214590739157146,0.5169668684886654,0.4463836604259594,0.08128646688247962,0.07220961147165854 -2006-02-13,0.39670693046367317,0.2664763344676194,0.5072845446824386,0.4400156155214512,0.07778949693060956,0.06952932033646296 -2006-02-21,0.38509004231028626,0.2612734077707015,0.49797950102453636,0.4339919495747513,0.07434882341875969,0.06678876070578674 -2006-03-01,0.3751105908929529,0.25653453288420824,0.4891985161393981,0.4283355561120068,0.071092025643499,0.06413106773472436 -2006-03-09,0.3669310275141391,0.2521494182382644,0.4810341816744554,0.4230153975939589,0.06813981506175125,0.06171682916695892 -2006-03-17,0.36064992788127853,0.2479792523876185,0.47355957190724846,0.41799743286495267,0.06559722035855835,0.05970301049692295 -2006-03-25,0.3562927903939437,0.2438615758473108,0.46682895022888626,0.4132493451313755,0.06355018694626388,0.05823943149090066 -2006-04-02,0.3538054043223342,0.23968955653370644,0.4609016710827867,0.40877501207270617,0.06205336237525494,0.05743163642116151 -2006-04-10,0.3530424366618055,0.23542675106119693,0.4558445474280578,0.40462207012595397,0.06112448947610209,0.05733084743489446 -2006-04-18,0.3537859985393515,0.23112568019114746,0.45172754726017456,0.4008745841450187,0.06074098563265792,0.057913105995296844 -2006-04-26,0.3557391798943752,0.2269399643220189,0.44862313146081767,0.39765626690961153,0.060834621423571614,0.05906803992825435 -2006-05-04,0.3585967694501978,0.22306778332993055,0.4465859237283171,0.3950912475581202,0.06131255403333023,0.060632234766964103 -2006-05-12,0.3620527394395767,0.21975082822186642,0.44564812606868426,0.39329975912827525,0.062056819327051684,0.06238698116120726 -2006-05-20,0.3658878065837976,0.2171898367576635,0.4458076672573679,0.3923613777094711,0.06296296726994331,0.06413469088779022 -2006-05-28,0.3699895157501927,0.21553070666793095,0.4470226570983512,0.3923061617461689,0.06394586075658526,0.06570977527673505 -2006-06-05,0.37437460306757664,0.21482362896931476,0.4492197009755341,0.3931096483323835,0.06496349686971328,0.06703150385676987 -2006-06-13,0.379206121009636,0.21500754638471709,0.4522907409964312,0.3946856765474993,0.06602440203176019,0.0681192436416738 -2006-06-21,0.3847065046600898,0.21594722949514578,0.4561134442697633,0.39691523502054876,0.06717022559242389,0.06906443512142761 -2006-06-29,0.39115440455932693,0.21742843671692896,0.46055135096548583,0.3996453838975763,0.06847747343581022,0.07003554310638914 -2006-07-07,0.3987488266761583,0.21923899647741846,0.4654718040123056,0.40272868504943227,0.07001188754768936,0.07119010354574357 -2006-07-15,0.4075851593863967,0.22117726179576847,0.4707481840148303,0.4060278084692984,0.07182236468736589,0.07266420413158735 -2006-07-23,0.4175999206227614,0.2231049364083778,0.47626891521028636,0.4094374478993898,0.07391337816224744,0.0745135369176905 -2006-07-31,0.42854765538332823,0.22495943614207597,0.4819404106563865,0.41289030710692554,0.07623741776691409,0.07669878284278438 -2006-08-08,0.44010978603312056,0.22673267823274734,0.48768957567940413,0.41634914175410187,0.07872058768260476,0.07912525296594616 diff --git a/inst/lucc_MT/patterns/fallow_cotton.csv b/inst/lucc_MT/patterns/fallow_cotton.csv deleted file mode 100644 index 661d7f5..0000000 --- a/inst/lucc_MT/patterns/fallow_cotton.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2014-08-29,0.32458615045236644,0.04483111222981985,0.27321659910476637,0.1551470341887938,0.1583254776412885,0.2623471682136184 -2014-09-06,0.3106759881639702,0.046990202845903134,0.2713149412945277,0.1503118441994078,0.16752266274276356,0.28003559786476884 -2014-09-14,0.2984662888558146,0.04884060255082182,0.2695560963000277,0.14582047958662542,0.1759782965852301,0.2961631108272869 -2014-09-22,0.28929272405563733,0.050122808829574345,0.26813651278010153,0.14193875055425664,0.1832361969534218,0.30957336331254304 -2014-09-30,0.2832820111326636,0.05074097420463183,0.26742795180309764,0.13867416926227608,0.18978214159669332,0.32044883235357435 -2014-10-08,0.2797454438835776,0.05075207353046599,0.26783912278400923,0.13589056259485677,0.19657371907072457,0.3297584146423624 -2014-10-16,0.2761765952149055,0.05059876599788874,0.26977698706618874,0.13316422253393573,0.20544902977884422,0.34010057132080396 -2014-10-24,0.26991685436813484,0.050855163103788625,0.2735421817097318,0.13013425393654007,0.21802171760053335,0.3539347367106866 -2014-11-01,0.2590792204424169,0.052256916335086355,0.2790876810381221,0.12684679619083758,0.23455343144456056,0.3721249903620214 -2014-11-09,0.24286089206123418,0.05547480831891833,0.2861564659893394,0.12365806804924583,0.25424802633486915,0.39410586219272603 -2014-11-17,0.22346814639352852,0.06085417793046075,0.2940462364350221,0.12165418649495799,0.2737917944464205,0.41574910357213035 -2014-11-25,0.20421615871820517,0.06846879734221127,0.30195913563771387,0.12203960437952124,0.28930414540247246,0.4320439365069533 -2014-12-03,0.18994984643706675,0.07772908196628821,0.30909529969244864,0.12589505075864232,0.2968088854940324,0.4375815720370841 -2014-12-11,0.18509484528864234,0.08775347644224649,0.31484617474534765,0.1338605135059722,0.29345378029753444,0.42838646847219564 -2014-12-19,0.19163759560721524,0.09716578022096907,0.31923144990544716,0.14521807726341895,0.28013565687365277,0.40545559741339143 -2014-12-27,0.20976754624159252,0.10452887098591955,0.3226229101850908,0.15855779382103066,0.25980607859124794,0.3725465202918141 -2015-01-04,0.235448050657969,0.10856551771205888,0.3261113362602031,0.1711508325413684,0.23954249216852946,0.33900240940170606 -2015-01-12,0.2624576952925817,0.10848077692581672,0.33104799362605997,0.1800058998110107,0.2276271003574432,0.3157599766336674 -2015-01-20,0.28312658619457065,0.1040313724303001,0.3388920123177188,0.18220126880273693,0.23260173782761523,0.3140616256419024 -2015-01-28,0.291681150592117,0.09558276070467521,0.350579252015439,0.1761260281195507,0.2595629541953228,0.3404521978421484 -2015-02-05,0.28481172037056257,0.08418624884804131,0.36639271906079574,0.1617613772945141,0.30937411754066024,0.3957181001150306 -2015-02-13,0.26268212573465305,0.07117444202274317,0.3858095103248504,0.1406011324149456,0.37848123141225976,0.4748951373809844 -2015-02-21,0.22941459388211002,0.05815094546501096,0.40739169336650544,0.11580926429793209,0.4584051185760658,0.5666446773601178 -2015-03-01,0.19062347885556416,0.04645544379200177,0.42928925257776934,0.09076982893556246,0.5395520582457403,0.6586910919111917 -2015-03-09,0.15332363621681655,0.03707753932514057,0.44924462537986304,0.06895856750780459,0.6114877414026436,0.7382900591826032 -2015-03-17,0.12279527944405015,0.03046208937915167,0.46527604703817127,0.05263344941433339,0.6667502392711206,0.797207684905126 -2015-03-25,0.10207970479612498,0.026427035196916415,0.47577928347495785,0.04256692048779066,0.7015656550086297,0.8327176855811949 -2015-04-02,0.09157650712572407,0.02446214233897985,0.47979351107615653,0.03825073730019449,0.7157547599197336,0.846728435347984 -2015-04-10,0.08872305821227383,0.023722041345228408,0.4771047932578324,0.0377790037986532,0.713116430122617,0.8461847132440975 -2015-04-18,0.09038403781808071,0.02349366543417893,0.46799529042921034,0.03922110043125131,0.6980659533580041,0.8378081089133941 -2015-04-26,0.09302954941309102,0.02325711540897861,0.4532821913091139,0.040757066013680406,0.6753767623532415,0.827544309153739 -2015-05-04,0.0947947694008833,0.022850146167779258,0.4339570116407338,0.041579475141877806,0.6477868337403597,0.8174309522908844 -2015-05-12,0.0958515355045177,0.0225221147568088,0.4111793713441183,0.04208693329589662,0.6155700989364364,0.8048965290385366 -2015-05-20,0.09774425185898966,0.022667343808525035,0.3861625876187858,0.04332746462616274,0.5777618197706642,0.7854150309486815 -2015-05-28,0.10345110783217765,0.02381202751703951,0.36017445948286764,0.04700328232848572,0.532181308308805,0.7526115914256568 -2015-06-05,0.11514261845310277,0.02625543278514033,0.33456814553084635,0.05425780245076179,0.47822811978285173,0.7030861941017337 -2015-06-13,0.13386400077489072,0.030008677342680718,0.3107914273671095,0.06549070340270738,0.41733749789828056,0.6372351224742523 -2015-06-21,0.15875937916761151,0.03476609196888986,0.2902342536171245,0.0800542327187306,0.35332040650982377,0.5598482876091457 -2015-06-29,0.18680396121032317,0.03987005814379809,0.2741952121194609,0.09611715516772,0.29261996488353736,0.48060601563554783 -2015-07-07,0.21452486501556195,0.04464578173359458,0.2635260576547996,0.1116685553151274,0.24136120601714087,0.40927119124934724 -2015-07-15,0.23814312270811827,0.048435950510942724,0.2585344236590764,0.12460871358894038,0.20497843993156117,0.3551639425917047 -2015-07-23,0.2555357721660967,0.050904276862920025,0.2588461351991495,0.13379822752396764,0.18557770821377076,0.3227179979005068 -2015-07-31,0.26660629321012436,0.05210058924219635,0.26331386210932706,0.13926498792422018,0.18131090548346243,0.31050943815502524 -2015-08-08,0.272788607452014,0.05234211762941006,0.2704966243353613,0.1418885196714222,0.18788818094858287,0.3132306812349024 diff --git a/inst/lucc_MT/patterns/millet_cotton.csv b/inst/lucc_MT/patterns/millet_cotton.csv deleted file mode 100644 index 480f407..0000000 --- a/inst/lucc_MT/patterns/millet_cotton.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2015-08-29,0.2604275634122556,0.050060081608309395,0.28126382183351173,0.13182369745021683,0.22899003728988915,0.3680519927867411 -2015-09-06,0.2810434198541945,0.05012558656890456,0.26780769902782814,0.13589170041093826,0.19447975104770093,0.3286470813471555 -2015-09-14,0.29862543341631603,0.05001973888196762,0.2569463537055084,0.13913760510954246,0.16682914875164329,0.29674623256519544 -2015-09-22,0.31028581335018574,0.049587114484218264,0.25103917943065246,0.14075632542738648,0.15235896767311669,0.2792826594649984 -2015-09-30,0.3137209774328378,0.04873600364938768,0.25150403409039584,0.14001082381971253,0.15523415838478394,0.2809061765369023 -2015-10-08,0.30730278360237423,0.04744895968319134,0.25904077392111724,0.13634577044762414,0.17778053734941823,0.3042653181980064 -2015-10-16,0.29152064191566696,0.045941873682259736,0.27241878423318455,0.12986461772175784,0.21712317128910774,0.3462868971063959 -2015-10-24,0.2679090022955916,0.04455994144444047,0.28966030423776107,0.12106995174071336,0.26815614989038905,0.4013943008223243 -2015-11-01,0.2400649220103271,0.043932239781534585,0.30772945967510973,0.11140211110542081,0.3220504044381008,0.45971893827357674 -2015-11-09,0.21219139382974095,0.04475424688932005,0.3235117992013837,0.10262359147716503,0.36931475650200496,0.5105798647835645 -2015-11-17,0.18885254518515804,0.047703165485517576,0.33463667660719865,0.09684776425630623,0.401532747193266,0.5443386998906731 -2015-11-25,0.1742149703238965,0.0533077337302309,0.33941166240119464,0.09603454007021989,0.41197893352876636,0.5532392464028868 -2015-12-03,0.17049399609572613,0.06152110703114184,0.3381133888297191,0.10117911071524124,0.3996172166070366,0.5360114759917782 -2015-12-11,0.1786525690642799,0.07194507933089608,0.33192133733982804,0.11257438316445739,0.3663980746587201,0.4948355627219174 -2015-12-19,0.1965953820461945,0.08335158146442176,0.32365732444370876,0.12866884392901293,0.32052895036126905,0.439206321366278 -2015-12-27,0.22101008291702215,0.09422742486629447,0.3165681063204915,0.1471627579826973,0.27222432592597473,0.38094389142804796 -2016-01-04,0.24678042592467944,0.10275562740222922,0.31394115826204483,0.1646137052642056,0.23342521135432942,0.33379896639812684 -2016-01-12,0.26861092872353887,0.10735051594687034,0.31832165993646,0.17755399143339504,0.21449221530361826,0.3094775907327749 -2016-01-20,0.28162281377345244,0.10687162096860113,0.33117992194601603,0.18290812813385965,0.22288464714884454,0.31604607345574653 -2016-01-28,0.28266646694939823,0.10095725158917122,0.35277828112454934,0.17893560415486406,0.2615987697456122,0.35617071234575043 -2016-02-05,0.27099079467258974,0.09029473406843232,0.3817950761134452,0.16573857389257288,0.3276577633654201,0.4253118642600916 -2016-02-13,0.24773686681815169,0.07620143297535682,0.4161065367567193,0.14479510170592877,0.41460508727315737,0.5150775307583125 -2016-02-21,0.21643846540402617,0.06075987982008719,0.45263875872482884,0.11930895013355111,0.5117686087816707,0.6124860787593174 -2016-02-29,0.18118865532387465,0.04596932095476696,0.4883609769780947,0.09270398792771281,0.6085792421626195,0.705429251466004 -2016-03-08,0.1465905746180397,0.03362009890270033,0.5204130443998819,0.06851098522280798,0.6951068045473999,0.7835732996492498 -2016-03-16,0.11627719229543335,0.024817052413168627,0.5465530202914952,0.04925820912460085,0.7643677247095945,0.8408214693247961 -2016-03-24,0.09245343954680027,0.01972185368830144,0.5653924900973403,0.03604687599723355,0.8134611461517876,0.876844927859886 -2016-04-01,0.07597779341733699,0.017892257846798692,0.5761390919334778,0.028842174625304695,0.8421759536436555,0.8944926325590805 -2016-04-09,0.06597574356231503,0.018144830935381637,0.578745997939311,0.026174225928916175,0.85362131739118,0.9004174441415054 -2016-04-17,0.06109071963117799,0.01931381445101922,0.573424805981414,0.026364612365294235,0.8511492779511751,0.9005483167084483 -2016-04-25,0.05948567157594012,0.02034448279474634,0.5606717680242436,0.027607112807597284,0.8380784027673932,0.8993537442743169 -2016-05-03,0.05983487291116413,0.020669124019616188,0.541085923277366,0.02876967298530645,0.8163178917614603,0.8982583234848963 -2016-05-11,0.061584744020291016,0.020391832376260837,0.5153857337745991,0.02968203250828371,0.7858120418281354,0.8945815905473193 -2016-05-19,0.06479502279895843,0.0199513660106216,0.48455654752473515,0.03077531293160924,0.7459287107514592,0.8842522448977497 -2016-05-27,0.0702962509026121,0.02017810170539111,0.449950472620899,0.03321298925915213,0.6954664956498116,0.8617336048777176 -2016-06-04,0.07886013808467196,0.02171725835799939,0.41331109629061746,0.03798621582322532,0.6345720861762769,0.8235244098209198 -2016-06-12,0.09108567077778273,0.024907671173560696,0.37688927698110564,0.04575815236953901,0.5653575420644534,0.7691652443744512 -2016-06-20,0.10701433861293741,0.029660694067143847,0.34300615540562357,0.05655621720999828,0.49159333023306595,0.7012668694039189 -2016-06-28,0.1259252794263854,0.035323167371147976,0.3140152370485871,0.06955349392278272,0.41914170890423036,0.6263982648846804 -2016-07-06,0.14680597965982006,0.04112110460545988,0.2916671389548497,0.08366389980838322,0.35348030781792394,0.5515374327402671 -2016-07-14,0.16832440153926548,0.04617783679281651,0.27685053720255776,0.09753828780528895,0.29928342513606804,0.48376255190647255 -2016-07-22,0.18947666752933917,0.0499462582514787,0.2694538662435771,0.11024549665892805,0.2586962452146369,0.4273067505324507 -2016-07-30,0.20975709765003048,0.05235432763622559,0.2680546498752702,0.12147198894272819,0.23039477099832806,0.38238810521524813 -2016-08-07,0.2291217585261297,0.05365961319110376,0.2707275855450678,0.13139122656377292,0.21115429489743398,0.3466382208257265 diff --git a/inst/lucc_MT/patterns/pasture.csv b/inst/lucc_MT/patterns/pasture.csv deleted file mode 100644 index 81db957..0000000 --- a/inst/lucc_MT/patterns/pasture.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2016-09-13,0.23920966863652407,0.05324099157934362,0.278130092609651,0.139916573962837,0.20285798112947495,0.3384492561508481 -2016-09-21,0.24136294926734936,0.05231182883929825,0.27974747193945737,0.13777853812831484,0.20930228838965234,0.3479968763672378 -2016-09-29,0.24220828329379546,0.05123655069631376,0.2815233187223889,0.1350779187378099,0.21725988746251143,0.3597817755394912 -2016-10-07,0.24047562059970928,0.04989640061261194,0.2836258452790298,0.13129415652440576,0.22817839703259934,0.37591705045767365 -2016-10-15,0.23546450933040886,0.048343881070807046,0.28622082681936245,0.12627134798483614,0.24274473072991956,0.39726535134113905 -2016-10-23,0.2268661857558188,0.046716932885645025,0.2894564263810738,0.12007133144052648,0.26115147259026317,0.4239047208205409 -2016-10-31,0.21539389469644535,0.04530823752237337,0.2934014869573877,0.11325461851591416,0.282362008115655,0.45402428493129243 -2016-11-08,0.20215681369138985,0.04443387143143766,0.29807984571718277,0.10654068748708932,0.3048941296889134,0.4851570662572863 -2016-11-16,0.18865767775529996,0.04429679944493434,0.30341727351998476,0.10067935212574486,0.3269368210452842,0.5144661534675093 -2016-11-24,0.17640512451525367,0.045002512919759784,0.30929366960662175,0.09632873746316788,0.3467510875493968,0.5393091059974018 -2016-12-02,0.16637038058663628,0.04636056417730973,0.3155135634822195,0.09367204790105688,0.3634086102844346,0.5584894421605889 -2016-12-10,0.1591773219514188,0.0480590615810741,0.32185084321116747,0.09264708810560122,0.37644339401607363,0.5716067415858185 -2016-12-18,0.15458840400014706,0.04964070387247814,0.32804318354519907,0.09272018306868839,0.386398274194833,0.5799261204336876 -2016-12-26,0.1520231552567197,0.05064283605121726,0.3338197494046638,0.09321341984741563,0.39416809138548686,0.5852768624879567 -2017-01-03,0.15052714555599478,0.05078949692508546,0.3389221552274749,0.09345091727261577,0.400849724040565,0.5897411937036444 -2017-01-11,0.14914728968402197,0.050026964576759904,0.3431254536827483,0.0929453727831405,0.40733017818939365,0.594980678082152 -2017-01-19,0.14712562865187537,0.048572907210943456,0.346254166577096,0.09152068413680625,0.4140451421218741,0.6018302718272516 -2017-01-27,0.1440336159250892,0.04673822276300134,0.34818582695175204,0.08922771302567682,0.421016582323149,0.6104098851840846 -2017-02-04,0.13995349766194276,0.04494082670755104,0.34887754548192906,0.08643872953985121,0.42766090903682835,0.6198005106126236 -2017-02-12,0.13519997168119918,0.04350174033999206,0.3483316077617575,0.08355989257182267,0.43326411853012103,0.6288751277850562 -2017-02-20,0.13038512326367674,0.04256355771119805,0.3466115574849073,0.08099860990127514,0.436999807172541,0.6363396224950634 -2017-02-28,0.1260833085517302,0.04209955155394222,0.34378871087292223,0.07900989589054438,0.43825422908260975,0.6412698037778047 -2017-03-08,0.12275672079678385,0.041834406742096805,0.3399369471957433,0.07759758951789368,0.4367829164494665,0.6433784245460198 -2017-03-16,0.12068517356649491,0.04145443978476237,0.33510967649041556,0.07663924323159235,0.4325831987929313,0.6427525085453476 -2017-03-24,0.11986694190674602,0.040625742832416004,0.32932657682681715,0.07584749032370458,0.4259882575299682,0.6399864355727191 -2017-04-01,0.12019764084786926,0.039157019420767825,0.3226097078638397,0.07498982532499404,0.4173453495195048,0.6355937113086831 -2017-04-09,0.1214511468113373,0.037096043229307145,0.31499174820157655,0.07395268837803472,0.406968893891707,0.6298681339548056 -2017-04-17,0.12345160908053812,0.034652417183704765,0.3065715886746145,0.07278596351346056,0.3950462084378469,0.622742749292191 -2017-04-25,0.1261261746446738,0.03226166173720521,0.2975565546408707,0.0717886884346916,0.38156263585967365,0.6135951174593004 -2017-05-03,0.1294919219413313,0.030355843113526196,0.28825472277107544,0.07133424753064999,0.3664959750837806,0.6016929516485673 -2017-05-11,0.13369606118850857,0.029322227011711277,0.2791222421857781,0.07187559004288807,0.3498477560925234,0.5862204937001947 -2017-05-19,0.13888179473579512,0.029383504015831247,0.2706572419451025,0.07375959569507021,0.3318047838331018,0.5667241169149208 -2017-05-27,0.1451630179493611,0.030498145875826808,0.26340109960571373,0.07714386188881729,0.31285584278542666,0.5433642875606252 -2017-06-04,0.15254995769892407,0.03247912252536433,0.2578070742653863,0.08200346743551917,0.2936623823974551,0.5168032011769416 -2017-06-12,0.1608864154903476,0.03494436640003517,0.254174125221138,0.08804398201484492,0.2751183445397487,0.48842726209557424 -2017-06-20,0.16991668601093768,0.03752867149309665,0.25261918303913544,0.09487911152408905,0.2580597716684922,0.459784182197853 -2017-06-28,0.17924905495563856,0.039927843917723235,0.2529837229012903,0.10201648917358751,0.243186242392901,0.43256415591308267 -2017-07-06,0.18848740534866124,0.041967410284604934,0.2549642781567556,0.10901758091892909,0.2309435786759792,0.4081636849137259 -2017-07-14,0.19725180921457278,0.0436762902313502,0.25806459781556984,0.1155507566494574,0.2213809238743919,0.3874686785387726 -2017-07-22,0.20525502971310947,0.04516244674320053,0.26179107548953995,0.12140394009728257,0.21432569427157389,0.3709068441244177 -2017-07-30,0.21236090106297717,0.046629523904123324,0.26568754328331695,0.12654475442533222,0.20930566492785635,0.35823358375402564 -2017-08-07,0.218556298527685,0.04823593920916203,0.2694281783659404,0.13102625688909492,0.20580553739336377,0.3489089890380152 -2017-08-15,0.2240062182896015,0.05004969678149681,0.272896145950505,0.13501276105229512,0.20328867360589803,0.3420173413660942 -2017-08-23,0.22893888464952597,0.05206403649641401,0.27609065767645774,0.13868372713747368,0.20130002644000372,0.3366292496869119 diff --git a/inst/lucc_MT/patterns/patt.rds b/inst/lucc_MT/patterns/patt.rds deleted file mode 100644 index 8ac5dd4..0000000 Binary files a/inst/lucc_MT/patterns/patt.rds and /dev/null differ diff --git a/inst/lucc_MT/patterns/soy_corn.csv b/inst/lucc_MT/patterns/soy_corn.csv deleted file mode 100644 index 78b4e22..0000000 --- a/inst/lucc_MT/patterns/soy_corn.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2015-08-29,0.2735236483084057,0.05639014052471263,0.29604342200768263,0.1646929078191135,0.18310394183559475,0.29310474229849093 -2015-09-06,0.28262967159658825,0.056608753539870334,0.28942838517603964,0.15957650440452187,0.1795791615009424,0.289918668874306 -2015-09-14,0.29094571644023703,0.05649555574004139,0.281846397247991,0.15484371386793155,0.17480143829133904,0.2866124552943123 -2015-09-22,0.2975205469574395,0.05573642103138438,0.2726817956272337,0.15077590823306314,0.16817492343812473,0.2836071755822302 -2015-09-30,0.30088990478237365,0.05419933429743339,0.26324460228696067,0.147063328410155,0.1625390063460346,0.2838801069775159 -2015-10-08,0.29947526773429356,0.051878672751977585,0.2557792479239479,0.14312764849672371,0.16230073352327654,0.2914601847330206 -2015-10-16,0.29190226674785824,0.04910162904888683,0.2541488135755081,0.13798126159487295,0.1742844627558268,0.31163231621920606 -2015-10-24,0.27714341102775797,0.04633301424365767,0.26242035279016207,0.13064672026394436,0.20535200428796785,0.34934726801892135 -2015-11-01,0.2555158690802326,0.04420663689868266,0.2832830964008524,0.12078631177116314,0.2590448687814893,0.40615889315326076 -2015-11-09,0.22810505184712734,0.04336557725277614,0.31816840028964744,0.108515965124543,0.33630637760629956,0.4812762546056717 -2015-11-17,0.1976808385173236,0.044252750193087456,0.3645630943313336,0.09516612013844777,0.4306183801566687,0.5677518717788571 -2015-11-25,0.16762382169172557,0.047145253548286636,0.41810146347683214,0.08251975797434233,0.5322150206512717,0.6561396196926902 -2015-12-03,0.14176475903301508,0.05182648667619364,0.4713717053979476,0.07273996824170342,0.6267810536028149,0.7340942919825874 -2015-12-11,0.12374662383909395,0.057837960104533494,0.516245187507782,0.06786443579623291,0.6993838415827822,0.7892704036486423 -2015-12-19,0.1156774422999326,0.0642920385371914,0.545604575692131,0.06879162950451645,0.7387386297128659,0.8135015887330183 -2015-12-27,0.11867426367118494,0.07017236296217254,0.5536253348301244,0.07571706402073346,0.7366786350810165,0.8016624625592101 -2016-01-04,0.1313796012085246,0.07447763142633639,0.5391594952709827,0.0871591566047357,0.6945442646037941,0.7569919902971272 -2016-01-12,0.1507783510469768,0.07643722932419093,0.5050386219248133,0.10068964068268144,0.6208981281615015,0.6886299090358838 -2016-01-20,0.17238656836495714,0.07565704724545788,0.458278345926067,0.1131812845017596,0.531362364343309,0.6111729586061826 -2016-01-28,0.19185130287867988,0.07210894339507738,0.4069328271239512,0.1219607397141028,0.44223317574854354,0.5391165308356789 -2016-02-05,0.20528969104440345,0.06630641632845394,0.36005478160679155,0.12522637378169946,0.36970627345748597,0.48584791922345594 -2016-02-13,0.21002227042129257,0.058976155972095684,0.3250358226180195,0.12220655705852407,0.32584751086843977,0.46079596122041955 -2016-02-21,0.205283119070888,0.051104803243390375,0.30634216880993437,0.113683888285296,0.315897524282285,0.46702168361632085 -2016-02-29,0.1915741923314455,0.043577542147811316,0.3055315911296879,0.10106876484557346,0.3398043127106176,0.5032053099825267 -2016-03-08,0.1711721069334221,0.0370672860539242,0.3198439285306078,0.0865452012917523,0.3899141731209043,0.5618551792524062 -2016-03-16,0.14683295110791714,0.03196438935663536,0.3446584429791989,0.07209713632468864,0.45607872499850455,0.6337366171405298 -2016-03-24,0.12180745367749851,0.028228562418791323,0.3730100748249182,0.05924894431159485,0.5253083888691018,0.7077703234599168 -2016-04-01,0.09899451797148237,0.025645423908416537,0.39830276916980745,0.0489711974144854,0.5859089791893901,0.7740312673514503 -2016-04-09,0.08062886983035777,0.023786522218997942,0.4149509958234092,0.04140027715999729,0.6288088035055733,0.8248111964845751 -2016-04-17,0.06830016410842245,0.02230256302565046,0.4192945922701785,0.036449682115088644,0.648069655811909,0.8546595428496521 -2016-04-25,0.06264187594532661,0.021011084808404974,0.41068967557661507,0.033805685741531956,0.6424979098671011,0.8614919771266916 -2016-05-03,0.06380478085618668,0.019921413710707866,0.39032027798741786,0.03334376441135852,0.6134279121902457,0.8450460042908713 -2016-05-11,0.07129816830265479,0.019330132254782986,0.3618853620066119,0.03528580603771342,0.5655855742780191,0.8075217810734723 -2016-05-19,0.08437201395438806,0.019585167061066756,0.32962358226487,0.03996528140404197,0.5046126543887965,0.7522498077966961 -2016-05-27,0.10195131458864273,0.021069241628143415,0.29824188788613154,0.047836269921278166,0.4371033377806099,0.6840296173457919 -2016-06-04,0.12282622884944508,0.02400178447299361,0.2717064687166513,0.059067779239664026,0.36928015245055107,0.6082562562708915 -2016-06-12,0.145618340592824,0.028326465392972316,0.2526788704117258,0.07333673370550782,0.30658922347896844,0.5310599062276113 -2016-06-20,0.16897669215058547,0.03379391329515136,0.24265687129253527,0.08994025724917933,0.2533919521172921,0.4582752772746646 -2016-06-28,0.19162755574463788,0.039883711260810056,0.24143226964869244,0.10763498307381841,0.21237548743690846,0.3951078067891884 -2016-07-06,0.21258400102614347,0.04607621064021891,0.24806479866397008,0.12521483608783068,0.1848352722399511,0.3452339843908175 -2016-07-14,0.2313233415983227,0.051899862074783996,0.26068933609491773,0.14162985627144736,0.17003822380179084,0.309895392350061 -2016-07-22,0.24773533118918856,0.057071585617900725,0.2774207833368457,0.15628784674234686,0.16615679976910075,0.28839893672742695 -2016-07-30,0.2623119369233548,0.06161059353267953,0.2964956739409139,0.1693241153382048,0.16990012924497339,0.2772722792361487 -2016-08-07,0.2756769673912236,0.06565975500383818,0.31658122379809783,0.181146543345654,0.17809051954117497,0.27266755186816016 diff --git a/inst/lucc_MT/patterns/soy_cotton.csv b/inst/lucc_MT/patterns/soy_cotton.csv deleted file mode 100644 index ef26a5e..0000000 --- a/inst/lucc_MT/patterns/soy_cotton.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2015-08-29,0.2772897583831439,0.05454503338749375,0.3052856008661409,0.14929503012122824,0.22350534437878933,0.33526267292731177 -2015-09-06,0.28786169883510315,0.05423765733740149,0.28457253823735645,0.14974872377309117,0.190069939435147,0.3080072820489781 -2015-09-14,0.2970593908543575,0.053866283236323045,0.2663139504228946,0.14992824312834016,0.16141389080247875,0.2848620432834411 -2015-09-22,0.3032836681139122,0.053343583196842403,0.2532051131318433,0.1494248188184882,0.1428988705984537,0.2705663215089862 -2015-09-30,0.3043033541548691,0.05249307180001413,0.2483243261503678,0.14728661659774303,0.14115845738193633,0.27154937656741085 -2015-10-08,0.2978519335917524,0.051129857591360885,0.25455990640136367,0.14242153194517215,0.16260235930327782,0.2942305525231539 -2015-10-16,0.2824137616984803,0.049159489618560544,0.2733471067191262,0.13386462616353373,0.2107349383419308,0.34252602811360977 -2015-10-24,0.2572319314609842,0.046592914088215254,0.3050143562237466,0.12098389755204508,0.28660934842219965,0.41796909802011206 -2015-11-01,0.22404931739615544,0.04380811884976997,0.3467600276928535,0.1045364285718336,0.3838958787693493,0.5144525523852844 -2015-11-09,0.18587404067189142,0.041353540400044304,0.39429317599715286,0.08602201282463426,0.49270067344362206,0.6221172985854935 -2015-11-17,0.14807412718599566,0.04005073795085719,0.4408136962019949,0.06843668944345059,0.5969740576677565,0.7244383237051774 -2015-11-25,0.11662672514666408,0.040737358520585076,0.4789607945056968,0.0551625299582584,0.679329018780638,0.8034322861936737 -2015-12-03,0.09697432956905477,0.04393148425794854,0.5023514756696208,0.04924038297699146,0.7249608617395765,0.8438662623840822 -2015-12-11,0.09358303444689842,0.049875339906399814,0.5058728075172608,0.05307553572135785,0.7223165662041107,0.8340147573152474 -2015-12-19,0.10708470807769097,0.057970997007026875,0.4889925027129889,0.06656804256465197,0.6714558354980691,0.7747277335706186 -2015-12-27,0.136082901285209,0.06723901437069033,0.45355030594606444,0.08832160578155988,0.5783494199027392,0.673198669926701 -2016-01-04,0.17513540144054454,0.07611478858321748,0.40597896004706485,0.11443742295218404,0.46011625251009186,0.5484126480625148 -2016-01-12,0.21689709103565466,0.0829541143917919,0.3548498828746523,0.14000097029350558,0.3385934056405868,0.42400618515731464 -2016-01-20,0.25296874955950965,0.0862871316345567,0.30989719659977766,0.15971603453151403,0.23764858289381763,0.32519578901603835 -2016-01-28,0.27659290087896415,0.08523718215932985,0.2790363201752771,0.16962480654647158,0.1758227361263221,0.27089676851971023 -2016-02-05,0.28391330833071327,0.07988021266327597,0.26696985284246544,0.1680520549440671,0.16248476162831335,0.26931714571914034 -2016-02-13,0.2738893233551201,0.07091564707850012,0.27528556857001374,0.15526036231741933,0.19892068827272102,0.31974781280404774 -2016-02-21,0.24964238484048615,0.059910383469762414,0.3009652933542765,0.13432598103271662,0.27461387219061173,0.4085651155774305 -2016-02-29,0.21576781261437458,0.04850557810853929,0.33925599801501993,0.10906275799634052,0.3754085564157979,0.5187795178314346 -2016-03-08,0.17858913163249787,0.03831626901009452,0.38331664280765576,0.08402301495067359,0.48313996953484567,0.630040169912891 -2016-03-16,0.14360621381381203,0.030470640580542607,0.4269198355784993,0.06282516863628114,0.5824892387157816,0.7261085856840055 -2016-03-24,0.11466711660193757,0.025297146762229637,0.4652003440635295,0.047420411672982915,0.6635081559696268,0.7981017911088364 -2016-04-01,0.0938416225590017,0.022628176495718296,0.49489391600417815,0.03838495262029371,0.7210951655157681,0.8430565931646528 -2016-04-09,0.08054578907012808,0.02161365654335625,0.5151460449593478,0.0343098842579312,0.75723042872291,0.8665237839676352 -2016-04-17,0.07337859479113706,0.02138710677772709,0.5259127398575216,0.03341106765856336,0.7754948985691275,0.8754791677296643 -2016-04-25,0.07003861959065581,0.021170808932013084,0.5280862044902822,0.033621507830038755,0.7809198637607058,0.877809872941293 -2016-05-03,0.06879335040311742,0.020553547373542258,0.5223856635222693,0.03359070441430608,0.7766769387111078,0.8782814554414025 -2016-05-11,0.06899540964817501,0.019712339139089898,0.5091016397284822,0.033194439446647725,0.7626525297416482,0.8763779888736223 -2016-05-19,0.07077642987150679,0.01906493917754294,0.48853235627684893,0.03298801275764185,0.7374378556611965,0.8691135689948968 -2016-05-27,0.07536299857227396,0.01932465322266882,0.4610751859664993,0.03443532695664912,0.6980734416681912,0.8503493115742164 -2016-06-04,0.08387426330951252,0.021030019367795624,0.4279364531621779,0.038785584461259945,0.6433813882512582,0.8156061876962883 -2016-06-12,0.09706537571695457,0.0243808144749889,0.3916412589668125,0.046779047500914155,0.5753897467354421,0.7636825440797135 -2016-06-20,0.11495877301962823,0.029237830221416804,0.3554058728938727,0.05843404267704357,0.4987279106604723,0.6968015439464128 -2016-06-28,0.1364340593040581,0.034961678269016805,0.3233507971211298,0.07266783742993721,0.421749198279428,0.6222597021559757 -2016-07-06,0.1600631512076564,0.040850503336804664,0.29893304563862677,0.08815335867122595,0.3523849458113096,0.5478414516832701 -2016-07-14,0.18410779405492433,0.04618567074530517,0.2843485605211492,0.10335013516437597,0.29721452794155134,0.4812561287959476 -2016-07-22,0.20730953547056857,0.0505242102602116,0.2801034346724192,0.11720915584890468,0.25934353428554524,0.4270730629299113 -2016-07-30,0.22925137486163066,0.053879462839976244,0.28413294286636315,0.1295254676894861,0.2363863858307801,0.38470756490150926 -2016-08-07,0.24999702584650102,0.05646438227358441,0.2937100496578344,0.14052409562056478,0.2240682096905262,0.351399756130339 diff --git a/inst/lucc_MT/patterns/soy_fallow.csv b/inst/lucc_MT/patterns/soy_fallow.csv deleted file mode 100644 index 7cb37b7..0000000 --- a/inst/lucc_MT/patterns/soy_fallow.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2007-08-29,0.31271536153606305,0.07991050188865595,0.3141044128659724,0.1996358149628646,0.1493434750927691,0.21782121255651327 -2007-09-06,0.3175608291777792,0.07967701777539567,0.30970236034202203,0.1987927200529474,0.1458651348252059,0.21902936900796455 -2007-09-14,0.32145911989488973,0.07951843795104174,0.30609160951868575,0.1976191937930884,0.1440148107378278,0.22090086918949986 -2007-09-22,0.32351450083860384,0.07950014004332509,0.3039263762924854,0.1957964868887932,0.1452196209682896,0.22403580076634003 -2007-09-30,0.32335400724412555,0.07962031863167618,0.3029878934989139,0.1931587988616228,0.1494880623515311,0.22854644568346402 -2007-10-08,0.3208965956464389,0.07983234896680783,0.3026787646943981,0.18962333226006925,0.15618878496432378,0.23438027189295033 -2007-10-16,0.3166791960897768,0.07997524498399214,0.3018832806605781,0.1852764115694328,0.163732869669302,0.24141787627483718 -2007-10-24,0.31138340031010336,0.0798407874009964,0.2995611551645492,0.1802355925521817,0.17058847895630067,0.24974353461926127 -2007-11-01,0.30548114024706524,0.0791629226465313,0.29577725253850495,0.17452164965814865,0.17687682252675654,0.26049989279277347 -2007-11-09,0.2992043517016313,0.07767207463783139,0.29130570201847295,0.16807459783532414,0.18380589065583322,0.27539743991679183 -2007-11-17,0.2919584098930772,0.07516552155735114,0.28866667007849917,0.1605944343103609,0.1953161145116755,0.29729354517818335 -2007-11-25,0.28282500243923747,0.07149819191269852,0.29092242704845905,0.15171108691975593,0.21618303426495852,0.3292534411622333 -2007-12-03,0.27055468762160717,0.06670076174598738,0.301167639556833,0.1410717292071978,0.2511690281024328,0.3736615465089182 -2007-12-11,0.25398208922078047,0.060898249360280834,0.32196162410393075,0.12842668969055804,0.3041239322570127,0.4320478620283726 -2007-12-19,0.23281680112249517,0.0544076216557853,0.3532909285853052,0.11398858843675551,0.37460960807555826,0.5028899338621586 -2007-12-27,0.2073427273451832,0.04762423883118836,0.3936967252303257,0.09823313001863014,0.4597656522271798,0.5830191109158057 -2008-01-04,0.1792961493691858,0.04103646628121533,0.4386072048308252,0.08221163775221059,0.5514788181605805,0.6658939298090727 -2008-01-12,0.1512907139880771,0.035138596194084534,0.4819185942033948,0.0672662660137844,0.6390016537005668,0.7433139553702076 -2008-01-20,0.12663938731970079,0.03038600689389058,0.5166076158564208,0.054921259524816066,0.7099542975022277,0.8060852919270844 -2008-01-28,0.10826212566931426,0.027122901838583086,0.5370030764405032,0.046436102320865806,0.7542186876447473,0.8465847831710493 -2008-02-05,0.09830415376411755,0.02551516160853584,0.5398364837137136,0.042598378432701436,0.7657154028709476,0.8599751363418174 -2008-02-13,0.09777784398293388,0.025572978022000297,0.5243945052145247,0.043643415660724535,0.7428540442113705,0.8447018508760472 -2008-02-21,0.1059509245153828,0.02707848702658494,0.49378301101249955,0.048971897475594085,0.6907883625747196,0.8041774221630501 -2008-02-29,0.1211243565287065,0.029701334730431406,0.4526991813502885,0.05756316823928611,0.6177232139488585,0.7444333361177938 -2008-03-08,0.1402971641783749,0.03296367466516881,0.4078197827959242,0.06784190445888429,0.5357587792003454,0.6749871923631731 -2008-03-16,0.16040921887760196,0.036386776507991346,0.36526076559464515,0.07825924490930908,0.45630908602945675,0.6053470687472063 -2008-03-24,0.17852540772455197,0.03951597535149236,0.32987252308301035,0.0874055479352703,0.38898576817842057,0.5444039509873903 -2008-04-01,0.1924894735434821,0.04200691975898444,0.3046132149447342,0.09428797555531741,0.3400385290842428,0.4985412758184928 -2008-04-09,0.20141331312131283,0.04369429748704873,0.28952753572389817,0.09858047703745663,0.31046106463237505,0.47009183715860253 -2008-04-17,0.2053017308813999,0.04455807675562841,0.28324549233684815,0.10041021703174557,0.2983010933625765,0.4584654282804456 -2008-04-25,0.20542180564754614,0.044788483811698365,0.2825522542276557,0.10054130090587551,0.2976256637003866,0.45892315897810027 -2008-05-03,0.20337763663890274,0.04465525826453531,0.28428350813217995,0.09990050883900244,0.30210029255334153,0.46562332075313556 -2008-05-11,0.20111732416670972,0.04452389729272346,0.28566849896177476,0.09957575694659952,0.305398972661703,0.4716072734662273 -2008-05-19,0.20028188516279238,0.044725133390019184,0.2849250064370655,0.10047912434227171,0.3029119029885051,0.47119971077450856 -2008-05-27,0.20194676679244322,0.04551572706719878,0.281835039607101,0.10320449152456636,0.29286237460175596,0.4610399847432806 -2008-06-04,0.20665255994228743,0.04704548702083617,0.27694598210007787,0.10804148849041866,0.27532055764974955,0.44003363586451627 -2008-06-12,0.21415117232302894,0.04929981482781337,0.2718049376840811,0.11483029764516844,0.2529295513054011,0.4104401927156921 -2008-06-20,0.22388954882104634,0.05217965506780273,0.2679334306386278,0.12322785874200642,0.22882653330124358,0.37574720559530667 -2008-06-28,0.23494335190273377,0.055474250795226984,0.2666434058526208,0.13266395405254683,0.20653142801232355,0.34090014664468066 -2008-07-06,0.24645700888910935,0.05897382986115893,0.2687436642488757,0.14260004012919192,0.18874280261863527,0.3101430597305859 -2008-07-14,0.2577508622239063,0.062489166724132864,0.27424824537074066,0.15258813200122373,0.17673810107315255,0.2863114066638118 -2008-07-22,0.2684074424212344,0.06589448651575996,0.282773248510148,0.16233361907747368,0.17070276304495616,0.27054856645163305 -2008-07-30,0.2784119322567452,0.06916506739297254,0.2934119699170484,0.17177955458002486,0.16926090580975797,0.2614474183910035 -2008-08-07,0.2878983445953214,0.07231842188787968,0.3052806633630193,0.18096447755815045,0.17077480914234944,0.2569192160984753 diff --git a/inst/lucc_MT/patterns/soy_millet.csv b/inst/lucc_MT/patterns/soy_millet.csv deleted file mode 100644 index 3a8ef6b..0000000 --- a/inst/lucc_MT/patterns/soy_millet.csv +++ /dev/null @@ -1,45 +0,0 @@ -date,mir,blue,nir,red,evi,ndvi -2014-08-29,0.25801036009202993,0.04654557089021999,0.2711427648877147,0.14478424194996978,0.1841063410337807,0.30762946525541257 -2014-09-06,0.2591508622908625,0.04686739580622259,0.2721686933374464,0.13790192849638838,0.19678260054906715,0.33201488241527816 -2014-09-14,0.2598540555466566,0.04724641244252613,0.27358005522576934,0.13162550783114263,0.20959210296062192,0.3553108568876747 -2014-09-22,0.259711168594624,0.0477351227849719,0.27568104913051783,0.12652486247006575,0.22252377040419433,0.3763799751318044 -2014-09-30,0.2584095918337763,0.04837022860940616,0.2785072047951641,0.12304743843302654,0.2350930794844956,0.3939359284527848 -2014-10-08,0.25573300814847705,0.0491665094609217,0.28198781605878975,0.12143228062696727,0.24676844885033455,0.4070424755172314 -2014-10-16,0.2517112553253624,0.05007808671782103,0.28589575809121515,0.12132260249019038,0.25720048209273216,0.41601523299222115 -2014-10-24,0.24640355664270414,0.05103371963337279,0.29018712156036525,0.12210862069544276,0.26659152387920154,0.4220379646606709 -2014-11-01,0.23979376294721155,0.05192134188051462,0.2955278441254339,0.1227671080790764,0.2768256565960276,0.42826656154709264 -2014-11-09,0.23177938864829856,0.05262496673605761,0.30295798132302837,0.12224870842069412,0.2905519888216652,0.4384570624038791 -2014-11-17,0.22202955819961262,0.05304414427288477,0.3142440959722114,0.11970634909871596,0.31174678945417883,0.4569531052659349 -2014-11-25,0.2102307737470608,0.05309916099731016,0.33102640732998984,0.11455398206713033,0.34406615677867564,0.4874239517366965 -2014-12-03,0.19628774033966218,0.052767384798209124,0.354020069029411,0.10691978995307355,0.38918412572983097,0.5309145568184324 -2014-12-11,0.18039810908230675,0.05205457513686981,0.3829623750349041,0.09728838673252892,0.4468236230290082,0.5864950654725254 -2014-12-19,0.1635675955406397,0.05101980047577779,0.41503152426990664,0.08681490452551179,0.5116958326682622,0.6485930552774678 -2014-12-27,0.14716695979387556,0.04973814745888745,0.4463632925452314,0.07681882448413833,0.57646180390035,0.7099176152339122 -2015-01-04,0.13319291441757475,0.048299815270835884,0.47149016475759037,0.06868041718605511,0.6306063512707716,0.7608321483845738 -2015-01-12,0.12360059844308971,0.04678747436825057,0.48536593367408043,0.06345805477873405,0.6642453292076048,0.7926363135679714 -2015-01-20,0.12010195183783265,0.04527026397882511,0.4839516421854803,0.06179316425516848,0.6692353435009427,0.7984978226087743 -2015-01-28,0.12340616627323969,0.043793547330879555,0.46611422754378007,0.0637348763699178,0.6428122943488305,0.7763102741667809 -2015-02-05,0.13302022922892107,0.04237499340050736,0.4341588705818192,0.06865124801913729,0.5885956104713781,0.7295457595720422 -2015-02-13,0.14734580020397842,0.041012229309557766,0.3929558104441157,0.07549956871913982,0.5152210327407428,0.6659564918807609 -2015-02-21,0.16356615075480216,0.03968113612975411,0.3501232228711259,0.08280647619858876,0.43674290105944313,0.597912853455777 -2015-03-01,0.17868934620293944,0.03835217923299197,0.3128226999076722,0.08917130369898753,0.36685524285589854,0.5375902043302571 -2015-03-09,0.1896432160968607,0.03699145888877982,0.28735921522298535,0.09332786449669105,0.3182232662748925,0.4964237786451656 -2015-03-17,0.1943555149612242,0.03557670381400324,0.2768142033143582,0.09449099207170432,0.2978356314574223,0.481132309279619 -2015-03-25,0.19198472979516099,0.03410050152495415,0.28043677532260414,0.09245028960574496,0.30587525299425217,0.49275699833026765 -2015-04-02,0.18302767647548324,0.032579258665774644,0.2946289887143264,0.08751520682337746,0.3369063848513542,0.5272741943443697 -2015-04-10,0.1694854721768942,0.031057179375492926,0.3127361104815966,0.08057198846786265,0.379370629766679,0.5750802488513563 -2015-04-18,0.15392233396695904,0.029602194114430783,0.32838027529219294,0.07271199118785253,0.42131730112379623,0.6254831868677837 -2015-04-26,0.1394343260208111,0.028308380431839277,0.33583756349060634,0.06521892445777942,0.4509489232668029,0.6670385393646447 -2015-05-04,0.12856264312333993,0.02727710124578175,0.3322329600434486,0.059191312462347975,0.46100474493672833,0.6915566336850021 -2015-05-12,0.12311960061662011,0.026615698993593756,0.3180778340913917,0.05548255696938244,0.44970946278067075,0.6948791033849775 -2015-05-20,0.12388599116717344,0.02641282541044073,0.2961968301058672,0.05457864453807875,0.4197391172790269,0.6769433226206019 -2015-05-28,0.13046039208497315,0.026733372234951494,0.27187120700041434,0.05654156364877784,0.37861141488988426,0.6422588804162348 -2015-06-05,0.14181668508802037,0.02760552364912903,0.24997782748551836,0.06116086041331577,0.33412846540167973,0.5967596093714802 -2015-06-13,0.1562839449191345,0.029014852843211994,0.23465857905535636,0.06793490544018219,0.29396210903180087,0.5476638827637333 -2015-06-21,0.17228068671064556,0.030914764299435606,0.22788551562290565,0.07631324202764014,0.26277701134919584,0.5006678317656541 -2015-06-29,0.18839401202860445,0.03322363586724639,0.22908608055661206,0.0857131982616931,0.24162811415211463,0.45950838395559224 -2015-07-07,0.20366451833290444,0.03585142465299536,0.23638063254639072,0.09566240260382727,0.22925120015430694,0.4258266176408447 -2015-07-15,0.2176661243505024,0.03870127680601536,0.24656497478876466,0.10582507454606294,0.221939189950495,0.39891760655574815 -2015-07-23,0.23032221431769195,0.04169101800750874,0.2571183789653998,0.11599298561114568,0.21647446985802343,0.3773536784538034 -2015-07-31,0.24193076036645372,0.04475695703271916,0.266464968557068,0.12609989032740213,0.21045017874596672,0.359045615578601 -2015-08-08,0.25286335834660056,0.04785545603930973,0.2742330924643236,0.13613518556408244,0.20299348135381698,0.3423586062850941 diff --git a/inst/mato_grosso_brazil/samples.gpkg b/inst/mato_grosso_brazil/samples.gpkg index 2fa2e37..888cbee 100644 Binary files a/inst/mato_grosso_brazil/samples.gpkg and b/inst/mato_grosso_brazil/samples.gpkg differ diff --git a/inst/reduce_time/patterns/cotton_MODIS13Q1.csv b/inst/reduce_time/patterns/cotton_MODIS13Q1.csv deleted file mode 100644 index ad8d67e..0000000 --- a/inst/reduce_time/patterns/cotton_MODIS13Q1.csv +++ /dev/null @@ -1,24 +0,0 @@ -date,ndvi,evi,red,nir,blue,mir -2001-01-03,0.444967950410115,0.374766005056835,0.127977524042588,0.348767034507488,0.0931126028333047,0.180567671675996 -2001-01-11,0.44311958982467486,0.3727575844039585,0.12828512170265896,0.34236421082028473,0.09247093833676046,0.17711373392112353 -2001-01-20,0.4528943047935253,0.3868140721218037,0.12639918520131932,0.3443528289553863,0.09042335422512508,0.17207438633500732 -2001-01-29,0.47954905797260644,0.42212182637897777,0.12160437335098868,0.35909517081606984,0.08683246838578486,0.1652165788651114 -2001-02-06,0.5252049716567443,0.4778213439933965,0.11350006074692098,0.386574055833113,0.0816772100061517,0.1563289273654477 -2001-02-15,0.587834621504687,0.5476113457525791,0.10222679044592942,0.4227993528028412,0.07510951840631076,0.14532466065506566 -2001-02-24,0.6613212569723727,0.6217219888548638,0.08856983717149404,0.4613187749972771,0.06747995454365656,0.1323994194224809 -2001-03-05,0.7367569672306624,0.6896728241770063,0.07386595308255885,0.4952951110615866,0.05930890889238878,0.11815991321888507 -2001-03-13,0.8045355419376814,0.7426512648511678,0.0597273678232791,0.5192354316834878,0.051206337631188155,0.10365107280850691 -2001-03-22,0.856542712330896,0.7748863054369326,0.04767535195379727,0.5299100336467649,0.04376268987117231,0.09025314846387622 -2001-03-31,0.8877056650068627,0.7838838969284608,0.03881748877657129,0.5264160188133048,0.03744683667641247,0.07947347683852109 -2001-04-09,0.896453505484532,0.76987088139346,0.0336817637739788,0.509700014461407,0.0325435655971044,0.0726942971625978 -2001-04-17,0.8841062361480679,0.7349802022449392,0.0322478419886479,0.48195033295890116,0.02914604377846018,0.07094465402010884 -2001-04-26,0.8536429475591288,0.6826004176614213,0.03412781418705604,0.44613079323896665,0.027195803379201054,0.07474466024730077 -2001-05-05,0.8084992100289524,0.6170278325356408,0.03878916991312941,0.40568961444019197,0.02654527670177526,0.08404061149786808 -2001-05-13,0.7519171467993048,0.5432913448150851,0.04571106331588187,0.36430574237228874,0.02701358372420419,0.09822477900988834 -2001-05-22,0.6869897137680677,0.46690921061751117,0.054420475908401586,0.32552583545229186,0.028416797268485273,0.11622585610897798 -2001-05-31,0.6171208018429609,0.39341389531868143,0.06443232207379304,0.29227589823430467,0.03057166663067295,0.13665919920275993 -2001-06-09,0.546377975615634,0.32766032672225087,0.07517484547920743,0.2663889217463392,0.03328709494925107,0.15802876530309223 -2001-06-17,0.4792889084682046,0.2730983287499943,0.08598559951325613,0.24835434047903487,0.036362152924044155,0.17896248398070724 -2001-06-26,0.4199404774833143,0.23123039005386997,0.09621951669697684,0.23742069178045377,0.03960288137971567,0.1984488751274781 -2001-07-05,0.37079339589385363,0.20146793921652217,0.10541987409401181,0.23200753301277785,0.04285219878988699,0.216002115582154 -2001-07-14,0.331634772021193,0.181355567059534,0.113487387628814,0.230252327695981,0.0460231258415498,0.231730119203362 diff --git a/inst/reduce_time/patterns/maize_MODIS13Q1.csv b/inst/reduce_time/patterns/maize_MODIS13Q1.csv deleted file mode 100644 index 0194889..0000000 --- a/inst/reduce_time/patterns/maize_MODIS13Q1.csv +++ /dev/null @@ -1,24 +0,0 @@ -date,ndvi,evi,red,nir,blue,mir -2001-02-20,0.462427363481723,0.311468097169514,0.105724165720738,0.290319002761205,0.0745369829577744,0.165140392344011 -2001-02-25,0.47125355746184183,0.3166962192844462,0.10232763667644283,0.2862547222377701,0.07214226078879223,0.16360998446659547 -2001-03-03,0.4972289300185806,0.34102884775668374,0.09662078043423508,0.29245070057858347,0.0682655583727576,0.15870292458399596 -2001-03-09,0.5369907259096383,0.38046439139402516,0.08924651776801391,0.30722744777757144,0.06323956827324381,0.15087062580076602 -2001-03-15,0.5867214167544808,0.43047036469167754,0.08087479447887254,0.32855068088668127,0.05742141789650735,0.1406675360099526 -2001-03-21,0.6414169658432533,0.48513635909404146,0.07223198710798212,0.35343613807482027,0.05122663756078658,0.12892811253905226 -2001-03-26,0.6959045710106386,0.538399876868478,0.06396902405167029,0.37862086601586503,0.04504742027908625,0.11659348191541159 -2001-04-01,0.7453138462468007,0.5846347775737952,0.05661391250352754,0.4009111936295438,0.03922135911747632,0.10463090866781981 -2001-04-07,0.7854351361685876,0.6191045872515706,0.05054659534620658,0.4174709580497008,0.03401267167553239,0.09396835658787607 -2001-04-13,0.8131434223858617,0.6384977615490575,0.04598052737152547,0.426182858359981,0.02959492395872,0.08541021636094934 -2001-04-19,0.82659042681737,0.6411661464499067,0.04298589353564932,0.42586821630366717,0.0260568591614169,0.07957967337873814 -2001-04-25,0.825203244437692,0.627106366988472,0.0415257644187671,0.416341899963063,0.0234178073510757,0.0768927489458227 -2001-04-30,0.8096374955828157,0.597875393454268,0.04149083972964175,0.3984159684112245,0.021643436672267827,0.07754078109180956 -2001-05-06,0.78161593080625,0.5563387827464932,0.04274115347140408,0.3737942077002905,0.02066658294362743,0.08148721467833592 -2001-05-12,0.7435931357409921,0.5061845470737054,0.0451394520404325,0.3447819285360683,0.02040730209424049,0.0885024719492727 -2001-05-18,0.6984644741201409,0.4514978242821856,0.04856717076719004,0.31400864739514367,0.02078433503364764,0.09820023812715498 -2001-05-24,0.6492713149060283,0.3963228316303556,0.05293226789061162,0.28412757935509697,0.021722788559789038,0.11007633333987754 -2001-05-29,0.598885448652584,0.3442052781007675,0.058162979440454246,0.2574768022946922,0.02315571627412572,0.12356147996306348 -2001-06-04,0.5497883829978768,0.2979049689707265,0.06419246981112835,0.23584772668606419,0.02502037112115918,0.1380733558875628 -2001-06-10,0.5039482785686681,0.2592497375435074,0.07094599764249833,0.22035917843902117,0.02725424318349895,0.15305196703428736 -2001-06-16,0.46272651829787836,0.22904944912951808,0.07832748234428298,0.21137481948518425,0.029790424272281894,0.16799515559785874 -2001-06-22,0.4267243653338272,0.2069942562969657,0.08620804077681277,0.20842023271014157,0.03255456263890864,0.18251981601564243 -2001-06-28,0.396283797982961,0.19240034761769,0.0944473109047571,0.210740708687293,0.0354696773659327,0.196285828408661 diff --git a/inst/reduce_time/patterns/soybean_MODIS13Q1.csv b/inst/reduce_time/patterns/soybean_MODIS13Q1.csv deleted file mode 100644 index ff318b1..0000000 --- a/inst/reduce_time/patterns/soybean_MODIS13Q1.csv +++ /dev/null @@ -1,24 +0,0 @@ -date,ndvi,evi,red,nir,blue,mir -2000-10-15,0.365295460859194,0.220981257645001,0.111894694872652,0.233946822206083,0.064580106907113,0.226400944649906 -2000-10-21,0.37409688589362905,0.23453494879036565,0.11237163767261167,0.24496255484650878,0.07054532406216199,0.22122346282336844 -2000-10-27,0.3938875377482366,0.2598725937146358,0.1112409622693289,0.2619579300264877,0.07528668134557456,0.21380645982131496 -2000-11-02,0.42648589601610293,0.29965534831520635,0.10831237295698622,0.2861079489636935,0.07823484428595609,0.20395123649009067 -2000-11-08,0.47260218350032934,0.35515320768270603,0.10352369179839634,0.3178942017741009,0.07895274882437292,0.19167965873869341 -2000-11-14,0.5306187407641799,0.4246706832048077,0.09707878024916958,0.35633199881310096,0.07731072482801925,0.17745937053634037 -2000-11-21,0.5967147913812151,0.5036085510254732,0.0894276182259448,0.39902925177664844,0.07354234415034622,0.1621362638207481 -2000-11-27,0.6651962194758281,0.5848864270733645,0.08123667566286298,0.44240903154619177,0.06823856979330911,0.14684401141048103 -2000-12-03,0.7293591244805184,0.6600923599575469,0.07330310627809536,0.4822820290602561,0.062249542822957,0.13282486581413197 -2000-12-09,0.7825236021829289,0.7208710488305187,0.06644911287378416,0.5145277732264604,0.05654351598796634,0.12123265956334842 -2000-12-15,0.8191417229705451,0.7604157265547651,0.06140457065523071,0.5358175466476227,0.05203539141279634,0.11293800987508722 -2000-12-22,0.8358009939775857,0.7748042970794984,0.05869112513050826,0.5442501761637007,0.0494069924916686,0.10837434631542692 -2000-12-28,0.8316618160020145,0.7635457647477194,0.058554995883485214,0.5395958782169152,0.04899363057568795,0.10750465636986029 -2001-01-03,0.8084111026604361,0.7294529223308609,0.06094622296149683,0.5232148244169779,0.050744546995138175,0.10987763870544076 -2001-01-09,0.7698309684738202,0.6779802627848236,0.06553702184742027,0.4977217525077458,0.05424726692631404,0.11475018170579324 -2001-01-15,0.7209948935606467,0.6160446949920498,0.071780765101067,0.4664111994139621,0.05882121268548136,0.12126749107487562 -2001-01-21,0.6672663837096635,0.5506062476024197,0.07900789918405153,0.43258806888389373,0.06367218203613158,0.12864642020175796 -2001-01-28,0.613526308348935,0.48760312814105455,0.08651599768288241,0.39908537311110615,0.06803839968170146,0.13628721627675403 -2001-02-03,0.563679076152823,0.43131183073733287,0.09364614584980943,0.36799528550707156,0.07131042009790925,0.14381057133491432 -2001-02-09,0.5204597183703794,0.38416380927837396,0.09984164860871392,0.3406237653474741,0.07311594365924114,0.1510204185968703 -2001-02-15,0.4856355520518031,0.34714344610003467,0.10466975283599628,0.31770823866116277,0.0733338621068568,0.15779799754644994 -2001-02-21,0.46039099687944673,0.32057208845961077,0.10785302650604524,0.29979996151419147,0.07207619381460215,0.16391451487020098 -2001-02-28,0.445711446516946,0.304612431025669,0.109192892489924,0.287436407004264,0.0695601132749626,0.169053416214151 diff --git a/inst/reduce_time/ts_MODIS13Q1.csv b/inst/reduce_time/ts_MODIS13Q1.csv deleted file mode 100644 index d2e8725..0000000 --- a/inst/reduce_time/ts_MODIS13Q1.csv +++ /dev/null @@ -1,94 +0,0 @@ -date,ndvi,evi,red,nir,blue,mir -2009-08-05,0.3169,0.1687,0.1167,0.225,0.0427,0.2193 -2009-08-28,0.2609,0.1385,0.1168,0.1993,0.0548,0.2657 -2009-09-13,0.1794,0.0809,0.1219,0.1752,0.0348,0.276 -2009-09-28,0.2785,0.1617,0.1308,0.2318,0.0608,0.2306 -2009-10-01,0.1769,0.0993,0.1698,0.2428,0.0565,0.315 -2009-10-31,0.6701,0.4491,0.0674,0.3413,0.0295,0.1442 -2009-11-16,0.7259,0.7358,0.0861,0.5422,0.0679,0.1245 -2009-12-02,0.9388,0.9227,0.0198,0.628,0.045,0.1292 -2009-12-14,0.7379,0.6383,0.0693,0.4597,0.1409,0.1007 -2009-12-25,0.9125,0.748,0.0219,0.4788,0.0111,0.08 -2010-01-05,0.6576,0.672,0.1183,0.5729,0.2465,0.263 -2010-01-18,0.3707,0.4547,0.1868,0.4069,0.1757,0.1892 -2010-02-03,0.2685,0.1622,0.1494,0.2591,0.062,0.2978 -2010-03-05,0.5466,0.3843,0.1022,0.3487,0.0478,0.1781 -2010-03-17,0.7121,0.5757,0.0737,0.4384,0.0396,0.1678 -2010-03-29,0.8159,0.6076,0.0452,0.446,0.0091,0.0684 -2010-04-09,0.8796,0.73740000000000006,0.0306,0.478,0.0193,0.0754 -2010-04-25,0.9009,0.8169,0.0283,0.5429,0.0184,0.062 -2010-05-20,0.8231,0.6455,0.0414,0.4267,0.0244,0.0814 -2010-05-27,0.7524,0.4739,0.0454,0.3214,0.0184,0.0805 -2010-06-12,0.6657,0.4061,0.0591,0.2945,0.0267,0.1183 -2010-06-28,0.49,0.301,0.0911,0.2662,0.0478,0.1538 -2010-07-21,0.389,0.2342,0.1075,0.2444,0.0571,0.1692 -2010-07-30,0.339,0.2106,0.1239,0.251,0.0648,0.1896 -2010-08-22,0.2862,0.1853,0.1591,0.2867,0.0693,0.2756 -2010-09-09,0.273,0.1873,0.1641,0.2874,0.0835,0.3027 -2010-09-16,0.2788,0.1765,0.1399,0.2481,0.074,0.276 -2010-10-09,0.2755,0.1686,0.1391,0.2449,0.0682,0.2755 -2010-10-20,0.211,0.1157,0.1428,0.2192,0.0567,0.2105 -2010-11-03,0.2755,0.1467,0.1308,0.2303,0.0427,0.2735 -2010-11-24,0.86,0.581,0.0265,0.3523,0.0146,0.0845 -2010-12-16,0.92,0.9011,0.026,0.6241,0.0161,0.1069 -2010-12-25,0.9351,0.9139,0.0212,0.6328,0.0116,0.0943 -2011-01-10,0.882,0.8581,0.0378,0.6031,0.0244,0.1223 -2011-01-31,0.3246,0.1792,0.1146,0.2248,0.0501,0.2465 -2011-02-14,0.2693,0.2263,0.185,0.3214,0.1362,0.221 -2011-02-23,0.2653,0.1654,0.1221,0.2103,0.0764,0.1717 -2011-03-18,0.3634,0.2891,0.1485,0.3181,0.0863,0.218 -2011-03-26,0.6881,0.5454,0.0724,0.3919,0.0407,0.1707 -2011-04-19,0.9062,0.6599,0.0197,0.4004,0.0102,0.0546 -2011-05-05,0.9096,0.8403,0.0265,0.5599,0.0176,0.0588 -2011-05-14,0.9063,0.8416,0.0286,0.5825,0.0145,0.065 -2011-05-30,0.878,0.7637,0.0329,0.5068,0.0204,0.0589 -2011-06-24,0.7817,0.5605,0.0466,0.3805,0.0228,0.0813 -2011-06-29,0.7457,0.5012,0.048,0.3296,0.0284,0.1029 -2011-07-15,0.5194,0.321,0.0833,0.2634,0.0481,0.1526 -2011-08-02,0.4033,0.2544,0.1114,0.262,0.0601,0.1709 -2011-08-18,0.3519,0.2277,0.1273,0.2656,0.0682,0.1958 -2011-09-03,0.3623,0.2233,0.1276,0.2726,0.0554,0.2733 -2011-09-19,0.3244,0.2097,0.1415,0.2774,0.0676,0.272 -2011-10-04,0.3296,0.2657,0.1595,0.3164,0.0752,0.334 -2011-10-27,0.2434,0.1737,0.1512,0.2485,0.1521,0.2697 -2011-11-13,0.4236,0.3227,0.1263,0.312,0.0674,0.2237 -2011-11-29,0.8641,0.6987,0.0337,0.4625,0.0174,0.1058 -2011-12-03,0.8935,0.8458,0.0334,0.5941,0.0183,0.1347 -2011-12-27,0.9041,0.8743,0.0252,0.5006,0.039,0.0672 -2012-01-09,0.8229,0.8991,0.0687,0.7075,0.0858,0.1228 -2012-01-28,0.5187,0.3247,0.0804,0.2537,0.0243,0.1249 -2012-02-07,0.2889,0.2382,0.175,0.3172,0.0639,0.2641 -2012-03-04,0.3884,0.2,0.0966,0.2193,0.0354,0.2113 -2012-03-13,0.3159,0.2237,0.18,0.3463,0.0758,0.165 -2012-03-22,0.8541,0.6418,0.0322,0.4094,0.0178,0.0609 -2012-04-07,0.8614,0.6827,0.0335,0.45,0.0168,0.0536 -2012-04-23,0.7968,0.5966,0.0461,0.4078,0.0225,0.0694 -2012-05-09,0.8309,0.6055,0.0361,0.391,0.019,0.0607 -2012-06-01,0.7016,0.3898,0.0469,0.2675,0.0179,0.0842 -2012-06-10,0.528,0.2432,0.0602,0.1949,0.0229,0.127 -2012-07-03,0.3303,0.1096,0.0589,0.117,0.0194,0.1637 -2012-07-19,0.3243,0.1171,0.0675,0.1323,0.0206,0.1701 -2012-08-04,0.322,0.1365,0.08,0.156,0.0326,0.1939 -2012-08-20,0.2579,0.1838,0.2229,0.3779,0.081,0.3034 -2012-09-05,0.2489,0.184,0.2055,0.3417,0.0966,0.2673 -2012-09-14,0.2675,0.2052,0.2092,0.362,0.1008,0.3167 -2012-10-06,0.3047,0.227,0.1476,0.277,0.0777,0.2942 -2012-10-29,0.3206,0.2337,0.1398,0.2718,0.0722,0.2951 -2012-11-01,0.276,0.2122,0.1608,0.2834,0.0726,0.3471 -2012-11-30,0.1941,0.193,0.2661,0.3943,0.2402,0.1671 -2012-12-07,0.5567,0.4362,0.1012,0.3554,0.0831,0.209 -2012-12-24,0.8884,0.7403,0.0287,0.486,0.0152,0.0868 -2013-01-09,0.9436,0.7865,0.0147,0.5075,0.0039,0.0715 -2013-01-17,0.942,0.9185,0.0185,0.6205,0.0124,0.108 -2013-02-15,0.5129,0.2975,0.0736,0.2286,0.0589,0.1131 -2013-02-19,0.3488,0.1924,0.0922,0.191,0.0337,0.2226 -2013-03-12,0.3518,0.1964,0.0932,0.1944,0.0439,0.1962 -2013-04-01,0.8148,0.6866,0.0445,0.4363,0.0369,0.0766 -2013-04-17,0.8964,0.6713,0.0226,0.4137,0.0124,0.0495 -2013-05-08,0.8669,0.5752,0.0244,0.3425,0.0142,0.0599 -2013-05-19,0.8378,0.5342,0.0285,0.3231,0.0154,0.0654 -2013-05-26,0.8295,0.5434,0.0307,0.3296,0.0185,0.0646 -2013-06-13,0.7172,0.4156,0.0457,0.2775,0.021,0.1003 -2013-06-29,0.3458,0.122,0.0627,0.129,0.0196,0.1603 -2013-07-15,0.3208,0.1163,0.0671,0.1305,0.0227,0.1745 -2013-07-31,0.317,0.1309,0.0827,0.1595,0.0253,0.2039 diff --git a/man/create_patterns.Rd b/man/create_patterns.Rd deleted file mode 100644 index 0ff9bd5..0000000 --- a/man/create_patterns.Rd +++ /dev/null @@ -1,41 +0,0 @@ -% Generated by roxygen2: do not edit by hand -% Please edit documentation in R/create_patterns.R -\name{create_patterns} -\alias{create_patterns} -\title{Create a Pattern Using GAM} -\usage{ -create_patterns( - x, - y, - formula = band ~ s(time), - start_column = "start_date", - end_column = "end_date", - label_colum = "label", - sampling_freq = NULL, - ... -) -} -\arguments{ -\item{x}{A three dimensions stars object (x, y, time) with the satellite image time series.} - -\item{y}{An sf object with the coordinates of the training points.} - -\item{formula}{A formula for the GAM. Default is \code{band ~ \link[mgcv]{s}(time)}.} - -\item{start_column}{Name of the column in y that indicates the start date. Default is 'start_date'.} - -\item{end_column}{Name of the column in y that indicates the end date. Default is 'end_date'.} - -\item{label_colum}{Name of the column in y that contains land use labels. Default is 'label'.} - -\item{sampling_freq}{The time frequency for sampling including unit, e.g '16 day'. If NULL, the function will infer it.} - -\item{...}{Additional arguments passed to the GAM function.} -} -\value{ -A list containing the predicted values for each label. -} -\description{ -This function creates a pattern based on Generalized Additive Models (GAM). -It uses the specified formula to fit the model and predict values. -} diff --git a/man/extract_time_series.Rd b/man/extract_time_series.Rd deleted file mode 100644 index c4aebe3..0000000 --- a/man/extract_time_series.Rd +++ /dev/null @@ -1,21 +0,0 @@ -% Generated by roxygen2: do not edit by hand -% Please edit documentation in R/create_patterns.R -\name{extract_time_series} -\alias{extract_time_series} -\title{Extract Time Series from a Stars Object for Specified Points} -\usage{ -extract_time_series(x, y) -} -\arguments{ -\item{x}{A stars object containing the raster time series data.} - -\item{y}{An sf object containing the point geometries and their associated labels.} -} -\value{ -A data.frame with the extracted time series for each point in the sf object, -with additional columns for the ID and label of each sample. -} -\description{ -This function extracts a time series from a stars object for each specified point in the sf object. -Each extracted sample is then labeled with an ID and the label from the sf object. -} diff --git a/man/get_stars_time_freq.Rd b/man/get_stars_time_freq.Rd deleted file mode 100644 index 47d90c4..0000000 --- a/man/get_stars_time_freq.Rd +++ /dev/null @@ -1,18 +0,0 @@ -% Generated by roxygen2: do not edit by hand -% Please edit documentation in R/create_patterns.R -\name{get_stars_time_freq} -\alias{get_stars_time_freq} -\title{Compute the Most Common Sampling Frequency in a Stars Object} -\usage{ -get_stars_time_freq(x) -} -\arguments{ -\item{x}{A stars object containing time series data.} -} -\value{ -A difftime object representing the most common time difference between consecutive samples. -} -\description{ -This function calculates the most common difference between consecutive time points in a stars object. -This can be useful for determining the sampling frequency of the time series data. -} diff --git a/man/get_time_series_freq.Rd b/man/get_time_series_freq.Rd new file mode 100644 index 0000000..02dea06 --- /dev/null +++ b/man/get_time_series_freq.Rd @@ -0,0 +1,18 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/train.R +\name{get_time_series_freq} +\alias{get_time_series_freq} +\title{Compute the Most Common Sampling Frequency across all observations} +\usage{ +get_time_series_freq(x) +} +\arguments{ +\item{x}{A data frame including a column called `observations`` with the time series} +} +\value{ +A difftime object representing the most common time difference between consecutive samples. +} +\description{ +This function calculates the most common difference between consecutive time points. +This can be useful for determining the aproximate sampling frequency of the time series data. +} diff --git a/man/plot.twdtw_knn1.Rd b/man/plot.twdtw_knn1.Rd new file mode 100644 index 0000000..866d3d7 --- /dev/null +++ b/man/plot.twdtw_knn1.Rd @@ -0,0 +1,75 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/plot.R +\name{plot.twdtw_knn1} +\alias{plot.twdtw_knn1} +\title{Plot Patterns from twdtw-knn1 model} +\usage{ +\method{plot}{twdtw_knn1}(x, bands = NULL, ...) +} +\arguments{ +\item{x}{A model of class \code{"twdtw_knn1"}.} + +\item{bands}{A character vector specifying the bands or indices to plot. +If NULL (default), all available bands or indices in the data will be plotted.} + +\item{...}{Additional arguments passed to \code{\link[ggplot2]{ggplot}}. Currently not used.} +} +\value{ +A \code{\link[ggplot2]{ggplot}} object displaying the time series patterns. +} +\description{ +This function visualizes time series patterns from the \code{"twdtw_knn1"} model. +It produces a multi-faceted plot, where each facet represents a different time series +label from the model's data. Within each facet, different bands or indices (attributes) +are plotted as distinct lines, differentiated by color. +} +\examples{ +\dontrun{ + +# Read training samples +samples <- + system.file("mato_grosso_brazil/samples.gpkg", package = "dtwSat") |> + st_read(quiet = TRUE) + +# Get satellite image time sereis files +tif_files <- system.file("mato_grosso_brazil", package = "dtwSat") |> + dir(pattern = "\\\\.tif$", full.names = TRUE) + +# Get acquisition dates +acquisition_date <- regmatches(tif_files, regexpr("[0-9]{8}", tif_files)) |> + as.Date(format = "\%Y\%m\%d") + +# Create a 3D datacube +dc <- read_stars(tif_files, + proxy = FALSE, + along = list(time = acquisition_date), + RasterIO = list(bands = 1:6)) |> + st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")) |> + split(c("band")) + +# Create a knn1-twdtw model +m <- twdtw_knn1(x = dc, + y = samples, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50), + formula = band ~ s(time)) + +print(m) + +# Visualize model patterns +plot(m) + +# Classify satellite images +system.time(lu <- predict(dc, model = m)) + +# Visualise land use classification +ggplot() + + geom_stars(data = lu) + + theme_minimal() + +} +} +\seealso{ +twdtw_knn1 +} diff --git a/man/plot_patterns.Rd b/man/plot_patterns.Rd deleted file mode 100644 index 97e24d6..0000000 --- a/man/plot_patterns.Rd +++ /dev/null @@ -1,26 +0,0 @@ -% Generated by roxygen2: do not edit by hand -% Please edit documentation in R/plot_patterns.R -\name{plot_patterns} -\alias{plot_patterns} -\title{Plot Patterns from Time Series Data} -\usage{ -plot_patterns(x, ...) -} -\arguments{ -\item{x}{A list where each element is a data.frame representing a time series. -Each data.frame should have the same number of rows and columns, -with columns representing different attributes (e.g., bands or indices) -and rows representing time points. -The name of each element in the list will be used as the facet title.} - -\item{...}{Not used.} -} -\value{ -A ggplot object displaying the time series patterns. -} -\description{ -This function takes a list of time series data and creates a multi-faceted plot -where each facet corresponds to a different time series from the list. -Within each facet, different attributes (columns of the time series) are -plotted as lines with different colors. -} diff --git a/man/predict.twdtw_knn1.Rd b/man/predict.twdtw_knn1.Rd new file mode 100644 index 0000000..d3d0083 --- /dev/null +++ b/man/predict.twdtw_knn1.Rd @@ -0,0 +1,75 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/predict.R +\name{predict.twdtw_knn1} +\alias{predict.twdtw_knn1} +\title{Predict using the twdtw_knn1 model} +\usage{ +\method{predict}{twdtw_knn1}(object, newdata, ...) +} +\arguments{ +\item{object}{A \code{twdtw_knn1} model object generated by the \code{twdtw_knn1} function.} + +\item{newdata}{A data frame or similar object containing the new observations +(time series data) to be predicted.} + +\item{...}{Additional arguments passed to the \link[twdtw]{twdtw} function. +If provided, they will overwrite twdtw arguments previously passed to \link[dtwSat]{twdtw_knn1}.} +} +\value{ +A vector of predicted classes for the \code{newdata}. +} +\description{ +This function predicts the classes of new data using the Time Warped Dynamic Time Warping (TWDTW) +method with a 1-nearest neighbor approach. The prediction is based on the minimum TWDTW distance +to the known patterns stored in the \code{twdtw_knn1} model. +} +\examples{ +\dontrun{ + +# Read training samples +samples <- + system.file("mato_grosso_brazil/samples.gpkg", package = "dtwSat") |> + st_read(quiet = TRUE) + +# Get satellite image time sereis files +tif_files <- system.file("mato_grosso_brazil", package = "dtwSat") |> + dir(pattern = "\\\\.tif$", full.names = TRUE) + +# Get acquisition dates +acquisition_date <- regmatches(tif_files, regexpr("[0-9]{8}", tif_files)) |> + as.Date(format = "\%Y\%m\%d") + +# Create a 3D datacube +dc <- read_stars(tif_files, + proxy = FALSE, + along = list(time = acquisition_date), + RasterIO = list(bands = 1:6)) |> + st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")) |> + split(c("band")) + +# Create a knn1-twdtw model +m <- twdtw_knn1(x = dc, + y = samples, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50), + formula = band ~ s(time)) + +print(m) + +# Visualize model patterns +plot(m) + +# Classify satellite images +system.time(lu <- predict(dc, model = m)) + +# Visualise land use classification +ggplot() + + geom_stars(data = lu) + + theme_minimal() + +} +} +\seealso{ +twdtw_knn1 +} diff --git a/man/prepare_time_series.Rd b/man/prepare_time_series.Rd new file mode 100644 index 0000000..dfce70e --- /dev/null +++ b/man/prepare_time_series.Rd @@ -0,0 +1,20 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/prepare_time_series.R +\name{prepare_time_series} +\alias{prepare_time_series} +\title{Prepare a Time Series Tibble from a 2D stars Object with Bands and Time Attributes} +\usage{ +prepare_time_series(x) +} +\arguments{ +\item{x}{A data frame derived from a stars object containing time series data in wide format. +The column names should adhere to one of the following formats: "band.YYYY.MM.DD", "XYYYY.MM.DD.band", or "YYYY.MM.DD.band".} +} +\value{ +A nested tibble in wide format. Each row of the tibble corresponds to a unique 'ts_id' that maintains the order from the original stars object. +The nested structure contains observations (time series) for each 'ts_id', including the 'time' of each observation, and individual bands are presented as separate columns. +} +\description{ +This function reshapes a data frame, which has been converted from a stars object, into a nested wide tibble format. +The stars object conversion often results in columns named in formats like "band.YYYY.MM.DD", "XYYYY.MM.DD.band", or "YYYY.MM.DD.band". +} diff --git a/man/pretty_arguments.Rd b/man/pretty_arguments.Rd new file mode 100644 index 0000000..a620918 --- /dev/null +++ b/man/pretty_arguments.Rd @@ -0,0 +1,23 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/train.R +\name{pretty_arguments} +\alias{pretty_arguments} +\title{Print Pretty Arguments} +\usage{ +pretty_arguments(args) +} +\arguments{ +\item{args}{A list of named arguments to display.} +} +\value{ +Invisible NULL. The function is mainly used for its side effect of printing. +} +\description{ +Display a list of arguments of a given function in a human-readable format. +} +\examples{ +\dontrun{ +pretty_arguments(formals(twdtw_knn1)) +} + +} diff --git a/man/print.twdtw_knn1.Rd b/man/print.twdtw_knn1.Rd new file mode 100644 index 0000000..c909f7a --- /dev/null +++ b/man/print.twdtw_knn1.Rd @@ -0,0 +1,20 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/train.R +\name{print.twdtw_knn1} +\alias{print.twdtw_knn1} +\title{Print method for objects of class twdtw_knn1} +\usage{ +\method{print}{twdtw_knn1}(x, ...) +} +\arguments{ +\item{x}{An object of class \code{twdtw_knn1}.} + +\item{...}{ignored} +} +\value{ +Invisible \code{twdtw_knn1} object. +} +\description{ +This method provides a structured printout of the important components +of a \code{twdtw_knn1} object. +} diff --git a/man/shift_dates.Rd b/man/shift_dates.Rd index 87d58c1..9af6b7e 100644 --- a/man/shift_dates.Rd +++ b/man/shift_dates.Rd @@ -1,5 +1,5 @@ % Generated by roxygen2: do not edit by hand -% Please edit documentation in R/utils.R +% Please edit documentation in R/shift_dates.R \name{shift_dates} \alias{shift_dates} \title{Shift Dates to Start on a Specified Origin Year} diff --git a/man/twdtw_knn1.Rd b/man/twdtw_knn1.Rd new file mode 100644 index 0000000..569963e --- /dev/null +++ b/man/twdtw_knn1.Rd @@ -0,0 +1,105 @@ +% Generated by roxygen2: do not edit by hand +% Please edit documentation in R/train.R +\name{twdtw_knn1} +\alias{twdtw_knn1} +\title{Train a KNN-1 TWDTW model with optional GAM resampling} +\usage{ +twdtw_knn1( + x, + y, + time_weight, + cycle_length, + time_scale, + formula = NULL, + start_column = "start_date", + end_column = "end_date", + label_colum = "label", + sampling_freq = NULL, + ... +) +} +\arguments{ +\item{x}{A three-dimensional stars object (x, y, time) with bands as attributes.} + +\item{y}{An sf object with the coordinates of the training points.} + +\item{time_weight}{A numeric vector with length two (steepness and midpoint of logistic weight) or a function. +See details in \link[twdtw]{twdtw}.} + +\item{cycle_length}{The length of the cycle, e.g. phenological cycles. Details in \link[twdtw]{twdtw}.} + +\item{time_scale}{Specifies the time scale for the observations. Details in \link[twdtw]{twdtw}.} + +\item{formula}{Either NULL or a formula to reduce samples of the same label using Generalized Additive Models (GAM). +Default is \code{band ~ s(time)}. See details.} + +\item{start_column}{Name of the column in y that indicates the start date. Default is 'start_date'.} + +\item{end_column}{Name of the column in y that indicates the end date. Default is 'end_date'.} + +\item{label_colum}{Name of the column in y containing land use labels. Default is 'label'.} + +\item{sampling_freq}{The time frequency for sampling, including the unit (e.g., '16 day'). +If NULL, the function will infer the frequency. This parameter is only used if a formula is provided.} + +\item{...}{Additional arguments passed to the \link[mgcv]{gam} function and to \link[twdtw]{twdtw} function.} +} +\value{ +A 'twdtw_knn1' model containing the trained model information and the data used. +} +\description{ +This function prepares a KNN-1 model with the Time Warp Dynamic Time Warping (TWDTW) algorithm. +If a formula is provided, the training samples are resampled using Generalized Additive Models (GAM). +} +\details{ +If \code{formula} is NULL, the KNN-1 model will retain all training samples. If a formula is passed (e.g., \code{band ~ \link[mgcv]{s}(time)}), +then samples of the same label (land cover class) will be resampled using GAM. +Resampling can significantly reduce prediction processing time. +} +\examples{ +\dontrun{ + +# Read training samples +samples <- + system.file("mato_grosso_brazil/samples.gpkg", package = "dtwSat") |> + st_read(quiet = TRUE) + +# Get satellite image time sereis files +tif_files <- system.file("mato_grosso_brazil", package = "dtwSat") |> + dir(pattern = "\\\\.tif$", full.names = TRUE) + +# Get acquisition dates +acquisition_date <- regmatches(tif_files, regexpr("[0-9]{8}", tif_files)) |> + as.Date(format = "\%Y\%m\%d") + +# Create a 3D datacube +dc <- read_stars(tif_files, + proxy = FALSE, + along = list(time = acquisition_date), + RasterIO = list(bands = 1:6)) |> + st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")) |> + split(c("band")) + +# Create a knn1-twdtw model +m <- twdtw_knn1(x = dc, + y = samples, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50), + formula = band ~ s(time)) + +print(m) + +# Visualize model patterns +plot(m) + +# Classify satellite images +system.time(lu <- predict(dc, model = m)) + +# Visualise land use classification +ggplot() + + geom_stars(data = lu) + + theme_minimal() + +} +} diff --git a/tests/testthat/test-twdtw_classify.R b/tests/testthat/test-twdtw_classify.R index f5c0c04..8af1da7 100644 --- a/tests/testthat/test-twdtw_classify.R +++ b/tests/testthat/test-twdtw_classify.R @@ -1,27 +1,89 @@ -library(stars) -library(stringr) - # Read training samples -samples <- st_read(system.file("mato_grosso_brazil/samples.gpkg", package = "dtwSat")) +samples <- system.file("mato_grosso_brazil/samples.gpkg", package = "dtwSat") |> + st_read(quiet = TRUE) # Satellite image time sereis files tif_files <- system.file("mato_grosso_brazil", package = "dtwSat") |> dir(pattern = "\\.tif$", full.names = TRUE) -# The acquisition date is in the file name are not the true acquisition date of each pixel -# MOD13Q1 is a 16-day composite product, so the acquisition date is the first day of the 16-day period -acquisition_date <- as.Date(str_extract(tif_files, "[0-9]{8}"), format = "%Y%m%d") +# The acquisition date is in the file name are not the true acquisition date +# of each pixel. MOD13Q1 is a 16-day composite product, so the acquisition date +# is the first day of the 16-day period +acquisition_date <- regmatches(tif_files, regexpr("[0-9]{8}", tif_files)) |> + as.Date(format = "%Y%m%d") + +# Read the data as a stars object setting the time/date for each observation +# using along. This will prodcue a 4D array (data-cube) which will then be converted +# to a 3D array by spliting the 'band' dimension +dc <- read_stars(tif_files, + proxy = FALSE, + along = list(time = acquisition_date), + RasterIO = list(bands = 1:6)) |> + st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")) |> + split(c("band")) + +# Create a knn1-twdtw model +system.time( + m <- twdtw_knn1(x = dc, + y = samples, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50), + formula = band ~ s(time)) +) + +print(m) + +# Visualize model patterns +plot(m) + +# Classify satellite image time series +system.time(lu <- predict(dc, model = m)) + +# Visualise land use classification +ggplot() + + geom_stars(data = lu) + + theme_minimal() + +### OTHER TESTS +# split time first +dc <- read_stars(tif_files, + proxy = FALSE, + along = list(time = acquisition_date), + RasterIO = list(bands = 1:6)) |> + st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")) |> + split(c("band")) + +m <- twdtw_knn1(x = dc, + y = samples, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50), + formula = band ~ s(time), + sampling_freq = 60) + +plot(m) + +system.time( + lu <- predict(dc, + model = m, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50)) +) -# Read the data as a stars object, setting time as a dimension and band as attribute -dc <- read_stars(tif_files, proxy = FALSE, along = list(time = acquisition_date)) |> - st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR", "DOY")) |> - split("band") +# Visualise land use classification +ggplot() + + geom_stars(data = lu) + + theme_minimal() -# Remove the DOY band - this will be supported int the future -dc <- dc[c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")] +# Test model without samples reduction +m <- twdtw_knn1(x = dc, + y = samples, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50)) -# Get temporal patters -ts_patterns <- create_patterns(x = dc, y = samples) +plot(m) -# Visualize patterns -plot_patterns(ts_patterns) +plot(m, bands = c('EVI', 'NDVI')) diff --git a/vignettes.awk b/vignettes.awk new file mode 100644 index 0000000..0517265 --- /dev/null +++ b/vignettes.awk @@ -0,0 +1,13 @@ +{ + if (NR == 4) { + print("output: rmarkdown::html_vignette") + } else if (NR > 4 && NR <= 10) { + } else if (NR == 16) { + print("\n**For a better version of the dtwSat vignettes see** https://vwmaus.github.io/dtwSat/articles/\n") + } else if (NR == 19) { + print + print("knitr::opts_chunk$set(fig.height = 4.5)") + print("knitr::opts_chunk$set(fig.width = 6)") + } else + print +} diff --git a/vignettes/landuse-mapping.Rmd b/vignettes/landuse-mapping.Rmd new file mode 100644 index 0000000..66d8af0 --- /dev/null +++ b/vignettes/landuse-mapping.Rmd @@ -0,0 +1,137 @@ +--- +title: "1. Land use mapping using TWDTW-1NN and stars" +author: "Victor Maus" +output: + html_document: + toc: true + toc_float: + collapsed: false + smooth_scroll: false + toc_depth: 2 +vignette: > + %\VignetteIndexEntry{1. Land use mapping using TWDTW-1NN and stars} + %\VignetteEncoding{UTF-8} + %\VignetteEngine{knitr::rmarkdown} +bibliography: ./../inst/REFERENCES.bib +--- + + +```{r setup, include=FALSE} +knitr::opts_chunk$set(echo = TRUE, collapse = TRUE, dev = "png") +``` + +This vignette offers a concise guide for using version 1.0.0 or higher of the `dtwSat` package to generate a land-use map. +The package utilizes Time-Weighted Dynamic Time Warping (TWDTW) along with a 1-Nearest Neighbor (1-NN) classifier. +The subsequent sections will walk you through the process of creating a land-use map based on a set of training samples +and a multi-band satellite image time series. + +## Reading Training Samples + +First, let's read a set of training samples that come with the `dtwSat` package installation. + +```{r , echo = TRUE, eval = TRUE, warning = FALSE, message = FALSE} +library(dtwSat) + +samples <- system.file("mato_grosso_brazil/samples.gpkg", package = "dtwSat") |> + st_read(quiet = TRUE) +``` + +## Preparing the Satellite Image Time Series + +The `dtwSat` package supports satellite images read into R using the `stars` package. +The instalation comes with a set of MOD13Q1 images for a region within the Brazilian Amazon. +Note that timing is crucial for the TWDTW distance metric. To create a consistent image time series, +we start by extracting the date of acquisition from the MODIS file names [@Didan:2015]. + +```{r , echo = TRUE, eval = TRUE, warning = FALSE, message = FALSE} +tif_files <- system.file("mato_grosso_brazil", package = "dtwSat") |> + dir(pattern = "\\.tif$", full.names = TRUE) + +acquisition_date <- regmatches(tif_files, regexpr("[0-9]{8}", tif_files)) |> + as.Date(format = "%Y%m%d") + +print(acquisition_date) +``` + +*Side note:* The date in the file name is not the true acquisition date for each pixel. +MOD13Q1 is a 16-day composite product, and the date in the file name is the first day of this 16-day period. + +With the files and dates in hand, we can construct a stars satellite image time series for `dtwSat`. + +```{r , echo = TRUE, eval = TRUE, warning = FALSE, message = FALSE} +dc <- read_stars(tif_files, + proxy = FALSE, + along = list(time = acquisition_date), + RasterIO = list(bands = 1:6)) |> + st_set_dimensions(3, c("EVI", "NDVI", "RED", "BLUE", "NIR", "MIR")) |> + split(c("band")) + +print(dc) +``` + +Note that it's important to set the date for each observation using the parameter `along`. +This will produce a 4D array (data-cube) that will be collapsed into a 3D array by converting +the 'band' dimension into attributes. This prepares the data for training the TWDTW-1NN model. + +## Create TWDTW-KNN1 model + +```{r , echo = TRUE, eval = TRUE, warning = FALSE, message = FALSE} +twdtw_model <- twdtw_knn1(x = dc, + y = samples, + cycle_length = 'year', + time_scale = 'day', + time_weight = c(steepness = 0.1, midpoint = 50), + formula = band ~ s(time)) + +print(twdtw_model) +``` + +In addition to the mandatory arguments `x` (satellite data-cube) and `y` (training samples), +the TWDTW distance calculation also requires setting `cycle_length`, `time_scale`, and `time_weight`. +For more details, refer to the documentation using `?twdtw`. The argument `formula = band ~ s(time)` is optional. +If provided, training samples time sereis are resampled using Generalized Additive Models (GAMs), +collapsing all samples with the same land-use label into a single sample. This reduces computational +demands. The sample in the model can be visualized as follows: + +```{r , echo = TRUE, eval = TRUE, warning = FALSE, message = FALSE} +plot(twdtw_model) +``` + + +## Land Use Prediction + +Finally, we predict the land-use classes for each pixel location in the data-cube: + +```{r , echo = TRUE, eval = TRUE, warning = FALSE, message = FALSE} +lu_map <- predict(dc, model = twdtw_model) +print(lu_map) +``` + +The 'time' dimension was reduced to a single map. We can now visualize it using `ggplot`: + +```{r , echo = TRUE, eval = TRUE, warning = FALSE, message = FALSE} +ggplot() + + geom_stars(data = lu_map) + + theme_minimal() +``` + +Note that some pixels (in a 3x3 box) in the northeast part of the map have `NA` values due +to a missing value in the blue band recorded on 2011-11-17. This limitation will be addressed +in future versions of the dtwSat package. + +Ultimately, we can write the map to a TIFF file we can use: + +```{r , echo = TRUE, eval = FALSE, warning = FALSE, message = FALSE} +write_stars(lu_map, "lu_map.tif") +``` + +## Further Reading + +This introduction outlined the use of `dtwSat` for land-use mapping. +For more in-depth information, refer to the papers by @Maus:2016 and @Maus:2019 and the +[`twdtw` R package documentation](https://cran.r-project.org/web/packages/twdtw). + +For additional details on how to manage input and output satellite images, +[`check` the stars documentation](https://cran.r-project.org/web/packages/stars). + +## References