The general V(λ) mismatch index, , of a photometer quantifies the deviation between its determined relative spectral responsivity, , and its nominal responsivity, the spectral luminous efficiency function for photopic vision, V(λ). LED-based light sources now dominate in general lighting, displacing incandescent and fluorescent lamps. Furthermore, the calibration of photometers will very likely see the replacement of the reference CIE standard illuminant A with CIE reference spectrum L41. The article evaluates the consequences. Should we also change the definition of the general V(λ) mismatch index? Based on performance criteria, the individual indices are evaluated using various datasets with the help of statistical analyses. As a result, the authors can conclude that the current definition works very well even under changed calibration and application conditions and does not need to be changed. Other evaluated indices for a new index definition perform only slightly better in some cases, but do not generate generally better properties.
Physical measurements of photometric quantities such as luminance or illuminance are typically performed with photometers incorporating a single sensor, with a filter adapted to a spectral luminous efficiency function, or with a spectroradiometer, where the weighting is performed by subsequent calculation [1].
In general, the task of a photometer is to determine an integral photometric quantity, Y, from the spectral distribution (SD) of a light source, , by weighted integration, according to:
where Km ≅ 683lmW−1, and weighting is carried out with the spectral luminous efficiency function for photopic vision, [1]. Other spectral luminous efficiency functions such as for scotopic vision or for the CIE 10° photopic photometric observer can also be used, but they are not the focus of this article.
Using the spectral responsivity of the photometer, s , the luminous responsivity, sv, is usually measured during the calibration using a calibration light source, denoted hereafter with the subscript C according to [2]:
where is the SD of the calibration light source. The handling of the integral limits is explained in appendix
The term `spectral distribution' (`SD'), with the symbol S(λ), is used here in a general form for the concentration of any radiometric quantity such as irradiance or radiance as a function of wavelength [3, entry 17-21-029]. In most cases, only the relative SD is used. Currently, usually a source with a relative SD similar to the one of CIE standard illuminant A, SA(λ), is used for calibration [4].
Since the relative spectral responsivity of the photometer [3, entry 17-25-064], , can never be perfectly matched to the spectral luminous efficiency function for photopic vision, V(λ), the luminous responsivity, sv, of a photometer depends on the SDto be measured, and in practice it is important to evaluate the effect of the spectral mismatch on the measurement, and correct it, if possible.
Figure1 shows an example for the measured relative spectral responsivity, , of a photometer (, photometer #73 in the CIE S 025 dataset [5]), the V(λ) function and the difference between both curves. The slight difference of the function or mismatch to the function leads to errors in photometric measurements when the SDin the actual measurement is different from that used during calibration.
1.1.Quality indices
Quality indices used in photometric measurement are defined by the International Commission on Illumination (CIE) [2]. They are physical quantities with unit one, determined for specific conditions to describe specific properties of photometers. They are positive numbers usually expressed as a percentage (symbol %). For an ideal photometer, any of the defined indices is zero, meaning the selected property is perfectly realised.
For this article, the general V(λ) mismatch index, , is of interest. quantifies the mismatch between the relative spectral responsivity of a photometer, , and the spectral luminous efficiency function for photopic vision, V(λ). is used for a description of the photometric performance of photometers under general lighting conditions and for selecting photometers requiring minimal corrections under light sources with SDs different from the reference source in photometric measurements.
A perfect match between the relative spectral responsivity and V(λ) implies that the luminous responsivity of the photometer does not depend on the relative SDof the light source to be measured. Therefore, it can be said that the general V(λ) mismatch index indicates the sensitivity of the luminous responsivity of the photometer with respect to the variation of the relative SDof the measured light source.
is one of the most important indices for the classification of photometers [6]. In the communication between customers and manufactures, is often the only relevant quality index, and in several other publications, e.g. [7, 8], is used to define selected required properties of photometers for applications to achieve a small error due to spectral mismatch.
1.2.Original definition
Introducing the normalised spectral responsivity,
calculated using a weighting with the SD of CIE standard illuminant A, SA(λ) (see figure2), is defined as:
A detailed review of the historical development of and further aspects of indices in general can be found in Krüger etal [9].
In previous internal studies regarding the measurement of white phosphor-based LEDs, it has been observed that calibrating a photometer according to SA(λ) usually results in much larger mismatch errors measuring white LED-based light sources, than calibrating according to the CIE reference spectrum L41, SL41(λ). This was reflected in a better correlation of the general V(λ) mismatch index with the expected spectral mismatch when using CIE reference spectrum L41 for calibration and weighting.
2.Objective
Today, in the daily measurement practice of electrical lighting, almost only LED-based light sources are evaluated. Therefore, the calibration of photometers will very likely see the replacement of CIE standard illuminant A with CIE reference spectrum L41 [10]. Figure2 shows the SDs of the two illuminants, SA(λ) and SL41(λ) (also referred to as SL(λ)), respectively.
This replacement raises the question whether an update of the definition of the general V(λ) mismatch index is necessary. A modification can be performed either by using a different normalisation than SA(λ) in for the normalised spectral responsivity of the photometer, (see equation (3)), or by introducing a different type of function for assessing this mismatch.
In this article, different versions for the general V(λ) mismatch index are studied to assess if one of them might be more adequate when using CIE reference spectrum L41 as calibration SD. In addition, different evaluation functions are computed to rate the relationship between evaluated indices and their ability to describe the photometric performance of photometers, in terms of minimising the errors due to the spectral mismatch under general lighting conditions. The computation is carried out for all the evaluated indices based on large datasets with relative spectral responsivities of photometers and SDs of different light sources. Finally, the results will be summarised in the conclusion section.
After preliminary considerations, it was concluded that any index adequate as general V(λ) mismatch index should have at least the following characteristics:
- To be well correlated with the expected measurement deviation from the true value when measuring an arbitrary white LED-based light source spectrally different from the calibration light source. This deviation is called hereafter 'expected spectral mismatch deviation'.
- To be independent of a specific set of test light sources, i.e. its definition must not include a specific set of SDs.
- To be well defined with minimal requirements for the spectral bandwidth and spectral resolution.
- To be easy to implement and to be readily understood by manufacturers and users.
Definitions of the evaluated indices studied in this work, as alternative of , are given below. Note that they are not defined as percentages for readability.
3.1.Classic indices
Classic indices only differ from as defined in equation (4) by a normalisation factor, αX, of the relative spectral responsivity, which is defined as:
where is the normalised spectral responsivity of the photometer using a weighting with the SD of a general illuminant X, SX(λ). The following cases of a general V(λ) mismatch index are considered:
, in the case of weighting with the SDof CIE standard illuminant A, , the normalisation factor is:
In this specific case, the normalised spectral responsivity is denoted as :
, in the case of no SDweighting, the normalisation factor is:
In this case, the normalised spectral responsivity is denoted as :
, in the case of weighting with the SDof CIE reference spectrum L41, , the normalisation factor is:
In this case, the normalised spectral responsivity is denoted as :
, in this case of calculating the normalisation factor αMin in a way that is minimised with the help of an optimisation algorithm [11]:
The normalised spectral responsivity is then denoted as :
3.2.Fourier indices
In the following, evaluated indices based on the Fourier transform will be introduced. For this purpose, we define a difference function, :
Using this difference function and αX = αA according to equation(6), the current definition (see equation (4)) can be rewritten as:
For the following index definitions, the normalisation factor for the case of no SD weighting, αE, given by equation(8)will be used. Based on the difference function and its Fourier transform, the following indices can be defined:
, this index was introduced by Ferrero etal [12] as:
In this expression, is the Fourier transform of , νλ is the spectral frequency, and νλ,c is the cut-off spectral frequency. See appendix
, this index is defined using the concept associated to , but with a difference function, , based on the bandwidth-limited signal that is used in equation(16).
The inverse of the Fourier transform, , yields:
With the bandwidth-limited version of , we can calculate the index according to equation(15), leading to:
The values of are comparable with the usual values and, at the same time, are fully correlated to the values defined by equation (16).
, this index makes use of bandwidth limitation for based on the assumption that coloured LEDs will be measured. In this case, very sharp structures in do not affect the measurement result. A simple symmetric model for the SDof a coloured LED, , (with centre wavelength λ0 and the full-width-at-half-maximum Δλ), is used:
The index is calculated as:
where * denotes the symbol of the convolution operation.
4.Comparison of the evaluated indices
This section provides some basic information about the evaluated indices and the datasets used for evaluation.
4.1.Spectral data for the evaluation
The relative spectral responsivities of 120 photometers from CIE S 025 are used for the evaluation.
For the white phosphor-based LEDs, the database from the PhotoLED project [10, 13], corresponding to 1496 SDs, is used. In addition, the database with 97 white RGB-based LEDs from CIE S 025 is used for the RGB-type calculations.
Figure 3 shows the mean normalised spectral responsivity , the min/max band for the CIES025 dataset and the standard deviation, , of and . Figure 4 shows a histogram of the values for the photometer CIES 025 dataset.
4.2.Classic indices
The relative differences between the values of the classic indices and are small, similar over the whole interval of the dataset and usually less than 15 %. It can also be observed that the histograms of the different classic indices are very similar.
4.3.Fourier indices
In comparison with the classic indices, it can be seen that the Fourier indices generate bigger relative difference values with respect to . For , the mean relative difference lies approximately around 15 % and, for , around 60 % with relative differences of about 100 % in the maximum.
For the photometers used in this work, the histograms of the and values are shown in figure5. The values are significantly smaller than the values and there is a significant compression of the range used for the photometers, whereas the values are shifted to smaller values, too, but not in the same magnitude. The distribution is not shown in this figure, since its values are not directly comparable with the usual values, but its distribution is very similar to values.
Figure6 shows the behaviour of the function defined in equation(14) for the different approaches.
The variation of the cut-off frequency, νλ,C, results in significant differences of the index value. The more frequencies are used in the Fourier transform, the more equivalent is the back-transferred delta function to the original one, as shown in figure7. In practice, the definition of a quality index should not have parameters. Therefore, in the case of , the parameter νλ,C needs to be fixed to avoid a dependency of the index value on this parameter. The cut-off frequency used here is a good choice for phosphor-type white LEDs and needs to be higher for RGB-type white LEDs.
4.4.Correlation between the evaluated indices
The linear coefficients of correlation [14] between the different values of the evaluated indices for the photometer database were calculated for comparison. Figure8 shows that the correlation between the classic indices is very high. Only the Fourier indices and (used with the parameters above) have a lower correlation to the classic indices and could therefore have significantly different features compared to the current index.
The previously defined indices are evaluated by comparing them to the spectral mismatch correction factors (SMCF), defined in appendix
The SMCF values, Fi,j , are calculated using the relative spectral responsivities and the SDs in the datasets introduced in section4. Using a specific selection from the database, we have relative spectral responsivities of D photometers and SDs of L light sources.
According to the derivation in appendix
where is the SD of the calibration light source C and is the relative spectral responsivity normalised with respect to that calibration source.
To make the evaluations more concise, only the absolute value of the spectral mismatch correction factor minus one, , will be used for further evaluation. This is probably the most significant value for a specific measurement and will be referred to hereafter as 'absolute deviation':
5.1.Evaluation based on quantiles of the SMCF
When evaluating a lot of SDs in the dataset and using the absolute deviation, a quantile-based approach can be used to get one stable number for each photometer:
where q is the quantile (e.g. q = 0.05) used for evaluation and q+ is the short form for the upper qualtile 1-q.
To get one specific number summarising the information of all photometers, a linear regression is made
The expression describes a linear regression based on a least square approach for the data points calculating the slope m and the intersection n at the origin.
The slope can be used to estimate the possible contribution of the SMCF for a specific photometer if no correction is possible because the relative spectral responsivity of the photometer or the specific SDof the light source to be measured are not known. With the value of the photometer and the type of the SD(either of a phosphor-type LED or an RGB-type LED) one can estimate that with a probability of about 95 %.
5.2.Evaluation based on the standard deviation of the SMCF
As demonstrated earlier in Mantela etal [15], the standard deviation of the spectral mismatch correction factor of each photometer over all light sources can also be used as a performance indicator:
In case the photometer calibration takes place with CIE standard illuminant A, the subscript C is replaced by A. In case the photometer calibration takes place with CIE reference spectrum L41, the subscript C is replaced by L [16]. With the database from the PhotoLED project, we get figure9.
In figure9, the black data points represent the upper quantile of the absolute deviation for each photometer according to equation(24) with the linear regression according to equation(25) illustrated by the black dashed line. The slope 0.48 %/% means that we get 0.48 % SMCF for every per cent of our index value as the upper 95 % quantile.
The red points represent the standard deviation of the SMCF calculated for each photometer as defined in equation(26). The red line is the result of the linear regression with these standard deviations defined by equation(27). The slope of 0.15 1/% means that we get a 0.15 standard deviation of the spectral mismatch correction factor for every percent of the mismatch index.
The coefficient of correlation, r, is also indicated in the legend of the figure (figures9 and 10).
Using RGB-type white LEDs with SDs used in CIE S025 and applying the same analysis shown before, we get results shown in figure 10.
5.3.Summary of the evaluations
From the evaluations for both calibrations with CIE standard illuminant A (C = A) and with CIE reference spectrum L41 (C = L) and all evaluated indices, we get the coefficients of correlation shown in figures 11 and 12.
Selected results can be found in tables1 and 2.
Table 1.Selected properties of mismatch indices for the case calibration with CIE standard illuminant A (C = A).
SD set | Index | ||
---|---|---|---|
PhotoLED | 0.85 | 0.91 | |
0.85 | 0.92 | ||
0.84 | 0.97 | ||
RGB | 0.91 | 0.88 | |
0.93 | 0.87 | ||
0.80 | 0.68 |
Table 2.Selected properties of mismatch indices for the case calibration with CIE referemce source L41 (C = L).
SD set | Index | ||
---|---|---|---|
PhotoLED | 0.95 | 0.91 | |
0.96 | 0.92 | ||
0.96 | 0.97 | ||
RGB | 0.86 | 0.87 | |
0.86 | 0.87 | ||
0.67 | 0.68 |
Summarising all the data, we conclude that:
- It makes sense to focus the calculation on the determination of an upper limit using the absolute deviation and the regression line of the quantiles using equation(25); all evaluated indices make a meaningful and stable regression for phosphor-type white indicated by high r values.
- There is a significant difference in the r value for the current index definition of when changing the SDof the calibration light source from SC = SA to SC = SL whereas the behaviour of and is very similar. The advantage of is negligible for phosphor-type LEDs.
- For RGB-type white LEDs, the prediction based on is no longer useful. There is a decrease of the coefficient of correlation independent from the SDof the calibration light source. This can only be changed by using a higher cut-off frequency.
5.4.Dependency on the SDand photometer dataset
Furthermore, we need to analyse whether our values (e.g. in figures 11 and 12) are stable or depend on selecting the dataset for photometers and/or light sources we use for the calculation. To this end, we introduce the bootstrap method [17], varying the photometer relative spectral responsivities and SDs used in the evaluation, implemented with a Monte Carlo simulation. The bootstrap method is used here only to estimate the standard deviation for the coefficients of correlation to verify whether valid statements can be made with respect to the calculated coefficients.
Figures13 and 14 show that the coefficients of correlation of the evaluated indices are very stable. In this case, this means the estimation of the coefficients of correlation did not depend on the concrete selection of photometers or SDs from the dataset and are therefore reliable to evaluate the quality of the evaluated indices.
It can be seen, that the standard deviations for the coefficients of correlation using SDs of RGB LEDs are a little bit higher. But the small difference has no practical impact.
6.Measurement uncertainty evaluation
Another important point is the determination of the measurement uncertainty for the characteristic values of the proposed indices. This is important because, during the classification of a photometer [6], the index value and its measurement uncertainty are needed. In the case of the general V(λ) mismatch index, the expected value of depends on the measurement uncertainty of the spectral measurements for well-matched photometers, due to the absolute value function [18, 19]. For this purpose, a Monte Carlo simulation (MC) has to be carried out in principle.
In a first step, the sensitivities are only determined for the amplitude noise and the shift of the wavelength scale separately.
6.1.Impact of the mismatch at single wavelength positions
First, we analyse the effect of the mismatch at single wavelength positions for all the evaluated indices. To do this, the index values are calculated for an ideal photometer with a small change Δ = 0.01 in the relative spectral responsivity, , at only one wavelength position, λi
Based on the data shown in figure15, one can evaluate what influence a signal difference at a single wavelength position has for the different evaluated indices. The behaviour for , and is predictable, one can see the weighting function very clearly, but the characteristics for and are somehow unexpected.
6.2.Wavelength shift
The next step is a rough estimation of the wavelength sensitivity of the quality indices. This can be analysed by shifts of Δλ of the ideal photometer's relative spectral responsivity, . In the model used inside the Monte Carlo simulation, this is represented by a correlated wavelength uncertainty contribution
Using the data shown in figure16, the change of the evaluated indices is about 0.02for every nm correlated wavelength shift. This is a very critical value, which in practice means that great care must be taken to realise a very accurate wavelength scale when determining the relative spectral responsivities of photometers.
6.3.Measurement uncertainty contribution in photometric measurements due to the SMCF
Based on the investigations presented here, one can estimate the contribution of the SMCF in a measurement uncertainty budget for photometric measurements of white LEDs (phosphor-type and RGB-type) if only the value of the photometer is known and no further information about the relative spectral responsivity of the photometer and the SDof the light source is available.
The presented data and regression results in figure17 show that about 95 % of the SMCFs (absolute deviations and photometer calibration with CIE reference source L41) are smaller than for phosphor-type white LEDs and smaller than for RGB-type white LEDs with very high coefficients of correlation.
The calculations used in this article can be found in the Python package empir19nrm02 available on GitHub [20]. There, the Python source code based on the LuxPy package [21] for all the specific related calculations can be found. It contains all datasets used here for demonstration purpose. Furthermore, the package contains a Jupyter Notebook, generating all figures used in this article.
8.Discussion
This article describes and evaluates the performance of several indices for general mismatch under general lighting conditions based on LED light sources. They can be classified as classic indices (where the only modification with respect to the current is the normalisation factor for the relative spectral responsivity), and those indices derived by a Fourier transform called Fourier indices.
It was found that the linear correlation between the classic indices (, , and ) is very high, which would imply that any reasonable normalisation of the relative spectral responsivity should not impact the index performance. There is also a high correlation between the two Fourier indices ( and . This is expected, since the latter is based on the former.
The coefficients of correlation between the different indices and statistical parameters of the absolute value of the spectral mismatch correction factor minus one (the absolute deviation), , namely the quantiles and standard deviations, have been studied using experimental data of relative spectral responsivities and SDs of white LED light sources. The aim was to classify the different indices by performance, under the assumption that the statistical parameters of are related with an expected deviation from the true value due to the mismatch.
The initial observation that the SMCF is much better correlated to than to (see figure18) is only due to the different calibration conditions.
When CIE reference spectrum L41 is used for the calibration, we might conclude that the Fourier indices seem to have a better performance measuring phosphor-based LEDs (broadband SDs) compared to classic indices. However, this is not the case when evaluating general lighting based on RGB LEDs (narrowband SDs). Here, classic indices have significantly better performance compared to the Fourier indices using a predefined cut-off frequency.
When the CIE standard illuminant A is used for the calibration, and general lighting based on LEDs is evaluated, both types of indices have a bad performance.
This means we can make the following more or less equal proposals for the time after changing the calibration illuminant from CIE standard illuminant A to CIE reference spectrum L41:
- Since, in practice, the value of does not change significantly with other reasonable normalisations of the relative spectral responsivities, one option would be to keep the current definition without any change; this is the simplest option for customers and manufacturers;
- Change the quality index to ; this would be in line with the original idea, which is to include the weighting with the SDof the calibration light source (which will probably be CIE reference spectrum L41 in the future) into the calculation of the normalised relative spectral responsivity, ;
- An interesting option is to change the quality index to or ; it would make the index independent of the SDof CIE standard illuminant A; this would be coherent with the idea that the spectral matching of a photometer should be independent of any form of light source SD;
- Use as an alternative index for general lighting under phosphor-based LEDs; the advantage is that it correlates better with the in those conditions; the drawback is its more complex implementation. Furthermore, one needs to change the cut-off frequency to get a good correlation to RGB-type LEDs, too.
The community is expected to decide based on the proposals and the basis laid out in this study. The following overview may be helpful using '+/++' for pro arguments, '−' for contra arguments and 'o' for neutral (table3):
Table 3.Summarising the index properties for selection.
Index-type | ||
---|---|---|
Feature | Classic | Fourier |
Current use | + | − |
Correlation with SMCF | ||
Phosphor-type | + | ++ |
RGB-type | + | − |
Implementation | + | o |
Understanding | + | o |
It should be noted that the spectral mismatch correction factor allows the correction of the spectral mismatch deviation in the specific case of measuring a specific light source with a specific SD. In contrast, the general V(λ) mismatch index gives an indication of the mismatch but allows no correction. It is valid for describing the photometer's expected performance when measuring an arbitrary and unknown SDof a white light source. This information can be used for a first approximation in a measurement uncertainty budget. Using a photometer calibrated with CIE reference source L41, 95 % of all SMCFs (absolute deviation) are smaller than for phosphor-type white LEDs and smaller than for RGB-type white LEDs with a very high coefficient of correlation.
Furthermore, the relative spectral responsivity normalised with respect to the calibration light source C, (see equation(22)) should be provided by the manufacturer as additional information with the calibration certificates used to easily calculate the spectral mismatch correction factor in the application.
Acknowledgments
This project 19NRM02 RevStdLED has received funding from the EMPIR program co-financed by the participating states and from the European Union's Horizon 2020 research and innovation program. Ville Mantela acknowledges partial support by the Academy of Finland Flagship Programme, Photonics Research and Innovation (PREIN), decision number: 320167. Furthermore, we thank Kevin G Smet for successfully implementing all the photometric and colorimetric functions in the LuxPy Python package, which we used to implement the calculations presented here.
Special thanks from the authors go also to the Reviewer and Peter Zwick for their valuable advice in the preparation of the article.
The absolute luminous responsivity [3, entry 17-25-061], sv,C, of a non-ideal photometer with a spectral responsivity depends on the SDof the calibration light source, SC(λ), according to:
where Km ≅ 683lmW−1 and s0 is the normalisation factor for the absolute spectral responsivity . A light source Z, with a different SD, SZ(λ), has a different absolute luminous responsivity, , accordingly:
Note that:
- For all integrals using V(λ), except for the original definition of , the lower and upper integration limits are 360nm and 830nm because V(λ) is null outside this range.
- For all integrals using the photometer's spectral responsivity, the lower and upper integration limits are λminandλmax, the spectral limits for which the spectral responsivity of the photometer is not null. However, the integration limits to be considered should at least cover the wavelength range in which V(λ) is defined, so that λmin ⩽ 360 nm and λmax ⩾ 830 nm.
There are two possibilities here:
- The integration limits are left at according to the current definition of . Since especially the measurement in the interval comprises a large measurement uncertainty and the two intervals and hardly contribute to the SMCFs because V(λ) is very small in these intervals.
- The integration limits are extended to the range 360nm to 830nm, since is defined in this interval.
The absolute luminous responsivity, sv,Z, must be corrected using the spectral mismatch correction factor (SMCF, symbol ), which is the ratio of the two different luminous responsivities [2]:
Defining a more general normalisation of the relative spectral responsivity with respect to a calibration source C, , using the SD of an arbitrary calibration source C, SC(λ), as weighting function instead of the SD of CIE standard illuminant A, SA(λ), and by expanding the integration limits to the spectral ranges where both, V(λ) and srel(λ) , are not zero, one gets a modified expression of equation (3):
With this definition, it is possible to rearrange the spectral mismatch correction factor as:
Using equation(34) instead of equation(32) would simplify the spectral mismatch correction factor calculation if the manufacturer provides from the calibration. Furthermore, it shows the relevance of the normalisation in (33).
Appendix B.: Implementation details for the Fourier indices
: The discrete Fourier transform (DFT) 8 , of the difference function, , defined by equation(14), is used to calculate its representation in the frequency domain, , according to:
With N spectral values in the wavelength domain, δs,i , the DFT results in N values in the frequency domain, . The mean value and the variance with the basic rules for the DFT 9 are obtained through:
For the discrete (positive) frequencies νλ,i , one gets νλ,i = i/(N ⋅ Δλ) for i = 0,...,N/2. That means, for both the standard configurations where with 5nm and 1nm steps, the frequency step is 2.47 × 10−3 nm−1. In this article, only 5nm steps are used for the evaluation. However, comparisons of the results of the authors showed that the differences for the values using 1nm or 5nm steps are negligible.
Using a cut-off frequency νλ,c as parameter and the associated index ic = ⌊νλ,c ⋅ N ⋅ Δλ⌋, is calculated as:
This means that is the variance of a bandwidth-limited version of . Using the variance (without bandwidth limitation) as an index was first published by Erb [22].