Publication Cover
LEUKOS
The Journal of the Illuminating Engineering Society
Volume 19, 2023 - Issue 2
2,126
Views
1
CrossRef citations to date
0
Altmetric
Research Article

Reconstruction of Daylight Spectral Power Distribution Based on Correlated Color Temperature: A Comparative Study between the CIE Approach and Localized Procedures in Assessing Non-image Forming Effects

ORCID Icon, ORCID Icon & ORCID Icon
Pages 118-145 | Received 26 May 2021, Accepted 01 Oct 2021, Published online: 15 Feb 2022

Figures & data

Fig. 1. Schematic representation of the CIE (Citation2018b) reconstruction procedure approach to reconstructing daylight spectral power distribution (Sλ) as a function of either chromaticity coordinates (x,y) or correlated color temperature (Tcp).

Fig. 1. Schematic representation of the CIE (Citation2018b) reconstruction procedure approach to reconstructing daylight spectral power distribution (Sλ) as a function of either chromaticity coordinates (x,y) or correlated color temperature (Tcp).

Table 1. Compilation of the reconstruction methods of daylight

Table 2. Summary of the CIE and locally adjusted reconstruction approaches

Table 3. Percentage of colorimetrically “accurate” (GFC0.9950), “good” (GFC0.9990), and “exact” (GFC0.9999) reconstructions depending on the reconstruction method

Table 4. MeanGFC values for three wavelength ranges and the complete wavelength range depending on the selected reconstruction method, color coded as follows: GFC<0.9950, GFC0.9950, GFC0.9990, GFC0.9999

Table 5. Mean (GFC) for sky conditions

Table 6. MeanGFC for summer and winter measurements

Fig. 2. Daylight measuring site at Technische Universität Berlin.

Fig. 2. Daylight measuring site at Technische Universität Berlin.

Fig. 3. Chromaticities of daylight in Berlin in the CIE 1931 x,y chromaticity diagram compared to the CIE daylight locus. The diagram on the left shows the modeling data set from the 597,998 SPD measurements gathered in 2016 (left); the one on the right shows the validation data set from the 538,205 SPD measurements gathered in 2015 according to the procedure laid out by the CIE (Citation2018b), limiting the xD from 0.25 to 0.38 (right).

Fig. 3. Chromaticities of daylight in Berlin in the CIE 1931 x,y chromaticity diagram compared to the CIE daylight locus. The diagram on the left shows the modeling data set from the 597,998 SPD measurements gathered in 2016 (left); the one on the right shows the validation data set from the 538,205 SPD measurements gathered in 2015 according to the procedure laid out by the CIE (Citation2018b), limiting the xD from 0.25 to 0.38 (right).

Fig. 4. Sky patches as defined in Tregenza (Citation1987, Citation2004) and the CIE (Citation1994). Patches shaded in gray were used for the accuracy analysis and sensitivity study; the patches used to examine atmospheric influences are circled in black.

Fig. 4. Sky patches as defined in Tregenza (Citation1987, Citation2004) and the CIE (Citation1994). Patches shaded in gray were used for the accuracy analysis and sensitivity study; the patches used to examine atmospheric influences are circled in black.

Fig. 5. Diagram of the approach to deriving the Berlin (BLN) daylight chromaticities xD,BLN,yD,BLN.

Fig. 5. Diagram of the approach to deriving the Berlin (BLN) daylight chromaticities xD,BLN,yD,BLN.

Fig. 6. Non-linear relation between chromaticities of daylight in Berlin (left) and the correlation between the daylight chromaticity coordinate xD and the inverse correlated color temperature Tcp (MK−1) (right).

Fig. 6. Non-linear relation between chromaticities of daylight in Berlin (left) and the correlation between the daylight chromaticity coordinate xD and the inverse correlated color temperature Tcp (MK−1) (right).

Fig. 7. Loci of the chromaticities of daylight obtained for different sites along with the CIE daylight locus and the Berlin daylight locus: Granada, Spain (Hernández-Andrés et al. Citation2001a, Citation2001b) (upper left); Vaulx-en-Velin, France (Chain Citation2004) (upper right); Bratislava, Slovakia (Rusnák Citation2014) (middle left); Berlin, Germany (Diakite et al. Citation2018) (middle right); Coburg and Bendigo, Australia (Dixon Citation1978) (lower left); Atsugi and Nagaoka, Japan (Kobayashi et al. Citation1997, Citation1996) (lower right).

Fig. 7. Loci of the chromaticities of daylight obtained for different sites along with the CIE daylight locus and the Berlin daylight locus: Granada, Spain (Hernández-Andrés et al. Citation2001a, Citation2001b) (upper left); Vaulx-en-Velin, France (Chain Citation2004) (upper right); Bratislava, Slovakia (Rusnák Citation2014) (middle left); Berlin, Germany (Diakite et al. Citation2018) (middle right); Coburg and Bendigo, Australia (Dixon Citation1978) (lower left); Atsugi and Nagaoka, Japan (Kobayashi et al. Citation1997, Citation1996) (lower right).

Fig. 8. Diagram of the procedure of deriving the Berlin components and Berlin weighting factors for SPD reconstruction.

Fig. 8. Diagram of the procedure of deriving the Berlin components and Berlin weighting factors for SPD reconstruction.

Fig. 9. Components of the daylight distribution S0(λ), S1(λ), S2(λ), CIE (left); Berlin (right).

Fig. 9. Components of the daylight distribution S0(λ), S1(λ), S2(λ), CIE (left); Berlin (right).

Fig. 10. Diagram of the reconstruction procedure.

Fig. 10. Diagram of the reconstruction procedure.

Fig. 11. Diagram of the approach to determine the Goodness-of-Fit Coefficient (GFC) of the CIE reconstruction procedure and locally adjusted procedures (BLN1 and BLN2).

Fig. 11. Diagram of the approach to determine the Goodness-of-Fit Coefficient (GFC) of the CIE reconstruction procedure and locally adjusted procedures (BLN1 and BLN2).

Table 7. Mean absolute percentage error MAPEα [%] for α-opic quantities

Table 8. Root mean square percentage error RMSPEα [%] for α-opic quantities

Table 9. Mean percentage error MPEα [%] for the α-opic quantities

Fig. 12. Diagram of the approaches for the sensitivity study.

Fig. 12. Diagram of the approaches for the sensitivity study.

Fig. 13. Normalized measured spectral irradiance (-, gray, SMeλ) with the best achieved maxGFCBLN2=0.99998 conducted on 10 September 2015, at 14:58:13 CEST, with a solar altitude of 37.17°, for sky patch 114, with a Tcp equal to 6016 K (overcast sky) compared to the corresponding reconstructions performed using: the CIE method (-, SCIE,Reλ); the BLN1 method, using xD,BLN,yD,BLN(--,SBLN1,Reλ); and the BLN2 method, using xD,BLN,yD,BLN and components of the daylight distribution S0,BLN, S1,BLN, and S2,BLN (:, SBLN2,Reλ).

Fig. 13. Normalized measured spectral irradiance (-, gray, SMeλ) with the best achieved maxGFCBLN2=0.99998 conducted on 10 September 2015, at 14:58:13 CEST, with a solar altitude of 37.17°, for sky patch 114, with a Tcp equal to 6016 K (overcast sky) compared to the corresponding reconstructions performed using: the CIE method (-, SCIE,Reλ); the BLN1 method, using xD,BLN,yD,BLN(--,SBLN1,Reλ); and the BLN2 method, using xD,BLN,yD,BLN and components of the daylight distribution S0,BLN, S1,BLN, and S2,BLN (:, SBLN2,Reλ).

Fig. 14. Goodness-of-fit of the reconstructions as expressed by GFC0.9950, as a function of the inverse correlated color temperature Tcp. Upper left, calculated using the CIE method; upper right, the BLN1 method, using xD,BLN,yD,BLN; and lower left, the BLN2 method, using xD,BLN,yD,BLN and the components S0,BLN, S1,BLN, and S2,BLN. The numerical data of the percentages of colorimetrically accurate reconstructions can be found in the supplemental material.

Fig. 14. Goodness-of-fit of the reconstructions as expressed by GFC≥0.9950, as a function of the inverse correlated color temperature Tcp. Upper left, calculated using the CIE method; upper right, the BLN1 method, using xD,BLN,yD,BLN; and lower left, the BLN2 method, using xD,BLN,yD,BLN and the components S0,BLN, S1,BLN, and S2,BLN. The numerical data of the percentages of colorimetrically accurate reconstructions can be found in the supplemental material.

Fig. 15. Percentage of colorimetrically “accurate” (GFC0.9950) reconstructions according to sky conditions for each of the three methods. Upper left, the CIE method; upper right, the BLN1 method (xD,BLN,yD,BLN); and lower left, the BLN2 method (xD,BLN,yD,BLN and components S0,BLN, S1,BLN, and S2,BLN). The values above the bars indicate the percentage of colorimetrically “accurate” (GFC0.9950) reconstructions.

Fig. 15. Percentage of colorimetrically “accurate” (GFC≥0.9950) reconstructions according to sky conditions for each of the three methods. Upper left, the CIE method; upper right, the BLN1 method (xD,BLN,yD,BLN); and lower left, the BLN2 method (xD,BLN,yD,BLN and components S0,BLN, S1,BLN, and S2,BLN). The values above the bars indicate the percentage of colorimetrically “accurate” (GFC≥0.9950) reconstructions.

Fig. 16. Percentage of colorimetrically “accurate” (GFC0.9950) reconstructions according to seven sky patches (north and zenith direction, nos. 45, 73, 96, 118, 132, 142, 145). Upper left, the CIE method; upper right, the BLN1 method (using xD,BLN,yD,BLN); and lower left, the BLN2 method (using xD,BLN,yD,BLN and components S0,BLN, S1,BLN, and S2,BLN). The values above the bars indicate the percentage of colorimetrically “accurate” (GFC0.9950) reconstructions.

Fig. 16. Percentage of colorimetrically “accurate” (GFC≥0.9950) reconstructions according to seven sky patches (north and zenith direction, nos. 45, 73, 96, 118, 132, 142, 145). Upper left, the CIE method; upper right, the BLN1 method (using xD,BLN,yD,BLN); and lower left, the BLN2 method (using xD,BLN,yD,BLN and components S0,BLN, S1,BLN, and S2,BLN). The values above the bars indicate the percentage of colorimetrically “accurate” (GFC≥0.9950) reconstructions.

Fig. 17. An overall (left) and detailed (right) display of a normalized measured spectral irradiance (-, gray, SMeλ) retrieved on 30 March 2015, at 17:29:01 CEST, with a solar altitude of 18.09°, for sky patch 75 with Tcp=76,294 K as well as the corresponding reconstructions performed using the CIE method (-, SCIE,Reλ); the BLN1 method, using xD,BLN,yD,BLN (--,SBLN1,Reλ); and the BLN2 method, using xD,BLN,yD,BLN and components of the daylight distribution S0,BLN, S1,BLN, and S2,BLN (:, SBLN2,Reλ).

Fig. 17. An overall (left) and detailed (right) display of a normalized measured spectral irradiance (-, gray, SMeλ) retrieved on 30 March 2015, at 17:29:01 CEST, with a solar altitude of 18.09°, for sky patch 75 with Tcp=76,294 K as well as the corresponding reconstructions performed using the CIE method (-, SCIE,Reλ); the BLN1 method, using xD,BLN,yD,BLN (--,SBLN1,Reλ); and the BLN2 method, using xD,BLN,yD,BLN and components of the daylight distribution S0,BLN, S1,BLN, and S2,BLN (:, SBLN2,Reλ).
Supplemental material

Supplemental Material

Download MS Word (329 KB)