Awesome
Score-Based Priors for Bayesian Inverse Imaging
This project turns score-based diffusion models into explicit priors for Bayesian inverse problems in imaging. A "score-based prior" allows us to model complex, data-driven posterior distributions while using established sampling and optimization methods (e.g., MCMC, variational inference) that require an explicit log-density function. In our work, we propose a variational-inference approach to posterior sampling, in addition to empirically validating the accuracy and variance of a score-based prior.
UPDATE (2/5/24): The repo has been updated to support the surrogate score-based prior based on an evidence lower bound (ELBO). The ELBO-based prior is much more time- and memory-efficient than the ODE-based exact score-based prior.
@inproceedings{feng2023score,
title={Score-Based Diffusion Models as Principled Priors for Inverse Imaging},
author={Feng, Berthy T and Smith, Jamie and Rubinstein, Michael and Chang, Huiwen and Bouman, Katherine L and Freeman, William T},
booktitle={International Conference on Computer Vision (ICCV)},
year={2023},
organization={IEEE}
}
@article{feng2024variational,
title={Variational Bayesian Imaging with an Efficient Surrogate Score-based Prior},
author={Berthy Feng and Katherine Bouman},
journal={Transactions on Machine Learning Research},
issn={2835-8856},
year={2024},
url={https://openreview.net/forum?id=db2pFKVcm1}
}
Getting started
Clone the repo:
git clone https://github.com/berthyf96/score_prior
Create Conda environment:
conda update -n base -c defaults conda
conda create -n score_prior -y python=3.10
conda activate score_prior
Install dependencies (please open an issue if this does not work out-of-the-box):
sh conda.sh
The demos are a great way to start!
demos/2d_posterior_sampling.ipynb
introduces our approach through a toy inverse problem with a known analytical prior. This is the best way to gain intuition for DPI + score-based priors.demos/interferometry.ipynb
applies DPI + score-based priors to the nonconvex problem of imaging from interferometric measurements. It uses a score-based prior trained on black-hole simulations to image a synthetic black hole. This demo focuses on the surrogate ELBO-based prior.demos/denoising.ipynb
demonstrates DPI + score-based priors on a simple image denoising task. It focuses on the exact ODE-based prior. As this is computationally heavy, it is recommended to start with the above two demos.
To run demos/denoising.ipynb
and demos/interferometry.ipynb
with pre-trained score models, please download the score_checkpoints
folder from Box.
Example workflow
Train score model
To train a score-based prior, you first need to train a score-based diffusion model on the desired dataset:
python train.py \
--config configs/score_config.py \
--workdir score_checkpoints/CELEBA_32 \
--config.data.dataset CELEBA \
--config.data.image_size 32 \
--config.data.tfds_dir /tmp/tensorflow_datasets \
--config.training.n_iters 1000000
Once trained, the score-based diffusion model with parameters $\theta$ represents the image prior $p_\theta$.
Optimize DPI for posterior sampling
Often, our goal is to sample from a posterior $$p_\theta(\mathbf{x}\mid\mathbf{y})\propto p(\mathbf{y}\mid\mathbf{x})p_\theta(\mathbf{x}).$$ We use Deep Probabilistic Imaging (DPI) for posterior sampling, which optimizes the parameters of a RealNVP to approximate the target posterior. Here is an example command to perform DPI optimization:
python train_dpi.py \
--score_model_config configs/score_config.py \
--config configs/dpi_config.py \
--workdir dpi_checkpoints/Denoising_CELEBAxCELEBA_32 \
--config.prob_flow.score_model_dir score_checkpoints/CELEBA_32/checkpoints/checkpoint_1000000 \
--config.data.dataset CELEBA \
--config.data.image_size 32 \
--config.data.num_channels 3 \
--config.likelihood.likelihood Denoising \
--config.likelihood.noise_scale 0.2 \
--config.training.batch_size 32 \
--config.training.n_iters 10000 \
--config.training.snapshot_freq 100 \
--config.optim.prior ode \
--config.optim.learning_rate 2e-4 \
--config.optim.grad_clip 1. \
--score_model_config.training.sde vpsde
Replace config.prob_flow.score_model_dir
with the path to your trained score model.
You can also change the forward model by changing config.likelihood
. To use the ELBO-based prior, set config.optim.prior
to 'dsm'
.
Customization
Datasets can be added in score_flow/datasets.py
.
Forward models can be added in forward_models.py
.