data-raw/Rosco/read.rosco.files.r

In aphalo/photobiologyFilters: Spectral Transmittance and Spectral Reflectance Data

# read all available filter data
# assumes that all files with ".txt" are filter-data files
# assumes that all files with exactly four columns have two replicate measurements that need
# to be averaged
# STEPS
# 1) clear the workspace
# 2) read .txt files
# 3) save all the data.frames created into individual Rda file in data folder
library(photobiology)
library(dplyr)
rm(list = ls())
setwd("data-raw/Rosco")
file.list <- list.files(pattern = '*.txt$')
rosco.lst <- list()
for (file.name in file.list) {
  name <- sub(pattern=".txt", replacement="", x=file.name)
  name <- sub("Supergel", "Supergel_no", name)
  name <- sub("EColour", "EColour_no", name)
  used <- grepl("used", name)
  tmp.df <- read.table(file.name, header=TRUE)
  if (ncol(tmp.df) == 4) {
    tmp.df <- transmute(tmp.df, w.length = w.length.1,
                        Tfr = (transmittance.1 + transmittance.2) / 200)
  } else if (ncol(tmp.df) == 2) {
    tmp.df <- transmute(tmp.df, w.length = w.length, Tfr = transmittance / 100)
  } else {
    stop("ncol != 2 || ncol!=4")
  }
  tmp.df <- mutate(tmp.df, Tfr = ifelse(Tfr < 1e-5, 1e-5, Tfr))
  setFilterSpct(tmp.df, Tfr.type = "total")
  tmp.df <- setFilterProperties(tmp.df,
                                Rfr.constant = 0.08,
                                thickness = ifelse(grepl("Supergel", name),
                                                   90e-6, 75e-6),
                                attenuation.mode = "absorption") %>%
    #    clip_wl(c(240, NA)) %>%
    smooth_spct(method = "supsmu", strength = 0.5)
  setWhatMeasured(tmp.df, paste("Theatrical 'gel', type '",
                                gsub("_", " ", sub("_used", "", name)),
                                ifelse(used, "; used", "; new"),
                                "; from Rosco, UK",
                                sep = ""))
  setFilterProperties(tmp.df) # Missing metadata
  rosco.lst[[name]] <- tmp.df
}
rosco.mspct <- filter_mspct(rosco.lst)
names(rosco.mspct) <- paste("Rosco", names(rosco.mspct), sep = "_")
setwd("../..")

save(rosco.mspct, file = "data-raw/rda/rosco.mspct.rda")