Exact Rebinning Research Articles

Exact and Approximate Fourier Rebinning Algorithms for the Solution of the Data Truncation Problem in 3-D PET

This paper presents an extended 3-D exact rebinning formula in the Fourier space that leads to an iterative reprojection algorithm (iterative FOREPROJ), which enables the estimation of unmeasured oblique projection data on the basis of the whole set of measured data. In first approximation, this analytical formula also leads to an extended Fourier rebinning equation that is the basis for an approximate reprojection algorithm (extended FORE). These algorithms were evaluated on numerically simulated 3-D positron emission tomography (PET) data for the solution of the truncation problem, i.e., the estimation of the missing portions in the oblique projection data, before the application of algorithms that require complete projection data such as some rebinning methods (FOREX) or 3-D reconstruction algorithms (3DRP or direct Fourier methods). By taking advantage of all the 3-D data statistics, the iterative FOREPROJ reprojection provides a reliable alternative to the classical FOREPROJ method, which only exploits the low-statistics nonoblique data. It significantly improves the quality of the external reconstructed slices without loss of spatial resolution. As for the approximate extended FORE algorithm, it clearly exhibits limitations due to axial interpolations, but will require clinical studies with more realistic measured data in order to decide on its pertinence.

Iterative Image Reconstruction Using Inverse Fourier Rebinning for Fully 3-D PET

We describe a fast forward and back projector pair based on inverse Fourier rebinning for use in iterative image reconstruction for fully 3-D positron emission tomography (PET). The projector pair is used as part of a factored system matrix that takes into account detector-pair response by using shift-variant sinogram blur kernels, thereby combining the computational advantages of Fourier rebinning with iterative reconstruction using accurate system models. The forward projector consists of a 2-D projector, which maps 3-D images into 2-D direct sinograms, followed by exact inverse rebinning which maps the 2-D into fully 3-D sinograms. The back projector is implemented as the transpose of the forward projector and differs from the true exact rebinning operator in the sense that it does not require reprojection to compute missing lines of response (LORs). We compensate for two types of inaccuracies that arise in a cylindrical PET scanner when using inverse Fourier rebinning: 1) nonuniform radial sampling and 2) nonconstant oblique angles in the radial direction in a single oblique sinogram. We examine the effects of these corrections on sinogram accuracy and reconstructed image quality. We evaluate performance of the new projector pair for maximum a posteriori (MAP) reconstruction of simulated and in vivo data. The new projector results in only a small loss in resolution towards the edge of the field-of-view when compared to the fully 3-D geometric projector and requires an order of magnitude less computation.

Fourier rebinning of time-of-flight PET data

This paper explores fast reconstruction strategies for 3D time-of-flight (TOF) positron emission tomography (PET), based on 2D data rebinning. Starting from pre-corrected 3D TOF data, a rebinning algorithm estimates for each transaxial slice the 2D TOF sinogram that would have been acquired by a single-ring scanner. The rebinned sinograms can then be reconstructed using any algorithm for 2D TOF reconstruction. We introduce TOF-FORE, an approximate rebinning algorithm obtained by extending the Fourier rebinning method for non-TOF data. In addition, we identify two partial differential equations that must be satisfied by consistent 3D TOF data, and use them to derive exact rebinning algorithms and to characterize the degree of the approximation in TOF-FORE. Numerical simulations demonstrate that TOF-FORE is more accurate than two different TOF extensions of the single-slice rebinning method, and suggest that TOF-FORE will be a valuable tool for practical TOF PET in the range of axial apertures and time resolutions typical of current scanners.

Improved two-dimensional rebinning of helical cone-beam computerized tomography datausing John’s equation

This paper deals with three-dimensional (3D) image reconstruction for helical computerizedtomography with multi-row detectors. We describe a method to improve the accuracy ofthe rebinning algorithms, which separate 3D reconstruction into independent 2Dreconstructions for a set of oblique slices. Each oblique slice is reconstructed from anestimate of its 2D Radon transform obtained from the cone-beam projections acquiredwhile the x-ray source moves along a segment of the helix. Because this helixsegment is not contained within the oblique slice, the estimated 2D Radon transformis approximate. In theory, exact rebinning could be achieved by solving John’spartial differential equation to virtually move the x-ray source within the obliqueslice. In contrast with previous work by Patch (2002 IEEE Trans. Med. Imaging 21801–13), we do not attempt to solve John’s equation exactly. Instead, we use John’sequation to compute a first order correction to the rebinning algorithm. Testswith simulated data demonstrate a significant improvement of image quality,obtained with a negligible increase of the computation time and of the sensitivity tonoise.

A fast rebinning algorithm for 3D positron emission tomography using John's equation

Volume imaging in positron emission tomography (PET) requires the inversion of the three-dimensional (3D) x-ray transform. The usual solution to this problem is based on 3D filtered-backprojection (FBP), but is slow. Alternative methods have been proposed which factor the 3D data into independent 2D data sets corresponding to the 2D Radon transforms of a stack of parallel slices. Each slice is then reconstructed using 2D FBP. These so-called rebinning methods are numerically efficient but are approximate. In this paper a new exact rebinning method is derived by exploiting the fact that the 3D x-ray transform of a function is the solution to the second-order partial differential equation first studied by John. The method is proposed for two sampling schemes, one corresponding to a pair of infinite plane detectors and another one corresponding to a cylindrical multi-ring PET scanner. The new FORE-J algorithm has been implemented for this latter geometry and was compared with the approximate Fourier rebinning algorithm FORE and with another exact rebinning algorithm, FOREX. Results with simulated data demonstrate a significant improvement in accuracy compared to FORE, while the reconstruction time is doubled. Compared to FOREX, the FORE-J algorithm is slightly less accurate but more than three times faster.

Exact rebinning methods for three-dimensional PET.

The high computational cost of data processing in volume PET imaging is still hindering the routine application of this successful technique, especially in the case of dynamic studies. This paper describes two new algorithms based on an exact rebinning equation, which can be applied to accelerate the processing of three-dimensional (3-D) PET data. The first algorithm, FOREPROJ, is a fast-forward projection algorithm that allows calculation of the 3-D attenuation correction factors (ACF's) directly from a two-dimensional (2-D) transmission scan, without first reconstructing the attenuation map and then performing a 3-D forward projection. The use of FOREPROJ speeds up the estimation of the 3-D ACF's by more than a factor five. The second algorithm, FOREX, is a rebinning algorithm that is also more than five times faster, compared to the standard reprojection algorithm (3DRP) and does not suffer from the image distortions generated by the even faster approximate Fourier rebinning (FORE) method at large axial apertures. However, FOREX is probably not required by most existing scanners, as the axial apertures are not large enough to show improvements over FORE with clinical data. Both algorithms have been implemented and applied to data simulated for a scanner with a large axial aperture (30 degrees), and also to data acquired with the ECAT HR and the ECAT HR+ scanners. Results demonstrate the excellent accuracy achieved by these algorithms and the important speedup when the sinogram sizes are powers of two.