Umeå University's logo

umu.sePublications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Bayesian non-linear regression with spatial priors for noise reduction and error estimation in quantitative MRI with an application in T1 estimation
Umeå University, Faculty of Medicine, Department of Radiation Sciences, Radiation Physics. Umeå University, Faculty of Science and Technology, Department of Computing Science.ORCID iD: 0000-0001-7119-7646
Umeå University, Faculty of Medicine, Department of Radiation Sciences, Radiation Physics.ORCID iD: 0000-0002-0200-6567
Umeå University, Faculty of Medicine, Department of Radiation Sciences, Radiation Physics.ORCID iD: 0000-0001-7539-2262
Umeå University, Faculty of Medicine, Department of Radiation Sciences, Radiation Physics.ORCID iD: 0000-0002-0532-232X
2020 (English)In: Physics in Medicine and Biology, ISSN 0031-9155, E-ISSN 1361-6560, Vol. 65, no 22, article id 225036Article in journal (Refereed) Published
Abstract [en]

Purpose. To develop a method that can reduce and estimate uncertainty in quantitative MR parameter maps without the need for hand-tuning of any hyperparameters.

Methods. We present an estimation method where uncertainties are reduced by incorporating information on spatial correlations between neighbouring voxels. The method is based on a Bayesian hierarchical non-linear regression model, where the parameters of interest are sampled, using Markov chain Monte Carlo (MCMC), from a high-dimensional posterior distribution with a spatial prior. The degree to which the prior affects the model is determined by an automatic hyperparameter search using an information criterion and is, therefore, free from manual user-dependent tuning. The samples obtained further provide a convenient means to obtain uncertainties in both voxels and regions. The developed method was evaluated on T1 estimations based on the variable flip angle method.

Results. The proposed method delivers noise-reduced T1 parameter maps with associated error estimates by combining MCMC sampling, the widely applicable information criterion, and total variation-based denoising. The proposed method results in an overall decrease in estimation error when compared to conventional voxel-wise maximum likelihood estimation. However, this comes with an increased bias in some regions, predominately at tissue interfaces, as well as an increase in computational time.

Conclusions. This study provides a method that generates more precise estimates compared to the conventional method, without incorporating user subjectivity, and with the added benefit of uncertainty estimation.

Place, publisher, year, edition, pages
Institute of Physics (IOP), 2020. Vol. 65, no 22, article id 225036
Keywords [en]
Bayesian statistics, quantitative MRI, noise reduction, tissue parameter estimation, WAIC
National Category
Radiology, Nuclear Medicine and Medical Imaging Computer graphics and computer vision Probability Theory and Statistics
Identifiers
URN: urn:nbn:se:umu:diva-177258DOI: 10.1088/1361-6560/abb9f5ISI: 000591796100001PubMedID: 32947277Scopus ID: 2-s2.0-85097228803OAI: oai:DiVA.org:umu-177258DiVA, id: diva2:1506404
Funder
Swedish Research Council, 2019-0432Region Västerbotten, RV-738491Available from: 2020-12-03 Created: 2020-12-03 Last updated: 2025-09-09Bibliographically approved
In thesis
1. Noise reduction and uncertainty estimation in quantitative magnetic resonance imaging
Open this publication in new window or tab >>Noise reduction and uncertainty estimation in quantitative magnetic resonance imaging
2025 (English)Doctoral thesis, comprehensive summary (Other academic)
Alternative title[sv]
Brusreducering och osäkerhetsuppskattning i kvantitativ magnetresonanstomografi
Abstract [en]

Quantitative magnetic resonance imaging (qMRI) is a subset of MRI that moves beyond qualitative image interpretation towards more objective and reproducible quantification of tissue parameters, a process referred to as parameter mapping. In cancer imaging, this capability has the potential to assist with characterizing tumors, monitoring treatment response, and comparing scans both over time and across different scanners or imaging sites. This thesis addresses the problem of uncertainty in parameter mapping. This is done by developing methods that both reduce uncertainty and estimate the uncertainty that remains, two important components for advancing the clinical reliability of qMRI. The first component, reducing uncertainty, focuses on reducing noise in the parameter maps. In this thesis, we propose methods that apply denoising priors to stabilize the parameter mapping process. The second component, estimating the remaining uncertainty, aims to quantify what uncertainty remains after denoising has been applied. This is achieved through uncertainty estimation, providing a more complete picture of the results.

Paper I presents a method for T1 mapping that combines an explicit denoising prior, in the form of total variation (TV) regularization, together with uncertainty estimation using Markov chain Monte Carlo (MCMC) sampling. The strength of the TV prior is selected in a data-driven way which reduces the requirement of manual tuning. This approach produces denoised parameter maps with reduced uncertainty and meaningful uncertainty estimates. Remaining challenges include the computational burden and the need to carefully manage prior-induced biases.

Paper II introduces an alternative to the method in Paper I. Instead of using an explicit prior, it exploits Deep Image Prior (DIP), which utilizes the structure of an untrained convolutional neural network to impose an implicit prior. This is implemented with Monte Carlo dropout, a simple approach to approximate Bayesian inference for uncertainty estimation. The method is easy to implement and denoises the tissue parameters T1, T2, and ADC. However, the results reveal that uncertainty calibration becomes more challenging, and the computation time remains too long for practical clinical use.

Paper III improves the DIP method from Paper II by addressing long computation times and difficulties in uncertainty calibration. To reduce computation times, we introduce warm-start initialization, which leverages information from both neighboring image slices and different subjects to accelerate parameter mapping. To improve calibration, we systematically tune MC dropout to reduce miscalibration. In addition, a data-driven early stopping criterion is proposed to automatically set the denoising level, removing the need for manual tuning. Together, these changes make the DIP method faster, better calibrated, and more clinically usable.

Paper IV investigates further improvements of the DIP method towards more advanced qMRI models, specifically evaluated with Patlak-based estimation of pharmacokinetic parameters in DCE-MRI. This study evaluates the feasibility of this approach by comparing parameter maps with and without the DIP method applied. Preliminary results show substantial noise reduction and improved feature representation. Remaining challenges include prior-induced biases that require further refinement.

In summary, this thesis addresses uncertainty in qMRI by focusing on two key components: reducing uncertainty through denoising and estimating the remaining uncertainty. These aspects are investigated using two fundamentally different strategies: one based on explicit denoising with robust Bayesian inference using MCMC (Paper I), and one based on implicit denoising combined with approximate Bayesian inference (Papers II-IV). The four included studies develop and refine these strategies to overcome practical limitations and enhance applicability. By evaluating their respective strengths and weaknesses, the work provides insight into how methodological choices affect accuracy, uncertainty, and computational feasibility. These contributions aim to support the development of qMRI as a more trustworthy and reproducible tool in cancer imaging and beyond.

Place, publisher, year, edition, pages
Umeå: Umeå University, 2025. p. 78
Series
Umeå University medical dissertations, ISSN 0346-6612 ; 2376
Keywords
Quantitative MRI, Denoising, Uncertainty estimation, Markov chain Monte Carlo, Deep Image Prior
National Category
Radiology and Medical Imaging
Research subject
radiation physics
Identifiers
urn:nbn:se:umu:diva-244020 (URN)978-91-8070-760-2 (ISBN)978-91-8070-759-6 (ISBN)
Public defence
2025-10-03, Triple Helix, Universitetsledningshuset (plan 3), Umeå, 09:00 (English)
Opponent
Supervisors
Available from: 2025-09-12 Created: 2025-09-09 Last updated: 2025-09-11Bibliographically approved

Open Access in DiVA

fulltext(3367 kB)401 downloads
File information
File name FULLTEXT01.pdfFile size 3367 kBChecksum SHA-512
6c1b99c2006e46caa0c387d8666e23122cb8c831331aa6b2e2b83c639a0e131034b3a28e07fcbfa83a05dd344cb677cf711ac6194e94d31e7cccb103843eb99a
Type fulltextMimetype application/pdf

Other links

Publisher's full textPubMedScopus

Authority records

Löfstedt, TommyHellström, MaxBylund, MikaelGarpebring, Anders

Search in DiVA

By author/editor
Löfstedt, TommyHellström, MaxBylund, MikaelGarpebring, Anders
By organisation
Radiation PhysicsDepartment of Computing Science
In the same journal
Physics in Medicine and Biology
Radiology, Nuclear Medicine and Medical ImagingComputer graphics and computer visionProbability Theory and Statistics

Search outside of DiVA

GoogleGoogle Scholar
Total: 401 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

doi
pubmed
urn-nbn

Altmetric score

doi
pubmed
urn-nbn
Total: 821 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf