DOI QR코드

DOI QR Code

Methods to Improve Convergence Rate of Statistical Reconstruction Algorithm in Transmission CT

투과형 CT에서 통계적 재구성 알고리즘의 수렴률 향상 방안

  • Min-Gu Song (Faculty of Liberal Arts, Yewon Arts University)
  • 송민구 (예원예술대학교 교양학부)
  • Received : 2024.04.27
  • Accepted : 2024.06.07
  • Published : 2024.06.30

Abstract

In tomographic image reconstruction, the focus is on developing CT image reconstruction methods that can maintain high image quality while reducing patient radiation exposure. Typically, statistical image reconstruction methods have the ability to generate high-quality and accurate images while significantly reducing patient radiation exposure. However, in cases like CT image reconstruction, which involve multi-dimensional parameter estimation, the degree of the Hessian matrix of the penalty function is very large, making it impossible to calculate. To solve this problem, the author proposed the PEMG-1 algorithm. However, the PEMG-1 algorithm has issues with the convergence speed, which is typical of statistical image reconstruction methods, and increasing the penalty log-likelihood. In this study, we propose a reconstruction algorithm that ensures fast convergence speed and monotonic increase in likelihood. The basic structure of this algorithm involves sequentially updating groups of pixels instead of updating all parameters simultaneously with each iteration.

토머그래피 영상재구성에서 초점은 높은 이미지 품질을 유지하면서 환자의 방사선 노출을 줄일 수 있는 CT 영상재구성 방법을 개발하는 것이다. 일반적으로 통계적 영상재구성 방법은 고품질 및 정확한 이미지를 생성할 수 있는 능력을 개선하면서 환자의 방사선 노출을 크게 줄일 수 있다. 그런데 CT 영상재구성과 같은 다차원의 모수 추정인 경우에서는 그것의 페널티 함수의 헤이지안 행렬의 역행렬 차수가 매우 크기 때문에 구할 수가 없다. 이러한 문제점을 해결하기 위하여 저자는 PEMG-1 알고리즘을 제안하였다. 그러나 PEMG-1 알고리즘은 일반 통계적 영상재구성 방법처럼 페널티 로그우도를 증가시키는 수렴속도에 문제점이 있다. 이에 본 연구에서는 수렴속도가 빠르고 우도의 단조 증가성을 보장하는 재구성 알고리즘을 제안한다. 이 알고리즘의 기본 구조는 반복마다 모수들을 동시에 갱신하지 않고 몇 개의 픽셀로 이루어진 그룹들을 순차적으로 갱신하는 방법이다.

Keywords

References

  1. Qiaoqiao Ding, Yong Long, and Jeffrey A. Fessler "Statistical Image Reconstruction Using Mixed Poisson Gaussian Noise Model for X-Ray CT" IEEE Transactions on Medical Imaging, Vol. 10, No. 3, pp.1-11, 2018.
  2. K. S. Lange. "A Gradient Algorithm Locally Equivalent to the EM Algorithm", Journal of Royal Statistical Society B, Vol.57, pp.425-437, 1995. https://doi.org/10.1111/j.2517-6161.1995.tb02037.x
  3. Tateishi, K, Yamaguchi, Y., Abou Al-Ola, O.M. and Yoshinaga, T. "Continuous Analog of Accelerated OS-EM Algorithm for Computed Tomography". Math. Probl. Eng. pp.34-42, 2017.
  4. Cho Z. W, John J. P. and Singh M.,"Foundation of Medical Imaging", John Wiley and Sons, 1993.
  5. A. P. Dempster. and N. M. Laird, and D. B. Rubin, "Maximum Likelihood from Incomplete Data from via the EM algorithm." Journal of Royal Statistical Society B, Vol.39, No.2. pp.1-38, 1977. https://doi.org/10.1111/j.2517-6161.1977.tb01600.x
  6. P. J. Green,."On use of the EM Algorithm for Penalized Likelihood Estimation", Journal of Royal Statistical SocietyB, Vol.52, pp.443-452, 1990. https://doi.org/10.1111/j.2517-6161.1990.tb01798.x
  7. Ing-Tsung Hsiao and Anand Rangarajan, "An Overview of Fast Convergent Ordered-Subsets Reconstruction Methods for Emission Tomography Based on the Incremental EM Algorithm" Nucl Instrum Methods Phys Res A. Vol.569, No.2, pp.429-422, 2006.
  8. M. K. Song. "A Study of Digital Image Restoration for Modified PEM Gradient Algorithm", The Transactions of the Korea Information Processing Society, Vol.7, No.1, pp.153-158, 2000.
  9. M. K. Song., "A Study of Penalty EM Gradient Algorithm in Digital Image Restoration," Doctoral Thesis., pp.1-81, 1998.
  10. J. A. Fessler. and H. Erdogan."A Paraboloidal Surrogates Algorithm for Convergent Penalized -likelihood Emission Image Reconstruction", IEEE, pp.1132-1135, 1999.
  11. J. A. Fessler. and E. P. Ficaro. & N. H. Clinthorne. & K. S. Lange, "Grouped-Coordinate Ascent Algorithms for Penalized-Likelihood Transmission Image Reconstruction", IEEE Transactions on Medical Imaging, Vol.16, No.2, pp.166-175, 1997. https://doi.org/10.1109/42.563662
  12. Omar M. Abou Al-Ola, Ryosuke and Kasai,Yusaku Yamaguchi, "Image Reconstruction Algorithm Using Weighted Mean of Ordered -Subsets EM and MART for Computed Tomography", Mathematics, Vol.10, No.22, pp.1-17, 2022. https://doi.org/10.3390/math10224277
  13. Sathish Ramani and Jeffrey A. Fessler, "A Splitting-Based Iterative Algorithm for Accelerated Statistical X-Ray CT Reconstruction", IEEE Transactions on Medical Imaging, Vol..31, No.3, pp. 677-688. ,2012. https://doi.org/10.1109/TMI.2011.2175233
  14. Sathish Ramani and Jeffrey A. Fessler,"A Hybrid Regularizer Combining Orthonormal Wavelets and Finite Differences for Statistical Reconstruction in 3-D CT, The second international conference on image formation in X-ray computed tomography, 2017.
  15. Sathish Ramani and Jeffrey A. Fessler, "Statistical X-Ray CT Reconstruction Using A Splitting-Based Iterative Algorithm with Orthonormal Wavelets, IEEE, pp.1008-1011, 2012.