“…It can be determined in a two-step procedure, i.e., image reconstruction for spectral images and then material decomposition from these spectral images to material compositions [3], [6]–[16], or alternatively material-specific sinogram decomposition and then material reconstruction [4], [17]–[19]. Various iterative reconstruction models have been developed [20], with energy-by-energy reconstruction [3], [4], [9], [11], [17]–[19] and joint reconstruction [7], [10], [15], [16], such as total variation (TV) sparsity [14], [16], HYPR algorithm [8], tight frame sparsity [3], [11], bilateral filtration [12], [13], patch-based low-rank model [15], rank-and-sparsity decomposition model [7] and its tensor version [10]. In order to fully utilize the image similarity in the spectral dimension, the joint reconstruction is a natural formulation [7], [10], [15], [16].…”