I want to deconvolve a signal with a point spread function to retrieve a "ground truth" signal. I know that the ground truth is roughly Gaussian and should be 7x7 pixels.
This works when I do deconvolution without noise. However, when I add noise I get this weird periodic semi-gaussian signal that is the same length as the point spread function (47 pixels).
Is it possible to restrict the size of the output to be a 7x7 image?
Here is the code. The variables are attached: