data-raw-old/save-irrad-data-2019.r

In aphalo/photobiologyLEDs: Spectral Data for Light-Emitting-Diodes

library(photobiology)
library(dplyr)

rm(list = ls(pattern = "*"))

path <- "./data-raw/Rda/"

# load collections of spectra

files2read <- list.files(path = path, pattern = "*[-.]mspct.rda$")
files2read <- setdiff(files2read, "leds-mspct.rda")
for (f in files2read) {
  load(paste(path, f, sep = ""))
}

# concatenate collections

collections2bind <- ls(pattern = "*.mspct$")
leds.mspct <- source_mspct()
for (mspct in mget(collections2bind)) {
  leds.mspct <- c(leds.mspct, mspct)
}
rm(list = c(collections2bind, "mspct"))

names(leds.mspct)

# read individual spectra

files2read <- list.files(path = path, pattern = "*[-.]spct.rda$")
# files2read <- list.files(path = path, pattern = "*.spct.rda$")
for (f in files2read) {
  load(paste(path, f, sep = ""))
}

# add spectra to the collection

spectra2bind <- ls(pattern = "*[.]spct$")
spectra.ls <- mget(spectra2bind)
names(spectra.ls) <- gsub(".spct$", "", names(spectra.ls))
leds.mspct <- c(leds.mspct, source_mspct(spectra.ls))

leds.mspct <- leds.mspct[order(names(leds.mspct))]

rm(list = spectra2bind)

names(leds.mspct)

# normalize all spectra to their maxima

leds.mspct <- normalize(leds.mspct)

# group into a single spectrum the three channels

NHXRGB090.spct <- rbindspct(list(R = leds.mspct$NHXRGB090_R,
                                 G = leds.mspct$NHXRGB090_G,
                                 B = leds.mspct$NHXRGB090_B),
                            idfactor = "channel")

leds.mspct[["NHXRGB090"]] <- NHXRGB090.spct

# save collection
leds.mspct <- leds.mspct[-which(names(leds.mspct) %in% c("spct_1", "UVA", "FR_OLD"))]
names(leds.mspct)[names(leds.mspct) == "BLUE"] <- "HPR40E_48K30BI"

save(leds.mspct, file = "leds-mspct.rda")

# metadata

how_measured.ls <- list()
for (s in names(leds.mspct)) {
  how_measured.ls[[s]] <- getHowMeasured(leds.mspct[[s]])
}

normalized.ls <- list()
for (s in names(leds.mspct)) {
  normalized.ls[[s]] <- getNormalised(leds.mspct[[s]])
}

led.whats <- list(#"BS436" = "Blue LED Roithner-Laser B5-436, 5mm",
                  "LY5436" = "Amber LED Osram",
                  "CB30" = "Blue LED Roithner-Laser CB30",
                  "LED405" = "Violet LED (InGaN) Roithner-Laser LED405",
                  "LED740" = "NIR LED Roithner-Laser LED740",
                  "UV395" = "UVA LED Roithner-Laser UV395",
                  "XSL365" = "UVA LED Roithner-Laser XSL-365-TF, TO46",
                  "XSL370" = "UVA LED Roithner-Laser XSL-370-TF, TO46",
                  "XSL375" = "UVA LED Roithner-Laser XSL-375-TF, TO46",
                  "UVMAX340" = "UVA LED 340nm Roithner-Laser (SeTi UVCLEAN) type UVMAX340-HL-15 (TO3 hemispherical lens), measured at 100 mm distance from cosine difuser",
                  "UVMAX305" = "UVA LED 310nm Roithner-Laser (SeTi UVCLEAN) type UVMAX305-HL-15 (TO3 hemispherical lens), measured at 100 mm distance from cosine difuser",
                  "TY_UV310nm" = "UVB LED TaoYuan 310nm",
                  "white" = "White LED from hardware store (Clas Ohlsson, Finland)",
                  "Q36_4000K" = "White LED Lumitronix Q36 4000K (neutralwhite) array, SKU 53681, 36 Nichia 757 LEDs, 39W",
                  "QDDH66002" = "Red LED 660nm Quantum Devices, 5mm epoxi package, 0.1W",
                  "QDDH68002" = "Red LED 680nm Quantum Devices, 5mm epoxi package, 0.1W",
                  "QDDH70002" = "Far-red LED 700nm Quantum Devices, 5mm epoxi package, 0.1W",
                  "QDDH73002" = "Far-red LED 7300nm Quantum Devices, 5mm epoxi package, 0.1W"  #,
                  # "G_P30R140A1_XT" = "",
                  # "weili430nm" = "",
                  # "weili540nm" = "",
                  # "weili665nm" = "",
                  # "weili740nm" = "",
                  # "HPR40E_48K30BG" = "",
                  # "NHXRGB090_R" = "",
                  # "NHXRGB090_G" = "",
                  # "NHXRGB090_B" = ""
                  )

