FIRST: Fast Iterative Reconstruction Software for (PET) tomography
Open Access
- 30 August 2006
- journal article
- Published by IOP Publishing in Physics in Medicine & Biology
- Vol. 51 (18), 4547-4565
- https://doi.org/10.1088/0031-9155/51/18/007
Abstract
Small animal PET scanners require high spatial resolution and good sensitivity. To reconstruct high-resolution images in 3D-PET, iterative methods, such as OSEM, are superior to analytical reconstruction algorithms, although their high computational cost is still a serious drawback. The higher performance of modern computers could make iterative image reconstruction fast enough to be viable, provided we are able to deal with the large number of probability coefficients for the system response matrix in high-resolution PET scanners, which is a difficult task that prevents the algorithms from reaching peak computing performance. Considering all possible axial and in-plane symmetries, as well as certain quasi-symmetries, we have been able to reduce the memory requirements to store the system response matrix (SRM) well below 1 GB, which allows us to keep the whole response matrix of the system inside RAM of ordinary industry-standard computers, so that the reconstruction algorithm can achieve near peak performance. The elements of the SRM are stored as cubic spline profiles and matched to voxel size during reconstruction. In this way, the advantages of 'on-the-fly' calculation and of fully stored SRM are combined. The on-the-fly part of the calculation (matching the profile functions to voxel size) of the SRM accounts for 10-30% of the reconstruction time, depending on the number of voxels chosen. We tested our approach with real data from a commercial small animal PET scanner. The results (image quality and reconstruction time) show that the proposed technique is a feasible solution.Keywords
This publication has 24 references indexed in Scilit:
- Overview of methods for image reconstruction from projections in emission computed tomographyProceedings of the IEEE, 2003
- A penalized-likelihood image reconstruction method for emission tomography, compared to postsmoothed maximum-likelihood with matched spatial resolutionIEEE Transactions on Medical Imaging, 2003
- Edge-preserving tomographic reconstruction with nonlocal regularizationIEEE Transactions on Medical Imaging, 2002
- SS3D-Fast fully 3D PET iterative reconstruction using stochastic samplingIEEE Transactions on Nuclear Science, 2002
- Comparison of post-filtering and filtering between iterations for SPECT reconstructionIEEE Transactions on Nuclear Science, 1999
- Multiresolution expectation maximization reconstruction algorithm for positron emission tomography using wavelet processingIEEE Transactions on Nuclear Science, 1999
- A row-action alternative to the EM algorithm for maximizing likelihood in emission tomographyIEEE Transactions on Medical Imaging, 1996
- Practical tradeoffs between noise, quantitation, and number of iterations for maximum likelihood-based reconstructionsIEEE Transactions on Medical Imaging, 1991
- Iterative methods for image deblurringProceedings of the IEEE, 1990
- Bayesian reconstructions from emission tomography data using a modified EM algorithmIEEE Transactions on Medical Imaging, 1990