R/mosaic_fun_sen2.R

In spatialstatisticsupna/rsat: Dealing with Multiplatform Satellite Images

mosaic_fun_sen2 <- function(mfiles) {
  allfun <- NULL
  allfun$readfromscratch <- function(m, ...) {
    file.path("/vsizip", m, utils::unzip(m, list = TRUE)$Name)
  }
  allfun$filterchunks <- function(allfiles, bnds) {
    allfiles[grepl(paste0(bnds, ".jp2"), allfiles, ignore.case = TRUE)]
  }
  bands <- allfun$readfromscratch(mfiles[1])
  bands <- bands[grepl(".jp2", bands, ignore.case = TRUE)]
  if (any(grepl("60m", bands))) {
    bands <- unlist(lapply(strsplit(gsub(".*/",
                                         "",
                                         gsub(".*\\\\", "", bands)),
                                    "_"),
                           function(x) {
      do.call(paste, as.list(c(x[c(length(x) - 1, length(x))], sep = "_")))
    }))
    bands <- gsub(".*\\d{8}T\\d{6}_", "", bands)
  } else {
    bands <- gsub(".*_", "", bands)
  }
  allfun$bands <- gsub(".jp2", "", bands)

  allfun$defineNodata <- function(chunks, bnds) {
    NULL
  }
  allfun
}