A new framework of image reconstruction from fan beam projections

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.

Image reconstruction from truncated fan-beam projections along a circular trajectory using the virtual fan-beam method

Objective.In this paper, we investigate how the virtual fan-beam (VFB) method can be used to perform mathematically correct 2D reconstruction in a region-of-interest (ROI), using truncated fan-beam projections acquired on a circular scan, for truncation that only occurs on one side of the object.Approach.We start by choosing a virtual fan-beam trajectory and specifying how to obtain the corresponding virtual projections. Then, three VFB formulas are obtained by applying known super-short-scan (SSS) formulas to this virtual trajectory. Two of them perform the backprojection in a virtual parallel geometry and the third in the virtual fan-beam geometry. Next, we develop two VFB formulas that perform the backprojection step in the fan-beam acquisition geometry.Main results.We present five VFB reconstruction formulas for this truncation setting. To our knowledge, the two VFB formulas performing the backprojection in the fan-beam acquisition geometry are new. Moreover, the five VFB formulas presented here obtain accurate reconstruction in a larger ROI than what has been previously reported in the literature in the same setting. A complete mathematical derivation of these five VFB formulas is given, and their implementation is described step by step. Numerical simulations, using the Forbild head and thorax phantoms, demonstrate the efficacy of these formulas. A spatial resolution analysis and a variance study indicate minor differences between these five VFB formulas.Significance.This work shows that many different VFB formulas can be applied to perform mathematically correct 2D reconstruction in a ROI, in case of truncated fan-beam projections acquired on a circular scan. Moreover, the two new VFB formulas, with backprojection in the acquisition geometry, may open the path for an extension of the VFB method to 3D reconstruction from transversely truncated cone-beam projection acquired on a circular scan.

A new weighting scheme for arc based circle cone-beam CT reconstruction

In this paper, we present an arc based fan-beam computed tomography (CT) reconstruction algorithm by applying Katsevich's helical CT image reconstruction formula to 2D fan-beam CT scanning data. Specifically, we propose a new weighting function to deal with the redundant data. Our weighting function ϖ(x_,λ) is an average of two characteristic functions, where each characteristic function indicates whether the projection data of the scanning angle contributes to the intensity of the pixel x_. In fact, for every pixel x_, our method uses the projection data of two scanning angle intervals to reconstruct its intensity, where one interval contains the starting angle and another contains the end angle. Each interval corresponds to a characteristic function. By extending the fan-beam algorithm to the circle cone-beam geometry, we also obtain a new circle cone-beam CT reconstruction algorithm. To verify the effectiveness of our method, the simulated experiments are performed for 2D fan-beam geometry with straight line detectors and 3D circle cone-beam geometry with flat-plan detectors, where the simulated sinograms are generated by the open-source software "ASTRA toolbox." We compare our method with the other existing algorithms. Our experimental results show that our new method yields the lowest root-mean-square-error (RMSE) and the highest structural-similarity (SSIM) for both reconstructed 2D and 3D fan-beam CT images.

Learning to Reconstruct Computed Tomography Images Directly From Sinogram Data Under A Variety of Data Acquisition Conditions

Computed tomography (CT) is widely used in medical diagnosis and non-destructive detection. Image reconstruction in CT aims to accurately recover pixel values from measured line integrals, i.e., the summed pixel values along straight lines. Provided that the acquired data satisfy the data sufficiency condition as well as other conditions regarding the view angle sampling interval and the severity of transverse data truncation, researchers have discovered many solutions to accurately reconstruct the image. However, if these conditions are violated, accurate image reconstruction from line integrals remains an intellectual challenge. In this paper, a deep learning method with a common network architecture, termed iCT-Net, was developed and trained to accurately reconstruct images for previously solved and unsolved CT reconstruction problems with high quantitative accuracy. Particularly, accurate reconstructions were achieved for the case when the sparse view reconstruction problem (i.e., compressed sensing problem) is entangled with the classical interior tomographic problems.

On the data acquisition, image reconstruction, cone beam artifacts, and their suppression in axial MDCT and CBCT - A review

PURPOSE

In the clinic, computed tomography (CT) has evolved into an essential modality for diagnostic imaging by multidetector row CT (MDCT) and image guided intervention by cone beam CT (CBCT). Recognizing the increasing importance of axial MDCT/CBCT in clinical and preclinical applications, and the existence of CB artifacts in MDCT/CBCT images, we provide a review of CB artifacts' root causes, rendering mechanisms and morphology, and possible solutions for elimination and/or reduction of the artifacts.

METHODS

By examining the null space in Radon and Fourier domain, the root cause of CB artifacts (i.e., data insufficiency) in axial MDCT/CBCT is analytically investigated, followed by a review of the data sufficiency conditions and the "circle +" source trajectories. The rendering mechanisms and morphology of CB artifacts in axial MDCT/CBCT and their special cases (e.g., half/short scan and full scan with latitudinally displaced detector) are then analyzed, followed by a survey of the potential solutions to suppress the artifacts. The phenomenon of imaged zone indention and its variation over FBP, BPF/DBPF, two-pass and iterative CB reconstruction algorithms and/or schemes are discussed in detail.

RESULTS

An interdomain examination of the null space provides an insightful understanding of the root cause of CB artifacts in axial MDCT/CBCT. The decomposition of CB artifacts rendering mechanisms facilitates understanding of the artifacts' behavior under different conditions and the potential solutions to suppress them. An inspection of the imaged zone intention phenomenon provides guidance on the design and implementation of CB image reconstruction algorithms and schemes for CB artifacts suppression in axial MDCT/CBCT.

CONCLUSIONS

With increasing importance of axial MDCT/CBCT in clinical and preclinical applications, this review article can update the community with in-depth information and clarification on the latest progress in dealing with CB artifacts and thus increase clinical/preclinical confidence.

Synchronized multiartifact reduction with tomographic reconstruction (SMART-RECON): A statistical model based iterative image reconstruction method to eliminate limited-view artifacts and to mitigate the temporal-average artifacts in time-resolved CT

PURPOSE

In x-ray computed tomography (CT), a violation of the Tuy data sufficiency condition leads to limited-view artifacts. In some applications, it is desirable to use data corresponding to a narrow temporal window to reconstruct images with reduced temporal-average artifacts. However, the need to reduce temporal-average artifacts in practice may result in a violation of the Tuy condition and thus undesirable limited-view artifacts. In this paper, the authors present a new iterative reconstruction method, synchronized multiartifact reduction with tomographic reconstruction (SMART-RECON), to eliminate limited-view artifacts using data acquired within an ultranarrow temporal window that severely violates the Tuy condition.

METHODS

In time-resolved contrast enhanced CT acquisitions, image contrast dynamically changes during data acquisition. Each image reconstructed from data acquired in a given temporal window represents one time frame and can be denoted as an image vector. Conventionally, each individual time frame is reconstructed independently. In this paper, all image frames are grouped into a spatial-temporal image matrix and are reconstructed together. Rather than the spatial and/or temporal smoothing regularizers commonly used in iterative image reconstruction, the nuclear norm of the spatial-temporal image matrix is used in SMART-RECON to regularize the reconstruction of all image time frames. This regularizer exploits the low-dimensional structure of the spatial-temporal image matrix to mitigate limited-view artifacts when an ultranarrow temporal window is desired in some applications to reduce temporal-average artifacts. Both numerical simulations in two dimensional image slices with known ground truth and in vivo human subject data acquired in a contrast enhanced cone beam CT exam have been used to validate the proposed SMART-RECON algorithm and to demonstrate the initial performance of the algorithm. Reconstruction errors and temporal fidelity of the reconstructed images were quantified using the relative root mean square error (rRMSE) and the universal quality index (UQI) in numerical simulations. The performance of the SMART-RECON algorithm was compared with that of the prior image constrained compressed sensing (PICCS) reconstruction quantitatively in simulations and qualitatively in human subject exam.

RESULTS

In numerical simulations, the 240(∘) short scan angular span was divided into four consecutive 60(∘) angular subsectors. SMART-RECON enables four high temporal fidelity images without limited-view artifacts. The average rRMSE is 16% and UQIs are 0.96 and 0.95 for the two local regions of interest, respectively. In contrast, the corresponding average rRMSE and UQIs are 25%, 0.78, and 0.81, respectively, for the PICCS reconstruction. Note that only one filtered backprojection image can be reconstructed from the same data set with an average rRMSE and UQIs are 45%, 0.71, and 0.79, respectively, to benchmark reconstruction accuracies. For in vivo contrast enhanced cone beam CT data acquired from a short scan angular span of 200(∘), three 66(∘) angular subsectors were used in SMART-RECON. The results demonstrated clear contrast difference in three SMART-RECON reconstructed image volumes without limited-view artifacts. In contrast, for the same angular sectors, PICCS cannot reconstruct images without limited-view artifacts and with clear contrast difference in three reconstructed image volumes.

