Next Article in Journal
Relationship between EMG and fNIRS during Dynamic Movements
Next Article in Special Issue
Spectral Filter Selection Based on Human Color Vision for Spectral Reflectance Recovery
Previous Article in Journal
Field-Deployable Determinations of Peroxide Index and Total Phenolic Content in Olive Oil Using a Promising Portable Sensor System
Previous Article in Special Issue
ABANICCO: A New Color Space for Multi-Label Pixel Classification and Color Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination and Measurement of Melanopic Equivalent Daylight (D65) Illuminance (mEDI) in the Context of Smart and Integrative Lighting

1
Laboratory of Adaptive Lighting Systems and Visual Processing, Technical University of Darmstadt, Hochschulstr. 4a, 64289 Darmstadt, Germany
2
ERCO GmbH, Brockhauser Weg 80-82, 58507 Lüdenscheid, Germany
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(11), 5000; https://doi.org/10.3390/s23115000
Submission received: 5 April 2023 / Revised: 7 May 2023 / Accepted: 17 May 2023 / Published: 23 May 2023
(This article belongs to the Special Issue Recent Trends and Advances in Color and Spectral Sensors)

Abstract

:
In the context of intelligent and integrative lighting, in addition to the need for color quality and brightness, the non-visual effect is essential. This refers to the retinal ganglion cells ( i p R G C s ) and their function, which were first proposed in 1927. The melanopsin action spectrum has been published in CIE S 026/E: 2018 with the corresponding melanopic equivalent daylight (D65) illuminance ( m E D I ), melanopic daylight (D65) efficacy ratio ( m D E R ), and four other parameters. Due to the importance of m E D I and m D E R , this work synthesizes a simple computational model of m D E R as the main research objective, based on a database of 4214 practical spectral power distributions (SPDs) of daylight, conventional, LED, and mixed light sources. In addition to the high correlation coefficient R 2 of 0.96795 and the 97% confidence offset of 0.0067802, the feasibility of the m D E R model in intelligent and integrated lighting applications has been extensively tested and validated. The uncertainty between the m E D I calculated directly from the spectra and that obtained by processing the RGB sensor and applying the m D E R model reached ± 3.3% after matrix transformation and illuminance processing combined with the successful m D E R calculation model. This result opens the potential for low-cost RGB sensors for applications in intelligent and integrative lighting systems to optimize and compensate for the non-visual effective parameter m E D I using daylight and artificial light in indoor spaces. The goal of the research on R G B sensors and the corresponding processing method are also presented and their feasibility is methodically demonstrated. A comprehensive investigation with a huge amount of color sensor sensitivities is necessary in a future work of other research.

1. Introduction

The development of lighting technology from the beginning of the 20th century to the present can be roughly divided into three stages, some of which overlap. Until the end of the 20th century, with thermal radiators and discharge lamps as light sources, the first component of integrative HCL (“Human Centric Lighting”) lighting technology was about “visual performance” (e.g., reaction time, contrast perception, reading speed, visual acuity) to enable visibility, improve work performance, and minimize error rates in industrial, educational, and healthcare facilities and offices, among others [1,2,3,4,5,6,7,8,9]. Since the 1970s and even more so since the mid-1990s, with the growing importance of education and the information society, the psychological effects of light (e.g., scene preference, room preference, user satisfaction, spatial perception, attractiveness, and color preference) have become more and more important. As a result of this research process, international and national standards (e.g., EN 12464 standard [9]) have been formulated for lighting design and the development of interior lighting systems. For special applications (e.g., exhibitions, meeting rooms in prestigious buildings, and treatment rooms in medical centers and hospitals) the following values in Table 1 [10,11,12,13,14,15,16,17] from research studies have been used.
In [18], De Boer and Fischer stated in 1978 that, for a horizontal illuminance of 1000 lx (which roughly corresponds to the preferred illuminance according to the above-mentioned study by Moosmann [13]), a luminance of 200 cd/m2 is preferred for the ceiling and 100 cd/m2 for the walls. At a horizontal illuminance of 500 lx, the above values can be estimated to be about 210 cd/m2 (ceiling) and 70 cd/m2 (walls). Another specification of the preferred luminance at 75 cd/m2 on the wall at a horizontal illuminance of E v = 500 lx was found in Balder’s 1957 publication [19]. Since the beginning of the 21st century, with the quantitative discovery of intrinsically photosensitive retinal ganglion cells ( i p R G C s ), other non-visual effects of light have been considered, such as circadian rhythm, hormone production and suppression, sleep quality, alertness, and mood. The overall picture of light effects on humans in the three contexts of “visual performance”, “psychological light effects”, and “non-visual light effects” is shown in Figure 1, which attempts to describe a chain of signal processing from the optical and temporal input parameters to the light users with influencing parameters to the output parameters expressing the physiological, behavioral, and human biological effects of light and lighting.
In recent years, there has been a large amount of research and publications on non-visual lighting effects that have been analyzed and reviewed in books and literature reviews [20,21,22,23]. In the majority of scientific publications cited in these reviews, the input parameters describing lighting conditions were photometric metrics (illuminance, luminance) and/or color temperature. In addition to the optimal lighting design of buildings and the development of modern luminaires according to the criteria of visual performance and psychological lighting effects, it is necessary to define suitable input parameters for the description of non-visual lighting effects to determine their optimal values and to measure these parameters in a wide variety of applications in the laboratory and in the field. To achieve the above goal, the following three research questions for lighting research must be formulated in general terms:
  • Which input parameters can be used to describe the non-visual effects of light in their variety of manifestations (alertness, sleep quality, hormone production and suppression, and phase shift)?
  • What values of these input variables are currently considered in the literature to be minimal, maximal, or optimal?
  • Which measuring devices, sensor systems, and measuring methods can be used to measure the input quantities for non-visual lighting effects and process them in the context of smart lighting in the course of the control and regulation of LED or OLED luminaires on the basis of the definition of personal and room-specific applications?
To define the parameters of non-visual lighting effects (Research Question 1), we must consider that, in the last 20 years, with more activities in the time range between 2016 and 2022, numerous international scientific discussions and analysis tasks have been carried out to specify the right metric for non-visual effects using datasets from experiments of different research groups on nocturnal melatonin suppression as a validation basis. Rea et al. defined and improved the C S model (Version 2018) in the two years 2020 and 2021, wherein they took into account the exposure time t (in hours) and the visual field, as well as modeled the contributions of the i p R G C channel, the S cones, the rods, and the ( L + M ) channel [24,25,26]. This improved formulation was validated using melatonin suppression datasets [26]. The measurement method for the C S parameter in the 2018 version was published by Truong et al. [27]. As well, the International Commission on Illumination (CIE) [28] and scientists in the fields of neurophysiology, sleep research, and lighting technology have recently made efforts to find effective parameters based on the evaluation of the five photoreceptor signals ( L M S cones, rods, and i p R G C cells) and the corresponding calculation tool [29]. This scientific process resulted in the definitions of “Melanopic Equivalent Daylight (D65) Illuminance (mEDI)” and the “Melanopic Daylight (D65) Efficacy Ratio (mDER)”. These definitions have been recognized and used by international experts for several years and form the basis of Section 2 of this paper. Additionally, similar issues have been addressed in the WELL Building Standard v2, Q3 2020 version [30,31], which includes recommendations for melanopic equivalent daylight illuminance ( m E D I 136 , 109, 218, or 180 depending on the type of light sources). To answer Research Question 2, we recognize that international sleep researchers and neuroscientists made some recommendations in 2022 regarding non-visual light effects, physiological aspects, sleep, and wakefulness based on their literature reviews [32]. The physiological aspects included hormonal regulation cycles, heart rate, core body temperature, and certain brain activities. The main findings were as follows:
(a)
Daytime light recommendations for indoor environments: Throughout the daytime, the recommended minimum melanopic EDI is 250 lux at the eye measured in the vertical plane at approximately 1.2 m height (i.e., the vertical illuminance at eye level when seated)”.
(b)
Evening light recommendations for residential and other indoor environments: During the evening, starting at least 3 h before bedtime, the recommended maximum melanopic EDI is 10 lux measured at the eye in the vertical plane approximately 1.2 m height. To help achieve this, where possible, white light should have a spectrum depleted in short wavelengths close to the peak of the melanopic action spectrum”.
(c)
Nighttime light recommendations for the sleep environment: The sleep environment should be as dark as possible. The recommended maximum ambient melanopic EDI is 1 lux measured at the eye”.
The first recommendation, with a minimum value of “melanopic EDI = Melanopic Equivalent Daylight (D65) Illuminance, mEDI” of 250 lx measured vertically at the observer’s eye, is relevant and of great interest for the professional sector during the day (offices, industrial halls, educational and health facilities). The other two recommendations are for the dark hours in the home. For the solution to Research Question 3 “Determination and measurement of input parameters for non-visual lighting effects”, which builds the focus of this present paper, two application areas can be targeted:
  • The determination and measurement of non-visual parameters after the completion of new lighting installations and comparison with the specifications of the previous lighting design; or verification of the results of the development of new luminaires for HCL lighting in the lighting laboratory or in the field. For this purpose, absolute spectroradiometers for spectral radiance or spectral irradiance are used to calculate the parameters m E D I and m D E R at different locations in the building. These two non-visual parameters cannot be measured directly with integral colorimeters and illuminance–luminance meters.
  • The control or regulation of modern semiconductor-based lights (LED-OLED) and window systems (daylight systems) with the help of sensors in the room (e.g., presence sensors, position sensors, light and color sensors) [33]. In order to achieve a predefined value of non-visual parameters such as m E D I at a specific location in the room (e.g., in the center of the room or at locations further away from the windows), while taking into account dynamically changing weather conditions and the whereabouts of the room users, in practice, relatively inexpensive but sufficiently accurate optical sensors ( R G B sensors) are required. The goal is to obtain not only the target photometric and colorimetric parameters such as illuminance E v (in lx), chromaticity coordinates (x, y, z), or correlated color temperature ( C C T in K), but also the non-visual parameters m E D I and m D E R . The principle of this Smart Lighting concept using R G B color sensors is illustrated in Figure 2. The methods for measuring non-visual parameters with low-cost but well-qualified R G B color sensors are the content of the Section 3 and Section 4 and the focus of this paper.
The notations “integrative lighting” and “human-centric-lighting” used in this publication have been defined by the I S O / C I E publication [34] and mean a lighting concept and practice that integrates both visual and non-visual effects and produces physiological as well as psychological benefits for light users. The terms “smart lighting” or “intelligent lighting” are used in the same way and mean a framework that combines the above-mentioned integrative lighting as a goal and content of lighting practice with the technological aspects of lighting (signal communication, Internet of Things, sensor systems, control electronics, software that includes the methods of artificial intelligence, L E D luminaires) while taking into account the individual needs of light users. In order to achieve the above objectives, the content of this paper is structured as follows:
  • Section 2 defines the non-visual parameters according to the C I E publication [28].
  • Section 3 describes the characterization and signal transformation mathematics for converting the R G B signals of the R G B sensors to the chromaticities x, y, and z according to the C I E publication [35] for a viewing angle of 2°.
  • In Section 4, a computational model is established and analysed that converts the chromaticity z into the non-visual parameter m D E R and the combination of z and illuminance E v into m E D I . Section 5 summarizes the results with some outlook.

2. Definition of Non-Visual Input Parameters [28]

The extensive mathematical treatment and definition of the non-visual input parameters are described in detail in the C I E publication S 026/E:2018CIE System for Metrology of Optical Radiation for ipRGC-Influenced Responses to Light” [28], and those contents are summarized in this paper. The purpose of the C I E publication is to define spectral response curves, quantities, and metrics to describe optical radiation for each of the five photoreceptors in the human eye that may contribute to non-visual processes. The spectral sensitivity functions of the five receptor types are specified as follows (see Figure 3):
  • S-cone-opic, s 10 ( λ ) ;
  • M-cone-opic, m 10 ( λ ) ;
  • L-cone-opic, l 10 ( λ ) ;
  • Rhodopic, V ( λ ) , rods;
  • Melanopic, m e l ( λ ) , ipRGCs.
The C I E defined the α -opic quantities in Table 2 with the aid of the spectral sensitivities in Figure 3, and these quantities include the spectral radian flux Φ e , λ ( λ ) , the spectral radiance L e , λ ( λ ) , and the spectral irradiance E e , λ ( λ ) ( α represents one of the five photoreceptor types).
From this point, further parameters can be derived for assessing the non-visual effect of light, where S α ( λ ) is the spectral sensitivity function of one of the five receptor types, K m is the known photometric radiation equivalent (=683.002 lm/W), and Φ D 65 , λ ( λ ) is the spectrum of standard daylight type D65 (see the Equations (4) and (5)). The quantity “ K α , v ” is the α -opic efficacy of the luminous radiation.
K α , v = Φ e , α ( λ ) · S α ( λ ) · d λ K m · Φ e , α ( λ ) · V α ( λ ) · d λ
Equation (4) is similar to Equation (3.4) on Page 4 of [28].
K α , v D 65 = Φ D 65 , α ( λ ) · S α ( λ ) · d λ K m · Φ D 65 , α ( λ ) · V α ( λ ) · d λ
Equation (5) is similar to Equation (3.7) on Page 5 of [28]. For melanopic efficacy ( α = mel), the parameter value K m e l , V D 65 equals 1.3262 mW/lm. Consequently, the parameters m E D I and m D E R can be set up as shown in Table 3. If one wants to interpret the meaning of m E D I or m D E R for lighting engineering, there are two main aspects:
  • m E D I is the illuminance of the standard daylight illuminant D65 that has as much melanopic efficacy as the test light source for a given illuminance E v (lx) caused by the test light source. See Equation (6) in Table 3.
  • m D E R is the ratio of the illuminance of the standard illuminant D65 ( m E D I ) to the illuminance of the test illuminant E v (in lx) when the absolute melanopic efficacy of both illuminants is set equal. See Equation (9) in Table 3.
The melanopically effective parameters m D E R and m E D I can be calculated from the Equations (3), (5), (6), and (9) if the spectral irradiance E e , λ ( λ ) is known from a spectroradiometric measurement. The key question in this publication is whether these two parameters can also be determined with sufficient accuracy in practical lighting technology using a well-characterized and inexpensive R G B sensor in the sense of intelligent lighting technology (smart lighting). Section 3 deals with the description of the R G B sensor and how the R G B sensor signals can be transformed into tristimulus values ( R G B ) and into chromaticity coordinates ( x , y , z ) by means of a comprehensive spectral analysis.

3. RGB Color Sensors: Characterization and Signal Transformation

After the definition of the non-visual input parameters [28] in Section 2, we can understand the main concepts and mathematical forms of m E D I and m E D R . In Section 4, the prediction of a simple computational model for the non-visual quantities m E D I and m D E R was carried out. As well, the verification of the feasibility of a R G B sensor using this model was implemented to check the synthesized prediction model for lighting applications. This was also performed to confirm that the model not only conformed to the mathematical calculations, but also to determine whether it can be applied in the lighting systems with low-cost R G B color sensors. Therefore, Section 3 must attempt to describe R G B color sensors, their characterization, and appropriate signal transformation techniques. This section serves as a link between Section 2 and Section 4, as well as to balance the essential material for the verification of the prediction model using a color sensor as an example in Section 4 later. As a demonstration of the methodology for processing R G B sensors, it is not necessary to collect all color sensor sensitivities, synthesize, and compare them, but the most important thing here is to prove that the methodology can work well in this approach. Consequently, the future work can be implemented more comprehensively for different color sensor families through other research.

3.1. Characterisation of R G B Color Sensors

