Per-Channel Regularization for Regression-Based Spectral Reconstruction

Abstract

Spectral reconstruction algorithms seek to recover spectra from RGB images. This estimation problem is often formulated as leastsquares regression, and a Tikhonov regularization is generally incorporated, both to support stable estimation in the presence of noise and to prevent over- tting. The degree of regularization is controlled by a single penalty-term parameter, which is often selected using the cross validation experimental methodology. In this paper, we generalize the simple regularization approach to admit a per-spectral-channel optimization setting, and a modi ed cross-validation procedure is developed. Experiments validate our method. Compared to the conventional regularization, our perchannel approach signi cantly improves the reconstruction accuracy at multiple spectral channels, by up to 17% increments for all the considered models.

Keywords

• Spectral reconstruction spectral imaging

Introduction

The light spectrum is a continuous intensity distribution across wavelengths. This spectral information is commonly used to help determine and/or discriminate the physical properties of object surfaces, for example in remote sensing [ 25, 8, 14, 23 ] and medical imaging [ 28, 29 ]. Also, in various practical applications, the devices (sensors or displays), light sources and object surfaces are characterized by spectral measurements [ 1, 9, 16, 27 ].

Despite the advantages of spectral capture, almost all images that we record contain just 3 measurements - the 3 weighted spectral averages over the Red, Green and Blue spectral regions. Perforce, much spectral information is therefore lost in this RGB image formation process. Indeed, it is a classical result in color science, that there are many spectra - called metamers [ 12 ] - which integrate to the same RGB, and of course, given only one RGB measurement we cannot tell which physical spectrum induced it. Still, by adopting learning approaches we can estimate the spectrum that likely corresponds to a given RGB. Estimating spectra from RGBs is called spectral reconstruction (SR).

In Fig. 1 we illustrate RGB image formation and the SR process. In the top-left panel we see a single radiance spectrum measured at one location in the hyperspectral image (bottom left). This spectrum is sampled by 3 sensors, resulting in the 3-value RGB response (top-right). Repeating this process for all image locations, the corresponding RGB image in the bottom right is derived from the hyperspectral image. Then, the spectral reconstruction algorithms attempt to recover the hyperspectral image back from the RGB image (or an approximation thereof).

Historically, this SR problem is e ectively solved by least-squares regression [ 15, 11, 17, 18, 2 ], where the map from RGBs (or the non-linear RGB features) to spectra is modelled as a simple linear transform. More recently, deep learning approaches [ 4, 22, 6 ] have been developed that provide even better SR performance. E ectively, this performance increment is achieved by regressing an RGB in the context of its neighborhood to its corresponding spectra. Clearly, this patch-based idea has merit. For example, if the algorithm can identify a patch in the scene as a `skin region' then spectral recovery is plausibly easier to solve, since skin spectra have characteristic spectral shapes [ 7, 19 ].

Despite clear rationale behind the deep-learning approach, Aeschbacher et al. [ 2 ] show that the regression-based A+ algorithm provides very competitive performance. Moreover, Lin and Finlayson [ 17 ] shows that several regression methods actually generalize better than the leading deep-learning models when the scene exposure changes.

The main concern of this paper is in the regularization step of regressionbased SR algorithms. The classical (multivariate) regression problem from statistics is written as

MA

B ; where A is an m N matrix as the table of measured data (m is the dimension of the measured data and N is the number of data samples, with N m), and B is the corresponding target data matrix, of dimension k N (k is the dimension of the target data). The aim is to nd the k m linear mapping M that makes the approximation as good as possible.

Now, let us suppose small uctuations in the target data, denoted as a matrix E of very small numbers (all entries are close to 0). The following regression is almost identical as Equation (1): (1) M0A

B + E : (2) And yet, often we nd that the best solution for M and M0 are very di erent from one another. The reason for this is that some dimensions of the measured data (the rows of the data matrix A) could be highly correlated, such that there can be very di erent M's that t B equally well.

Regularization theory [ 24 ] is a way of dealing with this kind of non-robustness. Given the example of Equation (1) and (2), we may ask that among all plausible (near optimized) solutions, which one is more likely to generalize to the `unseen data' better. Typically, the principle of regularization follows the idea that the best tting function (i.e. M) should be the simplest possible solution that can still ts the data well.

