#' Read the ECO-VSF data
#'
#' @param eco.file is the file name of the ECO-VSF3 data
#' @param ECO.bands is a character string indicating the bands (BLUE, GREEN and RED),
#' i.e. either "B", "G", "R" or "BG", or "BGR".
#' (Default is "BG")
#' @return A dataframe with the raw data
#'
#' @author Simon Belanger
#' @export
#'
read.ECOVSF<- function(eco.file, ECO.bands="BG"){
if (ECO.bands == "BG") {
id <- file(eco.file, "r")
line<- unlist(strsplit(readLines(con=id, n =1), " ") )
if (line[1] == "Created")
{
print("Data recorded using ECOHost")
date<- as.POSIXct(line[3], format="%m-%d-%Y")
line<- unlist(strsplit(readLines(con=id, n =1), " ") )
line<- unlist(strsplit(readLines(con=id, n =1), " ") )
# Start counting the number of record
nrec=0
while (line[1] != "EOF") {
line<- unlist(strsplit(readLines(con=id, n =1), " ") )
nrec=nrec+1
}
nrec = nrec - 1
print(nrec)
close(id)
raw = read.table(file=eco.file, skip=3, nrows = nrec)
names(raw) <- c("SN", "time", "B100", "B125", "B150", "Bref",
"G100", "G125", "G150", "Gref", "Vin")
}
else
{
print("Data recorded using ECOview")
date<-NA
raw = read.table(file=eco.file)
names(raw) <- c("SN", "time", "B100", "B125", "B150", "Bref",
"G100", "G125", "G150", "Gref", "Vin")
nrec = length(raw$time)
}
return(list(date=date, time=raw$time, raw=raw,
nrec=nrec))
}
if (ECO.bands == "B") {
print("Data recorded using ECOview")
date<-NA
raw = read.table(file=eco.file, skip=5)
names(raw) <- c("SN", "time", "B100", "B125", "B150", "Bref")
nrec = length(raw$time)
return(list(date=date, time=raw$time, raw=raw,
nrec=nrec))
}
if (ECO.bands == "G") {
print("Data recorded using ECOview")
date<-NA
raw = read.table(file=eco.file, skip=5)
names(raw) <- c("SN", "time", "G100", "G125", "G150", "Bref")
nrec = length(raw$time)
return(list(date=date, time=raw$time, raw=raw,
nrec=nrec))
}
if (ECO.bands == "R") {
print("Data recorded using ECOview")
date<-NA
raw = read.table(file=eco.file, skip=5)
names(raw) <- c("SN", "time", "R100", "R125", "R150", "Bref")
nrec = length(raw$time)
return(list(date=date, time=raw$time, raw=raw,
nrec=nrec))
}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.