CONCLUSIONS

In time-resolved CT, the proposed SMART-RECON method provides a new method to eliminate limited-view artifacts using data acquired in an ultranarrow temporal window, which corresponds to approximately 60(∘) angular subsectors.

Time-resolved cardiac interventional cone-beam CT reconstruction from fully truncated projections using the prior image constrained compressed sensing (PICCS) algorithm

C-arm cone-beam CT could replace preoperative multi-detector CT scans in the cardiac interventional setting. However, cardiac gating results in view angle undersampling and the small size of the detector results in projection data truncation. These problems are incompatible with conventional tomographic reconstruction algorithms. In this paper, the prior image constrained compressed sensing (PICCS) reconstruction method was adapted to solve these issues. The performance of the proposed method was compared to that of FDK, FDK with extrapolated projection data (E-FDK), and total variation-based compressed sensing. A canine projection dataset acquired using a clinical C-arm imaging system supplied realistic cardiac motion and anatomy for this evaluation. Three different levels of truncation were simulated. The relative root mean squared error and the universal image quality index were used to quantify the reconstruction accuracy. Three main conclusions were reached. (1) The adapted version of the PICCS algorithm offered the highest image quality and reconstruction accuracy. (2) No meaningful variation in performance was observed when the amount of truncation was changed. (3) This study showed evidence that accurate interior tomography with an undersampled acquisition is possible for realistic objects if a prior image with minimal artifacts is available.

Region-of-interest image reconstruction in circular cone-beam microCT

Cone-beam microcomputed tomography (microCT) is one of the most popular choices for small animal imaging which is becoming an important tool for studying animal models with transplanted diseases. Region-of-interest (ROI) imaging techniques in CT, which can reconstruct an ROI image from the projection data set of the ROI, can be used not only for reducing imaging-radiation exposure to the subject and scatters to the detector but also for potentially increasing spatial resolution of the reconstructed images. Increasing spatial resolution in microCT images can facilitate improved accuracy in many assessment tasks. A method proposed previously for increasing CT image spatial resolution entails the exploitation of the geometric magnification in cone-beam CT. Due to finite detector size, however, this method can lead to data truncation for a large geometric magnification. The Feldkamp-Davis-Kress (FDK) algorithm yields images with artifacts when truncated data are used, whereas the recently developed backprojection filtration (BPF) algorithm is capable of reconstructing ROI images without truncation artifacts from truncated cone-beam data. We apply the BPF algorithm to reconstructing ROI images from truncated data of three different objects acquired by our circular cone-beam microCT system. Reconstructed images by use of the FDK and BPF algorithms from both truncated and nontruncated cone-beam data are compared. The results of the experimental studies demonstrate that, from certain truncated data, the BPF algorithm can reconstruct ROI images with quality comparable to that reconstructed from nontruncated data. In contrast, the FDK algorithm yields ROI images with truncation artifacts. Therefore, an implication of the studies is that, when truncated data are acquired with a configuration of a large geometric magnification, the BPF algorithm can be used for effective enhancement of the spatial resolution of a ROI image.

Temporally targeted imaging method applied to ECG-gated computed tomography: preliminary phantom and in vivo experience

RATIONALE AND OBJECTIVES

Existing cardiac imaging methods do not allow for improved temporal resolution when considering a targeted region of interest (ROI). The imaging method presented here enables improved temporal resolution for ROI imaging (namely, a reconstruction volume smaller than the complete field of view). Clinically, temporally targeted reconstruction would not change the primary means of reconstructing and evaluating images, but rather would enable the adjunct technique of ROI imaging, with improved temporal resolution compared with standard reconstruction ( approximately 20% smaller temporal scan window). In gated cardiac computed tomography (CT) scans improved temporal resolution directly translates into a reduction in motion artifacts for rapidly moving objects such as the coronary arteries.

MATERIALS AND METHODS

Retrospectively electrocardiogram gated coronary angiography data from a 64-slice CT system were used. A motion phantom simulating the motion profile of a coronary artery was constructed and scanned. Additionally, an in vivo study was performed using a porcine model. Comparisons between the new reconstruction technique and the standard reconstruction are given for an ROI centered on the right coronary artery, and a pulmonary ROI.

RESULTS

In both a well-controlled motion model and a porcine model results show a decrease in motion induced artifacts including motion blur and streak artifacts from contrast enhanced vessels within the targeted ROIs, as assessed through both qualitative and quantitative observations.

CONCLUSION

Temporally targeted reconstruction techniques demonstrate the potential to reduce motion artifacts in coronary CT. Further study is warranted to demonstrate the conditions under which this technique will offer direct clinical utility. Improvement in temporal resolution for gated cardiac scans has implications for improving: contrast enhanced CT angiography, calcium scoring, and assessment of the pulmonary anatomy.

A local region of interest image reconstruction via filtered backprojection for fan-beam differential phase-contrast computed tomography

Recently, x-ray differential phase contrast computed tomography (DPC-CT) has been experimentally implemented using a conventional source combined with several gratings. Images were reconstructed using a parallel-beam reconstruction formula. However, parallel-beam reconstruction formulae are not directly applicable for a large image object where the parallel-beam approximation fails. In this note, we present a new image reconstruction formula for fan-beam DPC-CT. There are two major features in this algorithm: (1) it enables the reconstruction of a local region of interest (ROI) using data acquired from an angular interval shorter than 180 degrees + fan angle and (2) it still preserves the filtered backprojection structure. Numerical simulations have been conducted to validate the image reconstruction algorithm.

A new scheme for view-dependent data differentiation in fan-beam and cone-beam computed tomography

In computed tomography, analytical fan-beam (FB) and cone-beam (CB) image reconstruction often involves a view-dependent data differentiation. The implementation of this differentiation step is critical in terms of resolution and image quality. In this work, we present a new differentiation scheme that is robust to changes in the data acquisition geometry and to coarse view sampling. Our scheme was compared to two previously suggested methods, which we call the direct scheme and the chain-rule scheme. Image reconstructions were performed from computer-simulated data of the Shepp-Logan phantom, the FORBILD thorax phantom and a modified FORBILD head phantom. For FB reconstruction, we investigated three acquisition geometries: a circular, an ellipse-shaped and a square-shaped trajectory. For CB reconstruction, the circle-plus-line trajectory was considered. Image comparison showed that the new scheme performs consistently well when varying the scenario, in both FB and CB geometry, unlike the other two schemes.

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.

Hilbert transform based FBP algorithm for fan-beam CT full and partial scans

This paper presents a new type of filtered backprojection (FBP) algorithm for fan-beam full- and partial-scans. The filtering is shift-invariant with respect to the angular variable. The backprojection does not include position-dependent weights through the Hilbert transform and the one-dimensional transformation between the fan- and parallel-beam coordinates. The strong symmetry of the filtered projections directly leads to an exact reconstruction for partial data. The use of the Hilbert transform avoids the approximation introduced by the nonuniform cutoff frequency required in the ramp filter-based FBP algorithm. Variance analysis indicates that the algorithm might lead to a better uniformity of resolution and noise in the reconstructed image. Numerical simulations are provided to evaluate the algorithm with noise-free and noisy projections. Our simulation results indicate that the algorithm does have better stability over the ramp-filter-based FBP and circular harmonic reconstruction algorithms. This may help improve the image quality for in place computed tomography scanners with single-row detectors.

Development and evaluation of an exact fan-beam reconstruction algorithm using an equal weighting scheme via locally compensated filtered backprojection (LCFBP)

A novel exact fan-beam image reconstruction formula is presented and validated using both phantom data and clinical data. This algorithm takes the form of the standard ramp filtered backprojection (FBP) algorithm plus local compensation terms. This algorithm will be referred to as a locally compensated filtered backprojection (LCFBP). An equal weighting scheme is utilized in this algorithm in order to properly account for redundantly measured projection data. The algorithm has the desirable property of maintaining a mathematically exact result for: the full scan mode (2pi), the short scan mode (pi+ full fan angle), and the supershort scan mode [less than (pi+ full fan angle)]. Another desirable feature of this algorithm is that it is derivative-free. This feature is beneficial in preserving the spatial resolution of the reconstructed images. The third feature is that an equal weighting scheme has been utilized in the algorithm, thus the new algorithm has better noise properties than the standard filtered backprojection image reconstruction with a smooth weighting function. Both phantom data and clinical data were utilized to validate the algorithm and demonstrate the superior noise properties of the new algorithm.

