Tehran University of Medical Sciences

Science Communicator Platform

Stay connected! Follow us on X network (Twitter):
Share this content! On (X network) By
Computationally Efficient System Matrix Calculation Techniques in Computed Tomography Iterative Reconstruction Publisher



Mahmoudi G1, 2 ; Ay M1, 2 ; Rahmim A3, 4 ; Ghadiri H1, 2
Authors
Show Affiliations
Authors Affiliations
  1. 1. Department of Medical Physics and Biomedical Engineering, Tehran University of Medical Sciences, Tehran, Iran
  2. 2. Research Center for Molecular and Cellular Imaging, Tehran University of Medical Sciences, Tehran, Iran
  3. 3. Department of Radiology and Physics, University of British Columbia, Vancouver, BC, Canada
  4. 4. Department of Integrative Oncology, BC Cancer Research Centre, Vancouver, BC, Canada

Source: Journal of Medical Signals and Sensors Published:2020


Abstract

Background: Relative to classical methods in computed tomography, iterative reconstruction techniques enable significantly improved image qualities and/or lowered patient doses. However, the computational speed is a major concern for these iterative techniques. In the present study, we present a method for fast system matrix calculation based on the line integral model (LIM) to speed up the computations without compromising the image quality. In addition, we develop a hybrid line-area integral model (AIM) that highlights the advantages of both LIM and AIMs. Methods: The contributing detectors for a given pixel and a given projection view, and the length of corresponding intersection lines with pixels, are calculated using our proposed algorithm. For the hybrid method, the respective narrow-angle fan beam was modeled by multiple equally spaced lines. The computed system matrix was evaluated in the context of reconstruction using the simultaneous algebraic reconstruction technique (SART) as well as maximum likelihood expectation maximization (MLEM). Results: The proposed LIM offers a considerable reduction in calculation times compared to the standard Siddon algorithm: 2.9 times faster. Differences in root mean square error and peak signal-to-noise ratio were not significant between the proposed LIM and the Siddon algorithm for both SART and MLEM reconstruction methods (P > 0.05). Meanwhile, the proposed hybrid method resulted in significantly improved image qualities relative to LIM and the Siddon algorithm (P < 0.05), though computations were 4.9 times more intensive than the proposed LIM. Conclusion: We have proposed two fast algorithms to calculate the system matrix. The first is based on LIM and was faster than the Siddon algorithm, with matched image quality, whereas the second method is a hybrid LIM-AIM that achieves significantly improved images though with its computational requirements. © 2020 Isfahan University of Medical Sciences(IUMS). All rights reserved.