Skip to main content

CUQIpy plugin for CIL

Project description

CUQIpy-CIL

CUQIpy-CIL is a plugin for the CUQIpy software package.

It adds a thin wrapper around Computed Tomography (CT) forward models from the Core Imaging Library (CIL).

Installation

First install CIL. Then install CUQIpy-CIL with pip:

pip install cuqipy-cil

If CUQIpy is not installed, it will be installed automatically.

Quickstart

import numpy as np
import matplotlib.pyplot as plt
import cuqi
import cuqipy_cil

# Load a CT forward model and data from testproblem library
A, y_data, info = cuqipy_cil.testproblem.ParallelBeam2D.get_components(
    im_size=(128, 128),
    det_count=128,
    angles=np.linspace(0, np.pi, 180),
    phantom="shepp-logan"
)

# Set up Bayesian model
x = cuqi.distribution.Gaussian(np.zeros(A.domain_dim), cov=1) # x ~ N(0, 1)
y = cuqi.distribution.Gaussian(A@x, cov=0.05**2)              # y ~ N(Ax, 0.05^2)

# Set up Bayesian Problem
BP = cuqi.problem.BayesianProblem(y, x).set_data(y=y_data)

# Sample from the posterior
samples = BP.sample_posterior(200)

# Analyze the samples
info.exactSolution.plot(); plt.title("Exact solution")
y_data.plot(); plt.title("Data")
samples.plot_mean(); plt.title("Posterior mean")
samples.plot_std(); plt.title("Posterior standard deviation")

For more examples, see the demos folder.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

CUQIpy-CIL-0.5.0.post0.dev2.tar.gz (14.3 kB view hashes)

Uploaded Source

Built Distribution

CUQIpy_CIL-0.5.0.post0.dev2-py3-none-any.whl (15.8 kB view hashes)

Uploaded Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page