General formula for fan-beam computed tomography

In this Letter, we derive a general reconstruction formula for fan-beam computed tomography (CT) utilizing the two-dimensional Dirac function, which is useful in CT imaging for moving objects.

Fan-beam and cone-beam image reconstruction via filtering the backprojection image of differentiated projection data

In this paper, a new image reconstruction scheme is presented based on Tuy's cone-beam inversion scheme and its fan-beam counterpart. It is demonstrated that Tuy's inversion scheme may be used to derive a new framework for fanbeam and cone-beam image reconstruction. In this new framework, images are reconstructed via filtering the backprojection image of differentiated projection data. The new framework is mathematically exact and is applicable to a general source trajectory provided the Tuy data sufficiency condition is satisfied. By choosing a piece-wise constant function for one of the components in the factorized weighting function, the filtering kernel is one dimensional, viz. the filtering process is along a straight line. Thus, the derived image reconstruction algorithm is mathematically exact and efficient. In the cone-beam case, the derived reconstruction algorithm is applicable to a large class of source trajectories where the pi-lines or the generalized pi-lines exist. In addition, the new reconstruction scheme survives the super-short scan mode in both the fan-beam and cone-beam cases provided the data are not transversely truncated. Numerical simulations were conducted to validate the new reconstruction scheme for the fan-beam case.

A unified framework for exact cone-beam reconstruction formulas

In this paper, we present concise proofs of several recently developed exact cone-beam reconstruction methods in the Tuy inversion framework, including both filtered-backprojection and backprojection-filtration formulas in the cases of standard spiral, nonstandard spiral, and more general scanning loci. While a similar proof of the Katsevich formula was previously reported, we present a new proof of the Zou and Pan backprojection-filtration formula. Our proof combines both odd and even data extensions so that only the cone-beam transform itself is utilized in the backprojection-filtration inversion. More importantly, our formulation is valid for general smooth scanning curves, in agreement with an earlier paper from our group [Ye, Zhao, Yu, and Wang, Proc. SPIE 5535, 293-300 (Aug. 6 2004)]. As a consequence of that proof, we obtain a new inversion formula, which is in a two-dimensional filtering backprojection format. A possibility for generalization of the Katsevich filtered-backprojection reconstruction method is also discussed from the viewpoint of this framework.

Exact fan-beam image reconstruction algorithm for truncated projection data acquired from an asymmetric half-size detector

In this paper, we present a new algorithm designed for a specific data truncation problem in fan-beam CT. We consider a scanning configuration in which the fan-beam projection data are acquired from an asymmetrically positioned half-sized detector. Namely, the asymmetric detector only covers one half of the scanning field of view. Thus, the acquired fan-beam projection data are truncated at every view angle. If an explicit data rebinning process is not invoked, this data acquisition configuration will reek havoc on many known fan-beam image reconstruction schemes including the standard filtered backprojection (FBP) algorithm and the super-short-scan FBP reconstruction algorithms. However, we demonstrate that a recently developed fan-beam image reconstruction algorithm which reconstructs an image via filtering a backprojection image of differentiated projection data (FBPD) survives the above fan-beam data truncation problem. Namely, we may exactly reconstruct the whole image object using the truncated data acquired in a full scan mode (2pi angular range). We may also exactly reconstruct a small region of interest (ROI) using the truncated projection data acquired in a short-scan mode (less than 2pi angular range). The most important characteristic of the proposed reconstruction scheme is that an explicit data rebinning process is not introduced. Numerical simulations were conducted to validate the new reconstruction algorithm.

A new data consistency condition for fan-beam projection data

The sum of all attenuation data acquired in one view of parallel-beam projections is a view angle independent constant. This fact is known as a data consistency condition on the two-dimensional Radon transforms. It plays an important role in tomographic image reconstruction and artifact correction. In this paper, a novel fan-beam data consistency condition (FDCC) is derived and presented. Using the FDCC, individual projection data in one view of fan-beam projections can be estimated from filtering all other projection data measured from different view angles. Numerical simulations are performed to validate the new FDCC in correcting ring artifacts caused by malfunctioning detector cells.

