State of the art: iterative CT reconstruction techniques

Owing to recent advances in computing power, iterative reconstruction (IR) algorithms have become a clinically viable option in computed tomographic (CT) imaging. Substantial evidence is accumulating about the advantages of IR algorithms over established analytical methods, such as filtered back projection. IR improves image quality through cyclic image processing. Although all available solutions share the common mechanism of artifact reduction and/or potential for radiation dose savings, chiefly due to image noise suppression, the magnitude of these effects depends on the specific IR algorithm. In the first section of this contribution, the technical bases of IR are briefly reviewed and the currently available algorithms released by the major CT manufacturers are described. In the second part, the current status of their clinical implementation is surveyed. Regardless of the applied IR algorithm, the available evidence attests to the substantial potential of IR algorithms for overcoming traditional limitations in CT imaging

Tomography of atomic number and density of materials using dual-energy imaging and the Alvarez and Macovski attenuation model

Dual-energy computed tomography and the Alvarez and Macovski [Phys. Med. Biol. 21, 733 (1976)] transmitted intensity (AMTI) model were used in this study to estimate the maps of density (ρ) and atomic number (Z) of mineralogical samples. In this method, the attenuation coefficients are represented [Alvarez and Macovski, Phys. Med. Biol. 21, 733 (1976)] in the form of the two most important interactions of X-rays with atoms that is, photoelectric absorption (PE) and Compton scattering (CS). This enables material discrimination as PE and CS are, respectively, dependent on the atomic number (Z) and density (ρ) of materials [Alvarez and Macovski, Phys. Med. Biol. 21, 733 (1976)]. Dual-energy imaging is able to identify sample materials even if the materials have similar attenuation coefficients at single-energy spectrum. We use the full model rather than applying one of several applied simplified forms [Alvarez and Macovski, Phys. Med. Biol. 21, 733 (1976); Siddiqui et al., SPE Annual Technical Conference and Exhibition (Society of Petroleum Engineers, 2004); Derzhi, U.S. patent application 13/527,660 (2012); Heismann et al., J. Appl. Phys. 94, 2073–2079 (2003); Park and Kim, J. Korean Phys. Soc. 59, 2709 (2011); Abudurexiti et al., Radiol. Phys. Technol. 3, 127–135 (2010); and Kaewkhao et al., J. Quant. Spectrosc. Radiat. Transfer 109, 1260–1265 (2008)]. This paper describes the tomographic reconstruction of ρ and Z maps of mineralogical samples using the AMTI model. The full model requires precise knowledge of the X-ray energy spectra and calibration of PE and CS constants and exponents of atomic number and energy that were estimated based on fits to simulations and calibration measurements. The estimated ρ and Z images of the samples used in this paper yield average relative errors of 2.62% and 1.19% and maximum relative errors of 2.64% and 7.85%, respectively. Furthermore, we demonstrate that the method accounts for the beam hardening effect in density (ρ) and atomic number (Z) reconstructions to a significant extent.S.J.L., G.R.M., and A.M.K. acknowledge funding through the DigiCore consortium and the support of a linkage grant (LP150101040) from the Australian Research Council and FEI Company

The Estimation and Correction of Rigid Motion in Helical Computed Tomography

X-ray CT is a tomographic imaging tool used in medicine and industry. Although technological developments have significantly improved the performance of CT systems, the accuracy of images produced by state-of-the-art scanners is still often limited by artefacts due to object motion. To tackle this problem, a number of motion estimation and compensation methods have been proposed. However, no methods with the demonstrated ability to correct for rigid motion in helical CT scans appear to exist. The primary aims of this thesis were to develop and evaluate effective methods for the estimation and correction of arbitrary six degree-of-freedom rigid motion in helical CT. As a first step, a method was developed to accurately estimate object motion during CT scanning with an optical tracking system, which provided sub-millimetre positional accuracy. Subsequently a motion correction method, which is analogous to a method previously developed for SPECT, was adapted to CT. The principle is to restore projection consistency by modifying the source-detector orbit in response to the measured object motion and reconstruct from the modified orbit with an iterative reconstruction algorithm. The feasibility of this method was demonstrated with a rapidly moving brain phantom, and the efficacy of correcting for a range of human head motions acquired from healthy volunteers was evaluated in simulations. The methods developed were found to provide accurate and artefact-free motion corrected images with most types of head motion likely to be encountered in clinical CT imaging, provided that the motion was accurately known. The method was also applied to CT data acquired on a hybrid PET/CT scanner demonstrating its versatility. Its clinical value may be significant by reducing the need for repeat scans (and repeat radiation doses), anesthesia and sedation in patient groups prone to motion, including young children

Calibration by differentiation – Self‐supervised calibration for X‐ray microscopy using a differentiable cone‐beam reconstruction operator

