inst/examples/build_brick_maskcloud.R

In albhasan/sits.starfm: Use the StarFM Fusion Model on Landsat 8 Images

#!/usr/bin/env Rscript
# Build a brick by piling and cloud-masking Landsat 8 images
setwd("/home/alber/Documents/data/experiments/l8mod-fusion/Rpackage/sits.starfm")
suppressMessages(suppressPackageStartupMessages(library(dplyr)))
suppressMessages(suppressPackageStartupMessages(library(optparse)))
suppressMessages(suppressPackageStartupMessages(library(sits.starfm)))

option_list = list(
    make_option(c("-s", "--scene"), type = "character", default = NULL, 
        help = "landsat scene e.g. 225063",        metavar = "character"),
    make_option(c("-f", "--from"),  type = "character", default = "NULL", 
        help = "brick start date e.g. 2015-08-01", metavar = "character"),
    make_option(c("-t", "--to"),    type = "character", default = "NULL", 
        help = "brick end date e.g. 2016-08-01",   metavar = "character")
)
opt_parser <- OptionParser(option_list = option_list);
opt <- parse_args(opt_parser)

if (any(is.null(opt$scene), is.null(opt$from), is.null(opt$to))) {
    print_help(opt_parser)
    stop("Invalid input!", call.=FALSE)
}

brick_scene <- opt$scene
brick_from  <- opt$from 
brick_to    <- opt$to

# configuration
#brick_scene  <- "225063"
#brick_scene  <- "226064"
#brick_scene  <- "233067"
#brick_from   <- "2016-08-01"
#brick_to     <- "2017-07-30"
#brick_from   <- "2015-08-01"
#brick_to     <- "2016-07-30"
#brick_from   <- "2014-08-01"
#brick_to     <- "2015-07-30"
#brick_from   <- "2013-08-01"
#brick_to     <- "2014-07-30"
brick_bands  <- c("sr_band2", "sr_band3", "sr_band4", "sr_band5", "sr_band6", "sr_band7") # include the bands required by the mixture model
brick_path   <- "/home/alber/shared/brick_maskcloud"
brick_n_img  <- 4
brick_prefix <- "LC8SR-MASKCLOUD"
landsat_path <- "/home/alber/landsat8"
modis_path   <- "/home/alber/MOD13Q1"
scene_shp    <- "/home/alber/Documents/data/experiments/l8mod-fusion/data/shp/wrs2_descending.shp"
tile_shp     <- "/home/alber/Documents/data/experiments/l8mod-fusion/data/shp/modis-tiles.shp"
temp_dir     <- "/home/alber/shared/tmp_maskcloud"
no_data      <- -9999
gdal_format  <- "GTiff"
gdal_options  <- "BIGTIFF=YES"

# handle temporal directories
tmp_dir <- temp_dir %>% file.path(paste("L8MOD", brick_scene, lubridate::year(brick_to), sep = "_"))
if (!dir.exists(tmp_dir)) dir.create(tmp_dir)

build_brick_maskcloud(landsat_path = landsat_path,
                      modis_path = modis_path,
                      scene_shp = scene_shp,
                      tile_shp = tile_shp,
                      brick_scene = brick_scene,
                      brick_from = brick_from,
                      brick_to = brick_to,
                      brick_bands = brick_bands,
                      brick_prefix = brick_prefix
                      brick_path = brick_path,
                      brick_n_img = brick_n_img,
                      gdal_options = gdal_options,
                      gdal_format = gdal_format,
                      no_data = no_data,
                      tmp_dir = tmp_dir) %>%
    print()