image = imresize(imread('flowersonih35.bmp'),[256,256]); 
h = 1-fspecial('gaussian',256,70)./max(max(fspecial('gaussian',256,70))); 
z = ones(256,256); 
h1(:,:,1) = h; 
h1(:,:,2) = h; 
h1(:,:,3) = h; 
subplot(2,3,1) 
imshow(uint8(h1.*double(image)))