{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T01:48:56Z","timestamp":1760147336750,"version":"build-2065373602"},"reference-count":24,"publisher":"MDPI AG","issue":"3","license":[{"start":{"date-parts":[[2023,1,26]],"date-time":"2023-01-26T00:00:00Z","timestamp":1674691200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000002","name":"NIH\/NCI","doi-asserted-by":"publisher","award":["#CA206171"],"award-info":[{"award-number":["#CA206171"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Sensors"],"abstract":"<jats:p>Most penalized maximum likelihood methods for tomographic image reconstruction based on Bayes\u2019 law include a freely adjustable hyperparameter to balance the data fidelity term and the prior\/penalty term for a specific noise\u2013resolution tradeoff. The hyperparameter is determined empirically via a trial-and-error fashion in many applications, which then selects the optimal result from multiple iterative reconstructions. These penalized methods are not only time-consuming by their iterative nature, but also require manual adjustment. This study aims to investigate a theory-based strategy for Bayesian image reconstruction without a freely adjustable hyperparameter, to substantially save time and computational resources. The Bayesian image reconstruction problem is formulated by two probability density functions (PDFs), one for the data fidelity term and the other for the prior term. When formulating these PDFs, we introduce two parameters. While these two parameters ensure the PDFs completely describe the data and prior terms, they cannot be determined by the acquired data; thus, they are called complete but unobservable parameters. Estimating these two parameters becomes possible under the conditional expectation and maximization for the image reconstruction, given the acquired data and the PDFs. This leads to an iterative algorithm, which jointly estimates the two parameters and computes the to-be reconstructed image by maximizing a posteriori probability, denoted as joint-parameter-Bayes. In addition to the theoretical formulation, comprehensive simulation experiments are performed to analyze the stopping criterion of the iterative joint-parameter-Bayes method. Finally, given the data, an optimal reconstruction is obtained without any freely adjustable hyperparameter by satisfying the PDF condition for both the data likelihood and the prior probability, and by satisfying the stopping criterion. Moreover, the stability of joint-parameter-Bayes is investigated through factors such as initialization, the PDF specification, and renormalization in an iterative manner. Both phantom simulation and clinical patient data results show that joint-parameter-Bayes can provide comparable reconstructed image quality compared to the conventional methods, but with much less reconstruction time. To see the response of the algorithm to different types of noise, three common noise models are introduced to the simulation data, including white Gaussian noise to post-log sinogram data, Poisson-like signal-dependent noise to post-log sinogram data and Poisson noise to the pre-log transmission data. The experimental outcomes of the white Gaussian noise reveal that the two parameters estimated by the joint-parameter-Bayes method agree well with simulations. It is observed that the parameter introduced to satisfy the prior\u2019s PDF is more sensitive to stopping the iteration process for all three noise models. A stability investigation showed that the initial image by filtered back projection is very robust. Clinical patient data demonstrated the effectiveness of the proposed joint-parameter-Bayes and stopping criterion.<\/jats:p>","DOI":"10.3390\/s23031374","type":"journal-article","created":{"date-parts":[[2023,1,27]],"date-time":"2023-01-27T01:27:58Z","timestamp":1674782878000},"page":"1374","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":4,"title":["A Joint-Parameter Estimation and Bayesian Reconstruction Approach to Low-Dose CT"],"prefix":"10.3390","volume":"23","author":[{"ORCID":"https:\/\/orcid.org\/0000-0001-6169-3478","authenticated-orcid":false,"given":"Yongfeng","family":"Gao","sequence":"first","affiliation":[{"name":"Department of Radiology, Stony Brook University, Stony Brook, NY 11794, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Siming","family":"Lu","sequence":"additional","affiliation":[{"name":"Department of Radiology, Stony Brook University, Stony Brook, NY 11794, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Yongyi","family":"Shi","sequence":"additional","affiliation":[{"name":"Department of Radiology, Stony Brook University, Stony Brook, NY 11794, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Shaojie","family":"Chang","sequence":"additional","affiliation":[{"name":"Department of Radiology, Stony Brook University, Stony Brook, NY 11794, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Hao","family":"Zhang","sequence":"additional","affiliation":[{"name":"Department of Medical Physics, Memorial Sloan Kettering Cancer Center, New York, NY 10065, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Wei","family":"Hou","sequence":"additional","affiliation":[{"name":"Department of Preventive Medicine, Stony Brook University, Stony Brook, NY 11794, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Lihong","family":"Li","sequence":"additional","affiliation":[{"name":"Department of Engineering Science and Physics, CUNY\/CSI, Staten Island, NY 10314, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Zhengrong","family":"Liang","sequence":"additional","affiliation":[{"name":"Department of Radiology, Stony Brook University, Stony Brook, NY 11794, USA"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"1968","published-online":{"date-parts":[[2023,1,26]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"433","DOI":"10.1016\/j.jhep.2009.05.023","article-title":"Non-invasive assessment and quantification of liver steatosis by ultrasound, computed tomography and magnetic resonance","volume":"51","author":"Schwenzer","year":"2009","journal-title":"J. Hepatol."},{"key":"ref_2","unstructured":"Deans, S.R. (2007). The Radon Transform and Some of Its Applications, Courier Corporation."},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"21","DOI":"10.1109\/TNS.1974.6499235","article-title":"The Fourier reconstruction of a head section","volume":"21","author":"Shepp","year":"1974","journal-title":"IEEE Trans. Nucl. Sci."},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"2409","DOI":"10.1109\/TMI.2017.2768978","article-title":"Low-Dose CT: What Has Been Done, and What Challenges Remain","volume":"36","author":"Liang","year":"2017","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"1714","DOI":"10.1093\/imanum\/drr051","article-title":"Iterative parameter choice by discrepancy principle","volume":"32","author":"Jin","year":"2012","journal-title":"IMA J. Numer. Anal."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"215","DOI":"10.1080\/00401706.1979.10489751","article-title":"Generalized cross-validation as a method for choosing a good ridge parameter","volume":"21","author":"Golub","year":"1979","journal-title":"Technometrics"},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"561","DOI":"10.1137\/1034115","article-title":"Analysis of discrete ill-posed problems by means of the L-curve","volume":"34","author":"Hansen","year":"1992","journal-title":"SIAM Rev."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"1415","DOI":"10.1137\/100790756","article-title":"A regularization parameter for nonsmooth Tikhonov regularization","volume":"33","author":"Kazufumi","year":"2011","journal-title":"SIAM J. Sci. Comput."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"2466","DOI":"10.1109\/TMI.2017.2759819","article-title":"Z-index parameterization for volumetric CT image reconstruction via 3-D dictionary learning","volume":"36","author":"Bai","year":"2017","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_10","first-page":"383","article-title":"Optimal Selection for Regularization Parameter in Iterative CT Reconstruction Based on the Property of Natural Image Statistics","volume":"13","author":"Mou","year":"2015","journal-title":"Fully 3D"},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"413","DOI":"10.1109\/34.134041","article-title":"Image restoration using Gibbs priors: Boundary modeling, treatment of blurring, and selection of hyperparameter","volume":"5","author":"Johnson","year":"1991","journal-title":"IEEE Trans. Pattern Anal. Mach. Intell."},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"7","DOI":"10.1117\/1.1526103","article-title":"Bayesian image reconstruction for transmission tomography using deterministic annealing","volume":"12","author":"Hsiao","year":"2003","journal-title":"J. Electron. Imaging"},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"111916","DOI":"10.1118\/1.4820478","article-title":"A full-spectral Bayesian reconstruction approach based on the material decomposition model applied in dual-energy computed tomography","volume":"40","author":"Cai","year":"2013","journal-title":"Med. Phys."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"225","DOI":"10.1109\/TMI.2016.2603843","article-title":"Cone beam X-ray luminescence computed tomography based on Bayesian method","volume":"36","author":"Zhang","year":"2016","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"313","DOI":"10.1109\/TMI.1987.4307849","article-title":"Stopping rule for the MLE algorithm based on statistical hypothesis testing","volume":"6","author":"Veklerov","year":"1987","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"1221","DOI":"10.1088\/0031-9155\/35\/9\/003","article-title":"Statistical stopping criteria for iterative maximum likelihood reconstruction of emission images","volume":"35","author":"Herbert","year":"1990","journal-title":"Phys. Med. Biol."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"131","DOI":"10.1016\/j.compmedimag.2009.07.006","article-title":"PET image reconstruction: A stopping rule for the MLEM algorithm based on properties of the updating coefficients","volume":"34","author":"Gaitanis","year":"2010","journal-title":"Comput. Med. Imaging Graph."},{"key":"ref_18","doi-asserted-by":"crossref","unstructured":"Gao, Y., Liang, Z., Lu, S., Shi, Y., Chang, S., and Hou, W. (November, January 31). Hyperparameter Selection for Bayesian Image Reconstruction by Mimicking Physical Crystallization. Proceedings of the 2020 IEEE Symposium on Nuclear Science (NSS\/MIC), Boston, MA, USA.","DOI":"10.1109\/NSS\/MIC42677.2020.9508084"},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"5032","DOI":"10.1002\/mp.14449","article-title":"Characterization of tissue-specific pre-log Bayesian CT reconstruction by texture\u2013dose relationship","volume":"47","author":"Gao","year":"2020","journal-title":"Med. Phys."},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"4051","DOI":"10.1118\/1.4722751","article-title":"Variance Analysis of X-ray CT Sinograms in the Presence of Electronic Noise Background","volume":"39","author":"Ma","year":"2012","journal-title":"Med. Phys."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"1272","DOI":"10.1109\/TMI.2006.882141","article-title":"Penalized Weighted Least-squares Approach to Sinogram Noise Reduction and Image Reconstruction for Low-dose X-ray CT","volume":"25","author":"Wang","year":"2006","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"041916","DOI":"10.1118\/1.4869160","article-title":"Deriving Adaptive MRF Coefficients from Previous Normal-dose CT Scan for Low-dose Image Reconstruction via Penalized Weighted Least-squares Minimization","volume":"41","author":"Zhang","year":"2014","journal-title":"Med. Phys."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"860","DOI":"10.1109\/TMI.2015.2498148","article-title":"Extracting Information from Previous Full-dose CT Scan for Knowledge-based Bayesian Reconstruction of Current Low-dose CT Images","volume":"35","author":"Zhang","year":"2016","journal-title":"IEEE Trans. Med. Imaging"},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"1981","DOI":"10.1109\/TMI.2018.2890788","article-title":"A Feasibility Study of Extracting Tissue Textures from a Previous Full-Dose CT Database as Prior Knowledge for Bayesian Reconstruction of Current Low-dose CT Images","volume":"38","author":"Gao","year":"2019","journal-title":"IEEE Trans. Med. Imaging"}],"container-title":["Sensors"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/1424-8220\/23\/3\/1374\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,10]],"date-time":"2025-10-10T18:16:04Z","timestamp":1760120164000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/1424-8220\/23\/3\/1374"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,1,26]]},"references-count":24,"journal-issue":{"issue":"3","published-online":{"date-parts":[[2023,2]]}},"alternative-id":["s23031374"],"URL":"https:\/\/doi.org\/10.3390\/s23031374","relation":{},"ISSN":["1424-8220"],"issn-type":[{"type":"electronic","value":"1424-8220"}],"subject":[],"published":{"date-parts":[[2023,1,26]]}}}