Iterative Data Completion for Limited Angle Tomography using Filtered Backprojection

각도 제한 단층영상재구성을 위한 여현 역투사 기반 반복적 데이터 완결 기법

  • 이남용 (인제대학교 컴퓨터응용과학부)
  • Published : 2009.03.30

Abstract

When the range of projection angles is limited, tomographic reconstruction suffers from artifacts caused by incomplete data. One can consider a data completion technique, which estimates projection data at unobserved angles using a prior knowledge or mathematical exploration, but the result is often not improved; the improvement by the data completion often undermined by the artifacts by inaccurate estimation, In this paper, we propose an iterative method, which computes projection data at unobserved angles by using the current estimate on the image, links the computed projection data to the observed ones by using the consistence condition of Radon transform, and reconstruct the next estimate on the image by filtered backprojection. The proposed method does not require a prior knowledge on the image, and has much faster approximation rate than the expectation maximization method. The performance of the proposed method was tested through several simulation studies.

단층영상재구성에 있어서, 프로젝션 데이터를 제한된 각도에서만 관측할 수 있는 경우에는 불충분한 데이터에 의한 왜곡현상이 발생한다. 이러한 왜곡에 대한 해결책으로 선험적 정보 또는 수학적 외삽법 등을 이용하여 관측하지 않은 각도에서의 프로젝션 데이터를 완결하는 기법이 있으나, 부정확한 추정에 의한 새로운 왜곡이 데이터 완결에 의한 효과를 종종 상쇄한다. 본 논문에서는 각도 제한 단층영상재구성 방법으로 여현 역투사 기반 반복적 데이터 완결방법을 제안한다. 제안된 방법은 반복법으로, 현재의 영상으로부터 관측하지 않은 각도에서의 프로젝션 데이터를 계산하고 이를 관측한 각도에서의 프로젝션 데이터와 Radon 변환의 일관성 조건을 이용하여 결합한 후에, 여현 역투사에 적용하여 다음 단계 영상을 재구성한다. 제안된 방법은 선험적 정보를 필요로 하지 않는다는 것과 함께 기대치 최대화에 비해 그 수렴속도가 현저히 빠르다는 장점이 있다. 제안된 방법의 성능은 모의실험을 통해 확인하였다.

Keywords

References

  1. A.C. Kak and M. Slaney, Principles of Computerized Tomography Imaging, IEEE Press, New York, 1987.
  2. L.A. Shepp and B.F. Logan, "The Fourier reconstruction of a head section," IEEE Trans. Nuclear Sci., Vol.21, pp. 21-43, 1974.
  3. L.A. Shepp and Y. Vardi, "Maximum likel-ihood reconstruction for emission tomography," IEEE Trans. Medical Imaging, Vol.1, pp. 113-122, 1982. https://doi.org/10.1109/TMI.1982.4307558
  4. M.E. Davision, "The ill-conditioned nature of the limited angle tomography problem," SIAM J. Appl. Math. Vol.43, pp. 428-448, 1983. https://doi.org/10.1137/0143028
  5. A. Ramrn, "Inversion of limited-angle tomo-graphic data," Comput. Math Appl., Vol.22, pp, 101-111, 1992.
  6. F. Natterer and F. Wubbeling, Mathematical methods in image reconstruction, SIAM, Philadelphia, 2001.
  7. R.M. Lewitt, "Image reconstruction from projections III: Projection completion meth-ods(theory)," Optik, Vol.50, pp, 180-205, 1978.
  8. F. Natterer, The mathematics of computerized tomography, John Wiley, B.G. Teubner, New York, 1986.
  9. A.K. Louis and A. Reider, "Incomplete data problems in x-ray computerized tomography II: Truncated projections and region-of-interest tomography," Numer. Math. Vol.56, pp.371-383, 1989. https://doi.org/10.1007/BF01396611
  10. G.L. Zeng, G.T. Gullberg, B.M.W. Tsui, and J.A. Terry, "Three-dimensional iterative reconstruction algorithms with attenuation and geometric point response correction," IEEE Trans. Nuclear Science, Vol.38, No.2, pp. 693-702, 1991. https://doi.org/10.1109/23.289376
  11. J. Y. Hesterman and M. A. Kupinski, "The multi-module, multi-resolution system (M3R): A novel small-animal SPECT system," Med. Phys., Vol.34, Issue 3, pp. 987-993, 2007. https://doi.org/10.1118/1.2432071
  12. G. Wang, L. Fu, and J. Qi, "Maximum a posteriori reconstruction of the Patlak parametric image from sinograms in dynamic PET," Phys. Med. Biol., Vol.53, pp. 593-604, 2008. https://doi.org/10.1088/0031-9155/53/3/006
  13. W.H. Press, S.A. Teukolsky, W.T. Vettering, and B.P. Flannery, Numerical Recipes in C, University Press, Cambridge, 1992.