A question about PSF deconvolution

Could I obtain the deconvolved image if I have an original image and its PSF? And if I have a deconvolved image and its PSF, how can I obtain the convolved image?

1 Like

Convolution, see for images

is basically. a matrix multiplication, so your second part is relatively easy.
The inverse operation is not uniquely determined, you are facing an inverse problem, so that can not be done in a stable way (without “tricks” like regularisation).

1 Like

Tikhonov regularization is pretty common for this kind of problem, a version of which (if you have some estimate of the spectral density) is called Wiener deconvolution. (These can be implemented with FFTs.) This works pretty well in my experience if the noise is relatively small and the PSF is not too blurry.

2 Likes

Tikhonov sometimes gets a bit blurry, but you can also do TV-regularization (Total variation denoising - Wikipedia) as a regulariser. This usually keeps edges – and is more often used in (pure) denoising and not necessarily deconvolution.

2 Likes

Have you looked at the Julia package: DeconvOptim.jl?

2 Likes

Thank you, this works well :grinning:

1 Like

Author here, if you need some assistance, feel free to ping me :slight_smile:

1 Like

Thank you. The tutorial is clear on website and youtube :grinning: