Static coded aperture in robotic X-ray tomography systems.

Coded aperture X-ray computed tomography is a computational imaging technique capable of reconstructing inner structures of an object from a reduced set of X-ray projection measurements. Coded apertures are placed in front of the X-ray sources from different views and thus significantly reduce the radiation dose. This paper introduces coded aperture X-ray computed tomography for robotic X-ray systems which offer positioning flexibility. While single coded-aperture 3D tomography was recently introduced for standard trajectory CT scanning, it is shown that significant gains in imaging performance can be attained by simple modifications in the CT scanning trajectories enabled by emerging dual robotic CT systems. In particular, the subject is fixed on a plane and the CT system uniformly rotates around the r -axis which is misaligned with the coordinate axes. A single stationary coded aperture is placed on front of the robotic X-ray source above the plane and the corresponding X-ray projections are measured by a two-dimensional detector on the second arm of the robotic system. The compressive measurements with misalignment enable the reconstruction of high-resolution three-dimensional volumetric images from the low-resolution coded projections on the detector at a sub-sampling rate. An efficient algorithm is proposed to generate the rotation matrix with two basic sub-matrices and thus the forward model is formulated. The stationary coded aperture is designed based on the Pearson product-moment correlation coefficient analysis and the direct binary search algorithm is used to obtain the optimized coded aperture. Simulations using simulated datasets show significant gains in reconstruction performance compared to conventional coded aperture CT systems.

Static coded illumination strategies for low-dose x-ray material decomposition.

Static coded aperture x-ray tomography was introduced recently where a static illumination pattern is used to interrogate an object with a low radiation dose, from which an accurate 3D reconstruction of the object can be attained computationally. Rather than continuously switching the pattern of illumination with each view angle, as traditionally done, static code computed tomography (CT) places a single pattern for all views. The advantages are many, including the feasibility of practical implementation. This paper generalizes this powerful framework to develop single-scan dual-energy coded aperture spectral tomography that enables material characterization at a significantly reduced exposure level. Two sensing strategies are explored: rapid kV switching with a single-static block/unblock coded aperture, and coded apertures with non-uniform thickness. Both systems rely on coded illumination with a plurality of x-ray spectra created by kV switching or 3D coded apertures. The structured x-ray illumination is projected through the objects of interest and measured with standard x-ray energy integrating detectors. Then, based on the tensor representation of projection data, we develop an algorithm to estimate a full set of synthesized measurements that can be used with standard reconstruction algorithms to accurately recover the object in each energy channel. Simulation and experimental results demonstrate the effectiveness of the proposed cost-effective solution to attain material characterization in low-dose dual-energy CT.

Experimental demonstration and optimization of X-ray StaticCodeCT.

As the use of X-ray computed tomography (CT) grows in medical diagnosis, so does the concern for the harm a radiation dose can cause and the biological risks it represents. StaticCodeCT is a new low-dose imaging architecture that uses a single-static coded aperture (CA) in a CT gantry. It exploits the highly correlated data in the projection domain to estimate the unobserved measurements on the detector. We previously analyzed the StaticCodeCT system by emulating the effect of the coded mask on experimental CT data. In contrast, this manuscript presents test-bed reconstructions using an experimental cone-beam X-ray CT system with a CA holder. We analyzed the reconstruction quality using three different techniques to manufacture the CAs: metal additive manufacturing, cold casting, and ceramic additive manufacturing. Furthermore, we propose an optimization method to design the CA pattern based on the algorithm developed for the measurement estimation. The obtained results point to the possibility of the real deployment of StaticCodeCT systems in practice.

Nonlinear reconstruction of coded spectral X-ray CT based on material decomposition.

Coded spectral X-ray computed tomography (CT) based on K-edge filtered illumination is a cost-effective approach to acquire both 3-dimensional structure of objects and their material composition. This approach allows sets of incomplete rays from sparse views or sparse rays with both spatial and spectral encoding to effectively reduce the inspection duration or radiation dose, which is of significance in biological imaging and medical diagnostics. However, reconstruction of spectral CT images from compressed measurements is a nonlinear and ill-posed problem. This paper proposes a material-decomposition-based approach to directly solve the reconstruction problem, without estimating the energy-binned sinograms. This approach assumes that the linear attenuation coefficient map of objects can be decomposed into a few basis materials that are separable in the spectral and space domains. The nonlinear problem is then converted to the reconstruction of the mass density maps of the basis materials. The dimensionality of the optimization variables is thus effectively reduced to overcome the ill-posedness. An alternating minimization scheme is used to solve the reconstruction with regularizations of weighted nuclear norm and total variation. Compared to the state-of-the-art reconstruction method for coded spectral CT, the proposed method can significantly improve the reconstruction quality. It is also capable of reconstructing the spectral CT images at two additional energy bins from the same set of measurements, thus providing more spectral information of the object.

