Abstract Spectral CT with an energy-discriminating photon
counting detector collects each projection in different energy
channels simultaneously. How to fully utilize such a spectral
dataset for image reconstruction attracts an increasing attention.
Our previous work has shown that the conventional vector based
dictionary learning for low-dose CT can reduce image noise
remarkably. In this paper, we develop a tensor-based
spatio-spectral dictionary learning approach for spectral CT. The
proposed method is compared with two other algorithms using an
in vivo mouse scanned on the MARS micro-CT scanner with a
Medipix MXR CdTe detector. The results demonstrate that the
proposed method outperforms the competing algorithms.
Index Terms Spectral CT, dictionary learning, tensor, CPD
I. I
NTRODUCTION
PECTRALCT with a photon counting detector can
discriminate x-ray energy and collect multiple projections
from different energy channels simultaneously. How to fully
use the acquired information to reconstruct a multichannel
image is a hot topic. A direct way is to treat multichannel data
as a group of conventional CT data, and iteratively reconstruct
each channel image independently with regularization [1].
Because data in different channels are collected from the same
object, the corresponding images are highly correlated and low
rank priori can be applied to encourage correlation among
channels [2, 3]. Furthermore, the x-ray attenuation coefficients
for all channels can be reconstructed as a whole simultaneously.
For example, Semerci et al. treated x-ray attenuation
coefficients at different energy channels as a tensor and the
tensor nuclear norm regularization was designed for spectral
CT reconstruction [4].
In our previous work, vector based dictionary learning was
applied to conventional low dose CT, which demonstrated the
This work was supported in part by the National Natural Science Foundation
of China (NSFC) under Grant Nos. 61172163, 90920003, 61401349 and
61302136, in part by NIH/NIBIB U01 grant EB017140, in part by Natural
Science Basic Research Plan in Shaanxi Province of China (Program No.
2014JQ8317), in part by China Postdoctoral Science Foundation (Program No.
2012M521775).
Yanbo Zhang, Xuanqin Mou and Qiong Xu are with the Institute of Image
Processing and Pattern Recognition, Xi'an Jiaotong University, Xi'an, Shaanxi
710049, China. (e-mail:yanbozhang007@gmail.com,
xqmou@mail.xjtu.edu.cn, xuqiong@mail.xjtu.edu.cn).
Hengyong Yu is with Department of Electrical and Computer Engineering,
University of Massachusetts Lowell, Lowell, MA, 91854 (e-mail:
hengyong-yu@ieee.org).
Ge Wang is with Biomedical Imaging Center/Cluster, CBIS, Rensselaer
Polytechnic Institute, Troy, NY, 12180 (e-mail: wangg6@rpi.edu).
superior reconstruction performance [5]. Later on, it was
extended to spectral CT reconstruction [6, 7]. Very recently, we
proposed a tensor based spatio-temporal dictionary learning
method for dynamic CT, and our preliminary results
demonstrated that tensor based dictionary outperforms the
vector based dictionary [8]. Based on our previous work, in this
paper we will develop a tensor based spatio-spectral dictionary
learning approach for spectral CT reconstruction.
II. T
ENSOR
-B
ASED
D
ICTIONARY
L
EARNING
A tensor is a multidimensional array. An n
th
-order (n-way)
tensor is defined as , whose element in
is , , .
Specifically, for and , tensors are vectors and
matrices, respectively. For , is a
third-order tensor. By an element reordering transformation, a
tensor can be unfolded to a matrix. The mode-k unfolding of
is denoted by , , i.e. ,
and . A third-order tensor can
be factorized into a sum of weighted rank-one tensors using
CANDECOMP/PARAFAC decomposition (CPD) [9]
, (1)
where R is a positive integer,
1
I
r
a ,
2
I
r
b and
3
I
r
c
are normalized vectors for , and is the weight. The
Given a set of third-order training tensors ,
, the tensor based dictionary learning can be
formulated as the following optimization problem,
(2)
where L represents the sparsity level and
1 2 3
( )
{ }
I I I K
k
is the learned tensor dictionary, K is the
number of atoms, is the representation coefficients, the
-mode product of a tensor with a
vector resulting in a third-order tensor, represents L
0
-norm,
namely, the number of nonzero elements, and
F
represents
Frobenius norm. is the k
th
atom which is a
rank-1 tensor, so the atom can be rewritten as
Tensor-based Dictionary Learning for Spectral
CT Reconstruction
Yanbo Zhang, Xuanqin Mou, Qiong Xu, Hengyong Yu, and Ge Wang
S