The color sensors (see Figure 4) include arrays of individual sensors or a group of sensors based on semiconductors with silicon as the semiconductor material, which are most commonly used for the visible spectral range between 380 nm and 780 nm.
In order to generate the corresponding spectral sensitivity curves of the individual color channels in the red, green, and blue range (so-called R G B channels), optical color filters based on thin-film technology (interference filters), absorption glass filters or color varnishes were microstructured and applied to the respective silicon sensor. For the correct spatial evaluation of the optical radiation according to the cosine law, a so-called cos prefix (usually a small diffuser plate made of optical scattering materials) was applied to the R G B sensor. The photons of each wavelength are absorbed by the R G B sensor, and photon currents are generated, which are then converted into voltages in an amplifier circuit. These voltages are then digitized in different bit depths (8 bits, 12 bits, 16 bits) by an A/D converter (analog–digital). The spectral sensitivity of each color channel RGB is, therefore, made up of the components shown in the Equations (10)–(12).
R ( λ ) = S ( s i l i c o n , λ ) · τ ( d i f f u s o r , λ ) · τ ( color lacquer R , λ ) · K R = C V R ( λ ) E ( λ )
G ( λ ) = S ( s i l i c o n , λ ) · τ ( d i f f u s o r , λ ) · τ ( color lacquer G , λ ) · K G = C V G ( λ ) E ( λ )
B ( λ ) = S ( s i l i c o n , λ ) · τ ( d i f f u s o r , λ ) · τ ( color lacquer B , λ ) · K B = C V B ( λ ) E ( λ )
where:
  • S ( s i l i c o n , λ ) is the spectral sensitivity of the silicon sensor;
  • τ ( d i f f u s o r , λ ) is the spectral transmittance of the cos diffusor;
  • τ ( color lacquer i , λ ) is the spectral transmittance of the color filter layer for the respective color channel i = R, G, B;
  • K R , K G , and K B are absolute factors to account for current-to-voltage conversion and voltage digitization;
  • C V i ( λ ) values are output signals (analog or digital) of the respective RGB color sensor at wavelength λ .
If a certain spectral irradiance E ( λ ) is present on the R G B sensor, three corresponding output signals R, G, and B are generated in the three color channels R, G, and B, respectively (see the Equations (13)–(15)).
R = E ( λ ) · R ( λ ) · d λ
G = E ( λ ) · G ( λ ) · d λ
B = E ( λ ) · B ( λ ) · d λ
For the calculation of the R G B channel signals, it is, therefore, necessary to know or to determine the spectral sensitivity function of each individual color channel R, G, B by laboratory measurements. Consequently, the spectral apparatus for determining the spectral sensitivity of the semiconductor sensors in the authors’ light laboratory (see Figure 5) consists of a high-intensity xenon ultrahigh-pressure lamp, a grating monochromator (spectral half-width Δ λ = 2 nm, spectral measuring steps for the entire spectrum between 380 nm and 780 nm, Δ λ = 2 nm), and an integrating sphere for homogenizing the quasi-monochromatic radiation coming out of the monochromator. The R G B color sensor and a known calibrated reference sensor are located at two different locations on the inner wall of the sphere behind a shutter.
During spectral measurement of the R G B sensor, one can determine the spectral sensitivity for the red color channel, for example, according to Equation (16).
R ( λ ) = C V R ( λ ) E ( S p h e r e , λ )
with C V R ( λ ) as the output signal of the red color channel and E ( sphere , λ ) as the spectral irradiance at the inner wall of the sphere when the monochromator is set to the wavelength λ , which in turn can be determined from the measured output photocurrent i ( reference , λ ) and the known spectral sensitivity of the reference sensor S ( reference , λ ) . As a result, we can write Equation (17).
E ( sphere , λ ) = i ( reference , λ ) S ( reference , λ )
The absolute spectral sensitivity of each color channel R, G, B can be determined from the Equations (16) and (17). Figure 6 and Figure 7 show examples of the spectral response curves of some R G B color sensors measured in the authors’ light laboratory. In Figure 6, the two datasets of the same type, SeS1 and SeS2, differ by the different peak heights of the R channel and by the different slopes of the G and B curves, because the two sensors S e S 1 and S e S 2 came from different production sets. The S e S 1 and S e S 2 spectral sensitivity curves in Figure 6 are fundamentally different from the R G B spectral sensitivity curves of the color sensor type in Figure 7.
The most accurate characterization possible for R G B sensors and the best possible similarity between C I E color matching functions C I E x y z ( λ ) for a field of view of 2° and R G B ( λ ) sensitivities are very important and have a direct influence on the later determination of non-visual parameters such as the results in Table 13. The smaller the difference between R G B ( λ ) sensitivities and C I E color matching functions C I E x y z ( λ ) for a field of view of 2°, the closer the distance between the values R G B S e S . , the chromaticity coordinates x y z C I E , and later the values x y z S e S . after matrix transformation of the chromaticity coordinates x y z C I E . Finally, the difference between the determined z S e S . = ( 1 x S e S . y S e S . ) and the C I E chromaticity coordinates z C I E . = ( 1 x C I E . y C I E . ) must be as small as possible, and this is understood as the primary uncertainty of the non-visual parameter calculation. Then the non-visual parameter model m E D R = f ( z S e S . ) in Section 4 and the parameter m E D I = [ m E D R ( z S e S . ) · E v ] determine the secondary uncertainty of the non-visual parameter calculation. The abbreviation “SeS.” indicates that it is the value of the sensors.

3.2. Method of Signal Transformation from R G B to X Y Z

Figure 6 and Figure 7 show that the R G B curves of the real color sensors differ more or less strongly from the x y z curves of the C I E color matching function for a field of view of 2° [35]. In order to obtain the X Y Z tristimulus values, the generated R G B signals must be transformed into digital form by matrixing (Equation (18)).
X SeS . , i Y SeS . , i Z SeS . , i = m 1 , 1 m 1 , 2 m 1 , n m 2 , 1 m 2 , 2 m 2 , n m 3 , 1 m 3 , 2 m 3 , n · R 1 , i G 1 , i B 1 , i .
the sensor technology, the amplification electronics, and the A/D conversion often have a more or less pronounced non-linear behavior, it may be necessary to use different matrixing types from 3 × 3 to 3 × 22. See Table 4.
If the sensor electronics are linear and the spectral response curves have a similar relative shape to the x y z color matching functions or L M S sensitivity spectra of the retinal photoreceptors, the 3 × 3 to 3 × 8 matrices will yield contributions from R G B signals as a function of the first term. The more the sensor electronics deviate from linear behavior, the more the shape of the spectral curves of the real sensors deviate from the x y z color matching functions. Therefore, matrices with R G B contributions in quadratic or cubic functions (see Table 4) should be set up. The procedure for finding the optimal matrix based on the chromaticity difference Δ u v is shown in Figure 8.
The goal of the optimization is to get the chromaticity (x, y, z) of the R G B sensor as close as possible to that derived from the C I E calculation. Instead of working with the non-uniform x y color diagram, it is better to work with the more uniform u v diagram (CIE, 1964). Bieske [37] showed in her dissertation that, when Δ u v is less than 10 3 , the color difference is not perceived by the subjects. From 10 3 to 3 × 10 3 , the color difference can be perceived but is still acceptable. If this value is higher than 5 × 10 3 , the color difference is unacceptable to the subjects. Based on this scientific contribution, the optimization with u v not only gives the best z-value, but can also be directly checked and compared with the obtained perception thresholds.

3.3. Matrix Transformation in Practice and Verification with a Real R G B Color Sensor

