 Research
 Open Access
 Published:
Energy balance in Spectral Filter Array camera design
Journal of the European Optical SocietyRapid Publications volume 13, Article number: 1 (2017)
Abstract
Background
Multispectral imaging permits to capture more spectral information on object surface properties than color imaging. This is useful for machine vision applications. Transmittance spectral filter arrays combined with a solid state sensor form an emerging technology used for snapshot acquisition. In spectral filter arrays technology, the sensitivities of the camera have critical consequences, not only on applications, but also in the viability of the system. We discuss how to balance the energy of each channel in single exposure multispectral imaging.
Methods
We propose a methodology to design filters that permits to reach an optimal balance of energy. We apply this method on practical illuminations combined with a Gaussian model of transmittance filters.
Results and discussion
Our results demonstrate that we can optimize energetically the global camera response with few efforts for several cases of illumination environments, for a given sensor and a number of spectral channels. This methodology can be embedded in an applicationoriented optimization framework.
Conclusions
This methodology enables achieving a great range of optical design, and it can be embedded in an applicationoriented optimization framework.
Background
Multispectral imaging (MSI) is a sensing technique that is employed in different spectrum of wavelengths. When used in the context of typical silicon sensors, it may span the electromagnetic spectrum from ultraviolet to nearinfrared (NIR) range, approx. [3201100] nm. Up to now, multiple sensors or exposures were necessary to record such images [1–3]. The generalization of color filter arrays (CFA) into spectral filter arrays (SFA) [3] permits the use of a single sensor and a single exposure to record a multispectral image. This solution can be encapsulated into a compact and affordable MSI solution, easy to handle, which complies with a typical camera pipeline. The principle of SFA imaging is to combine a single sensor with a filter set. Each spatial position of the array captures only one spectral channel similarly to the Bayer filters [4], and the filter transmission characteristics determine the sensor spectral properties. Several realizations of such setup have been shown recently in the literature [5–7]. An example is provided in Fig. 1. Various research directions have been explored: filter characteristics [8–11], adequate full resolution image processing, e.g. demosaicing [12–17], and other aspects [18, 19].
In silicon based solidstate sensor image acquisition, when a photodiode receives too much radiant energy, the maximum number of electrons per pixel is reached [20]. This refers to saturation, and extra information is lost. In addition, the overloaded photosite distributes charges into adjacent pixels, generating blooming (especially for CCD technology) [21, 22]. On the other hand, if a photodiode does not receive sufficient photons, the signal to noise ratio (SNR) of the electrical signal generated is too low compared to the global noise level [23]. Thus, when integration time is fixed, the sensor performance is optimal only within a given range of intensities according to its sensitivity. In a graylevel intensity sensor, pixels have a similar sensitivity and it is relatively easy to set an optimal integration time [24, 25].
However, in the case of a SFA setup, different sensitivities are present on the same chip and may lead to over and underexposed classes of pixels, depending on the scene and on the illumination. This is illustrated by Thomas et al. [7], where the sensitivities of the lesssensitive band and the mostsensitive band (see Fig. 1) show a relative efficiency ratio in the interval [11%, 99%] for a set of practical filters, an ideal viewing illuminant and a perfect diffuser. In such cases, it is more likely that a single optimal integration time shall generate either saturated values or underexposed pixels, as indicated by results reported in their article. Figure 1(e) shows a raw image captured with an SFA camera. On a zoom on Fig. 1(e), we can observe a severe discrepancy between pixel intensity values for a uniform patch.
One might consider several solutions to solve this problem. One solution would be to set an optimal integration time for each pixel. This would be possible by using a multiple layer technology, see e.g. [26]. To the best of our knowledge, this approach is still in development and remains an expensive solution, which does not yet meet the requirements of a low price SFA. Another solution is to use high dynamic range imaging (HDR) based on multiple exposure, see e.g. [27, 28]. Although the solution may be embedded in a realtime framework, it precludes a singleshot SFA capability. Another solution would be to apply gains to rebalance the underexposed channel, but the SNR would still remain at a similar level despite denoising potential facilities. One might consider this problem in the sensor design and optimize the camera sensitivities for a particular condition. Such an optimization may be constrained by a specific application or by manufacturing constraints in filter characteristics. Energy balance of SFA has been previously considered by Péguillet et al. [29]. This paper develops such a solution within a constrained optimization framework.
Optimization of spectral filters have been extensively studied in the color imaging field for either accurate color acquisition or reproduction, by e.g. Vhrel [30, 31]. Sharma et al. [32] addressed the problem of obtaining a good perceived colorimetric information, under several viewing illuminants. A set of 3 to 7 filters have been designed using their optimization framework, and a linear minimum mean squared error estimator. Simulation shows that their technique performs better and with an improved SNR than the solution initiated by Vhrel [31]. They concluded that the use of 4 bands seems to be the optimal solution to get a good accuracy when recording color under multiilluminants. This conclusion is shared by Sadeghipoor et al. [33] who proposed an algorithm to optimize the spectral sensitivity function of the camera by finding the best demosaicing matrix for color acquisition. The constraint imposed is that the resulting function contains smooth practical filters, which are assumed to be physically realizable. Filters have been also optimized by Monno et al. [34] in the context of SFAs in order to provide the best demosaicing results. In their study, they used Gaussian models to approximate sensor sensitivities.
In this work, we propose to base on these previous works and to develop a methodology that provides the best filter bank for a singleshot SFA camera, while considering energy balance. The main contribution is in the definition of what is energy balance and how this may be handled within a more comprehensive optimization process that includes higher level information processing (such as demosaicing, spectral signature indicators or spectral reconstruction).
We use a Gaussian model to approximate transmittance filters and to provide simulated examples. The methodology can be extended to any optical model that can describe the transmittance of filters in a few meaningful parameters, e.g. hyperGaussians or multilobe single cavity FabryPérot, as suggested by one of the anonymous reviewer. The proposed methodology is based on an imaging model and on an empirical set of constraints such as:

1.
We must avoid overexposed and underexposed areas when shooting a scene with the multispectral camera.

2.
We must guarantee the best postprocessing efficiency (e.g demosaicing or spectral reconstruction) by homogenizing independent channel noise and preserving the design of the filter according to the application.

3.
We want to guarantee system versatility over a set of viewing illuminants.
Methods
Acquisition model and balanced sensor
We consider a SFA camera that consists of N spectral channels characterized by F, a set of spectral transmittance filters, such as F={f _{1},…,f _{ N }}. The spectral sensitivity of the sensor without a SFA is represented by S(λ). For a given illuminant l(λ), and an object of reflectance O(λ), we obtain a set of values Π={ρ _{1},…,ρ _{ N }}. Every value ρ _{ i } is given by filter f _{ i }, i∈⟦1,N ⟧ such as in Eq. 1.
This equation could also include the transmission of the medium in which the system is used, e.g. underwater imaging.
We define energy balance of a sensor for one illuminant when, for a perfect diffuser, the variance of Π is equal to 0, in other words, when ρ _{1}=ρ _{2}=…=ρ _{ N }. The flat reflectance object may be replaced by a specific reflectance property characterizing the application. For compactness, we denote Π=X(F), a function where X define a vector Π from a set of filters F, according to the imaging model of Eq. 1.
Definition 1
In order to guarantee energy balance, we want to find F among the general set of possible transmittance curves T, such as it minimizes V(Π), the variation among components of Π. If we note Y(F)=V(X(F)), this can be summed up as finding F such as min∀F∈T Y(F).
To generalize this definition to represent practical illuminant conditions for a range of applications, we introduce a set of illuminants L=[l _{1},…,l _{ M }] assumed to be representative of lightning conditions. This set may be dependent on the application. Then we rewrite Eq. 1 as Eq. 2, with j∈⟦1,M ⟧.
Definition 2
Again, for a perfect diffuser, we want to find an optimal set F, which minimize the variation among components of Π ^{′}. Π ^{′} is a set of Π _{ j }, such as Π ^{′}={Π _{1},…,Π _{ M }}. Similarly, we note Y ^{′}(F)=V(X ^{′}(F)), where X ^{′} is a set of functions according to Eq. 2 and we can write the problem as finding F such as min∀F∈T Y ^{′}(F).
The solution where this converges might not be suitable for multispectral acquisition. Indeed, an obvious solution to this problem is when F is composed of the same filters (i.e. f _{1}=f _{2}=…f _{ N }). To ensure the practicality of the filter set, we constrain the set F to a reduced subset of T. The reduced set is comprised of Gaussian filters distributed over the considered part of the spectrum. Then, F is characterized by its set of peak sensitivity locations Ω=μ _{1},…,μ _{ N }, which is fixed by the application. We propose to give a freedom to two parameters in order to find an optimal F using the standard deviation of each Gaussian σ _{1},…,σ _{ N } and the maximum transmittance B _{1},…,B _{ N }. This model could be adapted to several models of transmittance, but for clarity of communication, it is easier to work with Gaussian filters.
According to Definition 1, one may still find a minimum convergent solution. However, the set might not be viable in practice, due to limited efficiency or unpractical width. According to Definition 2, we can not converge to zero due to the variability in illuminations. But we can give a degree of freedom, ε, such as Y ^{′}(F)+ε=0.
In this context, we must put constraints on the model, in order to ensure the practicability of the filters set (and its realization), a reasonable variation among the set and relatively good maximum transmittance. These set of constraints are to be defined in next Sections.
Sensor considerations
A major element in this system is the shape of the sensor spectral response S(λ). The silicon response is a fixed parameter in our methodology because we assume the choice of sensor to be made a priori, based on the application. An example of such a characteristic is presented in Fig. 2.
A photon is characterized by a wavelength λ which is equivalent to an energy E. The inverse relationship between the energy of a photon and the wavelength is defined by:
where c is the speed of the light in vacuum and h the Planck’s constant.
Photons with an energy bellow 1.1e V will go through silicon without interaction [35]. This corresponds, according to Eq. 3, to a theoretical wavelength larger than 1125 n m. The structure of the sensor will, in fact, reduce this wavelength prior to reaching the theoretical limit [35]. Our model does not include the optical effects. Indeed, the layer of silicon can act as a FabryPérot interferometer [36], which modifies the absorption properties of the substrate. Without the optical effects, we can use directly the spectral response as defined in Fig. 2.
Filter considerations
The filters considered in SFAs are most often a realization of a FabryPérot interferometers today [5, 7, 37]. Future technical development may give rise to nanotube or nanohole SFA realizations [10]. Filters can be metallic filters [38, 39] or nano structure filters [40], and the substrate is often silicon or doped silicones, germanium, AsGa, indium, etc.
We use a simple Gaussian model to approximate FabryPérot interferometric filters, such as realized by SILIOS [37] and IMEC [8]. This approximation has been extensively done in the literature, in Fig. 3, we verify the impact of this approximation in term of energy inaccuracy before to push further our description. This figure shows the difference between a Gaussian approximation (G), a model of FabryPérot (FPM), and measurement of a practical FabryPérot designed by SILIOS (FPR). If we compute the differences between these curves as integral of function difference, we obtain the following numbers: d(F P R,G)=36.74%, d(F P M,F P R)=9.96% et d(F P M,G)=46.70%. Although the difference between G and the FabryPérot instances seems to indicate that the approximation may not be adequate, the difference between the FPM and the FPR indicates that depending on the technology, general models would not be adequate either due to noise and specific secret tuning of the manufacturer. We then make the statement that for illustration purpose it is more interesting to develop a methodology with Gaussian and to let the user of this methodology to tune it toward his own model of filters. In addition, if we only compute the difference in peak values, the Gaussian model is closer to the FabryPérot model than the measurement, due to the realization process. If we consider the differences between these curves within the standard deviation of the Gaussian, i.e. full width at half maximum (FWHM), the difference between FPM and FPG is greater than the difference with G and validate the use of the model by the assumption that the difference comes from the lobes of the FabryPérot filter, which generate only an offset spectrally consistent. This simple computation may also serve in the analysis of literature results obtained with Gaussian simulations.
The shape of the filter according to the Gaussian model depends on three parameters, namely the standard deviation σ, the scaling factor (the peak efficiency) B, and the average of the distribution (the location of the filter in the electromagnetic range) μ. Example of filters are given in Fig. 4 where we show a set of three Gaussian filters uniformly distributed. Within this model, we can write f _{ i } as in Eq. 4.
where μ _{ i } are fixed parameters according to the application and σ _{ i } and B _{ i } correspond to the given degree of freedom. The distribution of μ may be given a degree of freedom, but according to our knowledge, users usually optimize it for a specific application. Here, we consider that μ are provided in a priori design and as example, we consider that it is an equidistributed set or a list of selected peaks, such as shown in the Results section. In some cases, where μ are not provided or in a more global framework, it may be embedded into the optimization process with little extra work and restrained so that all channels do not end with the same transmittance.
As for example purpose, we define an equidistributed μ set, within the global applicable range of wavelengths Δ Λ=λ _{ max }−λ _{ min } and the number of channels N, so that the central positions μ _{ i } of each Gaussian are fixed before the optimization process. Consequently, in this case, the position μ _{ i } can be defined like in Eq. 5:
In a practical realization of an SFA system, we may consider two situations, as shown in Figs. 5(a) and (b).
We can consider two sensor configurations. One is very specific to some applications, such as in medical imaging [41, 42], where a very narrow spectral band of about 10 nm or ratio computed between specific narrow bands are required. In this case, Fig. 5(b) is a typical setup, and peaks cannot be given freedom. Another situation, where very wide bands are considered, does not to provide an optimal results either for demosaicing or for spectral reconstruction. After looking at the results of Sadeghipoor et al. [33] and Wang et al. [43], it seems that rather a relatively wide band is adequate for general demosaicing and spectral reconstruction. Indeed, we want to preserve separability between wavelengths if possible, e.g. NIR channels must not be polluted by visible information for computer vision, while maximizing spectral correlation for demosaicing. In this last case, we may authorize a little freedom on μ.
Singleparameter optimization
In this section, we derive the optimization process in the last section for a singleparameter optimization. We optimize the filters by modifying the parameters Bandσ of the Gaussian for each channel i, following Eq. 4. The aim of the single parameter optimization is to obtain the same camera response (i.e energy balance) for the N filters, while taking into account the substrate absorption and the spectral properties of the illumination. So we search for good settings of B _{ i } or σ _{ i } in the transmission characteristics of filters. We fix the sigma for all bands in the amplitude optimization, and we fix the amplitude in the standard deviation optimization. The response ρ _{ i } is calculated for a set of filter transmissions, f(B _{1},…,B _{ N }) and f(σ _{1},…,σ _{ N }), when searching respectively for the amplitudes and sigmas. In these two cases, optimal parameters are found such as \(\rho _{1}=\rho _{2}=\ldots =\rho _{N} = \tilde {C}\), a constant value.
For the analysis, we will consider the case of a camera based on N=[3,5,8,12] bands. We chose a threeband sensor to represent the typical color case, five bands are chosen for multispectral color sensors or colorimeters [44]. The eight bands are selected for comparison with an existing filter realization [3, 7]. Twelve bands are chosen, as they are expected to provide the best spectral reconstruction used in multispectral sensors [43, 45].
Energy balance by changing the magnitude
We change sequentially the magnitude of B _{ i } to obtain energetically balanced filters, whereas the value of σ is set as a constant for all bands during the optimization process.
We use a rescale algorithm. First, initial values of B _{ i,0} are set to 1. Equation 1 is applied to calculate initial response values ρ _{ i,0} for which channel responses are not balanced (uncorrected filters). Then, the optimal amplitude parameters are found by dividing the highest value among the set by ρ _{ i,0}. The equation could be written as this:
After normalization between [0,1], we finally find filter parameters B _{1},…,B _{ N }, for which each camera integrated responses ρ _{ i } share the same energy value \(\rho _{1}=\rho _{2}=\ldots =\rho _{N} = \tilde {C}\).
Figures 6(a), (b), (c) and (d) show the results after optimization under the illumination E for a perfect diffuser (O(λ)=1). In the range of wavelengths considered, σ is fixed and is equal to 48.0 for three filters, 28.8 for five filters, 18.0 for eight filters and 12.0 for twelve filters. For visualization and interpretation, we also show the response curve S(λ) of silicon contained in Fig. 2.
We observe that the magnitude of each filter depends directly on the response of silicon. These results highlight the fact that the silicon response is not uniform. In addition, if we consider a linear model of noise, the difference in amplitude is related to the difference of channel noise. The difference among filters increases as the bandwidth of the filters becomes narrow.
Energy balance by changing the standard deviation
This time, B is set as a constant for all bands. As for the amplitude optimization process, Eq. 1 is applied to calculate ρ _{ i,0} initial values for each filters. In order to obtain the Gaussian values of σ _{ i } (on a fixed interval) for which the filters share the same energy, we search for the minimum of the singlevariable function defined as this:
where \(C=\frac {\sum {\rho _{i,0}}}{N}\) and \(\tilde {C}\) is the integration result of camera response when changing sigma in Eq. 4.
The obtained energetically balanced filters are presented in Fig. 7(a), (b), (c) and (d), for a flat reflectance object and illuminant E. The amplitude B _{ i } is fixed and is equal to 1.
We observe that the standard deviation of each filter depends directly on the response of silicon. These results highlight the fact that the silicon response is not uniform. In addition, if we consider a linear model of noise, the difference in standard deviation is related to the difference of channel noise. The difference among filters increases as the bandwidth of the filters becomes narrow.
Limitations
When we optimize the magnitude or the standard deviation of the filters, we reach an exact solution in a few iterations, where standard deviation of the values {ρ _{1},ρ _{2},…,ρ _{ N }} is near or equal to zero. However, we observe an important reduction of efficiency in most of the optimal camera responses. Indeed, typically, 70% of the energy should be filtered out to guarantee the balance. This result is illustrated in Fig. 6(d). In Figs. 6(c) and (d), we can see that we severely reduced the transmittance efficiency of most of the filters. This appears to be an inefficient situation for practical implementation, considering the noise that might appear with low sensitivity filters. Such optimization is not adequate in practice. Spectral overlapping and spectral gap problems among bands are observed when we perform the optimization on σ (see Figs. 7(c) and (d)). Moreover, when illumination changes, filters change accordingly to the coupled pair {sensor, illuminant}. In a practical case, several illuminations may be considered, and an exact optimization may not give good results in several configurations.
These results suggest that a constrained multiparameter optimization may be a better fit to realize practical filter banks, which is the subject of the following section.
Multiparameter constrained optimization
In the last section, the optimization of parameters σ and B were proposed independently, such as presented by Péguillet et al. [29]. In this section, we present a new framework for finding optimal spectral filter responses according to the global energy balance of the camera system. To this aim, we define a set of constraints that should guarantee the following principles:

1.
There should be no gap in the spectral domain, otherwise some spectral information will be lost for some wavelengths.

2.
There should not be too much spectral overlapping of a filter.

3.
The sensor must globally be sensitive. Peak attenuation should be controlled.

4.
Energy balance must be achieved for an illumination set.
The framework tunes both parameters B and σ, in front of a set of illuminations given by the application, and constrained by the above principles. μ is presented as constant here, but may be embedded easily in the optimization parameters if the application allow freedom on the peak sensitivities. A block diagram of the method is presented in Fig. 8.
The silicon response S(λ) is a specific parameter as it corresponds to the selected imaging device. The illumination l _{ j }(λ) is also given as parameter. The global problem formula is as in Eq. 8:
given a set of linear inequalities that will be defined in the following sections. STD stands for the standard deviation of the set of computed sensor responses.
Constraints on amplitude
The value of the amplitude B must be as high as possible, otherwise the signal will be too weak compared to noise. In the same time, B _{ max } is not necessary 100% transmittance since a limitation may be imposed by the technology used to manufacture the filters. So, we take B _{ min } and B _{ max } in Eq. 9 as a minimum and maximum amplitude parameters. In any case, the maximum value can not be over 100%, and we notice that it may be beneficial for one band to be less efficient in order to be balanced with the others.
Constraints on overlapping
On Fig. 5(a) we observed too much overlapping. To avoid overlapping, according to the Fig. 9(b), we must respect the conditions of Eq. 10. We use the Gaussian percentiles to define this constraint, expressed as an integer weight of the standard deviation.
Equation 11 puts the Eq. 10 in the format A.x≤b in order to prepare a mathematical model for the optimization tool. A is a matrix and b is a vector. The constraint set of equation can be written as this:
where the values P and Q are user parameters. These are defined as 0<P<Q. The value of P can more likely be 1 and the value of Q, 3, according to the energy distribution of a Gaussian, presented on Fig. 9(a).
Constraints on gaps
On Fig. 5(b) we observed gaps. In the same manner as for the overlapping constraints, according to the Fig. 9(c), avoiding gap is done by respecting the conditions of Eq. 12. We can then define linear constraints as follow:
where the values P and Q are the same than the values used with the constraints to avoid the overlapping issue.
Resolution of the inverse problem
The optimization framework was implemented in Matlab using the fmincon routine to solve the function 8. The algorithm used is the activeset. In order to select the optimization solving method, we followed the MATLAB recommendation of use of fmincon function [46]. Indeed, the size of the problem is not very large (we can store all data during optimization, and linear algebra can be easily computed on them), the problem does not seem to show local minimum issues, although several approximated solutions may be possible. So this solution would provide the fastest optimization of parameters. In test, we noticed that other existing algorithms would reach the same results.
Results and discussion
General results
The code can be downloaded online along with different sensor responses to be tested. The set of constraints was defined as above. The range of wavelength chosen is the visible spectrum (380−780 n m), but the acquisition model described in the methods section can be generalized to NIR wavelengths, as long as the imaging model used is valid, i.e. without fluorescent effects and emission from objects. The set of constraints are linear functions. In practice, the multispectral device can be subject to a limitation of illuminant intensity, and to a limited range of exposure times, which is not addressed here. We also consider the sensor to have a linear response. If it was not the case, a linearization function may be used to correct the data. Figure 10 and Table 1 show a list of standard illuminants used. To perform simulations, we consider the case of a commercial sensor by SONY [47], and use the spectral sensor response presented in Fig. 2. A similar process can be used for most of the sensors, for which the absorption by wavelength is known or measurable.
Table 2 shows several iterations of the optimization process, for different illuminants and number of channels. By analyzing the results, we distinguish 3 specific behaviors:

1.
The optimization process converges to a feasible solution s t d(Π ^{′})≃0 (highlighted in blue in Table 2). See the model at the beggining of the document.

2.
The optimization process converges to an acceptable solution with a certain degree of freedom/spectral response deviation (highlighted in green in Table 2).

3.
The solution does not converge to an acceptable solution in term of standard deviation among the filter set (highlighted in red in Table 2).
These three cases are illustrated in Fig. 11. The first column shows the resulting set of filters for N=3,5,8 and 12, for which the optimization process converges successfully. So the integration for each area curve is quite constant along the wavelengths, and decreases with the number of channels N. We can clearly distinguish, for example, the shape contribution of the sensor sensitivity in Fig. 2 when using for example 12 filters and illuminant E.
The middle column is for the case where the optimization converges and where a minimum standard deviation among the filters is acceptable. Here, the D series of the illuminants, which have a globally homogeneous shape, are used. These are not very far from each other in term of emissivity, and the optimization succeeded.
The third column shows the case where the optimization process fails to recover proper filter shapes within this set of constraints. It is due to the fact that we use too many different viewing illuminants from which the shapes of emissivity are too far from each others. The transmission peaks remain at the minimal value (0.3) and the filters are not optimized since the standard deviation is more than 50% of the energy mean among the bands. Nevertheless, it is important to note that a good convergence could be found if constraints like B _{ min }, P or Q are adjusted and relaxed.
For comparison with a practical case, the response of the real camera design presented in Fig. 1 has a response mean of 3.8255 and a standard deviation of 1.71 through illuminant E. It typically corresponds to the case where the filters are not well optimized energetically.
Specific scenarios
To illustrate our method in practice, we propose to study three specific scenarios.
We study first the case of underwater imaging, where absorption of ocean water [48] is added in the model presented in the methods section. We use a visible range of wavelength (380−780 n m), one projector illuminant (tungsten), and 3 bands. A result of optimization with multiple parameters is shown in Fig. 12(a). We obtain an average score of 2.29 and a standard deviation of 6.74e ^{−6}. The optimization process succeeds in finding a minimum for the objective function in this case. The ocean water has a noticeable absorption in the red wavelength (between 700 and 780 n m), but we notice that the third "red" filter has a reduced transmission value than the first one. It could be explain by the fact that the tungsten illuminant has a emission even lower in the range of 380−500 n m, which explains why the optimized filters are compensated.
Another option is proposed on Fig. 12(b), where the peak sensitivities have been chosen according to the optimized sensitivities proposed by Parmar et al. [49]. We observe that the sensitivities proposed in this last case is more efficient. Index of quality between two filter banks, such as in Figs. 12(a) and (b), may be to choose the best sensitive option coupled with the best results in resulting color image.
A third case is a typical environment and camera system for colorimetric imaging purpose in a natural environment. This case uses 5 bands and performs well under all D illuminations. By doing the optimization process, we find an average of 3.07 and a standard deviation of 0.42. The corresponding filters are shown in Fig. 12(c). Here, the optimization process converges to an acceptable solution with a certain degree of freedom.
A fourth case example in Fig. 12(d) is the typical multispectral camera scenario for computer vision, where 8 bands in the visible and NIR (380−830 n m) are used. This scenario works with a relatively large variety of illuminants (A, F and D65). We obtain an average score of 0.82 with a standard deviation of 0.60, which is very high and not acceptable. Here, the optimization process fails to find a solution, due to the overconstrained environment. This is due to the selected set of illuminants that shows very different behaviors in term of spectral emissivity.
Also, for these last two cases, we observe low transmittance peak. Even if this order of efficiency magnitude is reasonable in practice according to the manufacturing status, this may be taken into account in releasing the constraints. Also, a larger standard deviation between the energies of the filters may be considered acceptable. It is though not too easy to decide on that before to incorporate noise influence into the simulations.
Discussion
The main contribution of this work is to show that the energy balance of SFA sensors is reachable in several illumination conditions by using an optimization process constrained by practical issues.
What we have not addressed here, is how this influences the actual imaging pipeline in term of noise reduction or in term of application accuracy. Indeed, we demonstrated here that we could manage to avoid saturated values while others would be very noisy low signal. In order to infer the effect on the complete imaging pipeline and fix/relax the constraints more accurately, one may either incorporate noise evaluation or incorporate this methodology into a more comprehensive optimization toward spectral reconstruction, demosaicing or color reproduction. Further work includes then the objective definition of the constraints.
Also, the optical model that we use is rather limited in practice since filter transmittance are not Gaussian. In order to demonstrate our purpose, it was not critical since the FWHM part of the curves are rather similar to real filters, however, when it comes to noise and correlation between channels, it makes no doubt that a more accurate optical transmittance model must be used. We also did not address multilobe sensitivities, indeed it could be beneficial for increasing the sensitivity of the sensor, however, the spectral community tends to prefer singlelobes in order to separate better between wavelengths. Such concept would be nevertheless beneficial for color imaging based on less than 4 bands. However, with the increasing of bands, even color will not gain much from multilobe curves.
Our main message is that this problem of energy balance, which could be reasonably ignored in the case of CFA or SFA cameras based on large interference filters, may not be ignored in the SFA case.
Conclusions
We present a framework that permits filter optimization that achieves energy balance in SFA device configurations. This method can be useful for a great range of applications, where multiple captures, noise or dynamic range are crucial aspects in the final response of the camera. The implementation permits an investigation on how constraints imposed on the system can affect the efficiency and the reliability of the final product under several illumination and capture conditions. Various simulations have been shown to test robustness of the proposed method. Further work envisions the incorporation of this methodology into the application dependent or general purpose filter optimization, its implication into the imaging pipeline and practical realizations of filter and practical evaluation of SFA cameras.
Abbreviations
 CFA:

Color filter array
 FPM:

Fabry Pérot model
 FPR:

Fabry Pérot real
 FWHM:

Full width at half maximum
 MSI:

MultiSpectral imaging
 NIR:

Near infrared
 SFA:

Spectral filter array
References
 1
Hagen, N, Kudenov, MW: Review of snapshot spectral imaging technologies. Opt Eng. 52(9), 090901–090901 (2013).
 2
Vagni, F: Survey of hyperspectral and multispectral imaging technologies (etude sur les technologies d’imagerie hyperspectrale et multispectrale). Technical report, DTIC Document (2007). http://oai.dtic.mil/oai/oai?verb=getRecord&metadataPrefix=html&identifier=ADA473675.
 3
Lapray, PJ, Wang, X, Thomas, JB, Gouton, P: Multispectral filter arrays: Recent advances and practical implementation. Sensors. 14(11), 21626–21659 (2014). doi:10.3390/s141121626.
 4
Bayer, BE: Color imaging array. Google Patents. US Patent 3, 971,065 (1976). http://www.google.com/patents/US3971065. Accessed 6 Jan 2017.
 5
Geelen, B, Tack, N, Lambrechts, A: A snapshot multispectral imager with integrated tiled filters and optical duplication. In: Advanced Fabrication Technologies for Micro/Nano Optics and Photonics VI. Proc. SPIE, p. 861314 (2013), doi:10.1117/12.2004072. http://dx.doi.org/10.1117/12.2004072.
 6
Monno, Y, Kikuchi, S, Tanaka, M, Okutomi, M: A practical oneshot multispectral imaging system using a single image sensor. IEEE Trans Image Process. 24(10), 3048–3059 (2015). doi:10.1109/TIP.2015.2436342.
 7
Thomas, JB, Lapray, PJ, Gouton, P, Clerc, C: Spectral characterization of a prototype sfa camera for joint visible and nir acquisition. Sensors. 16(7), 993 (2016). doi:10.3390/s16070993.
 8
Geelen, B, Tack, N, Lambrechts, A: A compact snapshot multispectral imager with a monolithically integrated perpixel filter mosaic. In: Advanced Fabrication Technologies for Micro/Nano Optics and Photonics VII. Proc. SPIE, p. 89740 (2014), doi:10.1117/12.2037607. http://dx.doi.org/10.1117/12.2037607.
 9
Walls, K, Chen, Q, Grant, J, Collins, S, Cumming, DRS, Drysdale, TD: Narrowband multispectral filter set for visible band. Opt. Express. 20(20), 21917–21923 (2012). doi:10.1364/OE.20.021917.
 10
Park, H, Crozier, KB: Multispectral imaging with vertical silicon nanowires. Sci. Rep. 3 (2013).
 11
Frey, L, Masarotto, L, Armand, M, Charles, ML, Lartigue, O: Multispectral interference filter arrays with compensation of angular dependence or extended spectral range. Opt. express. 23(9), 11799–11812 (2015).
 12
