{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,14]],"date-time":"2026-03-14T05:25:26Z","timestamp":1773465926085,"version":"3.50.1"},"reference-count":38,"publisher":"MDPI AG","issue":"2","license":[{"start":{"date-parts":[[2021,2,20]],"date-time":"2021-02-20T00:00:00Z","timestamp":1613779200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100002341","name":"Academy of Finland","doi-asserted-by":"publisher","award":["312343"],"award-info":[{"award-number":["312343"]}],"id":[{"id":"10.13039\/501100002341","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100002341","name":"Academy of Finland","doi-asserted-by":"publisher","award":["336791"],"award-info":[{"award-number":["336791"]}],"id":[{"id":"10.13039\/501100002341","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100016032","name":"V\u00e4is\u00e4l\u00e4n Rahasto","doi-asserted-by":"publisher","award":["-"],"award-info":[{"award-number":["-"]}],"id":[{"id":"10.13039\/501100016032","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100004012","name":"Jane ja Aatos Erkon S\u00e4\u00e4ti\u00f6","doi-asserted-by":"publisher","award":["-"],"award-info":[{"award-number":["-"]}],"id":[{"id":"10.13039\/501100004012","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["J. Imaging"],"abstract":"<jats:p>In dynamic MRI, sufficient temporal resolution can often only be obtained using imaging protocols which produce undersampled data for each image in the time series. This has led to the popularity of compressed sensing (CS) based reconstructions. One problem in CS approaches is determining the regularization parameters, which control the balance between data fidelity and regularization. We propose a data-driven approach for the total variation regularization parameter selection, where reconstructions yield expected sparsity levels in the regularization domains. The expected sparsity levels are obtained from the measurement data for temporal regularization and from a reference image for spatial regularization. Two formulations are proposed. Simultaneous search for a parameter pair yielding expected sparsity in both domains (S-surface), and a sequential parameter selection using the S-curve method (Sequential S-curve). The approaches are evaluated using simulated and experimental DCE-MRI. In the simulated test case, both methods produce a parameter pair and reconstruction that is close to the root mean square error (RMSE) optimal pair and reconstruction. In the experimental test case, the methods produce almost equal parameter selection, and the reconstructions are of high perceived quality. Both methods lead to a highly feasible selection of the regularization parameters in both test cases while the sequential method is computationally more efficient.<\/jats:p>","DOI":"10.3390\/jimaging7020038","type":"journal-article","created":{"date-parts":[[2021,2,21]],"date-time":"2021-02-21T21:15:01Z","timestamp":1613942101000},"page":"38","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":4,"title":["Data-Driven Regularization Parameter Selection in Dynamic MRI"],"prefix":"10.3390","volume":"7","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-5044-3203","authenticated-orcid":false,"given":"Matti","family":"Hanhela","sequence":"first","affiliation":[{"name":"Department of Applied Physics, University of Eastern Finland, 70211 Kuopio, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-1372-1651","authenticated-orcid":false,"given":"Olli","family":"Gr\u00f6hn","sequence":"additional","affiliation":[{"name":"A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, 70211 Kuopio, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2004-660X","authenticated-orcid":false,"given":"Mikko","family":"Kettunen","sequence":"additional","affiliation":[{"name":"A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, 70211 Kuopio, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Kati","family":"Niinim\u00e4ki","sequence":"additional","affiliation":[{"name":"Xray Division, Planmeca Oy, Asentajankatu 6, 00880 Helsinki, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6454-8524","authenticated-orcid":false,"given":"Marko","family":"Vauhkonen","sequence":"additional","affiliation":[{"name":"Department of Applied Physics, University of Eastern Finland, 70211 Kuopio, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5621-795X","authenticated-orcid":false,"given":"Ville","family":"Kolehmainen","sequence":"additional","affiliation":[{"name":"Department of Applied Physics, University of Eastern Finland, 70211 Kuopio, Finland"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2021,2,20]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"473","DOI":"10.1109\/42.97598","article-title":"Selection of a convolution function for Fourier inversion using gridding (computerised tomography application)","volume":"10","author":"Jackson","year":"1991","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"489","DOI":"10.1109\/TIT.2005.862083","article-title":"Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information","volume":"52","author":"Candes","year":"2006","journal-title":"IEEE Trans. Inf. Theory"},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"1289","DOI":"10.1109\/TIT.2006.871582","article-title":"Compressed sensing","volume":"52","author":"Donoho","year":"2006","journal-title":"IEEE Trans. Inf. Theory"},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"1182","DOI":"10.1002\/mrm.21391","article-title":"Sparse MRI: The application of compressed sensing for rapid MR imaging","volume":"58","author":"Lustig","year":"2007","journal-title":"Magn. Reson. Med."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1259\/bjr.20150487","article-title":"Compressed sensing MRI: A review of the clinical literature","volume":"88","author":"Jaspan","year":"2015","journal-title":"Br. J. Radiol."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"1217","DOI":"10.1088\/0266-5611\/10\/6\/003","article-title":"Analysis of bounded variation penalty methods for ill-posed problems","volume":"10","author":"Acar","year":"1994","journal-title":"Inverse Probl."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"466","DOI":"10.1002\/jmri.21585","article-title":"Acquisition and reconstruction of undersampled radial data for myocardial perfusion magnetic resonance imaging","volume":"29","author":"Adluru","year":"2009","journal-title":"J. Magn. Reson. Imaging"},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"63","DOI":"10.1016\/0041-5553(68)90034-7","article-title":"The error principle in the solution of operational equations by the regularization method","volume":"8","author":"Morozov","year":"1968","journal-title":"USSR Comput. Math. Math. Phys."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"1487","DOI":"10.1137\/0914086","article-title":"The use of the L-curve in the regularization of discrete ill-posed problems","volume":"14","author":"Hansen","year":"1993","journal-title":"SIAM J. Sci. Comput."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"1540","DOI":"10.1109\/TIP.2008.2001404","article-title":"Monte-Carlo Sure: A Black-Box Optimization of Regularization Parameters for General Denoising Algorithms","volume":"17","author":"Ramani","year":"2008","journal-title":"IEEE Trans. Image Process."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"1411","DOI":"10.1109\/TMI.2013.2257829","article-title":"Non-cartesian MRI reconstruction with automatic regularization via Monte-Carlo SURE","volume":"32","author":"Ramani","year":"2013","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"1760","DOI":"10.1002\/mrm.24840","article-title":"Monte-Carlo SURE-based parameter selection for parallel magnetic resonance imaging reconstruction","volume":"71","author":"Weller","year":"2014","journal-title":"Magn. Reson. Med."},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"025005","DOI":"10.1088\/0266-5611\/28\/2\/025005","article-title":"Sparsity-promoting Bayesian inversion","volume":"28","author":"Kolehmainen","year":"2012","journal-title":"Inverse Probl."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"B644","DOI":"10.1137\/120876277","article-title":"Sparse tomography","volume":"35","author":"Kallonen","year":"2013","journal-title":"SIAM J. Comput."},{"key":"ref_15","unstructured":"Niinim\u00e4ki, K. (2013). Computational Optimization Methods for Large-Scale Inverse Problems. [Ph.D. Thesis, University of Eastern Finland]."},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"938","DOI":"10.1137\/15M1034076","article-title":"Multiresolution parameter choice method for total variation regularised tomography","volume":"9","author":"Lassas","year":"2016","journal-title":"SIAM J. Imaging Sci."},{"key":"ref_17","doi-asserted-by":"crossref","unstructured":"Niinim\u00e4ki, K., Hanhela, M., and Kolehmainen, V. (2019). Parameter selection in dynamic contrast-enhanced magnetic resonance tomography. Centre International de Rencontres Math\u00e9matiques, Springer.","DOI":"10.1007\/978-3-030-48634-1_6"},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"67","DOI":"10.1023\/B:BREA.0000010700.11092.f4","article-title":"Monitoring response to primary chemotherapy in breast cancer using dynamic contrast-enhanced magnetic resonance imaging","volume":"83","author":"Martincich","year":"2004","journal-title":"Breast Cancer Res. Treat."},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1007\/s10549-004-5819-2","article-title":"Role of dynamic contrast enhanced MRI in monitoring early response of locally advanced breast cancer to neoadjuvant chemotherapy","volume":"91","author":"Pickles","year":"2005","journal-title":"Breast Cancer Res. Treat."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"1269","DOI":"10.1007\/s00234-015-1582-9","article-title":"Early biomarkers from dynamic contrast-enhanced magnetic resonance imaging to predict the response to antiangiogenic therapy in high-grade gliomas","volume":"57","author":"Piludu","year":"2015","journal-title":"Neuroradiology"},{"key":"ref_21","doi-asserted-by":"crossref","unstructured":"Merali, Z., Huang, K., Mikulis, D., Silver, F., and Kassner, A. (2017). Evolution of blood-brain-barrier permeability after acute ischemic stroke. PLoS ONE, 12.","DOI":"10.1371\/journal.pone.0171558"},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"433","DOI":"10.1212\/WNL.0000000000003566","article-title":"DCE-MRI blood\u2013brain barrier assessment in acute ischemic stroke","volume":"88","author":"Villringer","year":"2017","journal-title":"Neurology"},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"560","DOI":"10.1109\/TSP.2002.807005","article-title":"Nonuniform fast Fourier transforms using min-max interpolation","volume":"51","author":"Fessler","year":"2003","journal-title":"IEEE Trans. Signal Process."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"1135","DOI":"10.1214\/aos\/1176345632","article-title":"Estimation of the Mean of a Multivariate Normal Distribution","volume":"9","author":"Stein","year":"1981","journal-title":"Ann. Stat."},{"key":"ref_25","first-page":"26","article-title":"An in vivo MRI template set for morphometry, tissue segmentation, and fMRI localization in rats","volume":"5","author":"Sumiyoshi","year":"2011","journal-title":"Front. Neuroinform."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"1420","DOI":"10.1016\/j.mri.2010.06.017","article-title":"Reproducibility of the gadolinium concentration measurements and of the fitting parameters of the vascular input function in the superior sagittal sinus in a patient population","volume":"28","author":"Lavini","year":"2010","journal-title":"Magn. Reson. Imaging"},{"key":"ref_27","first-page":"223","article-title":"Estimating kinetic parameters from dynamic contrast-enhanced T1-weighted MRI of a diffusable tracer: Standardized quantities and symbols","volume":"10","author":"Tofts","year":"1999","journal-title":"J. Magn. Reson."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"1334","DOI":"10.1007\/s10851-020-00985-2","article-title":"Temporal regularization methods in DCE-MRI","volume":"62","author":"Hanhela","year":"2020","journal-title":"J. Math Imaging Vis."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"074001","DOI":"10.1088\/1361-6420\/aac3af","article-title":"Dynamic MRI reconstruction from undersampled data with an anatomical prescan","volume":"34","author":"Rasch","year":"2018","journal-title":"Inverse Probl."},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"68","DOI":"10.1109\/TMI.2006.885337","article-title":"An optimal radial profile order based on the golden ratio for time-resolved MRI","volume":"26","author":"Winkelmann","year":"2007","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"145","DOI":"10.1016\/j.acha.2005.11.003","article-title":"Fast and accurate Polar Fourier transform","volume":"21","author":"Averbuch","year":"2006","journal-title":"Appl. Comput. Harmon. Anal."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"816","DOI":"10.1109\/TBME.2016.2578930","article-title":"Pseudo-polar Fourier transform based compressed sensing MRI","volume":"64","author":"Yang","year":"2016","journal-title":"IEEE Trans. Biomed. Eng."},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"e4250","DOI":"10.1002\/nbm.4250","article-title":"Detection of lentiviral suicide gene therapy in C6 rat glioma using hyperpolarised [1-13C]pyruvate","volume":"33","author":"Olsson","year":"2020","journal-title":"NMR Biomed."},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"120","DOI":"10.1007\/s10851-010-0251-1","article-title":"A first-order primal-dual algorithm for convex problems with applications to imaging","volume":"40","author":"Chambolle","year":"2011","journal-title":"J. Math. Imaging Vis."},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"3065","DOI":"10.1088\/0031-9155\/57\/10\/3065","article-title":"Convex optimization problem prototyping for image reconstruction in computed tomography with the Chambolle-Pock algorithm","volume":"57","author":"Sidky","year":"2012","journal-title":"Phys. Med. Biol."},{"key":"ref_36","doi-asserted-by":"crossref","unstructured":"Pock, T., and Chambolle, A. (2011, January 6\u201313). Diagonal preconditioning for first order primal-dual algorithms in convex optimization. Proceedings of the 2011 International Conference on Computer Vision, Barcelona, Spain.","DOI":"10.1109\/ICCV.2011.6126441"},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"55012","DOI":"10.1088\/0266-5611\/30\/5\/055012","article-title":"A primal\u2013dual hybrid gradient method for nonlinear operators with applications to MRI","volume":"30","author":"Valkonen","year":"2014","journal-title":"Inverse Probl."},{"key":"ref_38","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1109\/TMI.2013.2268898","article-title":"Retrospective Rigid Motion Correction in k-Space for Segmented Radial MRI","volume":"33","author":"Vaillant","year":"2014","journal-title":"IEEE Trans. Med. Imaging"}],"container-title":["Journal of Imaging"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2313-433X\/7\/2\/38\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T05:26:49Z","timestamp":1760160409000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2313-433X\/7\/2\/38"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,2,20]]},"references-count":38,"journal-issue":{"issue":"2","published-online":{"date-parts":[[2021,2]]}},"alternative-id":["jimaging7020038"],"URL":"https:\/\/doi.org\/10.3390\/jimaging7020038","relation":{},"ISSN":["2313-433X"],"issn-type":[{"value":"2313-433X","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,2,20]]}}}