From a6e20fdeb48f0e437a21ae524fed5a7f94e33796 Mon Sep 17 00:00:00 2001 From: Emma Mendelsohn Date: Thu, 12 Oct 2023 13:45:52 -0400 Subject: [PATCH] modis ndvi dataset remove processing of quality files --- R/create_modis_ndvi_dataset.R | 14 ++++++++------ 1 file changed, 8 insertions(+), 6 deletions(-) diff --git a/R/create_modis_ndvi_dataset.R b/R/create_modis_ndvi_dataset.R index d3f0ee2..0668232 100644 --- a/R/create_modis_ndvi_dataset.R +++ b/R/create_modis_ndvi_dataset.R @@ -10,27 +10,29 @@ #' @return #' @author Emma Mendelsohn #' @export -create_modis_ndvi_dataset <- function(modis_ndvi_downloaded, +create_modis_ndvi_dataset <- function(modis_ndvi_downloaded_subset, continent_raster_template, modis_ndvi_directory_dataset, overwrite = FALSE) { - filename <- basename(modis_ndvi_downloaded) - + # Extract start and end dates from the raw downloaded file name + filename <- basename(modis_ndvi_downloaded_subset) year_doy <- sub(".*doy(\\d+).*", "\\1", filename) start_date <- as.Date(year_doy, format = "%Y%j") # confirmed this is start date through manual download tests end_date <- start_date + 16 - save_filename <- glue::glue("transformed_modis_NDVI_{start_date}_to_{end_date}.gz.parquet") - - existing_files <- list.files(modis_ndvi_directory_dataset) + # Set filename for saving + save_filename <- glue::glue("transformed_modis_NDVI_{start_date}_to_{end_date}.gz.parquet") message(paste0("Transforming ", save_filename)) + # Check if file already exists + existing_files <- list.files(modis_ndvi_directory_dataset) if(save_filename %in% existing_files & !overwrite){ message("file already exists, skipping transform") return(file.path(modis_ndvi_directory_dataset, save_filename)) } + # Transform with template raster transformed_raster <- transform_raster(raw_raster = rast(modis_ndvi_downloaded), template = rast(continent_raster_template))