If I have a PSF and a unconvolve image, Does there exist a expert PSF package to convolve the image?
If you can load your PSF and your image as matrices, then you can compute the convolution with ImageFiltering.jl. More information is given in the docs, but a short demo is given below
using ImageFiltering: reflect, imfilter, centered
psf = centered(rand(5,5)) # Use `centered` to properly configure matrix indices
psf ./= sum(psf) # Normalize psf to unity sum
# ImageFiltering.jl performs image correlation, not convolution, so we reflect
# the PSF before passing into `imfilter`
kernel = reflect(psf)
image = rand(100, 100)
conv_image = imfilter(image, kernel)
# `imfilter` supports many other options here (e.g., boundary conditions), see docs
1 Like
And here’s another example with doggies!
1 Like