R/mosaic_fun_mod.R

In spatialstatisticsupna/rsat: Dealing with Multiplatform Satellite Images

mosaic_fun_mod <- function(mfiles) {
  allfun <- NULL
  allfun$readfromscratch <- function(m, ...) {
    image.data <- unlist(strsplit(gdal_utils(
      util = "info",
      source = m, quiet = TRUE
    ), "\n"),
    recursive = FALSE
    )
    bands.names <- image.data[grepl(".*SUBDATASET_.*_NAME=", image.data)]
    bands.names <- gsub(".*SUBDATASET_.*_NAME=", "", bands.names)
    bands.names
  }
  allfun$filterchunks <- function(allfiles, bnds) {
    allfiles[grepl(pattern = bnds, x = allfiles, perl = TRUE)]
  }
  bands <- allfun$readfromscratch(mfiles[1])
  bands <- gsub(".*SUBDATASET_.*_NAME=", "", bands)
  allfun$bands <- gsub('.*":', "", bands)
  allfun$defineNodata <- function(chunks, bnds) {
    NULL
  }

  allfun
}