R/getSeries_from_node.R

In moritzlindner/PatchR: An environment for working with ephys (PatchClamp) recordings

#' getSeries_from_node
#'
#' function inherited from ephys2 (/github/tdanker/ephys2/). Undocumented.
#' @importFrom stringr str_replace
#' @noRd
getSeries_from_node<-function (node, sweeps = 1:length(node), trace = 1, read_data = T,encoding=getOption("encoding"))
{
  path <- attr(node, "path")
  filename <- path[1]
  if (!file.exists(filename)) {
    filename <- paste("../", filename, sep = "")
  }
  con <- file(filename, "rb", encoding=encoding)
  signature <- readChar(con, 8)
  version <- readChar(con, 32)
  time <- readBin(con, "double")
  time <- time- 1580970496.0 # JanFirst1990
  if (time < 0.0){time<-time+4294967296} #HIGH_DWORD
  time <- time+ 9561652096 #MAC_BASE
  time<-as.Date(as.POSIXct(time, origin="1601-01-01"))
  if (signature == "DAT1") {
    mypul <- str_replace(filename, ".dat$",
                                  ".pul")
    if (!file.exists(mypul))
      stop("while trying to read from *.dat file, detected the old 'DAT1' format (no bundle file), therefore expected but could not find a correspoding *.pul file")
    con_pul <- file(mypul, "rb")
  }
  else {
    con_pul <- con
  }
  series <- lapply(node[sweeps], function(sweep) {
    getTrace_(con_pul, sweep[[trace]], read_data = read_data,
              name = names(sweep)[[trace]], con_dat = con)
  })

  SwTimer <- unlist(lapply(node[sweeps], function(sweep) {
    readAny(attr(sweep, "dataptr"), con_pul, 56, "double",
            8)
  }))

  close(con)
  if (signature == "DAT1") {
    close(con_pul)
  }
  maxLength <- max(unlist(lapply(series, length)))
  series <- lapply(series, function(s) {
    length(s) <- maxLength
    s
  })
  series_Y <- do.call(cbind, series)
  Xinterval <- attr(series[[1]], "Xinterval")
  nDatapoints_ <- attr(series[[1]], "nDatapoints_")
  RecMode <- attr(series[[1]], "RecMode")
  RPip <- attr(series[[1]], "RPip")
  Urest <- attr(series[[1]], "Urest")
  RSeal <- attr(series[[1]], "RSeal")
  Cs <- attr(series[[1]], "Cs")
  Rs <- attr(series[[1]], "Rs")

  XUnit <- attr(series[[1]], "XUnit")
  YUnit <- attr(series[[1]], "YUnit")
  
  DataScaler <- attr(series[[1]], "DataScaler")
  Unit_from_file <- attr(series[[1]], "Unit_from_file")   
  TrYOffset <- attr(series[[1]], "TrYOffset")   
  TrYRange <- attr(series[[1]], "TrYRange")   

  series_X1 <- 1:nDatapoints_ * Xinterval
  series_X <- matrix(rep(series_X1, length(series_Y[1, ])),
                     nDatapoints_, length(series_Y[1, ]))
  series <- list(y = series_Y,
                 x = series_X,
                 filename = filename,
                 exp = path[2],
                 ser = path[3],
                 trace = trace,
                 RecMode = RecMode,
                 RPip = RPip,
                 Urest = Urest,
                 RSeal = RSeal,
                 Cs = Cs,
                 Rs = Rs,
                 time = time,
                 version = version,
                 XUnit = XUnit,
                 YUnit = YUnit,
                 DataScaler = DataScaler,
                 Unit_from_file = Unit_from_file,
                 TrYOffset = TrYOffset,
                 TrYRange = TrYRange,
                 tracename = attr(series[[1]],
                                  "name"),
                 Trace_Time = SwTimer,
                 Serieslabel = getStimName_from_unique_seriesName(path[3]),
                 Stimulus = attr(node, "StimulusName"))
  class(series) <- "HEKAseries"
  series
}