Determination of tailored filter sets to create rayfiles including spatial and angular resolved spectral information

To simulate and optimize optical designs regarding perceived color and homogeneity in commercial ray tracing software, realistic light source models are needed. Spectral rayfiles provide angular and spatial varying spectral information. We propose a spectral reconstruction method with a minimum of time consuming goniophotometric near field measurements with optical filters for the purpose of creating spectral rayfiles. Our discussion focuses on the selection of the ideal optical filter combination for any arbitrary spectrum out of a given filter set by considering measurement uncertainties with Monte Carlo simulations. We minimize the simulation time by a preselection of all filter combinations, which bases on factorial design. © 2015 Optical Society of America OCIS codes: (120.4820) Optical systems; (350.4600) Optical engineering; (120.5240) Photometry; (120.5630) Radiometry; (220.2945) Illumination design. References and links 1. I. Ashdown, “Near-field photometry: a new approach,” J. Illum. Eng. Soc. 22(1), 163–180(1993). 2. I. Ashdown, and M. Salsbury, “A near-field goniospectroradiometer for LED measurements,” Proc. SPIE 6342, 634215(2007). 3. M. López, K. Bredemeier, N. Rohrbeck, C. Véron, F. Schmidt, and A. Sperling “LED near-field goniophotometer at PTB,” Metrologia 49(2), 141–145(2012). 4. R. Rykowski, “Spectral ray tracing from near field goniophotometer measurements,” Light. Eng. 19(1), 23– 29(2011). 5. OSRAM Opto Semiconductors application note, “Importing rayfiles of LEDs from OSRAM Opto Semiconductors”, (OSRAM Opto Semiconductors, 2013), http://www.osram-os.com/Graphics/XPic5/00165120 0.pdf, (accessed, July, 2015). 6. V. Jacobs, J. Audenaert, J. Bleumers, G. Durinck, P. Rombauts, and P. Hanselaer, “Rayfiles including spectral and colorimetric information,” Opt. Express 23(7), A361–A370(2015). 7. V. Jacobs, J. Audenaert, J. Bleumers, G. Durinck, P. Rombauts, and P. Hanselaer, “On spectral ray files of light sources using principal component analysis,” in Proceedings of 28th CIE Session 2015 (CIE, 2015), pp. 543–547. 8. IES,“Ray File Format for the Description of the Emission Property of Light Sources,” TM25-13(2013). 9. F. Reifergeiste and J. Lienig, “Modelling of the temperature and current dependence of LED spectra,” J. Light Vis. Environ. 32(3), 288–294(2008). 10. CVI Laser Optics, “Datasheet selection guide,” pp. 232–233, https://marketplace.idexop.com/store/Support Documents/12 FiltersandEtalons 2013 final-e.pdf, (accessed, July, 2015). #248854 Received 27 Aug 2015; revised 29 Sep 2015; accepted 30 Sep 2015; published 3 Nov 2015 (C) 2015 OSA 16 Nov 2015 | Vol. 23, No. 23 | DOI:10.1364/OE.23.029543 | OPTICS EXPRESS 29543 11. Newport, “Newport resource e-catalog,” pp. 655–658, http://www.nxtbook.com/nxtbooks/newportcorp/resource 2011/#/684, (accessed, July, 2015). 12. Schott Glass Technologies catalog, Catalog optical filters, (Schott Glass Technologies, 2015), http://www.schott.com/advanced optics/german/download/schott-optical-filters-2015-catalog-complete-en.pdf, (accessed, July, 2015). 13. U. Krüger, and F. Schmidt, “The impact of cooling on CCD-based camera systems in the field of image luminance measuring devices,” Metrologia 46(4), 252–259(2009). 14. TechnoTeam, “LMK 50 color datasheet,” (TechnoTeam, 2012), http://www.technoteam.de/product overview /lmk/products/lmk 50 color/index eng.html, (accessed, July, 2015). 15. H. M. Antia, Numerical Methods for Scientists and Engineers (Springer Science & Business Media, 2002). 16. Schott Glass Technologies Technical Information, “TIE-35: Transmittance of optical glass,” (Schott Glass Technologies, 2005), http://fp.optics.arizona.edu/optomech/references/glass/Schott/tie-35 transmittance us.pdf, (accessed, July, 2015). 17. Z. M. Zhang, T. R. Gentile, A. L. Migdall, and R. U. Datla, “Transmittance measurements for filters of optical density between one and ten,” Appl. Opt. 36(34), 8889–8895(1997). 18. Oriel Instruments, “Filter characteristics”, pp. 10-30–10-32, http://ecee.colorado.edu/ mcleod/pdfs/AOL/labs /10030.pdf, (accessed, July, 2015). 19. I. H. Blifford, “Factors affecting the performance of commercial interference filters,” Appl. Opt. 5(1), 105– 111(1966). 20. CIE Technical Report, CIE 198:2011: Determination of Measurement Uncertainties in Photometry (Commission internationale de l’Eclairage, 2011). 21. C. Hoelen, J. Ansems, P. Deurenberg, T. Treurniet, E. van Lier, O. Chao, V. Mercier, G. Calon, K. van Os, G. Lijten, and J. Sondag-Huethorst, “Multi-chip color variable LED spot modules,” Proc. SPIE 5941, 59410A(2005). 22. C.-C. Sun, I. Moreno, Y.-C. Lo, B.-C. Chiu, and W.-T. Chien, “Collimating lamp with well color mixing of red/green/blue LEDs,” Opt. Express 20(S1), A75–A84(2012). 23. D. C. Montgomery, Design and Analysis of Experiments, 7th Edition (Wiley, 2009). 24. R. F. Gunst, and R. L. Mason, “Fractional factorial design,” Wiley Interdiscip Rev: Comput Stat 1(2), 234–244 (2009).