For the color sensor type whose spectral sensitivity curves are shown in Figure 7, the optimization processes are carried out according to the scheme in Figure 8 with a series of matrix types from 3 × 3 to 3 × 22 with nine different lamp spectra (as a training set). These nine lamp spectra are tabulated in Table 5 and shown in Figure 9. The illuminance for the optimization of the matrix was chosen to be E v = 750 lx. In this training set, the following lamp types were selected to reflect the variety of types used in practice today: fluorescent lamps between 2640 K and 4423 K with many spectral lines (see Figure 9), a halogen incandescent lamp, two phosphor-converted white L E D s, and a combination of R G B L E D s and warm white L E D s (RGBWW4500K).
The calculation with different matrix types and with individual lamp spectra yields color differences Δ u v of varying magnitudes, with the 3 × 3 matrix form yielding the smallest color difference (see Table 6). Optimization with more lamp types did not bring any significant improvement in this case. The final 3 × 3 optimized matrix for the transformation from R G B to X Y Z , as well as the formula for the calculation of the illuminance from the R G B signals, are shown in Table 7.
To verify the prediction quality of the formula in Table 7, eight spectra were selected in the next step. These were white phosphor converted L E D s (pc-LEDs) and a mixture of daylight spectra with white phosphor converted L E D s (TL-pc-LEDs), as is often the case in practical indoor lighting (see Figure 10). The illuminance was again set to 750 lx. The results of the check are shown in Table 8. From Table 8 it can be seen that
(a)
The deviation of the illuminances, calculated directly from the lamp spectra via the R G B color sensor signals and via the formula in Table 7, was below 0.65% in percentage terms;
(b)
For the majority of pc-LEDs and for the combinations daylight–white L E D s, the chromaticity difference Δ u v was in a small or moderately small range from the point of view of practical lighting technology, although the matrix transformation was synthesized only for the general case of many lamp types. No specific procedure was optimized for the specific case of only mixed light between daylight and conventional light, but the color and illuminance difference was still very small when using the achieved matrix transformation. Therefore, the authors did not build a separate case of mixed light with specific matrix transformation. Also, pc-LEDs and their mixing with daylight are dominant nowadays, so this verification in this way makes sense for applications. Therefore, the extreme cases of very complicated spectral shapes are not necessary to investigate their feasibility, because they are very rarely applied in today’s life. And an exception is the spectrum R12GB12-5000 K (green curve in Figure 10) with three distinct peaks of the three R G B L E D s (B around 455 nm, G around 525 nm, R around 660 nm), which is a very good demonstration example of a complicated spectral form. In this case, the color difference was Δ u v = 10.1 × 10 3 .
However, for the case of pure daylight spectra, a special matrix was found and described in Table 9 as a specific case. For the investigated R G B color sensor type, 185 daylight spectra were measured absolutely on a sunny summer day (19 August 2020 in Darmstadt, Central Europe) from 6:32:00 in the early morning to 20:35:32 in the late evening (see Figure 11). The chromaticity and illuminance values were calculated and used as a training set for the matrix. The largest color difference with a 3 × 3 matrix was found to be Δ u v = 1.2 × 10 3 , which was recorded in the last minutes of the evening before sunset when the correlated color temperature was very high, measuring approximately in the range of 17,000 K on the day of measurement. Most of the color differences were well below this value.
To verify the 3 × 3 matrix for daylight, 24 daylight spectra measured on a overcast day (23 September 2020 in Darmstadt, Germany) were considered. The verification results are shown in Table 10 with 9 out of 24 spectra as examples. There, the chromaticities and illuminances determined from the measured spectra are compared with the data obtained by matrixing the R G B color sensor (processed data). The calculated color difference Δ u v was very small and lay in the range of 10 4 .
The determination of the optimal matrix for the transformation of R G B sensor signals into X Y Z values and the illuminance E v described so far, as well as the verification results with an actual R G B color sensor type, prove that it is possible to obtain the colorimetric and photometric values X Y Z , E v , and C C T with relatively good results in the sense of a reliable, relatively accurate and adaptive lighting technology using inexpensive and commercially available color sensors. The next section deals with the accurate determination of the non-visual input quantities m E D I and m D E R indoors (daylight and artificial light combined or separated) and outdoors during the day (with daylight only) from the tristimulus values X Y Z and the illuminance E v of a qualified R G B color sensor, because, in practical lighting technology, a spectroradiometer is often too expensive and too cumbersome to handle.

4. Prediction of a Simple Computational Model for the Non-Visual Quantity mDER and Verification of the Feasibility of a RGB Sensor Using this Model

Once the matrices for electric light sources and for daylight spectra have been found, the tristimulus values X Y Z , the chromaticity coordinates x y z , the color temperature C C T , and the illuminance E v can be obtained from the R G B color sensors used (see the Table 7, Table 8, Table 9 and Table 10). To obtain the melanopic equivalent daylight illuminance m E D I , one only needs to know the value of m D E R and the illuminance E v (in lx) according to Equation (7), where the illuminance can be formed from the R G B signals (see the formula in Table 7 and Table 9) for determination. Consequently, the most important task is to calculate the dimensionless input quantity m D E R (see Equation (9)) from the chromaticity coordinates x, y, z (z turned out to be the most predictive, so it was the coordinate used).
For this transformation from z to m D E R , a large number of measured and simulated light source spectra were analyzed (see Table 11). The first four light source groups were real measured light sources from thermal radiators (28 light sources), compact and linear fluorescent lamps (252 light sources), different LED configurations (419 light sources), and 185 measured daylight spectra. These 884 light source spectra are shown in Figure 12. In order to simulate the combination of daylight with white LEDs or daylight with white fluorescent lamps, which can correspond to the lighting conditions of interiors with white LEDs or white fluorescent lamps with windows over a longer period of use, the authors of this paper took 185 measured daylight spectra and mixed them with a market-typical white LED spectrum or with a typical fluorescent lamp spectrum in nine different mixing ratios from 10% to 90%. This can be seen in Figure 13 for the case of combining daylight spectra with white LEDs. For the fluorescent lamps, it was similar, so a graphical representation of the mixing spectra is omitted here. The 185 daylight spectra with nine mixing ratios each resulted in 1665 simulated spectra each for the LED and fluorescent lamps (185 × 9 = 1665).
Using the 4214 light source spectra summarized in Table 11, 4214 m D E R values and 4214 z values were calculated, from which their relation was determined. The formula for m D E R is given in Equation (25).
m D E R = a · e b · z c · e d · z
The parameters a, b, c, and d for the fitting function and the correlation coefficients R 2 are given in Table 12. The course of the correlation between m D E R (ordinate) and the chromaticity coordinate z (abscissa) is shown in Figure 14, where the fitting function follows the green course of the curve.
To check the predictive quality of Equation (25) and the transformation from z to m D E R , the 4214 spectra in Table 11 were first used to calculate the 4214 m D E R values from the spectra themselves (denoted as m E D R [ C I E S 026 / E : 2018 ] ) and then another 4214 m D E R values indirectly via Equation (25), with z as the independent variable (denoted as m E D R = f ( z ) ) (see Figure 15). The goodness of fit ( R 2 ) was 0.97844, and the linear constants were 1.0214 and −0.05712. Compared to the ideal constants of 1 and 0, the quality of the fit is very good.
To verify the quality of the characterization of the R G B color sensor under test (see Figure 7), the quality of the “matrixing” in the transformation from R G B to X Y Z and E v (see Table 7) and the quality of the transformation from z to m D E R (see Equation (25)) and finally the calculation of m E D I (see the Equations (6) and (9)), we set each of the eight light source spectra in Figure 10 to 750 lx. From these 8 spectra, the chromaticity coordinates x and y and the quantities C C T , C R I R a , E v , m D E R o r i g i n a l , and m E D I o r i g i n a l were calculated. They are listed in Table 13. For these eight spectra, the chromaticity coordinates x n e w and y n e w , the value of E v , n e w , and the values m D E R n e w and m E D I n e w were also calculated from the R G B color sensor values by matrixing and transforming z to m D E R and m E D I , respectively. The maximum relative deviations Δ m D E R (in %) and Δ m E D I (in %) were found to be in the range of ±3.3%.

5. Conclusions and Discussion