A novel extension of the parallel-beam projection-slice theorem to divergent fan-beam and cone-beam projections

The general goal of this paper is to extend the parallel-beam projection-slice theorem to divergent fan-beam and cone-beam projections without rebinning the divergent fan-beam and cone-beam projections into parallel-beam projections directly. The basic idea is to establish a novel link between the local Fourier transform of the projection data and the Fourier transform of the image object. Analogous to the two- and three-dimensional parallel-beam cases, the measured projection data are backprojected along the projection direction and then a local Fourier transform is taken for the backprojected data array. However, due to the loss of the shift invariance of the image object in a single view of the divergent-beam projections, the measured projection data is weighted by a distance dependent weight w(r) before the local Fourier transform is performed. The variable r in the weighting function w(r) is the distance from the backprojected point to the x-ray source position. It is shown that a special choice of the weighting function, w(r)=1/r, will facilitate the calculations and a simple relation can be established between the Fourier transform of the image function and the local Fourier transform of the 1/r-weighted backprojection data array. Unlike the parallel-beam cases, a one-to-one correspondence does not exist for a local Fourier transform of the backprojected data array and a single line in the two-dimensional (2D) case or a single slice in the 3D case of the Fourier transform of the image function. However, the Fourier space of the image object can be built up after the local Fourier transforms of the 1/r-weighted backprojection data arrays are shifted and then summed in a laboratory frame. Thus the established relations Eq. (27) and Eq. (29) between the Fourier space of the image object and the Fourier transforms of the backprojected data arrays can be viewed as a generalized projection-slice theorem for divergent fan-beam and cone-beam projections. Once the Fourier space of the image function is built up, an inverse Fourier transform could be performed to reconstruct tomographic images from the divergent beam projections. Due to the linearity of the Fourier transform, an image reconstruction step can be performed either when the complete Fourier space is available or in parallel with the building of the Fourier space. Numerical simulations are performed to verify the generalized projection-slice theorem by using a disc phantom in the fan-beam case.

A two-step Hilbert transform method for 2D image reconstruction

The paper describes a new accurate two-dimensional (2D) image reconstruction method consisting of two steps. In the first step, the backprojected image is formed after taking the derivative of the parallel projection data. In the second step, a Hilbert filtering is applied along certain lines in the differentiated backprojection (DBP) image. Formulae for performing the DBP step in fanbeam geometry are also presented. The advantage of this two-step Hilbert transform approach is that in certain situations, regions of interest (ROIs) can be reconstructed from truncated projection data. Simulation results are presented that illustrate very similar reconstructed image quality using the new method compared to standard filtered backprojection, and that show the capability to correctly handle truncated projections. In particular, a simulation is presented of a wide patient whose projections are truncated laterally yet for which highly accurate ROI reconstruction is obtained.

An alternative derivation of Katsevich's cone-beam reconstruction formula

In this paper an alternative derivation of Katsevich's cone-beam image reconstruction algorithm is presented. The starting point is the classical Tuy's inversion formula. After (i) using the hidden symmetries of the intermediate functions, (ii) handling the redundant data by weighting them, (iii) changing the weighted average into an integral over the source trajectory parameter, and (iv) imposing an additional constraint on the weighting function, a filtered backprojection reconstruction formula from cone beam projections is derived. The following features are emphasized in the present paper: First, the nontangential condition in Tuy's original data sufficiency conditions has been relaxed. Second, a practical regularization scheme to handle the singularity is proposed. Third, the derivation in the cone beam case is in the same fashion as that in the fan-beam case. Our final cone-beam reconstruction formula is the same as the one discovered by Katsevich in his most recent paper. However, the data sufficiency conditions and the regularization scheme of singularities are different. A detailed comparison between these two methods is presented.

Nonuniform noise propagation by using the ramp filter in fan-beam computed tomography

It is observed that when the homogeneity property of the ramp filter is used to derive a filtered backprojection algorithm in fan-beam tomography, the reconstructed images have nonstationary frequency components and nonstationary noise. When a short focal-length is used, higher frequency components are amplified more at the edge of the image than at the center of the image, resulting in higher noise at the edge of the image.