Two-dimensional iterative region-of-interest (ROI) reconstruction from truncated projection data

Comparison of iterative reconstruction implementations for multislice helical CT

The most mature image reconstruction algorithms in multislice helical computed tomography are based on analytical and iterative methods. Over the past decades, several methods have been developed for iterative reconstructions that improve image quality by reducing noise and artifacts. In the regularization step of iterative reconstruction, noise can be significantly reduced, thereby making low-dose CT. The quality of the reconstructed image can be further improved by using model-based reconstructions. In these reconstructions, the main focus is on modeling the data acquisition process, including the behavior of the photon beams, the geometry of the system, etc. In this article, we propose two model-based reconstruction algorithms using a virtual detector for multislice helical CT. The aim of this study is to compare the effect of using a virtual detector on image quality for the two proposed algorithms with a model-based iterative reconstruction using the original detector model. Since the algorithms are implemented using multiple GPUs, the merging of separately reconstructed volumes can significantly affect image quality. This issue is often referred to as the "long object" problem, for which we also present a solution that plays an important role in the proposed reconstruction processes. The algorithms were evaluated using mathematical and physical phantoms, as well as patient cases. The SSIM, MS-SSIM and L1 metrics were utilized to evaluate the image quality of the mathematical phantom case. To demonstrate the effectiveness of the algorithms, we used the CatPhan 600 phantom. Additionally, anonymized patient scans were used to showcase the improvements in image quality on real scan data.

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.

Towards high-resolution synchrotron radiation imaging with statistical iterative reconstruction

Synchrotron radiation (SR) X-ray micro-computed tomography (CT) is an effective imaging modality for high-resolution investigation of small objects, with several applications in medicine, biology and industry. However, the limited size of the detector field of view (FOV) restricts the sample dimensions to only a few millimeters. When the sample size is larger than the FOV, images reconstructed using conventional methods suffer from DC-shift and low-frequency artifacts. This classical problem is known as the local tomography or the interior problem. In this paper, a statistical iterative reconstruction method is introduced to eliminate image artifacts resulting from the local tomography. The proposed method, which can be used in several SR imaging applications, enables high-resolution SR imaging with superior image quality compared with conventional methods. Real data obtained from different SR micro-CT applications are used to evaluate the proposed method. Results indicate a noteworthy quality improvement in the image reconstructed from the local tomography measurements.

Relationship between x-ray illumination field size and flat field intensity and its impacts on x-ray imaging

PURPOSE

X-ray cone-beam CT (CBCT) is being increasingly used for various clinical applications, while its performance is still hindered by image artifacts. This work investigates a new source of reconstruction error, which is often overlooked in the current CBCT imaging. The authors find that the x-ray flat field intensity (I(0)) varies significantly as the illumination volume size changes at different collimator settings. A wrong I(0) value leads to inaccurate CT numbers of reconstructed images as well as wrong scatter measurements in the CBCT research.

METHODS

The authors argue that the finite size of x-ray focal spot together with the detector glare effect cause the I(0) variation at different illumination sizes. Although the focal spot of commercial x-ray tubes typically has a nominal size of less than 1 mm, the off-focal-spot radiation covers an area of several millimeters on the tungsten target. Due to the large magnification factor from the field collimator to the detector, the penumbra effects of the collimator blades result in different I(0) values for different illumination field sizes. Detector glare further increases the variation, since one pencil beam of incident x-ray is scattered into an area of several centimeters on the detector. In this paper, the authors study these two effects by measuring the focal spot distribution with a pinhole assembly and the detector point spread function (PSF) with an edge-spread function method. The authors then derive a formula to estimate the I(0) value for different illumination field sizes, using the measured focal spot distribution and the detector PSF. Phantom studies are carried out to investigate the accuracy of scatter measurements and CT images with and without considering the I(0) variation effects.

RESULTS

On our tabletop system with a Varian Paxscan 4030CB flat-panel detector and a Varian RAD-94 x-ray tube as used on a clinical CBCT system, the focal spot distribution has a measured full-width-at-half-maximum (FWHM) of around 0.4 mm, while non-negligible off-focal-spot radiation is observed at a distance of over 2 mm from the center. The measured detector PSF has an FWHM of 0.510 mm, with a shape close to Gaussian. From these two distributions, the author calculate the estimated I(0) values at different collimator settings. The I(0) variation mainly comes from the focal spot effect. The estimation matches well with the measurements at different collimator widths in both horizontal and vertical directions, with an average error of less than 3%. Our method improves the accuracy of conventional scatter measurements, where the scatter is measured as the difference between fan-beam and cone-beam projections. On a uniform water cylinder phantom, more accurate I(0) suppresses the unfaithful high-frequency signals at the object boundaries of the measured scatter, and the SPR estimation error is reduced from 0.158 to 0.014. The proposed I(0) estimation also reduces the reconstruction error from about 20 HU on the Catphan©600 phantom in the selected regions of interest to less than 4 HU.

