1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
import numpy as np, os
import crl, crl.data
from matplotlib import pyplot as plt
from shared import *
rea2008 = [
# CCT Ra GAI FSCI
["CW1", 5137, 75, 65, 60],
["CW2", 6682, 78, 97, 74],
["CW3", 6126, 71, 81, 64],
["CW4", 5854, 94, 99, 74],
["WW1", 3174, 95, 55, 67],
["WW2", 3443, 80, 64, 68],
["WW3b", 3732, 81, 74, 72],
["WW4", 3261, 91, 93, 77]
]
with TexWriter(os.path.join(os.path.dirname(__file__), "../build/rea2008.tex")) as fd:
for name, CCT_ref, Ra_ref, GAI_ref, FSCI_ref in rea2008:
path = os.path.join(os.path.dirname(__file__), "../lab/Rea2008/", name + ".csv")
data = np.sort(crl.data.load(path))
data = data[np.lexsort((data[:, 0], data[:, 1]))]
spec = crl.Spectrum(data[:, 1], data[:, 0])
res = crl.cri_Ra(spec)
GAI, _ = crl.cri_GAI(spec)
FSCI = crl.cri_FSCI(spec)
fd.write("%s & %g & %.4g & %g & %.1f & %g & %.1f & %g & %.1f \\\\\n"
% (name, CCT_ref, res.CCT, Ra_ref, res.CRI,
GAI_ref, GAI, FSCI_ref, FSCI))
|