Deconvolve Implementation

A few weeks ago @jakirkham, myself, and other NVIDIA folks with @tlambert03, @sofroniewn, and other napari/czi folks we chatted about having better support for deconvolve (with Richardson Lucy) within the dask ecosystem. During that discussion we demonstrated a nascent effort to perform RL in fft space with a modified version from [AIAPY's deconvolve](https://aiapy.readthedocs.io/en/v0.2.0/_modules/aiapy/psf/deconvolve.html) . Current version of the demo lives here:https://gist.github.com/quasiben/fd6bc00ae72ef5ae15e7c1e49f92b1ae#file-rl-w-fft-w-dask-ipynb
Some things which @tlambert03 brought which would make a significant impact for image scientists would include options/features to handle the following:

- PSF rotational averaging within the implementation
- Interpolation of the PSF to handle different step sizes between PSF and chunks

We have briefly explored modifications to scikit-image but we have two issues:
- NEP-18/35 work which needs to be handled
- New FFT implementation of Richardson-lucy

What I would like the dask-image to consider is a deconvolution function which supports both host and device implementations.