Regarding deconvoultion using a point-spread function

I am deconvoluting an image with a simulated point-spread function (PSF), and the top pixels of my image obviously have their intensity decreased. The intensity of the top pixels should, in theory, be the highest. How do I adress this in a reasonable way?

load('data');
image2=deconvwnr(image1,PSF,0.5);

figure(), hold on
subplot(2,1,1), imagesc(image1), axis off, title('Original Image')
subplot(2,1,2), imagesc(image2), axis off, title('Deconvoluted Image')

[Walter Roberson]

I notice the line in the documentation,

• The output image J could exhibit ringing introduced by the discrete Fourier transform used in the algorithm. To reduce the ringing, use I = edgetaper(I,psf) before calling deconvwnr.

Perhaps you could try that?

[Sebastian Daneli]

@Walter Roberson, I'll try it.

Answer 1

[Catalytic]

load('data');
tmp=[flipud(image1);image1];
tmp=deconvwnr(tmp,PSF,0.5);
image2=tmp(end/2+1:end,:);

figure(), hold on
subplot(2,1,1), imagesc(image1), axis off, title('Original Image')
subplot(2,1,2), imagesc(image2), axis off, title('Deconvoluted Image')

●