Introduction
A lot of effort is put into enhancing the illumination quality, for instance color homogeneity, color mixture or other parameters related to actinic action spectra, in modern lighting technology.Especially color fringes and color mixture have become a problem due to spatially and angular varying spectral distributions of the light sources LED and OLED.Besides reducing the physical origins of those varying spectra, simulations of their interactions in an optical system to validate their impact have become more important.A typical case is the prediction of the color homogeneity appearance of a lighting system during a product developing process.Predicting those spectral related information requires much simulation time and realistic light sources models.The perfect light source model or rayfile would describe spatially and angular varying luminance and spectra.This information can be created by combining spectral and goniophotometric measurements [1][2][3].Rykowski [4] proposed different methods to combine goniophotometric measurements and spectral measurements.A state of the art technique, which describes conventional phosphor converted white LEDs, uses two spectrally sharp separated rayfiles for the LED and the phosphor [5].However, this separation does not work for any arbitrary spectrum.Jacobs et al. [6,7] combined a goniophotometric measurements and some spectral measurements at different angles and performed a principal component analysis (PCA) to describe varying spectra over angle.The idea is that a weighted sum of basis functions is capable of reconstructing the spectrum.This approach works for any arbitrary spectra and minimizes the amount of storage.Except for a few spectral basis functions just the angular varying amplitudes need to be stored.The basis function approach also fulfills the requirements to be included in the ray format TM25 [8].However, the measurement approach from [6,7] is not capable of describing spatially varying spectra if no hyperspectral camera is available.Additionally, this requires both additional hardware and in case of a tunable filter based hyperspectral camera a lot of measurement time [4].
In this work we extend the basis spectra approach to create light source models, which include spatially and angular varying luminance and spectra.The main focus of the article is a robust reconstruction of spectra described by a set of basis functions taking into account measurement uncertainties.The main idea is to use a minimum of time consuming goniophotometric measurements with different optical filters.We describe a method consisting of a filter preselection, which bases on factorial design, succeeded by a Monte Carlo simulation to determine the optimal filters with minimal error out of a given filter set.

