R/efit2file.R

In TIMP: Fitting Separable Nonlinear Models in Spectroscopy and Microscopy

efit2file <-
  function(filename, skip = 2, numcol, nrows = vector()) {
    xx <- as.matrix(scan(filename, skip = skip, na.strings = c(
      "NA", "NAN",
      "NaNQ"
    )))
    if (length(nrows) == 0) {
      nrows <- length(xx) / numcol
    }
    dim(xx) <- c(numcol, sum(nrows))
    xx <- t(xx)
    cnt <- 1
    for (i in 1:length(nrows)) {
      write.table(xx[cnt:(cnt + nrows[i] - 1), ],
        row.names = FALSE,
        col.names = FALSE, file = paste("plainmat", i, filename,
          sep = "_"
        )
      )
      cnt <- cnt + nrows[i]
      cat("Wrote the file", paste("plainmat", i, filename,
        sep = "_"
      ), "\n")
    }
  }

# either copy this file into R or use
# source("efit2file.R")
# in R, when the file is in your working directory
#
# call the function with for example:
#
# efit2file("SADSrc685nm5S22nm1010.efit", numcol=8, nrows=c(256,256,256))
#
# where the file "SADSrc685nm5S22nm1010.efit" contains 3 sets of spectra,
# with each set having 256 rows and 8 columns.
#
# if nrows is not given, then just 1 matrix will be written to file,
# that includes everything in the efit file, for example with:
# efit2file("SADSrc685nm5S22nm1010.efit", numcol=8)