In indoor lighting technology, including daylight components during the day, the daily task is to evaluate the lighting systems after the completion of the building or after reconstruction to assess the photometric and colorimetric quality of the newly developed interior luminaires and, in the context of intelligent lighting (smart lighting, H C L lighting, integrative lighting), to adaptively control the lighting systems in order to provide the room users with the best visual and non-visual conditions, as well as room atmosphere, at all times. In addition to the criteria of visual performance according to current national and international standards (e.g., EN DIN 12464 [9] or WELL [31]), aspects of psychological and emotional lighting effects, as well as the evaluation and adaptive control of lighting systems according to non-visual quality characteristics, should be considered. For this purpose, non-visual input parameters such as m E D I and m D E R need to be measured and processed with sufficient accuracy, reliability, and reasonable effort. Existing illuminance meters, luminance meters, color and luminance cameras, and small portable colorimeters are currently only capable of measuring photometric quantities such as illuminance E v , luminance L v , colorimetric parameters such as chromaticity coordinates x y z , and correlated color temperature ( C C T ). To measure the non-visual parameters, a transformation is needed to convert the tristimulus values to the non-visual parameters, and a sensor platform is needed to convert the R G B sensor signals to the tristimulus values X Y Z and E v (in lx). The authors of this paper have characterized some exemplary R G B color sensors on the basis of laboratory measurements, calculations, and optimizations by creating one or more matrices for the optimal transformation of R G B sensor signals to the tristimulus values. For different spectral sensitivities of the R G B color sensors, different matrices with linear or higher order RGB-signals can be optimized. As an example, for a certain R G B sensor, a 3 × 3-matrix was found by means of optimization with nine different spectra to be the best transformation that delivered a maximal chromaticity difference of Δ u v = 0.0083 (see Table 6). This matrix was verified with eight white phosphor-converted L E D s (pc-LEDs) and mixture of daylight spectra with white phosphor-converted L E D s (TL-pcLEDs) (see Table 8). For the spectra of daylight only, a special matrix was found that allowed for a much lower chromaticity difference (maximal Δ u v = 0.000857 , see Table 10). The reconstruction of the illuminance (in lux) from the R G B signals delivered a difference in the order of 0.64% (see Table 8). As next steps, on the basis of extensive calculations, a transformation from the chromaticity coordinate z to m D E R and, via E v , also to m E D I , was achieved with good accuracy. Table 13 has shown the verification results with eight light source spectra, which were in the range of Δ m E D I = ± 3.3 %, which can be accepted from a practical point of view of lighting engineers, lighting designers, and luminaire system developers. This work and the methodology described therein enable an accurate and financially justifiable assessment of lighting installations according to non-visual criteria, which will become increasingly important in the coming decades. This paper followed the aim to deliver, firstly, a framework and a general method on how to use commercially available R G B sensors to qualify them for measuring the metrics for non-visual effects. The results had some limitations, because only one R G B sensor type has been used in the methodical frame of this paper. Other R G B sensor types may have other spectral, optical, and electronic characteristics with different dark currents, signal-to-noise ratios, and cross-talk effects, which should require other forms of matrices (e.g., matrices with higher orders of R G B signals). Further studies will be planned by the authors of this present paper, in which the accuracies of several commercial RGB sensor types will be compared to the measurement results of an absolute measuring spectroradiometer for different test conditions (e.g., outdoor and indoor lighting conditions with different mixture ratios of daylight and artificial light sources).

Author Contributions

