I = imread('cameraman.tif');
I_noise = imnoise(I, 'gaussian', 0.01);

[I2, noise] = wiener2(I_noise,[3 3]);
imshow(I_noise), title('Original Image')
figure, imshow(I2)
title(['Filtered Image: noise variance =',...
        num2str(noise)])