Modeling LED spectra and reconstructing the amplitudes
The main assumption is that each spectrum R(x, y, z, φ , θ , λ ) which generally changes over the angles φ and θ as well as over the spatial dimensions x,y and z can be written as weighted sum of basis spectra S i (λ ).The relative intensity at each wavelength λ is assumed to be constant for each basis spectrum S i (λ ).Just their amplitudes A i (x, y, z, φ , θ ) vary as a function of angle or spatial dimension according to Eq. (1).
It is important to note some limitations of the constant basis spectra assumption.In the case of LEDs an unstable thermal environment or changing electrical parameters might cause a change in the relative basis spectra.Also a phosphor spectrum might change due to self absorption.It has to be ensured that all measurement data apply for the same steady state of the light source.Since our analysis focuses on LED based lighting we can access a large set of basis functions to model the individual emission spectra of the LEDs.The model functions used in this analysis are summed up in [9].Note that according to [9] some of them are capable of predicting spectral changes as a function of temperature or forward current, which could be included in the optical simulation and that the models are always greater or equal to zero, which may help to detect errors during the reconstruction process.
In LED based lighting the set of physical basis spectra S i generally consists of semiconductor emission spectra and phosphor emission spectra.We describe the phosphor distribution by a smoothed spline.The sum of all basis spectra is optimized with a constrained simulated annealing or a simplex algorithm.Therefore an optimization variable is necessary.Our optimization variable is defined as an optionally spectral weighted residual sum of squares (RSS) as described in Eq. ( 2).The weighting factors w j (λ ) may be set to w j (λ ) = 1 if all wavelengths have the same importance.Depending on the case they may be customized to account varying relative importance of individual wavelengths.In this case we aim to minimize chromaticity coordinate distances between the test spectra R(λ ) and the modeled sum of basis spectra R(λ ).Therefore we used the k = 3 color matching functions w 1 (λ ) = x(λ ), w 2 (λ ) = ȳ(λ ) and w 3 (λ ) = z(λ ) from the CIE 1931 2°Standard Observer.
Four different typical LED based lighting spectra R(λ ) are shown in Fig. 1.Note that all following analyses could be applied for any other arbitrary set of basis spectra.
In general, a spectral rayfile described by n basis spectra S i (λ ) can be created with n goniophotometric measurements and any arbitrary set of n different optical filters τ i (λ ) for example those shown in Fig. 2. By solving the system of linear equations according to Eq. (3) all n amplitudes A i , which change over angle φ , θ and over the spatial region x, y, z are reconstructed based on the measurement values M i .The variable τ sys (λ ) describes the transmission profile of the remaining measurement system such as objective, neutral density filters and CCD.
According to Eq. (3) both filter combinations visualized in Fig. 3 would be able to reconstruct the spectrum.However, after a real measurement the obtained results A i will differ by ΔA i due to the uncertainties η i of the measurement values M i .To account those uncertainties Eq. (3) has to be extended according to Eq. ( 4) .The matrix M Sτ is the interaction matrix of the basis spectra and the filters according to Eq. ( 3).
Therefore it is necessary to consider those uncertainties in the general reconstruction workflow by a filter selection as shown in Fig. 4.This filter selection process will be discussed in the following sections.

Modeling the uncertainty
The uncertainty η in M + η = (M + η τ ) • (1 + η N ) from Eq. ( 4) has two different physical origins.The first factor consists of uncertainties η τ regarding the filter functions τ i (λ ) and their interactions with the set of basis functions S i (λ ).The factor η τ depends on the individual filter and its technology such as optical glasses and interference filters.The evenly distributed filter set used in our analysis consists of 63 interference filters (three with 1.5 nm, 47 with 10 nm, seven with 40 nm and six with 70 nm FWHM) from [10,11] and 23 glass filters in the visible spectral range.The glass filters consist of the 18 edge filters GG380-RG715 and five band pass filters from [12].Typical internal transmission profiles, temperature coefficients and refractive indices were provided by the manufacturers.
The second term η N describes uncertainties, which apply for all filters independent on the specific filter function τ(λ ) at each step of the goniophotometric measurements.Examples are noise from the CCD chip and uncertainties of the remaining optical system τ sys (λ ).The noise from the CCD chip scales linear with the signal level of an individual pixel, if the dark signal noise, which is just in the order of a few least significant bits (LSB), is neglected [13].We assume the term η N to be Gaussian distributed in the order of the repeatability of a commercial luminance/color camera [14].
The uncertainties propagate while solving the system of linear equations.This error propagation can be estimated with the condition number of the matrix M Sτ according to Eq. ( 5) [15].Note that Eq. ( 5) is an approximation, which is not capable to distinguish between the different amplitude values ΔA i and is therefore not suited to judge the performance of a reconstruction as detailed as Eq. ( 2).However, it is suited as an rough estimation, especially if cond 2 (M Sτ ) is very large.In the following section we analyze empirical models for optical glass and common interference filters as they are provided by the manufacturers to approximate the physical caused measurement uncertainty η τ .

