Fourier-based reconstruction for fully 3-D PET: optimization of interpolation parameters

Fourier-based approaches for three-dimensional (3-D) reconstruction are based on the relationship between the 3-D Fourier transform (FT) of the volume and the two-dimensional (2-D) FT of a parallel-ray projection of the volume. The critical step in the Fourier-based methods is the estimation of the samples of the 3-D transform of the image from the samples of the 2-D transforms of the projections on the planes through the origin of Fourier space, and vice versa for forward-projection (reprojection). The Fourier-based approaches have the potential for very fast reconstruction, but their straightforward implementation might lead to unsatisfactory results if careful attention is not paid to interpolation and weighting functions. In our previous work, we have investigated optimal interpolation parameters for the Fourier-based forward and back-projectors for iterative image reconstruction. The optimized interpolation kernels were shown to provide excellent quality comparable to the ideal sinc interpolator. This work presents an optimization of interpolation parameters of the 3-D direct Fourier method with Fourier reprojection (3D-FRP) for fully 3-D positron emission tomography (PET) data with incomplete oblique projections. The reprojection step is needed for the estimation (from an initial image) of the missing portions of the oblique data. In the 3D-FRP implementation, we use the gridding interpolation strategy, combined with proper weighting approaches in the transform and image domains. We have found that while the 3-D reprojection step requires similar optimal interpolation parameters as found in our previous studies on Fourier-based iterative approaches, the optimal interpolation parameters for the main 3D-FRP reconstruction stage are quite different. Our experimental results confirm that for the optimal interpolation parameters a very good image accuracy can be achieved even without any extra spectral oversampling, which is a common practice to decrease errors caused by interpolation in Fourier reconstruction.     

An Iterative Image Space Reconstruction Algorthm Suitable for Volume ECT

The trend in the design of scanners for positron emission computed tomography has traditionally been to improve the transverse spatial resolution to several millimeters while maintaining relatively coarse axial resolution (1-2 cm). Several scanners are being built with fine sampling in the axial as well as transverse directions, leading to the possibility of the true volume imaging. The number of possible coincidence pairs in these scanners is quite large. The usual methods of image reconstruction cannot handle these data without making approximations. It is computationally most efficient to reduce the size of this large, sparsely populated array by back-projecting the coincidence data prior to reconstruction. While analytic reconstruction techniques exist for back-projected data, an iterative algorithm may be necessary for those cases where the point spread function is spatially variant. A modification of the maximum likelihood algorithm is proposed to reconstruct these back-projected data. The method, the iterative image space reconstruction algorithm (ISRA), is able to reconstruct data from a scanner with a spatially variant point spread function in less time than other proposed algorithms. Results are presented for single-slice data, simulated and actual, from the PENN-PET scanner.     

Rotate-and-slant projector for fast LOR-based fully-3-D iterative PET reconstruction

One of the greatest challenges facing iterative fully-3-D positron emission tomography (PET) reconstruction is the issue of long reconstruction times due to the large number of measurements for 3-D mode as compared to 2-D mode. A rotate-and-slant projector has been developed that takes advantage of symmetries in the geometry to compute volumetric projections to multiple oblique sinograms in a computationally efficient manner. It is based upon the 2-D rotation-based projector using the three-pass method of shears, and it conserves the 2-D rotator computations for multiple projections to each oblique sinogram set. The projector is equally applicable to both conventional evenly-spaced projections and unevenly-spaced line-of-response (LOR) data. The LOR-based version models the location and orientation of the individual LORs (i.e., the arc-correction), providing an ordinary Poisson reconstruction framework. The projector was implemented in C with several optimizations for speed, exploiting the vertical symmetry of the oblique projection process, depth compression, and array indexing schemes which maximize serial memory access. The new projector was evaluated and compared to ray-driven and distance-driven projectors using both analytical and experimental phantoms, and fully-3-D iterative reconstructions with each projector were also compared to Fourier rebinning with 2-D iterative reconstruction. In terms of spatial resolution, contrast, and background noise measures, 3-D LOR-based iterative reconstruction with the rotate-and-slant projector performed as well as or better than the other methods. Total processing times, measured on a single cpu Linux workstation, were approximately 10x faster for the rotate-and-slant projector than for the other 3-D projectors studied. The new projector provided four iterations fully-3-D ordered-subsets reconstruction in as little as 15 s--approximately the same time as FORE + 2-D reconstruction. We conclude that the rotate-and-slant projector is a viable option for fully-3-D PET, offering quality statistical reconstruction in times only marginally slower than 2-D or rebinning methods.     

