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.     

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.     

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 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.     

Simultaneous MR-Compatible Emission and Transmission Imaging for PET Using Time-of-Flight Information

Quantitative positron emission tomography (PET) imaging relies on accurate attenuation correction. Predicting attenuation values from magnetic resonance (MR) images is difficult because MR signals are related to proton density and relaxation properties of tissues. Here, we propose a method to derive the attenuation map from a transmission scan. An annulus transmission source is positioned inside the field-of-view of the PET scanner. First a blank scan is acquired. The patient is injected with FDG and placed inside the scanner. 511-keV photons coming from the patient and the transmission source are acquired simultaneously. Time-of-flight information is used to extract the coincident photons originating from the annulus. The blank and transmission data are compared in an iterative reconstruction method to derive the attenuation map. Simulations with a digital phantom were performed to validate the method. The reconstructed attenuation coefficients differ less than 5% in volumes of interest inside the lungs, bone, and soft tissue. When applying attenuation correction in the reconstruction of the emission data a standardized uptake value error smaller than 9% was obtained for all tissues. In conclusion, our method can reconstruct the attenuation map and the emission data from a simultaneous scan without prior knowledge about the anatomy or the attenuation coefficients of the tissues.     

A general class of preconditioners for statistical iterativereconstruction of emission computed tomography

A major drawback of statistical iterative image reconstruction for emission computed tomography is its high computational cost. The ill-posed nature of tomography leads to slow convergence for standard gradient-based iterative approaches such as the steepest descent or the conjugate gradient algorithm. Here, new theory and methods for a class of preconditioners are developed for accelerating the convergence rate of iterative reconstruction. To demonstrate the potential of this class of preconditioners, a preconditioned conjugate gradient (PCG) iterative algorithm for weighted least squares reconstruction (WLS) was formulated for emission tomography. Using simulated positron emission tomography (PET) data of the Hoffman brain phantom, it was shown that the convergence rate of the PCG can reduce the number of iterations of the standard conjugate gradient algorithm by a factor of 2-8 times depending on the convergence criterion     

Simultaneous maximum a posteriori reconstruction of attenuation andactivity distributions from emission sinograms

In order to perform attenuation correction in emission tomography an attenuation map is required. We propose a new method to compute this map directly from the emission sinogram, eliminating the transmission scan from the acquisition protocol. The problem is formulated as an optimization task where the objective function is a combination of the likelihood and an a priori probability. The latter uses a Gibbs prior distribution to encourage local smoothness and a multimodal distribution for the attenuation coefficients. Since the attenuation process is different in positron emission tomography (PET) and single photon emission tomography (SPECT), a separate algorithm for each case is derived. The method has been tested on mathematical phantoms and on a few clinical studies. For PET, good agreement was found between the images obtained with transmission measurements and those produced by the new algorithm in an abdominal study. For SPECT, promising simulation results have been obtained for nonhomogeneous attenuation due to the presence of the lungs.     

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.     

Efficient Monte Carlo based scatter artifact reduction in cone-beam micro-CT

Cupping and streak artifacts caused by the detection of scattered photons may severely degrade the quantitative accuracy of cone-beam X-ray computed tomography (CT) images. In order to overcome this problem, we propose and validate the following iterative scatter artifact reduction scheme. First, an initial image is reconstructed from the scatter-contaminated projections. Next, the scatter component of the projections is estimated from the initial reconstruction by a Monte Carlo (MC) simulation. The estimate obtained is then utilized during the reconstruction of a scatter-corrected image. The last two steps are repeated until an adequate correction is obtained. The estimation of the noise-free scatter projections in this scheme is accelerated in the following way: first, a rapid (i.e., based on a low number of simulated photon tracks) MC simulation is executed. The noisy result of this simulation is de-noised by a three-dimensional fitting of Gaussian basis functions. We demonstrate that, compared to plain MC, this method shortens the required simulation time by three to four orders of magnitude. Using simulated projections of a small animal phantom, we show that one cycle of the scatter correction scheme is sufficient to produce reconstructed images that barely differ from the reconstructions of scatter-free projections. The reconstructions of data acquired with a charge-coupled device based micro-CT scanner demonstrate a nearly complete removal of the scatter-induced cupping artifact. Quantitative errors in a water phantom are reduced from around 12% for reconstructions without the scatter correction to 1% after the proposed scatter correction has been applied. In conclusion, a general, accurate, and efficient scatter correction algorithm is developed that requires no mechanical modifications of the scanning equipment and results in only a moderate increase in the total reconstruction time.     

Fully three-dimensional reconstruction for a PET camera with retractable septa

A fully 3-D reconstruction algorithm has been developed to reconstruct data from a 16 ring PET camera (a Siemens/CTI 953B) with automatically retractable septa. The tomograph is able to acquire coincidences between any pair of detector rings and septa retraction increases the total system count rate by a factor of 7.8 (including scatter) and 4.7 (scatter subtracted) for a uniform, 20 cm diameter cylinder. The reconstruction algorithm is based on 3-D filtered backprojection, expressed in a form suitable for the multi-angle sinogram data. Sinograms which are not measured due to the truncated cylindrical geometry of the tomograph, but which are required for a spatially invariant response function, are obtained by forward projection. After filtering, the complete set of sinograms is backprojected into a 3-D volume of 128x128x31 voxels using a voxel-driven procedure. The algorithm has been validated with simulation, and tested with both phantom and clinical data from the 953B.     

2.5-D simultaneous multislice reconstruction by series expansion methods from Fourier-rebinned PET data

True three-dimensional (3-D) volume reconstruction from fully 3-D data in positron emission tomography (PET) has only a limited clinical use because of its large computational burden. Fourier rebinning (FORE) of the fully 3-D data into a set of 2-D sinogram data decomposes the 3-D reconstruction process into multiple 2-D reconstructions of decoupled 2-D image slices, thus substantially decreasing the computational burden even in the case when the 2-D reconstructions are performed by an iterative reconstruction algorithm. On the other hand, the approximations involved in the rebinning combined with the decoupling of the image slices cause a certain reduction of image quality, especially when the signal-to-noise ratio of the data is low. We propose a 2.5-D Simultaneous Multislice Reconstruction approach, based on the series expansion principle, where the volume is represented by the superposition of 3-D spherically symmetric bell-shaped basis functions. It takes advantage of the time reduction due to the use of the FORE (2-D) data, instead of the original fully 3-D data, but at the same time uses a 3-D iterative reconstruction approach with 3-D basis functions. The same general approach can be applied to any reconstruction algorithm belonging to the class of series expansion methods (iterative or noniterative) using 3-D basis functions that span multiple slices, and can be used for any multislice sinogram or list mode data whether obtained by a special rebinning scheme or acquired directly by a PET scanner in the 2-D mode using septa. Our studies confirm that the proposed 2.5-D approach provides a considerable improvement in reconstruction quality, as compared to the standard 2-D reconstruction approach, while the reconstruction time is of the same order as that of the 2-D approach and is clinically practical even on a general-purpose computer.     

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.     

A projector/backprojector with slice-to-slice blurring forefficient three-dimensional scatter modeling

Scatter correction is an important factor in single photon emission computed tomography (SPECT). Many scatter correction techniques, such as multiple-window subtraction and intrinsic modeling with iterative algorithms, have been under study for many years. Previously, we developed an efficient slice-to-slice blurring technique to model attenuation and system geometric response in a projector/backprojector pair, which was used in an ML-EM algorithm to reconstruct SPECT data. This paper proposes a projector/backprojector that models the three-dimensional (3-D) first-order scatter in SPECT, also using an efficient slice-to-slice blurring technique. The scatter response is estimated from a known nonuniform attenuation distribution map. It is assumed that the probability of detection of a first-order scattered photon from a photon that is emitted in a given source voxel and scattered in a given scatter voxel is proportional to the attenuation coefficient value at that voxel. Monte Carlo simulations of point sources and an MCAT torso phantom were used to verify the accuracy of the proposed projector/backprojector model. An experimental Jaszczak torso/cardiac phantom SPECT study was also performed. For a 64 x 64 x 64 image volume, it took 8.7 s to perform each iteration per slice on a Sun ULTRA Enterprise 3000 (167 MHz, 1 Gbyte RAM) computer, when modeling 3-D scatter, attenuation, and system geometric response functions. The main advantage of the proposed method is its easy implementation and the possibility of performing reconstruction in clinically acceptable time.     

Non-Gaussian smoothing of low-count transmission scans for PET whole-body studies

A non-Gaussian smoothing (NGS) technique is developed for filtering low count transmission (TR) data to be used for attenuation correction (AC) of positron emission tomography (PET) studies. The method is based on a statistical technique known as the generalized linear mixed model that allows an inverse link function that avoids the inversion of the observed transmission data. The NGS technique has been implemented in the sinogram domain in one-dimensional mode as angle-by-angle computation. To make it adaptive as a function of the TR count statistics we also develop and validate an objective procedure to choose an optimal smoothing parameter. The technique is assessed using experimental phantoms, simulating PET whole-body studies, and applied to real patient data. Different experimental conditions, in terms of TR scan time (from 1 h to 1 min), covering a wide range of TR counting statistic are considered. The method is evaluated, in terms of mean squared error (MSE), by comparing pixel by pixel the distribution for high counts statistics TR scan (1 h) with the corresponding counts distribution for low count statistics TR scans (e.g., 1 min). The smoothing parameter selection is shown to have high efficiency, meaning that it tends to choose values close to the unknown best value. Furthermore, the counts distribution of emission (EM) images, reconstructed with AC generated using low count TR data (1 min), are within 5% of the corresponding EM images reconstructed with AC generated using the high count statistics TR data (1 h). An application to a real patient whole-body PET study shows the promise of the technique for routine use.     

Accurate event-driven motion compensation in high-resolution PET incorporating scattered and random events

With continuing improvements in spatial resolution of positron emission tomography (PET) scanners, small patient movements during PET imaging become a significant source of resolution degradation. This work develops and investigates a comprehensive formalism for accurate motion-compensated reconstruction which at the same time is very feasible in the context of high-resolution PET. In particular, this paper proposes an effective method to incorporate presence of scattered and random coincidences in the context of motion (which is similarly applicable to various other motion correction schemes). The overall reconstruction framework takes into consideration missing projection data which are not detected due to motion, and additionally, incorporates information from all detected events, including those which fall outside the field-of-view following motion correction. The proposed approach has been extensively validated using phantom experiments as well as realistic simulations of a new mathematical brain phantom developed in this work, and the results for a dynamic patient study are also presented.      

Performance characteristics of the 3-D OSEM algorithm in the reconstruction of small animal PET images. Ordered-subsets expectation-maximixation

Rat brain images acquired with a small animal positron emission tomography (PET) camera and reconstructed with the three-dimensional (3-D) ordered-subsets expectation-maximization (OSEM) algorithm with resolution recovery have better quality when the brain is imaged by itself than when inside the head with surrounding background activity. The purpose of this study was to characterize the dependence of this effect on the level of background activity, attenuation, and scatter. Monte Carlo simulations of the imaging system were performed. The coefficient of variation from replicate images, full-width at half-maximum (FWHM) from point sources and image profile fitting, and image contrast and uniformity were used to evaluate algorithm performance. A rat head with the typical levels of five and ten times the brain activity in the surrounding background requires additional iterations to achieve the same resolution as the brain-only case at a cost of 24% and 64% additional noise, respectively. For the same phantoms, object scatter reduced contrast by 3%-5%. However, attenuation degraded resolution by 0.2 mm and was responsible for up to 12% nonuniformity in the brain images suggesting that attenuation correction is useful. Given the effects of emission and attenuation distribution on both resolution and noise, simulations or phantom studies should be used for each imaging situation to select the appropriate number of OSEM iterations to achieve the desired resolution-noise levels.     

Emission image reconstruction for randoms-precorrected PET allowing negative sinogram values

Most positron emission tomography (PET) emission scans are corrected for accidental coincidence (AC) events by real-time subtraction of delayed-window coincidences, leaving only the randoms-precorrected data available for image reconstruction. The real-time randoms precorrection compensates in mean for AC events but destroys the Poisson statistics. The exact log-likelihood for randoms-precorrected data is inconvenient, so practical approximations are needed for maximum likelihood or penalized-likelihood image reconstruction. Conventional approximations involve setting negative sinogram values to zero, which can induce positive systematic biases, particularly for scans with low counts per ray. We propose new likelihood approximations that allow negative sinogram values without requiring zero-thresholding. With negative sinogram values, the log-likelihood functions can be nonconcave, complicating maximization; nevertheless, we develop monotonic algorithms for the new models by modifying the separable paraboloidal surrogates and the maximum-likelihood expectation-maximization (ML-EM) methods. These algorithms ascend to local maximizers of the objective function. Analysis and simulation results show that the new shifted Poisson (SP) model is nearly free of systematic bias yet keeps low variance. Despite its simpler implementation, the new SP performs comparably to the saddle-point model which has shown the best performance (as to systematic bias and variance) in randoms-precorrected PET emission reconstruction.     