High‐resolution X‐ray microscopy (XRM) is gaining interest for biological investigations of extremely small‐scale structures. XRM imaging of bones in living mice could provide new insights into the emergence and treatment of osteoporosis by observing osteocyte lacunae, which are holes in the bone of few micrometres in size. Imaging living animals at that resolution, however, is extremely challenging and requires very sophisticated data processing converting the raw XRM detector output into reconstructed images. This paper presents an open‐source, differentiable reconstruction pipeline for XRM data which analytically computes the final image from the raw measurements. In contrast to most proprietary reconstruction software, it offers the user full control over each processing step and, additionally, makes the entire pipeline deep learning compatible by ensuring differentiability. This allows fitting trainable modules both before and after the actual reconstruction step in a purely data‐driven way using the gradient‐based optimizers of common deep learning frameworks. The value of such differentiability is demonstrated by calibrating the parameters of a simple cupping correction module operating on the raw projection images using only a self‐supervisory quality metric based on the reconstructed volume and no further calibration measurements. The retrospective calibration directly improves image quality as it avoids cupping artefacts and decreases the difference in grey values between outer and inner bone by 68–94%. Furthermore, it makes the reconstruction process entirely independent of the XRM manufacturer and paves the way to explore modern deep learning reconstruction methods for arbitrary XRM and, potentially, other flat‐panel computed tomography systems. This exemplifies how differentiable reconstruction can be leveraged in the context of XRM and, hence, is an important step towards the goal of reducing the resolution limit of in vivo bone imaging to the single micrometre domain

Reduced and coded sensing methods for x-ray based security

Current x-ray technologies provide security personnel with non-invasive sub-surface imaging and contraband detection in various portal screening applications such as checked and carry-on baggage as well as cargo. Computed tomography (CT) scanners generate detailed 3D imagery in checked bags; however, these scanners often require significant power, cost, and space. These tomography machines are impractical for many applications where space and power are often limited such as checkpoint areas. Reducing the amount of data acquired would help reduce the physical demands of these systems. Unfortunately this leads to the formation of artifacts in various applications, thus presenting significant challenges in reconstruction and classification. As a result, the goal is to maintain a certain level of image quality but reduce the amount of data gathered. For the security domain this would allow for faster and cheaper screening in existing systems or allow for previously infeasible screening options due to other operational constraints. While our focus is predominantly on security applications, many of the techniques can be extended to other fields such as the medical domain where a reduction of dose can allow for safer and more frequent examinations. This dissertation aims to advance data reduction algorithms for security motivated x-ray imaging in three main areas: (i) development of a sensing aware dimensionality reduction framework, (ii) creation of linear motion tomographic method of object scanning and associated reconstruction algorithms for carry-on baggage screening, and (iii) the application of coded aperture techniques to improve and extend imaging performance of nuclear resonance fluorescence in cargo screening. The sensing aware dimensionality reduction framework extends existing dimensionality reduction methods to include knowledge of an underlying sensing mechanism of a latent variable. This method provides an improved classification rate over classical methods on both a synthetic case and a popular face classification dataset. The linear tomographic method is based on non-rotational scanning of baggage moved by a conveyor belt, and can thus be simpler, smaller, and more reliable than existing rotational tomography systems at the expense of more challenging image formation problems that require special model-based methods. The reconstructions for this approach are comparable to existing tomographic systems. Finally our coded aperture extension of existing nuclear resonance fluorescence cargo scanning provides improved observation signal-to-noise ratios. We analyze, discuss, and demonstrate the strengths and challenges of using coded aperture techniques in this application and provide guidance on regimes where these methods can yield gains over conventional methods

Interferometric Optical Tomography

Embodiments of the present disclosure provide systems and methods for constructing a profile of sample object. Briefly described, in architecture, one embodiment of the system, among others, can be implemented as follows. An interferometer device is used to collect interference images of a sample object at a sequence of angles around the sample object. Accordingly, a controller device rotates the sample object to enable acquisition of the interference images; and a projection generator produces projections of the sample object from the interference images at the sequence of angles. Further, a tomographic device constructs the profile of the optical device from the projections of the interference images. The profile is capable of characterizing small index variations of less than 1x10?4. Other systems and methods are also included.Georgia Tech Research Corporatio

Algorithm for the reconstruction of dynamic objects in CT-scanning using optical flow

Computed Tomography is a powerful imaging technique that allows non-destructive visualization of the interior of physical objects in different scientific areas. In traditional reconstruction techniques the object of interest is mostly considered to be static, which gives artefacts if the object is moving during the data acquisition. In this paper we present a method that, given only scan results of multiple successive scans, can estimate the motion and correct the CT-images for this motion assuming that the motion field is smooth over the complete domain using optical flow. The proposed method is validated on simulated scan data. The main contribution is that we show we can use the optical flow technique from imaging to correct CT-scan images for motion