This paper analyzes in detail the process of tomographic image reconstruction by pseudo-inversion of the blurring matrix of a PET imaging system. Eigenvector and eigenvalue decomposition is used as a method to evaluate the physical reasons for the ill-conditioned nature of the problem. It is shown that finding an accurate pseudo-inverse for even a modest PET array of 8 × 8 pixels is a difficult task for a computer with 48-bit mantissa. The problem is caused by the strong ambiguity with which the detector system measures the activity at each pixel. For a problem in which imaging with a complete detector ring is not possible, and in which invariance of the point response function cannot be maintained, the pseudo-inverse method of reconstruction is, however, shown to be very useful. Advantage is taken of the fact that the activitiy to be measured is localized in a single plane, without over-or underlying activity. A planar camera configuration yields very well conditioned matrices that are separable for a large number of useful cases. It is even possible to define pixel sizes which are considerably smaller than the detector size and solve the problem without a substantial increase in the noise magnification factor. Recognizing that the above application is equivalent to a case of very well defined time-of-flight (TOF) measurement, the simple initial PET study is reevaluated by inclusion of TOF information.
[1]
G. G. Stokes.
"J."
,
1890,
The New Yale Book of Quotations.
[2]
B. R. Hunt,et al.
Digital Image Restoration
,
1977
.
[3]
E. Paul.
Nuclear Science
,
1957,
Nature.
[4]
Jorge Llacer.
Theory of Imaging with a Very Limited Number of Projections
,
1979,
IEEE Transactions on Nuclear Science.
[5]
Z. Cho,et al.
General views on 3-D image reconstruction and computerized transverse axial tomography
,
1974
.
[6]
C A Tobias,et al.
High energy beams of radioactive nuclei and their biomedical applications.
,
1981,
International journal of radiation oncology, biology, physics.