
Computed tomography (CT) examinations are commonly used to predict lung
nodule malignancy in patients, which are shown to improve noninvasive early
diagnosis of lung cancer. It remains challenging for computational approaches
to achieve performance comparable to experienced radiologists. Here we present
NoduleX, a systematic approach to predict lung nodule malignancy from CT data,
based on deep learning convolutional neural networks (CNN). For training and
validation, we analyze >1000 lung nodules in images from the LIDC/IDRI cohort.
All nodules were identified and classified by four experienced thoracic
radiologists who participated in the LIDC project. NoduleX achieves high
accuracy for nodule malignancy classification, with an AUC of ~0.99. This is
commensurate with the analysis of the dataset by experienced radiologists. Our
approach, NoduleX, provides an effective framework for highly accurate nodule
malignancy prediction with the model trained on a large patient population. Our
results are replicable with software available at
http://bioinformatics.astate.edu/NoduleX.

Image reconstruction in Xray transmission tomography has been an important
research field for decades. In light of data volume increasing faster than
processor speeds, one needs accelerated iterative algorithms to solve the
optimization problem in the Xray CT application. Incremental methods, in which
a subset of data is being used at each iteration to accelerate the
computations, have been getting more popular lately in the machine learning and
mathematical optimization fields. The most popular member of this family of
algorithms in the Xray CT field is orderedsubsets. Even though it performs
well in earlier iterations, the lack of convergence in later iterations is a
known phenomenon. In this paper, we propose two incremental methods that use
Jensen surrogates for the Xray CT application, one stochastic and one
orderedsubsets type. Using measured data, we show that the stochastic variant
we propose outperforms other algorithms, including the gradient descent
counterparts.

In xray coherent scatter tomography, tomographic measurements of the forward
scatter distribution are used to infer scatter densities within a volume. A
radiopaque 2D pattern placed between the object and the detector array enables
the disambiguation between different scatter events. The use of a fan beam
source illumination to speed up data acquisition relative to a pencil beam
presents computational challenges. To facilitate the use of iterative
algorithms based on a penalized Poisson loglikelihood function, efficient
computational implementation of the forward and backward models are needed. Our
proposed implementation exploits physical symmetries and structural properties
of the system and suggests a joint systemalgorithm design, where the system
design choices are influenced by computational considerations, and in turn lead
to reduced reconstruction time. Computationaltime speedups of approximately
146 and 32 are achieved in the computation of the forward and backward models,
respectively. Results validating the forward model and reconstruction algorithm
are presented on simulated analytic and Monte Carlo data.

We consider Xray coherent scatter imaging, where the goal is to reconstruct
momentum transfer profiles (spectral distributions) at each spatial location
from multiplexed measurements of scatter. Each material is characterized by a
unique momentum transfer profile (MTP) which can be used to discriminate
between different materials. We propose an iterative image reconstruction
algorithm based on a Poisson noise model that can account for photonlimited
measurements as well as various second order statistics of the data. To improve
image quality, previous approaches use edgepreserving regularizers to promote
piecewise constancy of the image in the spatial domain while treating each
spectral bin separately. Instead, we propose spectrally grouped regularization
that promotes piecewise constant images along the spatial directions but also
ensures that the MTPs of neighboring spatial bins are similar, if they contain
the same material. We demonstrate that this group regularization results in
improvement of both spectral and spatial image quality. We pursue an
optimization transfer approach where convex decompositions are used to lift the
problem such that all hypervoxels can be updated in parallel and in
closedform. The group penalty introduces a challenge since it is not directly
amendable to these decompositions. We use the alternating directions method of
multipliers (ADMM) to replace the original problem with an equivalent sequence
of subproblems that are amendable to convex decompositions, leading to a
highly parallel algorithm. We demonstrate the performance on real data.

We propose a globally convergent alternating minimization (AM) algorithm for
image reconstruction in transmission tomography, which extends automatic
relevance determination (ARD) to Poisson noise models with Beer's law. The
algorithm promotes solutions that are sparse in the pixel/voxeldifferences
domain by introducing additional latent variables, one for each pixel/voxel,
and then learning these variables from the data using a hierarchical Bayesian
model. Importantly, the proposed AM algorithm is free of any tuning parameters
with image quality comparable to standard penalized likelihood methods. Our
algorithm exploits optimization transfer principles which reduce the problem
into parallel 1D optimization tasks (one for each pixel/voxel), making the
algorithm feasible for largescale problems. This approach considerably reduces
the computational bottleneck of ARD associated with the posterior variances.
Positivity constraints inherent in transmission tomography problems are also
enforced. We demonstrate the performance of the proposed algorithm for xray
computed tomography using synthetic and realworld datasets. The algorithm is
shown to have much better performance than prior ARD algorithms based on
approximate Gaussian noise models, even for high photon flux.

Threedimensional xray CT image reconstruction in baggage scanning in
security applications is an important research field. The variety of materials
to be reconstructed is broader than medical xray imaging. Presence of high
attenuating materials such as metal may cause artifacts if analytical
reconstruction methods are used. Statistical modeling and the resultant
iterative algorithms are known to reduce these artifacts and present good
quantitative accuracy in estimates of linear attenuation coefficients. However,
iterative algorithms may require computations in order to achieve
quantitatively accurate results. For the case of baggage scanning, in order to
provide fast accurate inspection throughput, they must be accelerated
drastically. There are many approaches proposed in the literature to increase
speed of convergence. This paper presents a new method that estimates the
wavelet coefficients of the images in the discrete wavelet transform domain
instead of the image space itself. Initially, surrogate functions are created
around approximation coefficients only. As the iterations proceed, the wavelet
tree on which the updates are made is expanded based on a criterion and detail
coefficients at each level are updated and the tree is expanded this way. For
example, in the smooth regions of the image the detail coefficients are not
updated while the coefficients that represent the highfrequency component
around edges are being updated, thus saving time by focusing computations where
they are needed. This approach is implemented on real data from a SureScan (TM)
x1000 Explosive Detection System and compared to straightforward implementation
of the unregularized alternating minimization of O'Sullivan and Benac [1].