StaticCodeCT: single coded aperture tensorial X-ray CT.

Coded aperture X-ray CT (CAXCT) is a new low-dose imaging technology that promises far-reaching benefits in industrial and clinical applications. It places various coded apertures (CA) at a time in front of the X-ray source to partially block the radiation. The ill-posed inverse reconstruction problem is then solved using l1-norm-based iterative reconstruction methods. Unfortunately, to attain high-quality reconstructions, the CA patterns must change in concert with the view-angles making the implementation impractical. This paper proposes a simple yet radically different approach to CAXCT, which is coined StaticCodeCT, that uses a single-static CA in the CT gantry, thus making the imaging system amenable for practical implementations. Rather than using conventional compressed sensing algorithms for recovery, we introduce a new reconstruction framework for StaticCodeCT. Namely, we synthesize the missing measurements using low-rank tensor completion principles that exploit the multi-dimensional data correlation and low-rank nature of a 3-way tensor formed by stacking the 2D coded CT projections. Then, we use the FDK algorithm to recover the 3D object. Computational experiments using experimental projection measurements exhibit up to 10% gains in the normalized root mean square distance of the reconstruction using the proposed method compared with those attained by alternative low-dose systems.

Fast optimization of coded apertures in X-ray computed tomography.

Coded aperture X-ray computed tomography (CAXCT) is a novel X-ray imaging system capable of reconstructing high quality images from a reduced set of measurements. Coded apertures are placed in front of the X-ray source in CAXCT so as to obtain patterned projections onto a detector array. Then, compressive sensing (CS) reconstruction algorithms are used to reconstruct the linear attenuation coefficients. The coded aperture is an important factor that influences the point spread function (PSF), which in turn determines the capability to sample the linear attenuation coefficients of the object. A coded aperture optimization approach was recently proposed based on the coherence of the system matrix; however, this algorithm is memory intensive and it is not able to optimize the coded apertures for large image sizes required in many applications. This paper introduces a significantly more efficient approach for coded aperture optimization that reduces the memory requirements and the execution time by orders of magnitude. The features are defined as the inner product of the vectors representing the geometric paths of the X-rays with the sparse basis representation of the object; therefore, the algorithm aims to find a subset of features that minimizes the information loss compared to the complete set of projections. This subset corresponds to the unblocking elements in the optimized coded apertures. The proposed approach solves the memory and runtime limitations of the previously proposed algorithm and provides a significant gain in the reconstruction image quality compared to that attained by random coded apertures in both simulated datasets and real datasets.

Coded aperture optimization in compressive X-ray tomography: a gradient descent approach.

Coded aperture X-ray computed tomography (CT) has the potential to revolutionize X-ray tomography systems in medical imaging and air and rail transit security - both areas of global importance. It allows either a reduced set of measurements in X-ray CT without degradation in image reconstruction, or measure multiplexed X-rays to simplify the sensing geometry. Measurement reduction is of particular interest in medical imaging to reduce radiation, and airport security often imposes practical constraints leading to limited angle geometries. Coded aperture compressive X-ray CT places a coded aperture pattern in front of the X-ray source in order to obtain patterned projections onto a detector. Compressive sensing (CS) reconstruction algorithms are then used to recover the image. To date, the coded illumination patterns used in conventional CT systems have been random. This paper addresses the code optimization problem for general tomography imaging based on the point spread function (PSF) of the system, which is used as a measure of the sensing matrix quality which connects to the restricted isometry property (RIP) and coherence of the sensing matrix. The methods presented are general, simple to use, and can be easily extended to other imaging systems. Simulations are presented where the peak signal to noise ratios (PSNR) of the reconstructed images using optimized coded apertures exhibit significant gain over those attained by random coded apertures. Additionally, results using real X-ray tomography projections are presented.

Coded aperture optimization for compressive X-ray tomosynthesis.

Radiation dose is a concern in X-ray tomographic imaging; coded aperture compressive X-ray tomosynthesis is an approach used to reduce radiation. It places a coded aperture in front of an X-ray source in order to obtain 2D patterned projections of a three-dimensional object onto a detector plane. By using different coded apertures in a multiple source system, multiplexed projections can be obtained instead of sequential projections as in conventional tomosynthesis systems. Compressed sensing (CS) reconstruction algorithms are then used to recover the three-dimensional data cube. An optimization approach to design the structure of the coded apertures in a multiple source compressive X-ray tomosynthesis imaging system is presented. A uniform energy criteria on the voxels and detector elements is used so that the object is uniformly sensed and the elements of the detector plane uniformly sense the information. Simulations and experimental results for optimized coded apertures are shown, and their performance is compared to the use of random coded apertures.