Optical glass filters
An optical filter glass is characterized by its internal transmission spectrum τ in (λ ) given at a reference thickness d 0 .The ideal transmission function τ(λ ) is the product of the resulting internal transmission τ in for a thickness d and the Fresnel losses τ Fr at the air/glass and glass/air interfaces, which depend on the angle of incidence α and the refractive index n r .The ideal relation is given in Eq. ( 6) [12].We added the uncertainties from the room temperature u(T ) and temperature shifts u(T abs ) due to absorption as well as spatial roughness u(d), angle of incidence u(α), the uncertainty of the transmission curve measurement u(τ in ) and its wavelength precision u(λ ) [16] as absolute values to the model.All uncertainties are assumed to be correlated and are therefore the same for all wavelengths.Note that all uncorrelated uncertainties of τ(λ ) are averaged out due to the integration on the CCD.The model of the transmission function therefore depends on the angle of incidence (AOI) α, the refractive index n r , the thickness of the glass d and the temperature coefficient dλ /dT .Hence the model is defined as ) using β = arcsin {sin[α + u(α)]/n r }.We do not model multiple reflections since the absolute difference between multiple reflection and a two boundary approach is small.Of course multiple reflections occur during the measurement of the final transmission τ(λ ) and therefore the absolute difference will be calibrated in the system for normal incidence.By neglecting the effect in our model we do not assign this difference to τ Fr (α, n r ) but treat it as a constant, which will not change as a function of the uncertainty terms.If the refractive index of the optical glass is 1.56, the difference of the transmission at normal incidence would be smaller than 0.25%.Therefore by not considering the effect as a functional relation we just neglect the change of those 0.25% difference as a function of the uncertainty terms.For the same reason we also neglect dispersion effects and thermal expansion [17].Also note that the linear temperature dependence in Eq. ( 7) is limited to edge filters [12].

Interference filters
The transmission spectrum τ(λ ) of interference filters relies on constructive and destructive interference at a stack of specially designed interfaces.They can be designed with different FWHMs at different center wavelengths.We apply the same uncertainty as described for the glass filters despite heating due to absorption.The shift of the transmission spectrum can be estimated for small angles (α < 15°), where n e is the effective refraction index of the spacer layer and λ 0 the wavelength of the shifted spectral feature [10,12,18,19].Decreasing transmission, the broadening of the transmission spectrum and a separation of s-and p polarization can be neglected for the angle of incidence used in this simulation [18].The model of the interference filters for small angles α < 15°can therefore be written as

Filter selection and Monte Carlo simulation
Given a set of m different filter functions τ(λ ) one can determine the ideal subset of n filters for each specific R(x, y, z, φ , θ , λ ) consisting of n basis spectra S i (λ ) with a Monte Carlo simulation according to Eqs. ( 4), (7) and (8).The uncertainty terms are modeled as described in Fig. 5.The shape of the individual distributions is chosen in agreement with [20].The angular distribution is modeled based on geometrical assumptions regarding the diameter of the spherical light sources (small: 20 mm, large: 70 mm), the circular filters diameter (40 mm) and the measurement distance (300 mm).We assume the filter to be in front of the objective.Both extended sources result in an individual angular distribution visualized in Fig. 5.To demonstrate our geometrical assumptions regarding the derivation of the angular distributions we added the resulting distribution of a point source.In the Monte Carlo simulations just the small and large angular distributions are used.During the Monte Carlo simulation 10 5 random amplitude sets modeled as rectangular distribution within a high dynamic range of 140 dB [14] are generated and reconstructed with all filter combinations from the 86 filters.
We performed four simulations for each test spectrum of Fig. 1 and varied both the uncertainty angular distribution and the signal noise of the CCD as shown in Fig. 5.All 10 5 random amplitude sets were reconstructed according to Eq. ( 4).As figure of merit for an individual spectrum we used the weighted RSS according to Eq. ( 2).We again weighted the RSS with the CIE 1931 2°Standard Observer color matching functions x(λ ), ȳ(λ ) and z(λ ) for the same reasons as described in section 2. As figure of merit for the whole filter combination the mean value of all RSS RSS mean is used.The filter combination minimizing RSS mean would be the preferred one for the goniophotometric measurements and the spectral reconstruction.
The histograms in Fig. 6 visualize the results of the Monte Carlo simulation for the reconstruction of test spectrum RGBY from Fig. 1 for three plausible filter combinations using u(α) = large and η N = 0.4%.The first filter combination performs best according to RSS mean .The histograms visualize the error as Δu v since the absolute values of RSS mean do not allow an intuitive judgment of the colorimetric precision.Those color coordinate differences are recently used for describing illumination color uniformity [21,22].Note that Δu v alone is not suitable to evaluate the performance of the reconstruction, if there are more than three basis spectra.There are in general infinite possibilities to adjust a color coordinate with more than three spectral sources.

