I=imread('008.tif');
I1=imnoise(I,'gaussian',0,0.005);%添加高斯噪声
I2=imnoise(I,'gaussian',0,0.005);
I3=imnoise(I,'gaussian',0,0.005);
I4=imnoise(I,'gaussian',0,0.005);
K=imlincomb(0.25,I1,0.25,I3,0.25,I4);%线性组合
subplot(2,2,1);imshow(I);
subplot(2,2,2);imshow(I1);
subplot(2,2,3);imshow(K,[]);
K=double(K);%
m=max(max(K));%
n=min(min(K));
K=uint8((K-n).*255/(m-n));
imwrite(I,'008.bmp','bmp');
imwrite(I1,'008noised.bmp','bmp');
imwrite(K,'008denoised.bmp','bmp');