Time-efficient Fourier domain evaluation of pharmacokinetic model in dynamic contrast-enhanced magnetic resonance imaging

0490728 - ÚTIA 2020 RIV SG eng C - Conference Paper (international conference)
Bartoš, Michal - Šorel, Michal - Jiřík, Radovan
Time-efficient Fourier domain evaluation of pharmacokinetic model in dynamic contrast-enhanced magnetic resonance imaging.
IFMBE Proceedings, Volume 68, Issue 1. Singapore: Springer, 2019, s. 777-781. ISBN 978-981-10-9034-9. ISSN 1680-0737.
[World Congress on Medical Physics and Biomedical Engineering. Praha (CZ), 03.06.2018-08.06.2018]
R&D Projects: GA ČR(CZ) GA16-13830S
Grant - others:AV ČR(CZ) MSM100751802
Program: Program na podporu mezinárodní spolupráce začínajících výzkumných pracovníků
Institutional support: RVO:67985556 ; RVO:68081731
Keywords : DCE-MRI * Tissue homogeneity model * Tracer kinetic modelling
OECD category: Computer sciences, information science, bioinformathics (hardware development to be 2.2, social aspect to be 5.8); Computer sciences, information science, bioinformathics (hardware development to be 2.2, social aspect to be 5.8) (UPT-D)
http://library.utia.cas.cz/separaty/2018/ZOI/bartos-0490728.pdf

Dynamic contrast-enhanced magnetic resonance imaging obtains information about tissue perfusion and permeability. Following the administration of a contrast agent, concentration-time curves measured in each voxel are fitted by a pharmacokinetic model formulated as a time-domain convolution of an arterial input function (AIF) and an impulse residue function (IRF). Since the measurement window contains hundreds of time samples, the discrete convolution is demanding, even when it is performed via discrete Fourier transform (DFT). Additionally, its discretization causes convergence complications in the curve fitting and it is not applicable to functions without a closed-form expression in the time domain, e.g. tissue homogeneity model IRF. Both issues can be solved by formulating the functions in a closed form in the Fourier domain. In the Fourier domain, the model transforms to multiplication of IRF and AIF, followed by the inverse DFT. To avoid time-domain aliasing, the number of samples in the Fourier domain must be higher than the sum of supports of the functions in the time domain. If the functions are slowly decaying exponentials, the support is theoretically infinite, which dramatically reduces the computational performance. In this contribution, we propose a modification of IRF in the Fourier domain to consider the measurement window. Our solution reduces the required number of samples to three times the measurement window compared to dozens needed without the modification and reduces the number of DFTs. This provides faster evaluation of the pharmacokinetic model and its derivatives for each voxel in each iteration of the curve fitting.
Permanent Link: http://hdl.handle.net/11104/0285273