Exact BPF and FBP algorithms for nonstandard saddle curves

Completeness map evaluation demonstrated with candidate next-generation cardiac CT architectures

PURPOSE

In this report, the authors introduce the general concept of the completeness map, as a means to evaluate the completeness of data acquired by a given CT system design (architecture and scan mode). They illustrate the utility of completeness map by applying the completeness map concept to a number of candidate CT system designs, as part of a study to advance the state-of-the-art in cardiac CT.

METHODS

In order to optimally reconstruct a point within a volume of interest (VOI), the Radon transform on all possible planes through that point should be measured. The authors quantified the extent to which this ideal condition is satisfied for the entire image volume. They first determined a Radon completeness number for each point in the VOI, as the percentage of possible planes that is actually measured. A completeness map is then defined as a 3D matrix of the completeness numbers for the entire VOI. The authors proposed algorithms to analyze the projection datasets in Radon space and compute the completeness number for a fixed point and apply these algorithms to various architectures and scan modes that they are evaluating. In this report, the authors consider four selected candidate architectures, operating with different scan modes, for a total of five system design alternatives. Each of these alternatives is evaluated using completeness map.

RESULTS

If the detector size and cone angle are large enough to cover the entire cardiac VOI, a single-source circular scan can have ≥99% completeness over the entire VOI. However, only the central z-slice can be exactly reconstructed, which corresponds to 100% completeness. For a typical single-source architecture, if the detector is limited to an axial dimension of 40 mm, a helical scan needs about five rotations to form an exact reconstruction region covering the cardiac VOI, while a triple-source helical scan only requires two rotations, leading to a 2.5x improvement in temporal resolution. If the source and detector of an inverse-geometry (IGCT) system have the same axial extent, and the spacing of source points in the axial and transaxial directions is sufficiently small, the IGCT can also form an exact reconstruction region for the cardiac VOI. If the VOI can be covered by the x-ray beam in any view, a composite-circling scan can generate an exact reconstruction region covering the VOI.

CONCLUSIONS

The completeness map evaluation provides useful information for selecting the next-generation cardiac CT system design. The proposed completeness map method provides a practical tool for analyzing complex scanning trajectories, where the theoretical image quality for some complex system designs is impossible to predict, without yet-undeveloped reconstruction algorithms.

Cone-Beam CT with a Flat-Panel Detector: From Image Science to Image-Guided Surgery

The development of large-area flat-panel x-ray detectors (FPDs) has spurred investigation in a spectrum of advanced medical imaging applications, including tomosynthesis and cone-beam CT (CBCT). Recent research has extended image quality metrics and theoretical models to such applications, providing a quantitative foundation for the assessment of imaging performance as well as a general framework for the design, optimization, and translation of such technologies to new applications. For example, cascaded systems models of Fourier domain metrics, such as noise-equivalent quanta (NEQ), have been extended to these modalities to describe the propagation of signal and noise through the image acquisition and reconstruction chain and to quantify the factors that govern spatial resolution, image noise, and detectability. Moreover, such models have demonstrated basic agreement with human observer performance for a broad range of imaging conditions and imaging tasks. These developments in image science have formed a foundation for the knowledgeable development and translation of CBCT to new applications in image-guided interventions - for example, CBCT implemented on a mobile surgical C-arm for intraoperative 3D imaging. The ability to acquire high-quality 3D images on demand during surgical intervention overcomes conventional limitations of surgical guidance in the context of preoperative images alone. A prototype mobile C-arm developed in academic-industry partnership demonstrates CBCT with low radiation dose, sub-mm spatial resolution, and soft-tissue visibility potentially approaching that of diagnostic CT. Integration of the 3D imaging system with real-time tracking, deformable registration, endoscopic video, and 3D visualization offers a promising addition to the surgical arsenal in interventions ranging from head-and-neck / skull base surgery to spine, orthopaedic, thoracic, and abdominal surgeries. Cadaver studies show the potential for significant boosts in surgical performance under CBCT guidance, and early clinical trials demonstrate feasibility, workflow, and image quality within the surgical theatre.

Cone-beam composite-circling scan and exact image reconstruction for a quasi-short object

Here we propose a cone-beam composite-circling mode to solve the quasi-short object problem, which is to reconstruct a short portion of a long object from longitudinally truncated cone-beam data involving the short object. In contrast to the saddle curve cone-beam scanning, the proposed scanning mode requires that the X-ray focal spot undergoes a circular motion in a plane facing the short object, while the X-ray source is rotated in the gantry main plane. Because of the symmetry of the proposed mechanical rotations and the compatibility with the physiological conditions, this new mode has significant advantages over the saddle curve from perspectives of both engineering implementation and clinical applications. As a feasibility study, a backprojection filtration (BPF) algorithm is developed to reconstruct images from data collected along a composite-circling trajectory. The initial simulation results demonstrate the correctness of the proposed exact reconstruction method and the merits of the proposed mode.

Exact interior reconstruction from truncated limited-angle projection data

Using filtered backprojection (FBP) and an analytic continuation approach, we prove that exact interior reconstruction is possible and unique from truncated limited-angle projection data, if we assume a prior knowledge on a subregion or subvolume within an object to be reconstructed. Our results show that (i) the interior region-of-interest (ROI) problem and interior volume-of-interest (VOI) problem can be exactly reconstructed from a limited-angle scan of the ROI/VOI and a 180 degree PI-scan of the subregion or subvolume and (ii) the whole object function can be exactly reconstructed from nontruncated projections from a limited-angle scan. These results improve the classical theory of Hamaker et al. (1980).

Knowledge-based dynamic volumetric cardiac computed tomography with saddle curve trajectory

Motion artifact is still a major issue in cardiac computed tomography because the current motion correction and electrocardiogram gating techniques have not fully addressed this problem. The image quality can be significantly improved by using information about the actual state of the heart and an exact reconstruction algorithm. We propose to extend a cardiac computed tomographic technique, using the knowledge of the volume and the relation between the state and the phase of the heart, to a saddle curve trajectory. This will optimize the image quality by reducing the artifacts resulting from approximate reconstruction and solve the long-object problem. Necessary background is provided, and the effectiveness of the algorithms is demonstrated in numerical simulations with the dynamic thorax phantom.

Anniversary paper. Development of x-ray computed tomography: the role of medical physics and AAPM from the 1970s to present

The AAPM, through its members, meetings, and its flagship journal Medical Physics, has played an important role in the development and growth of x-ray tomography in the last 50 years. From a spate of early articles in the 1970s characterizing the first commercial computed tomography (CT) scanners through the "slice wars" of the 1990s and 2000s, the history of CT and related techniques such as tomosynthesis can readily be traced through the pages of Medical Physics and the annals of the AAPM and RSNA/AAPM Annual Meetings. In this article, the authors intend to give a brief review of the role of Medical Physics and the AAPM in CT and tomosynthesis imaging over the last few decades.

Lambda tomography with discontinuous scanning trajectories

Lambda tomography (LT) is a well-known local reconstruction technology to reduce the radiation dose or accommodate a limited imaging geometry. After a theoretical analysis of the so-called Calderon operator (CO), the necessary conditions for exact LT reconstruction are presented in terms of the 2D and 3D COs. Based on our previous results on LT, a general scheme is proposed to construct exact LT formulae in terms of the 2D CO with multiple segment trajectories. Every 2D formula has a corresponding 3D cone-beam formula in the Feldkamp framework in terms of the 2D CO which was illustrated in a triple-segment case. Our simulation results verify the correctness and demonstrate the merits of the proposed scheme.

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.

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

It is well known that a rigid in-plane motion can be decomposed into a translation and a rotation around an origin. Based on our previous work, we first extend the Helgason-Ludwig consistency condition (HLCC) to cover a general rigid motion in fan-beam geometry. Then, we model the general motion by several parameters, and develop an iterative scheme for estimation of the in-plane motion parameters. This scheme determines the motion parameters by numerically minimizing an objective function constructed based on the HLCC. After the motion parameters are estimated, image reconstruction can be performed to compensate for the motion effects. Finally, we implement the algorithm and evaluate its performance in numerical simulations.

Practical cone-beam lambda tomography