Miao, L, Qi, H, Ramanath, R, Snyder, WE: Binary treebased generic demosaicking algorithm for multispectral filter arrays. IEEE Trans. Image Process. 15(11), 3550–3558 (2006). doi:10.1109/TIP.2006.877476.
 13
Wang, X, Thomas, JB, Hardeberg, JY, Gouton, P: Discrete wavelet transform based multispectral filter array demosaicking. In: Colour and Visual Computing Symposium (CVCS), 2013, pp. 1–6 (2013), doi:10.1109/CVCS.2013.6626274.
 14
Wang, X, Thomas, JB, Hardeberg, JY, Gouton, P: Median Filtering in Multispectral Filter Array Demosaicking (2013). doi:10.1117/12.2005256. 10.1117/12.2005256.
 15
Monno, Y, Kiku, D, Kikuchi, S, Tanaka, M, Okutomi, M: Multispectral demosaicking with novel guide image generation and residual interpolation. In: 2014 IEEE International Conference on Image Processing, ICIP 2014, Paris, France, October 2730, 2014, pp. 645–649 (2014), doi:10.1109/ICIP.2014.7025129. http://dx.doi.org/10.1109/ICIP.2014.7025129.
 16
Aggarwal, HK, Majumdar, A: Compressive sensing multispectral demosaicing from single sensor architecture. In: Signal and Information Processing (ChinaSIP), 2014 IEEE China Summit & International Conference on. IEEE, pp. 334–338 (2014).
 17
Mihoubi, S, Losson, O, Mathon, B, Macaire, L: Multispectral demosaicing using intensitybased spectral correlation. In: Image Processing Theory, Tools and Applications (IPTA), 2015 International Conference On, pp. 461–466 (2015), doi:10.1109/IPTA.2015.7367188.
 18
Wang, X, Pedersena, M, Thomas, JB: The influence of chromatic aberration on demosaicking. In: Visual Information Processing (EUVIP), 2014 5th European Workshop On, pp. 1–6 (2014), doi:10.1109/EUVIP.2014.7018410.
 19
Wang, X, Green, PJ, Thomas, J, Hardeberg, JY, Gouton, P: Evaluation of the colorimetric performance of singlesensor image acquisition systems employing colour and multispectral filter array. In: Trémeau, A, Schettini, R, Tominaga, S (eds.)Computational Color Imaging  5th International Workshop, CCIW 2015, Saint Etienne, France, March 2426, 2015, Proceedings. Lecture Notes in Computer Science, pp. 181–191. Springer (2015), doi:10.1007/9783319159799. http://dx.doi.org/10.1007/9783319159799.
 20
Holst, GC, Lomheim, TS: CMOS/CCD Sensors and Camera Systems, SPIEInternational Society for Optical Engine (2007). http://spie.org/Publications/Book/740657.
 21
Theuwissen, AJ: Solidstate Imaging with Chargecoupled Devices, Vol. 1. Springer (2006). doi:10.1007/0306471191. http://link.springer.com/book/10.1007\%2F0306471191.
 22
Fellers, TJ, Davidson, MW: Concepts in Digital Imaging Technology–CCD Saturation and Blooming (2005). http://micro.magnet.fsu.edu/primer/digitalimaging/concepts/ccdsatandblooming.html. Accessed 6 Jan 2017.
 23
