Pre-computed system matrix calculation based on a piece-wise method for PET

Abdella M. Ahmed, Yohei Kikuchi, Shigeo Matsuyama, Atsuki Terakawa, Sodai Takyu, Hiroyuki Sugai, Keizo Ishii

Research output: Contribution to journalArticlepeer-review

1 Citation (Scopus)


A system matrix (SM) is the basic component of iterative image reconstruction algorithms. Calculation of the SM needs a considerable amount of time due to an enormous number of lines of response (LORs) being modeled. In this study, we developed a technique based on a piece-wise calculation method in which symmetry and further division of the voxels are applied. The detector response function for all detectable pairs of photons along certain LORs originating from each voxel is calculated analytically. The total number of LORs in 300 × 300 × 120 voxels (with 2 × 2 × 2 mm3) is ~44 billion, and the SM was calculated by the use of three different computers independently; the calculation time was 5 h. The SM took 5 days when calculated by the use of the conventional method (where symmetry and the piece-wise method are not used). The sensitivity correction factor was stored; it had a size of 42 MB in a four-byte computer memory.

Original languageEnglish
Pages (from-to)88-96
Number of pages9
JournalRadiological Physics and Technology
Issue number1
Publication statusPublished - 2015 Jan


  • Image reconstruction algorithm
  • PET
  • System matrix

ASJC Scopus subject areas

  • Radiation
  • Physical Therapy, Sports Therapy and Rehabilitation
  • Radiology Nuclear Medicine and imaging


Dive into the research topics of 'Pre-computed system matrix calculation based on a piece-wise method for PET'. Together they form a unique fingerprint.

Cite this