A light fork of sigpy, a signal processing library for Python.
Project description
NOTE: This is a fork of the repository below. Changes have been made to lighten it’s dependency on Numba and PyWavelets. If you are using a function that used Numba either install with the [full] extra or expect much slower results.
SigPy
Source Code | Documentation | MRI Recon Tutorial | MRI Pulse Design Tutorial
SigPy is a package for signal processing, with emphasis on iterative methods. It is built to operate directly on NumPy arrays on CPU and CuPy arrays on GPU. SigPy also provides several domain-specific submodules: sigpy.plot for multi-dimensional array plotting, sigpy.mri for MRI reconstruction, and sigpy.mri.rf for MRI pulse design.
Installation
SigPy requires Python version >= 3.5. The core module depends on numba, numpy, PyWavelets, scipy, and tqdm.
Additional features can be unlocked by installing the appropriate packages. To enable the plotting functions, you will need to install matplotlib. To enable CUDA support, you will need to install cupy. And to enable MPI support, you will need to install mpi4py.
Via conda
We recommend installing SigPy through conda:
conda install -c frankong sigpy # (optional for plot support) conda install matplotlib # (optional for CUDA support) conda install cupy # (optional for MPI support) conda install mpi4py
Via pip
SigPy can also be installed through pip:
pip install sigpy # (optional for plot support) pip install matplotlib # (optional for CUDA support) pip install cupy # (optional for MPI support) pip install mpi4py
Installation for Developers
If you want to contribute to the SigPy source code, we recommend you install it with pip in editable mode:
cd /path/to/sigpy pip install -e .
To run tests and contribute, we recommend installing the following packages:
pip install coverage ruff sphinx sphinx_rtd_theme black isort
and run the script run_tests.sh.
Features
CPU/GPU Signal Processing Functions
SigPy provides signal processing functions with a unified CPU/GPU interface. For example, the same code can perform a CPU or GPU convolution on the input array device:
# CPU convolve
x = numpy.array([1, 2, 3, 4, 5])
y = numpy.array([1, 1, 1])
z = sigpy.convolve(x, y)
# GPU convolve
x = cupy.array([1, 2, 3, 4, 5])
y = cupy.array([1, 1, 1])
z = sigpy.convolve(x, y)
Iterative Algorithms
SigPy also provides convenient abstractions and classes for iterative algorithms. A compressed sensing experiment can be implemented in four lines using SigPy:
# Given some observation vector y, and measurement matrix mat
A = sigpy.linop.MatMul([n, 1], mat) # define forward linear operator
proxg = sigpy.prox.L1Reg([n, 1], lamda=0.001) # define proximal operator
x_hat = sigpy.app.LinearLeastSquares(A, y, proxg=proxg).run() # run iterative algorithm
PyTorch Interoperability
Want to do machine learning without giving up signal processing? SigPy has convenient functions to convert arrays and linear operators into PyTorch Tensors and Functions. For example, given a cupy array x, and a Linop A, we can convert them to Pytorch:
x_torch = sigpy.to_pytorch(x)
A_torch = sigpy.to_pytorch_function(A)
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for sigpy_lite-0.1.26-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 62fa81d4c86b7470f3a2c8f065b76dd9b5b40318e18453c7d6fac05870534cdf |
|
MD5 | c6fc01756b65dc036ad8821a695bde33 |
|
BLAKE2b-256 | 6fc8466c32b34c56002b8d597465f3f33b6c66680a99ee1068d88a6ec9746aa8 |