Clark, RN: Digital camera sensor performance summary (2008). Available at: (http://www.clarkvision.com/articles/digital.sensor.performance.summary/). Accessed 15 Oct 2010.
 24
Richards, WA: Lightness scale from image intensity distributions. Applied Optics. 21(14), 2569–2582 (1982).
 25
Johnson, BK: Photographic exposure control system and method. Google Patents. US Patent 4,423,936 (1984). https://www.google.com/patents/US4423936.
 26
Ieong, M, Guarini, KW, Chan, V, Bernstein, K, Joshi, R, Kedzierski, J, Haensch, W: Three dimensional cmos devices and integrated circuits. In: Custom Integrated Circuits Conference, 2003. Proceedings of the IEEE 2003, pp. 207–213. IEEE (2003). doi:10.1109/CICC.2003.1249391. http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=1249391&isnumber=27953.
 27
Lapray, PJ, Heyrman, B, Ginhac, D: Hardwarebased smart camera for recovering high dynamic range video from multiple exposures. Optical Engineering. 53(10), 102110 (2014).doi:10.1117/1.OE.53.10.102110.
 28
Lapray, PJ, Heyrman, B, Ginhac, D: Hdrartist: an adaptive realtime smart camera for high dynamic range imaging. J RealTime Image Process. 12, 1–16 (2014).
 29
Péguillet, H, Thomas, JB, Ruicheck, Y, Gouton, P: Energy Balance in Single Exposure Multispectral Sensors. In: Colour and Visual Computing Symposium (CVCS). IEEE (2013). doi:10.1109/CVCS.2013.6626277.
 30
Vrhel, MJ, Trussell, HJ: Filter considerations in color correction. Image Process. IEEE Trans. 3(2), 147–161 (1994).
 31
Vhrel, M, Trussell, H, Bosch, J: Design and realization of optimal color filters for multiilluminant correction. J. Electronic Imaging. 4(6), 147–161 (1995).
 32
Sharma, G, Trussell, HJ: Optimal filter design for multiilluminant color correction. In: Proc. IS&T/OSA’s Optics and Imaging in the Information Age, pp. 83–86 (1996).
 33
Sadeghipoor Kermani, Z, Lu, Y, Sússtrunk, S: Optimum Spectral Sensitivity Functions for Single Sensor Color Imaging. In: Proceedings of IS&T/SPIE EI: Digital Photography VIII. Proceedings of SPIE, p. 829904. SpieInt Soc Optical Engineering, Po Box 10, Bellingham, Wa 982270010 Usa (2012), doi:10.1117/12.907904.
 34
Monno, Y, Kitao, T, Tanaka, M, Okutomi, M: Optimal spectral sensitivity functions for a singlecamera oneshot multispectral imaging system. In: 2012 19th IEEE International Conference on Image Processing, pp. 2137–2140 (2012), doi:10.1109/ICIP.2012.6467315.
 35
Darmont, A: Spectral response of silicon image sensors. Technical report, Aphesa, April 2009. White paper.
 36
Wenjie, L, Bockrath, M, Bozovic, D, Hafner, JH, Tinkham, M, Hongkun, P: FabryPerot interference in a nanotube electron waveguide. Nature. 1 (2001).
 37
TECHNOLOGIES, S: MICROOPTICS Supplier. http://www.silios.com/. Accessed 1 July 2016.
 38
Wang, Y: Lowcost thinmetalfilm interference filters. Google Patents. US Patent 6,031,653 (2000). https://www.google.com/patents/US6031653. Accessed 6 Jan 2017.
 39
AJ, H, GR, M: Infrared interference filter. Google Patents. US Patent 3,682,528 (1972). https://www.google.ch/patents/US3682528. Accessed 6 Jan 2017.
 40
NguyenHuu, N, Lo, YL, Chen, YB, Yang, TY: Subwavelength metallic gratings as an integrated device: polarized color filter. In: Proceedings of SPIE  The International Society for Optical Engineering, pp. 79340–79340. Proc. SPIE 7934, Optical Components and Materials VIII, 79340U (2011). doi:10.1117/12.871600.
 41
Kong, L, Yi, D, Sprigle, S, Wang, F, Wang, C, Liu, F, Adibi, A, Tummala, R: Single sensor that outputs narrowband multispectral images. J. Biomed. Optics. 15(1), 010502–0105023 (2010). doi:10.1117/1.3277669.
 42
Flux Data Incorporated, Availableonline. http://www.fluxdata.com. Accessed on 1 July 2016.
 43
Wang, X, Thomas, JB, Hardeberg, JY, Gouton, P: Multispectral imaging: narrow or wide band filters?. JAICJ. Int. Colour Assoc. 12 (2014).
 44
Monno, Y, Tanaka, M, Okutomi, M: Multispectral demosaicking using adaptive kernel upsampling. In: 2011 18th IEEE International Conference on Image Processing, pp. 3157–3160 (2011), doi:10.1109/ICIP.2011.6116337.
 45
Wang, X, Thomas, JB, Hardeberg, JY, Gouton, P: A Study on the Impact of Spectral Characteristics of Filters on Multispectral Image Acquisition. In: Lindsay MacDonald, SW, Stephen Westland (eds.)12th Congress of the International Colour Association, pp. 1765–1768, Gateshead, United Kingdom (2013). https://hal.archivesouvertes.fr/hal00847828. Accessed 6 Jan 2017.
 46
Choosing the Algorithm. Available online: https://fr.mathworks.com/help/optim/ug/choosingthealgorithm.html. Accessed 1 Oct 2016.
 47
Sony Pregius Global Shutter CMOS Imaging Performance. Available online: http://www.ptgrey.com/support/downloads/10414. Accessed 1 July 2016.
 48
Shifrin, K: Physical Optics of Ocean Water. Springer (1988).
 49
Parmar, M, Reeves, SJ: Selection of optimal spectral sensitivity functions for color filter arrays. In: 2006 International Conference on Image Processing, pp. 1005–1008 (2006), doi:10.1109/ICIP.2006.312669.
 50
Protassov, K: Analyse Statistique de Données Expérimentales. In: SOFEDIS, p. 29. EDP Sciences ISBN 2868835902 (2012).
Funding
The research project is funded by the EU and the H2020 EXIST project No. 662222.
Availability of data and materials
The datasets used during the current study are available in the Github repository, https://github.com/pjlapray/Energy_balance_in_SFA following this link. A password to extract the compressed content is available from the corresponding author on reasonable request.
Authors’ contributions
Concept and methodology: JBT (35%) PJL (35%) PG (15%) YR (15%) Mathematics: JBT (50%) PJL (40%) PG (5%) YR (5%) Implementation: JBT (30%) PJL (70%) Writing: JBT (30%) PJL (50%) PG (10%) YR (10%). All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Author information
Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Lapray, PJ., Thomas, JB., Gouton, P. et al. Energy balance in Spectral Filter Array camera design. J. Eur. Opt. Soc.Rapid Publ. 13, 1 (2017). https://doi.org/10.1186/s4147601600317
Received:
Accepted:
Published:
Keywords
 Spectral filter arrays
 Multispectral imaging
 Optical design