As a potentially important technology for medical x-ray computed tomography (CT), lambda tomography (LT) is to reconstruct a gradient-like image only from local projection data. Based on our recently derived exact fan-beam LT formula, [H. Y. Gu and G. Wang, Int. J. Biomed. Imaging 2006(1), 1-9 (2006)] here we propose a practical cone-beam LT algorithm for LT reconstruction from local data collected along an arbitrary smooth three-dimensional curve. A key step in our algorithm is to determine an appropriate vector perpendicular to the line connecting the x-ray source and an image point. The algorithm is implemented assuming an equispatial planar detector and a nonstandard spiral trajectory. The numerical simulation results demonstrate the merits of our method.

Cone-beam mammo-computed tomography from data along two tilting arcs

Over the past several years there has been an increasing interest in cone-beam computed tomography (CT) for breast imaging. In this article, we propose a new scheme for theoretically exact cone-beam mammo-CT and develop a corresponding Katsevich-type reconstruction algorithm. In our scheme, cone-beam scans are performed along two tilting arcs to collect a sufficient amount of information for exact reconstruction. In our algorithm, cone-beam data are filtered in a shift-invariant fashion and then weighted backprojected into the three-dimensional space for the final reconstruction. Our approach has several desirable features, including tolerance of axial data truncation, efficiency in sequential/parallel implementation, and accuracy for quantitative analysis. We also demonstrate the system performance and clinical utility of the proposed technique in numerical simulations.

View-independent reconstruction algorithms for cone beam CT with general saddle trajectory

In Yang et al (2006 Phys. Med. Biol. 51 1157-72), an exact filtered backprojection (FBP) reconstruction algorithm was proposed for cone beam tomography with saddle trajectory based on the seminal works of Pack and Noo (2005a Inverse Problems 21 1105-20; 2005b 8th Int. Meeting on Fully 3D Reconstruction in Radiology and Nuclear Medicine (Salt Lake City) ed F Noo, H Kudo and L G Zeng pp 287-90). However, the artefacts due to discretization and/or sampling errors in the reconstructed images by this method were still visible, especially when the pitch is large. In this paper, two view-independent (VI) algorithms, which are similar to the FDK-type algorithms (Feldkamp et al 1984 J. Opt. Soc. Am. A 1 612-19), are proposed for planar detector geometry. The first VI algorithm involves two filtered projections and a small additional term (two-dimensional (2D) Radon transform term). One of the filtered projections is obtained by ramp filtering (as in the FDK algorithm for circular trajectory) and the other one is obtained by Hilbert transform. The 2D Radon transform term is just like the term which was first derived by Hu (1996 Scanning 18 572-81) for a circular trajectory. The second VI algorithm involves only one filtered projection term, which is obtained by differentiation followed by Hilbert transform and the 2D Radon transform term. Both algorithms involve only one backprojection step with a weighting factor as in the FDK algorithm. The simulation studies show that the pixel values of the reconstructed images by the VI algorithms are more accurate than those by the original view differencing (VD) algorithm, the streak artefacts are also reduced, and their computational times are comparable to that of the original VD algorithm. We also generalize the concept of saddle trajectory and the corresponding reconstruction algorithm. The generalized algorithm is also theoretically exact, has a shift-invariant FBP structure, and does not depend on the concept of pi-line.

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.

A differentiable Shepp–Logan phantom and its applications in exact cone-beam CT

Recently, several exact cone-beam reconstruction algorithms, such as the generalized filtered-backprojection (FBP) and backprojection-filtration (BPF) methods, have been developed to solve the long object problem. Although the well-known 3D Shepp-Logan phantom (SLP) is often used to validate these algorithms, it is deficient due to the discontinuity of the SLP. In this paper, we first construct a differentiable polynomial function to approximate the unit rectangular function on [-1, 1]. Then, we use this function to obtain a differentiable ellipsoid phantom, whose x-ray transform is differentiable for any smooth scanning trajectory. Finally, we propose a differentiable Shepp-Logan phantom (DSLP) for numerical simulation of the exact cone-beam CT algorithms. Our numerical simulation shows that the reconstructed DSLP has a better image quality than the reconstructed SLP, and is complementary to the traditional SLP for evaluation of the exact cone-beam CT algorithms.