R/import.R

In mendivilg/eisr: eisR - Fitting EIS

openDTA <- function(arch) {

    conn <- file(arch, open = "r")
    linn <- readLines(conn)
    for (i in 1:length(linn)) {
        if (linn[i] == "ZCURVE\tTABLE")
            comienzo <- i + 2
    }
    close(conn)
    tmp <- read.delim2(arch, encoding = "UTF8", header = FALSE, skip = comienzo)
    tmp <- tmp[, c(3, 4, 5, 6, 10,  11)]
    colnames(tmp) <- c("time", "freq", "real", "imag", "Idc", "Vdc") # OJO el
    tmp
}

openZb <- function(arch) {

    conn <- file(arch, open = "r")
    linn <- readLines(conn)
    for (i in 1:length(linn)) {
        if (linn[i] == "\"Frequency\"")
            comienzo <- i + 3
    }
    close(conn)
    tmp <- read.delim2(arch, encoding = "UTF8", header = FALSE, sep = ",", skip = comienzo, stringsAsFactors = FALSE)
    tmp <- tmp[, c(1, 5, 6)]
    colnames(tmp) <- c("freq", "real", "imag")
    tmp$freq <- as.double(tmp$freq)
    tmp$real <- as.double(tmp$real)
    tmp$imag <- as.double(tmp$imag)
    tmp
}

openZ <- function(arch) {

    conn <- file(arch, open = "r")
    linn <- readLines(conn)
    for (i in 1:length(linn)) {
        if (linn[i] == "End Comments")
            comienzo <- i
    }
    close(conn)
    tmp <- read.delim2(arch, encoding = "UTF8", header = FALSE, skip = comienzo, stringsAsFactors = FALSE)
    tmp <- tmp[, c(1, 5, 6)]
    colnames(tmp) <- c("freq", "real", "imag")
    tmp$freq <- as.double(tmp$freq)
    tmp$real <- as.double(tmp$real)
    tmp$imag <- as.double(tmp$imag)
    tmp
}

openNOVA <- function(arch) {

    conn <- file(arch, open = "r")
    tmp <- read.delim2(arch, encoding = "UTF8", header = FALSE, sep = ";", skip = 1)
    tmp <- tmp[, c(2, 3, 4)]
    colnames(tmp) <- c("freq", "real", "imag")
    tmp$imag <- tmp$imag * -1
    tmp
}


#'\code{openEIS()}
#'
#' Importa datos de impedancia de distintos softwares adquisidores
#'
#'
#' Esta función extrae la frecuencia , parte imaginaria y parte real de los archivos generados por:
#'
#'\itemize{
#'   \item *.DTA file,  Gamry Framework 7.06 version
#'   \item *.z file,  ZPlot software 2.6 and 2.6b version
#'   \item *.txt ascii AUTOLAB NOVA 2.1.3 version
#'}
#' @return dataframe con variables freq, real, imag
#' @export
#'
#' @examples
#' mydata <- openEIS()
openEIS <- function() {
    arch <- path.expand(file.choose())
    conn <- file(arch, open = "r")
    linn <- readLines(conn)
    for (i in 1:length(linn)) {
        # .DTA file, Gamry Framework 7.06 version
        if (linn[i] == "EXPLAIN") {
            close(conn)
            df <- openDTA(arch)
            break
        }
        # .z file, ZPlot software 2.6 version
        if (linn[i] == "ZPLOT2 ASCII") {
            close(conn)
            df <- openZ(arch)
            break
        }
        # .z file, ZPlot software 2.6b version
        if (linn[i] == "\"ZPlotW Data File: Version 2.6b\"") {
            close(conn)
            df <- openZb(arch)
            break
        }
        # .txt ascii NOVA 2.1.3 version

        if (linn[i] == "Index;Frequency (Hz);Z' (Ω);-Z'' (Ω);Z (Ω);-Phase (°);Time (s)") {
            close(conn)
            df <- openNOVA(arch)
            break
        }

    }  #for
    df
}