Generalized true three-dimensional reconstruction algorithm

A generalized true three-dimensional reconstruction (TTR) algorithm applicable to a sphere or cylindrical geometry is described. A new composite filter concept which is applicable to all forms of three-dimensional geometries, especially suitable for a truncated sphere is introduced and incorporated into image reconstruction of physically realizable detector array system such as the truncated spherical positron detector array.     

Exact and approximate rebinning algorithms for 3-D PET data

This paper presents two new rebinning algorithms for the reconstruction of three-dimensional (3-D) positron emission tomography (PET) data. A rebinning algorithm is one that first sorts the 3-D data into an ordinary two-dimensional (2-D) data set containing one sinogram for each transaxial slice to be reconstructed; the 3-D image is then recovered by applying to each slice a 2-D reconstruction method such as filtered-backprojection. This approach allows a significant speedup of 3-D reconstruction, which is particularly useful for applications involving dynamic acquisitions or whole-body imaging. The first new algorithm is obtained by discretizing an exact analytical inversion formula. The second algorithm, called the Fourier rebinning algorithm (FORE), is approximate but allows an efficient implementation based on taking 2-D Fourier transforms of the data. This second algorithm was implemented and applied to data acquired with the new generation of PET systems and also to simulated data for a scanner with an 18° axial aperture. The reconstructed images were compared to those obtained with the 3-D reprojection algorithm (3DRP) which is the standard “exact” 3-D filtered-backprojection method. Results demonstrate that FORE provides a reliable alternative to 3DRP, while at the same time achieving an order of magnitude reduction in processing time     

Direct reconstruction of kinetic parameter images from dynamic PET data

Our goal in this paper is the estimation of kinetic model parameters for each voxel corresponding to a dense three-dimensional (3-D) positron emission tomography (PET) image. Typically, the activity images are first reconstructed from PET sinogram frames at each measurement time, and then the kinetic parameters are estimated by fitting a model to the reconstructed time-activity response of each voxel. However, this "indirect" approach to kinetic parameter estimation tends to reduce signal-to-noise ratio (SNR) because of the requirement that the sinogram data be divided into individual time frames. In 1985, Carson and Lange proposed, but did not implement, a method based on the expectation-maximization (EM) algorithm for direct parametric reconstruction. The approach is "direct" because it estimates the optimal kinetic parameters directly from the sinogram data, without an intermediate reconstruction step. However, direct voxel-wise parametric reconstruction remained a challenge due to the unsolved complexities of inversion and spatial regularization. In this paper, we demonstrate and evaluate a new and efficient method for direct voxel-wise reconstruction of kinetic parameter images using all frames of the PET data. The direct parametric image reconstruction is formulated in a Bayesian framework, and uses the parametric iterative coordinate descent (PICD) algorithm to solve the resulting optimization problem. The PICD algorithm is computationally efficient and is implemented with spatial regularization in the domain of the physiologically relevant parameters. Our experimental simulations of a rat head imaged in a working small animal scanner indicate that direct parametric reconstruction can substantially reduce root-mean-squared error (RMSE) in the estimation of kinetic parameters, as compared to indirect methods, without appreciably increasing computation.     

Weighted backprojection approach to cone beam 3D projection reconstruction for truncated spherical detection geometry

A new analytical three-dimensional cone beam reconstruction algorithm is presented for truncated spherical detection geometry. The basic idea of the proposed algorithm is the formation of spatially invariant 3D blurred back-projected volumetric image by the use of the weighted backprojection of cone beam projection data and subsequent 3D filtering using an acceptance angle dependent rho filter. The backprojection weighting function is calculated on the basis of each given geometrical condition, i.e. detection geometry or degree of truncation, position of cone beam apex, and backprojection point. The proposed algorithm is derived analytically and is computationally efficient. Performance of the algorithm is evaluated by the reconstruction of 3D volumetric images using simulated data from arbitrarily truncated spherical detector geometries.     

A multigrid expectation maximization reconstruction algorithm for positron emission tomography

The problem of reconstruction in positron emission tomography (PET) is basically estimating the number of photon pairs emitted from the source. Using the concept of the maximum-likelihood (ML) algorithm, the problem of reconstruction is reduced to determining an estimate of the emitter density that maximizes the probability of observing the actual detector count data over all possible emitter density distributions. A solution using this type of expectation maximization (EM) algorithm with a fixed grid size is severely handicapped by the slow convergence rate, the large computation time, and the nonuniform correction efficiency of each iteration, which makes the algorithm very sensitive to the image pattern. An efficient knowledge-based multigrid reconstruction algorithm based on the ML approach is presented to overcome these problems.     

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 three-dimensional (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 two-dimensional (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 line of responses (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.     

A 3D reconstruction of vascular structures from two X-ray angiograms using an adapted simulated annealing algorithm

A three-dimensional (3D) reconstruction of the vessel lumen from two angiographic views, based on the reconstruction of a series of cross-sections, is proposed. Assuming uniform mixing of contrast medium and background subtraction, the cross-section of each vessel is reconstructed through a binary representation. A priori information about both the slice to be reconstructed and the relationships between adjacent slices are incorporated to lessen ambiguities on the reconstruction. Taking into account the knowledge of normal vessel geometry, an initial solution of each slice is created using an elliptic model-based method. This initial solution is then deformed to be made consistent with projection data while being constrained into a connected realistic shape. For that purpose, properties on the expected optimal solution are described through a Markov random field. To find an optimal solution, a specific optimization algorithm based on simulated annealing is used. The method performs well both on single vessels and on branching vessels possessing an additional inherent ambiguity when viewed at oblique angles. Results on 2D slice independent reconstruction and 3D reconstruction of a stack of spatially continuous 2D slices are presented for single vessels and bifurcations.     

Efficient interpolators and filter banks using multiplier blocks

We examine the use of efficient shift-and-add multiplier structures and multiplier blocks to reduce computational complexity in filter banks. This is more efficient than treating each bank filter separately. We also examine the Farrow (see Proc. Int. Symp. Circuits Syst. (ISCAS), p. 2641-2645) structure, which is used in interpolators. Applying multiplier blocks makes this structure cheaper than the more recognized Lagrange interpolator     

A methodology for specifying PET VOIs using multimodalitytechniques

Volume-of-interest (VOI) extraction for radionuclide and anatomical measurements requires correct identification and delineation of the anatomical feature being studied. The authors have developed a toolset for specifying three dimensional (3-D) VOIs on a multislice positron emission tomography (PET) dataset. The software is particularly suited for specifying cerebral cortex VOIs which represent a particular gyrus or deep brain structure. A registered 3-D magnetic resonance image (MRI) dataset is used to provide high-resolution anatomical information, both as oblique two-dimensional (2-D) sections and as volume renderings of a segmented cortical surface. VOIs are specified indirectly in two dimensions by drawing a stack of 2-D regions on the MRI data. The regions are tiled together to form closed triangular mesh surface models, which are subsequently transformed into the observation space of the PET scanner. Quantification by this method allows calculation of radionuclide activity in the VOIs, as well as their statistical uncertainties and correlations. The methodology for this type of analysis and validation results are presented     

基于稀疏重构的全极化SAR联合多维重建

各极化通道独立处理和三维分步成像会忽视数据之间的关联性,造成散射中心的失配以及极化散射矩阵获取的不准确。鉴于此,该文提出一种基于稀疏重构的全极化联合多维重建方法。该方法通过设置联合稀疏约束对所有极化通道及所有维度进行联合,将全极化多维重建建模为多通道联合稀疏重构问题。通过数据插值对模型简化后,结合三维快速傅里叶变换、共轭梯度法和牛顿迭代法给出一种高效的模型求解方法,可以同时得到极化散射矩阵和目标三维信息。该文方法保证了不同极化通道、不同维度的稀疏支撑集一致,且充分利用了数据之间的关联性带来的额外信息。基于仿真数据和电磁计算数据的实验结果表明,该方法的性能不受目标类型影响,具有一定的抗噪性,能有效地获取目标的多维重建结果,得到的三维成像结果分辨率高且极化散射矩阵估计精度高。      

Inverse scattering method for one-dimensional inhomogeneous layered media

An inverse scattering method to reconstruct simultaneously the permittivity profile and the conductivity profile of one-dimensional inhomogeneons medium which makes use of the transverse electric (TE) wave and/or transverse magnetic (TM) wave, is proposed. The medium is illuminated by the TE and/or TM plane wave at oblique incidence, and the data are taken as the reflection coefficients for a set of discrete frequencies and/or a finite number of incident angles. Furthermore, the reflection coefficient data contain the Gaussian noise. The nonlinear integral equation relating the unknown constitutive parameter of the medium to the reflection coefficient for TE wave and/or TM wave is solved by the Newton iteration method. The inverse operator in the Newton method is determined by the regularization method. It is demonstrated in terms of the numerical examples that this method utilizing both polarizations and the incident angle of the incident plane wave is very effective even if the reflection coefficient contains the practical measurement error, or the phase of the reflection coefficient is unknown. Moreover, the relationships between the errors of reconstructed profile and the measured reflection coefficient are also discussed.     

Novel scatter compensation of list-mode PET data using spatial and energy dependent corrections

With the widespread use of positron emission tomography (PET) crystals with greatly improved energy resolution (e.g., 11.5% with LYSO as compared to 20% with BGO) and of list-mode acquisitions, the use of the energy of individual events in scatter correction schemes becomes feasible. We propose a novel scatter approach that incorporates the energy of individual photons in the scatter correction and reconstruction of list-mode PET data in addition to the spatial information presently used in clinical scanners. First, we rewrite the Poisson likelihood function of list-mode PET data including the energy distributions of primary and scatter coincidences and show that this expression yields an MLEM reconstruction algorithm containing both energy and spatial dependent corrections. To estimate the spatial distribution of scatter coincidences we use the single scatter simulation (SSS). Next, we derive two new formulae which allow estimation of the 2-D (coincidences) energy probability density functions (E-PDF) of primary and scatter coincidences from the 1-D (photons) E-PDFs associated with each photon. We also describe an accurate and robust object-specific method for estimating these 1-D E-PDFs based on a decomposition of the total energy spectra detected across the scanner into primary and scattered components. Finally, we show that the energy information can be used to accurately normalize the scatter sinogram to the data. We compared the performance of this novel scatter correction incorporating both the position and energy of detected coincidences to that of the traditional approach modeling only the spatial distribution of scatter coincidences in 3-D Monte Carlo simulations of a medium cylindrical phantom and a large, nonuniform NCAT phantom. Incorporating the energy information in the scatter correction decreased bias in the activity distribution estimation by ~20% and ~40% in the cold regions of the large NCAT phantom at energy resolutions 11.5% and 20% at 511 keV, respectively, compared to when using the spatial information alone.     

Performance evaluation of the PC-2048: a new 15-slice encoded-crystal PET scanner for neurological studies

Initial experience is reported with the Scanditronix PC 2048-15B, a 15-slice positron emission tomography (PET) system using multicrystal/multiphoto-multiplier modules to obtain high spatial resolution. Random and scattered events are reduced using an orbiting (68)Ge rod source for transmission scans by only accepting coincidence lines which intersect the instantaneous position of the source. Scatter correction of the emission data is removed with a deconvolution kernel, random and dead-time correction by the use of observed singles rates. The peak count rates are 11.7/20.0 Kcps for the direct cross slices at concentrations of 4.5/5.1 muCi/cc. respectively. Over the radial range 0-9 cm from the ring center, radial transverse resolution is 4.6-6.4 mm, aid tangential transverse resolution is 4.6-5.1 mm using a Hanning filter. Over the same range, axial resolution varies from 6.1-6.2 mm in direct slices and from 5.4-7.1 mm in cross slices. This near-isotropic resolution allows collection of image volume data with no preferred direction for signal averaging errors.     

Tomographic reconstruction using an adaptive tetrahedral mesh defined by a point cloud

Medical images in nuclear medicine are commonly represented in three dimensions as a stack of two-dimensional images that are reconstructed from tomographic projections. Although natural and straightforward, this may not be an optimal visual representation for performing various diagnostic tasks. A method for three-dimensional (3-D) tomographic reconstruction is developed using a point cloud image representation. A point cloud is a set of points (nodes) in space, where each node of the point cloud is characterized by its position and intensity. The density of the nodes determines the local resolution allowing for the modeling of different parts of the image with different resolution. The reconstructed volume, which in general could be of any resolution, size, shape, and topology, is represented by a set of nonoverlapping tetrahedra defined by the nodes. The intensity at any point within the volume is defined by linearly interpolating inside a tetrahedron from the values at the four nodes that define the tetrahedron. This approach creates a continuous piecewise linear intensity over the reconstruction domain. The reconstruction provides a distinct multiresolution representation, which is designed to accurately and efficiently represent the 3-D image. The method is applicable to the acquisition of any tomographic geometry, such as parallel-, fan-, and cone-beam; and the reconstruction procedure can also model the physics of the image detection process. An efficient method for evaluating the system projection matrix is presented. The system matrix is used in an iterative algorithm to reconstruct both the intensity and location of the distribution of points in the point cloud. Examples of the reconstruction of projection data generated by computer simulations and projection data experimentally acquired using a Jaszczak cardiac torso phantom are presented. This work creates a framework for voxel-less multiresolution representation of images in nuclear medicine.     

Mean and covariance properties of dynamic PET reconstructions from list-mode data

We derive computationally efficient methods for the estimation of the mean and variance properties of penalized likelihood dynamic positron emission tomography (PET) images. This allows us to predict the accuracy of reconstructed activity estimates and to compare reconstruction algorithms theoretically. We combine a bin-mode approach in which data is modeled as a collection of independent Poisson random variables at each spatiotemporal bin with the space-time separabilities in the imaging equation and penalties to derive rapidly computable analytic mean and variance approximations. We use these approximations to compare bias/variance properties of our dynamic PET image reconstruction algorithm with those of multiframe static PET reconstructions.     

Elimination of nonpivotal plane images from x-ray motion tomograms

A method for eliminating the contribution of nonpivotal plane images into a conventional X-ray motion tomogram is described. A set of tomograms, each focused on one of distinct parallel planes, are combined to form a blur-free tomogram. From the mathematical analysis of conventional tomography, it is found that the multisection tomography system reduces to a three-dimensional incomplete linear system with angularly limited frequency response. An iterative constrained reconstruction algorithm incorporating the object boundary constraint is used for reconstruction of a blur-free tomogram without a contribution from nonpivotal planes. Some computer simulations were made for verification of the algorithm. Although the algorithm accomplishes the three-dimensional reconstruction of volume rather than image restoration of a single slice plane, the optimal high-pass spatial filter applying to a single tomogram is derived on the way to solving the inverse problem of the multisection tomography system.