In Fig. 2 we show a 1-D example. In the least-squares sense, the wiggly red curve is found to best t the training data (black data points). Yet intuitively, this is not the correct t, as the data points appear to follow a much simpler distribution. In contrast, the regularized t (blue curve) seems to model the data better.

Returning to our spectral reconstruction problem, taking linear regression as an example, the matrix A corresponds to a set of image RGBs (m = 3), and B refers to the spectra we are trying to recover (k is the number of spectral channels we have measured). The insight that we explore in this paper is that there is no reason why the ts for di erent spectral channels should be regularized altogether. Rather, we seek to consider the per-channel regression problem, where each spectral channel is recovered in turn, and correspondingly, each row of M is solved in turn. This simple modi cation allows us to carry out a per-channel regularization that ensures individual optimizations for all spectral channels in the spectral reconstruction problem.

Either for the conventional global regularization or for our per-channel approach, care must be taken not to overly tune the terms in M to the data at hand. This led us to develop a modi ed cross validation procedure. Our method separates the data on hand into three subsets, respectively for training, regularization and testing, which is a novel adjustment from the standard methodology [ 20 ] and is another contribution of this paper. 2 2.1

Background

Hyperspectral Images and RGB Simulation

In a hyperspectral image, spectra are measured discretely at some sampled wavelengths. Suppose the visible spectrum runs from 400 to 700 nanometers and the spectral sampling is every 10 nanometers, we get a 31-dimensional discrete representation of spectra, denoted as r 2 R31.

Correspondingly, the spectral sensitivities of the R, G and B camera sensors can also be represented in discrete vector form (i.e. as 31-dimensional vectors), respectively denoted as sR, sG and sB. Then, as per the illustration in Fig. 1 we can write image formation as [ 26 ]: 0R1 0sT 1

R x = BGC = B@sTGCA r ; @ A

B sTB where x = (R; G; B)T is the 3-value RGB camera response.

In the SR problem (the bottom of Fig. 1) we seek to recover hyperspectral images from the RGB images. Denote an SR algorithm as : R3 7! R31, 2.2

Regression-Based Spectral Reconstruction

The general regression-based formulation of the SR problem is written as (x)

r : (x) = M'(x) ; where '( ) is a feature transformation which maps each RGB to a corresponding p-term feature vector, and in turn it is mapped by the regression matrix (3) (4) (5) M. Each of the various regression-based models [ 15, 11, 17, 18, 2 ] adopts a bespoke de nition of '(x). For details of the considered models, including Linear, Root-Polynomial and Adjusted Anchored Neighborhood Regression [ 15, 17, 2 ], see Appendix A.

Least-Squares Optimization The most common least-squares optimization seeks to minimize the sum of squared errors between the ground-truth training spectral data and the reconstruction from their corresponding RGBs: (x). Given the formulation of (x) in Equation (5), the least-squares optimization of M is formulated as:

M = arg min

M

N X i=1 jjri

M'(xi)jj22 ; where N is the number of data points in the training set and i indexes an individual spectrum.

