A unified statistical framework for material decomposition using multienergy photon counting x-ray detectors

Estimation of Basis Line-Integrals in a Spectral Distortion-Modeled Photon Counting Detector Using Low-Order Polynomial Approximation of X-ray Transmittance

Photon counting detector (PCD)-based computed tomography exploits spectral information from a transmitted x-ray spectrum to estimate basis line-integrals. The recorded spectrum, however, is distorted and deviates from the transmitted spectrum due to spectral response effect (SRE). Therefore, the SRE needs to be compensated for when estimating basis line-integrals. One approach is to incorporate the SRE model with an incident spectrum into the PCD measurement model and the other approach is to perform a calibration process that inherently includes both the SRE and the incident spectrum. A maximum likelihood estimator can be used to the former approach, which guarantees asymptotic optimality; however, a heavy computational burden is a concern. Calibration-based estimators are a form of the latter approach. They can be very efficient; however, a heuristic calibration process needs to be addressed. In this paper, we propose a computationally efficient three-step estimator for the former approach using a low-order polynomial approximation of x-ray transmittance. The low-order polynomial approximation can change the original non-linear estimation method to a two-step linearized approach followed by an iterative bias correction step. We show that the calibration process is required only for the bias correction step and prove that it converges to the unbiased solution under practical assumptions. Extensive simulation studies validate the proposed method and show that the estimation results are comparable to those of the ML estimator while the computational time is reduced substantially.

Sparse-view spectral CT reconstruction using spectral patch-based low-rank penalty

Spectral computed tomography (CT) is a promising technique with the potential for improving lesion detection, tissue characterization, and material decomposition. In this paper, we are interested in kVp switching-based spectral CT that alternates distinct kVp X-ray transmissions during gantry rotation. This system can acquire multiple X-ray energy transmissions without additional radiation dose. However, only sparse views are generated for each spectral measurement; and the spectra themselves are limited in number. To address these limitations, we propose a penalized maximum likelihood method using spectral patch-based low-rank penalty, which exploits the self-similarity of patches that are collected at the same position in spectral images. The main advantage is that the relatively small number of materials within each patch allows us to employ the low-rank penalty that is less sensitive to intensity changes while preserving edge directions. In our optimization formulation, the cost function consists of the Poisson log-likelihood for X-ray transmission and the nonconvex patch-based low-rank penalty. Since the original cost function is difficult to minimize directly, we propose an optimization method using separable quadratic surrogate and concave convex procedure algorithms for the log-likelihood and penalty terms, which results in an alternating minimization that provides a computational advantage because each subproblem can be solved independently. We performed computer simulations and a real experiment using a kVp switching-based spectral CT with sparse-view measurements, and compared the proposed method with conventional algorithms. We confirmed that the proposed method improves spectral images both qualitatively and quantitatively. Furthermore, our GPU implementation significantly reduces the computational cost.