Noise characteristics of 3-D and 2-D PET images

The authors analyzed the noise characteristics of two-dimensional (2-D) and three-dimensional (3-D) images obtained from the GE Advance positron emission tomography (PET) scanner. Three phantoms were used: a uniform 20-cm phantom, a 3-D Hoffman brain phantom, and a chest phantom with heart and lung inserts. Using gated acquisition, the authors acquired 20 statistically equivalent scans of each phantom in 2-D and 3-D modes at several activity levels. From these data, they calculated pixel normalized standard deviations (NSD's), scaled to phantom mean, across the replicate scans, which allowed them to characterize the radial and axial distributions of pixel noise. The authors also performed sequential measurements of the phantoms in 2-D and 3-D modes to measure noise (from interpixel standard deviations) as a function of activity. To compensate for the difference in axial slice width between 2-D and 3-D images (due to the septa and reconstruction effects), they developed a smoothing kernel to apply to the 2-D data. After matching the resolution, the ratio of image-derived NSD values (NSD_2D/NSD_3D)² averaged throughout the uniform phantom was in good agreement with the noise equivalent count (NEC) ratio (NEC_3D/NEC_2D). By comparing different phantoms, the authors showed that the attenuation and emission distributions influence the spatial noise distribution. The estimates of pixel noise for 2-D and 3-D images produced here can be applied in the weighting of PET kinetic data and may be useful in the design of optimal dose and scanning requirements for PET studies. The accuracy of these phantom-based noise formulas should be validated for any given imaging situation, particularly in 3-D, if there is significant activity outside the scanner field of view     

Quantitative PET with positron emitters that emit prompt gamma rays

The purpose of this work was to determine the feasibility of using positron emitting isotopes that emit prompt gammas to acquire quantitative positron emission tomography (PET) data using standard PET instrumentation. Prompt gammas can contaminate PET data by increasing dead time, converting singles into invalid coincidences, and producing multiple coincidences which can lead to the replacement of valid coincidences by invalid coincidences. The measurements in this work were made by scanning point sources containing F-18, Na-22, and Co-60 and studying the effects of the prompt gammas on the PET data, We found that for the Na-22 point source, the annihilation photon coincidence rate was about 25 times the prompt gamma-annihilation photon coincidence rate in the entire active volume of the scanner. With scatter, the Na-22 prompt gamma-annihilation photon coincidence rate was 1.3 times higher than the F-18 scatter coincidence rate. The most significant effect of the prompt gamma was to increase dead time; the dead time correction factor for Cu-60 was 2.4 times higher than the correction factor for N-13 for the same source activity. We conclude that, in many cases, quantitative PET data can be readily obtained with isotopes that emit prompt gammas, using standard PET 2-D instrumentation. However there are some cases, such as 3-D PET, where prompt gammas could significantly contaminate the PET data.     

Optimal rebinning of time-of-flight PET data

Time-of-flight (TOF) positron emission tomography (PET) scanners offer the potential for significantly improved signal-to-noise ratio (SNR) and lesion detectability in clinical PET. However, fully 3D TOF PET image reconstruction is a challenging task due to the huge data size. One solution to this problem is to rebin TOF data into a lower dimensional format. We have recently developed Fourier rebinning methods for mapping TOF data into non-TOF formats that retain substantial SNR advantages relative to sinograms acquired without TOF information. However, mappings for rebinning into non-TOF formats are not unique and optimization of rebinning methods has not been widely investigated. In this paper we address the question of optimal rebinning in order to make full use of TOF information. We focus on FORET-3D, which approximately rebins 3D TOF data into 3D non-TOF sinogram formats without requiring a Fourier transform in the axial direction. We optimize the weighting for FORET-3D to minimize the variance, resulting in H(2)-weighted FORET-3D, which turns out to be the best linear unbiased estimator (BLUE) under reasonable approximations and furthermore the uniformly minimum variance unbiased (UMVU) estimator under Gaussian noise assumptions. This implies that any information loss due to optimal rebinning is as a result only of the approximations used in deriving the rebinning equation and developing the optimal weighting. We demonstrate using simulated and real phantom TOF data that the optimal rebinning method achieves variance reduction and contrast recovery improvement compared to nonoptimized rebinning weightings. In our preliminary study using a simplified simulation setup, the performance of the optimal rebinning method was comparable to that of fully 3D TOF MAP.