Technical Note: Rapid multiexponential curve fitting algorithm for voxel-based targeted radionuclide dosimetry

Jackson, Price, McIntosh, Lachlan, Hofman, Michael S., Kong, Grace, and Hicks, Rodney J. (2020) Technical Note: Rapid multiexponential curve fitting algorithm for voxel-based targeted radionuclide dosimetry. Medical Physics, 47 (9). pp. 4332-4339.

PDF (Accepted Author Version) - Accepted Version
Download (533kB) | Preview
View at Publisher Website:


Background: Dosimetry in nuclear medicine often relies on estimating pharmacokinetics based on sparse temporal data. As analysis methods move toward image-based three-dimensional computation, it becomes important to interpolate and extrapolate these data without requiring manual intervention; that is, in a manner that is highly efficient and reproducible. Iterative least-squares solvers are poorly suited to this task because of the computational overhead and potential to optimize to local minima without applying tight constraints at the outset.

Methodology: This work describes a fully analytical method for solving three-phase exponential time-activity curves based on three measured time points in a manner that may be readily employed by image-based dosimetry tools. The methodology uses a series of conditional statements and a piecewise approach for solving exponential slope directly through measured values in most instances. The proposed algorithm is tested against a purpose-designed iterative fitting technique and linear piecewise method followed by single exponential in a cohort of ten patients receiving 177Lu-DOTA-Octreotate therapy.

Results: Tri-exponential time-integrated values are shown to be comparable to previously published methods with an average difference between organs when computed at the voxel level of 9.8 ± 14.2% and −3.6 ± 10.4% compared to iterative and interpolated methods, respectively. Of the three methods, the proposed tri-exponential algorithm was most consistent when regional time-integrated activity was evaluated at both voxel- and whole-organ levels. For whole-body SPECT imaging, it is possible to compute 3D time-integrated activity maps in '5 min processing time. Furthermore, the technique is able to predictably and reproducibly handle artefactual measurements due to noise or spatial misalignment over multiple image times.

Conclusions: An efficient, analytical algorithm for solving multiphase exponential pharmacokinetics is reported. The method may be readily incorporated into voxel-dose routines by combining with widely available image registration and radiation transport tools.

Item ID: 75274
Item Type: Article (Research - C1)
ISSN: 2473-4209
Keywords: image processing, pharmacokinetics, radionuclide dosimetry
Copyright Information: Published Version: © Wiley-Blackwell. Author Accepted Version may be made open access in an Institutional Repository after a 12 month embargo.
Funders: National Health and Medical Research Council of Australia (NHMRC)
Projects and Grants: NHMRC 1108050
Date Deposited: 17 Aug 2022 00:51
FoR Codes: 32 BIOMEDICAL AND CLINICAL SCIENCES > 3202 Clinical sciences > 320215 Nuclear medicine @ 50%
51 PHYSICAL SCIENCES > 5105 Medical and biological physics > 510502 Medical physics @ 50%
Downloads: Total: 677
Last 12 Months: 82
More Statistics

Actions (Repository Staff Only)

Item Control Page Item Control Page