A new data consistency condition for fan-beam projection data

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.

Optimization based beam-hardening correction in CT under data integral invariant constraint

In computed tomography (CT), the polychromatic characteristics of x-ray photons, which are emitted from a source, interact with materials and are absorbed by a detector, may lead to beam-hardening effect in signal detection and image formation, especially in situations where materials of high attenuation (e.g. the bone or metal implants) are in the x-ray beam. Usually, a beam-hardening correction (BHC) method is used to suppress the artifacts induced by bone or other objects of high attenuation, in which a calibration-oriented iterative operation is carried out to determine a set of parameters for all situations. Based on the Helgasson-Ludwig consistency condition (HLCC), an optimization based method has been proposed by turning the calibration-oriented iterative operation of BHC into solving an optimization problem sustained by projection data. However, the optimization based HLCC-BHC method demands the engagement of a large number of neighboring projection views acquired at relatively high and uniform angular sampling rate, hindering its application in situations where the angular sampling in projection view is sparse or non-uniform. By defining an objective function based on the data integral invariant constraint (DIIC), we again turn BHC into solving an optimization problem sustained by projection data. As it only needs a pair of projection views at any view angle, the proposed BHC method can be applicable in the challenging scenarios mentioned above. Using the projection data simulated by computer, we evaluate and verify the proposed optimization based DIIC-BHC method's performance. Moreover, with the projection data of a head scan by a multi-detector row MDCT, we show the proposed DIIC-BHC method's utility in clinical applications.

Geometry and energy constrained projection extension

BACKGROUND

In clinical computed tomography (CT) applications, when a patient is obese or improperly positioned, the final tomographic scan is often partially truncated. Images directly reconstructed by the conventional reconstruction algorithms suffer from severe cupping and direct current bias artifacts. Moreover, the current methods for projection extension have limitations that preclude incorporation from clinical workflows, such as prohibitive computational time for iterative reconstruction, extra radiation dose, hardware modification, etc.METHOD:In this study, we first established a geometrical constraint and estimated the patient habitus using a modified scout configuration. Then, we established an energy constraint using the integral invariance of fan-beam projections. Two constraints were extracted from the existing CT scan process with minimal modification to the clinical workflows. Finally, we developed a novel dual-constraint based optimization model that can be rapidly solved for projection extrapolation and accurate local reconstruction.

RESULTS

Both numerical phantom and realistic patient image simulations were performed, and the results confirmed the effectiveness of our proposed approach.

CONCLUSION

We establish a dual-constraint-based optimization model and correspondingly develop an accurate extrapolation method for partially truncated projections. The proposed method can be readily integrated into the clinical workflow and efficiently solved by using a one-dimensional optimization algorithm. Moreover, it is robust for noisy cases with various truncations and can be further accelerated by GPU based parallel computing.

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.

Volumetric CT with sparse detector arrays (and application to Si-strip photon counters)

Novel x-ray medical imaging sensors, such as photon counting detectors (PCDs) and large area CCD and CMOS cameras can involve irregular and/or sparse sampling of the detector plane. Application of such detectors to CT involves undersampling that is markedly different from the commonly considered case of sparse angular sampling. This work investigates volumetric sampling in CT systems incorporating sparsely sampled detectors with axial and helical scan orbits and evaluates performance of model-based image reconstruction (MBIR) with spatially varying regularization in mitigating artifacts due to sparse detector sampling. Volumetric metrics of sampling density and uniformity were introduced. Penalized-likelihood MBIR with a spatially varying penalty that homogenized resolution by accounting for variations in local sampling density (i.e. detector gaps) was evaluated. The proposed methodology was tested in simulations and on an imaging bench based on a Si-strip PCD (total area 5 cm  ×  25 cm) consisting of an arrangement of line sensors separated by gaps of up to 2.5 mm. The bench was equipped with translation/rotation stages allowing a variety of scanning trajectories, ranging from a simple axial acquisition to helical scans with variable pitch. Statistical (spherical clutter) and anthropomorphic (hand) phantoms were considered. Image quality was compared to that obtained with a conventional uniform penalty in terms of structural similarity index (SSIM), image uniformity, spatial resolution, contrast, and noise. Scan trajectories with intermediate helical width (~10 mm longitudinal distance per 360° rotation) demonstrated optimal tradeoff between the average sampling density and the homogeneity of sampling throughout the volume. For a scan trajectory with 10.8 mm helical width, the spatially varying penalty resulted in significant visual reduction of sampling artifacts, confirmed by a 10% reduction in minimum SSIM (from 0.88 to 0.8) and a 40% reduction in the dispersion of SSIM in the volume compared to the constant penalty (both penalties applied at optimal regularization strength). Images of the spherical clutter and wrist phantoms confirmed the advantages of the spatially varying penalty, showing a 25% improvement in image uniformity and 1.8  ×  higher CNR (at matched spatial resolution) compared to the constant penalty. The studies elucidate the relationship between sampling in the detector plane, acquisition orbit, sampling of the reconstructed volume, and the resulting image quality. They also demonstrate the benefit of spatially varying regularization in MBIR for scenarios with irregular sampling patterns. Such findings are important and integral to the incorporation of a sparsely sampled Si-strip PCD in CT imaging.

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

Cone-beam consistency conditions (also known as range conditions) are mathematical relationships between different cone-beam projections, and they therefore describe the redundancy or overlap of information between projections. These redundancies have often been exploited for applications in image reconstruction. In this work we describe new consistency conditions for cone-beam projections whose source positions lie on a plane. A further restriction is that the target object must not intersect this plane. The conditions require that moments of the cone-beam projections be polynomial functions of the source positions, with some additional constraints on the coefficients of the polynomials. A precise description of the consistency conditions is that the four parameters of the cone-beam projections (two for the detector, two for the source position) can be expressed with just three variables, using a certain formulation involving homogeneous polynomials. The main contribution of this work is our demonstration that these conditions are not only necessary, but also sufficient. Thus the consistency conditions completely characterize all redundancies, so no other independent conditions are possible and in this sense the conditions are full. The idea of the proof is to use the known consistency conditions for 3D parallel projections, and to then apply a 1996 theorem of Edholm and Danielsson that links parallel to cone-beam projections. The consistency conditions are illustrated with a simulation example.

Software architecture for multi-bed FDK-based reconstruction in X-ray CT scanners

Most small-animal X-ray computed tomography (CT) scanners are based on cone-beam geometry with a flat-panel detector orbiting in a circular trajectory. Image reconstruction in these systems is usually performed by approximate methods based on the algorithm proposed by Feldkamp et al. (FDK). Besides the implementation of the reconstruction algorithm itself, in order to design a real system it is necessary to take into account numerous issues so as to obtain the best quality images from the acquired data. This work presents a comprehensive, novel software architecture for small-animal CT scanners based on cone-beam geometry with circular scanning trajectory. The proposed architecture covers all the steps from the system calibration to the volume reconstruction and conversion into Hounsfield units. It includes an efficient implementation of an FDK-based reconstruction algorithm that takes advantage of system symmetries and allows for parallel reconstruction using a multiprocessor computer. Strategies for calibration and artifact correction are discussed to justify the strategies adopted. New procedures for multi-bed misalignment, beam-hardening, and Housfield units calibration are proposed. Experiments with phantoms and real data showed the suitability of the proposed software architecture for an X-ray small animal CT based on cone-beam geometry.

Statistical projection completion in X-ray CT using consistency conditions

Projection data incompleteness arises in many situations relevant to X-ray computed tomography (CT) imaging. We propose a penalized maximum likelihood statistical sinogram restoration approach that incorporates the Helgason-Ludwig (HL) consistency conditions to accommodate projection data incompleteness. Image reconstruction is performed by the filtered-backprojection (FBP) in a second step. In our problem formulation, the objective function consists of the log-likelihood of the X-ray CT data and a penalty term; the HL condition poses a linear constraint on the restored sinogram and can be implemented efficiently via fast Fourier transform (FFT) and inverse FFT. We derive an iterative algorithm that increases the objective function monotonically. The proposed algorithm is applied to both computer simulated data and real patient data. We study different factors in the problem formulation that affect the properties of the final FBP reconstructed images, including the data truncation level, the amount of prior knowledge on the object support, as well as different approximations of the statistical distribution of the available projection data. We also compare its performance with an analytical truncation artifacts reduction method. The proposed method greatly improves both the accuracy and the precision of the reconstructed images within the scan field-of-view, and to a certain extent recovers the truncated peripheral region of the object. The proposed method may also be applied in areas such as limited angle tomography, metal artifacts reduction, and sparse sampling imaging.

Fourier properties of the fan-beam sinogram

PURPOSE

P. R. Edholm, R. M. Lewitt, and B. Lindholm, "Novel properties of the Fourier decomposition of the sinogram," in Proceedings of the International Workshop on Physics and Engineering of Computerized Multidimensional Imaging and Processing [Proc. SPIE 671, 8-18 (1986)] described properties of a parallel beam projection sinogram with respect to its radial and angular frequencies. The purpose is to perform a similar derivation to arrive at corresponding properties of a fan-beam projection sinogram for both the equal-angle and equal-spaced detector sampling scenarios.

METHODS

One of the derived properties is an approximately zero-energy region in the two-dimensional Fourier transform of the full fan-beam sinogram. This region is in the form of a double-wedge, similar to the parallel beam case, but different in that it is asymmetric with respect to the frequency axes. The authors characterize this region for a point object and validate the derived properties in both a simulation and a head CT data set. The authors apply these results in an application using algebraic reconstruction.

RESULTS

In the equal-angle case, the domain of the zero region is (q,k) for which / k/(k-q) / > R/L, where q and k are the frequency variables associated with the detector and view angular positions, respectively, R is the radial support of the object, and L is the source-to-isocenter distance. A filter was designed to retain only sinogram frequencies corresponding to a specified radial support. The filtered sinogram was used to reconstruct the same radial support of the head CT data. As an example application of this concept, the double-wedge filter was used to computationally improve region of interest iterative reconstruction.

CONCLUSIONS

Interesting properties of the fan-beam sinogram exist and may be exploited in some applications.

Approximate and exact cone-beam reconstruction with standard and non-standard spiral scanning

The long object problem is practically important and theoretically challenging. To solve the long object problem, spiral cone-beam CT was first proposed in 1991, and has been extensively studied since then. As a main feature of the next generation medical CT, spiral cone-beam CT has been greatly improved over the past several years, especially in terms of exact image reconstruction methods. Now, it is well established that volumetric images can be exactly and efficiently reconstructed from longitudinally truncated data collected along a rather general scanning trajectory. Here we present an overview of some key results in this area.

Statistical reconstruction for x-ray CT systems with non-continuous detectors

We analyse the performance of statistical reconstruction (SR) methods when applied to non-continuous x-ray detectors. Robustness to projection gaps is required in x-ray CT systems with multiple detector modules or with defective detector pixels. In such situations, the advantage of statistical reconstruction is that it is able to ignore missing or faulty pixels and that it makes optimal use of the remaining line integrals. This potentially obviates the need to fill the sinogram discontinuities by interpolation or any other approximative pre-processing techniques. In this paper, we apply SR to cone beam projections of (i) a hypothetical modular detector micro-CT scanner and of (ii) a system with randomly located defective detector elements. For the modular-detector system, SR produces reconstruction volumes free of noticeable gap-induced artefacts as long as the location of detector gaps and selection of the scanning range provide complete object sampling in the central imaging plane. When applied to randomly located faulty detector elements, SR produces images free of substantial ring artefacts even for cases where defective pixels cover as much as 3% of the detector area.

Data consistency based translational motion artifact reduction in fan-beam CT

A basic assumption in the classic computed tomography (CT) theory is that an object remains stationary in an entire scan. In biomedical CT/micro-CT, this assumption is often violated. To produce high-resolution images, such as for our recently proposed clinical micro-CT (CMCT) prototype, it is desirable to develop a precise motion estimation and image reconstruction scheme. In this paper, we first extend the Helgason-Ludwig consistency condition (HLCC) from parallel-beam to fan-beam geometry when an object is subject to a translation. Then, we propose a novel method to estimate the motion parameters only from sinograms based on the HLCC. To reconstruct the moving object, we formulate two generalized fan-beam reconstruction methods, which are in filtered backprojection and backprojection filtering formats, respectively. Furthermore, we present numerical simulation results to show that our approach is accurate and robust.