Hierarchical Bayesian myocardial perfusion quantification

Cian M. Scannell (Corresponding author), Amedeo Chiribiri, Adriana D. M. Villa, Marcel Breeuwer, Jack Lee

Research output: Contribution to journalArticleAcademicpeer-review

12 Downloads (Pure)

Abstract

Purpose: Tracer-kinetic models can be used for the quantitative assessment of contrast-enhanced MRI data. However, the model-fitting can produce unreliable results due to the limited data acquired and the high noise levels. Such problems are especially prevalent in myocardial perfusion MRI leading to the compromise of constrained numerical deconvolutions and segmental signal averaging being commonly used as alternatives to the more complex tracer-kinetic models. Methods: In this work, the use of hierarchical Bayesian inference for the parameter estimation is explored. It is shown that with Bayesian inference it is possible to reliably fit the two-compartment exchange model to perfusion data. The use of prior knowledge on the ranges of kinetic parameters and the fact that neighbouring voxels are likely to have similar kinetic properties combined with a Markov chain Monte Carlo based fitting procedure significantly improves the reliability of the perfusion estimates with compared to the traditional least-squares approach. The method is assessed using both simulated and patient data. Results: The average (standard deviation) normalised mean square error for the distinct noise realisations of a simulation phantom falls from 0.32 (0.55) with the least-squares fitting to 0.13 (0.2) using Bayesian inference. The assessment of the presence of coronary artery disease based purely on the quantitative MBF maps obtained using Bayesian inference matches the visual assessment in all 24 slices. When using the maps obtained by the least-squares fitting, a corresponding assessment is only achieved in 16/24 slices. Conclusion: Bayesian inference allows a reliable, fully automated and user-independent assessment of myocardial perfusion on a voxel-wise level using the two-compartment exchange model.
Original languageEnglish
Article number1906.02540
Number of pages37
JournalarXiv
Volume2019
Publication statusPublished - 6 Jun 2019

Fingerprint

Perfusion
Least-Squares Analysis
Noise
Markov Chains
Coronary Artery Disease

Cite this

Scannell, C. M., Chiribiri, A., Villa, A. D. M., Breeuwer, M., & Lee, J. (2019). Hierarchical Bayesian myocardial perfusion quantification. arXiv, 2019, [1906.02540].
Scannell, Cian M. ; Chiribiri, Amedeo ; Villa, Adriana D. M. ; Breeuwer, Marcel ; Lee, Jack. / Hierarchical Bayesian myocardial perfusion quantification. In: arXiv. 2019 ; Vol. 2019.
@article{822a766f9e4f4057b5e736728e0ec6de,
title = "Hierarchical Bayesian myocardial perfusion quantification",
abstract = "Purpose: Tracer-kinetic models can be used for the quantitative assessment of contrast-enhanced MRI data. However, the model-fitting can produce unreliable results due to the limited data acquired and the high noise levels. Such problems are especially prevalent in myocardial perfusion MRI leading to the compromise of constrained numerical deconvolutions and segmental signal averaging being commonly used as alternatives to the more complex tracer-kinetic models. Methods: In this work, the use of hierarchical Bayesian inference for the parameter estimation is explored. It is shown that with Bayesian inference it is possible to reliably fit the two-compartment exchange model to perfusion data. The use of prior knowledge on the ranges of kinetic parameters and the fact that neighbouring voxels are likely to have similar kinetic properties combined with a Markov chain Monte Carlo based fitting procedure significantly improves the reliability of the perfusion estimates with compared to the traditional least-squares approach. The method is assessed using both simulated and patient data. Results: The average (standard deviation) normalised mean square error for the distinct noise realisations of a simulation phantom falls from 0.32 (0.55) with the least-squares fitting to 0.13 (0.2) using Bayesian inference. The assessment of the presence of coronary artery disease based purely on the quantitative MBF maps obtained using Bayesian inference matches the visual assessment in all 24 slices. When using the maps obtained by the least-squares fitting, a corresponding assessment is only achieved in 16/24 slices. Conclusion: Bayesian inference allows a reliable, fully automated and user-independent assessment of myocardial perfusion on a voxel-wise level using the two-compartment exchange model.",
keywords = "eess.IV, cs.CV, physics.med-ph",
author = "Scannell, {Cian M.} and Amedeo Chiribiri and Villa, {Adriana D. M.} and Marcel Breeuwer and Jack Lee",
year = "2019",
month = "6",
day = "6",
language = "English",
volume = "2019",
journal = "arXiv",
publisher = "Cornell University Library",

}

