R/soundscapes_second.R

In audioblast/audioBlastAnalyse: Analysis service for audioBlast!

#' Analyses soundscapes in one second chunks
#'
#' Analyses performed (three second chunks):
#' - TDSC
#'
#' @param db database connector
#' @param source Source
#' @param id id (unique within source)
#' @param file url to file
#' @param type MIME type
#' @param duration audio duration
#' @param tmp Location to download temp file
#' @param force If TRUE recalculates all values
#' @param verbose If TRUE outputs debugging information
#' @export
#' @importFrom sonicscrewdriver readAudio rainfallDetection allChannels channels_se
#' @importFrom rjson toJSON
#' @importFrom tdsc tdsc

soundscapes_by_second <- function(db, source, id, file, type, duration, tmp, force=FALSE, verbose=FALSE) {
  duration <- av_media_info(tmp)$duration
  n <- ceiling(duration/3)
  if (force==TRUE) {
    deleteAllAnalyses(db, source, id, justR=TRUE)
  }

  if (duration == 0) {
    if (verbose) print("Provided duration is zero -- skipping")
    return()
  }

  for (i in (1:n)) {
    from <- (i-1)*3
    to <- (i)*3
    complete <- 1

    if (duration < 3 || (duration - (i-1)*3) < 3) {
      complete <- 0
    }

    if (complete == 0 || i*3 == duration) {
      w <- readAudio(tmp, from=from, units="seconds")
    } else {
      w <- readAudio(tmp, from=from, to=to, units="seconds")
    }

    if (is.logical(w)) return()

    #If only one value can't calculate TDSC (a failure mode for recordings)
    if (length(unique(w@left)) == 1) return()

    if (verbose) print(paste("tdsc startTime:", from))
    v <- allChannels(w, tdsc, max_D=14, channel.param=NULL, output.FUN = channels_tdsc)
    insertAnalysis(db, "analysis_3sec-tdsc", source, id, from, v, complete)

  }
  sql = paste0("UPDATE `recordings-calculated` SET `soundscapes_second` = 1 ",
               "WHERE `source` = ", dbQuoteString(db, source),
               " AND `id` = ", dbQuoteString(db, id), ";")
  dbeq <- abdbExecute(db, sql)
  if (verbose) print(dbeq)

}

channels_tdsc <- function(...) {
  params = list(...)
  return(toJSON(params[[1]]@a_matrix))
}