Conceptualization, V.Q.T., P.B. and T.Q.K.; Data curation, V.Q.T.; Formal analysis, V.Q.T. and P.B.; Methodology, V.Q.T. and T.Q.K.; Software, V.Q.T.; Supervision, T.Q.K.; Validation, V.Q.T. and P.B.; Visualization, V.Q.T.; Writing—original draft, V.Q.T. and P.B.; Writing—review and editing, V.Q.T., P.B. and T.Q.K.; Project administration, V.Q.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors. The publication of the manuscript is supported by the Open Access Publishing Fund of the Technical University of Darmstadt.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed to support the findings of the present study are included this article. The raw data can be obtained from the authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rea, M.S.; Ouellette, M.J. Relative Visual Performance: A Basis for Application; SAGE Publications Sage UK: London, UK, 1991; Volume 23, pp. 135–144. [Google Scholar] [CrossRef]
  2. Rea, M.S.; Ouellette, M.J. Visual performance using reaction times. Light. Res. Technol. 1988, 20, 139–153. [Google Scholar] [CrossRef]
  3. Weston, H.C. The Relation between Illumination and Visual Efficiency-The Effect of Brightness Contrast. In The Relation between Illumination and Visual Efficiency-The Effect of Brightness Contrast; H.M. Stationery Office: London, UK, 1945. [Google Scholar]
  4. Boyce, P.R. Human Factors in Lighting; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar] [CrossRef]
  5. Lindner, H. Beleuchtungsstärke und Arbeitsleistung–Systematik experimenteller Grundlagen (Illumination levels and work performance–systematic experimental principles). Zeitschrift für Die Gesamte Hygiene und Ihre Grenzgebiete (J. Hyg. Relat. Discip.) 1975, 22, 101–107. [Google Scholar]
  6. Schmidt-Clausen, H.J.; Finsterer, H. Beleuchtung eines Arbeitsplatzes mit erhöhten Anforderungen im Bereich der Elektronik und Feinmechanik, 2nd ed.; Fachverlag NW in Carl Ed. Schünemann KG, 1989; ISBN-10: 388314925X, ISBN-13: 978-3883149257. [Google Scholar]
  7. Blackwell, H.R. Contrast thresholds of the human eye. J. Opt. Soc. Am. 1946, 36, 624–643. [Google Scholar] [CrossRef] [PubMed]
  8. Commission International de l’Eclairage (CIE). An Analytical Model for Describing the Influence of Lighting Parameters upon Visual Performance. 1981. Available online: https://cie.co.at/publications/analytic-model-describing-influence-lighting-parameters-upon-visual-performance-2nd-0 (accessed on 21 April 2023).
  9. DIN Deutsches Institut für Normung. DIN EN 12464-1 Light and Lighting—Lighting of Work Places—Part 1: Indoor Work Places. (2019). Available online: https://www.din.de/en/getting-involved/standards-committees/fnl/drafts/wdc-beuth:din21:302583817?destinationLanguage=&sourceLanguage (accessed on 21 April 2023).
  10. Houser, K.W.; Tiller, D.K.; Bernecker, C.A.; Mistrick, R. The subjective response to linear fluorescent direct/indirect lighting systems. Light. Res. Technol. 2002, 34, 243–260. [Google Scholar] [CrossRef]
  11. Tops, M.; Tenner, A.; Van Den Beld, G.; Begemann, S. The Effect of the Length of Continuous Presence on the Preferred Illuminance in Offices. In Proceedings CIBSE Lighting Conference; 1998; Available online: https://research.tue.nl/en/publications/the-effect-of-the-length-of-continuous-presence-on-the-preferred- (accessed on 21 April 2023).
  12. Juslén, H. Lighting, Productivity and Preferred Illuminances: Field Studies in the Industrial Environment; Helsinki University of Technology Helsinki: Espoo, Finland, 2007. [Google Scholar]
  13. Moosmann, C. Visueller Komfort und Tageslicht am Bueroarbeitsplatz: Eine Felduntersuchung in neun Gebaeuden; KIT Scientific Publishing: Karlsruhe, Deutschand, 2015. [Google Scholar]
  14. Park, B.C.; Chang, J.H.; Kim, Y.S.; Jeong, J.W.; Choi, A.S. A study on the subjective response for corrected colour temperature conditions in a specific space. Indoor Built Environ. 2010, 19, 623–637. [Google Scholar] [CrossRef]
  15. Lee, C.W.; Kim, J.H. Effect of LED lighting illuminance and correlated color temperature on working memory. Int. J. Opt. 2020, 2020, 3250364. [Google Scholar] [CrossRef]
  16. Fleischer, S.E. Die Psychologische Wirkung Veränderlicher Kunstlichtsituationen auf den Menschen. Ph.D. Thesis, ETH Zurich, Zürich, Switzerland, 2001. [Google Scholar]
  17. Bodrogi, P.; Brückner, S.; Krause, N.; Khanh, T.Q. Semantic interpretation of color differences and color-rendering indices. Color Res. Appl. 2014, 39, 252–262. [Google Scholar] [CrossRef]
  18. de Boer, J.B.D.F. Interior Lighting; Macmillan: London, UK, 1978. [Google Scholar]
  19. Balder, J. Erwunschte Leuchtdichten in Buroraumen (Preferred Luminace in Offices). Lichttechnik 1957, 9, 455–461. [Google Scholar]
  20. Khanh, T.Q.; Bodrogi, P.; Trinh, Q.V. Beleuchtung in Innenräumen–Human Centric Integrative Lighting Technologie, Wahrnehmung, Nichtvisuelle Effekte; Wiley-VCH: Berlin, Germany, 2022. [Google Scholar]
  21. Lok, R.; Smolders, K.C.; Beersma, D.G.; de Kort, Y.A. Light, alertness, and alerting effects of white light: A literature overview. J. Biol. Rhythm. 2018, 33, 589–601. [Google Scholar] [CrossRef] [PubMed]
  22. Souman, J.L.; Tinga, A.M.; Te Pas, S.F.; Van Ee, R.; Vlaskamp, B.N. Acute alerting effects of light: A systematic literature review. Behav. Brain Res. 2018, 337, 228–239. [Google Scholar] [CrossRef] [PubMed]
  23. Dautovich, N.D.; Schreiber, D.R.; Imel, J.L.; Tighe, C.A.; Shoji, K.D.; Cyrus, J.; Bryant, N.; Lisech, A.; O’Brien, C.; Dzierzewski, J.M. A systematic review of the amount and timing of light in association with objective and subjective sleep outcomes in community-dwelling adults. Sleep Health 2019, 5, 31–48. [Google Scholar] [CrossRef] [PubMed]
  24. Rea, M.S.; Figueiro, M. Light as a circadian stimulus for architectural lighting. Light. Res. Technol. 2018, 50, 497–510. [Google Scholar] [CrossRef]
  25. Rea, M.S.; Nagare, R.; Figueiro, M.G. Predictions of melatonin suppression during the early biological night and their implications for residential light exposures prior to sleeping. Sci. Rep. 2020, 10, 14114. [Google Scholar] [CrossRef] [PubMed]
  26. Rea, M.S.; Nagare, R.; Figueiro, M.G. Modeling circadian phototransduction: Quantitative predictions of psychophysical data. Front. Neurosci. 2021, 15, 615322. [Google Scholar] [CrossRef] [PubMed]
  27. Truong, W.; Trinh, V.; Khanh, T. Circadian stimulus–A computation model with photometric and colorimetric quantities. Light. Res. Technol. 2020, 52, 751–762. [Google Scholar] [CrossRef]
  28. Commission International de l’Eclairage (CIE). CIE S 026/E: 2018: CIE System for Metrology of Optical Radiation for ipRGC-Influenced Responses to Light; CIE Central Bureau: Vienna, Austria, 2018. [Google Scholar] [CrossRef]
  29. Commission International de l’Eclairage (CIE). CIE S 026. 2020: User Guide to the α-Opic Toolbox for Implementing; CIE Central Bureau: Vienna, Austria, 2020. [Google Scholar] [CrossRef]
  30. Houser, K.W.; Esposito, T. Human-centric lighting: Foundational considerations and a five-step design process. Front. Neurol. 2021, 12, 630553. [Google Scholar] [CrossRef] [PubMed]
  31. International WELL Building Institute pbc. The WELL Building Standard, Version 2. New York, NY, USA, 2020. Available online: https://v2.wellcertified.com/wellv2/en/overview (accessed on 25 April 2023).
  32. Brown, T.M.; Brainard, G.C.; Cajochen, C.; Czeisler, C.A.; Hanifin, J.P.; Lockley, S.W.; Lucas, R.J.; Münch, M.; O’Hagan, J.B.; Peirson, S.N.; et al. Recommendations for daytime, evening, and nighttime indoor light exposure to best support physiology, sleep, and wakefulness in healthy adults. PLoS Biol. 2022, 20, e3001571. [Google Scholar] [CrossRef] [PubMed]
  33. Fan, B.; Zhao, X.; Zhang, J.; Sun, Y.; Yang, H.; Guo, L.J.; Zhou, S. Monolithically Integrating III-Nitride Quantum Structure for Full-Spectrum White LED via Bandgap Engineering Heteroepitaxial Growth. Laser Photonics Rev. 2023, 17, 2200455. [Google Scholar] [CrossRef]
  34. ISO/CIE TR 21783:2022|ISO/CIE TR 21783; Light and Lighting—Integrative Lighting—Non-Visual Effects. ISO: Geneva, Switzerland, 2022.
  35. Commission Internationale de l’Éclairage (CIE). Colourimetry, 2018. Available online: https://onlinelibrary.wiley.com/doi/10.1002/col.22387 (accessed on 4 April 2023).
  36. Westland, S.; Ripamonti, C.; Cheung, V. Computational Colour Science Using MATLAB; John Wiley & Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  37. Bieske, K. Über die Wahrnehmung von Lichtfarbenänderungen zur Entwicklung Dynamischer Beleuchtungssysteme. Ph.D. Thesis, Technische Universität Ilmenau, 2010. Available online: https://d-nb.info/1002583519/04 (accessed on 25 April 2023).
