Full data consistency conditions for cone-beam projections with sources on a plane

AIRPORT: A Data Consistency Constrained Deep Temporal Extrapolation Method To Improve Temporal Resolution In Contrast Enhanced CT Imaging

Typical tomographic image reconstruction methods require that the imaged object is static and stationary during the time window to acquire a minimally complete data set. The violation of this requirement leads to temporal-averaging errors in the reconstructed images. For a fixed gantry rotation speed, to reduce the errors, it is desired to reconstruct images using data acquired over a narrower angular range, i.e., with a higher temporal resolution. However, image reconstruction with a narrower angular range violates the data sufficiency condition, resulting in severe data-insufficiency-induced errors. The purpose of this work is to decouple the trade-off between these two types of errors in contrast-enhanced computed tomography (CT) imaging. We demonstrated that using the developed data consistency constrained deep temporal extrapolation method (AIRPORT), the entire time-varying imaged object can be accurately reconstructed with 40 frames-per-second temporal resolution, the time window needed to acquire a single projection view data using a typical C-arm cone-beam CT system. AIRPORT is applicable to general non-sparse imaging tasks using a single short-scan data acquisition.

Equal parallel and cone-beam projections: a curious property of D-symmetric object functions

In tomographic image reconstruction, the object density function is the unknown quantity whose projections are measured by the scanner. In the three-dimensional case, we define the D-reflection of such a density function as the object obtained by a particular weighted reflection about the planez=D, and a D-symmetric function as one whose D-reflection is equal to itself. D-symmetric object functions have the curious property that their parallel projection onto the detector planez=Dis equal to their cone-beam projection onto the same detector with x-ray source location at the origin. Much more remarkable is the additional fact that for any fixed D-symmetric object,everyoblique parallel projection onto this same detector plane equals the cone-beam projection for a corresponding source location. The mathematical proof is straight forward but not particularly enlightening, and we also provide here an alternative physical demonstration that explains the various weighting terms in the context of classical tomosynthesis. Furthermore, we clarify the distinction between the new formulation presented here, and the original formulation of Edholm and co-workers who obtained similar properties but for a pair of objects whose divergent and parallel projections matched, but with no D-symmetry. We do not claim any immediate imaging application or useful physics from these notions, but we briefly comment on consequences for methods that apply data consistency conditions in image reconstruction.

Motion estimation and correction in SPECT, PET and CT

Patient motion impacts single photon emission computed tomography (SPECT), positron emission tomography (PET) and X-ray computed tomography (CT) by giving rise to projection data inconsistencies that can manifest as reconstruction artifacts, thereby degrading image quality and compromising accurate image interpretation and quantification. Methods to estimate and correct for patient motion in SPECT, PET and CT have attracted considerable research effort over several decades. The aims of this effort have been two-fold: to estimate relevant motion fields characterizing the various forms of voluntary and involuntary motion; and to apply these motion fields within a modified reconstruction framework to obtain motion-corrected images. The aims of this review are to outline the motion problem in medical imaging and to critically review published methods for estimating and correcting for the relevant motion fields in clinical and preclinical SPECT, PET and CT. Despite many similarities in how motion is handled between these modalities, utility and applications vary based on differences in temporal and spatial resolution. Technical feasibility has been demonstrated in each modality for both rigid and non-rigid motion, but clinical feasibility remains an important target. There is considerable scope for further developments in motion estimation and correction, and particularly in data-driven methods that will aid clinical utility. State-of-the-art machine learning methods may have a unique role to play in this context.

An Exact and Fast CBCT Reconstruction via Pseudo-Polar Fourier Transform based Discrete Grangeat's Formula

The recent application of Fourier Based Iterative Reconstruction Method (FIRM) has made it possible to achieve high-quality 2D images from a fan beam Computed Tomography (CT) scan with a limited number of projections in a fast manner. The proposed methodology in this article is designed to provide 3D Radon space in linogram fashion to facilitate the use of FIRM with cone beam projections (CBP) for the reconstruction of 3D images in a sparse view angles Cone Beam CT (CBCT). For this reason, in the first phase, the 3D Radon space is generated using CBP data after discretization and optimization of the famous Grangeat's formula. The method used in this process involves fast Pseudo Polar Fourier transform (PPFT) based on 2D and 3D Discrete Radon Transformation (DRT) algorithms with no wraparound effects. In the second phase, we describe reconstruction of the objects with available Radon values, using direct inverse of 3D PPFT. The method presented in this section eliminates noises caused by interpolation from polar to Cartesian space and exhibits no thorn, V-shaped and wrinkle artifacts. This method reduces the complexity to for images of size n × n × n The Cone to Radon conversion (Cone2Radon) Toolbox in the first phase and MATLAB/ Python toolbox in the second phase were tested on three digital phantoms and experiments demonstrate fast and accurate cone beam image reconstruction due to proposed.

Spectrum Estimation-Guided Iterative Reconstruction Algorithm for Dual Energy CT