CONCLUSIONS

The I(0) variation is identified as one additional error source in x-ray imaging. By measuring the focal-spot distribution and detector PSF, the authors propose an accurate method of estimating the I(0) value for different illumination field sizes. The method obtains more accurate scatter measurements and therefore facilitates scatter correction algorithm designs. As correction methods for other CBCT artifacts become more successful, our research is significant in further improving the CBCT imaging accuracy.

Comparison of extended field-of-view reconstructions in C-arm flat-detector CT using patient size, shape or attenuation information

In C-arm-based flat-detector computed tomography (FDCT) it frequently happens that the patient exceeds the scan field of view (SFOV) in the transaxial direction because of the limited detector size. This results in data truncation and CT image artefacts. In this work three truncation correction approaches for extended field-of-view (EFOV) reconstructions have been implemented and evaluated. An FDCT-based method estimates the patient size and shape from the truncated projections by fitting an elliptical model to the raw data in order to apply an extrapolation. In a camera-based approach the patient is sampled with an optical tracking system and this information is used to apply an extrapolation. In a CT-based method the projections are completed by artificial projection data obtained from the CT data acquired in an earlier exam. For all methods the extended projections are filtered and backprojected with a standard Feldkamp-type algorithm. Quantitative evaluations have been performed by simulations of voxelized phantoms on the basis of the root mean square deviation and a quality factor Q (Q = 1 represents the ideal correction). Measurements with a C-arm FDCT system have been used to validate the simulations and to investigate the practical applicability using anthropomorphic phantoms which caused truncation in all projections. The proposed approaches enlarged the FOV to cover wider patient cross-sections. Thus, image quality inside and outside the SFOV has been improved. Best results have been obtained using the CT-based method, followed by the camera-based and the FDCT-based truncation correction. For simulations, quality factors up to 0.98 have been achieved. Truncation-induced cupping artefacts have been reduced, e.g., from 218% to less than 1% for the measurements. The proposed truncation correction approaches for EFOV reconstructions are an effective way to ensure accurate CT values inside the SFOV and to recover peripheral information outside the SFOV.

Comparison of analytic and algebraic methods for motion-compensated cone-beam CT reconstruction of the thorax

Respiratory motion is a major concern in cone-beam (CB) computed tomography (CT) of the thorax. It causes artifacts such as blur, streaks, and bands, in particular when using slow-rotating scanners mounted on the gantry of linear accelerators. In this paper, we compare two approaches for motion-compensated CBCT reconstruction of the thorax. The first one is analytic; it is heuristically adapted from the method of Feldkamp, Davis, and Kress (FDK). The second one is algebraic: the system of linear equations is generated using a new algorithm for the projection of deformable volumes and solved using the Simultaneous Algebraic Reconstruction Technique (SART). For both methods, we propose to estimate the motion on patient data using a previously acquired 4-D CT image. The methods were tested on two digital and one mechanical motion-controlled phantoms and on a patient dataset. Our results indicate that the two methods correct most motion artifacts. However, the analytic method does not fully correct streaks and bands even if the motion is perfectly estimated due to the underlying approximation. In contrast, the algebraic method allows us full correction of respiratory-induced artifacts.

Exact iterative reconstruction for the interior problem

There is a trend in single photon emission computed tomography (SPECT) that small and dedicated imaging systems are becoming popular. For example, many companies are developing small dedicated cardiac SPECT systems with different designs. These dedicated systems have a smaller field of view (FOV) than a full-size clinical system. Thus data truncation has become the norm rather than the exception in these systems. Therefore, it is important to develop region of interest (ROI) reconstruction algorithms using truncated data. This paper is a stepping stone toward this direction. This paper shows that the common generic iterative image reconstruction algorithms are able to exactly reconstruct the ROI under the conditions that the convex ROI is fully sampled and the image value in a sub-region within the ROI is known. If the ROI includes a sub-region that is outside the patient body, then the conditions can be easily satisfied.

Tiny a priori knowledge solves the interior problem in computed tomography

Based on the concept of differentiated backprojection (DBP) (Noo et al 2004 Phys. Med. Biol. 49 3903, Pan et al 2005 Med. Phys. 32 673, Defrise et al 2006 Inverse Problems 22 1037), this paper shows that the solution to the interior problem in computed tomography is unique if a tiny a priori knowledge on the object f(x, y) is available in the form that f(x, y) is known on a small region located inside the region of interest. Furthermore, we advance the uniqueness result to obtain more general uniqueness results which can be applied to a wider class of imaging configurations. We also develop a reconstruction algorithm which can be considered an extension of the DBP-POCS (projection onto convex sets) method described by Defrise et al (2006 Inverse Problems 22 1037), where we not only extend this method to the interior problem but also introduce a new POCS algorithm to reduce computational cost. Finally, we present experimental results which show evidence that the inversion corresponding to each obtained uniqueness result is stable.