Preselection of filters and possible filter combinations
If the filter set consists of m filters and the number of spectral sources is n, there are possible filter combinations.In our example this leads to over 2 × 10 6 combinations for test spectrum RGBY (n = 4) and even more for larger n.An individual Monte Carlo simulation for each possible filter according to Eqs. ( 7) or ( 8) and each possible combination according to Eq. ( 4) results in an very high computation time for the filter selection process.Therefore we developed a preselection method to reduce the number of filters and filter combinations to be addressed by a Monte Carlo simulation.
Given a specific test spectrum we first use a two level full factorial design [23] to estimate the physical filter errors according to Eqs. ( 7) and (8).A two level full factorial design describes the permutation of chosen minimal and maximal values for all k input factors resulting in 2 k iterations.As minimal and maximal factors u min and u max we used every uncertainty Δu'v' mean =6.3×10 -4   Δu'v' mean =3.16×10 -3   Δu'v' mean =1.54×10 -3   Wavelength (nm) u min = u mean − u σ and u max = u mean + u σ , that is the difference of the uncertainty factors from their mean value and the standard deviation derived from Fig. 5.There are six factors resulting in 64 iterations for each optical filter glass and four factors resulting in 16 iterations for each interference filter according to Eqs. ( 7) and (8).If more detailed models, which account more factors than Eqs.( 7) and (8), are used, the number of iterations can be reduced by using a fractional factorial design [24].We estimate the measurements standard deviation η τ for each filter from the results.
After estimating the uncertainty of each filter, the error propagation resulting from the solution of Eq. ( 4) has to be considered.Therefore the second part of our preselection process estimates the mathematical error of each filter combination with a 3-level full factorial design [23].We use the results of the first simulation process and add the CCD signal noise η N resulting in a minimal, mean and maximal value for each filter: Note M represents the measurement value without any uncertainty.The test spectrum is reconstructed for all 3 n possible permutations for each filter combination.Again the number of iterations could be reduced by either using a two level factorial design, a fractional factorial design or both.As figure of merit we use the RSS mean from the original spectrum and the 3 n reconstructed spectra like in the Monte Carlo simulation.After selecting a certain amount of the m n possible combinations, which performed best in the preselection process, we perform the Monte Carlo simulation.