X-ray spectrum plays a very important role in dual energy computed tomography (DECT) reconstruction. Because it is difficult to measure x-ray spectrum directly in practice, efforts have been devoted into spectrum estimation by using transmission measurements. These measurement methods are independent of the image reconstruction, which bring extra cost and are time consuming. Furthermore, the estimated spectrum mismatch would degrade the quality of the reconstructed images. In this paper, we propose a spectrum estimation-guided iterative reconstruction algorithm for DECT which aims to simultaneously recover the spectrum and reconstruct the image. The proposed algorithm is formulated as an optimization framework combining spectrum estimation based on model spectra representation, image reconstruction, and regularization for noise suppression. To resolve the multi-variable optimization problem of simultaneously obtaining the spectra and images, we introduce the block coordinate descent (BCD) method into the optimization iteration. Both the numerical simulations and physical phantom experiments are performed to verify and evaluate the proposed method. The experimental results validate the accuracy of the estimated spectra and reconstructed images under different noise levels. The proposed method obtains a better image quality compared with the reconstructed images from the known exact spectra and is robust in noisy data applications.

Motion compensation for cone-beam CT using Fourier consistency conditions

In cone-beam CT, involuntary patient motion and inaccurate or irreproducible scanner motion substantially degrades image quality. To avoid artifacts this motion needs to be estimated and compensated during image reconstruction. In previous work we showed that Fourier consistency conditions (FCC) can be used in fan-beam CT to estimate motion in the sinogram domain. This work extends the FCC to [Formula: see text] cone-beam CT. We derive an efficient cost function to compensate for [Formula: see text] motion using [Formula: see text] detector translations. The extended FCC method have been tested with five translational motion patterns, using a challenging numerical phantom. We evaluated the root-mean-square-error and the structural-similarity-index between motion corrected and motion-free reconstructions. Additionally, we computed the mean-absolute-difference (MAD) between the estimated and the ground-truth motion. The practical applicability of the method is demonstrated by application to respiratory motion estimation in rotational angiography, but also to motion correction for weight-bearing imaging of knees. Where the latter makes use of a specifically modified FCC version which is robust to axial truncation. The results show a great reduction of motion artifacts. Accurate estimation results were achieved with a maximum MAD value of 708 μm and 1184 μm for motion along the vertical and horizontal detector direction, respectively. The image quality of reconstructions obtained with the proposed method is close to that of motion corrected reconstructions based on the ground-truth motion. Simulations using noise-free and noisy data demonstrate that FCC are robust to noise. Even high-frequency motion was accurately estimated leading to a considerable reduction of streaking artifacts. The method is purely image-based and therefore independent of any auxiliary data.

John's Equation-based Consistency Condition and Corrupted Projection Restoration in Circular Trajectory Cone Beam CT

In transmitted X-ray tomography imaging, the acquired projections may be corrupted for various reasons, such as defective detector cells and beam-stop array scatter correction problems. In this study, we derive a consistency condition for cone-beam projections and propose a method to restore lost data in corrupted projections. In particular, the relationship of the geometry parameters in circular trajectory cone-beam computed tomography (CBCT) is utilized to convert an ultra-hyperbolic partial differential equation (PDE) into a second-order PDE. The second-order PDE is then transformed into a first-order ordinary differential equation in the frequency domain. The left side of the equation for the newly derived consistency condition is the projection derivative of the current and adjacent views, whereas the right side is the projection derivative of the geometry parameters. A projection restoration method is established based on the newly derived equation to restore corrupted data in projections in circular trajectory CBCT. The proposed method is tested in beam-stop array scatter correction, metal artifact reduction, and abnormal pixel correction cases to evaluate the performance of the consistency condition and corrupted projection restoration method. Qualitative and quantitative results demonstrate that the present method has considerable potential in restoring lost data in corrupted projections.

Epipolar Consistency in Transmission Imaging

This paper presents the derivation of the Epipolar Consistency Conditions (ECC) between two X-ray images from the Beer-Lambert law of X-ray attenuation and the Epipolar Geometry of two pinhole cameras, using Grangeat's theorem. We motivate the use of Oriented Projective Geometry to express redundant line integrals in projection images and define a consistency metric, which can be used, for instance, to estimate patient motion directly from a set of X-ray images. We describe in detail the mathematical tools to implement an algorithm to compute the Epipolar Consistency Metric and investigate its properties with detailed random studies on both artificial and real FD-CT data. A set of six reference projections of the CT scan of a fish were used to evaluate accuracy and precision of compensating for random disturbances of the ground truth projection matrix using an optimization of the consistency metric. In addition, we use three X-ray images of a pumpkin to prove applicability to real data. We conclude, that the metric might have potential in applications related to the estimation of projection geometry. By expression of redundancy between two arbitrary projection views, we in fact support any device or acquisition trajectory which uses a cone-beam geometry. We discuss certain geometric situations, where the ECC provide the ability to correct 3D motion, without the need for 3D reconstruction.