Scannell, CM, Chiribiri, A, Villa, ADM, Breeuwer, M & Lee, J 2019, 'Hierarchical Bayesian myocardial perfusion quantification', arXiv, vol. 2019, 1906.02540.

Hierarchical Bayesian myocardial perfusion quantification. / Scannell, Cian M. (Corresponding author); Chiribiri, Amedeo; Villa, Adriana D. M.; Breeuwer, Marcel; Lee, Jack.

In: arXiv, Vol. 2019, 1906.02540, 06.06.2019.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Hierarchical Bayesian myocardial perfusion quantification

AU - Scannell, Cian M.

AU - Chiribiri, Amedeo

AU - Villa, Adriana D. M.

AU - Breeuwer, Marcel

AU - Lee, Jack

PY - 2019/6/6

Y1 - 2019/6/6

N2 - Purpose: Tracer-kinetic models can be used for the quantitative assessment of contrast-enhanced MRI data. However, the model-fitting can produce unreliable results due to the limited data acquired and the high noise levels. Such problems are especially prevalent in myocardial perfusion MRI leading to the compromise of constrained numerical deconvolutions and segmental signal averaging being commonly used as alternatives to the more complex tracer-kinetic models. Methods: In this work, the use of hierarchical Bayesian inference for the parameter estimation is explored. It is shown that with Bayesian inference it is possible to reliably fit the two-compartment exchange model to perfusion data. The use of prior knowledge on the ranges of kinetic parameters and the fact that neighbouring voxels are likely to have similar kinetic properties combined with a Markov chain Monte Carlo based fitting procedure significantly improves the reliability of the perfusion estimates with compared to the traditional least-squares approach. The method is assessed using both simulated and patient data. Results: The average (standard deviation) normalised mean square error for the distinct noise realisations of a simulation phantom falls from 0.32 (0.55) with the least-squares fitting to 0.13 (0.2) using Bayesian inference. The assessment of the presence of coronary artery disease based purely on the quantitative MBF maps obtained using Bayesian inference matches the visual assessment in all 24 slices. When using the maps obtained by the least-squares fitting, a corresponding assessment is only achieved in 16/24 slices. Conclusion: Bayesian inference allows a reliable, fully automated and user-independent assessment of myocardial perfusion on a voxel-wise level using the two-compartment exchange model.

AB - Purpose: Tracer-kinetic models can be used for the quantitative assessment of contrast-enhanced MRI data. However, the model-fitting can produce unreliable results due to the limited data acquired and the high noise levels. Such problems are especially prevalent in myocardial perfusion MRI leading to the compromise of constrained numerical deconvolutions and segmental signal averaging being commonly used as alternatives to the more complex tracer-kinetic models. Methods: In this work, the use of hierarchical Bayesian inference for the parameter estimation is explored. It is shown that with Bayesian inference it is possible to reliably fit the two-compartment exchange model to perfusion data. The use of prior knowledge on the ranges of kinetic parameters and the fact that neighbouring voxels are likely to have similar kinetic properties combined with a Markov chain Monte Carlo based fitting procedure significantly improves the reliability of the perfusion estimates with compared to the traditional least-squares approach. The method is assessed using both simulated and patient data. Results: The average (standard deviation) normalised mean square error for the distinct noise realisations of a simulation phantom falls from 0.32 (0.55) with the least-squares fitting to 0.13 (0.2) using Bayesian inference. The assessment of the presence of coronary artery disease based purely on the quantitative MBF maps obtained using Bayesian inference matches the visual assessment in all 24 slices. When using the maps obtained by the least-squares fitting, a corresponding assessment is only achieved in 16/24 slices. Conclusion: Bayesian inference allows a reliable, fully automated and user-independent assessment of myocardial perfusion on a voxel-wise level using the two-compartment exchange model.

KW - eess.IV

KW - cs.CV

KW - physics.med-ph

M3 - Article

VL - 2019

JO - arXiv

JF - arXiv

M1 - 1906.02540

ER -

Scannell CM, Chiribiri A, Villa ADM, Breeuwer M, Lee J. Hierarchical Bayesian myocardial perfusion quantification. arXiv. 2019 Jun 6;2019. 1906.02540.