S. Soltani A Matlab Package for DL Approach to CT
2.4.3 Lea_Algorithm_T_D2.m . . . . . . . . . . . . . . . . . . . 13
2.4.4 Lea_Algorithm_T_Dinf.m . . . . . . . . . . . . . . . . . . 13
2.5 The Approximation Error by the Tensor Dictionary . . . . . . . 14
2.5.1 MAET_demo.m . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.5.2 MAE_T.m . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.6 Tomographic Reconstruction with Tensor Dictionary . . . . . . . 16
2.6.1 RecT_demo.m . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.6.2 RecM_Algorithm.m . . . . . . . . . . . . . . . . . . . . . . 16
3 Alphabetic List of Package Functions 19
1 Introduction
1.1 Overall
In our recent works [5, 6] we have developed a two-stage algorithm using a dictio-
nary learning approach in a discrete tomographic reconstruction problem both
with a matrix and a tensor formulation. We first build a training data base from
given training images and then construct a dictionary purely from the training
data by looking for a number of basis elements in terms of the dictionary and the
coefficients/representations of the training data in the dictionary. In these two
works ([5, 6]) dictionary learning is formulated as a non-negative matrix/tensor
formulation with sparsity constraints on the representation, and then the recon-
struction problem is formulated as a sparse approximation problem in a convex
optimization framework. We optimized the dictionary learning problem locally
by using the Alternating Direction Method of Multipliers (ADMM), see e.g., [2].
The tomographic reconstruction problems are solved using the software package
TFOCS (Templates for First-Order Conic Solvers) [1].
This documentation describes the implementations details of the proposed
algorithms for solving the tomographic image reconstruction problem in Mat-
lab. This package is designed for a discretized computed tomography problem
formulation, i.e., we use a matrix formulation of the reconstruction problem
(Ax ≈ b), where b ∈ R
m
contains the noisy measurement data and A is the
system matrix. The vector x ∈ R
n
represents an M × N image, with n = MN,
of absorption coefficients. The total number of tomographic measurement is
m. The system matrix A ∈ R
m×n
is available. We use the AIR Tools [4] to
compute the system matrix A. The discretized tomographic problem is a large
sparse and often ill-posed system and incorporating a priori information about
the solution is a necessity to improve the reconstruction and regularize the solu-
tion. This package aims at providing a computational framework for the use of
training images as priors for the solution in tomographic image reconstruction
in the scheme described above. The algorithms work for any size of the sys-
tem matrix A; however, we are more concerned with underdetermined problems
where m < n, because the need for regularization is even more pronounced in
that scenarios.
This package of Matlab routines provides the user with easy-to-use routines
and demo scripts, based on formulation and numerical algorithms described in
[5, 6]. With the demo script files, the user can easily specifies a few values of
problem parameters to obtain a solution to the dictionary learning problem,
2
评论1
最新资源