Results and discussion
If possible, a robust filter selection tends to create a diagonal matrix M Sτ to minimize the error propagation of Eq. ( 4) according to Eq. ( 5).This is shown for example by filter combination 1 in Fig. 6.However, filter combination 2 and 3 show also that the relations of the spectral reconstruction, which bases on the interaction of different filter transmission profiles and their uncertainties with given basis spectra, are more complicated and not intuitive.Figure 6 also proves that a non-optimal filter selection may result in an unnecessary large uncertainty of the reconstructed spectrum.Therefore the Monte Carlo simulation should be performed.
The preselection is capable of reducing the overall amount of necessary Monte Carlo simulations.To validate the performance of our preselection we executed all simulations with a very large amount of possible filter combinations, that is 50% for test spectra BY, RBY and BBY and 10% for test spectrum RGBY since its large amount of possible combinations results in a high computation time.
Figure 7 shows four exemplary results out of the 16 simulations.The mean residual sum of squares RSS mean from each combinations' Monte Carlo simulation is plotted on a logarithmic scale against the preselection number N to check the validity of the preselection process.A perfect preselection would be a monotonically increasing graph.However, in practical implications, most important is that the ideal filter combination, or at least a combination which is nearly as good as the ideal one, is within the first part of the preselection.There must not be strong outliers towards small values at late preselection numbers since it would prohibit a more detailed analysis in the real Monte Carlo simulation.However, an overestimation of bad filter combinations is uncritical since it will be detected during the succeeding Monte Carlo simulation.The highlighted areas in Fig. 7 correspond to the first two percent of all filter combinations.All most promising combinations including the global minimum of RSS mean are this area.
We summed up the preselection number N of the ideal filter combination, its relative position out of all possible combinations in percent and its associated mean coordinate distance Δu v mean in Table 1.All relative positions are within the first few percent of the possible combinations.Note, that the results in percent are more critical for fewer possible combinations.However, one can for example save a factor 50 of computation time by just simulating the best 2% from the preselection.Instead of an ideal combination you may also search for a combination better than a given threshold which will save even more computation time if the simulation starts according to the preselection number.
It is important to note that the reconstruction error summed up in Table 1 results just from the uncertainties of the filter transmission profiles τ i (λ ) and the signal noise of the CCD.The absolute values of the uncertainty depend on other parameters as well.The tendency to create a diagonal matrix M Sτ increases if more uncertainties, which are basisally the same for the different filters τ i (λ ), are considered.The general concept of the preselection and the Monte Carlo simulation is not affected.
The importance of different uncertainty factors in Eqs. ( 7) and ( 8) depends on the interaction between the specific test spectrum and the specific filter functions as well as the filter technology and therefore differs for all 86 filters.In general the impact of uncertainties affecting the wavelength precision of the transmission function τ(λ ) strongly depends on the spectral position of steep slopes.If the filter and a basis spectrum both have a steep slope in the same spectral region, the impact on the measurement result and therefore the relative importance of those uncertainties increases.
We analyzed the main effects of the physical filter errors for the absolute input uncertainty values shown in Fig. 5 according to Eqs. ( 7) and ( 8) with the results of the full factorial design.It can be summarized that the most important effect for optical glass filters is the wavelength   precision of the transmission curve measurement u(λ ).The thermal induced wavelength shifts are less important but still in the same order of magnitude.This effect increases for edge filters with a higher edge wavelength since they have a higher temperature coefficient than the optical glass filters with a smaller edge wavelength.If the wavelength shifts overall have a small impact due to the relative positions of the slopes, the uncertainties of the angular distribution u(α) as well as the spatial roughness u(d) of 3 μm are as important as the given temperature uncertainties u(T ) and u(T abs ).The effect of the uncertainty from the amplitudes of the transmission curve measurement u(τ) is the least important since it is at an order of magnitude smaller.
In case of interference filters the angular uncertainty u(α) and again the wavelength precision of the transmission curve measurement u(λ ) are the most important factors.The impact of the temperature uncertainty u(T ) is approximately one order of magnitude smaller since the temperature coefficients of interference filters are at least an order of magnitude smaller than for optical glass filters.If the wavelength shifts overall have a small impact on the measurement results due to the relative positions of steep slopes, the uncertainty of the transmission curve measurement u(τ) is at least a few factors smaller but might still be in the same order of magnitude.It is important to note, that the relative importance changes if input uncertainties or other model parameters like the spatial roughness u(d) or the thickness of an optical glass filter change.

Conclusion and outlook
Based on the assumption that a test spectrum can be described as weighted sum of constant basis spectra, we developed a workflow for the fast measurement of angular and spatial resolved spectral rayfiles.The minimum number of goniophotometric measurements with different optical filters depends on the number of basis spectra used to model the source.Each goniophotometric measurement provides information about a spatially and angular varying spectral section.The whole spectrum is reconstructed as function of angle and position by solving the resulting system of linear equations.An ideal filter combination out of a given filter set or the ideal channels of a mounted hyperspectral camera based on tunable filters can be found by a Monte Carlo simulation of the filter combinations.If the Monte Carlo simulation becomes computationally expensive due to a large number of filter combinations, one can use a small preselection simulation which bases on a few iterations according to a factorial design to save computation time.
In a future work the simulation time might be further reduced by removing less important uncertainty factors to simplify the models.Therefore the significance of the main uncertainty factors and their interactions has to be obtained.Those significances can be estimated for the individual case by a further analysis of the results from the full factorial design [23].The absolute uncertainties in Table 1 just account for the spectral filters and the signal noise of the CCD.Therefore we will further investigate the different terms of uncertainty, for the creation of spectral rayfiles, including the spectral uncertainties of the remaining optical system and the basis spectra as well as spatial shifts resulting from different goniophotometric measurements and focus our studies on both experimental validation and a filter technology comparison.

Fig. 5 .
Fig. 5. Model parameters of the uncertainties.The shape of the probability density function (PDF) is chosen in agreement with [20].

Fig. 6 .
Fig. 6.Histograms visualizing the reconstruction error derived from 10 5 Monte Carlo simulations for test spectrum RGBY with the large angular distribution and η N = 0.4% for three different filter combinations.

Fig. 7 .
Fig. 7. Validation of preselection: The graphs visualize four of the performed Monte Carlo simulations.The reconstruction performance RSS mean of different filter combinations is plotted on a logarithmic scale against their preselection number N. The highlighted area corresponds to the first two percent of the preselection.The minimum equals the optimal combination. ) )

Table 1 .
Validation of preselection:The table shows the preselection number N of the optimal combination, its relative position in percent and its associated Δu v mean of all simulations.