Search

Ultra low radiation dose in X-ray Computed Tomography (CT) is an important clinical objective in order to minimize the risk of carcinogenesis. Compressed Sensing (CS) enables significant reductions in radiation dose to be achieved by producing diagnostic images from a limited number of CT projections. However, the excessive computation time that conventional CS-based CT reconstruction typically requires has limited clinical implementation. In this paper, we first demonstrate that a thorough analysis of CT reconstruction through a Maximum a Posteriori objective function results in a weighted compressive sensing problem. This analysis enables us to formulate a low dose fan beam and helical cone beam CT reconstruction. Subsequently, we provide an efficient solution to the formulated CS problem based on a Fast Composite Splitting Algorithm-Latent Expected Maximization (FCSA-LEM) algorithm. In the proposed method we use pseudo polar Fourier transform as the measurement matrix in order to decrease the computational complexity; and rebinning of the projections to parallel rays in order to extend its application to fan beam and helical cone beam scans. The weight involved in the proposed weighted CS model, denoted by Error Adaptation Weight (EAW), is calculated based on the statistical characteristics of CT reconstruction and is a function of Poisson measurement noise and rebinning interpolation error. Simulation results show that low computational complexity of the proposed method made the fast recovery of the CT images possible and using EAW reduces the reconstruction error by one order of magnitude. Recovery of a high quality 512$\times$ 512 image was achieved in less than 20 sec on a desktop computer without numerical optimizations.

A time domain system of equations is proposed to model elastic wave propagation in an unbounded two-dimensional anisotropic solid using perfectly matched layer (PML). Starting from a system of first-order frequency domain stress-velocity equations and using complex coordinate stretching approach with a two-parameter stretch function, a second-order formulation is obtained. The final system, which consists of just two second order equations along with four auxiliary equations, is smaller than existing formulations, thereby simplifying the problem and reducing the computational cost. The discrete stability of the solutions for a given mesh size is examined with the help of a plane-wave analysis of the corresponding continuous problem. It is shown that increasing the scaling parameter of the stretch function leads to significant stability improvements for certain anisotropic media that have known issues. Numerical computations for different isotropic and anisotropic media are used to illustrate the results.