a = imread('PSWebImageBitMap0001.bmp',  'bmp');
b = a(:,:,2);

b = imnoise(b,'gaussian', 0, 0.3);

PSF = fspecial('gaussian',15,3);
NSR = 0.1;

c = deconvwnr(b,PSF,NSR);

figure(4)
clf
imagesc(c)
colormap(gray)
title('Weiner filter')

