Papers › Variational Bayesian Imaging with an Efficient Surrogate Score-based Prior

Variational Bayesian Imaging with an Efficient Surrogate Score-based Prior

5 Sep 2023arXiv:2309.01949archive 2025-07-28

Berthy T. Feng, Katherine L. Bouman

We propose a surrogate function for efficient yet principled use of score-based priors in Bayesian imaging. We consider ill-posed inverse imaging problems in which one aims for a clean image posterior given incomplete or noisy measurements. Since the measurements do not uniquely determine a true image, a prior is needed to constrain the solution space. Recent work turned score-based diffusion models into principled priors for solving ill-posed imaging problems by appealing to an ODE-based log-probability function. However, evaluating the ODE is computationally inefficient and inhibits posterior estimation of high-dimensional images. Our proposed surrogate prior is based on the evidence lower bound of a score-based diffusion model. We demonstrate the surrogate prior on variational inference for efficient approximate posterior sampling of large images. Compared to the exact prior in previous work, our surrogate accelerates optimization of the variational image distribution by at least two orders of magnitude. We also find that our principled approach gives more accurate posterior estimation than non-variational diffusion-based approaches that involve hyperparameter-tuning at inference. Our work establishes a practical path forward for using score-based diffusion models as general-purpose image priors.

PaperPDFCodeCode Syntology ran

In Syntology Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="2309.01949")

Code

Syntology Ran 9 of 13 code samples harvested from 1 repository linked to this paper; 4 have no recorded run. Of those that ran: 1 ran · violated contract; 4 ran · our draft was wrong; 4 ran with no contract checked.

By repository: official repository: 13 samples from 1 repository, 9 ran. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

berthyf96/score_prior officialmentioned on GitHubjax report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

13 samples harvested; 9 ran; 0 honoured the contract we drafted; 4 have no recorded run. Read from Syntology's graph 2026-09-24; that is when this build read the record, not when the samples ran.

1ran · violated contract
4ran · our draft was wrong
4ran
4unverified

Licence: 13 of the 13 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from berthyf96/score_prior. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

batch_add berthyf96/score_prior/score_flow/utils.py official repository ran fingerprinted no licence file found · pointer only · ba575cda3781df1f · report
batch_mul berthyf96/score_prior/score_flow/sde_lib.py official repository ran · violated contract fingerprinted no licence file found · pointer only · 705aea0955c9050d · report
get_data_inverse_scaler berthyf96/score_prior/score_flow/datasets.py official repository ran · our draft was wrong no licence file found · pointer only · 6c419026e778dee2 · report
get_data_scaler berthyf96/score_prior/score_flow/datasets.py official repository ran · our draft was wrong no licence file found · pointer only · 346d2e6cc8b48a5d · report
get_optimizer berthyf96/score_prior/score_flow/losses.py official repository ran no licence file found · pointer only · 69394cd870ed9898 · report
optimization_manager berthyf96/score_prior/score_flow/losses.py official repository ran no licence file found · pointer only · 53f39e2a10eaea69 · report
plot_image_grid berthyf96/score_prior/score_flow/utils.py official repository ran no licence file found · pointer only · e2e76a11168d6b4a · report
register_corrector berthyf96/score_prior/score_flow/sampling.py official repository ran · our draft was wrong no licence file found · pointer only · a05b1d0cf3f508b4 · report
register_predictor berthyf96/score_prior/score_flow/sampling.py official repository ran · our draft was wrong no licence file found · pointer only · 473066ddf62a075a · report
classifier_fn_from_tfhub berthyf96/score_prior/score_flow/evaluation.py official repository unverified no licence file found · pointer only · 54a00db9aadf092f · report
crop_resize berthyf96/score_prior/score_flow/datasets.py official repository unverified no licence file found · pointer only · fd14d238ccd5f679 · report
get_predictor berthyf96/score_prior/score_flow/sampling.py official repository unverified no licence file found · pointer only · a58facb8fe71823a · report
load_dataset_stats berthyf96/score_prior/score_flow/evaluation.py official repository unverified no licence file found · pointer only · beaa2ce502047cbf · report

Tasks

Variational Inference

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Methods

DiffusionVariational Inference

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections