data-raw/maya-rda/Weili/12-channels-arrays/make-RDS-array-B.R

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

library(photobiology)
library(ggspectra)

energy_as_default()

# clear workspace
rm(list = ls(pattern = "*"))

load("~/R-pkgs-owned/photobiologyLEDs/data-raw/maya-rda/Weili/12-channels-arrays/collection.Custom.array.2.with.12.channel.x.10W.Rda")

autoplot(collection.Custom.array.2.with.12.channel.x.10W.irrad.mspct, facets = 4)

# autoplot(s_sum(collection.Array.1.12.channels.irrad.mspct[-c(1,6,7)]), span = 21)

Weili_120W_12ch_B.mspct <- collection.Custom.array.2.with.12.channel.x.10W.irrad.mspct

old.names <- names(Weili_120W_12ch_B.mspct)

names(Weili_120W_12ch_B.mspct) <-
  paste("ch", gsub(" ", "0", format(1:12, width = 2)), sep = "")

# summary(Weili_12ch_1.mspct, which.metadata = "when.measured")

how.measured <- "Array spectrometer, Ocean Optics Maya 2000 Pro; Bentham cosine diffuser D7H; distance 70 mm; LED maximum current 350 mA per channel, adjusted to avoid clipping."

comment.text <- paste("Custom assembled 120W 12 channel LED array",
                      "in a 40x56mm metal based package with 10 x 1W LED dies per channel.\n",
                      "The LED dies are located within 24.6x24.6mm area and encased in clear silicone.\n",
                      "Array supplied by Shezhen Weili (https://www.leds-global.com/), assembled with dies from Epileds/Epistar/Bridgelux")
what.measured <- "Custom assembled 120W 12 channel LED array"

for (s in names(Weili_120W_12ch_B.mspct)) {
  temp.spct <- Weili_120W_12ch_B.mspct[[s]]
  temp.spct <- smooth_spct(temp.spct)
  temp.spct <- thin_wl(temp.spct)
  temp.spct <- trim_wl(temp.spct, range = c(280, 900), fill = 0)
  # if (s == "ch10") {
  #   idx <- which(abs(temp.spct[["w.length"]] - 407.69) < 0.1)
  #   print(idx)
  #   temp.spct[idx, "s.e.irrad"] <- 
  #     (temp.spct[idx - 1, "s.e.irrad"] + temp.spct[idx + 1, "s.e.irrad"]) / 2
  # }
  setHowMeasured(temp.spct, how.measured)
  setWhatMeasured(temp.spct, 
                  paste("Custom assembled 120W 12 channel LED array,", s))
  comment(temp.spct) <- NULL
  trimInstrDesc(temp.spct)
  trimInstrSettings(temp.spct)
#  print(str(get_attributes(temp.spct)))
  print(autoplot(temp.spct, annotations = c("+", "title:what:when:comment")))
  Weili_120W_12ch_B.mspct[[s]] <- temp.spct
  readline("next:")
}

autoplot(Weili_120W_12ch_B.mspct)
# Weili_120W_12ch_B.mspct <- normalise(Weili_120W_12ch_B.mspct)
# Weili_120W_12ch_B.mspct <- setNormalised(Weili_120W_12ch_B.mspct)
autoplot(normalise(Weili_120W_12ch_B.mspct))

Weili_120W_12ch_B.spct <- rbindspct(Weili_120W_12ch_B.mspct, idfactor = "channel")
how_measured(Weili_120W_12ch_B.spct) <- how.measured
what_measured(Weili_120W_12ch_B.spct) <- "Custom assembled 120W 12 channel LED array"
comment(Weili_120W_12ch_B.spct) <- comment.text

names(Weili_120W_12ch_B.spct)

saveRDS(Weili_120W_12ch_B.spct, "data-raw/maya-rda/Weili/12-channels-arrays/weili_12ch_B.rds")