what_measured(leds.mspct[["HPR40E_48K30BI"]]) <- unname(getWhatMeasured(leds.mspct[["HPR40E_48K30BI"]]))

what_measured.ls <- list()
for (s in names(leds.mspct)) {
  what_measured.ls[[s]] <- getWhatMeasured(leds.mspct[[s]])
}
na_what <- names(what_measured.ls)[is.na(what_measured.ls)]
short_what <- names(what_measured.ls)[grepl("^LED, type  ", what_measured.ls)]
replace_what <- intersect(names(leds.mspct), names(led.whats))

for (s in replace_what) {
  what_measured(leds.mspct[[s]]) <- led.whats[[s]]
}

leds.mspct <- leds.mspct[order(names(leds.mspct))]

# Create category index vectors

roithner_laser <- 
  c("BS436", "LED405", "LED740", "UV395", "LED435_66_60",
    "LED740_01AV", "B5_436_30D", "UVMAX340", "UVMAX305", 
    "XSL365", "XSL370", "XSL375")
seti <- c("UVMAX340", "UVMAX305", "XSL365", "XSL370", "XSL375")
tao_yuan <- "TY_UV310nm"
# unknown <- c("white", "FR_OLD")
lumitronix <- "Q36_4000K"
hewlett_packard <- agilent <- c("HLMB_CB30", "HLMP_CB31", "HLMP_CM30", 
                                "HLMP_CM31", "HLMP_DJ32", "HLMP_DL32")
quantum_devices <- c("QDDH66002", "QDDH68002", "QDDH70002", "QDDH73502")
osram <- "LY5436"
norlux <- c("NHXRGB090_R", "NHXRGB090_G", "NHXRGB090_B", "NHXRGB090")
shenzhen_weili <- leds_global <- 
  c("G_P30R140A1_XT", "weili430nm", "weili540nm", "weili665nm", "weili740nm")
huey_jann <- "HPR40E_48K30BG"
led_engin <- c("LZ1_10DB00", "LZ1_10UA00_00U8", "LZ1_10UA00_00U4", 
               "LZ1_10UV00", "LZ1_10R302")

uv_leds <- c("UVMAX340", "UVMAX305", "XSL365", "XSL370", "XSL375",
             "G_P30R140A1_XT", "LZ1_10UA00_00U4", "LZ1_10UV00")
red_leds <- c("QDDH66002", "QDDH68002", "QDDH70002", "QDDH73502", 
              "weili740nm", "weili665nm",
              "NHXRGB090_R", "LED740_01AV", "FR_OLD", "LZ1_10R302")
amber_leds <- c("LY5436")
green_leds <- c("NHXRGB090_G", "weili540nm")
blue_leds <- c("NHXRGB090_B", "LED435_66_60", "weili430nm", "LZ1_10UA00_00U8")
white_leds <- c("white", "Q36_4000K")
multichannel_leds <- c("NHXRGB090")

oo_maya_leds <- names(leds.mspct)

save(leds.mspct, roithner_laser, seti, tao_yuan, unknown, lumitronix,
     hewlett_packard, quantum_devices, osram, norlux, shenzhen_weili,
     leds_global, huey_jann, led_engin, oo_maya_leds,
     uv_leds, red_leds, amber_leds, green_leds, blue_leds,
     white_leds, multichannel_leds,
     file = "data/leds-mspct.rda")

tools::resaveRdaFiles("data", compress="auto")
print(tools::checkRdaFiles("data"))