Numerical deconvolution by a Monte Carlo approach with application to dynamic cardiac perfusion Tc-99m SPECT

  • Authors:

  • Koon-Pong Wong;Steven R. Meikle;David Dagan Feng;Michael J. Fulham

  • Affiliations:

  • Dept. of PET and Nuclear Medicine, Royal Prince Alfred Hospital, Camperdown, Australia and Biomedical and Multimedia Info. Tech. Gp., Basser Dept. of Comp. Sci., Univ. Sydney, Australia;Dept. of PET and Nuclear Medicine, Royal Prince Alfred Hospital, Camperdown, Australia and Biomedical and Multimedia Info. Tech. Gp., Basser Dept. of Comp. Sci., Univ. Sydney, Australia;Biomedical and Multimedia Info. Tech. Gp., Basser Dept. of Comp. Sci., Univ. Sydney, Australia and Center for Multimedia Signal Proc., Dept. of Elec. and Info. Eng., Hong Kong Poly. Univ., Hong Ko ...;Dept. of PET and Nuclear Medicine, Royal Prince Alfred Hospital, Camperdown, Australia and Biomedical and Multimedia Info. Tech. Gp., Basser Dept. of Comp. Sci., Univ. of Sydney, Australia and Fac ...

  • Venue:
  • VIP '01 Proceedings of the Pan-Sydney area workshop on Visual information processing - Volume 11
  • Year:
  • 2001

Quantified Score

Hi-index 0.00

Visualization

Abstract

Accurate determination of the input function is essential for absolute quantification of physiological parameters in PET and SPECT imaging but it requires an invasive and tedious procedure of blood sampling that is impractical in clinical studies. We previously proposed a technique that simultaneously estimates kinetic parameters and the input function from the tissue impulse response functions and which, requires only two blood samples. A nonlinear least squares method was used to estimate all the parameters in the impulse response functions and the input function but it fails occasionally due to high noise levels in the data causing an ill-conditioned cost function. This study investigates the feasibility of applying a Monte Carlo method called simulated annealing to estimate kinetic parameters in the impulse response functions and the input function. Time-activity curves of teboroxime, which is very sensitive to changes in the input function, were simulated based on published data obtained from a canine model. The equations describing the tracer kinetics in different regions were minimised simultaneously by simulated annealing and nonlinear least squares. We found that the physiological parameters obtained with simulated annealing are more accurate and the estimated input function more closely resembled the simulated curve. We conclude that simulated annealing reduces bias in the estimation of physiological parameters and determination of the input function.