|
object = double(imread('cameraman.tif'));
DC = 0.5;
sigma = 40; [Idc, H] = gfilt(object.*DC,sigma); figure;
subplot(121); imshow(object.*DC,[]); title('Original');
subplot(122); imshow(Idc,[]); title('Blurred');
psf = abs(fftshift(ifft2(H)));
J = deconvwnr(Idc,psf,0.0001);
figure;
subplot(221); imshow(psf,[]); title('PSF');
subplot(222); imshow(Idc,[]); title('DC');
subplot(223); imshow(J,[]); title('Deconvolution');
subplot(224); imshow(object,[]); title('Original');
Reference:
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-26 23:47
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社