Collating all spectral training data in a data matrix R = (r1; r2; :::; rN ) and the corresponding feature vector matrix = ('(x1); '(x2); :::; '(xN )), Equation (6) can be written as:

M = arg min

M jjR

M jj2F : Here jj jj2F is the squared Frobenius norm, which is exactly the sum-of-squares of all entries of the enclosed matrix.

Tikhonov Regularization In regression-based SR, the most common method to regularize a model is Tikhonov Regularization [ 15, 24 ], which hypothesizes that a more natural t is obtained when the `magnitude' (or the `matrix norm') of M is bounded to some extent. Based on this assumption, the least-squares optimization in Equation (7) is extended to incorporate a regularization term: M = arg min

M jjR

M jj2F + jjMjj2F : Here, the jjMjj2F term (the regularization term, or penalty term) is controlled by a user-de ned regularization parameter 0, which is usually determined empirically [ 13, 21 ].

Equation (8) is solved in closed form [ 15 ]:

M = R

T(

T + I) 1 ; where I is the p vectors '(x)).

p identity matrix (recall that p is the dimension of the feature (6) (7) (8) (9)

Proposed Method

In spectral reconstruction, we wish to recover spectral measurements in the range from 400 to 700 nanometers (the visible spectrum). Suppose we know the intensity of light entering the camera at 400 nanometers. Given this knowledge if we wished to predict the value of the spectrum at 410 nanometers, it makes sense to assume a similar value as the one at 400 nanometers. Indeed, the fact that the intensity values at close-by wavelengths are similar is why we can represent the continuous visible spectrum at discrete wavelengths. Conversely, one could not use the knowledge of light at 400 nanometers to predict the spectral value at, say, 700 nanometers.

And yet, in the literature when we regularize the regression-based SR models, we are - in some sense - assuming that all wavelengths are related. Our new per-channel reformulation of Tikhonov regularization for spectral reconstruction e ectively allows the recovery of the spectral values at distant wavelengths to be considered more independently from one another. 3.1

Per-Channel Regularization

Let us split the regression matrix M by row: M = (m1; m2; :::; m31)T, such that the general form of regression-based SR formulated in Equation (5) becomes (10) (11) 0 m1T 1

0 r^1 1 (x) = M'(x) = BBBB m...2T CCCC '(x) = BBBB r^...2 CCCC ; where (r^1; r^2; :::; r^31)T = ^r is the reconstructed spectrum. For an arbitrary kth spectral channel, the estimated intensity value r^k is given by r^k = mkT'(x).

Note that as we represent the regression model by channel, we do not alter the original model. This says that the regression-based spectral reconstruction has always been in such a way that the reconstruction for each spectral channel depends exclusively on the corresponding row of M.

Given this fact, we might expect that each row of M would be optimized independently. However, this was not the case for the conventional regularized least-squares solution in Equation (9). Indeed, we see in Equation (8) the regularization term is only controlled by one single regularization parameter and the ts for all spectral channels (all rows of M) are regularized by the same . Regardless of how we optimize this parameter, this setting clearly contradicts to the inherent independence between the rows of M.

Let us split the spectral reconstruction problem into 31 independent problems, where the function k : R3 7! R reconstructs the kth-channel values of the reconstructed spectra by the kth row of M: k(x) = mkT'(x) : where kT includes the kth channel values of all training spectral data. mkT is then optimized following the regularized least-squares optimization: Then, we are to determine mkT, as the least-squares t for the kth channel data. Recall the training spectral data matrix R whose columns are individual training spectra, now we split R by spectral channel instead: (12) (13) (14) mkT = arg min mTk

T jj k mkT jj22 + kjjmkjj2 2 ; and solved in closed form: mkT = kT T(

T + kI) 1 : Here k represents the channel-wise regularization parameter that only controls the regularization for the kth channel.

Clearly, our per-channel regularization scheme (Equation (14)) solves the regression matrix M row-by-row, such that each row is ready to be regularized independently. The remaining question is then how we are going to optimize these regularization parameters. 3.2

Modi ed Cross Validation

Perforce, regularization parameters (conventionally the single and our perchannel k's) are determined empirically. In the literature, a grid-search approach is adopted, where di erent parameters are tried to regularize the model. These `intermediate models' are then used to recover spectra from a set of unseen RGB images, and the model that minimizes the evaluation criteria is selected. For example, see [ 11, 17, 2 ].

As usual we would like to train, regularize and test a model using the images from the same database, we must partition the database into several subsets for these di erent usages. All (to our knowledge) deep-learning models simply separate the image database into 3 subsets randomly (respectively for training, validation and testing, in the parlance of deep learning), see [ 22, 6 ]. However, this setting can potentially create so-called `unfair' separations, such that if the database is separated di erently, the results may vary.

A better practice is using a cross-validation process. In this paper we develop our own cross validation scheme, which is modi ed from the conventional Kfold cross validation [ 20 ]. This is because the conventional K-fold only seeks to separate a dataset into training and testing data, and here we need an additional partition as the regularization data.

In Fig. 3 we show the comparison between the conventional 4-fold cross validation (left) and our method (right). For both methods the same experiment is conducted 4 times. In each trial, the conventional method selects 3 out of 4 portions of data for training (marked in blue) and the remaining part is for testing (marked in orange). In our method, however, only 2 out of 4 portions of data are for training, which allows 1 portion of data (marked in green) used for regularization, that is to determine the and k parameters. Subject to these terms we solve for the best regression model for the training (blue) data. The performance statistics are calculated based on the recovery errors on the testing (orange) data and averaged over the 4 trials.

Notice for our cross validation method there actually exists more possible permutations than the presented 4-trial setting. To be exact, there should be 12 di erent permutations. We remark that according to our empirical study, experimenting with more trials (than the presented setting) does not make signi cant di erence in the performance statistics. 4 4.1

Experiments

Considered Models

In this paper we consider 3 regression-based models: { Linear Regression (LR) [ 15 ] { Root-Polynomial Regression (RPR) [ 17 ] { Adjusted Anchored Neighborhood Regression (A+) [ 2 ] .

For all the above models we adopt both the original regularization (as reported in their citations and as per Equation (9)) and our per-channel regularization (as per Equation (14)).

Database

We use the ICVL hyperspectral image database [ 3 ] (Fig. 4), which provides 201 hyperspectral images of spatial dimension 1392 1300 and with 31 spectral channels. The spectral channels represent narrow-band intensity measurements, respectively at every 10 nanometers (nm) between 400 and 700 nm.

The corresponding RGB images are simulated following the linear RGB simulation setting (Equation (3)), and the CIE 1964 color matching functions [ 10 ] are used as the spectral sensitivities. The selected evaluation metric is Relative Absolute Error (RAE) [ 4, 5 ], which is de ned per channel as:

RAE(rk; r^k) = rk rk r^k ; (15) where rk and r^k are respectively the kth-channel values of the ground-truth and reconstructed spectra. E ectively, this metric measures the percentage absolute error. RAE is the most common performance measure used in recent research, and the rationale of using this metric can be found in [ 4 ]. 5

Results and Discussion

In Table 1, we present the per-channel error statistics of LR (left table), RPR (middle table) and A+ (right table) under the original settings - where a single penalty term is used in the regularization [ 15, 17, 2 ] - and our per-channel regularization method. The spectral channels are represented by the wavelengths (nm). We also calculate the percentage `improvements' as:

Improve (%) = 100

RAEoriginal

RAEours ;

RAEoriginal

(16) which is presented in the rightmost column of each table. In the bottom of each table, the Mean RAE results (averaging over all spectral channels) are shown.

First, we see that for all considered models, our method improves the RAE in multiple channels by over 10% (marked in bold and with underlines), with maximal improvements around 16-17%.

Secondly, in terms of Mean RAE performance, our method improves the RPR model the most, by an 8.6% increment, compared to 3.2% for A+ and 3.1% for LR. Signi cantly, the A+ model is the leading sparse coding model, which is shown able to perform equally-well as some deep-learning solutions [ 2 ]. By improving the A+ model, we e ectively bring forward the shallow-learned baseline. Moreover, our method reduces the gap between RPR and A+. Relative to A+, RPR model is much simpler (with signi cantly less model parameters) [ 17 ], which ensures more e ective model re-training and shorter runtime.

Lastly, for the A+ model, it seems curious that the per-channel performances in the rst three channels (400, 410 and 420 nm) degrade by minute di erences. Indeed, this means the regularization parameters we chose for these channels are not actually optimized for the test-set data. We remark that this is most likely originated from the unequal separation of data subsets in cross validation, such that the best regularization parameter for the regularization-set data does not correspond to the best for the test-set data. We are investigating how to remedy this issue.

For one example image in the ICVL database [ 3 ], we visualize the spectral reconstruction errors as the Mean RAE error maps in Fig. 5. It is clear that for all models our method improves the Mean RAE in various parts of the image. For example, the tree stem for LR and RPR, and the grassy ground for A+. In the spectral reconstruction (SR) problem, hyperspectral images are reconstructed from RGB images. Many approaches are based on least-squares regression, where the tting function is modelled by a simple linear transformation, and a Tikhonov regularization process is applied to improve model generalizability. Conventionally, the ts for all spectral channels are jointly regularized. We demonstrate that the t for each spectral channel can be formulated independently, such that the t for each channel is regularized (therefore optimized) independently. We also provide a novel modi cation of K-fold cross validation so that the models can be fairly trained, regularized and tested with a single image database. Compared to the original models, our per-channel regularization method improves the accuracy of recovery for individual spectral channels by up to 17% increments, and by 3-9% in mean improvements over all spectral channels.

Appendix: Regression Models

A.1 Linear regression (LR) [ 15 ] assumes a linear map from RGB to spectra. The spectral estimate is written as (x) = Mx ; (17)

3 regression matrix. As a simple non-linear extension from LR, in Root-Polynomial Regression (RPR) [ 17 ] we expand a series of root-polynomial terms from each RGBs. Denote ' : R3 7! Rp as the -order root-polynomial transformation, the example transformations for the 2nd; 3rd and 4th order RPR are:

Order Root-Polynomials '2(x) R; G; B; pRG; pGB; pBR

R; G; B; pRG; pGB; pBR; '3(x) 3 RG2; p3GB2; p3BR2; 3 R2G; p3G2B; p3B2R; p3RGB p p R; G; B; pRG; pGB; pBR; p 3 RG2; p3GB2; p3BR2; p3R2G; p3G2B; p3B2R; p3RGB; '4(x) p4R3G; 4 R3B; 4 G3R; 4 G3B; 4 B3R; 4 B3G;

p p p p p p p p 4 R2GB; 4 G2RB; 4 B2RG The spectral reconstruction then seeks to linearly map these root-polynomial vectors to spectra: (x) = M' (x) : (18) In this paper we set

= 6, which is the 6th-order RPR.

A.3

Adjusted Anchored Neighborhood Regression (A+ Sparse

Coding)

The leading sparse coding method `A+' [ 2 ] assumes linear maps from RGBs to spectra (e ectively, operates LR in every neighborhood). Denote i(x) as the spectral reconstruction mapping for the data in the ith neighborhood. On input of an RGB x, the reconstruction is written as: i(x) = Mix : (19) See [ 2 ] for more details about the model.

References

1. Abrardo , A. , Alparone , L. , Cappellini , I. , Prosperi , A. : Color constancy from multispectral images . In: Proceedings of the International Conference on Image Processing . vol. 3 , pp. 570 { 574 . IEEE ( 1999 )

2. Aeschbacher , J. , Wu , J. , Timofte , R.: In defense of shallow learned spectral reconstruction from RGB images . In: Proceedings of the International Conference on Computer Vision . pp. 471 { 479 . IEEE ( 2017 )

3. Arad , B. , Ben-Shahar , O. : Sparse recovery of hyperspectral signal from natural RGB images . In: Proceedings of the European Conference on Computer Vision . pp. 19 { 34 . Springer ( 2016 )

4. Arad , B. , Ben-Shahar , O. , Timofte , R. , et al.: NTIRE 2018 challenge on spectral reconstruction from RGB images . In: Proceedings of the Conference on Computer Vision and Pattern Recognition Workshops . pp. 929 { 938 . IEEE ( 2018 )

5. Arad , B. , Timofte , R. , Ben-Shahar , O. , Lin , Y. , Finlayson , G.D. , et al.: NTIRE 2020 challenge on spectral reconstruction from an RGB image . In: Proceedings of the Conference on Computer Vision and Pattern Recognition Workshops . IEEE ( 2020 )

6. Arun , P.V. , Buddhiraju , K.M. , Porwal , A. , Chanussot , J.: CNN based spectral super-resolution of remote sensing images . Signal Processing 169 , 107394 ( 2020 )

7. Bashkatov , A.N. , Genina , E.A. , Kochubey , V.I. , Tuchin , V.V. : Optical properties of the subcutaneous adipose tissue in the spectral range 400{2500 nm . Optics and Spectroscopy 99 ( 5 ), 836 { 842 ( 2005 )

8. Chen , Y. , Zhao , X. , Jia , X. : Spectral{spatial classi cation of hyperspectral data based on deep belief network . IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing 8 ( 6 ), 2381 { 2392 ( 2015 )

9. Cheung , V. , Westland , S. , Li , C. , Hardeberg , J.Y. , Connah , D.R. : Characterization of trichromatic color cameras by using a new multispectral imaging technique . Journal of the Optical Society of America A 22 ( 7 ), 1231 { 1240 ( 2005 )

10. Commission Internationale de l' Eclairage: CIE proceedings 1964 Vienna session, committee report E-1.4. 1 ( 1964 )

11. Connah , D.R. , Hardeberg , J.Y.: Spectral recovery using polynomial models . In: Color Imaging X: Processing, Hardcopy, and Applications . vol. 5667 , pp. 65 { 75 . International Society for Optics and Photonics ( 2005 )

12. Finlayson , G.D. , Morovic , P. : Metamer sets . Journal of the Optical Society of America A 22 ( 5 ), 810 { 819 ( 2005 )

13. Galatsanos , N.P. , Katsaggelos , A.K. : Methods for choosing the regularization parameter and estimating the noise variance in image restoration and their relation . IEEE Transactions on Image Processing 1 ( 3 ), 322 { 336 ( 1992 )

14. Ghamisi , P. , Dalla Mura , M. , Benediktsson , J.A. : A survey on spectral{spatial classi cation techniques based on attribute pro les . IEEE Transactions on Geoscience and Remote Sensing 53 ( 5 ), 2335 { 2353 ( 2014 )

15. Heikkinen , V. , Lenz , R. , Jetsu , T. , Parkkinen , J. , Hauta-Kasari , M. , Jaaskelainen, T.: Evaluation and uni cation of some methods for estimating re ectance spectra from RGB images . Journal of the Optical Society of America A 25 ( 10 ), 2444 { 2458 ( 2008 )

16. Lam , A. , Sato , I. : Spectral modeling and relighting of re ective- uorescent scenes . In: Proceedings of the Conference on Computer Vision and Pattern Recognition . pp. 1452 { 1459 . IEEE ( 2013 )

17. Lin , Y. , Finlayson , G.D.: Exposure invariance in spectral reconstruction from RGB images . In: Proceedings of the Color and Imaging Conference . vol. 2019 , pp. 284 { 289 . Society for Imaging Science and Technology ( 2019 )

18. Nguyen , R.M.H. , Prasad , D.K. , Brown , M.S.: Training-based spectral reconstruction from a single RGB image . In: Proceedings of the European Conference on Computer Vision . pp. 186 { 201 . Springer ( 2014 )

19. Pan , Z. , Healey , G. , Prasad , M. , Tromberg , B. : Face recognition in hyperspectral images . IEEE Transactions on Pattern Analysis and Machine Intelligence 25 ( 12 ), 1552 { 1560 ( 2003 )

20. Refaeilzadeh , P. , Tang , L. , Liu , H.: Cross-Validation . In: Encyclopedia of Database Systems , pp. 532 { 538 . Springer US , Boston, MA ( 2009 )

21. Reginska , T.: A regularization parameter in discrete ill-posed problems . SIAM Journal on Scienti c Computing 17 ( 3 ), 740 { 749 ( 1996 )

22. Shi , Z. , Chen , C. , Xiong , Z. , Liu , D. , Wu , F. : HSCNN+: Advanced CNN-based hyperspectral recovery from RGB images . In: Proceedings of the Conference on Computer Vision and Pattern Recognition Workshops . pp. 939 { 947 . IEEE ( 2018 )

23. Tao , C. , Pan , H. , Li , Y. , Zou , Z. : Unsupervised spectral{spatial feature learning with stacked sparse autoencoder for hyperspectral imagery classi cation . IEEE Geoscience and Remote Sensing Letters 12 ( 12 ), 2438 { 2442 ( 2015 )

24. Tikhonov , A.N. , Goncharsky , A.V. , Stepanov , V.V. , Yagola , A.G. : Numerical Methods for the Solution of Ill-posed Problems , vol. 328 . Springer Science & Business Media ( 2013 )

25. Veganzones , M.A. , Tochon , G. , Dalla-Mura , M. , Plaza , A.J. , Chanussot , J.: Hyperspectral image segmentation using a new spectral unmixing-based binary partition tree representation . IEEE Transactions on Image Processing 23 ( 8 ), 3574 { 3589 ( 2014 )

26. Wandell , B.A. : The synthesis and analysis of color images . IEEE Transactions on Pattern Analysis and Machine Intelligence (1) , 2 { 13 ( 1987 )

27. Xu , P. , Xu , H. , Diao , C. , Ye , Z. : Self-training-based spectral image reconstruction for art paintings with multispectral imaging . Applied Optics 56 ( 30 ), 8461 { 8470 ( 2017 )

28. Zhang , Y. , Mou , X. , Wang , G. , Yu , H.: Tensor-based dictionary learning for spectral CT reconstruction . IEEE Transactions on Medical Imaging 36 ( 1 ), 142 { 154 ( 2016 )

29. Zhang , Y. , Xi , Y. , Yang , Q. , Cong , W. , Zhou , J. , Wang , G. : Spectral CT reconstruction with image sparsity and spectral mean . IEEE Transactions on Computational Imaging 2 ( 4 ), 510 { 523 ( 2016 )