Fast Richardson-Lucy deconvolution of 3D volume data using GPU or CPU with napari plugin.
Richardson-Lucy deconvolution for fishes, scientists and engineers.
This software is for filtering 3D data using the Richardson-Lucy deconvolution algorithm.
Richardson-Lucy is an iterative deconvolution algorithm that is used to remove point spread function (PSF) or optical transfer function (OTF) artifacts from experimental images.
The method was originally developed for astronomy to remove optical effects and simultaneously reduce poisson noise in 2D images.
Lucy, L. B. An iterative technique for the rectification of observed distributions. The Astronomical Journal 79, 745 (1974). DOI: 10.1086/111605
The method can also be applied to 3D data. Nowadays this filtering technique is also widely used by microscopists.
The Richardson-Lucy deconvolution algorigthm is iterative. Each iteration involves the calculation of 2 convolutions, one element-wise multiplication and one element-wise division.
When dealing with 3D data, the Richardson-Lucy algorithm is quite computional intensive primarly due to the calculation of the convolution, and can take a while to complete depending on the resources available. Convolution is significantly sped up using FFT compared to raw convolution.
This software was developed with the aim to make the R-L computation faster by exploiting GPU resources, and with the use of FFT convolution.
To make RedLionfish easily accessible, it is available through PyPi and anaconda (conda-forge channel). A useful plugin for Napari is also available.
Please note that this software only works with 3D data. For 2D data there are many alternatives such as the DeconvolutionLab2 in Fiji (ImageJ) and sckikit-image.
You can now use the Napari's plugin installation in Menu -> Plugins -> Install/Uninstall Plugins.... However, if you chose to use this method, GPU acceleration may not be available and it will use the CPU backend. Better check.
Alternatively, if you follow the installation instructions below, and install the napari in the same python environment then the plugin should be immediately available in the Menu -> Plugins -> RedLionfish.
Previously there was a problem in installing using
pip, because no PyOpenCL wheels for windows were avaiable. It is now avaialble. This package can be installed using pip or conda.
This package is available in conda-forge channel. It contains the precompiled libraries and it will install all the requirments for GPU-accelerated RL calculations.
conda install redlionfish -c conda-forge
Install from PyPi¶
pip install redlionfish
In Linux , the package
ocl-icd-system may also be useful.
conda install reikna pyopencl ocl-icd-system -c conda-forge
Manual installation using the conda package file.¶
Download the appropriate conda package .bz2 at https://github.com/rosalindfranklininstitute/RedLionfish/releases
In the command line, successively run:
conda install <filename.bz2> conda update --all -c conda-forge
The second line is needed because you are installing from a local file, conda installer will not install dependencies. Right after this you should run the update command given.
Manual installation (advanced and for developers)¶
Please note that in order to use OpenCL GPU accelerations, PyOpenCL must be installed. The best way to get it working is to install it under a conda environment.
The installation is similar to the previously described for PyPi.
conda install reikna pyopencl
conda install reikna pyopencl ocl-icd-system -c conda-forge (Linux)
Clone/download from source https://github.com/rosalindfranklininstitute/RedLionfish/
python setup.py install
If you want to test and modify the code then you should probably install in debug mode using:
python setup.py develop
pip install -e .
The software has algorithms for Richardson-Lucy deconvolution that use either CPU and GPU.
The CPU version is very similar to the skimage.restoration.richardson_lucy code, with improvments in speed. major differences are:
- the convolution steps use FFT only.
- PSF and PSF-flipped FFTs are precalculated before starting iterations.
The GPU version, was written in to use Reikna package, which does FFT using OpenCL, via PyOpenCL.
Unfortunately, a major limitation in RAM usage exists with PyOpenCL. Large 3D data volumes with cause out-of-memory error when trying to upload data to the GPU for FFT calculations. As such, to overcome this problem, a block algorithm is used, which splits data into blocks with padded data. The results are then combined together to give the final result. This affects the perfomance of the calculation rather significantly, but with the advantage of being possible to handle large data volumes.
If Richardson-Lucy deconvolution using the GPU method fails, RedLionfish will falls back to CPU calculation. Check console output for messages.
If you are using the RedLionfish in your code, note that, by default,
def doRLDeconvolutionFromNpArrays() method it uses the GPU OpenCL version.
Many examples can be found in `/scripts' folder.
A quick and benchmarking installation can be run from the proect root using the command:
or in windows
This will print out information about your GPU device (if available) and run some deconvolutions. It initially creates some data programatically, convolutes with a gaussian PSF, and add Poisson noise. Then it executes executes the Richardson-Lucy deconvolution calculation using CPU and GPU methods, for 10 iterations. During the calculation it will print some information to the console/terminal, including the time it takes to run the calculation.
Computer generated data and an experimental PSF can be found in
Testing Redlionfish in napari¶
Here is an example testing the Redlionfish plugin in napari:
- load data
scripts/testdata/gendata_psfconv_poiss_large.tif(can use draga and drop)
- load psf data
- In the RedLionfish side window ensure that 'gendata_psfconv_poiss_large' is selected in data dropdown widget, and
PSF_RFI_8bitis selected in psfdata widget.
- Choose number of iterations (default=10)
- Click 'Go' button and wait until result shows as a new data layer.
- Use controls of the left panel to compare before and after RL deconvolution: select 'RL-deconvolution' layer and set colormap to red. Hide PSF_RFI_8bit. Make sure that both 'RL-deconvolution' and 'gendata-psfconv' are visible. Now, hide/unhide RL-deconvolution layer to see before and after deconvolution. Adjust contrast limits of each layer as desired.
GPU vs CPU¶
You may notice that choosing GPU does not make RL-calculation much faster compared with CPU, and sometimes is slower.
Which method runs the R-L deconvolution faster. That depends on the computer configuration/architecture.
GPU calculations will be generally faster than CPU with bigger data volumes.
GPU calculation will be significantly faster if using a dedicated GPU card.
Please see benchmark values that highlights significant variability in calculation speeds.
Please feel free to browse the
/scripts folder for examples.
In your code, add the import.
in order to use the functions.
The most useful function is perhaps the following.
def doRLDeconvolutionFromNpArrays(data_np , psf_np ,*, niter=10, method='gpu', useBlockAlgorithm=False, callbkTickFunc=None, resAsUint8 = False)
This will do the Richardson-Lucy deconvolution on the data_np (numpy, 3 dimensional data volume) using the provided PSF data volume, for 10 iterations. GPU method is generally faster but it may fail. If it does fail, the program will automatically use the CPU version that uses the scipy fft package.
Manually building the conda package¶
For this installation, ensure that the conda-build package is installed
conda install conda-build
In windows, simply execute
Or, execute the following command-line to create the installation package.
conda-build --output-folder ./conda-built-packages -c conda-forge conda-recipe
and the conda package will be created in folder conda-built-packages.
Otherwise, navigate to
conda-recipe, and execute on the command-line
conda build .
It will take a while to complete.
Report issues and questions in project's github page, please. Please don't try to send emails as they may be igored or spam-filtered.
- 01 March 2023
- 19 November 2021
- Information not submitted
- Stars: 15
- Forks: 3
- Issues + PRs: 5