A backprojection kernel (KRNL3D) for very-wide-aperture 3D tomography applied to PET with Multigrid for precise use...

In 'KRNL3D' we derive a kernel functionK(y₁,y₂,φ) whose backprojections from all directions (θ,φ) in the spherical band∣φ∣<φ¯maxon the celestial sphere, when integrated with respect to solid angle, yieldρ, the 3D Gaussian point response function (PRF) of radius 1. ThisK, when convolved against line integral data from an unknown density functionf, yields an integral formula for the 'mollification'ff=ρ∗f, which is a slightly blurred version off, and which stabilizes the mild ill-posedness. Applied to positron emission tomography that backprojection reconstruction occurs stochastically and one emission event at a time, after needed data corrections. We describe Octave (≈Matlab) codes to tabulateKand to test its use with a large apertureφ¯max=π/3orπ/6. 'KRNL3D-TOF' truncates backprojection to a cylindrical patch about the TOF approximate location of each event. These 'backplacements' decrease the computational cost and limit noise and streaking in one region from contaminating the reconstruction in more distant regions. They also retain the ability to count emission events in an isolated blob despiteverylow event counts, a valuable feature fordynamicstudies of metabolic processes. 'Multigrid' allows further reduction in the radius and lengths of the cylinders, thereby enabling even moreprecise use of the TOF information. This precision should be especially important as researchers decrease the TOF uncertainty in newer generation scanners. Finally, we discuss 'further work' that needs to be done. Our codes are being made freely available athttps://github.com/keithmillerberkeley/PET-codes.