Figure 1. Input variables, influencing factors, and output variables in a comprehensive view of the effects of light on humans.
Figure 1. Input variables, influencing factors, and output variables in a comprehensive view of the effects of light on humans.
Sensors 23 05000 g001
Figure 2. Principle of a Smart Lighting concept with R G B color sensors.
Figure 2. Principle of a Smart Lighting concept with R G B color sensors.
Sensors 23 05000 g002
Figure 3. Five relative photoreceptor sensitivities in CIE S 026/E:2018 [28].
Figure 3. Five relative photoreceptor sensitivities in CIE S 026/E:2018 [28].
Sensors 23 05000 g003
Figure 4. 4 Examples of R G B color sensors. From left to right, the color sensor chips MTCS-CDCAF, MRGBiCS, and S11059 are shown. (Image Source: TU Darmstadt).
Figure 4. 4 Examples of R G B color sensors. From left to right, the color sensor chips MTCS-CDCAF, MRGBiCS, and S11059 are shown. (Image Source: TU Darmstadt).
Sensors 23 05000 g004
Figure 5. Schematic arrangement for measuring the spectral sensitivity of R G B sensors.
Figure 5. Schematic arrangement for measuring the spectral sensitivity of R G B sensors.
Sensors 23 05000 g005
Figure 6. Examples of the spectral sensitivity curves of some R G B color sensors.
Figure 6. Examples of the spectral sensitivity curves of some R G B color sensors.
Sensors 23 05000 g006
Figure 7. Spectral sensitivity curves of another R G B color sensor type.
Figure 7. Spectral sensitivity curves of another R G B color sensor type.
Sensors 23 05000 g007
Figure 8. Flowchart for the optimization of a transformation matrix from R G B to X Y Z based on the chromaticity difference Δ u v [35].
Figure 8. Flowchart for the optimization of a transformation matrix from R G B to X Y Z based on the chromaticity difference Δ u v [35].
Sensors 23 05000 g008
Figure 9. Spectra of the lamp types for matrix optimization.
Figure 9. Spectra of the lamp types for matrix optimization.
Sensors 23 05000 g009
Figure 10. Lamp spectra of 8 light source types for the validation of the 3 × 3 matrix in Table 7.
Figure 10. Lamp spectra of 8 light source types for the validation of the 3 × 3 matrix in Table 7.
Sensors 23 05000 g010
Figure 11. Some daylight spectra on a sunny day in Darmstadt, Germany (on 19 August 2020).
Figure 11. Some daylight spectra on a sunny day in Darmstadt, Germany (on 19 August 2020).
Sensors 23 05000 g011
Figure 12. Spectra of 884 measured light source spectra.
Figure 12. Spectra of 884 measured light source spectra.
Sensors 23 05000 g012
Figure 13. Spectra of 185 phases of daylight, mixed with white LEDs in 9 mixture ratios from 10% until 90%.
Figure 13. Spectra of 185 phases of daylight, mixed with white LEDs in 9 mixture ratios from 10% until 90%.
Sensors 23 05000 g013
Figure 14. Correlation between m D E R (ordinate) and the chromaticity coordinate z (abscissa).
Figure 14. Correlation between m D E R (ordinate) and the chromaticity coordinate z (abscissa).
Sensors 23 05000 g014
Figure 15. Comparison of the 4214 m D E R values on the ordinate, calculated from the spectra themselves (denoted as m E D R [ C I E S 026 / E : 2018 ] ), with the 4214 m D E R ( = f ( z ) ) values calculated using Equation (25). R 2 = 0.97844 . The linear constants a = 1.0214 , b = 0.015712 compared to the ideal constants 1 and 0 and the Δ m E D R 97% confidence interval of 0.067802 show a good quality.
Figure 15. Comparison of the 4214 m D E R values on the ordinate, calculated from the spectra themselves (denoted as m E D R [ C I E S 026 / E : 2018 ] ), with the 4214 m D E R ( = f ( z ) ) values calculated using Equation (25). R 2 = 0.97844 . The linear constants a = 1.0214 , b = 0.015712 compared to the ideal constants 1 and 0 and the Δ m E D R 97% confidence interval of 0.067802 show a good quality.
Sensors 23 05000 g015
Table 1. Parameters for interior lighting in [10,11,12,13,14,15,16,17].
Table 1. Parameters for interior lighting in [10,11,12,13,14,15,16,17].
No.ParameterPreferred Values
1Horizontal illuminance E v , h in lux [11,12,13]Greater than 850 lx;
Recommended range: 1300 lx–1500 lx
2Correlated color temperature C C T in K [14,15]4000 K–5000 K
3Color rendering index C I E C R I [17]>87
4Indirect to total illuminance ratio γ [10,16]>0.6–0.8
Table 2. α -opic quantities in CIE S 026/E:2018 [28].
Table 2. α -opic quantities in CIE S 026/E:2018 [28].
α -Opic Quantities
ParameterEquationEquation No.
α -opic-radian flux
Φ e , α or Φ α = Φ e , α ( λ ) · S α ( λ ) · d λ
similar to (3.1) Page 4 [28]
α -opic-radiance
L e , α or L α = L e , α ( λ ) · S α ( λ ) · d λ
similar to (3.5) Page 4 [28]
α -opic-irradiance
E e , α or E α = E e , α ( λ ) · S α ( λ ) · d λ
similar to (3.6) Page 5 [28]
Table 3. Melanopic equivalent D65 quantities of CIE S 026/E:2018 [28].
Table 3. Melanopic equivalent D65 quantities of CIE S 026/E:2018 [28].
α -Opic Quantities
ParameterEquationEquation No.
Melanopic Equivalent
Daylight (D65)
Illuminance (mEDI)
in lx
m E D I = 380 780 S m e l . ( λ ) · E e , λ ( λ ) · d λ K m e l . , V D 65
similar to (3.9) Page 6 [28]
with α =mel.
m D E R = K m e l . , V K m e l . , V D 65 ( * )
Melanopic
Daylight (D65)
Efficacy Ratio ( m D E R )
m D E R = 1 K m e l . , V D 65 · 380 780 S m e l . ( λ ) · E e , λ ( λ ) · d λ E v
m D E R = 1 K m e l . , V D 65 · 380 780 S m e l . ( λ ) · E e , λ ( λ ) · d λ E v
similar to (3.10) Page 7 [28]
with α =mel.
(*) Note: For the parameter K x y (indices x, y according to the corresponding definitions), see Equations (4) and (5)
when α = mel, and apply Equation (3) when the calculated parameter is the mel-opic irradiance.
Table 4. Different matrix types for color sensing according to [36].
Table 4. Different matrix types for color sensing according to [36].
Nr.SizeContent
13 × 3[R G B]
23 × 5[R G B R G B 1]
33 × 7[R G B RG RB GB 1]
43 × 8[R G B RG RB GB R G B 1]
53 × 10[R G B RG RB GB R 2   G 2   B 2 1]
63 × 11[R G B RG RB GB R 2   G 2   B 2   R G B 1]
83 × 14[R G B RG RB GB R 2   G 2   B 2   R G B   R 3   G 3   B 3 1]
93 × 16[R G B RG RB GB R 2   G 2   B 2   R G B   R 2 G G 2 B B 2 R R 3   G 3   B 3 ]
103 × 17[R G B RG RB GB R 2   G 2   B 2   R G B   R 2 G G 2 B B 2 R R 3   G 3   B 3 1]
113 × 19[R G B RG RB GB R 2   G 2   B 2   R G B   R 2 G G 2 B B 2 R R 2 B G 2 R B 2 G R 3   G 3   B 3 ]
123 × 20[R G B RG RB GB R 2   G 2   B 2   R G B   R 2 G G 2 B B 2 R R 2 B G 2 R B 2 G R 3   G 3   B 3 1]
133 × 22[R G B RG RB GB R 2   G 2   B 2   R G B   R 2 G G 2 B B 2 R R 2 B G 2 R B 2 G R 3   G 3   B 3
R 2 GB R G 2 B RG B 2 ]
Table 5. Properties of the lamp types for matrix optimization.
Table 5. Properties of the lamp types for matrix optimization.
Lamp TypeTungsten HalogenXenon-2CFL 3000 KCFL 5000 KFL 627FL 645LED C3LLED c3NLED RGBWW4500
C C T (in K)276241002640442327854423264045804500
D u v ( · 10 3 )37.20.612.21.82.26.31.3−1.0
CIE R 9 85−8648−60−72−60−28−3936
CIE R a 976790686468676990
Table 6. Chromaticity differences Δ u v for different matrix types for the lamp types of the training set.
Table 6. Chromaticity differences Δ u v for different matrix types for the lamp types of the training set.
NameXenon-2CFL-930CFL-5KFL627FL645C3L-pcLEDsC3N-pcLEDRGBWW4500 KMax
Δ u v 3 × 3 · 10 3 7.47.71.54.51.50.422.78.38.3
Δ u v 3 × 5 · 10 3 4.84812411251161651
Δ u v 3 × 7 · 10 3 8.55.3 3.1 × 10 4 1.6 1.8 × 10 3 7.0206.420
Δ u v 3 × 8 · 10 3 6.932114011487.01648
Δ u v 3 × 10 · 10 3 2.94.5 4.4 × 10 4 0.92 8.1 × 10 4 12 2.8 × 10 5 3.812
Δ u v 3 × 11 · 10 3 7.9 2.8 × 10 7 2.2372.242128.142
Δ u v 3 × 15 · 10 3 7.110 9.8 × 10 4 7.8 3.5 × 10 3 5.4 × 10 5 1.76.710
Δ u v 3 × 16 · 10 3 8.45.7 5.9 × 10 7 3.0 4.1 × 10 3 2.51.91010
Δ u v 3 × 17 · 10 3 8.45.7 5.9 × 10 9 3.0 4.1 × 10 3 2.51.91010
Δ u v 3 × 19 · 10 3 7.29.1 3.6 × 10 2 4.3 3.6 × 10 3 9.3 × 10 6 5.77.59.1
Δ u v 3 × 20 · 10 3 7.29.1 3.6 × 10 2 4.3 3.6 × 10 2 9.3 × 10 6 5.77.59.1
Δ u v 3 × 22 · 10 3 6.68.9 4.8 × 10 4 2.2 1.8 × 10 3 2.3 × 10 4 3.27.28.9
Table 7. Optimum matrix (3 × 3) for the transformation R G B to X Y Z , as well as the formula for the calculation of the illuminance E v from the R G B signals.
Table 7. Optimum matrix (3 × 3) for the transformation R G B to X Y Z , as well as the formula for the calculation of the illuminance E v from the R G B signals.
E v = 683 × ( 1.0001 · E v , R G B 0.0066 )
E v (lx) from R, G, B
E v , R G B = 0.6802 · R + 0.3651 · G + 0.1751 · B
R 2 = 1.0 ; R M S E = 0.15
Matrix 3 × 3 in case of
the 9 light sources of training set



5.73 × 10 5 4.17 × 10 5 2.27 × 10 5 3.64 × 10 5 4.50 × 10 4 3.26 × 10 5 6.61 × 10 4 2.04 × 10 5 1.13 × 10 6
Table 8. Verification of optimal matrix (3 × 3) in Table 7 with 8 LEDs plus mixed daylight spectra.
Table 8. Verification of optimal matrix (3 × 3) in Table 7 with 8 LEDs plus mixed daylight spectra.
NamepcLED-WWG3pcLED-WWVG1vpcLED- WWEX4pcLED-NWVG1pcLED-NWEXCV2pcLED-CWEX2R12GB12-5000 KTL-pcLED-30%
C C T (K)28013105296946143942505950014391
CRI R a 84.0685.6394.4190.9193.0995.9989.7389.85
x0.44340.42170.42520.35690.37640.34390.34490.3644
y0.39290.38410.37650.36080.35500.35570.34950.3650
E v , n e w 752.29752.54754.17751.40754.79750.01747.60751.94
x x e w 0.44340.42170.42520.35460.37790.33680.33190.3665
y n e w 0.39290.38410.37650.36430.36980.35710.35290.3696
Δ u v 1.36 × 10 10 1.20 × 10 10 1.68 × 10 10 3.15 × 10 3 8.76 × 10 3 5.32 × 10 3 10.1 × 10 3 2.42 × 10 3
Δ E v in %0.310.340.560.190.640.00-0.320.26
Table 9. Optimal matrix (3 × 3) for transforming R G B to X Y Z and formula for calculating illuminance E v from R G B signals in case of pure daylight.
Table 9. Optimal matrix (3 × 3) for transforming R G B to X Y Z and formula for calculating illuminance E v from R G B signals in case of pure daylight.
E v = 683 × ( 0.9987 × E v , R G B 0.1124 )
E v (lx) from R, G, B
E v , R G B = 0.6333 × R + 0.4804 × G
R 2 = 1.0 ; R M S E = 0.04
Matrix 3 × 3 in the case of the 9
light sources of the training set
1.8558 1.6603 1.4322 1.1084 0.23047 0.53993 0.53135 0.97596 6.0956
Table 10. Verification of the 3 × 3 matrix for daylight spectra (trained on a sunny day, verified on a overcast day).
Table 10. Verification of the 3 × 3 matrix for daylight spectra (trained on a sunny day, verified on a overcast day).
Sampling Time07:27:1710:03:3111:06:0112:03:1913:05:5014:07:4015:10:1016:12:4119:14:58
C C T in K12,464803363135651591458535470817416,066
E v , m e a s u r e in lx401939723,94058,21239,46343,15071,00615,058461
y m e a s u r e 0.28310.30660.32930.34130.33640.33760.34490.30810.2705
E v , p r o c e s s e d in lx400.609403.0223,944.0758,205.5739,465.2043,149.5970,987.1315,070.62 460.41
Δ u v m e a s u r e c a l c u l a t e × 10 4 4.75 2.81 × 10 2 1.741.201.821.960.5964.038.57
Table 11. Measured and simulated spectra (4214 in total).
Table 11. Measured and simulated spectra (4214 in total).
No.Type of light sources and their parameter ranges
C C T = 2201 – 17,815 K; −1.467× 10 2 < D u v < 1.529× 10 2
0.0797 < z < 0.4792; 80 < C I E R a < 100; 0.2784 < M D E R < 1.46
1.Conventional incandescent lamps + filtered incandescent lamps
2.Fluorescent tubes + Compact fluorescent lamps
3.LED lamps + LED luminaires
4.Daylight (CIE-model + measurements)
5.Mixtures (DL+LED) − [ mixture ratio = 10–90% ]
6.Mixtures (DL+FL) − [ mixture ratio = 10–90% ]
Table 12. Parameters a, b, c, and d for the fit function in Equation (25) and the correlation coefficient R 2 .
Table 12. Parameters a, b, c, and d for the fit function in Equation (25) and the correlation coefficient R 2 .
Parameterabcd R 2
Value0.538062.3212−0.37585−0.858080.968
Table 13. Verification of the m D E R and m E D I values with the color sensor under test on the basis of the 8 spectra in Figure 10 set to 750 lx.
Table 13. Verification of the m D E R and m E D I values with the color sensor under test on the basis of the 8 spectra in Figure 10 set to 750 lx.
NamepcLED-WWG3pcLED-WWVG1pcLED-WWEX4pcLED-NWVG1pcLED-NWEXCV2pcLED-CWEX2R12GB12-5000 KTL-pcLED-30%
C C T (K)28013105296946143942505950014391
CIE R a 84.0685.6394.4190.9193.0995.9989.7389.85
x0.44340.42170.42520.35690.37640.34390.34490.3644
y0.39290.38410.37650.36080.35500.35570.34950.3650
E v (lx)750750750750750750750750
m D E R o r i g i n a l 0.460.510.530.750.670.830.830.71
m E D I o r i g i n a l 347.67383.55395.93562.91500.06625.10625.79532.40
E v , n e w (lx)752.29752.54754.17751.40754.79750.01747.60751.94
x n e w 0.44340.42170.42520.35460.37790.33680.33190.3665
y n e w 0.39290.38410.37650.36430.36980.35710.35290.3696
m E D R n e w 0.460.530.540.740.660.810.830.69
m E D I n e w 346.28396.00403.79554.58500.85604.47621.78521.19
Δ u v 1.36 × 10 10 1.20 × 10 10 1.68 × 10 10 3.15 × 10 3 8.76× 10 3 5.32 × 10 3 1.01 × 10 2 2.42 × 10 3
Δ E v in %0.310.340.560.190.640.00−0.320.26
Δ m E D R in %−0.72.91.4−1.7−0.5−3.3−0.3−2.4
Δ m E D I in %0.403.251.991.480.163.300.642.10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Trinh, V.Q.; Bodrogi, P.; Khanh, T.Q. Determination and Measurement of Melanopic Equivalent Daylight (D65) Illuminance (mEDI) in the Context of Smart and Integrative Lighting. Sensors 2023, 23, 5000. https://doi.org/10.3390/s23115000

AMA Style

Trinh VQ, Bodrogi P, Khanh TQ. Determination and Measurement of Melanopic Equivalent Daylight (D65) Illuminance (mEDI) in the Context of Smart and Integrative Lighting. Sensors. 2023; 23(11):5000. https://doi.org/10.3390/s23115000

Chicago/Turabian Style

Trinh, Vinh Quang, Peter Bodrogi, and Tran Quoc Khanh. 2023. "Determination and Measurement of Melanopic Equivalent Daylight (D65) Illuminance (mEDI) in the Context of Smart and Integrative Lighting" Sensors 23, no. 11: 5000. https://doi.org/10.3390/s23115000

APA Style

Trinh, V. Q., Bodrogi, P., & Khanh, T. Q. (2023). Determination and Measurement of Melanopic Equivalent Daylight (D65) Illuminance (mEDI) in the Context of Smart and Integrative Lighting. Sensors, 23(11), 5000. https://doi.org/10.3390/s23115000

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop